Hi all, I have a list of rasters (DataSet objects) that I want to be treated as overviews for the BuildVRT function and not as a mosaic. I couldn't find how to make such a vrt programmatically. IIRC I could create such a vrt by manually editing the vrt file and inserting <Overview> sections.
Long story: Say I have a few rasters (with overviews) with different extents, assuming that the rasters with the larger extents have lower resolution. I'm trying to write a function that takes all these rasters/overviews and creates a vrt using the "best" overviews from all the rasters, so that I have a raster with the maximum extent for each overview level. "almost" working implementation: https://github.com/talos-gis/gdalos/blob/master/src/gdalos/make_vrt_with_multiple_extent_overviews.py For that I created a list of DataSets objects opened with OpenEx with OVERVIEW_LEVEL: I want to use these DataSets as a set of overviews for the base DataSet. I've tried using: vrt_ds = gdal.BuildVRT(str(vrt_filename), ds_list, options= gdal.BuildVRTOptions(outputBounds=extent)) But I think it tries to create a mosaic of the input DataSets instead of treating them as overviews. Kind regards, Idan
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
