[GRASS-user] Setting Region from Display

[Micha Silver]

On Wed, 2011-07-27 at 10:13 -0700, Rich Shepard wrote:

> On Wed, 27 Jul 2011, Rich Shepard wrote:
> 
> >  Searching my saved mail list threads I found the answer: v.in.region
> > out=<region_name>. However, I have not found the follow up to this: now
> > that I have defined my limited analytical region, how do I clip the large
> > maps (and their associated attribute data) to this area?
> 
>    Well, that didn't work as expected. When I kill the monitor and start the
> GUI I see only the bottom portion of the zoomed maps. Back to the drawing
> board.
> 


First you'll want to set your analytical region with g.region (use *all
four* parameters, n-= s= e= w=. If you don't specifiy any one of them,
it will stay as it was.)
Now, for rasters, just run:
r.mapcalc <analytical region raster>=<full region rater>
Since r.mapcalc honors the current region, the new rasters will be
limited to the new region

Next you can use the v.in.,region to get a rectangular vector of the
current region, and use that with v.overlay to clip your various line
and polygon vectors like:
v.overlay a_in=streams b_in=analytical_region_vector out=streams_clipped

v.overlay should merge all data columns into the resulting output
vector.
 
-- 
Micha


> 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/20110729/c0bac82b/attachment-0001.html