Le 30/05/2023 à 19:24, Nicolas Bellaiche a écrit :
Ok, i kind of fixed my issue even if that's not all very clear to me.
One of the things that I can't understand is that outdated "+init="
For instance if I try to use
echo 2d20\'11.5\"E46d47\'59.775\"N |./bin/projIGNF:LAMBE -Swithout the
+init= syntaxI get
proj_create: unrecognized format / unknown name
Rel. 9.1.1, December 1st, 2022
<proj>:
projection initialization failure
cause: Invalid PROJ string syntax
program abnormally terminated
Is it normal?
Maybe? the old proj utility works with coordinate operations, actually
only map projections, but not CRS definitions. In PROJ 4.x era, there
was no such distinction. IGNF:LAMBE is a CRS definition now. Furthermore
the proj utility has not been updated to resolve objects by code/names.
(with the +init= prefix, proj errors out with "can't initialize
operations that take non-angular input coordinates" for me. likely some
mismatch between radian vs degree input expectations). I presume it
could be enhanced to work with projected CRS definitions, although it is
not totally clear to me what should be the expectations regarding
input/output axis order.
Even
Sincerely,
N.B.
------------------------------------------------------------------------
*De: *"Even Rouault" <[email protected]>
*À: *"Nicolas Bellaiche" <[email protected]>
*Cc: *"proj" <[email protected]>
*Envoyé: *Mardi 30 Mai 2023 14:43:03
*Objet: *Re: [PROJ] Migrating to PROJ-9
Weird indeed:
and what about trying directly the pipeline with cct :
echo 652760.73700000 7033791.24300000 43.64200000 | cct +proj=pipeline
+step +inv +proj=lcc +lat_0=46.5 +lon_0=3 +lat_1=44 +lat_2=49
+x_0=700000 +y_0=6600000 +ellps=GRS80 +step +proj=cart +ellps=GRS80
you should get:
4070327.2514 166078.9321 4891298.6149 inf
Le 30/05/2023 à 14:35, Nicolas Bellaiche a écrit :
Bonjour M. Rouault,
*
*
Voilà ce que j'obtiens:
*
*
*PROJ9:*
*
*
*Release % PROJ_DEBUG=3 projinfo -s IGNF:LAMB93 -t IGNF:RGF93 -o PROJ*
pj_open_lib(proj.ini): call
fopen(/Users/nbellaiche/dev/sd-externs/macosx-clang/proj-9.1.1/nad/proj.ini)
- succeeded
pj_open_lib(proj.db): call
fopen(/Users/nbellaiche/dev/sd-externs/macosx-clang/proj-9.1.1/nad/proj.db)
- succeeded
Candidate operations found: 1
-------------------------------------
Operation No. 1:
unknown id, Inverse of LAMBERT-93 + Conversion from RGF93
geographiques (dms) to RGF93 cartesiennes geocentriques, 0 m, World
PROJ string:
+proj=pipeline
+step +inv +proj=lcc +lat_0=46.5 +lon_0=3 +lat_1=44 +lat_2=49
+x_0=700000
+y_0=6600000 +ellps=GRS80
+step +proj=cart +ellps=GRS80
* Release % echo 652760.73700000 7033791.24300000 43.64200000 |
cs2cs +init=IGNF:LAMB93 +to +init=IGNF:RGF93*
4070355.05166080.07 4891332.24
* Release % echo $PROJ_DATA*
/Users/nbellaiche/dev/sd-externs/macosx-clang/proj-9.1.1/nad
* Release % echo $PROJ_LIB*
/Users/nbellaiche/dev/sd-externs/macosx-clang/proj-9.1.1/nad
je ne vois pas de différence. Je n'ai pas le début d'une idée de
comment identifier le problème...
N.
------------------------------------------------------------------------
*De: *"Even Rouault" <[email protected]>
<mailto:[email protected]>
*À: *"proj" <[email protected]> <mailto:[email protected]>
*Envoyé: *Vendredi 26 Mai 2023 21:50:27
*Objet: *Re: [PROJ] Migrating to PROJ-9
Ah, and you should avoid using the "+init=" prefix. This is a
deprecated legacy syntax that doesn't enforce axis order of the
CRS definition (without practical consequence for IGNF: CRS since
they use "GIS friendly" axis order)
To help debugging this, please provide the output of:
PROJ_DEBUG=3 projinfo -s IGNF:LAMB93 -t IGNF:RGF93 -o PROJ
On my system I get:
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
Candidate operations found: 1
-------------------------------------
Operation No. 1:
unknown id, Inverse of LAMBERT-93 + Conversion from RGF93
geographiques (dms) to RGF93 cartesiennes geocentriques, 0 m, World
PROJ string:
+proj=pipeline
+step +inv +proj=lcc +lat_0=46.5 +lon_0=3 +lat_1=44 +lat_2=49
+x_0=700000
+y_0=6600000 +ellps=GRS80
+step +proj=cart +ellps=GRS80
Le 26/05/2023 à 21:46, Even Rouault a écrit :
Hi,
I get the expected 4070327.25 166078.93 4891298.61 result with
proj 6.0.0, 7.2.1, 9.2.0 and master
There's no grid involved in this operation, just inverse LCC
projection on the RGF93 datum
I also tried to reproduce your issue if the PROJ < 6 IGNF file
is found in the PROJ resource path(s) (non nominal situation -
this file is not needed any more, and will generally lead to
suboptimal results. Its content is included in the proj.db
file), but even with that I get the expected result. So no
idea why you get wrong results
Even
Le 26/05/2023 à 20:54, Nicolas Bellaiche a écrit :
Hi all,
I'm migrating my C++ code from PROJ4 to PROJ9 and I meet a
worrying issue with the computation. Some of my regress
tests don't pass so I've tested with cs2cs and gets the
same difference:
PROJ4
echo652760.737000007033791.24300000 43.64200000 | cs2cs
+init=IGNF:LAMB93 +to +init=IGNF:RGF93:
4070327.25166078.93 4891298.61
PROJ7 or PROJ9
echo652760.737000007033791.24300000 43.64200000 | cs2cs
+init=IGNF:LAMB93 +to +init=IGNF:RGF93:
4070355.05166080.07 4891332.24
I suspect that PROJ7 and PROJ9 don't read properly the
geoid file but I'm not sure. Someone can help me tracking
and fixing this problem?
Thanks,
Nicolas Bellaiche
Institut Géographique National France
_______________________________________________
PROJ mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/proj
--
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
--
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
--
http://www.spatialys.com
My software is free, but my time generally not.
--
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
