Learn R Programming

intensity.analysis (version 0.1.6)

CIA: Category level intensity analysis

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.

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.

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.

Examples

Run this code
# NOT RUN {
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)
# }

Run the code above in your browser using DataLab