Removes sinks in a digital elevation model by filling depressions

1 | ```
sinkfill(DEM,res,degree)
``` |

`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] |

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.

A matrix of the same size as the DEM.

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

See http://rwiki.sciviews.org/doku.php?id=guides:tutorials:hydrological_data_analysis for examples.

topmodel, topidx

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

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.