Discussion:
[otb-users] otbcli_Pansharpening claims "Inputs do not occupy the same physical space!" but it is not true
Agustin Lobo
2018-11-07 17:49:12 UTC
Permalink
I run
otbcli_Pansharpening -inp pan.tif -inxs reflhyper.tif -out
reflPansharp.tif float
and get
2018-11-07 18:29:10 (INFO): Default RAM limit for OTB is 128 MB
2018-11-07 18:29:10 (INFO): GDAL maximum cache size is 798 MB
2018-11-07 18:29:10 (INFO): OTB will use at most 8 threads
2018-11-07 18:29:10 (INFO): No kwl metadata found in file pan.tif
2018-11-07 18:29:10 (INFO): No kwl metadata found in file reflhyper.tif
2018-11-07 18:29:10 (FATAL): itk::ERROR:
SimpleRcsPanSharpeningFusionImageFilter(0x1db9d50): Inputs do not
occupy the same physical space!
InputImage Origin: [1.0000000e+01, -1.0000000e+01], InputImage_1
Origin: [5.0000000e-01, -5.0000000e-01]
Tolerance: 2.0000000e-05
InputImage Spacing: [2.0000000e+01, 2.0000000e+01], InputImage_1
Spacing: [1.0000000e+00, 1.0000000e+00]
Tolerance: 2.0000000e-05

But I can see them one on top of the other in QGIS and according to gdalinfo:
pan.tif:
Origin = (0.000000000000000,0.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=LZW
INTERLEAVE=BAND
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)

reflhyper.tif
Origin = (0.000000000000000,0.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=LZW
INTERLEAVE=BAND
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)

Test images
pan.tif
https://www.dropbox.com/s/cqb11j8ul17c9oq/pan.tif?dl=0
reflhyper.tif
https://www.dropbox.com/s/xz3kuejbrnqjvyb/reflhyper.tif?dl=0

Is it perhaps a wrong message and the proble is that the difference in
resolution is too high?

Thanks
Agus
--
--
Check the OTB FAQ at
http://www.orfeo-toolbox.org/FAQ.html

You received this message because you are subscribed to the Google
Groups "otb-users" group.
To post to this group, send email to otb-***@googlegroups.com
To unsubscribe from this group, send email to
otb-users+***@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/otb-users?hl=en
---
You received this message because you are subscribed to the Google Groups "otb-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to otb-users+***@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
Agustin Lobo
2018-11-08 10:32:27 UTC
Permalink
Just to add some more info: it works with gdal_pansharpen.py
While this alternative is ok for us by now, I think otb offers more
choices of methods
and I keep interested on a solution for this problem.
Agus
Post by Agustin Lobo
I run
otbcli_Pansharpening -inp pan.tif -inxs reflhyper.tif -out
reflPansharp.tif float
and get
2018-11-07 18:29:10 (INFO): Default RAM limit for OTB is 128 MB
2018-11-07 18:29:10 (INFO): GDAL maximum cache size is 798 MB
2018-11-07 18:29:10 (INFO): OTB will use at most 8 threads
2018-11-07 18:29:10 (INFO): No kwl metadata found in file pan.tif
2018-11-07 18:29:10 (INFO): No kwl metadata found in file reflhyper.tif
SimpleRcsPanSharpeningFusionImageFilter(0x1db9d50): Inputs do not
occupy the same physical space!
InputImage Origin: [1.0000000e+01, -1.0000000e+01], InputImage_1
Origin: [5.0000000e-01, -5.0000000e-01]
Tolerance: 2.0000000e-05
InputImage Spacing: [2.0000000e+01, 2.0000000e+01], InputImage_1
Spacing: [1.0000000e+00, 1.0000000e+00]
Tolerance: 2.0000000e-05
Origin = (0.000000000000000,0.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
AREA_OR_POINT=Area
COMPRESSION=LZW
INTERLEAVE=BAND
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)
reflhyper.tif
Origin = (0.000000000000000,0.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)
AREA_OR_POINT=Area
COMPRESSION=LZW
INTERLEAVE=BAND
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)
Test images
pan.tif
https://www.dropbox.com/s/cqb11j8ul17c9oq/pan.tif?dl=0
reflhyper.tif
https://www.dropbox.com/s/xz3kuejbrnqjvyb/reflhyper.tif?dl=0
Is it perhaps a wrong message and the proble is that the difference in
resolution is too high?
Thanks
Agus
--
--
Check the OTB FAQ at
http://www.orfeo-toolbox.org/FAQ.html

You received this message because you are subscribed to the Google
Groups "otb-users" group.
To post to this group, send email to otb-***@googlegroups.com
To unsubscribe from this group, send email to
otb-users+***@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/otb-users?hl=en
---
You received this message because you are subscribed to the Google Groups "otb-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to otb-users+***@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
Loading...