[GRASS-user] distance from extent to origin cost surfaces

[Dylan Beaudette]

On Sunday 30 March 2008, Hamish wrote:
> 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
>
>

As usual Hamish has a clever and concise solution!

-- 
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341