[PROJ] Strange behaviour using PROJ from conda

[Javier Jimenez Shaw]

Hi

I installed today proj via conda, that includes a PROJ 9.4.1.
conda install -c conda-forge proj
...
proj               conda-forge/linux-64::proj-9.4.1-h54d7996_1

I am working with a (time dependent) transformation from ITRF2020 to
ETRS89. (not modifying the db ;)

The problem is that projinfo is suggesting a transformation that makes
sense, but cs2cs is doing a ballpark.

However, if I compile PROJ 9.4.1 on my own, it works fine.

Here are some commands.

In the Conda environment:

$ which cs2cs
/home/jshaw/anaconda3/envs/testproj/bin/cs2cs

$ cs2cs
Rel. 9.4.1, June 1st, 2024
...

$ echo 46 6 0 2024.6 | PROJ_DEBUG=3  cs2cs EPSG:9989 EPSG:4937 -d 9
pj_open_lib(proj.ini): call
fopen(/home/jshaw/anaconda3/envs/testproj/share/proj/proj.ini) - succeeded
pj_open_lib(proj.db): call
fopen(/home/jshaw/anaconda3/envs/testproj/share/proj/proj.db) - succeeded
... (removed for clarity. I can send it if needed) ...
Using coordinate operation Ballpark geographic offset from ITRF2020 to
ETRS89
46.000000000 6.000000000 0.000000000 2024.6

$ projinfo -s EPSG:9989 -t EPSG:4937 -o proj --bbox 5.5,45.5,6.5,46.5
Candidate operations found: 4
-------------------------------------
Operation No. 1:

unknown id, Conversion from ITRF2020 (geog3D) to ITRF2020 (geocentric) +
ITRF2020 to ETRF2020 (2) + Conversion from ETRF2020 (geocentric) to
ETRF2020 (geog2D) + Inverse of ETRS89 to ETRF2020 + Null geographic offset
from ETRS89 (geog2D) to ETRS89 (geog3D), 0.1 m, Europe - onshore and
offshore: Albania; Andorra; Austria; Belgium; Bosnia and Herzegovina;
Bulgaria; Croatia; Cyprus; Czechia; Denmark; Estonia; Faroe Islands;
Finland; France; Germany; Gibraltar; Greece; Hungary; Ireland; Italy;
Kosovo; Latvia; Liechtenstein; Lithuania; Luxembourg; Malta; Moldova;
Monaco; Montenegro; Netherlands; North Macedonia; Norway including Svalbard
and Jan Mayen; Poland; Portugal; Romania; San Marino; Serbia; Slovakia;
Slovenia; Spain; Sweden; Switzerland; United Kingdom (UK) including Channel
Islands and Isle of Man; Vatican City State.

PROJ string:
+proj=pipeline
  +step +proj=axisswap +order=2,1
  +step +proj=unitconvert +xy_in=deg +z_in=m +xy_out=rad +z_out=m
  +step +proj=cart +ellps=GRS80
  +step +proj=helmert +x=0 +y=0 +z=0 +rx=0.002236 +ry=0.013494
+rz=-0.019578 +s=0
        +dx=0 +dy=0 +dz=0 +drx=8.6e-05 +dry=0.000519 +drz=-0.000753 +ds=0
        +t_epoch=2015 +convention=position_vector
  +step +inv +proj=cart +ellps=GRS80
  +step +proj=unitconvert +xy_in=rad +z_in=m +xy_out=deg +z_out=m
  +step +proj=axisswap +order=2,1

-------------------------------------
...

If I run it in my compiled PROJ 9.4.1 out of Conda:

$ echo 46 6 0 2024.6 | PROJ_DEBUG=ON ./cs2cs EPSG:9989 EPSG:4937 -d 9
pj_open_lib(proj.ini): call
fopen(/home/jshaw/work/projs/proj_9.4.1/share/proj/proj.ini) - succeeded
pj_open_lib(proj.db): call
fopen(/home/jshaw/work/projs/proj_9.4.1/share/proj/proj.db) - succeeded
pj_open_lib(eur_nkg_nkgrf17vel.tif): call fopen(eur_nkg_nkgrf17vel.tif) -
failed
pj_open_lib(no_kv_NKGETRF14_EPSG7922_2000.tif): call
fopen(no_kv_NKGETRF14_EPSG7922_2000.tif) - failed
pj_open_lib(eur_nkg_nkgrf03vel_realigned.tif): call
fopen(eur_nkg_nkgrf03vel_realigned.tif) - failed
Using coordinate operation Conversion from ITRF2020 (geog3D) to ITRF2020
(geocentric) + ITRF2020 to ETRF2020 (2) + Conversion from ETRF2020
(geocentric) to ETRF2020 (geog2D) + Inverse of ETRS89 to ETRF2020 + Null
geographic offset from ETRS89 (geog2D) to ETRS89 (geog3D)
45.999994986 5.999991132 -0.001870790 2024.6

Any clue what is happening there?

Thank you.
Javier
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20240904/cc798ecb/attachment.htm>