Tôi đã sử dụng gIntersecthàm từ rgeosgói để xác định đa giác nào giao nhau trong hai SpatialPolygonsDataFrames (spdfs). Kết quả là một ma trận logic với các tên gọi bao gồm các ID đa giác từ sfdf # 1 và đặt tên cho các ID đa giác từ spdf # 2. Tôi muốn sử dụng ma trận này để tập hợp hai spdf gốc, bỏ qua các đa giác không có sự trùng lặp trước khi sử dụng hàm gIntersection. Tôi có thể dễ dàng chuyển đổi ma trận logic thành hai vectơ ID (một cho mỗi spdf) đại diện cho các đa giác đó với ít nhất một số chồng chéo.
Câu hỏi của tôi là - tôi có thể sử dụng một vectơ ID đa giác để đặt một spdf không? Nếu vậy, nó sẽ tập hợp chính xác các giá trị trong @datakhe? Tôi biết tôi thường có thể tập hợp các spdfs bằng cú pháp này:
spdf[spdf@data$myVariable[someVector], ]Điều này sẽ không hoạt động đối với ID, vì chúng được chôn sâu hơn nhiều trong cấu trúc danh sách spdf.
spdf[r, ])