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

[Rich Shepard]

On Sun, 11 Mar 2012, Hermann Peifer wrote:

> r.mapcalc small_dem_east = dem_east

Herman,

   The map names have changed, but the issue remains unresolved.

   After setting g.region with both the raster dem and vector basins map (the
drainage basins are the analytical area), I ran,

 	r.mapcalc "dem_basin = dem"

yet the raster map display is not limited to the drainage basins, but to a
square area that just includes the basins. I assume the reason is that the
region is a square. What command could I pass to r.mapcalc that would clip
the raster map to the basins themselves (with some overlap so whole cells of
the raster map form the edge)? Screenshot attached.

Rich
-------------- next part --------------
A non-text attachment was scrubbed...
Name: screenshot.png
Type: image/png
Size: 269170 bytes
Desc: 
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20120314/bfcfd519/screenshot-0001.png