У меня есть шейп-файл земного покрова с двумя категориями полигонов: лес [1] и вырубка леса [0]. Я хочу создать гексагональную сетку только на полигонах обезлесения, чтобы вычислить метрики связности. Идея состоит в том, чтобы рассчитать, насколько добавление 1 га леса улучшит ландшафтную связь. Я использую функции "spsample" и "HexPoints2SpatialPolygon" из пакета sp в R. Однако я могу найти только один способ сделать это для всего шейп-файла, а не только для категории обезлесения. Кроме того, при выполнении этой функции для всего шейп-файла выходные данные не различают категории
>####cellsize parameter, which is the distance (cell_diameter; i.e., apothem) between centroids of neighbouring cells
>####defining cellsize by area in square meter (1ha == 10000m2)
>cell_area <- 10000
>cell_diameter <- sqrt(2 * cell_area / sqrt(3))
>candidateCT.shpHexpt <- spsample(candidateCT.shpUTM, type = "hexagonal", cellsize = cell_diameter, offset = c(0, 0))
>candidateCT.shpHexpol <- HexPoints2SpatialPolygons(candidateCT.shpHexpt)
Любая помощь будет принята с благодарностью