[postgis-tickets] [PostGIS] #4954: ST_Transform error involving NAD27

[PostGIS]

#4954: ST_Transform error involving NAD27
------------------------+--------------------------
 Reporter:  turntwo463  |      Owner:  pramsey
     Type:  defect      |     Status:  new
 Priority:  critical    |  Milestone:  PostGIS GDAL
Component:  postgis     |    Version:  3.1.x
 Keywords:              |
------------------------+--------------------------
 Transforms between epsg:26757 and epsg:4326 produce incorrect results
 which differ from results under earlier versions of PostGIS such as
 version 2.4.

 The query used to reproduce this is:

 SELECT
 ST_asText(ST_Transform(ST_SetSrid(ST_GeomFromText('POINT(454052.363184914
 367472.910555989)'), 26757), 4326))

 Using version 2.4 the result is:

 POINT(-75.5780007645488 39.0090346781401)

 using 3.1.1/3.1.2:

 POINT(-75.57798270944112 39.00894843606639)

 This is off by nearly 10 meters

 SELECT ST_Distance(ST_GeomFromText('POINT(-75.5780007645488
 39.0090346781401)')::geography,
 ST_GeomFromText('POINT(-75.57798270944112 39.00894843606638)')::geography)


 Windows 10 Environment
 GDAL_DATA=d:\PostgreSQL\13\gdal-data
 PROJ_LIB=d:\PostgreSQL\13\share\contrib\postgis-3.1\proj

 PostGIS version:

 POSTGIS="3.1.2 3.1.2" [EXTENSION] PGSQL="130" GEOS="3.9.1-CAPI-1.14.1"
 PROJ="7.1.1" LIBXML="2.9.9" LIBJSON="0.12" LIBPROTOBUF="1.2.1"
 WAGYU="0.5.0 (Internal)" (core procs from "3.1.1 3.1.1" need upgrade)

-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/4954>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.