Em R, como participo e participo o SpatialPolygonsDataFrame?
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?