[GRASS-dev] r.fillnulls failure (grass7 )

epy@epinux:~/dev/grass_trunk$ svn up
Updating '.':
At revision 55645.


GRASS 7.0.svn (utm19N_wgs84):~ > g.region rast=dtm1 -ap
projection: 1 (UTM)
zone:       19
datum:      wgs84
ellipsoid:  wgs84
north:      4720170
south:      4715050
west:       391340
east:       398580
nsres:      10
ewres:      10
rows:       512
cols:       724
cells:      370688

GRASS 7.0.svn (utm19N_wgs84):~ > r.fillnulls input=dtm1 output=dtm1 --overwrite

Using RST interpolation...
Locating and isolating NULL areas...
 100%
Growing NULL areas
Assigning IDs to NULL areas
 100%
Processing 1 map holes
Filling hole 2 of 1
 100%
 100%
 100%
Note: The following warnings about number of points for interpolation may
be ignored.
ERROR: Database connection not defined for layer 1
ERROR: Failed to fill hole 2
WARNING: Raster map <r_fillnulls_5883_hole_2_dem> not found
WARNING: <r_fillnulls_5883_hole_2_dem> nothing removed

GRASS 7.0.svn (utm19N_wgs84):~ > 

On Sat, Apr 6, 2013 at 4:25 PM, epi <massimodisasha@gmail.com> wrote:
...

GRASS 7.0.svn (utm19N_wgs84):~ > r.fillnulls input=dtm1 output=dtm1
--overwrite

You try to overwrite the source with the target, to my knowledge invalid
(sure, the module should say so).

Try
r.fillnulls input=dtm1 output=dtm1_filled

Markus

I used the --overwrite option

i had problems with r.fillnulls recently see :
http://osgeo-org.1560.n6.nabble.com/r-fillnulls-error-lt-z-gt-is-not-a-valid-flag-grass7-td5034768.html
after that fix the same syntax was working

changing the output to a new name i got the same error

GRASS 7.0.svn (utm19N_wgs84):/home/epy > r.fillnulls input=dtm1 output=dtm1fill --overwrite
Using RST interpolation...
Locating and isolating NULL areas...
100%
Growing NULL areas
Assigning IDs to NULL areas
100%
Processing 1 map holes
Filling hole 2 of 1
100%
100%
100%
Note: The following warnings about number of points for interpolation may
be ignored.
ERROR: Database connection not defined for layer 1
ERROR: Failed to fill hole 2
WARNING: Raster map <r_fillnulls_4517_hole_2_dem> not found
WARNING: <r_fillnulls_4517_hole_2_dem> nothing removed
GRASS 7.0.svn (utm19N_wgs84):/home/epy >

Il giorno 06/apr/2013, alle ore 11:08, Markus Neteler <neteler@osgeo.org> ha scritto:

On Sat, Apr 6, 2013 at 4:25 PM, epi <massimodisasha@gmail.com> wrote:
...

GRASS 7.0.svn (utm19N_wgs84):~ > r.fillnulls input=dtm1 output=dtm1
--overwrite

You try to overwrite the source with the target, to my knowledge invalid
(sure, the module should say so).

Try
r.fillnulls input=dtm1 output=dtm1_filled

Markus

On Sat, Apr 6, 2013 at 5:15 PM, epi <massimodisasha@gmail.com> wrote:
...

GRASS 7.0.svn (utm19N_wgs84):/home/epy > r.fillnulls input=dtm1 output=dtm1fill --overwrite
Using RST interpolation...
Locating and isolating NULL areas...
100%
Growing NULL areas
Assigning IDs to NULL areas
100%
Processing 1 map holes
Filling hole 2 of 1

^^^^^^^^^^^^^^^^^^^^^^^

... this looks odd

100%
100%
100%
Note: The following warnings about number of points for interpolation may
be ignored.
ERROR: Database connection not defined for layer 1

^^^- same problem here (I tried the example in the manual page
which fails - a bug (or two).

Markus

On Sun, Apr 7, 2013 at 10:47 PM, Markus Neteler <neteler@osgeo.org> wrote:

On Sat, Apr 6, 2013 at 5:15 PM, epi <massimodisasha@gmail.com> wrote:
...

GRASS 7.0.svn (utm19N_wgs84):/home/epy > r.fillnulls input=dtm1 output=dtm1fill --overwrite
Using RST interpolation...
Locating and isolating NULL areas...
100%
Growing NULL areas
Assigning IDs to NULL areas
100%
Processing 1 map holes
Filling hole 2 of 1

^^^^^^^^^^^^^^^^^^^^^^^

... this looks odd

100%
100%
100%
Note: The following warnings about number of points for interpolation may
be ignored.
ERROR: Database connection not defined for layer 1

^^^- same problem here (I tried the example in the manual page
which fails - a bug (or two).

Try to remove zcol = 'value' from the call to v.surf.rst in the
r.fillnulls script. The vector created by r.fillnulls with r.to.vect
is 3D, so no z column is needed (or present).

Markus M

On Sun, Apr 7, 2013 at 11:00 PM, Markus Metz
<markus.metz.giswork@gmail.com> wrote:
...

Try to remove zcol = 'value' from the call to v.surf.rst in the
r.fillnulls script. The vector created by r.fillnulls with r.to.vect
is 3D, so no z column is needed (or present).

Tested, this solves the problem. I guess it is an artifact of the -z usage
in the GRASS 6 version or the r.to.vect usage.

markusN

thanks this fixed my problem! :slight_smile:

Il giorno 07/apr/2013, alle ore 17:20, Markus Neteler <neteler@osgeo.org> ha scritto:

On Sun, Apr 7, 2013 at 11:00 PM, Markus Metz
<markus.metz.giswork@gmail.com> wrote:
...

Try to remove zcol = 'value' from the call to v.surf.rst in the
r.fillnulls script. The vector created by r.fillnulls with r.to.vect
is 3D, so no z column is needed (or present).

Tested, this solves the problem. I guess it is an artifact of the -z usage
in the GRASS 6 version or the r.to.vect usage.

markusN