Polígonos do diagrama de Voronoi entre fronteiras geográficas

Estou tentando criar polígonos de Voronoi (também conhecidos como mosaicos de Dirichlet ou polígonos de Thiessen) em uma região geográfica fixa para um conjunto de pontos. No entanto, estou tendo problemas para encontrar um método em R que limite os polígonos dentro das bordas do mapa. Meu principal objetivo é obter cálculos de área precisos (não apenas para produzir um gráfico visual). Por exemplo, o seguinte comunica visualmente o que estou tentando alcançar:

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 produz o seguinte:

Conceitualmente, quero cruzarcounties comvt que deve fornecer um conjunto de polígonos delimitados pelas fronteiras do condado e cálculos precisos de área para cada um. Agora mesmo,vt$summary fornece cálculos de área para cada polígono, mas eles são obviamente superestimados para todos, exceto para o polígono interno, edeldir() parece aceitar apenas gabinetes retangulares para suarw argumento. Sou novo nas capacidades geoespaciais de R, portanto estou aberto a outras abordagens além do que descrevi acima.