<div dir="ltr">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.<div><br></div><div>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: <a href="https://apps.epsg.org/api/v1/Change/2020.020/export?format=gml">https://apps.epsg.org/api/v1/Change/2020.020/export?format=gml</a></div><div><br></div><div>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?<div><br></div><div>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.</div><div><br></div><div>When I run projinfo EPSG:2163, I get this message:</div><div><font face="monospace">> projinfo EPSG:2163</font></div><font face="monospace">Warning: object is deprecated<br>Alternative non-deprecated CRS:<br>  EPSG:9311</font><div><br></div><div>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:<br><br></div><div><font face="monospace">>gdalinfo test.tif<br>Driver: GTiff/GeoTIFF<br>Files: test.tif<br>       test.tif.aux.xml<br>Size is 49, 1336<br>Coordinate System is:<br>PROJCRS["NAD27 / US National Atlas Equal Area",<br>    BASEGEOGCRS["NAD27",<br>        DATUM["North American Datum 1927",<br>            ELLIPSOID["Clarke 1866",6378206.4,294.978698213898,<br>                LENGTHUNIT["metre",1]]],<br>        PRIMEM["Greenwich",0,<br>            ANGLEUNIT["degree",0.0174532925199433]],<br>        ID["EPSG",4267]],<br>    CONVERSION["US National Atlas Equal Area",<br>        METHOD["Lambert Azimuthal Equal Area (Spherical)",<br>            ID["EPSG",1027]],<br>        PARAMETER["Latitude of natural origin",45,<br>            ANGLEUNIT["degree",0.0174532925199433],<br>            ID["EPSG",8801]],<br>        PARAMETER["Longitude of natural origin",-100,<br>            ANGLEUNIT["degree",0.0174532925199433],<br>            ID["EPSG",8802]],<br>        PARAMETER["False easting",0,<br>            LENGTHUNIT["metre",1],<br>            ID["EPSG",8806]],<br>        PARAMETER["False northing",0,<br>            LENGTHUNIT["metre",1],<br>            ID["EPSG",8807]]],<br>    CS[Cartesian,2],<br>        AXIS["easting (X)",east,<br>            ORDER[1],<br>            LENGTHUNIT["metre",1]],<br>        AXIS["northing (Y)",north,<br>            ORDER[2],<br>            LENGTHUNIT["metre",1]],<br>    USAGE[<br>        SCOPE["unknown"],<br>        AREA["USA"],<br>        BBOX[15.56,167.65,74.71,-65.69]],<br>    ID["EPSG",9311]]<br>Data axis to CRS axis mapping: 1,2</font><br></div><div><br></div><div>I put some more details on this <a href="https://gis.stackexchange.com/questions/377099/deprecated-crs-epsg2163-gets-reinterpreted-as-another-crs-epsg9311-by-gdal">GIS StackExchange question</a><br><br clear="all"><div>Thank you in advance for any advice!</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Henry Rodman<div>Forest Biometrician</div><div>SilviaTerra</div><div><a href="mailto:henry@silviaterra.com" target="_blank">henry@silviaterra.com</a></div></div></div></div></div></div>