Opérateurs spatiaux dans R?

Comment puis-je faire des sélections spatiales ou des filtres dans R? Je peux le faire dans PostGIS en utilisant (par exemple):

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

et c'est également un jeu d'enfant dans la plupart des autres logiciels SIG. Je m'attends donc à ce que quelque chose dans R ressemble à:

mypoint[within(mypoint, mypoly),]

mais je ne trouve aucune piste nulle part. Je cherche quelque chose qui utilise la terminologie commune de prédicat spatial utilisée dans PostGIS / GEOS / JTS / Shapely (à l'intérieur, contient, intersecte, se chevauche, se touche, etc.; DWithin serait un bonus). Cette fonctionnalité ou bibliothèque existe-t-elle pour R? Sinon, y a-t-il quelque chose de similaire?

Remarque: Dans R, j'utilise la spclasse pour mes objets ( SpatialPointsDataFrame, etc.).

Pour votre question spécifique, vous devriez pouvoir utiliser point.in.polygon( docs ) ou overlay( vignette ) pour effectuer l'analyse.

Un langage de prédicat spatial n'est pas trivial à implémenter, et tous les exemples que vous avez mentionnés se résument à des ports ou des wrappers de JTS à un moment donné. Il existe une version pré-alpha de rgeos , une entrée récente de Google SOC qui donne accès à GEOS directement à partir de R.

En général, il semble que les développeurs de R cherchent à maintenir une implémentation SIG complète dans R lui-même, ou plutôt de s'appuyer sur des packages facilitant l'interaction avec PostGIS, GRASS ( spgrass), SAGA ( rsaga) et similaires pour fournir des informations plus approfondies. savoir-faire spatial.