[GRASS-user] Converting from Lat/Long to Projection -- NOT SOLVED

Glynn Clements glynn at gclements.plus.com
Thu Dec 3 01:17:15 EST 2009


Rich Shepard wrote:

> On Wed, 2 Dec 2009, Rich Shepard wrote:
> 
> > cs2cs +proj=latlong +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=1312336 +y_0=0
> >      +nadgrids=WO
> > <<EOF
> > 45d19'19.49"	122d30'32.43"
> > 45d18'52.45"	122d30'17.92"
> > 45d18'47.16"	122d29'34.08"
> > EOF

1. The default input format has X (easting) followed by Y (northing). 
You need to use the -r switch if the northing comes first.

2. In the absence of a trailing E or W, positive eastings are east of
the prime meridian. You need either a leading "-" or a trailing "W"
for eastings in the western hemisphere.

E.g.:

	cs2cs -r +proj=latlong +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=1312336 +y_0=0 \
	     +nadgrids=WO \
	<<EOF
	45d19'19.49"N	122d30'32.43"W
	45d18'52.45"N	122d30'17.92"W
	45d18'47.16"N	122d29'34.08"W
	EOF

gives:

	1154847.51	398797.24 -0.00
	1155143.02	397955.05 -0.00
	1156093.62	397768.55 -0.00

-- 
Glynn Clements <glynn at gclements.plus.com>


More information about the grass-user mailing list