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 > [email protected] > http://lists.osgeo.org/mailman/listinfo/grass-user > > This mail was received via Mail-SeCure System. > >
_______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
