Skip to contents

Landscape division index (Aggregation metric)

Usage

lsm_c_division(landscape, directions = 8)

Arguments

landscape

A categorical raster object: SpatRaster; Raster* Layer, Stack, Brick; stars or a list of SpatRasters.

directions

The number of directions in which patches should be connected: 4 (rook's case) or 8 (queen's case).

Value

tibble

Details

$$DIVISON = (1 - \sum \limits_{j = 1}^{n} (\frac{a_{ij}} {A}) ^ 2)$$ where \(a_{ij}\) is the area in square meters and \(A\) is the total landscape area in square meters.

DIVISION is an 'Aggregation metric. It can be in as the probability that two randomly selected cells are not located in the same patch of class i. The landscape division index is negatively correlated with the effective mesh size (lsm_c_mesh).

Because the metric is based on distances or areas please make sure your data is valid using check_landscape.

Units

Proportion

Ranges

0 <= Division < 1

Behaviour

Equals DIVISION = 0 if only one patch is present. Approaches DIVISION = 1 if all patches of class i are single cells.

References

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

Jaeger, J. A. 2000. Landscape division, splitting index, and effective mesh size: new measures of landscape fragmentation. Landscape ecology, 15(2), 115-130.

Examples

landscape <- terra::rast(landscapemetrics::landscape)
lsm_c_division(landscape)
#> # A tibble: 3 × 6
#>   layer level class    id metric   value
#>   <int> <chr> <int> <int> <chr>    <dbl>
#> 1     1 class     1    NA division 0.985
#> 2     1 class     2    NA division 0.987
#> 3     1 class     3    NA division 0.924