[Proj] Lat-Long to Lambert Conformal Conic(2SP) transformation inaccuracy issue

[Gerald I. Evenden]

On Saturday 17 February 2007 2:05 pm, Tonu Lukk wrote:
> Hello all
>
> I'm trying to convert between lag-long and Lambert Conformal Conic (2SP)
> projection.
> The command line is:
> proj +proj=lcc +lat_1=58N +lat_2=59d20'N +lat_0=57d31'03.19415"N
> +lon_0=24E +x_0=500000 +y_0=6375000 +ellps=WGS84 -r
>
> The problem is that when I input the origin's coordinates,
> lat=57d31'03.19415"N lon=24E, I get x=500000 ja y=6375098.83, not x=500000
> y=6375000 as coordinate system's definition suggests - a 98.83m inaccuracy
> in northing value.
> What I'm doing wrong?

The ' and " in the runline are screwing up the shell so that it is not 
interpreting the DMS values properly.  A common error.

The following script gens the expected values:

lproj +proj=lcc +lat_1=58N +lat_2=59d20N +lat_0=57d31\'03.19415N +lon_0=24E 
+x_0=500000 +y_0=6375000 +ellps=WGS84  <<EOF
24e 57d31'03.19415"N
EOF

gives:

500000.00       6375000.00

Note: I removed the -r because I initially got whacko values because I am so 
used to putting longitude first.

-- 
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