【R】Braunschweigの交通事故

2020年11月4日

Braunschweigでの交通事故（unfall）のデータを取得して地図上にプロットしてみます。

データは、Statische Aemter des Bundes und der LaenderのUnfallatlasサイトから、２０１９年のデータ「Download Unfallorte 2019 – shapefile (zip)」をダウンロードして使用します。

Braunschweigのde:amtlicher_gemeindeschluesselは、03101000　です。OpenStreetMapで確認しました。

地図のシェープファイルは別途用意しておきます。

Alexandra Kapp氏のコードを参考にしました。

library(sf)

library(tidyverse)

library(osmdata)

library(mapview)

accidents <- st_read("./Shapefile/Unfallorte2019_LinRef.shp") %>%

st_zm()

BS <- opq(bbox = 'braunschweig germany') %>%

add_osm_feature(key = 'de:amtlicher_gemeindeschluessel', value = '03101000') %>%

osmdata_sf ()

BS <- BS$osm_multipolygons$geometry

accidents_BS_index <- st_intersects(accidents, st_transform(BS, 25832), sparse = F)

accidents_BS <- accidents[accidents_BS_index, ] %>% st_transform(4326)

coords <- st_coordinates(accidents_BS)

accidents_BS["lng"] <- coords[,1]

accidents_BS["lat"] <- coords[,2]

map <- read_sf("./Shapefile2/DEU_adm3.shp")

BS_map <- map %>%

filter(NAME_2 == "Braunschweig")

mapview(BS_map, legend = FALSE) +

mapview(accidents_BS, color = "red", size = 0.3, legend = FALSE)

library(sf) library(tidyverse) library(osmdata) library(mapview) accidents <- st_read("./Shapefile/Unfallorte2019_LinRef.shp") %>% st_zm() BS <- opq(bbox = 'braunschweig germany') %>% add_osm_feature(key = 'de:amtlicher_gemeindeschluessel', value = '03101000') %>% osmdata_sf () BS <- BS$osm_multipolygons$geometry accidents_BS_index <- st_intersects(accidents, st_transform(BS, 25832), sparse = F) accidents_BS <- accidents[accidents_BS_index, ] %>% st_transform(4326) coords <- st_coordinates(accidents_BS) accidents_BS["lng"] <- coords[,1] accidents_BS["lat"] <- coords[,2] map <- read_sf("./Shapefile2/DEU_adm3.shp") BS_map <- map %>% filter(NAME_2 == "Braunschweig") mapview(BS_map, legend = FALSE) + mapview(accidents_BS, color = "red", size = 0.3, legend = FALSE)

library(sf)
library(tidyverse)
library(osmdata)
library(mapview)

accidents <- st_read("./Shapefile/Unfallorte2019_LinRef.shp") %>% 
  st_zm()
BS <-   opq(bbox = 'braunschweig germany') %>%
  add_osm_feature(key = 'de:amtlicher_gemeindeschluessel', value = '03101000') %>% 
  osmdata_sf ()
BS <- BS$osm_multipolygons$geometry

accidents_BS_index <- st_intersects(accidents, st_transform(BS, 25832), sparse = F) 

accidents_BS <- accidents[accidents_BS_index, ] %>% st_transform(4326)
coords <- st_coordinates(accidents_BS)
accidents_BS["lng"] <- coords[,1]
accidents_BS["lat"] <- coords[,2]

map <- read_sf("./Shapefile2/DEU_adm3.shp")
BS_map <- map %>%
  filter(NAME_2 == "Braunschweig") 

mapview(BS_map, legend = FALSE) +
  mapview(accidents_BS, color = "red", size = 0.3, legend = FALSE)

mapviewでインタラクティブな地図を表示しました。