[GRASS-user] Correctly Converting LL to LCC

[Rich Shepard]

On Mon, 28 Dec 2009, Glynn Clements wrote:

>> cs2cs ... +x_0=1312336
> Versus:
>>     false east 1312335.958000 (International feet)
> Versus:
> 	If one uses the internal ellipsoid figure table, ellps=, then because
> 	all a= entries in the table are in meters the user must refer to the
> 	other linear unit items in meters.

   I still have errors in that command line. Changing to +x_0=400000
produces:

[rshepard at salmo /usr4/keypoints]$ cs2cs +proj=longlat +datum=NAD83 +to +proj=lcc +datum=NAD83 +ellps=GRS80 \
> +lat_1=43.0 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=400000 +y_0=0 \
> +nadgrids=WO <<EOF
> 122d30'32.43"W  45d19'19.49"N
> 122d30'17.92"W  45d18'52.45"N
> 122d29'34.08"W  45d18'47.16"N
> EOF
242511.51       398797.24 -0.00
242807.02       397955.05 -0.00
243757.62       397768.55 -0.00

and these are still off. I must still have parameters incorrect.

Rich