[Proj] cs2cs syntax help

Szilard Albert szilard.albert at gmail.com
Mon Oct 14 23:06:34 EST 2013


Dear list,

very new to this, I have a question that might be very easily answered by
someone familiar with cs2cs.
I have the following Lat/Lon coordinates pair (few thousand pairs actually)

262555.03S 1261004.91E
261455.02S 1261504.90E
262555.03S 1261004.91E
265555.04S 1261404.93E

with datum/spheroid/etc description as below. Zone 51, Australia.



H1400 GEODETIC DATUM SURVEYED   WGS84       WGS84        6378137.000
298.2572229
H1401 TRANSFORMATION PARAMETERS    0.0   0.0   0.0 0.000 0.000 0.000
0.0000000
H1500 GEODETIC DATUM AS PLOTTED WGS84       WGS84        6378137.000
298.2572229
H1501 TRANSFORMATION PARAMETERS    0.0   0.0   0.0 0.000 0.000 0.000
0.0000000
H1600 DATUM SHIFTS                 0.0   0.0   0.0 0.000 0.000 0.000
0.0000000
H1700 VERTICAL DATUM            N/A
H1800 PROJECTION                GEODETIC
H1900 ZONE                      N/A
H2002 ANGULAR UNITS             1    DEGREES
H2100 STANDARD PARALLELS          0 0 0.000N  0 0 0.000N
H2200 CENTRAL MERIDIAN            0 0 0.000E
H2301 GRID ORIGIN                 0 0 0.000N  0 0 0.000E
H2302 GRID COORDINATES AT ORIGIN       0.00W       0.00S


I need to transform these to XY coordinates.

Using

cs2cs +proj=latlong +datum=WGS84 +to  +proj=utm +zone=51 +datum=WGS84 -r
<<EOF
26d25'55" 126d15'04"
EOF

I get

824270.31       2927618.51 0.00


X looks OK, but for the Y values I am expecting values around 6,9xx,xxx ~
7,0xx,xxx

Would you please:

1. Point me towards using the command properly
2. Suggest best format of my text file for batch conversion. Can easily
manipulate/separate columns, would prefer not to have to inset d, ', "
symbols, but would be able if necessary.

Thanks for your help,
regards,
Szilard
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20131015/a4a169f3/attachment.htm 


More information about the Proj mailing list