[GRASS-user] Re: A Script that only selects raster inside a polygon

Marcello Gorini gorini at gmail.com
Fri Sep 24 18:09:56 EDT 2010

>>>Pedro Roma:

>>>I want to develop a script that crops a raster in order to only have the
>>>pixels that are located inside a vector polygon.

>>Marcello Gorini: 

>> There is probably a more elegant solution, but this will do:
>> g.region vect=your_vector
>> v.to.rast input=your_vector output=your_vector_rasterized use=val

> Nikos:

>(comment of minor importance I guess... the above needs "value=" to be fed 
with a value of course).

>> r.mask input=your_vector_rasterized
>> r.mapcalc "clipped_raster=original_raster"
>> r.mask -r


>honestly, it's very nice to see how you stepped up from 0 to 100 within
>zero time :-)

>(...watch out for curves and steep slopes :-p)

Thanks a lot for your comments Nikos! I am just trying to make up for lost
time due to my "Windows point-and-click" past life :-)


View this message in context: http://osgeo-org.1803224.n2.nabble.com/A-Script-that-only-selects-raster-inside-a-polygon-tp5567368p5568572.html
Sent from the Grass - Users mailing list archive at Nabble.com.

More information about the grass-user mailing list