Create an integrated co-occurrence matrix (incoma)

get_incoma(x, neighbourhood = 4, classes = NULL)

Arguments

x

A list object containing categorical matrices 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 the same length as x with the values of selected classes from all of the objects in x. It is used to calculate incoma only for selected classes.

Value

An integrated co-occurrence matrix

Examples

data(raster_x, package = "comat") data(raster_w, package = "comat") x = list(raster_x, raster_w, raster_x) get_incoma(x)
#> 1 2 3 2 4 6 8 9 1 2 3 #> 1 4 1 3 3 3 1 0 1 4 1 3 #> 2 1 2 2 0 2 1 1 1 1 2 2 #> 3 3 2 6 2 1 1 2 5 3 2 6 #> 2 3 0 2 2 1 1 0 1 3 0 2 #> 4 3 2 1 1 0 2 2 1 3 2 1 #> 6 1 1 1 1 2 0 0 0 1 1 1 #> 8 0 1 2 0 2 0 0 1 0 1 2 #> 9 1 1 5 1 1 0 1 4 1 1 5 #> 1 4 1 3 3 3 1 0 1 4 1 3 #> 2 1 2 2 0 2 1 1 1 1 2 2 #> 3 3 2 6 2 1 1 2 5 3 2 6 #> attr(,"no_unique") #> [1] 3 5 3
get_incoma(x, classes = list(1:2, 2:4, 1))
#> 1 2 2 3 4 1 #> 1 4 1 3 0 3 4 #> 2 1 2 0 0 2 1 #> 2 3 0 2 0 1 3 #> 3 0 0 0 0 0 0 #> 4 3 2 1 0 0 3 #> 1 4 1 3 0 3 4 #> attr(,"no_unique") #> [1] 2 3 1