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

Rich Shepard rshepard at appl-ecosys.com
Sun Mar 11 12:35:32 EDT 2012


   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



More information about the grass-user mailing list