[GRASS5] [bug #3023] (grass) r.fillnulls doesn't work

[Hamish]

> >> this bug's URL: http://intevation.de/rt/webrt?serial_num=3023
> >---------------------------------------------------------------------
> >> > Subject: r.fillnulls doesn't work
> >
..
> >> I'm reffering here to a single 3" SRTM tile loaded using r.in.srtm.
> >>
> >> r.fillnulls works very long and uses so much swap that mouse stops
> >> answering, temporarily, after several minutes.

...?

> > What is your region set to?
> > try 'g.region rast=(map)' first?
> 
> I did it of course:
> 
> GRASS 6.0.cvs:~ > g.region rast=N51E016 -ap
> projection: 3 (Latitude-Longitude)
> zone:       0
> datum:      wgs84
> ellipsoid:  wgs84
> north:      52:00:01.5N
> south:      50:59:58.5N
> west:       15:59:58.5E
> east:       17:00:01.5E
> nsres:      0:00:03
> ewres:      0:00:03
> rows:       1201
> cols:       1201

looks ok.


> > what is the output of "g.region -gm" ??
> 
> GRASS 6.0.cvs:~ > g.region -g
> n=52.00041667
> s=50.99958333
> w=15.99958333
> e=17.00041667
> nsres=0.00083333
> ewres=0.00083333

that looks ok too.  The buffer distance is supposed to be nsres+ewres.
This was fixed to work with floating point resolutions (i.e. lat/lon)
before 6beta2 was released.


> > I've just added some more error checking in CVS, can you update and
> > try again?  maybe uncomment "echo RES=.." line?
> > (It won't fix your problem but will stop the script if this step
> > fails.)
> 
> I can't access CVS at home, my expensive, 2kB/s connection doesn't let
> me to. But I'll grab the next CVS snapshot at my Uni.

It is possible from CVS snapshot to get just one file:

cd grass6/scripts/r.fillnulls
cvs -z3 update -dP

that will just update the files in that directory (and any
subdirectories), very little bandwidth.

or from the CVS web interface you can grab the files:

http://freegis.org/cgi-bin/viewcvs.cgi/*checkout*/grass6/scripts/r.in.srtm/r.in.srtm?rev=HEAD&content-type=text/plain

http://freegis.org/cgi-bin/viewcvs.cgi/*checkout*/grass6/scripts/r.fillnulls/r.fillnulls?rev=HEAD&content-type=text/plain

then run "make" from within the modules' directory. (+ install?)

both scripts have new bugfixes and extra checks. can you test them?


> > - In the script we should temporarily zoom to the raster limits with
> > g.region first?
> 
> Not really a good idea I guess. Anytime it should be up to the user
> what region he wants to patch.

ok.


> On thing more: my bug report was based r.fillnuls in GUI mode. After
> running  it from the command line I see there are two more errors. Let
> me quote them  all now:

.. pretty much everything will be broken after the r.buffer step fails.
Please try again with the two updated scripts.



Hamish