compile_pixel_states: Deriving qualitative pixel states from raster data

In Issoufou-Liman/SpatialProbs: Solving Spatial Problems under Uncertainity

Description

make_pixel_states takes a raster stack or brick to qualify pixel following scale either defined by the user or derived from the natural breaks of data following the following boxplot statistics and outliers.

Usage

compile_pixel_states(x, w = matrix(1, 3, 3), split_IQR = FALSE,
  custum_rclmat = NULL, sample = FALSE, size = 1000, ties = "NA",
  op = c("sampler", "proba"), as_states = FALSE, inparallel = NULL,
  ...)

Arguments

x	a raster stack or brick
w	matrix of weights (the moving window), e.g. a 3 by 3 matrix with values 1; see Details. The matrix does not need to be square, but the sides must be odd numbers. If you need even sides, you can add a column or row with weights of zero
split_IQR	logical. Should the inter-quartile range be split at the median to form different ranges? The default is TRUE
custum_rclmat	numeric vector indicating the limits based on which to reclassify x.
sample	Logical. Should the data be sampled? If TRUE, the result is based on a ramdon sampling of the original data.
size	Integer. Sample size (if sample is TRUE or if the data cannot be process in memory).
ties	character. Indicates how to treat ties. Either 'random', 'lowest', 'highest', 'first', or 'NA'
op	character string indicating the type of returned value: either raw states or their probabilities
as_states	Logical, are the states already made (e.g. x is an output of make_pixel_states function call)? If TRUE, the x is assumed to be a stack of ready made states and is directly processed.
inparallel	integer indicating the number of processes to run in parallel
...	Additional arguments as for writeRaster

Details

The ranges are bounded by the lower outlier (if any), extreme of the lower whisker, the lower ‘hinge’, the median, the upper ‘hinge’, the extreme of the upper whisker, and the upper outlier (if any).

Value

A raster brick where each layers correspond to one of the states.

Author(s)

Issoufou Liman

Examples

sample <- rnorm(100)
get_boxplot_range_1d (sample)

Issoufou-Liman/SpatialProbs documentation built on Oct. 30, 2019, 7:27 p.m.