Eu estou tentando descobrir o meu caminho em como executar (tão fácil em GIS) operações em R.

Vamos dar um exemplo de conjunto de dados de polígonospdep pacote

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

Eu consegui descobrir que eu posso escolher polígonos com instruções lógicas usandosubset. Por exemplo:

plotagem cc <- subconjunto (c, c @ data $ POLYID <5) (cc)

Agora, vamos supor que eu tenha outro quadro de dados que gostaria de unir aos meus dados espaciais:

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

Usando esses dois conjuntos de dados, como posso obter dois quadros de dados espaciais:

Primeiro, consistindo de polígonos que têm seu ID no segundo quadro de dadosSegundo, consistindo no conjunto oposto - polígonos que não existem no segundo quadro de dados.

Como eu poderia chegar a este ponto?