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

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


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

   This does not produce correct values. When it runs the output is:

-5056134.11	5898037.45 	-0.0
-5057005.89	5898228.79	-0.0
-5057306.14	5899149.81	-0.0

   The window values in other maps in the LCC projection are these:

proj:       99
zone:       0
north:      414531.46530409
south:      349974.06580951
east:       1201003.03577829
west:       1126515.54469129
cols:       74487
rows:       64557
e-w resol:  1.00000659
n-s resol:  1.00000619
top:        0.5
bottom:     -0.5
cols3:      74487
rows3:      64557
depths:     1
e-w resol3: 1.00000659
n-s resol3: 1.00000619
t-b resol:  1

   What might I have done incorrectly to produce both negative x values and
values that are so different from that of other maps in the same region and
projection?

Rich


More information about the grass-user mailing list