<div dir="ltr"><div><br></div>Hi all,<div><br></div><div>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. </div><div>I couldn't find how to make such a vrt programmatically. </div><div>IIRC I could create such a vrt by manually editing the vrt file and inserting <Overview> sections.</div><div><br></div><div>Long story:</div><div>Say I have a few rasters (with overviews) with different extents, assuming that the rasters with the larger extents have lower resolution.<br></div><div>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.</div><div><br></div><div>"almost" working implementation:</div><div><a href="https://github.com/talos-gis/gdalos/blob/master/src/gdalos/make_vrt_with_multiple_extent_overviews.py">https://github.com/talos-gis/gdalos/blob/master/src/gdalos/make_vrt_with_multiple_extent_overviews.py</a> </div><div> <br></div><div>For that I created a list of DataSets objects opened with OpenEx with 

OVERVIEW_LEVEL:</div><div>I want to use these DataSets as a set of overviews for the base DataSet.</div><div><br></div><div>I've tried using:</div><div>vrt_ds = gdal.BuildVRT(str(vrt_filename), ds_list, options=

gdal.BuildVRTOptions(outputBounds=extent))</div><div><br></div><div>But I think it tries to create a mosaic of the input DataSets instead of treating them as overviews.</div><div><br></div><div>Kind regards,</div><div>Idan</div><div><br></div></div>