Wie kann ich Straßen aus OSM-Daten in R verwenden?

Ich muss die Straßen verwenden, um GPS-Positionen zu vergleichen, um sicherzustellen, dass sie keine Fehler in den Daten sind.

Ich habe die map.osm von einer Site heruntergeladen, die Informationen über die Stadt Rom enthält, und sie in R in Spatial Lines umgewandelt.

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

Also habe ich mit dem Spatial Lines Object versucht, die distm-Funktion zu verwenden, aber es wird beschuldigt, dass das Spatial Object nicht verwendet werden kann (nicht in einem geeigneten Format).

Gibt es eine andere Möglichkeit, mit OSM nur Informationen zu Straßen abzurufen, die ich mit GPS-Daten in R vergleichen kann?

Ich habe mit dem osmarPaket einen Datensatz aus der Stadt Rom heruntergeladen . Danach folgte ich Ihrer Anfrage, um die gewünschten Autobahnen zu erhalten und dann ein zufälliges SpatialPoints- Objekt innerhalb des Rom-Begrenzungsrahmens zu erstellen , in dem Daten abgerufen wurden. Ich habe Entfernungen zwischen Punkten und Linien mit der Funktion dist2lineaus dem geospherePaket gemessen .

Bitte versuchen Sie den folgenden Code:

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