Hi,
It is quite uncommon that the input sensor geometry image has a negative
y spacing, and this sure confuses Orfeo ToolBox, which always assume
positive spacing for y spacing in sensor geometry.
Can you :
- Give insight on your Orfeo ToolBox version and package,
- post the result of gdalinfo command on the panchro image,
- Check that there are no tfw files associated with the panchro tif file
(usually we handle jpeg2000 products, and have troubles with the j2w
files ...). If there is a tfw file, try removing it and ReadImageInfo again.
There are two different solutions depending on the outcome of those tests :
- If the problem comes from the tfw file, removing it or using the
skipcarto extended filename option will do the trick,
- Or you could use gdal_edit.py to reset spacing to a correct value
Julien
Le 13/07/2016 à 04:37, [email protected] a écrit :
Hi everybody,
I have a small problem with the use of the OTB. I use the library
Orfeo in a professional environnement to process Pleiades image using
command line.
We use ADS products at PRIMARY BUNDLE brought to L2A (orthorectified
and atmo corrected). We never had any problems in our workflow until
this week:
- We ordered a picture "INSTANT TASKING" format PRIMARY
- We ran the process, but we had a problem with orthorectification:
the image is upside down (mirror effect) and just above the position
it should occupy on the map.
After going trought the OTB manual and the forum and the metadata of
the image, it seems that this problem could be because "yspacing" is a
negative value ... Another possibility is that the product is
delivered in the format SYSTEM_SENSOR and not SENSOR as usual.
My question is this: is it possible to overcome this problem? Is it
possible to set up a routine that automatically overcomes this problem
in a script, when necessary?
To light my problem, here is the command line used
set
PAN="C:\Users\huert\Desktop\Demo_Instant_Tasking_Event_SO16011807-2-01_DS_PHR1A_201606292324213_FR1_PX_E165S22_0216_01013\FCGC600379740\IMG_PHR1A_P_001\IMG_PHR1A_P_201606292324213_SEN_1847607101-001_R1C1.TIF"
set
MS="C:\Users\huert\Desktop\Demo_Instant_Tasking_Event_SO16011807-2-01_DS_PHR1A_201606292324213_FR1_PX_E165S22_0216_01013\FCGC600379740\IMG_PHR1A_MS_002\IMG_PHR1A_MS_201606292324213_SEN_1847607101-002_R1C1.TIF"
set
OUT="C:\Users\huert\Desktop\Demo_Instant_Tasking_Event_SO16011807-2-01_DS_PHR1A_201606292324213_FR1_PX_E165S22_0216_01013\FCGC600379740\"
set DEM="C:\Users\huert\Desktop\PROVISOIRE\SRTM\S22E165.hgt"
set GEOID="C:\Users\huert\Desktop\SCRIPT\Script
Pléiades\EGM96\OTB\egm96.grd"
set EPSG="4326"
otbcli_OpticalCalibration -in %PAN% -out %PAN%_TOC.TIF uint16 -level
toc -milli 1 -ram 15000
otbcli_OpticalCalibration -in %MS% -out %MS%_TOC.TIF uint16 -level toc
-milli 1 -ram 15000
otbcli_BundleToPerfectSensor -inp %PAN%_TOC.TIF -inxs %MS%_TOC.TIF
-mode phr -out %OUT%BUNDLE.TIF uint16 -ram 15000
otbcli_OrthoRectification -io.in %OUT%BUNDLE.TIF -io.out
%OUT%ORTHO.TIF uint16 -elev.dem %DEM% -elev.geoid %GEOID%
-map.epsg.code %EPSG%
As the result of a Read Image Picture information Panchromatic the bundle:
general image information:
Number of bands: 1
No data flags: Not found
Start index: [0,0]
Size: [21373.20912]
Origin: [1,1]
Spacing: [1, -1]
Estimated ground spacing (in meters): [0.540622,0.519378]
Image information acquisition:
Sensor: PHR 1A
Image identification number: 1847607101-001
Acquisition time: 2016-06-29T23: 24: 00
Country: New Caledonia
Town N├®pou
Image footprint coordinates:
Upper left corner (latitude, longitude) = [-21.3196,165.058]
Upper right corner (latitude, longitude) = [-21.3187,165.169]
Lower left corner (latitude, longitude) = [-21.4169,165.058]
Lower right corner (latitude, longitude) = [-21.4164,165.169]
Image default RGB composition:
[R, G, B] = [0,1,2]
Ground control information points:
Number of GCPs = 0
GCPs projection =
Output parameters value:
indexx: 0
indexy: 0
sizex: 21373
SizeY: 20912
spacingx: 1
spacingy: -1
originX: 1
originY: 1
estimatedgroundspacingx: 0.5406222343
estimatedgroundspacingy: 0.5193778872
numberbands: 1
Sensor: PHR 1A
id: 1847607101-001
time: 2016-06-29T23: 24: 00
ullat: -21.31964302
ullon: 165.057724
urlat: -21.31869698
urlon: 165.1693726
lrlat: -21.41637611
lrlon: 165.1694641
lllat: -21.41688919
lllon: 165.0575409
town: N├®pou
country: New Caledonia
rgb.r: 0
rgb.g: 1
rgb.b: 2
projectionref:
keyword:
gcp.count: 0
gcp.proj:
gcp.ids:
gcp.info:
gcp.imcoord:
gcp.geocoord:
--
--
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 [email protected]
To unsubscribe from this group, send email to
[email protected]
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 [email protected]
<mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout.
--
Julien MICHEL
CNES - DCT/SI/AP - BPI 1219
18, avenue Edouard Belin
31401 Toulouse Cedex 09 - France
Tel: +33 561 282 894 - Fax: +33 561 283 109
--
--
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 [email protected]
To unsubscribe from this group, send email to
[email protected]
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
