[PROJ] Question about transformation using cs2cs
José Luis García Pallero
jgpallero at gmail.com
Wed Sep 29 07:18:56 PDT 2021
Hello:
I'm using PROJ 7.2.1 from Debian repositories and I have some
questions about cs2cs. If I run
echo 1445748.66 5053500.03 | cs2cs +proj=utm +zone=10 +ellps=intl +to
+proj=utm +zone=10 +ellps=WGS84
I obtain
1445704.69 5053404.02 0.00
so the transformation of UTM coordinates between a point on the
international ellipsoid to WGS84 was performed.
But if I try to do the same with geographical coordinates as
echo 111W 45N | cs2cs +proj=latlong +ellps=intl +to +proj=latlong +ellps=WGS84
I obtain
111dW 45dN 0.000
so no transformation was performed. Why? Setting PROJ_DEBUG I obtain
for this case
pj_open_lib(proj.db): call fopen(/usr/share/proj/proj.db) - succeeded
pj_open_lib(proj.ini): call fopen(/usr/share/proj/proj.ini) - succeeded
pj_ellipsoid - final: a=6378137.000 f=1/298.257, errno=0
pj_ellipsoid - final: ellps=GRS80
Here I understand that only the second declaration of +ellps was taken
into account (althoug EGS84 and GRS80 ellipsoids are not exactly the
same)
Thanks
--
*****************************************
José Luis García Pallero
jgpallero at gmail.com
(o<
/ / \
V_/_
Use Debian GNU/Linux and enjoy!
*****************************************
More information about the PROJ
mailing list