Create a weighted co-occurrence matrix (wecoma)
get_wecoma( x, w, neighbourhood = 4, classes = NULL, fun = "mean", na_action = "replace" )
x | A matrix with categories |
---|---|
w | A matrix with weights |
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 vector or a list with the values of selected classes from the |
fun | Function to calculate values from adjacent cells to contribute to output matrix, |
na_action | Decides on how to behave in the presence of missing values in |
A weighted co-occurrence matrix
library(comat) data(raster_x, package = "comat") data(raster_w, package = "comat") wom = get_wecoma(raster_x, raster_w) wom#> 1 2 3 #> 1 12.0 5.0 13.5 #> 2 5.0 12.0 14.5 #> 3 13.5 14.5 49.0#> 1 3 #> 1 12.0 13.5 #> 3 13.5 49.0