Ich versuche, Voronoi-Polygone (auch Dirichlet-Tessellationen oder Thiessen-Polygone genannt) innerhalb einer festgelegten geografischen Region für eine Reihe von Punkten zu erstellen. Ich habe jedoch Probleme, eine Methode in R zu finden, die die Polygone innerhalb der Kartengrenzen begrenzt. Mein Hauptziel ist es, genaue Flächenberechnungen zu erhalten (nicht nur um ein visuelles Diagramm zu erstellen). Das folgende Beispiel vermittelt visuell, was ich erreichen möchte:

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)

was folgendes ergibt:

Konzeptionell möchte ich mich schneidencounties mitvt Dies sollte eine Reihe von Polygonen bereitstellen, die durch die Bezirksgrenzen begrenzt sind, und genaue Flächenberechnungen für jedes. Jetzt sofort,vt$summary liefert Flächenberechnungen für jedes Polygon, aber sie sind offensichtlich für alle außer dem einen inneren Polygon überbewertetdeldir() scheint nur rechteckige Gehäuse zu akzeptierenrw Streit. Ich bin neu in Rs geospacial Fähigkeiten, bin also offen für andere Ansätze als die, die ich oben skizziert habe.