On 04/05/13 20:17, BLANDENIER Lucien wrote:
Dear Grass users,
I'm working with rasters maps of large area and I want to calculate surface
related to specific colours. Before this, I have to replace the surfaces of the
rivers by the surrounding value (or the mean of the surrounding value).
I looked at the manual of r.mapcalc but I didn't found how to do this...
For instance, I want to replace all value between 4 and 7 by the mean of the
value of the 4 nearest pixels.
Is there any way of doing this?
Something like this:
newMap =if(oldMap >= 4 && oldMap <=7,
(oldMap[-1,0]+oldMap[0,-1]+oldMap[0,1]+oldMap[1,0])/4, oldMap)
Which command shall I use then for calculating the surface for each color?
When you say color, I suppose you mean cell values.
r.stats? or is there a command that directly computes the surface?
r.stats -a _does_ directly compute the surface.
Moritz
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
