Fill sinks in a digital elevation model

Share:

Description

Removes sinks in a digital elevation model by filling depressions

Usage

1
  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 http://rwiki.sciviews.org/doku.php?id=guides:tutorials:hydrological_data_analysis for examples.

See Also

topmodel, topidx

Examples

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