quantileAtSamplingLocation: Quantile at coordinates function

In cjcampbell/IsoModAT: Isotope Origin Clustering and Assignment Tools

Quantile at coordinates function

Description

Function estimates percentile of each non-NA value within a RasterLayer using the empirical cumulative distribution function, and extracts value at location specified. For more information, see help(ecdf).

Usage

quantileAtSamplingLocation(indivraster, Lat, Lon)

Arguments

indivraster	A RasterLayer representing normalized probability of origin surface
Lat	Integer latitude
Lon	Integer longitude

See Also

makeQuantileSurfaces

Examples

# Generate example probability surface.
myiso <- raster::rasterFromXYZ(isoscape)
myiso_sd <- rasterFromXYZ(isoscape_sd)
exampleSurface <- isotopeAssignmentModel(
         ID = "A",
         isotopeValue = -100,
         SD_indv = 5,
         precip_raster = myiso,
         precip_SD_raster = myiso_sd,
         nClusters = FALSE
         )
# Calculate quantile probability value at specific point.
set.seed(1)
x <- sample( which( !is.na(exampleSurface[]) ), size = 1)
pt <- raster::xyFromCell(exampleSurface, x)
quantileAtSamplingLocation(exampleSurface, Lat = pt[2], Lon = pt[1])

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