3.1-alpine builds projection issue

[cwygoda]

Using a projection like EPSG:5652, projection seems to be broken in 3.1-alpine builds:

select st_astext(st_transform(st_setsrid(st_point(32538051.270000093, 6007341.444999786), 5652), 4326));

yields POINT (Inf Inf), in the images listed as broken below, whereas POINT(9.583477583113051 54.2126738208946) is expected (confirmed in the images listed as working below).

broken images:

• postgis/postgis:13-3.1-alpine
• postgis/postgis:12-3.1-alpine

working images:

• postgis/postgis:12-2.5-alpine
• postgis/postgis:12-3.1
• postgis/postgis:13-3.1

[ImreSamu]

IMHO: This is probably a PROJ issue .. ( related to the Installed PROJ version )

The ~latest PROJ master has a better error message: postgis/postgis:13-master --> ERROR: transform: Point outside of projection domain (2050)

• PROJ_GIT_HASH='dcc2a94eb1c5eda5ede6275b58ae5873e78bb9c0'

psql (13.3 (Debian 13.3-1.pgdg100+1))
Type "help" for help.

postgres=# select st_astext(st_transform(st_setsrid(st_point(32538051.270000093, 6007341.444999786), 5652), 4326));
ERROR:  transform: Point outside of projection domain (2050)

postgres=#  SELECT PostGIS_Full_Version();
                                                                                        postgis_full_version                                                                                        
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 POSTGIS="3.2.0dev 3.1.0rc1-198-g429bf56aa" [EXTENSION] PGSQL="130" GEOS="3.10.0dev-CAPI-1.15.0" PROJ="8.1.0" LIBXML="2.9.4" LIBJSON="0.12.1" LIBPROTOBUF="1.3.1" WAGYU="0.5.0 (Internal)" TOPOLOGY
(1 row)