Contiguity index (Shape metric)
lsm_p_contig(landscape, directions = 8)tibble
A categorical raster object: SpatRaster; Raster* Layer, Stack, Brick; stars or a list of SpatRasters.
The number of directions in which patches should be connected: 4 (rook's case) or 8 (queen's case).
$$CONTIG = \frac{\Bigg[\frac{\sum\limits_{r=1}^z c_{ijr}}{a_{ij}}\Bigg] - 1 }{ v - 1} $$
where \(c_{ijr}\) is the contiguity value for pixel r in patch ij, \(a_{ij}\) the area of the respective patch (number of cells) and \(v\) is the size of the filter matrix (13 in this case).
CONTIG is a 'Shape metric'. It asses the spatial connectedness (contiguity) of cells in patches. CONTIG coerces patch values to a value of 1 and the background to NA. A nine cell focal filter matrix:
filter_matrix <- matrix(c(1, 2, 1,
2, 1, 2,
1, 2, 1), 3, 3, byrow = T)
... is then used to weight orthogonally contiguous pixels more heavily than diagonally contiguous pixels. Therefore, larger and more connections between patch cells in the rookie case result in larger contiguity index values.
None
0 >= CONTIG <= 1
Equals 0 for one-pixel patches and increases to a limit of 1 (fully connected patch).
McGarigal K., SA Cushman, and E Ene. 2023. FRAGSTATS v4: Spatial Pattern Analysis Program for Categorical Maps. Computer software program produced by the authors; available at the following web site: https://www.fragstats.org
LaGro, J. 1991. Assessing patch shape in landscape mosaics. Photogrammetric Engineering and Remote Sensing, 57(3), 285-293
lsm_c_contig_mn,
lsm_c_contig_sd,
lsm_c_contig_cv,
lsm_l_contig_mn,
lsm_l_contig_sd,
lsm_l_contig_cv
landscape <- terra::rast(landscapemetrics::landscape)
lsm_p_contig(landscape)
Run the code above in your browser using DataLab