layerize: Layerize

In raster: Geographic Data Analysis and Modeling

Description

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 signature 'RasterLayer,missing'
layerize(x, classes=NULL, falseNA=FALSE, filename='', ...)

## S4 method for signature '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)

Example output

Loading required package: sp

raster documentation built on Jan. 5, 2021, 3:01 a.m.