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 _______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
