Discussion:
[otb-users] otbcli_BundleToPerfectSensor is generating an invalid image
p***@deimos-space.com
2018-10-05 09:00:03 UTC
Permalink
Hi,
I am using otbcli_BundleToPerfectSensor to pan-sharpen a worldview3 image.
It is generating an image which seems to be invalid, when I try to load the
image on QGIS I get std::bad_alloc. I have ubuntu 14.04 server where I
installed otb-bin 5.4.0+dfsg-1~trusty2 and used the following command to
generate the image:
otbcli_BundleToPerfectSensor -ram 4096 -inp input_pan_image
-inxs input_multispectral_image -out pansharpened_image uint16

Am I missing something? Any help will be appreciated.
Pritimoy
--
--
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.
Michel Julien
2018-10-05 09:18:02 UTC
Permalink
Hi,

What file format are you using for input and outputs ? Can gdal read the inputs/output with gdalinfo ?

Regards,

Julien

De : otb-***@googlegroups.com <otb-***@googlegroups.com> De la part de ***@deimos-space.com
Envoyé : vendredi 5 octobre 2018 11:00
À : otb-users <otb-***@googlegroups.com>
Objet : [otb-users] otbcli_BundleToPerfectSensor is generating an invalid image

Hi,
I am using otbcli_BundleToPerfectSensor to pan-sharpen a worldview3 image. It is generating an image which seems to be invalid, when I try to load the image on QGIS I get std::bad_alloc. I have ubuntu 14.04 server where I installed otb-bin 5.4.0+dfsg-1~trusty2 and used the following command to generate the image:
otbcli_BundleToPerfectSensor -ram 4096 -inp input_pan_image -inxs input_multispectral_image -out pansharpened_image uint16

Am I missing something? Any help will be appreciated.
Pritimoy

