Hello,

I have a raster with different categories (adjacent patches)  but when I 
buffer this raster file I get a buffer for the whole area and not for each 
single patch.

I thought of a loop ( r.mapcalc "select patch cat x" - r.buffer x  - 
r.mapcalc  " set buffer of patch x to null()" ) and then r.patch.

I already did some experiences with 

COUNTER = 1

# while [  $COUNTER -lt $max ]; do
while [  $COUNTER < $max ]; do
        echo The counter is $COUNTER
        
        r.mapcalc "t1 = if(rast_in == $COUNTER,1,null())"
        r.buffer input=t1 output=b.$COUNTER distance=500

        let COUNTER=COUNTER+1 
                
   done

but failed because COUNTER was not recognized and I don't know at all, if 
these bash statements work at all using GRASS.

I welcome any other idea.

Martin
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to