[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