First post! I'm new to python and I'm getting into it because I'm constantly facing novel problems in gis, and scripts seem to be the only way to get around them.
I'm tasked with determining the total area of land affected by road salt in Northern New York. More or less 6 million acres. Going small area by small area is easy, but it would take months.
What I need my script to do:
The program must go cell by cell on a dem array. Once it hits a cell flagged as a ROAD a 8 neighbor window will be created. All cells less than the elevation of the road in the 8 neighbor window will be exported. This will run until the end of the array. The script will loop back and go cell by cell until it hits one of these "exported" cells. Once it hits these exported cells another 8 neighbor window will be created. All cells less than the elevation of the targeted (previously exported) cell will exported. This will loop and loop to completion. This process will be completed once there aren't any cells less than the targeted cells. When it hits a stream it will stop at the stream because the only cells less than the targeted cells will be downstream... Hopefully producing a raster representing all land affected by road salt (land lower in elevation, yet continuous on a down slope.)
So I'm interested in what particular module will help with this. And where should I start looking for this cell by cell encoding stuff. I know fortran would be perfect. But I need esri compatibility. Hence python