<div dir="ltr">Hi all, <br><br>We do stereo photogrammetry on Pleiades/PNeo images, and in one case we run our algorithms on pansharpened images. <br><br>Our workflow therefore was that we would pan-sharpen the images and then run our stereo processing (using the panchromatic RPC model for the pansharpened image). <br><br>However, with the change in this PR: <a href="https://github.com/OSGeo/gdal/pull/7373" target="_blank">https://github.com/OSGeo/gdal/pull/7373</a> which was integrated into release 3.7.0, this is no longer possible. gdal_pansharpen.py now gives: `<span style="color:rgb(31,35,40);font-family:ui-monospace,SFMono-Regular,"SF Mono",Menlo,Consolas,"Liberation Mono",monospace;font-size:11.9px">RuntimeError: Panchromatic band has no associated geotransform`. <br></span><br>What is the recommended workflow for this now? We can't orthorectify the images first, as we need the raw row/column information in the image in order for the photogrammetry processing to work. <br><br>Is there some other way to add a geotransform that does not require warping the image? <div><br></div><div><br></div><div>Cheers,</div><div><br></div><div>Ferdi </div></div>