Digital elevation model (DEM) data in a raster format can be used to automatically derive the drainage characteristics of an area. A procedure has been designed that is capable of operating on matrices of elevation data having no algorithmically imposed size limit, while performing within the resolution and accuracy tolerances of the DEM data.

Each cell is processed as the center of a 3- by 3-cell spatial window in the raster elevation data. If a cell is a local minimum in comparison with two of its non-adjacent neighbors, it is labeled as a drainage cell. The linkages of the drainage cells within user-specified distance and elevation thresholds are established in a separate process. The products of these processing steps are digital masks of the drainage cells and the watershed basins, both in raster format.

A drainage cell mask derived using this procedure is useful in computing slope values for a raster data base. Slope has traditionally been calculated for each cell by fitting a plane through the eight nearest cells. However, if the terrain represented by these cells is V-shaped, such as a gully, a plane does not fit well; in fact, the desired slope value is the slope along the bottom of the gully, regardless of the steepness of the gully sides. The automated drainage process will label such a cell as a drainage cell, and its slope can then be computed from the elevation values of neighboring drainage cells.