<div dir="ltr"><div><div><div><div><div><div>Dear List,<br><br></div>I would like some clarifications/assurances regarding my usage of cs2cs.<br></div>The task is to convert UTM Zone 53 XY coordinates into Lat/Lon coordinates, keeping datum.<br>
<br></div>I am using:<br><br><b><span style="color:rgb(76,17,48)">cs2cs -W2 +proj=utm +zone=53 +south +datum=WGS84 +to +proj=latlong +datum=WGS84<br>669734.6 7607027.2<br>136d38'24.07"E 21d37'54.60"S 0.000</span></b><br>
<br></div>getting a result that seems fine.<br><br></div>However, the following issues arise:<br><br></div>1. My input data is on the GDA94 datum. When I try +datum=GDA94 (or GDA), I get <br><div><div><i><cs2cs>: <br>
projection initialization failure<br>cause: unknown elliptical parameter name.</i><br><br></div><div>Am I safe using WGS84, since the difference between these 2 datums is less than 1 m? (<a href="http://www.geoproject.com.au/gda.faq.html">http://www.geoproject.com.au/gda.faq.html</a>).<br>
<br></div><div>2. I have heaps of other parameters of my input data - as below - should any of these be incorporated in my command, <br></div><div>or datum and zone are enough?<br></div><div><br></div><div>Datum GDA94<br>
</div><div>Spheroid GRS80<br></div><div>Semi-major axis <br>Projection UTM MGA Z 53 <br>Central meridian 135 E <br>False E/N, Scale factor 0.9996 etc<br><br></div><div>Thank you for your help.<br>regards,<br>Szilard<br><br>
</div><div><br><br></div><div><br><br><br><div><div><div> <br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 15 October 2013 14:06, Szilard Albert <span dir="ltr"><<a href="mailto:szilard.albert@gmail.com" target="_blank">szilard.albert@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Dear list, <br><br>very new to this, I have a question that might be very easily answered by someone familiar with cs2cs.<br>
</div>I have the following Lat/Lon coordinates pair (few thousand pairs actually)<br>
<br>262555.03S 1261004.91E<br>261455.02S 1261504.90E<br>262555.03S 1261004.91E<br>265555.04S 1261404.93E<br><br></div>with datum/spheroid/etc description as below. Zone 51, Australia.<br><div><br><br><br>H1400 GEODETIC DATUM SURVEYED WGS84 WGS84 6378137.000 298.2572229<br>
H1401 TRANSFORMATION PARAMETERS 0.0 0.0 0.0 0.000 0.000 0.000 0.0000000<br>H1500 GEODETIC DATUM AS PLOTTED WGS84 WGS84 6378137.000 298.2572229<br>H1501 TRANSFORMATION PARAMETERS 0.0 0.0 0.0 0.000 0.000 0.000 0.0000000<br>
H1600 DATUM SHIFTS 0.0 0.0 0.0 0.000 0.000 0.000 0.0000000<br>H1700 VERTICAL DATUM N/A<br>H1800 PROJECTION GEODETIC<br>H1900 ZONE N/A<br>H2002 ANGULAR UNITS 1 DEGREES<br>
H2100 STANDARD PARALLELS 0 0 0.000N 0 0 0.000N<br>H2200 CENTRAL MERIDIAN 0 0 0.000E<br>H2301 GRID ORIGIN 0 0 0.000N 0 0 0.000E<br>H2302 GRID COORDINATES AT ORIGIN 0.00W 0.00S<br>
<br><br></div><div>I need to transform these to XY coordinates.<br><br></div><div>Using <br><br><span style="color:rgb(0,0,255)">cs2cs +proj=latlong +datum=WGS84 +to +proj=utm +zone=51 +datum=WGS84 -r <<EOF</span><br style="color:rgb(0,0,255)">
<span style="color:rgb(0,0,255)">26d25'55" 126d15'04"</span><br style="color:rgb(0,0,255)"></div><div><span style="color:rgb(0,0,255)">EOF</span><br><br></div><div>I get<br><br><span style="color:rgb(0,0,255)">824270.31 2927618.51 0.00</span><br>
<br><br></div><div>X looks OK, but for the Y values I am expecting values around 6,9xx,xxx ~ 7,0xx,xxx<br><br></div><div>Would you please:<br><br></div><div>1. Point me towards using the command properly<br></div><div>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.<br>
<br></div><div>Thanks for your help,<br></div><div>regards,<br>Szilard<br><br></div><div><br></div></div>
</blockquote></div><br></div></div></div></div>