Saya telah menggunakan gIntersect
fungsi dari rgeos
paket untuk menentukan poligon yang bersilangan di dua SpatialPolygonsDataFrames (spdfs). Hasilnya adalah matriks logis dengan rownames yang terdiri dari ID poligon dari sfdf # 1 dan colnames ID poligon dari spdf # 2. Saya ingin menggunakan matriks ini untuk mengelompokkan dua spdf asli, menghilangkan poligon yang tidak memiliki tumpang tindih, sebelum menggunakan fungsi gIntersection
. Saya dapat dengan mudah mengkonversi matriks logis menjadi dua vektor ID (satu untuk setiap spdf) yang mewakili poligon dengan setidaknya beberapa tumpang tindih.
Pertanyaan saya adalah - dapatkah saya menggunakan vektor ID poligon untuk subset sebuah spdf? Jika demikian, apakah ini akan dengan benar mengelompokkan nilai dalam @data
slot? Saya tahu saya biasanya dapat mengelompokkan spdf menggunakan sintaks ini:
spdf[spdf@data$myVariable[someVector], ]
Ini tidak akan berfungsi untuk ID, karena mereka terkubur lebih dalam di dalam struktur daftar spdf.
spdf[r, ]
)