Hi Lisa, I assume you're meaning you want to eliminate contigious areas whose size is lower than 20 ha, right?
I'd suggest converting the raster into a vector with r.to.vect and then running v.to.db to compute the area of each polygon. The polygons will be have accepted the values of the cells they were made from and if several cells were adjacent with the same value they'll be made into one polygon. Then you can use v.edit to delete all the polygons that don't have a large enough area. After that, a quick trick would be to convert the vector map back into a raster map and use r.mapcalc to divide it by itself so that you've only got values of 1 or NULL and then multiplying that by your land use raster. The result should be all the land use classes you calculated, but the cells that were multiplied with NULL will now also have the value of NULL. Maybe that's not the most elegant way to do things, but it's the one that popped into my mind. Anybody have a better method? Best, Daniel -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Defining-using-MMU-tp5256926p5256983.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
