Ich versuche herauszufinden, wie man (in GIS so einfach) Operationen in R durchführt.

Nehmen wir ein Beispiel für einen Polygondatensatz ausspdep Paket

library("spdep")
c <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1])
plot(c)

Ich habe es geschafft herauszufinden, dass ich mit Polygonen mit logischen Anweisungen auswählen kannsubset. Zum Beispiel:

cc <- Teilmenge (c, c @ data $ POLYID <5) Diagramm (cc)

Angenommen, ich habe einen anderen Datenrahmen, den ich mit meinen Geodaten verknüpfen möchte:

POLYID=1:9
TO.LINK =101:109
link.data <- data.frame(POLYID=POLYID, TO.LINK=TO.LINK)

Wie kann ich mithilfe dieser beiden Datensätze zwei räumliche Datenrahmen abrufen:

Erstens, bestehend aus Polygonen, deren ID im zweiten Datenrahmen enthalten istZweitens, bestehend aus den entgegengesetzten Mengenpolygonen, die im zweiten Datenrahmen nicht existieren.

Wie könnte ich an diesen Punkt gelangen?