¿Cómo puedo usar carreteras de datos OSM en R?

Necesito usar las carreteras para comparar las posiciones de GPS para verificar que no haya errores en los datos.

Descargué el map.osm de un sitio que tiene información de la ciudad de Roma y lo transformé en Spatial Lines 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")  

Entonces, con el objeto Spatial Lines intenté usar la función distm pero se acusa de que el objeto espacial no se puede usar (no está en un formato adecuado).

¿Hay otra forma de usar OSM para obtener información de solo carreteras que puedo comparar con datos de GPS en R?

Descargué un conjunto de datos de la ciudad de Roma usando el osmarpaquete. Después de eso, seguí su consulta para obtener las carreteras deseadas y luego construí un objeto SpatialPoints aleatorio dentro de los datos recuperados del cuadro delimitador de Roma. Medí distancias entre puntos y líneas usando la función dist2linedel geospherepaquete.

Por favor, intente el siguiente código:

# 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)