Eu usei a gIntersect
função do rgeos
pacote para determinar quais polígonos se cruzam em dois SpatialPolygonsDataFrames (spdfs). O resultado é uma matriz lógica com nomes de nomes que compreendem os IDs de polígonos do sfdf # 1 e os nomes de nomes dos IDs de polígonos do spdf # 2. Eu gostaria de usar essa matriz para o subconjunto das duas spdfs originais, omitindo esses polígonos que têm nenhuma sobreposição, antes de usar a função gIntersection
. Posso facilmente converter a matriz lógica em dois vetores de IDs (um para cada spdf) que representam esses polígonos com pelo menos alguma sobreposição.
Minha pergunta é - posso usar um vetor de IDs de polígono para definir um spdf? Em caso afirmativo, ele subconecta corretamente os valores no @data
slot? Eu sei que geralmente posso agrupar spdfs usando esta sintaxe:
spdf[spdf@data$myVariable[someVector], ]
Isso não funcionará para IDs, pois eles estão enterrados muito mais profundamente na estrutura da lista spdf.
spdf[r, ]
)