r <- terra::rast(system.file("extdata/SICILY_1000m.tif", package="leastcostpath"))
slope_cs <- create_slope_cs(x = r, cost_function = "tobler", neighbours = 4)
locs <- sf::st_sf(geometry = sf::st_sfc(
sf::st_point(c(960745, 4166836)),
crs = terra::crs(r)))
locs <- sf::st_buffer(x = locs, dist = 25000)
slope_cs2 <- update_values(x = slope_cs, sf = locs,
FUN = function(j) { j * 0.6})
slope_cs3 <- update_values(x = slope_cs, sf = locs,
FUN = function(j) { j + 10})
slope_cs4 <- update_values(x = slope_cs, sf = locs,
FUN = function(j) { replace(x = j, values = 0)})
Run the code above in your browser using DataLab