otbcli_BundleToPerfectSensor shifts the image
Description
After using otbcli_BundleToPerfectSensor, the input location is shifted. No world files involved.
Steps to reproduce
This is what I do (OTB-6.6.0-Linux64):
gdalinfo multi.tif
gdalinfo KAISER2-3000S0_PAN2.tif
otbcli_BundleToPerfectSensor -inp KAISER2-3000S0_PAN2.tif -inxs multi.tif -out KAISER2-3001S0_REFpansharpOTB.tif
gdalinfo KAISER2-3001S0_REFpansharpOTB.tif
gdalinfo for input images:
Corner Coordinates:
Upper Left ( 0.0000000, 0.0000000) ( 1d29’19.48"W, 0d 0’ 0.01"N)
Lower Left ( 0.000, -1000.000) ( 1d29’19.48"W, 0d 0’32.47"S)
Upper Right ( 1000.000, 0.000) ( 1d28’47.23"W, 0d 0’ 0.01"N)
Lower Right ( 1000.000, -1000.000) ( 1d28’47.23"W, 0d 0’32.47"S)
Center ( 500.000, -500.000) ( 1d29’ 3.35"W, 0d 0’16.23"S)
gdalinfo for output:
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 1000.0)
Upper Right ( 1000.0, 0.0)
Lower Right ( 1000.0, 1000.0)
Center ( 500.0, 500.0)
Note that input and output images are displayed in the same region by QGIS, but not by Monteverdi and R
> m <- brick(file.path(imadir,"multi.tif"))
> p <- brick(file.path(imadir,"KAISER2-3000S0_PAN2.tif"))
> out <- brick(file.path(imadir,"KAISER2-3001S0_REFpansharpOTB.tif"))
> extent(m)
class : Extent
xmin : 0
xmax : 1000
ymin : -1000
ymax : 0
> extent(p)
class : Extent
xmin : 0
xmax : 1000
ymin : -1000
ymax : 0
> extent(out)
class : Extent
xmin : 0
xmax : 1000
ymin : 0
ymax : 1000
The whole thing here: https://www.dropbox.com/sh/rh0sf74axtzy6il/AABKe3oyZvyJGhnJVji_Vh2Ua?dl=0
Configuration information
OS OTB-6.6.0-Linux64, OTB-6.6.0-Linux64