패키지 의 gIntersect함수를 사용하여 rgeos두 개의 SpatialPolygonsDataFrames (spdfs)에서 교차하는 다각형을 결정했습니다. 결과는 sfdf # 1의 폴리곤 ID를 포함하는 행 이름과 spdf # 2의 폴리곤 ID를 colnames하는 논리 행렬입니다. 이 행렬을 사용하여 함수를 사용하기 전에 겹치지 않는 다각형을 생략하고 두 개의 원본 spdf를 하위 세트로 만들고 싶습니다 gIntersection. 논리 행렬을 적어도 일부 겹치는 다각형을 나타내는 두 개의 ID 벡터 (각 spdf마다 하나씩)로 쉽게 변환 할 수 있습니다.
내 질문은-다각형 ID 벡터를 사용하여 spdf의 하위 집합을 사용할 수 있습니까? 그렇다면 @data슬롯의 값을 올바르게 서브 세트 합니까? 이 구문을 사용하여 일반적으로 spdf를 하위 집합으로 만들 수 있다는 것을 알고 있습니다.
spdf[spdf@data$myVariable[someVector], ]spdf 목록 구조에 훨씬 더 깊이 묻혀 있기 때문에 ID에는 작동하지 않습니다.
spdf[r, ])