[gdal-dev] Reprojection issue that works in GDAL 2.4.0 but not 3.0+

[Even Rouault]

On mardi 21 janvier 2020 17:09:07 CET Even Rouault wrote:
> Brian,
> 
> > I am trying to go from (this is the .prj from an ESRI shapefile)
> > PROJCS["USA_Contiguous_Albers_Equal_Area_Conic_USGS_version",GEOGCS["GCS_N
> > or
> > th_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",63781
> > 37.
> > 0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.017453292519943
> > 3]]
> > ,PROJECTION["Albers"],PARAMETER["False_Easting",0.0],PARAMETER["False_Nor
> > thi
> > ng",0.0],PARAMETER["Central_Meridian",-96.0],PARAMETER["Standard_Parallel
> > _1"
> > ,29.5],PARAMETER["Standard_Parallel_2",45.5],PARAMETER["Latitude_Of_Origi
> > n", 23.0],UNIT["Foot_US",0.3048006096012192]]
> > 
> > to EPSG:4326
> 
> Can you share the shapefile, and the command line you use ?

OK, so with the test data you provided to me, I've determined that this is an 
issue in PROJ. Fix submitted per	https://github.com/OSGeo/PROJ/pull/1863

The issue was due to a misidentification to ESRI:102039 which has the same 
name and projection parameters, except that ESRI:102039 uses metre as the 
linear unit, and not foot_us...

As a workaround, you can edit the .prj to modify the PROJCS name to be
"USA_Contiguous_Albers_Equal_Area_Conic_USGS_version_footus"
(addition of "_footus" at the end. Anything would do)

Even

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