台風の情報をRで描く

その4、21号です。前回からの変更点は図の下に書いてます:

前回からの変更点

  • 色を変更。今回はgrDevice::heat.colors()を使用しました。
  • 透明度を、中心気圧で変化させるよう設定。

感想

Rの色については結構難しいですね。選択肢が多すぎるのも難です。あと新規機能に関するネタが尽きてきました。

以下説明。

この内容は、気象庁の台風情報を、ある時点において取得して、そのまま図示するための一連のコードを実行しています。基本は私の練習のためです。

詳しい情報、および速報については、上記リンクで必ず確認してください。

基本的な流れ

データの取得

{rvest}を用いて上記Webページより必要なデータをスクレイピングします。

データの処理

この後{leaflet}にて描画するために、データを処理します。未熟なため非常に汚く効率が悪いです。だれか教えてください。

描画

{leaflet}にてプロットします。また必要に応じて事前にいくつかデータを準備してます。

コード

必要なパッケージをインストールして実行すれば再現できます。ただし、実行時の時刻でデータを取得しますので、台風の中心位置は再現されません。

# 台風情報
library("dplyr")
library("rvest")
library("stringr")
library("leaflet")

trg <- html("http://www.jma.go.jp/jp/typh/") %>% html_table
tbl <- trg[[4]] # テーブルは4番目が台風情報
res <- data.frame(
  timing = tbl$X1 %>% str_subset("[0-9][0-9]日[0-9][0-9]時"),
  lat = tbl$X2 %>% 
    str_subset("北緯") %>%
    str_extract("[0-9][0-9]\\.[0-9]度") %>%
    str_replace("度","") %>%
    as.numeric,
  lon = tbl$X2 %>%
    str_subset("東経") %>%
    str_extract("[0-9][0-9][0-9]\\.[0-9]度") %>%
    str_replace("度","") %>%
    as.numeric,
  hpa = tbl$X2 %>% 
    str_subset("hPa") %>% 
    str_replace("hPa","") %>% 
    as.numeric,
  c_wspeed = tbl %>% 
    dplyr::filter(X1 == "中心付近の最大風速") %>% 
    dplyr::select(X2) %>% 
    dplyr::rename(c_wspeed=X2),
  m_wspeed = tbl %>% 
    dplyr::filter(X1 == "最大瞬間風速") %>% 
    dplyr::select(X2) %>% 
    dplyr::rename(m_wspeed=X2)
)

# 取得した台風20号のデータを表示
knitr::kable(res)
timing lat lon hpa c_wspeed m_wspeed
<24日15時の実況> 18.7 132.9 985 25m/s(50kt) 35m/s(70kt)
<25日03時の予報> 18.8 132.1 980 30m/s(55kt) 40m/s(80kt)
<25日15時の予報> 19.1 131.5 975 30m/s(60kt) 45m/s(85kt)
<26日15時の予報> 21.5 129.3 970 35m/s(65kt) 50m/s(95kt)
<27日15時の予報> 23.0 127.1 965 35m/s(70kt) 50m/s(100kt) 
# leafletで描写
# カラーパレット作成
colpal <- colorNumeric(palette = grDevices::heat.colors(n=80), domain = c(930,1000))
# ポップアップ作成
res <- dplyr::mutate(res, popup = paste(timing,paste("中心気圧",hpa, "hPa"),paste("最大風速",c_wspeed),paste("瞬間最大風速",m_wspeed),sep = "<br/>"))
# 説明用テキスト
text <- paste(paste(Sys.time(),"時点の情報です"),"点をクリックすると情報が出ます","詳細は<a href='http://www.jma.go.jp/jp/typh/'>気象庁の台風情報</a>を確認ください。",sep = "<br/>")

m <- leaflet::leaflet(res) %>% 
  addTiles() %>% 
  setView(lng = 139.0000, lat = 35.0000, zoom = 4) %>% 
  addPolylines(~lon,~lat, color="#039", weight = 3) %>% 
  addCircleMarkers(
    radius = 10, weight = 0,
    fillOpacity = ~(1000-hpa)/(1000-930), color = ~colpal(hpa),
    stroke = TRUE, popup = ~popup
  ) %>% 
  addPopups(lng = 125.0000, lat = 42.0000, text,
            option = popupOptions(closeButton = FALSE)) %>% 
  addLegend(position = "topright", pal = colpal, values = ~hpa,
            title = "中心気圧", labFormat=labelFormat(suffix = " hpa"))
## Assuming 'lon' and 'lat' are longitude and latitude, respectively
m

参考