Wouldn't the simpledt be to convert the vector to a raster named MASK, then
r.mapcalc analyticalrast=dem will give you a raster only in the analytical
region
?
--
Micha Silver
052-366-5918
-----Original message-----
From: Rich Shepard <[email protected]>
To: [email protected]
Sent: Sun, 11 Mar 2012, 18:35:50 GMT+02:00
Subject: [GRASS-user] r.region for intersection with vector map?
I've used v.overlay for its set intersection ability to save features
only
within the boundary of another vector map. Can I use r.region to do the same
operation on raster maps?
Specifically, I have two huge (about 51G each) DEM maps and I need only
the region within the project's analytical boundary. The analytical boundary
is a vector polygon map.
If I run
r.region map=dem_east 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 extraneous cells (and clip the raster to
the nearest cell boundaries of the vector map) is greatly appreciated.
TIA,
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
