Opened 3 years ago
#4954 new defect
ST_Transform error involving NAD27
| Reported by: | turntwo463 | Owned by: | pramsey |
|---|---|---|---|
| Priority: | critical | Milestone: | PostGIS GDAL |
| Component: | postgis | Version: | 3.1.x |
| Keywords: | Cc: |
Description
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)
