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