[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