<div dir="ltr"><div>Hi Even,</div><div><br></div><div>I've created ticket <a href="http://trac.osgeo.org/gdal/ticket/5960"><font color="#0066cc">#5960</font></a></div><div><br></div><div>Best Regards,</div><div>Dennis</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 13, 2015 at 2:58 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dennis,<br>
<br>
Could you open a ticket in GDAL Trac about that ?<br>
<br>
Even<br>
<br>
> Hi guys,<br>
><br>
><br>
><br>
> we have a Pleiades image here, the jp2 coming directly from the provider<br>
> with an GML header like this:<br>
><br>
><br>
><br>
> <gml:featureMember><br>
><br>
> <gml:FeatureCollection><br>
><br>
> <gml:featureMember><br>
><br>
> <gml:RectifiedGridCoverage dimension="2"<br>
> gml:id="COVERAGE_DS_PHR1A_201306090940372_FR1_PX_E017N43_0315_01434"><br>
><br>
> <gml:description>Ortho<br>
> DS_PHR1A_201306090940372_FR1_PX_E017N43_0315_01434<br>
> coverage</gml:description><br>
><br>
><br>
> <gml:name>DS_PHR1A_201306090940372_FR1_PX_E017N43_0315_01434</gml:name><br>
><br>
> <gml:rectifiedGridDomain><br>
><br>
> <gml:RectifiedGrid dimension="2"><br>
><br>
> <gml:limits><br>
><br>
><br>
><br>
> <!-- The PLEIADES convention is used: first pixel on the product<br>
> is pixel (1,1) --><br>
><br>
> <gml:GridEnvelope><br>
><br>
> <gml:low>1 1</gml:low><br>
><br>
> <gml:high>8528 7086</gml:high><br>
><br>
> </gml:GridEnvelope><br>
><br>
> </gml:limits><br>
><br>
><br>
><br>
> <!--first axis easting: urn:ogc:def:axis-name:EPSG::3035--><br>
><br>
> <gml:axisName>Easting</gml:axisName><br>
><br>
><br>
><br>
> <!--second axis northing: urn:ogc:def:axis-name:EPSG::3035--><br>
><br>
> <gml:axisName>Northing</gml:axisName><br>
><br>
> <gml:origin><br>
><br>
> <gml:Point srsName="urn:ogc:def:crs:EPSG::3035"><br>
><br>
> <gml:pos>4895767.000000000000000<br>
> 2296945.000000000000000</gml:pos><br>
><br>
> </gml:Point><br>
><br>
> </gml:origin><br>
><br>
> <gml:offsetVector<br>
> srsName="urn:ogc:def:crs:EPSG::3035">2.000000000000000 0</gml:offsetVector><br>
><br>
> <gml:offsetVector srsName="urn:ogc:def:crs:EPSG::3035">0<br>
> -2.000000000000000</gml:offsetVector><br>
><br>
> </gml:RectifiedGrid><br>
><br>
> </gml:rectifiedGridDomain><br>
><br>
> <gml:rangeSet><br>
><br>
> <gml:File><br>
><br>
> <gml:rangeParameters><br>
><br>
> <gml:QuantityExtent uom="int">0 255</gml:QuantityExtent><br>
><br>
> </gml:rangeParameters><br>
><br>
> <gml:fileName>gmljp2://codestream/0</gml:fileName><br>
><br>
> <gml:fileStructure>Record Interleaved</gml:fileStructure><br>
><br>
> </gml:File><br>
><br>
> </gml:rangeSet><br>
><br>
> </gml:RectifiedGridCoverage><br>
><br>
> </gml:featureMember><br>
><br>
> </gml:FeatureCollection><br>
><br>
> </gml:featureMember><br>
><br>
><br>
><br>
> </gml:FeatureCollection><br>
><br>
><br>
><br>
><br>
><br>
> The LAEA projection (EPSG::3035) is defined with having an axis order of Y,<br>
> X, which means GDAL swaps the axis of the geotransform. But in this case<br>
> this gives the wrong result. The JP2 GML defines an X,Y axis order by<br>
> putting in <gml:axisName> tags, which are not yet considered by GDAL,<br>
> neither GDAL 1.11 nor GDAL 2.0.0 beta1.<br>
><br>
> We can solve this issue by setting the GDAL_IGNORE_AXIS_ORIENTATION to<br>
> true, but maybe it would be better if GDAL would consider the<br>
> <gml:axisName><br>
><br>
><br>
> Best Regards,<br>
> Dennis<br>
><br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a><br>
</font></span></blockquote></div><br></div>