So-net無料ブログ作成

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

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

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

Rplot.png

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

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

タグ:R ggplot2
nice!(3)  コメント(0) 
共通テーマ:日記・雑感

nice! 3

コメント 0

コメントを書く

お名前:
URL:
コメント:
画像認証:
下の画像に表示されている文字を入力してください。

Facebook コメント