Description Usage Arguments Note Author(s) References Examples
Determine upslope contributing area based on an elevation raster and, optionally, compute the topographic wetness index.
1 | upslope.area(dem, log = TRUE, atb = FALSE, deg = 0.1, fill.sinks = TRUE)
|
dem |
raster Elevation raster (in m), using a projected coordinate system with identical x and y resolutions. |
log |
Boolean Return the natural log of the values. |
atb |
Boolean If TRUE, include both the upslope contributing area and the topographic wetness index ln(a/tan(beta)). Otherwise calculate just the upslope area. |
deg |
numeric Minimum intercell slope to identify with a sink (degrees). |
fill.sinks |
Fill sinks before calculation using the threshold angle given by deg. |
This is a wrapper to the function implemented in the TOPMODEL package by Wouter Buytaert.
Peter Metcalfe and Wouter Buytaert
Quinn, P. FALSE., Beven, K. J., & Lamb, R. (1995). The In (a/tan/beta) index: How to calculate it and how to use it within the Topmodel framework. Hydrological processes, 9(2), 161-182.
1 2 3 4 5 6 7 8 | ## Not run:
require(dynatopmodel)
data(brompton)
a.atb <- upslope.area(brompton$dem, atb=TRUE)
sp::plot(a.atb, main=c("Upslope area (log(m^2/m))", "TWI log(m^2/m)"))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.