Bagaimana saya bisa menggunakan jalan dari data OSM di R?

Saya perlu menggunakan jalan untuk membandingkan posisi GPS untuk memverifikasi bahwa tidak ada kesalahan dalam data.

Saya mengunduh map.osm dari situs yang memiliki informasi Kota Roma dan mengubahnya menjadi Garis Spasial di 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")  

Jadi dengan Objek Garis Spasial saya mencoba menggunakan fungsi distm tetapi dituduh bahwa objek spasial tidak dapat digunakan (tidak dalam format yang sesuai).

Ada cara lain menggunakan OSM untuk mendapatkan informasi hanya jalan yang bisa saya bandingkan dengan data GPS di R?

Saya mengunduh dataset dari Kota Roma menggunakan osmarpaket. Setelah itu, saya mengikuti permintaan Anda untuk mendapatkan jalan raya yang diinginkan dan kemudian membangun objek SpatialPoints acak di dalam kotak pembatas Roma mengambil data. Saya mengukur jarak antara titik dan garis menggunakan fungsi dist2linedari geospherepaket.

Silakan coba kode di bawah ini:

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