Landscape metrics of categorical rasters (FRAGSTATS, landscapemetrics conventions)

PatchMetrics: area (ha), perim (m, boundary included), shape \(0.25\,p/\sqrt{a}\), para, enn (nearest cell-centre distance to the class's other patches). ClassMetrics: ca, pland, np, pd, lpi, ed (boundary excluded), lsi, ai, clumpy (FRAGSTATS definition), cohesion, iji, area_mn, shape_mn, enn_mn. LandscapeMetrics: ta, np, pd, lpi, ed, lsi, ai, contag (adjacency proportions), iji, pr, shdi, shei, sidi, siei, dominance \(\ln m - SHDI\), shape_mn, area_mn, enn_mn (McGarigal and Marks 1995; Li and Reynolds 1993; O'Neill et al. 1988; Hesselbarth et al. 2019). Identical to the Python arm morie.fn.lsmets.

Usage

PatchMetrics(landscape, res = 1, directions = 8)

ClassMetrics(landscape, res = 1, directions = 8)

LandscapeMetrics(landscape, res = 1, directions = 8)

Arguments

landscape: Integer matrix of classes (NA allowed).
res: Cell size (m).
directions: 4 or 8 (patch connectivity).

Value

List of metrics.

References

McGarigal, K. and Marks, B. J. (1995). FRAGSTATS. USDA Forest Service General Technical Report PNW-351.

Hesselbarth, M. H. K., Sciaini, M., With, K. A., Wiegand, K. and Nowosad, J. (2019). landscapemetrics: an open-source R tool to calculate landscape metrics. Ecography 42, 1648-1657.

Examples

g <- rbind(c(1, 1, 2), c(2, 1, 2))
PatchMetrics(g, directions = 4)$shape
#> [1] 1.154701 1.000000 1.060660
ClassMetrics(g, directions = 4)$pland
#> [1] 50 50
LandscapeMetrics(g, directions = 4)$shdi
#> [1] 0.6931472