[postgis-tickets] [PostGIS] #4844: Wrong st_transform EPSG 3575 North Pole LAEA Europe

[PostGIS]

#4844: Wrong st_transform EPSG 3575 North Pole LAEA Europe
-------------------------------+---------------------------
 Reporter:  richan             |      Owner:  pramsey
     Type:  defect             |     Status:  new
 Priority:  medium             |  Milestone:  PostGIS 3.1.2
Component:  postgis            |    Version:  3.0.x
 Keywords:  st_transform 3575  |
-------------------------------+---------------------------
 Wrong transformation only happens in newer PostGIS versions.

 Correct result in PostGIS 2.5.0:

 {{{
 select version();
 --PostgreSQL 11.1 on x86_64-pc-mingw64, compiled by gcc.exe (Rev5, Built
 by MSYS2 project) 4.9.2, 64-bit
 select postgis_full_version();
 --POSTGIS="2.5.0 r16836" [EXTENSION] PGSQL="110" GEOS="3.7.0-CAPI-1.11.0
 3.7.0" PROJ="Rel. 4.9.3, 15 August 2016" GDAL="GDAL 2.2.4, released
 2018/03/19 GDAL_DATA not found" LIBXML="2.7.8" LIBJSON="0.12"
 LIBPROTOBUF="1.2.1" RASTER
 SELECT
 st_asewkt(st_transform(ST_GeomFromEWKT('SRID=3575;POINT(370182.35945313
 -2213980.8213281)'),4326));
 --SRID=4326;POINT(19.4921714866834 69.7902178774315)
 }}}

 Wrong result in PostGIS 3.0.3:

 {{{
 select version();
 --PostgreSQL 13.1 on x86_64-pc-linux-musl, compiled by gcc (Alpine 9.3.0)
 9.3.0, 64-bit
 select postgis_full_version();
 --POSTGIS="3.0.3 0" [EXTENSION] PGSQL="130" GEOS="3.8.1-CAPI-1.13.3"
 PROJ="7.0.1" LIBXML="2.9.10" LIBJSON="0.14" LIBPROTOBUF="1.3.3"
 WAGYU="0.4.3 (Internal)" TOPOLOGY
 SELECT
 st_asewkt(st_transform(ST_GeomFromEWKT('SRID=3575;POINT(370182.35945313
 -2213980.8213281)'),4326));
 --SRID=4326;POINT(-89.4921714866835 69.7902178774315)
 }}}


 Might be related to: #4842

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