[PROJ] EPSG:2163 is deprecated

[Henry Rodman]

Hello, I recently updated to new versions of PROJ (7.1.1) and GDAL (3.1.3)
on my Ubuntu machine and am experiencing some unexpected behavior when I
work with files projected in EPSG:2163 (US National Atlas Equal Area). I do
not experience the same issues in environments with PROJ 7.0.0 and GDAL
3.0.4 installed.

It seems as if EPSG:2163 has been "deprecated" or "orphaned" but I cannot
find any information about how or why this happened. I found this request
from Even Rouault with a recent search:
https://apps.epsg.org/api/v1/Change/2020.020/export?format=gml

I use this projection quite often and would like to know if there is
anything I can do to resuscitate EPSG:2163! Maybe there is some parameter I
can set locally that would override the automatic reassignment to
EPSG:9311? Or maybe the CRS needs updated documentation or something like
that?

Either GDAL or PROJ (I don't know which) is attempting to redefine the
projection of my files as EPSG:9311 which I believe has a different datum.

When I run projinfo EPSG:2163, I get this message:
> projinfo EPSG:2163
Warning: object is deprecated
Alternative non-deprecated CRS:
  EPSG:9311

When I run gdalinfo on a file projected in EPSG:2163 it looks as if the
projection has been re-interpreted as EPSG:9311 already:

>gdalinfo test.tif
Driver: GTiff/GeoTIFF
Files: test.tif
       test.tif.aux.xml
Size is 49, 1336
Coordinate System is:
PROJCRS["NAD27 / US National Atlas Equal Area",
    BASEGEOGCRS["NAD27",
        DATUM["North American Datum 1927",
            ELLIPSOID["Clarke 1866",6378206.4,294.978698213898,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4267]],
    CONVERSION["US National Atlas Equal Area",
        METHOD["Lambert Azimuthal Equal Area (Spherical)",
            ID["EPSG",1027]],
        PARAMETER["Latitude of natural origin",45,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",-100,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["unknown"],
        AREA["USA"],
        BBOX[15.56,167.65,74.71,-65.69]],
    ID["EPSG",9311]]
Data axis to CRS axis mapping: 1,2

I put some more details on this GIS StackExchange question
<https://gis.stackexchange.com/questions/377099/deprecated-crs-epsg2163-gets-reinterpreted-as-another-crs-epsg9311-by-gdal>

Thank you in advance for any advice!

-- 
Henry Rodman
Forest Biometrician
SilviaTerra
henry at silviaterra.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20201105/33e68c32/attachment.html>