【R】気象庁のデータから地震の震源地マップ

気象庁が様々なデータを公開しています。その多くは、2次的に利用しやすい形で提供されているものもあります。この中で、日本で多い地震のデータをスクレイピングして表示してみたいと思います。2020年4月30日の地震データを取得して地図上に表示することにします。

最初に結果を示すと、以下のようになります。マグニチュードの大きさに合わせて円の大きさを変化させています。

驚くほど多くの地震があることがわかります。多くは人が感じることはないのでしょうが。これを見ると日本は地震大国であることがわかりますね。

Rのコードは次の通りとなります。

library(leaflet)
library(knitr)
library(kableExtra)
library(dplyr)
library(tidyr)
library(rvest)
library(stringr)
library(clipr)

dat_jma <- html("https://www.data.jma.go.jp/svd/eqev/data/daily_map/20200430.html") #気象庁からデータ取得
dat<-dat_jma %>%
  html_nodes("pre") %>%
  html_text()

dat1<-as.data.frame(str_split(dat, pattern = "\n")) #一行ごとのデータフレームに
dat_lat<-as.data.frame(str_sub(dat1[,1], start=23, end=32)) #緯度抽出
dat_lon<-as.data.frame(str_sub(dat1[,1], start=34, end=45)) #経度抽出
dat_depth<-as.data.frame(str_sub(dat1[,1], start=47, end=51)) #深さ抽出
dat_M<-as.data.frame(str_sub(dat1[,1], start=53, end=56)) #M抽出

dat_lat_do<-as.data.frame(str_sub(dat_lat[,1], start=2, end=3))
dat_lat_fun<-as.data.frame(str_sub(dat_lat[,1], start=5, end=6))
dat_lat_byo<-as.data.frame(str_sub(dat_lat[,1], start=8, end=8))
dat_lat_N<-as.data.frame(as.numeric(as.character(dat_lat_do[,1]))+
  as.numeric(as.character(dat_lat_fun[,1]))/60+
  as.numeric(as.character(dat_lat_byo[,1]))/3600)

dat_lon_do<-as.data.frame(str_sub(dat_lon[,1], start=1, end=3))
dat_lon_fun<-as.data.frame(str_sub(dat_lon[,1], start=5, end=6))
dat_lon_byo<-as.data.frame(str_sub(dat_lon[,1], start=8, end=8))
dat_lon_E<-as.data.frame(as.numeric(as.character(dat_lon_do[,1]))+
  as.numeric(as.character(dat_lon_fun[,1]))/60+
  as.numeric(as.character(dat_lon_byo[,1]))/3600)

dat_earthquake<-data.frame("lat"=dat_lat_N, "lon"=dat_lon_E, "depth"=dat_depth, "M"=dat_M)
colnames(dat_earthquake)<-c("lat", "lon", "depth", "M")

dat_earthquake %>% 
  leaflet %>% 
  addTiles() %>% 
  addCircleMarkers(
    radius = ~2*as.numeric(as.character(M)),
    stroke = FALSE, fillOpacity = 0.5
  )

kable(dat_earthquake, align = "c") %>%
  kable_styling(full_width = F) %>%
  collapse_rows(columns = 1, valign = "middle")  %>%
  write_clip

長いですが、プロット用に修正されたデータフレームdat_earthquakeを次に示します。