Showing posts with label arcgis. Show all posts
Showing posts with label arcgis. Show all posts

October 09, 2012

digging out the edge of the hill

Recently I had a possibility and time to look practically at the boundary questions, the ones I always like to philosophy about...  I got a simple task: to calculate the volume of a small heap (much smaller then hill) having an image with elevation values (DEM) and so I just needed to automatize and to simplify the steps that were already determined by others, mostly some routine work with ArcGIS (ESRI) tools, so just a bit of patience learning to geoprocess with python.

Digging the boundary  of a heap
But before starting the calculations I had to find out where is this heap and where it's not anymore. I had to draw clear boundaries in order to start some calculations at all. Apparently there was no tool for a heap determination founded where I work, till now all was done manually, by digitizing: zooming out and in, searching out that boundary where the heap starts or ends and going around.
The human brain is great, fast recognizing the pattern it searches for. Maybe recognizing just approximately, but nicely and unconsciously ignoring all other: small depressions, pikes or elevation model noise due to calculated heights of some bushes or some artificial constructions  or photogrametrical noise...
But the problem is that when having a huge images with a great resolution it takes hours and hours to load it and while zooming back and forward trying to draw a line in order to determine "exactly" where the heap starts or/and ends. Besides, when you start looking at a greater and greater resolution the picture of  all the image disappears and for our brain it becomes hard to process just various tones of colours in pixels in front. So, I needed to find how to automatize this first step in order I could work on the calculations. And just looking rather at a simple task or indeed at the DEM (it's prety nice!),  I fastly ran into a dark dark forest searching the way how to determine where this heap finishes or starts in order to be able to define it and build the script.