<div dir="ltr">Hi Evan, <div><br></div><div>Thanks for the response, I'll give it a go. <br>As you are on the ticket I mentioned, I assume the rationale behind the change was to remove the implicit assumption of overlap? <br>Was there also a subpixel shift then when doing it the "old" way? </div><div><br></div><div>Cheers,</div><div><br></div><div>Ferdi</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 14, 2023 at 10:50 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
  
    
  
  <div>
    <p>Ferdinand,</p>
    <p>you can for example use "gdal_edit.py -ro -a_ullr X1 Y1 X2 Y2
      your.tif" to create a geotransform in a sidecar .aux.xml file.  If
      the pan and ms images have the same extent, you could possibly
      just use -a_ullr 0 0 1 -1 (untested though). This hypothesis is
      not totally true as there's some subpixel shift for those products
      between the pan and ms bands if I remember well, so you might need
      to tune that a bit to get optimal results, but generally assuming
      same extent should give you already a decent result.<br>
    </p>
    <p>Even<br>
    </p>
    <div>Le 14/07/2023 à 17:36, Ferdinand a
      écrit :<br>
    </div>
    <blockquote type="cite">
      
      <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>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>
      <br>
      <fieldset></fieldset>
      <pre>_______________________________________________
gdal-dev mailing list
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
    </blockquote>
    <pre cols="72">-- 
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
  </div>

</blockquote></div>