Removes sinks in a digital elevation model by filling depressions
Usage
sinkfill(DEM,res,degree)
Value
A matrix of the same size as the DEM.
Arguments
DEM
A matrix representing a digital elevation model [m] with equally sized pixels and equal x and y resolution
res
Resolution of the digital elevation model [m]
degree
Minimum slope to be kept between cells when filling [degrees]
Author
Wouter Buytaert, Imperial College London, based on an implementation from the Hydrology Group of Lancaster University
Details
For deep sinks or large maps, it may be possible that not all sinks are filled in one run. Then the function should be applied repeatedly over the same object.