[GRASS-user] r.contour fails to create contours on large dataset

[Eric Goddard]

Hi all,

 I'm working on a Grass 7 python module to create a local relief model
from lidar derived DTMs. I've tested the module on small and medium
size rasters with less than 200 million cells, and the script returns
results (whether or not the results are correct is another matter).
However, when I try to run the script on an entire county with 2.5
billion cells, all of the steps up to extracting the contours from the
difference between the DTM and the low pass filter work correctly, but
the resulting contours vector is empty, causing the remainder of the
script to fail.

I know I'm short on technical details, but any input into why this
might be happening would be appreciated. The full script is located
here: http://pastebin.com/DudYK0Se

The command used to run r.contours:
grass.run_command("r.contour", input=lp_subtract, output=LP_contour,
minlevel=0, maxlevel=0, step=10)


Thanks,
Eric