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 'DiscreteLulcRasterStack,ANY'
crossTabulate(x, y, times, ...)
|
x |
RasterLayer representing land use map from an earlier timestep or an LulcRasterStack 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 LulcRasterStack |
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.
LulcRasterStack
, raster::crosstab
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | ## Not run:
## Plum Island Ecosystems
## load observed land use maps
lu <- DiscreteLulcRasterStack(x=stack(pie[1:3]),
categories=c(1,2,3),
labels=c("Forest","Built","Other"),
t=c(0,6,14))
crossTabulate(x=lu, times=c(0,6))
crossTabulate(x=lu, times=c(0,14))
## RasterLayer input
crossTabulate(x=lu[[1]],
y=lu[[3]],
categories=c(1,2,3),
labels=c("forest","built","other"))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.