[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