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

[Rich Shepard]

On Sun, 11 Mar 2012, Helmut Kudrnovsky wrote:

> align=name
>    Adjust region cells to cleanly align with this raster map
> to set the region boundary to your analytical vector: i.e.
> g.region -a rast=dem_east vect=analytical_boundary align=dem_east

Helmut,

   I read about align on the man page but thought it required a raster map
different from the one being clipped. This is good to know!

> then r.mapcalc "dem_east_vect_analytical_boundary = dem_east"

   Ah-ha! I wondered if r.mapcalc was the appropriate tool.

Thanks,

Rich