[postgis-tickets] [PostGIS] #5389: ST_Transform broken with PROJ 9.2.0

[PostGIS]

#5389: ST_Transform broken with PROJ 9.2.0
----------------------+---------------------------
  Reporter:  tjay     |      Owner:  pramsey
      Type:  defect   |     Status:  new
  Priority:  medium   |  Milestone:  PostGIS 3.1.9
 Component:  postgis  |    Version:  3.1.x
Resolution:           |   Keywords:  proj
----------------------+---------------------------
Comment (by tjay):

 auditing the DEBUG-Logs it seems some missing (and not needed) Datumgrids
 are differently handled in PostGIS vs cs2cs:

 cs2cs 9.1.1
 {{{
 pj_open_lib(proj.ini): call fopen(/opt/giz/rhel8-pg12-postgis31-930402e-
 gitlab.it.nrw.de/share/proj/proj.ini) - succeeded
 pj_open_lib(proj.db): call fopen(/opt/giz/rhel8-pg12-postgis31-930402e-
 gitlab.it.nrw.de/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/opt/giz/rhel8-pg12-postgis31-930402e-
 gitlab.it.nrw.de/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/opt/giz/rhel8-pg12-postgis31-930402e-
 gitlab.it.nrw.de/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/opt/giz/rhel8-pg12-postgis31-930402e-
 gitlab.it.nrw.de/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/opt/giz/rhel8-pg12-postgis31-930402e-
 gitlab.it.nrw.de/share/proj/proj.db) - succeeded
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 pj_open_lib(de_lgvl_saarland_SeTa2016.tif): call
 fopen(de_lgvl_saarland_SeTa2016.tif) - failed
 pj_open_lib(SeTa2016.gsb): call fopen(SeTa2016.gsb) - failed
 Using coordinate operation Inverse of UTM zone 32N + Inverse of DHDN to
 ETRS89 (8) + 3-degree Gauss-Kruger zone 2
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 }}}

 postgis + proj 9.1.1
 {{{
 pj_open_lib(proj.ini): call fopen(/share/proj/proj.ini) - succeeded
 pj_open_lib(proj.db): call fopen(/share/proj/proj.db) - succeeded
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 pj_open_lib(de_lgvl_saarland_SeTa2016.tif): call
 fopen(de_lgvl_saarland_SeTa2016.tif) - failed
 pj_open_lib(SeTa2016.gsb): call fopen(SeTa2016.gsb) - failed
 Using coordinate operation Inverse of UTM zone 32N + Inverse of DHDN to
 ETRS89 (8) + 3-degree Gauss-Kruger zone 2
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 }}}

 cs2cs 9.2.0
 {{{
 pj_open_lib(proj.ini): call fopen(/share/proj/proj.ini) - succeeded
 pj_open_lib(proj.db): call fopen(/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/share/proj/proj.db) - succeeded
 pj_open_lib(epsg): call fopen(epsg) - failed
 pj_open_lib(proj.db): call fopen(/share/proj/proj.db) - succeeded
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 pj_open_lib(de_lgvl_saarland_SeTa2016.tif): call
 fopen(de_lgvl_saarland_SeTa2016.tif) - failed
 pj_open_lib(SeTa2016.gsb): call fopen(SeTa2016.gsb) - failed
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 Using coordinate operation Inverse of UTM zone 32N + Inverse of DHDN to
 ETRS89 (8) + 3-degree Gauss-Kruger zone 2
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 }}}


 postgis + proj 9.2.0
 {{{
 pj_open_lib(proj.ini): call fopen(/share/proj/proj.ini) - succeeded
 pj_open_lib(proj.db): call fopen(/share/proj/proj.db) - succeeded
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 pj_open_lib(de_lgvl_saarland_SeTa2016.tif): call
 fopen(de_lgvl_saarland_SeTa2016.tif) - failed
 pj_open_lib(SeTa2016.gsb): call fopen(SeTa2016.gsb) - failed
 pj_open_lib(de_adv_BETA2007.tif): call
 fopen(/share/proj/de_adv_BETA2007.tif) - succeeded
 }}}


 Transformation works correctly by adding the missing datumgrid
 (de_lgvl_saarland_SeTa2016.tif) -- short fix for now
 However, it would be important to know why this is happening...
-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/5389#comment:2>
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.