[gdal-dev] vrt inverse RPC warping

[thomas bonfort]

Hi gdal,
I have use case where I programmatically need to extract a window of pixels
from a raster, accompanied by the same window (number of pixels and ground
extents) of a third party dem dataset.
When the raster dataset has a projection and a geotransform this is
straightforward as a I can create a warped vrt of the DEM and read from
that.
However, I would like to know if it is possible to do this when my original
raster has no geotransform but instead has RPCs.
I know how to do this by materializing the warped dem (to a file or to
memory) via something like:

gdal_translate -b 1 raster.tif dem-for-raster.tif  #example only, this only
actually used to initialize the dem size and RPC metadata
gdalwarp dem.tif dem-for-raster.tif #this will use dem-for-raster's rpcs

but to avoid unnecessary computations and minimize memory consumption, if
there is a VRT solution which will result in pixels equivalent to what's
inside dem-for-raster.tif

thanks,
Thomas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20251211/ae5f02c4/attachment-0001.htm>