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 <[email protected]>
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
