Estoy tratando de crear polígonos de Voronoi (también conocidos como teselaciones de Dirichlet o polígonos de Thiessen) dentro de una región geográfica fija para un conjunto de puntos. Sin embargo, tengo problemas para encontrar un método en R que limite los polígonos dentro de los bordes del mapa. Mi objetivo principal es obtener cálculos de área precisos (no simplemente para producir un gráfico visual). Por ejemplo, lo siguiente comunica visualmente lo que estoy tratando de lograr:

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)

que produce lo siguiente:

Conceptualmente quiero cruzarcounties convt que debe proporcionar un conjunto de polígonos delimitados por las fronteras del condado y cálculos precisos del área para cada uno. Ahora mismo,vt$summary proporciona cálculos de área para cada polígono, pero obviamente están exagerados para todos menos un polígono interior, ydeldir() parece aceptar solo recintos rectangulares para surw argumento. Soy nuevo en las capacidades geoespaciales de R, así que estoy abierto a otros enfoques más allá de lo que describí anteriormente.