Description Usage Arguments Value References See Also Examples
Cross tabulate land use transitions using
raster::crosstab
. This step should form the basis of
further research into the processes driving the most important transitions in
the study region (Pontius et al., 2004).
1 2 3 4 5 6 7 8 | crossTabulate(x, y, ...)
## S4 method for signature 'RasterLayer,RasterLayer'
crossTabulate(x, y, categories,
labels = as.character(categories), ...)
## S4 method for signature 'ObsLulcRasterStack,ANY'
crossTabulate(x, y, times, ...)
|
x |
RasterLayer representing land use map from an earlier timestep or an ObsLulcRasterStack object containing at least two land use maps for different points in time |
y |
RasterLayer representing land use map from a later timestep. Not used
if |
... |
additional arguments to |
categories |
numeric vector containing land use categories to consider.
Not used if |
labels |
character vector (optional) with labels corresponding to
|
times |
numeric vector representing the time points of two land use maps from ObsLulcRasterStack |
A data.frame.
Pontius Jr, R.G., Shusas, E., McEachern, M. (2004). Detecting important categorical land changes while accounting for persistence. Agriculture, Ecosystems & Environment 101(2):251-268.
ObsLulcRasterStack
, raster::crosstab
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ## Plum Island Ecosystems
## Load observed land use maps
obs <- ObsLulcRasterStack(x=pie,
pattern="lu",
categories=c(1,2,3),
labels=c("forest","built","other"),
t=c(0,6,14))
crossTabulate(x=obs, times=c(0,14))
## RasterLayer input
crossTabulate(x=obs[[1]],
y=obs[[3]],
categories=c(1,2,3),
labels=c("forest","built","other"))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.