crosstab: Cross-tabulate
In raster: Geographic Data Analysis and Modeling

Description Usage Arguments Value See Also Examples

Cross-tabulate two RasterLayer objects, or mulitiple layers in a RasterStack or RasterBrick to create a contingency table.

## S4 method for signature 'Raster,Raster'
crosstab(x, y, digits=0, long=FALSE, useNA=FALSE, progress='', ...)

## S4 method for signature 'RasterStackBrick,missing'
crosstab(x, digits=0, long=FALSE, useNA=FALSE, progress='', ...)

`x`	Raster* object
`y`	Raster* object if `x` is a RasterLayer; Can be missing if `x` is a RasterStack or RasterBrick
`digits`	integer. The number of digits for rounding the values before cross-tabulation
`long`	logical. If `TRUE` the results are returned in 'long' format data.frame instead of a table
`useNA`	logical, indicting if the table should includes counts of `NA` values
`progress`	character. "text", "window", or "" (the default, no progress bar), only for large files that cannot be processed in one step
`...`	additional arguments. none implemented

A table or data.frame

freq, zonal

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)

Loading required package: sp
       layer.2
layer.1 0 1 2 3
      0 2 1 0 0
      1 1 6 6 1
      2 1 2 4 1
   layer.1 layer.2 layer.3 Freq
1        0       0       0    1
2        0       0       2    1
3        1       0       8    1
4        1       1       0    1
5        1       1       1    2
6        1       2       0    2
7        1       2       1    1
8        1       3       0    1
9        1      NA       0    1
10       1      NA       1    3
11       2       0     472    1
12       2       1       2    1
13       2       2       1    2
14       2      NA       1    2
15      NA       1       0    1
16      NA       1       2    1
17      NA       1       3    1
18      NA       2       0    1
19      NA       2       1    1