<div dir="ltr"><div>Hi gdal,</div><div>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.</div><div>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.</div><div>However, I would like to know if it is possible to do this when my original raster has no geotransform but instead has RPCs.</div><div>I know how to do this by materializing the warped dem (to a file or to memory) via something like:</div><div><br></div><div>gdal_translate -b 1 raster.tif dem-for-raster.tif  #example only, this only actually used to initialize the dem size and RPC metadata</div><div>gdalwarp dem.tif dem-for-raster.tif #this will use dem-for-raster's rpcs</div><div><br></div><div>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</div><div><br></div><div>thanks,</div><div>Thomas</div></div>