Просторові оператори в R?

Як я можу зробити просторові виділення або фільтри в R? Я можу це зробити в PostGIS, використовуючи (наприклад):

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

і це аналогічно шматок пирога в більшості інших програм ГІС. Тож я очікую, що щось у R виглядатиме так:

mypoint[within(mypoint, mypoly),]

але я ніде не можу знайти жодного відвідника. Я шукаю те, що використовує загальну просторову термінологію предикатів, що використовується в PostGIS / GEOS / JTS / Shapely (всередині, містить, перетинає, перекриває, торкається тощо; DWithin буде бонусом). Чи існує ця функціональність або бібліотека для R? Якщо ні, чи є щось подібне?

Примітка. У R я використовую spклас для своїх об'єктів ( SpatialPointsDataFrameтощо).

Для вашого конкретного запитання ви повинні мати можливість використовувати point.in.polygon( документи ) або overlay( віньєтку ) для проведення аналізу.

Просторова мова предикатів не тривіальна для реалізації, і всі згадані вами приклади зводяться до портів або обгортки JTS в якийсь момент. Існує попередня альфа-версія rgeos , нещодавній запис Google SOC, який забезпечує доступ до GEOS безпосередньо з Р.

Загалом, здається, що розробники R стикаються з тим, чи хочуть вони підтримувати повну реалізацію ГІС в межах самого R, або замість цього покладаються на полегшення пакетів, які дозволяють взаємодію з PostGIS, GRASS ( spgrass), SAGA ( rsaga) тощо надавати більш глибокі просторовий кмітливість.