[PROJ] Migrating to PROJ-9

[Even Rouault]

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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20230530/5b083cc4/attachment.htm>