Operator spasial di R?

Bagaimana saya bisa melakukan seleksi spasial atau filter dalam R? Saya dapat melakukan ini di PostGIS menggunakan (misalnya):

SELECT * FROM mypoint, mypoly
WHERE ST_Within(mypoint.geom, mypoly.geom);

dan juga sepotong kue di sebagian besar perangkat lunak GIS lainnya. Jadi saya berharap sesuatu dalam R terlihat seperti:

mypoint[within(mypoint, mypoly),]

tetapi saya tidak dapat menemukan petunjuk di mana pun. Saya mencari sesuatu yang menggunakan terminologi predikat spasial umum yang digunakan dalam PostGIS / GEOS / JTS / Shapely (di dalam, berisi, berpotongan, tumpang tindih, menyentuh, dll.; DWithin akan menjadi bonus). Apakah fungsi atau pustaka ini ada untuk R? Jika tidak, adakah yang serupa?

Catatan: Di R, saya menggunakan spkelas untuk objek saya ( SpatialPointsDataFrame, dll).

Untuk pertanyaan spesifik Anda, Anda harus dapat menggunakan point.in.polygon( dokumen ) atau overlay( sketsa ) untuk melakukan analisis.

Bahasa predikat spasial tidak sepele untuk diterapkan, dan semua contoh yang Anda sebutkan bermuara pada port atau pembungkus JTS di beberapa titik. Ada rilis pre-alpha dari rgeos , entri SOC Google baru-baru ini yang menyediakan akses ke GEOS langsung dari R.

Secara umum, tampaknya pengembang R bergulat dengan apakah mereka ingin mempertahankan implementasi GIS penuh dalam R itu sendiri, atau sebaliknya mengandalkan pada paket fasilitasi yang memungkinkan interaksi dengan PostGIS, GRASS ( spgrass), SAGA ( rsaga) dan sejenisnya untuk memberikan lebih dalam pengetahuan spasial.