Я строю график, используя данные osm для визуализации дорожной сети города Эр-Рияд, и пытаюсь рассматривать автомобильные аварии на нем как точки. Но северная часть дорожной сети отсутствует, и я не могу расширить площадь участка.
Это код:
library(readxl)
library(dplyr)
library(osmdata)
q1 <- opq('Riyadh', format_out = 'polygon') %>%
add_osm_feature(key = 'highway', value = 'motorway')
cway_sev1 <- osmdata_sp(q1)
sp::plot(cway_sev1$osm_lines)
par(new=TRUE)
q2 <- opq('Riyadh') %>%
add_osm_feature(key = 'highway', value = 'primary')
cway_sev2 <- osmdata_sp(q2)
lines(cway_sev2$osm_lines, col = "red")
q3 <- opq('Riyadh') %>%
add_osm_feature(key = 'highway', value = 'secondary')
cway_sev3 <- osmdata_sp(q3)
lines(cway_sev3$osm_lines, col = "blue")
d<-read_excel('AccDataset.xlsm')
d1 = d[,9:10]
points(d1, pch = 20)
Это выходные данные сети OSM.
Это вывод после визуализации аварий