[GRASS-user] r.fillnulls failing

GRASS GRASS user
1

I have a rather large map that I need to fill the null values in. I
have compiled grass 6.5. Any help would be greatly appreciated. Is
this the proper way to go about this?

r.fillnulls --overwrite input=ft_bragg_lidar_2009
output=ft_bragg_lidar_2009_no_nullsLocating and isolating NULL
areas...
100%
ERROR: G_malloc: unable to allocate 18446744072071437530 bytes of memory at
       read_map.c:38
ERROR: /usr/local/grass-6.5.svn/scripts/r.fillnulls abandoned. Removing
       temporary map, restoring user mask if needed:
Removing raster <MASK>
WARNING: Raster map <MASK> not found
WARNING: <MASK> nothing removed
Removing raster <r_fillnulls_6692>
Removing raster <r_fillnulls_6692.buf>
WARNING: Raster map <r_fillnulls_6692.buf> not found
WARNING: <r_fillnulls_6692.buf> nothing removed
Removing raster <r_fillnulls_6692_filled>
WARNING: Raster map <r_fillnulls_6692_filled> not found
WARNING: <r_fillnulls_6692_filled> nothing removed
WARNING: raster <usermask_mask.6692> not found

--
Stephen Sefick

Let's not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods. We are mammals, and have not exhausted the
annoying little problems of being mammals.

                -K. Mullis

2

stephen sefick wrote:

ERROR: G_malloc: unable to allocate 18446744072071437530
bytes of memory at
read_map.c:38

It wants 1.84e+10 gigabytes of RAM to do the job.

check that your region settings have not gone bigger than say
40000x40000 rows & columns with 'g.region -p'

probably running this before r.fillnulls helps:

g.region rast=inputmap

Hamish

3

stephen:

I have a rather large map that I need to fill the null values
in.

I believe this is why v.surf.bspline comes with the v.LIDAR
suite of modules, have a look there (+ r.to.vect).

the r.surf.nnbathy addon module might also help.

Hamish

4

projection: 99 (Lambert Conformal Conic)
zone: 0
datum: nad83
ellipsoid: a=6378137 es=0.006694380022900787
north: 570351.02
south: 440360.99
west: 1829042.99
east: 2049032.97
nsres: 3.28075387
ewres: 3.28073939
rows: 39622
cols: 67055
cells: 2656853210

I will have a look at that, thanks.

On Tue, Jan 26, 2010 at 12:18 AM, Hamish <hamish_b@yahoo.com> wrote:

stephen:

I have a rather large map that I need to fill the null values
in.

I believe this is why v.surf.bspline comes with the v.LIDAR
suite of modules, have a look there (+ r.to.vect).

the r.surf.nnbathy addon module might also help.

Hamish

--
Stephen Sefick

Let's not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods. We are mammals, and have not exhausted the
annoying little problems of being mammals.

                -K. Mullis

5

On Tuesday 26 January 2010 02:13:28 stephen sefick wrote:

I have a rather large map that I need to fill the null values in. I
have compiled grass 6.5. Any help would be greatly appreciated. Is
this the proper way to go about this?

r.fillnulls --overwrite input=ft_bragg_lidar_2009
output=ft_bragg_lidar_2009_no_nullsLocating and isolating NULL

v.surf.bspline, v.surf.idw and v.surf.rst are good for LiDAR and proven
through various journal articles. I have found that v.surf.idw has as good
results as these others with default parameters (only 2!) but is 10 times
faster in processing. So if you don't have the time for v.surf.rst, idw may be
the way to go for a quick or large map solution.

John