Create an integrated co-occurrence matrix (incoma)

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