'Re: [GRASS-user] Setting Region from Display'

[prev in list] [next in list] [prev in thread] [next in thread] 

List:       grass-user
Subject:    Re: [GRASS-user] Setting Region from Display
From:       Micha Silver <micha () arava ! co ! il>
Date:       2011-07-29 12:20:43
Message-ID: 1311942043.3177.13.camel () MS-SL6
[Download RAW message or body]

[Attachment #2 (multipart/alternative)]

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@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
> 
> This mail was received via Mail-SeCure System.
> 
> 

[Attachment #5 (text/html)]

<!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:

&gt;  Searching my saved mail list threads I found the answer: v.in.region
&gt; out=&lt;region_name&gt;. However, I have not found the follow up to this: now
&gt; that I have defined my limited analytical region, how do I clip the large
&gt; 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 &lt;analytical region raster&gt;=&lt;full region rater&gt;<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>

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[prev in list] [next in list] [prev in thread] [next in thread] 