[Proj] Reversible Calculation? (cs2cs)

Eric Miller EMiller at dfg.ca.gov
Tue Sep 12 08:28:59 PDT 2006


You've mixed up your X's and Y's in the reverse calculation.

C:\>cs2cs +proj=latlong +datum=WGS84 +to +proj=utm +datum=WGS84
+zone=47
98d36'27.58"E 3d30'17.47"N
456423.28       387405.51 0.00

C:\>cs2cs +proj=utm +datum=WGS84 +zone=47 +to +proj=latlong
+datum=WGS84
456423.28 387405.51
98d36'27.58"E   3d30'17.47"N 0.000
387405.51 456423.28
97d59'8.086"E   4d7'43.261"N 0.000


>>> oon at informatika.org 9/12/2006 5:27:40 AM >>>

Dear All Proj Users & Developers,

I've tried to convert from WGS84 to UTM (zone 47) with cs2cs:
=====
minnie$ ./cs2cs -v +proj=latlong +datum=WGS84 +to +proj=utm 
+datum=WGS84 +zone=47 -r +units=m << EOF
> 3d30'17.47"N 98d36'27.58"E
> EOF
# ---- From Coordinate System ----
#Lat/long (Geodetic)
#
# +proj=latlong +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
# ---- To Coordinate System ----
#Universal Transverse Mercator (UTM)
#       Cyl, Sph
#       zone= south
# +proj=utm +datum=WGS84 +zone=47 +units=m +ellps=WGS84 +towgs84=0,0,0
456423.28       387405.51 0.00
=====

point 3d30'17.47"N 98d36'27.58"E in WGS84 converted into UTM (zone47)
(X,Y,Z) = (456423.28,387405.51,0.00)

Then I want to reverse it back to WGS84.
I only replace the syntax of cs2cs, before "+to" replace with after
"+to".

=====
minnie$ ./cs2cs -v +proj=utm +datum=WGS84 +zone=47 -r +units=m +to 
+proj=latlong +datum=WGS84 << EOF
> 456423.28       387405.51 0.00
> EOF
# ---- From Coordinate System ----
#Universal Transverse Mercator (UTM)
#       Cyl, Sph
#       zone= south
# +proj=utm +datum=WGS84 +zone=47 +units=m +ellps=WGS84 +towgs84=0,0,0
# ---- To Coordinate System ----
#Lat/long (Geodetic)
#
# +proj=latlong +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
97d59'8.086"E   4d7'43.261"N 0.000
minnie:~/Desktop/WGS84/proj-4.4.9/src oon$
=====

The point UTM (zone 47) (X,Y,Z) = (456423.28,387405.51,0.00) no
reversible
into 3d30'17.47"N 98d36'27.58"E in WGS84.
But the result said "97d59'8.086"E   4d7'43.261"N 0.000"


What's wrong with my cs2cs syntax?

btw, why does the Proj library (or cs2cs command) didn't calculate the
zone automatically for UTM?
(the manual describe the UTM Zone: every 6 deg, start from 180W to
180E)

thanks.


-- 
regards,
Oon Arfiandwi
Software Engineer,
Bandung, Indonesia.



_______________________________________________
Proj mailing list
Proj at lists.maptools.org 
http://lists.maptools.org/mailman/listinfo/proj



More information about the Proj mailing list