[GRASS-user] execGRASS("r.diversity") Done. No rasters created (Large rasters)

[Markus Neteler]

Dear Sandra,

On Sat, Jun 14, 2014 at 9:21 AM, Sandra MacFadyen <sandramf at live.co.za> wrote:
> Dear GRASS Users,
>
> I am using r.diversity (GRASS GIS 7.0.0svn build 60785 win32) through R (R version 3.0.2 win32) on Windows 7 64bit.
>
> Testing from GRASS and North_Carolina dataset's lsat5_1987_10 raster (cells=29535), [r.diversity --o input=lsat5_1987_10 at landsat prefix=lsat5_1987_div_ alpha=0.5 size=5] reports as Done with the different diversity rasters added to the location fine (see # North_Carolina# code and results below).

Yes, same here on Linux. So far so nice.

BTW: I have changed r.diversity in SVN and removed the quiet flags to
print the output map names rather than only "Done".

> However, when running the same code on a larger image (cells=6746328) from my own location, although it reports Done, no rasters are created. If I subset the image (cells=1632830) and run it again its works (see # sub2Kruger # code and results below). So I'm guessing it is a memory issue?

ok I simply simulate higher resolution:

g.region -p rast=lsat5_1987_10 res=5 -a
...
cells:      8137540

and then execute the r.diversity call again. This then calls

r.li.simpson --o input=lsat5_1987_10 at landsat conf=conf_diversity_5.0
output=lsat5_1987_div__simpson_size_5.0

and it consumes a lot of memory... will check on a bigger machine,
perhaps a memory leak.

Markus