LS_alpine
LS_alpine.Rmd
library(Rsagacmd)
library(terra)
library(rusleR)
# 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)