Create a co-located co-occurrence matrix (cocoma)

get_cocoma(x, y, neighbourhood = 4, classes = NULL)

Arguments

x

A matrix with categories

y

A matrix with categories

neighbourhood

The number of directions in which cell adjacencies are considered as neighbours: 4 (rook's case) or 8 (queen's case). The default is 4.

classes

A list of length 2 with the values of selected classes from the x and y objects. It is used to calculate cocoma only for selected classes.

Value

A co-located co-occurrence matrix

Examples

library(comat) data(raster_x, package = "comat") data(raster_x_na, package = "comat") coom = get_cocoma(raster_x, raster_x_na) coom
#> 1 2 3 #> 1 3 1 3 #> 2 0 2 1 #> 3 2 2 2
get_cocoma(raster_x, raster_x_na, classes = list(c(1, 2), 3))
#> 3 #> 1 3 #> 2 1