<html><head><style type="text/css">body {word-wrap: break-word; background-color:#ffffff;}</style></head><body><div style="font-family: sans-serif; font-size: 16px">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<br>?<br><br>--<br>Micha Silver<br><a href="tel:0523665918">052-366-5918</a></div><br><br>-----Original message-----<br><blockquote style="; border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><div style="font-family: sans-serif; font-size: 14px"><b>From: </b>Rich Shepard <rshepard@appl-ecosys.com><b><br>To: </b>grass-users@lists.osgeo.org<b><br>Sent: </b>Sun, 11 Mar 2012, 18:35:50 GMT+02:00<b><br>Subject: </b>[GRASS-user] r.region for intersection with vector map?<br><br></div> I've used v.overlay for its set intersection ability to save features only<br>within the boundary of another vector map. Can I use r.region to do the same<br>operation on raster maps?<br><br> Specifically, I have two huge (about 51G each) DEM maps and I need only<br>the region within the project's analytical boundary. The analytical boundary<br>is a vector polygon map.<br><br> If I run<br><br> r.region map=dem_east vect=analytical_boundary<br><br>the raster map's region is that of the vector map. But, I don't believe that<br>this removes all cells beyond that region.<br><br> Advice on how to eliminate the extraneous cells (and clip the raster to<br>the nearest cell boundaries of the vector map) is greatly appreciated.<br><br>TIA,<br><br>Rich<br><br>_______________________________________________<br>grass-user mailing list<br>grass-user@lists.osgeo.org<br><a href="http://lists.osgeo.org/mailman/listinfo/grass-user">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br><br>This mail was received via Mail-SeCure System.<br><br><br></blockquote></body></html>