<div dir="ltr"><div>Your file prj.adf is</div><div><br></div><div><span style="font-family:monospace">Projection LAMBERT_AZIMUTHAL<br>Units METERS<br>Zunits NO<br>Xshift 0.0<br>Yshift 0.0<br>Parameters 6378137.0 6378137.0<br>6378137.0 /* radius of the sphere of reference<br> 20 0 0.0 /* longitude of center of projection<br> 5 0 0.0 /* latitude of center of projection<br>0.0 /* false easting (meters)<br>0.0 /* false northing (meters)</span></div><div><br></div><div>Looking at the code at<br></div><div></div><div><a href="https://github.com/OSGeo/gdal/blob/158e538d63458414a5348dae12ddf820f62cf30d/ogr/ogr_srs_esri.cpp#L421">https://github.com/OSGeo/gdal/blob/158e538d63458414a5348dae12ddf820f62cf30d/ogr/ogr_srs_esri.cpp#L421</a></div><div>it is reading the parameters in order 2, 1, 3 and 4.<br></div><div><br></div><div>That code is 13 years old, changed due to <a href="https://trac.osgeo.org/gdal/ticket/4302">https://trac.osgeo.org/gdal/ticket/4302</a></div><div><br></div><div>As you can see, the content of the prj.adf is different than the example in the ticket. Yours includes an extra one at the beginning, the radius of the sphere.</div><div><br></div><div>I do not know if there is any documentation from esri about that specific format for LAMBERT_AZIMUTHAL. It would be great.</div><div>And also I am surprised that in 13 years nobody complained.</div><div><br></div><div>If you remove the line with "<span style="font-family:monospace">6378137.0 /* radius of the sphere of reference"</span> it works, but probably it is not what you need.<br></div><div><br></div><div>Cheers,</div><div>Javier.<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 12 May 2024 at 01:07, Liyuneh Gebre via gdal-dev <<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@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"><div>Hi,<br></div><div>I've encountered a discrepancy in the outputs of GDAL and arcpy. While
using gdalinfo on an arcGrid file with a user-defined projection, the
reported latitude and longitude of natural origin and False easting seems to be incorrect. However,
when accessing the same file with arcpy and examining its
spatialReference it seems to align with the prj.adf file and is correct. How
can I resolve this discrepancy? <br></div><div>For your reference link to the file. <br></div><div><a href="https://u.pcloud.link/publink/show?code=XZvlVj0ZWHt5AMCJnjHl4pKMmozOJh8RYVE7" target="_blank">https://u.pcloud.link/publink/show?code=XZvlVj0ZWHt5AMCJnjHl4pKMmozOJh8RYVE7</a></div><div><br></div><div>best regards<br></div><br><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><i><b>Liyuneh Gebre</b></i><i>, </i><i><i><span><i>Associate Researcher</i></span></i></i></div><div><i>Ethiopian Institute of Agricultural Research,EIAR<br></i></div><div dir="ltr"><div><i>Cell Phone:+251 911858155<br></i></div><div><i>linkedIn: <a href="https://www.linkedin.com/in/liyuneh-gebre-b842b440/" target="_blank">https://www.linkedin.com/in/liyuneh-gebre-b842b440/</a><br></i></div><div><i><i>skype:liyenew_1</i></i></div><div><i></i></div><div><i>Addis Ababa, Ethioipa</i></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div>