[GRASS-user] Re: A Script that only selects raster inside a
nikos.alexandris at felis.uni-freiburg.de
Fri Sep 24 17:23:04 EDT 2010
> > >I want to develop a script that crops a raster in order to only have the
> > >pixels that are located inside a vector polygon.
> > 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