multicrosstab: Generate crosstabulation table(s)

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

Description

Generate crosstabulation table(s)

Usage

multicrosstab(raster.layers, time.points, categories)

Arguments

raster.layers	List of raster objects.
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

Read categorical raster objects, checks their dimentionality, resolution and extent are a perfect match. Then the function returns a list consist of a crosstabulation table for each consequtive pair of raster layers.

Value

The output is a list of crosstabulation table(s). The length of the list is equal to the number of intervals.

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)

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