Comment puis-je utiliser les routes à partir des données OSM dans R?

Je dois utiliser les routes pour comparer les positions GPS pour vérifier qu'il n'y a pas d'erreurs dans les données.

J'ai téléchargé le map.osm à partir d'un site contenant des informations sur la ville de Rome et je l'ai transformé en lignes spatiales en R.

src <- osmsource_file(file = "rome.osm")
bb <- center_bbox(41.9, 12.3, 1000, 1000)
rome <- get_osm(bb, src)
ways <- find(rome, way(tags(k == "highway")))
ways <- find_down(rome, way(ways))
ways <- subset(rome, ids = ways)
hw_lines <- as_sp(ways, "lines")  

Donc, avec l'objet Lignes spatiales, j'ai essayé d'utiliser la fonction distm mais on l'accuse que l'objet spatial ne peut pas être utilisé (il n'est pas dans un format approprié).

Il existe un autre moyen d'utiliser OSM pour obtenir des informations sur les routes uniquement que je peux comparer avec les données GPS dans R?

J'ai téléchargé un ensemble de données de la ville de Rome en utilisant le osmarpackage. Après cela, j'ai suivi votre requête pour obtenir les autoroutes souhaitées , puis j'ai construit un objet SpatialPoints aléatoire à l' intérieur des données récupérées de la boîte englobante de Rome. J'ai mesuré les distances entre les points et les lignes en utilisant la fonction dist2linedu geospherepackage.

Veuillez essayer le code ci-dessous:

# Load libraries
library('osmar')
library('geosphere')

# Define the spatial extend of the OSM data we want to retrieve
rome.box <- center_bbox(center_lon = 12.5450, center_lat = 41.8992, width =  2000, height = 2000)

# Download all osm data inside this area
api <- osmsource_api()
rome <- get_osm(rome.box, source = api)

# General plot
plot(rome)

# Find highways
ways <- find(rome, way(tags(k == "highway")))
ways <- find_down(rome, way(ways))
ways <- subset(rome, ids = ways)

# SpatialLinesDataFrame object
hw_lines <- as_sp(ways, "lines")  

# Plot
spplot(hw_lines, zcol = "uid")

# Interactive view
mapview::mapview(hw_lines) 

# Make a random points dataset (like GPS)
gpsPoints <- spsample(x = hw_lines, n = 100, type = "random")

# Plot points
plot(hw_lines, xlab = "Lon", ylab = "Lat")
plot(gpsPoints, add = TRUE, pch = 19, col = "red")
box()

# Distances between Higways and random points
distances <- dist2Line(p = gpsPoints, line = hw_lines)