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

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