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

[Patton, Eric]

It's ok, Glynn tested the script out and it worked for him, so I can't
explain what was happening. My script even works for me now. User error,
I guess ;)

I didn't know r.fillnulls was a script, as I've never used it before.
Maybe I could take a look at using other interpolation methods as you
suggest.

~ Eric.

--
Eric Patton 
email: epatton at nrcan.gc.ca


 

> -----Original Message-----
> From: Hamish [mailto:hamish_nospam at yahoo.com] 
> Sent: Tuesday, February 06, 2007 4:58 AM
> To: Patton, Eric
> Cc: grassuser at grass.itc.it
> Subject: Re: [GRASS-user] r.mapcalc median filters entire map 
> instead of only null values
> 
> 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
>