Layerize — layerize • raster

Create a RasterBrick with a Boolean layer for each class (value, or subset of the values) in a RasterLayer. For example, if the cell values of a RasterLayer indicate what vegetation type they are, this function will create a layer (presence/absence; dummy variable) for each of these classes. Classes and cell values are always truncated to integers.

You can supply a second spatially overlapping RasterLayer with larger cells (do not use smaller cells!). In this case the cell values are counts for each class. A similar result might be obtained more efficiently by using layerize with a single RasterLayer followed by aggregate(x, , sum).

Usage

# S4 method for class 'RasterLayer,missing'
layerize(x, classes=NULL, falseNA=FALSE, filename='', ...)

# S4 method for class 'RasterLayer,RasterLayer'
layerize(x, y, classes=NULL, filename='', ...)

Arguments

x: RasterLayer
y: RasterLayer or missing
classes: numeric. The values (classes) for which layers should be made. If NULL all classes are used
falseNA: logical. If TRUE, cells that are not of the class represented by a layer are NA rather then FALSE
filename: character. Output filename (optional)
...: Additional arguments as for writeRaster

Value

RasterBrick

Examples

r <- raster(nrow=20, ncol=20)
values(r) <- c(rep(NA, 50), rep(1:5, 70))
b <- layerize(r)

r2 <- raster(nrow=5, ncol=5)
b2 <- layerize(r, r2)