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
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
