<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
<META NAME="GENERATOR" CONTENT="GtkHTML/3.28.3">
</HEAD>
<BODY>
On Wed, 2011-07-27 at 10:13 -0700, Rich Shepard wrote:
<BLOCKQUOTE TYPE=CITE>
<PRE>
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.
</PRE>
</BLOCKQUOTE>
<BR>
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.)<BR>
Now, for rasters, just run:<BR>
r.mapcalc <analytical region raster>=<full region rater><BR>
Since r.mapcalc honors the current region, the new rasters will be limited to the new region<BR>
<BR>
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:<BR>
v.overlay a_in=streams b_in=analytical_region_vector out=streams_clipped<BR>
<BR>
v.overlay should merge all data columns into the resulting output vector.<BR>
<BR>
-- <BR>
Micha<BR>
<BR>
<BLOCKQUOTE TYPE=CITE>
<PRE>
Rich
_______________________________________________
grass-user mailing list
<A HREF="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</A>
<A HREF="http://lists.osgeo.org/mailman/listinfo/grass-user">http://lists.osgeo.org/mailman/listinfo/grass-user</A>
This mail was received via Mail-SeCure System.
</PRE>
</BLOCKQUOTE>
<BR>
</BODY>
</HTML>