[postgis-tickets] [PostGIS] #4922: Axis order wrong in PostGIS 3.1.2 (may releates to #4842)

[PostGIS]

#4922: Axis order wrong in PostGIS 3.1.2 (may releates to #4842)
---------------------+--------------------------------
 Reporter:  tjay     |      Owner:  pramsey
     Type:  defect   |     Status:  new
 Priority:  medium   |  Milestone:  PostGIS PostgreSQL
Component:  postgis  |    Version:  3.1.x
 Keywords:           |
---------------------+--------------------------------
 Since 3.1.2 ST_Transform swaps the coordinates of EPSG:31466 on output. In
 all cases Proj uses proj-datumgrid-1.8.zip.
 cs2cs does not show these swaps...

 Postgis 3.1.2 and Proj 8.0.1
 {{{
 POSTGIS="3.1.2 cbe925d" [EXTENSION] PGSQL="120" GEOS="3.9.1-CAPI-1.14.2"
 SFCGAL="1.3.8" PROJ="8.0.1" GDAL="GDAL 3.2.3, released 2021/04/27"
 LIBXML="2.9.1" LIBJSON="0.11"
 TOPOLOGY RASTER
 }}}

 Unexpected result
 {{{
 test=# select
 ST_AsEWKT(ST_Transform(ST_GeomFromEWKT('SRID=25832;POINT(343806.1686
 5676555.6896)'),31466),4);
                 st_asewkt
 -----------------------------------------
  SRID=31466;POINT(5676286.92 2553378.32)
 (1 Zeile)
 }}}

 Postgis 3.1.2 and Proj 7.2.1
 {{{
 POSTGIS="3.1.2 cbe925d" [EXTENSION] PGSQL="120" GEOS="3.9.1-CAPI-1.14.2"
 SFCGAL="1.3.8" PROJ="7.2.1" GDAL="GDAL 3.2.3, released 2021/04/27"
 LIBXML="2.9.1" LIBJSON="0.11"
 TOPOLOGY RASTER
 }}}

 Unexpected result
 {{{
 test=# select
 ST_AsEWKT(ST_Transform(ST_GeomFromEWKT('SRID=25832;POINT(343806.1686
 5676555.6896)'),31466),4);
                 st_asewkt
 -----------------------------------------
  SRID=31466;POINT(5676286.92 2553378.32)
 (1 Zeile)
 }}}

 Postgis 3.1.1 and Proj 7.2.1
 {{{
 POSTGIS="3.1.1 aaf4c79" [EXTENSION] PGSQL="120" GEOS="3.9.1-CAPI-1.14.2"
 SFCGAL="1.3.8" PROJ="7.2.1" GDAL="GDAL 3.2.3, released 2021/04/27"
 LIBXML="2.9.1"
 LIBJSON="0.11" TOPOLOGY RASTER
 }}}

 Expected result
 {{{
 select
 ST_AsEWKT(ST_Transform(ST_GeomFromEWKT('SRID=25832;POINT(343806.1686
 5676555.6896)'),31466),4);
                 st_asewkt
 -----------------------------------------
  SRID=31466;POINT(2553378.32 5676286.92)
 (1 Zeile)
 }}}

-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/4922>
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.