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

[Nikos Alexandris]

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).

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

Nikos

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