CIA: Category level intensity analysis

In intensity.analysis: Intensity of Change for Comparing Categorical Maps from Sequential Intervals

Description

Category level intensity analysis

Usage

CIA(crosstabulation, time.points, categories)

Arguments

crosstabulation	List of crosstabulation tables generated by multicrosstab function.
time.points	a charachter vector showing the time point of each raster layer in chronological order.
categories	A charachter vector showing the categories in the map. Order of categories decided bases on the equivalent IDs in the raster attribute table.

Details

Gets the list of crosstabulation tables, time points and categories vectors and returns a list of gain and loss metrics accompanied with relevant bar graphs.

Value

The output is a list of lists. Elements of the list include: gross.loss, gross.gain, loss.intensity, gain.intensity, Uniform.category.intensity, loss.behavior and gain.behavior.

Examples

raster_2005 <- raster::raster(system.file("external/RASTER_2005.RST", package="intensity.analysis"))
raster_2010 <- raster::raster(system.file("external/RASTER_2010.RST", package="intensity.analysis"))
raster_2012 <- raster::raster(system.file("external/RASTER_2012.RST", package="intensity.analysis"))
raster.layers <- list(raster_2005, raster_2010, raster_2012)
time.points <- c("2005","2010","2012")
categories <- c("Water","Trees","Impervious")
crosstabulation <- multicrosstab(raster.layers, time.points, categories)
CIA.output <- CIA(crosstabulation, time.points, categories)

intensity.analysis documentation built on May 2, 2019, 2:51 p.m.