Saya telah menggunakan gIntersectfungsi dari rgeospaket 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 @dataslot? 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, ])