千葉コロナマップ

千葉日報の 市町村別の感染者数（一覧表） | 千葉日報オンライン  が毎日更新されており、そのデータをCSVファイルとして Gist に保存してある。

　　ファイル名：covid_chiba.csv

これの最終日のデータをマップとしてプロットする。

使用する千葉県のシェープファイルは geojson 形式で、筆者の Dropbox に置いてある。

　　URL：https://www.dropbox.com/s/4310z8hktjijzol/chiba_admin.geojson?dl=0 

R スクリプト：

library(sf)
library(dplyr)
library(ggplot2)
library(ggthemes)

chiba <- st_read("https://www.dropbox.com/s/4310z8hktjijzol/chiba_admin.geojson?dl=1")
df <- read.csv("df <- read.csv("https://gist.githubusercontent.com/fusion0202/fdc1bfb603603b9dd8b42930483bbdac/raw/58e650e9087f2d121d19f1088335bd6c82fb867f/covid_chiba.csv", check.names = FALSE)
Date <- tail(colnames(df), 1)
df %>% rename(cases = Date) %>%
   select(sichoson, cases) %>%
   mutate(cases = replace(cases, cases == 0, NA)) -> dt

left_join(chiba, dt, by = c("SIKUCHOSON" = "sichoson")) %>%
   ggplot() +
   geom_sf(aes(fill = cases), alpha = 0.8, colour = 'grey5', size = 0.1) +
   scale_fill_gradient(low = "#fef9f9", high = "#cd0505", na.value = "white", name = "No. of Cases") +
   geom_sf_text(aes(label = cases)) +
   theme_map() +
   theme(legend.position = c(0.80, 0.05)) +
   labs(title = "COVID-19 Cases in Chiba",
           subtitle = paste("Data as of ", Date),
           caption = "Data Source: https://www.chibanippo.co.jp/news/national/681627")

 プロット例：