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 _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
