[GRASS-user] r.mapcalc median filters entire map instead of only null values

[Hamish]

Patton, Eric wrote:
> I'm trying to use r.neighbors and r.mapcalc in a script to fill holes
> in a raster where there are null values. A temporary raster is created
> with r.neighbors that contains a 3x3 median-interpolated surface of
> the primary DEM. r.mapcalc should only use the value in this temporary
> raster where a null value occurs in the primary DEM.


perhaps modify/extend the r.fillnulls script? A flag to use an alternate
interpolation technique might be nice.


Hamish