Ferdinand,
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.
Even
Le 14/07/2023 à 17:36, Ferdinand a écrit :
Hi all,
We do stereo photogrammetry on Pleiades/PNeo images, and in one case
we run our algorithms on pansharpened images.
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).
However, with the change in this PR:
https://github.com/OSGeo/gdal/pull/7373 which was integrated into
release 3.7.0, this is no longer possible. gdal_pansharpen.py now
gives: `RuntimeError: Panchromatic band has no associated geotransform`.
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.
Is there some other way to add a geotransform that does not require
warping the image?
Cheers,
Ferdi
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
http://www.spatialys.com
My software is free, but my time generally not.
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
