R: 北海道胆振東部地震による山腹崩壊地を地図に表示する

タイトルのとおりです。データは北海道水産林務部森林計画課森林計画関係資料ダウンロードページから入手できます(ランセンスはCC-BYです)。

Rコードです。geom_sfのおかげであっさりとプロットできます。背景をOpenStreetMapにしていますが、「メモ：ggplot2::geom_sf()の背景にOpenStreetMapの地図画像を使いたいときはggspatial::annotation_map_tile()」を参考にしました。

また、sum(st_area(houkai))で崩壊地面積の合計が計算できます。

> sum(st_area(houkai))
43819262 m^2

タグ：R ggplot2