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

Nikos Alexandris nikos.alexandris at felis.uni-freiburg.de
Fri Sep 24 17:23:04 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

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

...forgot (that the) default (is) "value=1" which will work just fine. :-P


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

More information about the grass-user mailing list