--
--
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<mailto:otb-***@googlegroups.com>
To unsubscribe from this group, send email to
otb-users+***@googlegroups.com<mailto: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<mailto:otb-users+***@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.
--
--
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.
p***@deimos-space.com
2018-10-05 09:51:28 UTC
Permalink
Hi Julien,
Many thanks for your quick reply. Both Input and Output file
formats are GeoTIFF.gdalinfo can read both inputs and output. gdalinfo for
the output image is as follows:
Driver: GTiff/GeoTIFF
Files: pansharpened_image
Size is 1129, 1257
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (116.442644579287958,40.183680965026390)
Pixel Size = (0.000003177881044,-0.000003178646589)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 116.4426446, 40.1836810) (116d26'33.52"E, 40d11' 1.25"N)
Lower Left ( 116.4426446, 40.1796854) (116d26'33.52"E, 40d10'46.87"N)
Upper Right ( 116.4462324, 40.1836810) (116d26'46.44"E, 40d11' 1.25"N)
Lower Right ( 116.4462324, 40.1796854) (116d26'46.44"E, 40d10'46.87"N)
Center ( 116.4444385, 40.1816832) (116d26'39.98"E, 40d10'54.06"N)
Band 1 Block=1129x1 Type=UInt16, ColorInterp=Gray
NoData Value=0
Band 2 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 3 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 4 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 5 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 6 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 7 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 8 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0

FYI, when I am loading the output on QGIS I am also getting following error
messages:
ERROR 1:pansharpened_image, band 1: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 1: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 2: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 2: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 3: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 3: Failed to compute statistics, no valid
pixels found in sampling.

Kind regards,
Pritimoy
Post by Michel Julien
Hi,
What file format are you using for input and outputs ? Can gdal read the
inputs/output with gdalinfo ?
Regards,
Julien
*Envoyé :* vendredi 5 octobre 2018 11:00
*Objet :* [otb-users] otbcli_BundleToPerfectSensor is generating an
invalid image
Hi,
I am using otbcli_BundleToPerfectSensor to pan-sharpen a worldview3 image.
It is generating an image which seems to be invalid, when I try to load the
image on QGIS I get std::bad_alloc. I have ubuntu 14.04 server where I
installed otb-bin 5.4.0+dfsg-1~trusty2 and used the following command to
otbcli_BundleToPerfectSensor -ram 4096 -inp input_pan_image
-inxs input_multispectral_image -out pansharpened_image uint16
Am I missing something? Any help will be appreciated.
Pritimoy
--
--
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.
<javascript:>
To unsubscribe from this group, send email to
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
For more options, visit https://groups.google.com/d/optout.
--
--
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.
Guillaume Pasero
2018-10-05 13:10:27 UTC
Permalink
p***@deimos-space.com
2018-10-05 14:29:16 UTC
Permalink
Hi,
They overlap correctly and yes you are right the input images already have
a projection. I tried gdal_pansharpen and it has generated a valid image, I
am a bit lost here!

Pritimoy
Hi,
it looks like your output image is empty (filled with 0). It also looks
like your input images already have a projection, when you open them in
QGIS do they overlap correctly ?
pansharpened_image.tif
Regards,
Guillaume
Hi Julien,
Many thanks for your quick reply. Both Input and Output file
formats are GeoTIFF.gdalinfo can read both inputs and output. gdalinfo for
Driver: GTiff/GeoTIFF
Files: pansharpened_image
Size is 1129, 1257
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (116.442644579287958,40.183680965026390)
Pixel Size = (0.000003177881044,-0.000003178646589)
AREA_OR_POINT=Area
INTERLEAVE=PIXEL
Upper Left ( 116.4426446, 40.1836810) (116d26'33.52"E, 40d11' 1.25"N)
Lower Left ( 116.4426446, 40.1796854) (116d26'33.52"E, 40d10'46.87"N)
Upper Right ( 116.4462324, 40.1836810) (116d26'46.44"E, 40d11' 1.25"N)
Lower Right ( 116.4462324, 40.1796854) (116d26'46.44"E, 40d10'46.87"N)
Center ( 116.4444385, 40.1816832) (116d26'39.98"E, 40d10'54.06"N)
Band 1 Block=1129x1 Type=UInt16, ColorInterp=Gray
NoData Value=0
Band 2 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 3 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 4 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 5 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 6 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 7 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 8 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
FYI, when I am loading the output on QGIS I am also getting following
ERROR 1:pansharpened_image, band 1: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 1: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 2: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 2: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 3: Failed to compute statistics, no valid
pixels found in sampling.
ERROR 1:pansharpened_image, band 3: Failed to compute statistics, no valid
pixels found in sampling.
Kind regards,
Pritimoy
Post by Michel Julien
Hi,
What file format are you using for input and outputs ? Can gdal read the
inputs/output with gdalinfo ?
Regards,
Julien
*Envoyé :* vendredi 5 octobre 2018 11:00
*Objet :* [otb-users] otbcli_BundleToPerfectSensor is generating an
invalid image
Hi,
I am using otbcli_BundleToPerfectSensor to pan-sharpen a worldview3
image. It is generating an image which seems to be invalid, when I try to
load the image on QGIS I get std::bad_alloc. I have ubuntu 14.04 server
where I installed otb-bin 5.4.0+dfsg-1~trusty2 and used the following
otbcli_BundleToPerfectSensor -ram 4096 -inp input_pan_image
-inxs input_multispectral_image -out pansharpened_image uint16
Am I missing something? Any help will be appreciated.
Pritimoy
--
--
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 unsubscribe from this group, send email to
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
For more options, visit https://groups.google.com/d/optout.
--
--
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.
<javascript:>
To unsubscribe from this group, send email to
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
For more options, visit https://groups.google.com/d/optout.
--
<http://www.c-s.fr> *Guillaume PASERO*
Responsable technique
*Business Unit ESPACE & GeoInformation - Département Payload Data &
Applications*
*CS SystÚmes d'Information*
Parc de la Grande Plaine - 5, Rue Brindejonc des Moulinais - BP 15872
31506 Toulouse Cedex 05 - FRANCE
--
--
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.
Michel Julien
2018-10-05 15:10:48 UTC
Permalink
Then it should work with OTB as well 
 Can you file a bug report ? And post the result of the ReadImageInfo app on the inputs, it helps to understand how otb sees your data.

De : otb-***@googlegroups.com <otb-***@googlegroups.com> De la part de ***@deimos-space.com
Envoyé : vendredi 5 octobre 2018 16:29
À : otb-users <otb-***@googlegroups.com>
Objet : Re: [otb-users] otbcli_BundleToPerfectSensor is generating an invalid image

Hi,
They overlap correctly and yes you are right the input images already have a projection. I tried gdal_pansharpen and it has generated a valid image, I am a bit lost here!

Pritimoy

On Friday, 5 October 2018 14:10:32 UTC+1, Guillaume Pasero wrote:

Hi,

it looks like your output image is empty (filled with 0). It also looks like your input images already have a projection, when you open them in QGIS do they overlap correctly ?

For clarity, you should also specify the file extension of output image : pansharpened_image.tif
Regards,
Guillaume
On 10/05/2018 11:51 AM, ***@deimos-space.com<javascript:> wrote:
Hi Julien,
Many thanks for your quick reply. Both Input and Output file formats are GeoTIFF.gdalinfo can read both inputs and output. gdalinfo for the output image is as follows:
Driver: GTiff/GeoTIFF
Files: pansharpened_image
Size is 1129, 1257
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (116.442644579287958,40.183680965026390)
Pixel Size = (0.000003177881044,-0.000003178646589)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 116.4426446, 40.1836810) (116d26'33.52"E, 40d11' 1.25"N)
Lower Left ( 116.4426446, 40.1796854) (116d26'33.52"E, 40d10'46.87"N)
Upper Right ( 116.4462324, 40.1836810) (116d26'46.44"E, 40d11' 1.25"N)
Lower Right ( 116.4462324, 40.1796854) (116d26'46.44"E, 40d10'46.87"N)
Center ( 116.4444385, 40.1816832) (116d26'39.98"E, 40d10'54.06"N)
Band 1 Block=1129x1 Type=UInt16, ColorInterp=Gray
NoData Value=0
Band 2 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 3 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 4 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 5 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 6 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 7 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0
Band 8 Block=1129x1 Type=UInt16, ColorInterp=Undefined
NoData Value=0

FYI, when I am loading the output on QGIS I am also getting following error messages:
ERROR 1:pansharpened_image, band 1: Failed to compute statistics, no valid pixels found in sampling.
ERROR 1:pansharpened_image, band 1: Failed to compute statistics, no valid pixels found in sampling.
ERROR 1:pansharpened_image, band 2: Failed to compute statistics, no valid pixels found in sampling.
ERROR 1:pansharpened_image, band 2: Failed to compute statistics, no valid pixels found in sampling.
ERROR 1:pansharpened_image, band 3: Failed to compute statistics, no valid pixels found in sampling.
ERROR 1:pansharpened_image, band 3: Failed to compute statistics, no valid pixels found in sampling.

Kind regards,
Pritimoy

On Friday, 5 October 2018 10:18:05 UTC+1, Julien Michel wrote:
Hi,

What file format are you using for input and outputs ? Can gdal read the inputs/output with gdalinfo ?

Regards,

Julien

De : otb-***@googlegroups.com <otb-***@googlegroups.com> De la part de ***@deimos-space.com
Envoyé : vendredi 5 octobre 2018 11:00
À : otb-users <otb-***@googlegroups.com>
Objet : [otb-users] otbcli_BundleToPerfectSensor is generating an invalid image

Hi,
I am using otbcli_BundleToPerfectSensor to pan-sharpen a worldview3 image. It is generating an image which seems to be invalid, when I try to load the image on QGIS I get std::bad_alloc. I have ubuntu 14.04 server where I installed otb-bin 5.4.0+dfsg-1~trusty2 and used the following command to generate the image:
otbcli_BundleToPerfectSensor -ram 4096 -inp input_pan_image -inxs input_multispectral_image -out pansharpened_image uint16

Am I missing something? Any help will be appreciated.
Pritimoy

--
--
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.
--
--
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<javascript:>
To unsubscribe from this group, send email to
otb-users+***@googlegroups.com<javascript:>
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<javascript:>.
For more options, visit https://groups.google.com/d/optout.

--
<http://www.c-s.fr>[Image supprimée par l'expéditeur.]<http://www.c-s.fr/><http://www.c-s.fr>

Guillaume PASERO
Responsable technique
Business Unit ESPACE & GeoInformation - Département Payload Data & Applications

CS SystÚmes d'Information
Parc de la Grande Plaine - 5, Rue Brindejonc des Moulinais - BP 15872
31506 Toulouse Cedex 05 - FRANCE
+33 561 17 64 21 - ***@c-s.fr<javascript:>

--
--
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<mailto:otb-users+***@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.
--
--
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...