sinkfill {topmodel}

Fill sinks in a digital elevation model

Description

Removes sinks in a digital elevation model by filling depressions

Usage

  sinkfill(DEM,res,degree) 

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]

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.

Value

A matrix of the same size as the DEM.

Author(s)

Wouter Buytaert, Imperial College London, based on an implementation from the Hydrology Group of Lancaster University

References

See https://github.com/ICHydro/topmodel for examples.

See Also

topmodel, topidx

Examples

data(huagrahuma.dem)
filled.dem <- sinkfill(huagrahuma.dem, 25, 0.1)