Dear Mxxxx,
I programmed some of my recent zonalstats/DOCELL/IFTHENELSE (AG slang)
work in numpy/scipys functionality for terrain analysis.
Have a look there, and in combination with gdal its rather
straightforward. The only problem you might run into is that at one
stage you need to implement tiling as your memory might not be able to
fit your dem into it.
Cheers hannes
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev