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

Rich Shepard rshepard at appl-ecosys.com
Wed Dec 2 20:28:38 EST 2009


   I've let the following script run for 3+ minutes before killing the
process. I don't understand why it's taking so long for only three points. I
must have an error here, despite following the man page and the example it
contains (except that I've not used the '-r' option).

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

   Suggestions, pointers, clue sticks welcome.

Rich


More information about the grass-user mailing list