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 sp
classe pour mes objets ( SpatialPointsDataFrame
, etc.).
point.in.polygon
est un bon début pour les polygones de base, bien qu'un peu grossier (par exemple):point.in.polygon(pnt@coords[,1], pnt@coords[,2], lake@polygons[[1]]@Polygons[[1]]@coords[,1], lake@polygons[[1]]@Polygons[[1]]@coords[,2])
je suis vraiment impatient de cergeos
que vous avez mentionné, et a tous les prédicats binaires et plus encore.