[GRASS-user] r.region for intersection with vector map?

[Micha Silver]

Wouldn't the simpledt be to convert the vector to a raster named MASK, then  
r.mapcalc analyticalrast=dem will give you a raster only in the analytical  
region
?

--
Micha Silver
052-366-5918

-----Original message-----
From: Rich Shepard <rshepard at appl-ecosys.com>
To: grass-users at lists.osgeo.org
Sent: Sun, 11 Mar 2012, 18:35:50 GMT+02:00
Subject: [GRASS-user] r.region for intersection with vector map?

   I've used v.overlay for its set intersection ability to save features  
only
within the boundary of another vector map. Can I use r.region to do the same
operation on raster maps?

   Specifically, I have two huge (about 51G each) DEM maps and I need only
the region within the project's analytical boundary. The analytical boundary
is a vector polygon map.

   If I run

 	r.region map=dem_east vect=analytical_boundary

the raster map's region is that of the vector map. But, I don't believe that
this removes all cells beyond that region.

   Advice on how to eliminate the extraneous cells (and clip the raster to
the nearest cell boundaries of the vector map) is greatly appreciated.

TIA,

Rich

_______________________________________________
grass-user mailing list
grass-user at lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

This mail was received via Mail-SeCure System.



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20120311/0381d9da/attachment.html