Cross-tabulate
Cross-tabulate the layers of a SpatRaster to create a contingency table. methods
## S4 method for signature 'SpatRaster,missing' crosstab(x, digits=0, long=FALSE, useNA=FALSE)
x
: SpatRasterdigits
: integer. The number of digits for rounding the values before cross-tabulationlong
: logical. If TRUE
the results are returned in 'long' format data.frame instead of a tableuseNA
: logical, indicting if the table should includes counts of NA
valuesA table or data.frame
freq
, zonal
r <- s <- rast(nc=5, nr=5) set.seed(1) values(r) <- runif(ncell(r)) * 2 values(s) <- runif(ncell(r)) * 3 x <- c(r, s) crosstab(x) rs <- r/s r[1:5] <- NA s[20:25] <- NA x <- c(r, s, rs) crosstab(x, useNA=TRUE, long=TRUE)
Useful links