[GRASS-user] Reprojecting Vector Map
Rich Shepard
rshepard at appl-ecosys.com
Tue Nov 24 14:20:21 EST 2009
From the US Census TIGER files I imported the 'edges' (streets) using
v.in.ogr and using the projection information in the *.prj file. This worked
well, except that 'g.region -p' shows lat/lon geographic coordinates:
GRASS 6.4.0svn (streets):/usr4/grassbase > g.region -p
projection: 3 (Latitude-Longitude)
zone: 0
datum: nad83
ellipsoid: grs80
north: 45:27:41.8788N
south: 44:53:05.766N
west: 122:52:05.2068W
east: 121:39:06.1704W
nsres: 0:34:36.1128
ewres: 1:12:59.0364
rows: 1
cols: 1
cells: 1
I changed pwd to PERMANENT and ran g.setproj. Specified lcc as the
projection and entered all the relevant data, then changed back to my mapset
and ran 'g.region -p'. No change. Running 'g.proj -w' tells me that I've not
made any difference; it's still in degrees:
GRASS 6.4.0svn (streets):/usr4/grassbase > g.proj -w
PROJCS["Lambert Conformal Conic",
GEOGCS["grs80",
DATUM["North_American_Datum_1983",
SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",43],
PARAMETER["standard_parallel_2",45.5],
PARAMETER["latitude_of_origin",41.75],
PARAMETER["central_meridian",-120.5],
PARAMETER["false_easting",1312336],
PARAMETER["false_northing",0],
UNIT["degree",1]]
My other vector layers are in LCC with boundaries in, I believe, metres.
What step have I missed with this streets layer?
Rich
More information about the grass-user
mailing list