[Proj] cs2cs syntax help

Szilard Albert szilard.albert at gmail.com
Wed Oct 23 21:46:44 PDT 2013


Dear List,

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,
keeping datum.

I am using:

*cs2cs -W2 +proj=utm +zone=53 +south  +datum=WGS84 +to  +proj=latlong
+datum=WGS84
669734.6 7607027.2
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),
I get
*<cs2cs>:
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?

Datum GDA94
Spheroid GRS80
Semi-major axis
Projection UTM MGA Z 53
Central meridian 135 E
False E/N, Scale factor 0.9996 etc

Thank you for your help.
regards,
Szilard









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
> 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.osgeo.org/pipermail/proj/attachments/20131024/afc67913/attachment.html>


More information about the Proj mailing list