>>>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 >Marcello, >honestly, it's very nice to see how you stepped up from 0 to 100 within almost >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 :-) Cheers, Marcello. -- 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. _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
