[Proj] ITRF96/UTM 36N conversion

[Gerald I. Evenden]

On Friday 07 August 2009 10:07:19 am fernandohoze wrote:
> Hello,
>
> I would like to convert coordinates in ITRF96 / UTM-36N format (measured
> with DGPS) to latitude and longitude values that I can use with Google
> Maps. I tried to do the conversion using proj but I couldn't do it. What
> command should I use for this?
>
> This is what my coordinates file looks like:
>
> 2   4415699.244828   482321.027219   950.950000   5 | 2 | 1.400
> 3   4415697.320372   482321.307750   951.050000   5 | 2 | 1.400
> 4   4415696.356924   482321.975236   951.050000   4 | 2 | 1.000
>
> This is what I want:
>
> 2   39.8911748   32.7932165
> 3   39.8911574   32.7932198
> 4   39.8911488   32.7932277
>
> I also would like to find the distance (in meters) between two points whose
> ITRF96 / UTM-36N coordinates are known. I googled it but I didn't find
> anything about it. Can you please help me with this?
>
> Thanks.

Using the script:

#!/bin/bash
lproj +proj=utm +ellps=WGS84 +lon_0=-39.1 -I -r -f %.7f <<EOF
4415699.244828   482321.027219   950.950000   5 | 2 | 1.400
4415697.320372   482321.307750   951.050000   5 | 2 | 1.400
4415696.356924   482321.975236
EOF

This is what I get:

-39.2067835     39.8911747   950.950000   5 | 2 | 1.400
-39.2067802     39.8911573   951.050000   5 | 2 | 1.400
-39.2067723     39.8911487

Note: the true CM gets automatically set in this case to 39W.  Note that the 
program had to be told that the x-y input was in reverse order (-r option).  
I also had to guess that the longitude you gave was west as the output would 
not have agreed.

The longitude seems to agree reasonably well but the latitude is real strange!

-- 
The whole religious complexion of the modern world is due
to the absence from Jerusalem of a lunatic asylum.
-- Havelock Ellis (1859-1939) British psychologist