Próbuję znaleźć sposób na wykonywanie (tak łatwo w GIS) operacji w R.

Weźmy przykładowy zestaw danych wielokąta zspdep pakiet

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

Udało mi się dowiedzieć, że mogę wybierać wielokąty za pomocą logicznych instrukcjisubset. Na przykład:

cc <- podzbiór (c, c @ dane $ POLYID <5) wykres (cc)

Załóżmy teraz, że mam inną ramkę danych, którą chciałbym dołączyć do moich danych przestrzennych:

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

Za pomocą tych dwóch zestawów danych mogę uzyskać dwie ramki danych przestrzennych:

Po pierwsze, składający się z wielokątów, które mają swój identyfikator w drugiej ramce danychPo drugie, składający się z przeciwnego zestawu - wielokątów, które nie istnieją w drugiej ramce danych.

Jak mogłem dojść do tego punktu?