Rで描いた地図にスケールバーを追加する方法

この記事では、Rで描いた地図へ、スケールバーを追加する方法について説明します。

スケールバー

記事の概要です。

  1. スケールバーを追加するためにggsnパッケージをインストールします。
  2. 茨城県の地図を作成します。
  3. ggsnパッケージのscalebar関数を使い、スケールバーを追加します。

スケールバーを追加する部分です。

## スケールーバーを追加して描画
library(ggsn)
ibaraki + scalebar(shpCsv,           # シェープファイル
                   dist      = 15,   # 1メモリが表す距離
                   dist_unit = "km", # 距離の単位
                   transform = TRUE, # 地図の座標が緯度・経度の場合TRUE
                   model     = "GRS80", # 楕円体モデル。シェープファイルに依存。
                   location  = "bottomleft") # "topright", "bottomright", "bottomleft", "topleft"

ソースコード全体です。

## シェープファイルとCSVファイルを結合

library(sf)
shp <- st_read("./data/N03-20150101_08_GML", options = "ENCODING=SHIFT-JIS")

d <- read.csv("./data/kekkahyo1_ibaraki_2015-2045.csv")

library(dplyr)
shpCsv <- left_join(shp, d, by = c("N03_004" = "市区町村"))

## 2020年の推計人口を、3グループ(増加、横ばい、減少)に分ける

x <- shpCsv$cy2020
xGrp <- case_when(x >= 102 ~ "増加",
                  x >= 98 & x < 102 ~ "横ばい",
                  x < 98 ~ "減少",
                  TRUE ~ NA_character_) # グループ化
xGrp <- factor(xGrp, levels = c("増加", "横ばい", "減少")) # 順序(≠大小関係)を設定
shpCsv$cy2020Grp <- xGrp # cy2020Grpカラムに格納

## 地図を ibaraki に代入

library(ggplot2)
ibaraki <- ggplot(shpCsv) + geom_sf(aes(fill = cy2020Grp)) +
    scale_fill_manual("2020年の推計人口", values = c("#fb9a99", "#ffffff", "#a6cee3")) +
    theme(legend.position = c(0.2, 0.8)) # 左下は 0, 0、右上は 1, 1

ibaraki                                  # 一旦、スケールバー無しで描画

## スケールーバーを追加して描画
library(ggsn)
ibaraki + scalebar(shpCsv,           # シェープファイル
                   dist      = 15,   # 1メモリが表す距離
                   dist_unit = "km", # 距離の単位
                   transform = TRUE, # 地図の座標が緯度・経度の場合TRUE
                   model     = "GRS80", # 楕円体モデル。シェープファイルに依存。
                   location  = "bottomleft") # "topright", "bottomright", "bottomleft", "topleft"


ミサキ
ミサキ

Rの地図にも、スケールバーは表示できます!
私は知りませんでした(- -;
また、この記事は、QGISのこちらの記事と対応しています。


この記事で使用したRとパッケージのバージョンを記載しておきます。
バージョンが異なると、動作や表示も若干異なる場合があるからです。

  • R 3.6.1
  • sf 0.7-6
  • dplyr 0.8.3
  • ggplot 3.2.0
  • ggsn 0.5.0



準備:ggsnパッケージのインストール

ggsnパッケージをインストールします。

このパッケージは、ggplot2などで描いた地図にスケールバー方位記号を追加するパッケージです。

方位記号を追加するときにも使いますので、ぜひインストールを!

インストールの仕方は、こちらを参考にして下さい。

準備:茨城県の地図を作成

スケールバーを追加する地図を作成します。

作成する地図は、こちらの記事で作成した地図をベースにし、凡例のタイトルと位置を変更したものです。
以下にソースコードを示します。

記事とこのソースコードを基に、地図をibarakiに代入して下さい。

## シェープファイルとCSVファイルを結合

library(sf)
shp <- st_read("./data/N03-20150101_08_GML", options = "ENCODING=SHIFT-JIS")

d <- read.csv("./data/kekkahyo1_ibaraki_2015-2045.csv")

library(dplyr)
shpCsv <- left_join(shp, d, by = c("N03_004" = "市区町村"))

## 2020年の推計人口を、3グループ(増加、横ばい、減少)に分ける

x <- shpCsv$cy2020
xGrp <- case_when(x >= 102 ~ "増加",
                  x >= 98 & x < 102 ~ "横ばい",
                  x < 98 ~ "減少",
                  TRUE ~ NA_character_) # グループ化
xGrp <- factor(xGrp, levels = c("増加", "横ばい", "減少")) # 順序(≠大小関係)を設定
shpCsv$cy2020Grp <- xGrp # cy2020Grpカラムに格納

## 地図を ibaraki に代入

library(ggplot2)
ibaraki <- ggplot(shpCsv) + geom_sf(aes(fill = cy2020Grp)) +
    scale_fill_manual("2020年の推計人口", values = c("#fb9a99", "#ffffff", "#a6cee3")) +
    theme(legend.position = c(0.2, 0.8)) # 左下は 0, 0、右上は 1, 1

一旦、スケールバーなしで描画してみましょう。

ibaraki                                  # 一旦、スケールバー無しで描画
茨城県の地図(スケールバーなし)

スケールバーの追加

では、スケールバーを追加します。

スケールバーの追加には、ggsnパッケージのscalebar関数を使います。

## スケールーバーを追加
library(ggsn)
ibaraki + scalebar(shpCsv,           # シェープファイル
                   dist      = 15,   # 1メモリが表す距離
                   dist_unit = "km", # 距離の単位
                   transform = TRUE, # 地図の座標が緯度・経度の場合TRUE
                   model     = "GRS80", # 楕円体モデル。シェープファイルに依存。
                   location  = "bottomleft") # "topright", "bottomright", "bottomleft", "topleft"

各引数の意味は、コメントから理解いただけると思います。
model引数はシェープファイルの仕様を確認する必要があります。

以下の地図が表示されれば、完了です。

茨城県の地図(スケールバーあり)

おわりに

お疲れさまでした。
地図さえ描けていれば、スケールバーの追加は簡単ですね!
次回は方位記号を追加します。
では。


YouTube: ミサキさん

問い合わせ先
Mail: caprico.aries@gmail.com
Twitter: https://twitter.com/caprico_aries(無言フォローもお気軽に)

コメント

  1. […] Rで描いた地図にスケールバーを追加する方法 […]

  2. […] 凡例の調整法 | R, GISRで描いた地図にスケールバーを追加する方法Rで描いた地図に方位記号を追加する方法 […]

タイトルとURLをコピーしました