| envChange {changeRangeR} | R Documentation |
Calculate change in suitable SDM area through time
Description
Calculate SDM area after masking for environmental variables through time
Usage
envChange(rStack, binaryRange, threshold, bound, correlation = FALSE)
Arguments
rStack |
rasterStack of environmental variable to measure within binary SDM through time |
binaryRange |
raster object or shapefile of binary range (SDM, AOO, EOO) with same projection as rStack |
threshold |
integer (or integers if bound = "both") of where rStack layers should be thresholded |
bound |
character string characterizing the way the threshold should happen. "upper" removes values above the threshold (e.g., maximum human footprint)."lower" removes values below the threshold (e.g., minimum forest cover). "neither" does not threshold at any point. "both" thresholds at both threshold values (if provided; e.g., minimum and maximum temperature). |
correlation |
boolean. If FALSE, environmental variable will be converted to a binary map and used as a mask. If TRUE, environmental variable is only thresholded by bounds, but left continuous. Then, Pearson's correlation coefficient with SDM will be computed for overlapping areas. |
Value
A list two objects. Area is a data.frame showing the total area of each masked raster. masks is a rasterStack of each masked raster.
Author(s)
pgalante@amnh.org
Examples
# create rStack
r1 <- raster::raster(nrows=108, ncols=108, xmn=-50, xmx=50)
raster::values(r1)<- runif(n = (108*108))
r2 <- raster::raster(nrows=108, ncols=108, xmn=-50, xmx=50)
raster::values(r2)<- runif(n = (108*108))
rStack <- raster::stack(r1,r2)
# create binaryRange
binaryRange <- raster::crop(r1, raster::extent(c(-50, 50, 0, 90)))
binaryRange <- raster::extend(binaryRange, r1)
binaryRange[!is.na(binaryRange)] <- 1
# set threshold
threshold <- 0.5
# set bound
bound <- "upper"
# Run function
envChange(rStack = rStack, binaryRange = binaryRange, threshold = threshold, bound = bound)