<div dir="ltr"><div>The definition of the CRS is in the .prj file of the shapefile, as plain text. In this case:<br></div><div><br></div><div>GEOGCS["WGS84(DD)", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]]</div><div><br></div><div>That formal definition of EPSG:4326 is northing-easting, and that WKT is saying easting-northing. So there is an inconsistency in the definition of the system, because at the end it says that the authority is EPSG:4326<br></div><div>You can either remove the axis or the final authority (the error message says that proj is removing the authority ;)<br></div><div><br></div><div>The axis order I think is not used to read the coordinates in a shp file, so the outputs should be the same.<br></div><div><br></div><div>The actual definition of EPSG:4326 in ESRI (from <a href="https://spatialreference.org/ref/epsg/4326/esriwkt.txt">https://spatialreference.org/ref/epsg/4326/esriwkt.txt</a>)</div><div>GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]</div><div><br></div><div>In case anybody needs an equivalent to EPSG:4326 but easting-northing in other context, there is <a href="https://spatialreference.org/ref/ogc/CRS84/">https://spatialreference.org/ref/ogc/CRS84/</a></div><div><br></div><div>Cheers,</div><div>Javier<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 28 Mar 2024 at 03:06, Simon Eves via PROJ <<a href="mailto:proj@lists.osgeo.org">proj@lists.osgeo.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">When trying to import this Shapefile bundle into our system with GDAL...<div><br clear="all"><div><a href="https://drive.google.com/file/d/13Mwhnugcy8HrnGlc81gyroxnp9s15Skt/view?usp=sharing" target="_blank">https://drive.google.com/file/d/13Mwhnugcy8HrnGlc81gyroxnp9s15Skt/view?usp=sharing</a><br></div><div><br></div><div>...we get an <b>OGRERR_CORRUPT_DATA</b> from <b>OGRSpatialReference::Validate()</b> on the SR obtained from the first feature.</div><div><br></div><div>Tracing through the GDAL and PROJ code, this appears to be because the PROJ import of that CRS results in the following WKT Import Error from PROJ...</div><div><br></div><div><i>Coordinate system of GeographicCRS in the WKT definition is different from the one of the authority. Unsetting the identifier to avoid confusion<br></i></div><div><br></div><div>...and any WKT Import Error makes GDAL report that error code.</div><div><br></div>The PROJ code throws this in <b>WKTParser::Private::buildGeodeticCRS</b> (io.cpp:3229 in 9.3.0) which is old code (although the second branch which can throw the same error was seemingly added between 8.2 and 9.3 in <a href="https://github.com/OSGeo/PROJ/pull/3274" target="_blank">https://github.com/OSGeo/PROJ/pull/3274</a> ...but it's the first branch that's triggering in this case).</div><div><br></div><div><b>ogrinfo</b> on the Shapefile reports no errors or warnings, although it does report a CRS ID of 9001 even though it reports WGS84 elsewhere in the CRS block.</div><div><br></div><div>I was able to trivially <b>ogr2ogr</b> it to GeoJSON, also with no errors or warnings, and that imports with no issue.</div><div><br></div><div><span class="gmail_signature_prefix">Is there an actual problem (a bad file) here, or are we doing something wrong?</span></div><div><span class="gmail_signature_prefix"><br></span></div><div><span class="gmail_signature_prefix">Thanks, as ever, in advance.</span></div><div><span class="gmail_signature_prefix"><br></span></div><div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr">Simon Eves<div>Senior Rendering Engineer</div><div>+1 (415) 902-1996<br><a href="mailto:simon.eves@heavy.ai" target="_blank">simon.eves@heavy.ai</a></div><div><span><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><a href="http://www.heavy.ai" target="_blank"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(17,85,204);background-color:transparent;vertical-align:baseline;white-space:pre-wrap"><span style="border:medium;display:inline-block;overflow:hidden;width:156px;height:29px"><img width="156" height="29" style="margin-left: 0px; margin-top: 0px;"></span></span></a></p></span></div></div></div></div></div>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</blockquote></div>