Flow

There are phenomena  that do not spread in all directions, but move or “flows” along a given, least-cost path, determined by characteristics of local terrain. The typical case arises when we want to determine drainage patterns in a catchment area: rain water “chooses” a way to leave the area. 

We can illustrate the principles involved in this typical case with the simple elevation raster provided in Figure a. For each cell in that raster, the steepest downward slope to a neighbour cell is computed and its direction is stored in a new raster (Figure b). This computation determines the elevation difference between the cell and the neighbour cell and it takes into account cell distance - 1 for neighbour cells in N–S or W–E direction, 2 for cells in a NE–SW or NW–SE direction. From among its eight neighbour cells, it picks the one with the steepest path to it. The directions thus obtained in raster (b) are encoded in integer values; we have “decoded” them for the sake of illustration. Raster (b) can be called the flow-direction raster . From raster (b), the GIS can compute the accumulated flow-count raster, a raster that for each cell indicates how many cells have their water flow into that cell.

Cells with a high accumulated flow count represent areas of concentrated flow and may, thus, belong to a stream. By using some appropriately chosen threshold value in a map algebra expression, we may decide whether they do or not. Cells with an accumulated flow count of zero are local topographic highs and can be used to identify ridges.

While computing flow, the GIS takes proper care of the path lengths. Obviously the flow from a cell C to its neighbour cell to the east C_e is shorter than to the cell that is its northeastern neighbour C_ne . The distance ratio between these two cases is 1/√2.