Я пытаюсь понять, как выполнять (так легко в ГИС) операции в R.

Давайте рассмотрим пример набора данных многоугольника изspdep пакет

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

Мне удалось выяснить, что я могу выбирать полигоны с логическими утверждениями, используяsubset, Например:

cc <- участок подмножества (c, c @ data $ POLYID <5) (cc)

Теперь предположим, что у меня есть еще один фрейм данных, который я хотел бы присоединить к своим пространственным данным:

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

Используя эти два набора данных, как я могу получить два кадра пространственных данных:

Во-первых, состоящий из многоугольников, которые имеют свой идентификатор во втором фрейме данныхВторое, состоящее из противоположного набора - многоугольников, которых нет во втором фрейме данных.

Как я мог добраться до этой точки?