|
The basis for the derivation of many hydrologic parameters is a drainage network, which can manually be digitised into a grid or extracted automatically. The latter is usually the case. You will learn how to do this in the following sections. Local drain directions (i.e. flow directions) in grids are often coded according to the numeric pad on a computer keyboard. A cell discharging southwest would therefore be allocated the value 1, a cell discharging west would be allocated the value 4 and so on.
A map assigning a local drain direction to each grid point is called a local drain direction net or ldd net (Burrough et al. 1998). Various algorithms exist for computing ldd nets, each based on different assumptions on the drainage of water in the terrain.
The D8 algorithm is based on the quite strict assumption, that within a window of 9 cells water
always flows in the direction of the steepest slope. I.e. the slope to each of the 8 neighbours of a
central source cell has to be computed using the following formula:
where d is the distance and
is the elevation difference between two grid points.
In a raster, where the cell width is 1 we have to distinguish between two cases:
Let’s have a look at the following example, where the source cell and its 8 neighbours have defined elevation values and the cell size is 1, 2 or 3. Where does the water drain to? Select the correct cell by clicking on it. (Attention: this little interactive demonstration is not perfect. Where two cells might serve as target cells equally, there is no special solution implemented - the application simply chooses the first of these two cells in a sequence!)
We now know how to calculate the local drain direction for a single cell in the center of a window of 9 cells. By moving the window stepwise through the entire map we can assign a local drain direction to every cell and get a "local drain direction net" or "ldd net" at the end.
The example illustration of the flow directions below shows a 25m raster DHM from the Tuerlersee area and the flow direction map derived from it. The area of the lake in the southern part of the map had to be excluded because the drain direction on a horizontal plain is undefined.
An ldd net can alternatively be visualised with a vector map as shown in the following image.