Unable to read projection from aux.xml file

[robolt]

It seems that I am not able to read projection information from aux.xml files.

Looking more in depth it looks like the xml parser is failing to parse the projection information into the tree. Ie. the <WKT> information

I am attaching is one zipped file that does not work. I have another that seems to parse some of the data (LOCAL_CS) but again not the WKT information in PROJCS.

This same method I have been using works fine with just an .aux file.

[robolt]

The sid.aux.xml

[robolt]

The smallest dataset I have at the moment is ~4M, I can zip and email as needed.

[robolt]

another sid.aux.xml

[robolt]

This is one file I have been testing with. Be warned it is large (800M+).

​http://gis3.oit.ohio.gov/ziparchives/IMAGERY/1ftsidmosaic/Lake.sid.zip

[warmerdam]

I have downloaded and review 754219_4.sid.aux.xml and it does not contain conventional .aux.xml coordinate system information. Instead it includes specially formatted XML metadata in the xml:ESRI metadata domain looking like this:

  <Metadata domain="xml:ESRI" format="xml">
    <GeodataXform xsi:type="typens:PolynomialXform" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:typens="http://www.esri.com/schemas/ArcGIS/9.2">
      <PolynomialOrder>1</PolynomialOrder>
      <SpatialReference xsi:type="typens:ProjectedCoordinateSystem">
        <WKT>PROJCS[&quot;NAD_1983_HARN_StatePlane_California_IV_FIPS_0404_Feet&quot;,GEOGCS[&quot;GCS_North_American_1983_HARN&quot;,DATUM[&quot;D_North_American_1983_HARN&quot;,SPHEROID[&quot;GRS_1980&quot;,6378137.0,298.257222101]],PRIMEM[&quot;Greenwich&quot;,0.0],UNIT[&quot;Degree&quot;,0.0174532925199433]],PROJECTION[&quot;Lambert_Conformal_Conic&quot;],PARAMETER[&quot;False_Easting&quot;,6561666.666666666],PARAMETER[&quot;False_Northing&quot;,1640416.666666667],PARAMETER[&quot;Central_Meridian&quot;,-119.0],PARAMETER[&quot;Standard_Parallel_1&quot;,36.0],PARAMETER[&quot;Standard_Parallel_2&quot;,37.25],PARAMETER[&quot;Latitude_Of_Origin&quot;,35.33333333333334],UNIT[&quot;Foot_US&quot;,0.3048006096012192]]</WKT>
        <HighPrecision>true</HighPrecision>
      </SpatialReference>
      <CoeffX xsi:type="typens:ArrayOfDouble">
        <Double>6422614.90434035</Double>
        <Double>-0.9998792879354</Double>
        <Double>-0.03347614294263</Double>
      </CoeffX>
...

I presume this is written by ArcGIS in some circumstances. It is not immediately obvious to me whether GDAL ought to be trying to interprete this metadata or not, but it isn't surprising that it isn't having the desired effect.

Adding Gao Peng in case he has any thoughts.

[warmerdam]

I don't have any plans to add support for this formulation. It is used to model specific things within ArcGIS.