Wrong st_transform EPSG 3575 North Pole LAEA Europe

[richan]

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

[pramsey]

I'm seeing correct answers in 3.1 and master branches, and wrong answer from 3.0 branch.

[Paul Ramsey <pramsey@…>]

In b9987885/git:

Handle axis order issues in polar stereographic / etc, closes #4844

[Paul Ramsey <pramsey@…>]

In fca92fa/git:

Add tests to support earlier projection axis issues, references #4844

[Paul Ramsey <pramsey@…>]

In b108f45/git:

Add test to back-up earlier projection axis fix, references #4844