schoenersD: Calculates Schoener's D-value between two RasterLayers.

In cjcampbell/IsoModAT: Isotope Origin Clustering and Assignment Tools

Calculates Schoener's D-value between two RasterLayers.

Description

RasterLayers must have identical resolutions and extents. The function will automatically normalize surfaces to sum to 1.

Usage

schoenersD(rast1, rast2)

Arguments

rast1	Input RasterLayer
rast2	Input RasterLayer 2

Details

Calculates similarity value of two RasterLayers using Schoener's D-metric.

Examples

# Generate example probability surfaces.
myiso <- raster::rasterFromXYZ(isoscape)
myiso_sd <- raster::rasterFromXYZ(isoscape_sd)
df <- data.frame(
        ID = c(-100, -80, -50),
        isotopeValue = c(-100, -80, -50),
        SD_indv = rep(5, 3)
        )
assignmentModels <- isotopeAssignmentModel(
        ID = df$ID,
        isotopeValue = df$isotopeValue,
        SD_indv = df$SD_indv,
        precip_raster = myiso,
        precip_SD_raster = myiso_sd,
        nClusters = FALSE
        )

# Calculate Schoener's D-metric of spatial similarity between two of the
# example probability surfaces.
schoenersD(assignmentModels[[1]], assignmentModels[[2]])
## 0.969156

cjcampbell/IsoModAT documentation built on Feb. 21, 2024, 10:27 a.m.