<div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">AFAICT, r.neighbors ignores NULL values and assigns the new value from</div>
the surrounding non-NULL values. You could then patch your original<br>
map with the output of r.neighbors.<br></blockquote><div><br></div><div>I have considered r.neighbors but I need to apply the filter only to some specific categories (border cells in my case). A MASK wouldn't solve my problem, because it would mask the neighbors...</div>
<div><br></div><div>giovanni</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
HTH,<br>
<br>
Markus M<br>
<div class="im"><br>
><br>
> giovanni<br>
><br>
><br>
> 2012/6/19 Marcello Gorini <<a href="mailto:gorini@gmail.com">gorini@gmail.com</a>><br>
>><br>
>><br>
>><br>
>> G. Allegri wrote:<br>
>><br>
>>> I need to assign values to a the cells on the "border" of a raster. The<br>
>>> inside and the outside are distinguished by having or not having null values<br>
>>> assigned.<br>
>>> I also need to keep the other cell values (internals) untouched.<br>
>>><br>
>><br>
>>  Hey,<br>
>><br>
>> Check out r.grow. I am pretty sure you can modify the example given in the<br>
>> manual to find the borders of your raster. Something like:<br>
>><br>
>> # creates an inverted raster from your raster<br>
>> > r.mapcalc "raster_inverted=if(isnull(raster,1,null())"<br>
>> # grow this inverted raster by one cell<br>
>> > r.grow in=raster_inverted out=raster_inverted_grown<br>
>> # now both rasters overlap at the border, so you can do whatever you want<br>
>> with it, for instance, extract it<br>
>> > r.mapcalc "border=if(raster_inverted_grown==1 &&<br>
>> > isnull(raster)==0,raster,null())"<br>
>><br>
>> Hope it helps.<br>
>><br>
>> Cheers,<br>
>> Marcello.<br>
>><br>
>><br>
>><br>
><br>
><br>
</div>> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
><br>
</blockquote></div><br>