Skip to contents

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")

Value

SpatRaster

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)
}