From earlier post:
“I’ve managed to pilot most of a fast high resolution landscape position workflow with PovRay as my magic tool. The final steps I hope to pipe through PostGIS Raster. In the meantime a screenshot and description: blues are riparian, raw ocre, etc upland categories, grey is mostly flat lake plain and mid slopes, all derived from just a high res DEM input (no hydro lines as input, so it works on areas where you don’t know where the streams may be). There will be more categories in the final, in the meantime, welcome to December.”
I didn’t use PostGIS Raster, but the incredibly useful gdal_calc.py to finish out my landscape position calculations. I will endeavor to post my code to github soon, but the parts and pieces include using gdal and PovRay. PovRay helps with the sampling (really!) of nearby neighbors in the raster DEM, and gdal does the averaging and differencing of those to get relative landscape position. I spent some time yesterday scratching my head over how to show all the new landscape position information on a readable and useable map, and after discussion with collegues, decided to use it to divide the world into two categories– riparian & lowland + headwaters & upland (not reflected yet in the labels). To find out more about landscape position, follow this link, or better yet this one. (BTW, the green is park land, blue is riparian/lowland/stream networks, purple is the basin boundary).