J'ai deux classes partageant le même CRS (Latitute et Longitude):
bolognaQuartieriMap
: unSpatialPolygonDataFrame
contenant des données d'un arrondissement de la ville.crashPoints
: unSpatialPointsDataFrame
contenant des données d'accidents.
Ils sont bien tracés en utilisant:
plot(bolognaQuartieriMap)
title("Crash per quartiere")
plot(crashPoints, col="red",add=TRUE)
Ce dont j'ai besoin, c'est d'obtenir le nombre de points ( crashPoints
) dans chaque polygone qui le constitue bolognaQuartieriMap
. On m'a suggéré d'utiliser over()
mais je n'ai pas réussi.