[gdal-dev] Discrepancy in reported Spatial Reference b/n gdalinfo and arcpy

[Javier Jimenez Shaw]

Your file prj.adf is

Projection    LAMBERT_AZIMUTHAL
Units         METERS
Zunits        NO
Xshift        0.0
Yshift        0.0
Parameters    6378137.0  6378137.0
6378137.0 /* radius of the sphere of reference
  20  0  0.0 /* longitude of center of projection
   5  0  0.0 /* latitude of center of projection
0.0 /* false easting (meters)
0.0 /* false northing (meters)

Looking at the code at
https://github.com/OSGeo/gdal/blob/158e538d63458414a5348dae12ddf820f62cf30d/ogr/ogr_srs_esri.cpp#L421
it is reading the parameters in order 2, 1, 3 and 4.

That code is 13 years old, changed due to
https://trac.osgeo.org/gdal/ticket/4302

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.

I do not know if there is any documentation from esri about that specific
format for LAMBERT_AZIMUTHAL. It would be great.
And also I am surprised that in 13 years nobody complained.

If you remove the line with "6378137.0 /* radius of the sphere of reference"
it works, but probably it is not what you need.

Cheers,
Javier.

On Sun, 12 May 2024 at 01:07, Liyuneh Gebre via gdal-dev <
gdal-dev at lists.osgeo.org> wrote:

> Hi,
> 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?
> For your reference link to the file.
>
> https://u.pcloud.link/publink/show?code=XZvlVj0ZWHt5AMCJnjHl4pKMmozOJh8RYVE7
>
> best regards
>
> --
> *Liyuneh Gebre**, **Associate Researcher*
>
> *Ethiopian Institute of Agricultural Research,EIAR*
>
> *Cell Phone:+251 911858155*
>
> *linkedIn: https://www.linkedin.com/in/liyuneh-gebre-b842b440/
> <https://www.linkedin.com/in/liyuneh-gebre-b842b440/>*
> *skype:liyenew_1*
> *Addis Ababa, Ethioipa*
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20240512/8dd62fdf/attachment.htm>