[gdal-dev] Importing ESRI (EPSG?) 102008 in GDAL 3.x/PROJ 6

[Even Rouault]

Daniel,
 
> While upgrading some Python software from using GDAL 2.3 to GDAL 3, I've
> come across this snippet of code:
> 
> srs = osr.SpatialReference()
> srs.ImportFromEPSG(102008)
> 
> This used to work, but now fails with GDAL 3 / Proj 6 with "crs not found".

yep, we're more pedantic now, and don't mix EPSG and ESRI codes together.

> Projection code 102008 corresponds to "ESRI:102008" described here, which
> the WKT output describes as "EPSG 102008": https://epsg.io/102008

epsg.io should not be taken as a reference website (even if more up-to-date 
than spatialreference.org). The WKT you see there has been built by some piece 
of software, and is in no way the original ESRI WKT (ESRI WKT has no AUTHORITY 
node)

> Is there an easy way to create an SRS from the code 102008 (perhaps with
> some prefix/suffix) in GDAL 3 / PROJ 6? The projection and its code are
> included in the PROJ source data
> (https://github.com/OSGeo/PROJ/blob/master/data/sql/esri.sql), but I can't
> figure out what combination of "ImportFrom..." function and projection code
> format to use.

srs.SetFromUserInput("ESRI:102008")

> (There is a possibility that the GDAL 2.3 build we have has 102008 added
> into the data files - the history of the custom fork of GDAL we're moving
> away from isn't well documented)

No, stock GDAL 2.4 will also output

$ gdalsrsinfo EPSG:102008

PROJ.4 : +proj=aea +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 
+datum=NAD83 +units=m +no_defs 

OGC WKT :
PROJCS["North_America_Albers_Equal_Area_Conic",
    GEOGCS["GCS_North_American_1983",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS_1980",6378137.0,298.257222101]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Albers_Conic_Equal_Area"],
    PARAMETER["False_Easting",0.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["longitude_of_center",-96.0],
    PARAMETER["Standard_Parallel_1",20.0],
    PARAMETER["Standard_Parallel_2",60.0],
    PARAMETER["latitude_of_center",40.0],
    UNIT["Meter",1.0],
    AUTHORITY["Esri","102008"]]

This is an arguable convenience that is no longer offered.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com