Update a DisturbanceMetrics object in order to avoid reprocessing parts of the data that have not changed. New data is supplied as a named list and the object is updated depending on the elements provided in the list.
updateDisturbance(distMet, newData, linBuffMethod = "raster", isPercent = TRUE)DisturbanceMetrics object
named list of objects to be used to update distMet. Potential names are: natDist, anthroDist, and linFeat.
character. The method used to buffer linear features if they are supplied as sf lines. The default is "raster" in which case they are rasterized then buffered using a moving window method. If "sf" then the lines are buffered with st_buffer and then rasterized. Either way points are included in the raster output.
logical. Should the results be returned as a percentage?
A DisturbanceMetrics Object see DisturbanceMetrics-class()
If newData contains only linFeat or anthroDist then the new data will
be combined with the existing data in the DisturbanceMetrics object to
determine the buffered anthropogenic disturbance
If newData contains only natDist then only Fire is updated.
DisturbanceMetrics-class() for information on the object
returned and disturbanceMetrics() for making a new
DisturbanceMetrics object.
Functions for calculating disturbance:
DisturbanceMetrics-class,
disturbanceMetrics(),
reclassDist(),
results()
# create example rasters
lc <- terra::rast(xmin = 0, xmax = 10, ymin = 0, ymax = 10,
ncols = 10, nrow = 10, crs = "EPSG:5070")
nd <- lc
nd[1:3, 1:3] <- 1
ad <- lc
ad[3:5, 3:5] <- 1
lc[] <- 1
# create sf objects
lf <- sf::st_as_sf(sf::st_sfc(list(sf::st_linestring(matrix(c(0, 0, 10, 10),
ncol = 2, byrow = TRUE))),
crs = 5070))
projPol <- sf::st_sf(sf::st_as_sfc(sf::st_bbox(ad)))
# calculate disturbance
dm <- disturbanceMetrics(landCover = lc,
linFeat = lf,
natDist = nd,
anthroDist = ad,
projectPoly = projPol,
padFocal = TRUE,
bufferWidth = 1)
#> buffering anthropogenic disturbance
#> calculating disturbance metrics
# new linear features
lf2 <- sf::st_as_sf(sf::st_sfc(list(sf::st_linestring(matrix(c(0, 0, 10, 10),
ncol = 2, byrow = TRUE)),
sf::st_linestring(matrix(c(0, 9, 9, 10),
ncol = 2, byrow = TRUE)),
sf::st_linestring(matrix(c(10, 0, 5, 5),
ncol = 2, byrow = TRUE))),
crs = 5070))
dm2 <- updateDisturbance(dm, newData = list(linFeat = lf2))
#> buffering anthropogenic disturbance
#> calculating disturbance metrics