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

2012-03-14 Thread Rich Shepard
On Wed, 14 Mar 2012, Hermann Peifer wrote: Define a MASK before running r.mapcalc http://grass.osgeo.org/grass64/manuals/html64_user/r.mask.html Hermann, Duh! Of course! Thanks, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http:/

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

2012-03-14 Thread Hermann Peifer
On 14/03/2012 19:02, Rich Shepard wrote: I assume the reason is that the region is a square. Indeed. What command could I pass to r.mapcalc that would clip the raster map to the basins themselves Define a MASK before running r.mapcalc http://grass.osgeo.org/grass64/manuals/html64_user/r.ma

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

2012-03-12 Thread Rich Shepard
On Sun, 11 Mar 2012, Helmut Kudrnovsky wrote: g.region sets the extent of the mapset (i.e. for raster calculations), but r.region defines the extent of the whole raster map. IMHO g.region is the choice ... Helmut, Thanks for the clarification. Rich ___

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

2012-03-12 Thread Hermann Peifer
On 12/03/2012 10:13, Moritz Lennert wrote: On 12/03/12 05:57, Hermann Peifer wrote: On 12/03/2012 00:16, Helmut Kudrnovsky wrote: but r.region defines the extent of the whole raster map. Why would one want to do this, compared to leaving the extent as detected during import? Use cases mig

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

2012-03-12 Thread Moritz Lennert
On 12/03/12 05:57, Hermann Peifer wrote: On 12/03/2012 00:16, Helmut Kudrnovsky wrote: but r.region defines the extent of the whole raster map. Why would one want to do this, compared to leaving the extent as detected during import? Use cases might be: - the file you import is not georefe

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

2012-03-11 Thread Hermann Peifer
On 12/03/2012 00:16, Helmut Kudrnovsky wrote: but r.region defines the extent of the whole raster map. Why would one want to do this, compared to leaving the extent as detected during import? Hermann ___ grass-user mailing list grass-user@lists.o

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

2012-03-11 Thread Helmut Kudrnovsky
>I suppose that both r.region and g.region will work. g.region and r.region are different tools: g.region - Manages the boundary definitions for the geographic region. (http://grass.osgeo.org/grass64/manuals/html64_user/g.region.html) r.region - Sets the boundary definitions for a raster m

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

2012-03-11 Thread Rich Shepard
On Sun, 11 Mar 2012, Hermann Peifer wrote: Try: g.region vect=analytical_boundary To be honest: until 5 minutes ago, I didn't even know that r.region existed. I thought you made a typo. As I don't understand what r.region is good for: I would use g.region as already said, then followed by:

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

2012-03-11 Thread 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 t

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

2012-03-11 Thread Helmut Kudrnovsky
>As I don't understand what r.region >is good for: r.region sets the extent/boundary definition of the whole raster, but no "clipping" to any region defined by g.region is done - best regards Helmut -- View this message in context: http://osgeo-org.1560.n6.nabble.com/r-region-for-intersectio

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

2012-03-11 Thread Hermann Peifer
On 11/03/2012 17:35, Rich Shepard wrote: If I run r.region map=dem_east vect=analytical_boundary Try: g.region vect=analytical_boundary the raster map's region is that of the vector map. But, I don't believe that this removes all cells beyond that region. Advice on how to eliminate the ex

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

2012-03-11 Thread Helmut Kudrnovsky
> Advice on how to eliminate the extraneous cells (and clip the raster to >the nearest cell boundaries of the vector map) is greatly appreciated. see http://grass.osgeo.org/grass64/manuals/html64_user/g.region.html align=name Adjust region cells to cleanly align with this raster map vect=