This function calculates an LS-factor from Universal Soil Loss Equation (USLE), as described in Schmidt et al. (2019).
Usage
ls_alpine(dem, threshold = 120, saga_obj = saga)
Arguments
- dem
SpatRaster, single layer with elevation values. Values should have the same unit as the map units, or in meters when the crs is longitude/latitude
- threshold
numeric/double, specifying a flow threshold in meters. Default is 120 m
- saga_obj
SAGA-GIS geoprocessor object. Run Rsagacmd::saga_gis(raster_backend = "terra")
References
Schmidt, Simon, Simon Tresch, and Katrin Meusburger. “Modification of the RUSLE Slope Length and Steepness Factor (LS-Factor) Based on Rainfall Experiments at Steep Alpine Grasslands.” MethodsX 6 (2019): 219–29. https://doi.org/10.1016/j.mex.2019.01.004.
Examples
if (FALSE) {
library(Rsagacmd)
# initiate a saga object
saga <- saga_gis(raster_backend = "terra")
# load DEM
f <- system.file("extdata/dem.tif", package="rusleR")
DEM <- rast(f)
# calculate LS-alpine
ls <- ls_alpine(dem = DEM)
plot(ls)
}