[gdal-dev] create a vrt with overviews from a raster list

[Idan Miara]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20200709/5aa84518/attachment.html>