Làm cách nào tôi có thể sử dụng đường từ dữ liệu OSM trong R?

Tôi cần sử dụng các con đường để so sánh các vị trí GPS để xác minh rằng chúng không có lỗi trong dữ liệu.

Tôi đã tải xuống map.osm từ một trang web có thông tin về Thành phố Rome và nó đã được chuyển đổi thành Spatial Lines ở 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")  

Vì vậy, với Đối tượng Đường không gian, tôi đã cố gắng sử dụng chức năng distm nhưng bị cáo buộc rằng đối tượng không gian không thể được sử dụng (không phải ở định dạng phù hợp).

Có một cách khác sử dụng OSM để lấy thông tin chỉ những con đường mà tôi có thể so sánh với dữ liệu GPS trong R?

Tôi đã tải xuống một bộ dữ liệu từ Thành phố Rome bằng cách sử dụng osmargói. Sau đó, tôi đã theo dõi truy vấn của bạn để có được các đường cao tốc mong muốn và sau đó xây dựng một đối tượng SpatialPoints ngẫu nhiên bên trong hộp giới hạn lấy dữ liệu ở Rome. Tôi đo khoảng cách giữa các điểm và đường bằng cách sử dụng hàm dist2linetừ geospheregói.

Vui lòng thử mã dưới đây:

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