Dylan:
> This is slightly incorrect -- as we are growing the cost surface to
> compute the perimeter, when we should be shrinking it. However it
> appears that r.grow cannot use a negative radius.

create a negative mask, grow that, then invert again:


r.mapcalc "inv_mask = if(isnull(mask), 1, null())"

r.grow in=inv_mask out=inv_mask.grown

r.mapcalc "mask_shrunk = if(isnull(inv_mask.grown), 1, null())"



Hamish




      
____________________________________________________________________________________
You rock. That's why Blockbuster's offering you one month of Blockbuster Total 
Access, No Cost.  
http://tc.deals.yahoo.com/tc/blockbuster/text5.com

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

Reply via email to