[Proj] cs2cs syntax help
szilard.albert at gmail.com
Wed Oct 23 23:46:44 EST 2013
I would like some clarifications/assurances regarding my usage of cs2cs.
The task is to convert UTM Zone 53 XY coordinates into Lat/Lon coordinates,
I am using:
*cs2cs -W2 +proj=utm +zone=53 +south +datum=WGS84 +to +proj=latlong
136d38'24.07"E 21d37'54.60"S 0.000*
getting a result that seems fine.
However, the following issues arise:
1. My input data is on the GDA94 datum. When I try +datum=GDA94 (or GDA),
projection initialization failure
cause: unknown elliptical parameter name.*
Am I safe using WGS84, since the difference between these 2 datums is less
than 1 m? (http://www.geoproject.com.au/gda.faq.html).
2. I have heaps of other parameters of my input data - as below - should
any of these be incorporated in my command,
or datum and zone are enough?
Projection UTM MGA Z 53
Central meridian 135 E
False E/N, Scale factor 0.9996 etc
Thank you for your help.
On 15 October 2013 14:06, Szilard Albert <szilard.albert at gmail.com> wrote:
> 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
> H1401 TRANSFORMATION PARAMETERS 0.0 0.0 0.0 0.000 0.000 0.000
> H1500 GEODETIC DATUM AS PLOTTED WGS84 WGS84 6378137.000
> H1501 TRANSFORMATION PARAMETERS 0.0 0.0 0.0 0.000 0.000 0.000
> H1600 DATUM SHIFTS 0.0 0.0 0.0 0.000 0.000 0.000
> 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.
> cs2cs +proj=latlong +datum=WGS84 +to +proj=utm +zone=51 +datum=WGS84 -r
> 26d25'55" 126d15'04"
> I get
> 824270.31 2927618.51 0.00
> X looks OK, but for the Y values I am expecting values around 6,9xx,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,
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Proj