projectSummaryMaxSurface: Project probability-of-origin surfaces into one summary...

In isocat: Isotope Origin Clustering and Assignment Tools

Project probability-of-origin surfaces into one summary surface.

Description

Create a summary surface showing which RasterLayer in a Stack has the highest value at a given location. For each cell in a RasterStack, this function returns the identity of the RasterLayer with the highest value at that cell. This surface is intended as a visual summary of common origins, not a basis for quantitative analysis.

Usage

projectSummaryMaxSurface(surfaces, nClust = FALSE)

Arguments

surfaces	SpatRaster where each layer represents a probability-of-origin surface
nClust	Depreciated. Formerly enabled multi-core processing; retained for back-compatibility. A non-'FALSE' value now issues a message and proceeds serially.

Examples

# Create and cluster example assignment surfaces.
myiso <- rast(isoscape, type="xyz")
myiso_sd <- rast(isoscape_sd, type="xyz")
assignmentModels <- isotopeAssignmentModel(
        ID = LETTERS[1:4],
        isotopeValue = seq(-120,-25,length.out = 4),
        SD_indv = rep(5,4),
        precip_raster = myiso,
        precip_SD_raster = myiso_sd,
        nClusters = FALSE
        )
plot(assignmentModels)

# Project mean aggregate surfaces into space.
summaryMap <- projectSummaryMaxSurface(
        surfaces = assignmentModels,
        nClust = FALSE
        )
plot(summaryMap)

isocat documentation built on July 9, 2026, 5:10 p.m.