I don't reproduce this (tried with PROJ master):
$ echo 3356123.5400 1303218.3090 5247430.6050 | PROJ_DEBUG=3 bin/cs2cs
IGNF:RGR92 +to IGNF:REUN47
pj_open_lib(proj.ini): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.ini) - succeeded
pj_open_lib(proj.db): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.db) - succeeded
Using coordinate operation RGR92 (RESEAU GEODESIQUE DE LA REUNION 1992) vers
REUNION-PITON DES NEIGES as a fallback due to lack of more appropriate
operations
3353420.95 1304075.02 5248935.14
$ echo 3356123.5400 1303218.3090 5247430.6050 | PROJ_DEBUG=3 bin/cs2cs
+init=IGNF:RGR92 +to +init=IGNF:REUN47
pj_open_lib(proj.ini): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.ini) - succeeded
pj_open_lib(proj.db): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.db) - succeeded
pj_open_lib(IGNF): call fopen(IGNF) - failed
pj_open_lib(proj.db): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.db) - succeeded
pj_open_lib(IGNF): call fopen(IGNF) - failed
pj_open_lib(proj.db): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.db) - succeeded
pj_open_lib(IGNF): call fopen(IGNF) - failed
pj_open_lib(proj.db): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.db) - succeeded
pj_open_lib(IGNF): call fopen(IGNF) - failed
pj_open_lib(proj.db): call
fopen(/home/even/proj/install-proj-master/share/proj/proj.db) - succeeded
Using coordinate operation RGR92 (RESEAU GEODESIQUE DE LA REUNION 1992) vers
REUNION-PITON DES NEIGES as a fallback due to lack of more appropriate
operations
3353420.95 1304075.02 5248935.14
$ bin/projinfo -s IGNF:RGR92 -t IGNF:REUN47 -o PROJ
Candidate operations found: 2
-------------------------------------
Operation No. 1:
IGNF:TSG355_RGR92_REUN47, RGR92 (RESEAU GEODESIQUE DE LA REUNION 1992) vers
REUNION-PITON DES NEIGES, unknown accuracy, ILE DE LA REUNION
PROJ string:
+proj=helmert +x=-789.99 +y=627.333 +z=89.685 +rx=-0.6072 +ry=-76.8019
+rz=10.568
+s=32.2083 +convention=position_vector
-------------------------------------
Operation No. 2:
INVERSE(IGNF):TSG601_REUN47_RGR92, Inverse of REUNION-PITON DES NEIGES vers
RGR92 (RESEAU GEODESIQUE DE LA REUNION 1992), unknown accuracy, ILE DE LA
REUNION
PROJ string:
+proj=pipeline
+step +inv +proj=helmert +x=789.524 +y=-626.486 +z=-89.904 +rx=0.6006
+ry=76.7946 +rz=-10.5788 +s=-32.3241 +convention=position_vector
$ echo 3356123.5400 1303218.3090 5247430.6050 | bin/cct +proj=helmert
+x=-789.99 +y=627.333 +z=89.685 +rx=-0.6072 +ry=-76.8019 +rz=10.568 +s=32.2083
+convention=position_vector
3353420.9490 1304075.0212 5248935.1441 inf
Le 30/05/2023 à 19:46, Nicolas Bellaiche a écrit :
Hi all,
Is there any reason why the deprecated syntax "+init=" doesn't give
the same result as the new syntax?
echo 3356123.5400 1303218.3090 5247430.6050 | ./cs2cs IGNF:RGR92 +to
IGNF:REUN47
3353420.951304075.02 5248935.14
echo 3356123.5400 1303218.3090 5247430.6050 | ./cs2cs +init=IGNF:RGR92
+to +init=IGNF:REUN47
3356123.541303218.31 5247430.61
Thanks in advance for your answers,
--
http://www.spatialys.com
My software is free, but my time generally not.
_______________________________________________
PROJ mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/proj
