Como posso usar estradas a partir de dados OSM no R?

Preciso usar as estradas para comparar as posições do GPS e verificar se não há erros nos dados.

Eu baixei o map.osm de um site que tinha informações da cidade de Roma e o transformou em Spatial Lines in 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")  

Então, com o Objeto de Linhas Espaciais, tentei usar a função distm, mas é acusado de que o objeto espacial não pode ser usado (não está em um formato adequado).

Existe outra maneira de usar o OSM para obter informações apenas de estradas que eu posso comparar com dados de GPS no R?

Eu baixei um conjunto de dados da cidade de Roma usando o osmarpacote. Depois disso, segui sua consulta para obter as rodovias desejadas e depois construir um objeto SpatialPoints aleatório dentro dos dados recuperados da caixa delimitadora de Roma. Eu medi distâncias entre pontos e linhas usando a função dist2linedo geospherepacote.

Por favor, tente o código abaixo:

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