Added the ability for GDAL overviews to be recursively searched.

plugins/input/gdal/gdal_featureset.cpp

@@ -131,12 +131,22 @@ void gdal_featureset::find_best_overview(int bandNumber,
                                          int & current_height) const
 {
     GDALRasterBand * band = dataset_.GetRasterBand(bandNumber);
+    find_best_overview(band, ideal_width, ideal_height, current_width, current_height);
+}
+
+void gdal_featureset::find_best_overview(GDALRasterBand * band,
+                                         int ideal_width,
+                                         int ideal_height,
+                                         int & current_width,
+                                         int & current_height) const
+{
     int band_overviews = band->GetOverviewCount();
     if (band_overviews > 0)
     {
         for (int b = 0; b < band_overviews; b++)
         {
             GDALRasterBand * overview = band->GetOverview(b);
+            find_best_overview(overview, ideal_width, ideal_height, current_width, current_height);
             int overview_width = overview->GetXSize();
             int overview_height = overview->GetYSize();
             if ((overview_width < current_width ||

plugins/input/gdal/gdal_featureset.hpp

@@ -78,6 +78,12 @@ private:
                             int & current_width,
                             int & current_height) const;
     
+    void find_best_overview(GDALRasterBand * band,
+                            int ideal_width,
+                            int ideal_height,
+                            int & current_width,
+                            int & current_height) const;
+
     mapnik::feature_ptr get_feature(mapnik::query const& q);
     mapnik::feature_ptr get_feature_at_point(mapnik::coord2d const& p);
     GDALDataset & dataset_;