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 _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
