Я пытаюсь создать полигоны Вороного (тесселяции Дирихле или полигоны Тиссена) в фиксированном географическом регионе для набора точек. Однако у меня возникают проблемы с поиском метода в R, который будет связывать полигоны в границах карты. Моя главная цель - получить точные расчеты площади (а не просто создать визуальный график). Например, следующее визуально сообщает, чего я пытаюсь достичь:

library(maps)
library(deldir)
data(countyMapEnv)
counties <- map('county', c('maryland,carroll','maryland,frederick', 'maryland,montgomery', 'maryland,howard'), interior=FALSE)
x <- c(-77.208703, -77.456582, -77.090600,  -77.035668, -77.197144)
y <- c(39.188603, 39.347019, 39.672818, 39.501898, 39.389203)
points(x,y)
vt <- deldir(x, y, rw=counties$range)
plot(vt, wlines="tess", lty="solid", add=TRUE)

который производит следующее:

Концептуально хочу пересечьcounties сvt который должен обеспечить набор полигонов, ограниченных границами округа, и точные расчеты площади для каждого. Прямо сейчас,vt$summary обеспечивает расчеты площади для каждого многоугольника, но они явно завышены для всех, кроме одного внутреннего многоугольника, иdeldir() кажется, принимает только прямоугольные оболочки для егоrw аргумент. Я новичок в геопространственных возможностях R, поэтому я открыт для других подходов, помимо того, что я изложил выше.