Cross-tabulate two RasterLayer objects, or mulitiple layers in a RasterStack or RasterBrick to create a contingency table.
# S4 method for Raster,Raster
crosstab(x, y, digits=0, long=FALSE, useNA=FALSE, progress='', ...)# S4 method for RasterStackBrick,missing
crosstab(x, digits=0, long=FALSE, useNA=FALSE, progress='', ...)
A table or data.frame
Raster* object
Raster* object if x is a RasterLayer; Can be missing if x is a RasterStack or RasterBrick
integer. The number of digits for rounding the values before cross-tabulation
logical. If TRUE the results are returned in 'long' format data.frame instead of a table
logical, indicting if the table should includes counts of NA values
character. "text", "window", or "" (the default, no progress bar), only for large files that cannot be processed in one step
additional arguments. none implemented
r <- raster(nc=5, nr=5)
values(r) <- runif(ncell(r)) * 2
s <- setValues(r, runif(ncell(r)) * 3)
crosstab(r,s)
rs <- r/s
r[1:5] <- NA
s[20:25] <- NA
x <- stack(r, s, rs)
crosstab(x, useNA=TRUE, long=TRUE)
Run the code above in your browser using DataLab