Hi Micha,

Are you sure that your WV2 products are not already rectified ? Because from the errors you describe (and the ReadImageInfo) it looks like they are.

Regards,

Julien

Le 08/12/2017 à 20:36, Micha Silver a écrit :
Hello:
I'm working with some WorldView 2 images, supplied as Ortho-Ready products. When I attempt OrthoRectification I get this error:

2017-12-08 21:25:02 (FATAL): itk::ERROR: ImageToGenericRSOutputParameters(0x2583330): No information in the metadata, please set an image with non empty metadata

Here's my command line:
otbcli_OrthoRectification -io.in " 17FEB19083841-M3DM-057103427020_01_P001_ITM.TIF?&skipcarto=true " -io.out rectified2.tif -elev.geoid egm96.grd -elev.dem ./dem/

Running the above without the "?&skipcarto=true" extensions works (no errors) but the rectified images are off by 150 meters, or so.
Any suggestions? Further image information below

Thanks, Micha


Here's the ImageInfo output :

$ otbcli_ReadImageInfo -in "17FEB19083841-M3DM-057103427020_01_P001.TIF" -keywordlist 1
2017-12-08 21:30:01 (INFO):
Image general information:
    Number of bands : 8
    No data flags : Not found
    Start index :  [0,0]
    Size :  [8577,5599]
    Origin :  [710127,3.44276e+06]
    Spacing :  [2,-2]
    Estimated ground spacing (in meters): [1.99552,2.0055]

Image acquisition information:
    Sensor :
    Image identification number:
    Image projection : PROJCS["WGS 84 / UTM zone 36N",
    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"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",33],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32636"]]

Image default RGB composition:
    [R, G, B] = [0,1,2]

Ground control points information:
    Number of GCPs = 0
    GCPs projection =

Image OSSIM keywordlist (optional):
       Ossim Keyword list:



Output parameters value:
indexx: 0
indexy: 0
sizex: 8577
sizey: 5599
spacingx: 2
spacingy: -2
originx: 710127
originy: 3442757
estimatedgroundspacingx: 1.995521903
estimatedgroundspacingy: 2.005495071
numberbands: 8
sensor:
id:
time:
ullat: 0
ullon: 0
urlat: 0
urlon: 0
lrlat: 0
lrlon: 0
lllat: 0
lllon: 0
town:
country:
rgb.r: 0
rgb.g: 1
rgb.b: 2
projectionref: PROJCS["WGS 84 / UTM zone 36N",
    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"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",33],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32636"]]
keyword:    Ossim Keyword list:


gcp.count: 0
gcp.proj:
gcp.ids:
gcp.info:
gcp.imcoord:
gcp.geocoord:


and finally, ImageInfo with the skipcarto parameter:

$ otbcli_ReadImageInfo -in "17FEB19083841-M3DM-057103427020_01_P001.TIF?&skipcarto=true" -keywordlist 1
2017-12-08 21:32:21 (INFO):
Image general information:
    Number of bands : 8
    No data flags : Not found
    Start index :  [0,0]
    Size :  [8577,5599]
    Origin :  [0.5,0.5]
    Spacing :  [1,1]
    Estimated ground spacing (in meters): [297.377,647.235]

Image acquisition information:
    Sensor :
    Image identification number:

Image default RGB composition:
    [R, G, B] = [0,1,2]

Ground control points information:
    Number of GCPs = 0
    GCPs projection =

Image OSSIM keywordlist (optional):
       Ossim Keyword list:



Output parameters value:
indexx: 0
indexy: 0
sizex: 8577
sizey: 5599
spacingx: 1
spacingy: 1
originx: 0.5
originy: 0.5
estimatedgroundspacingx: 297.3768005
estimatedgroundspacingy: 647.2349854
numberbands: 8
sensor:
id:
time:
ullat: 0
ullon: 0
urlat: 0
urlon: 0
lrlat: 0
lrlon: 0
lllat: 0
lllon: 0
town:
country:
rgb.r: 0
rgb.g: 1
rgb.b: 2
projectionref:
keyword:    Ossim Keyword list:


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 - DSO/SI/2A

--
--
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.

Reply via email to