読者です 読者をやめる 読者になる 読者になる

疫学と医療統計学と遺伝学と時々、大学院生活

疫学を専門とする大学院生の研究に関する備忘録的ページ。

R package 'ggmap'を用いたマッピング

今回は'ggmap'と呼ばれるパッケージを使用して、愛知県にあるナゴヤドームへの経路(最寄駅(大曽根駅)から)を図示することにする。

1. 早速、パッケージのインストールと読み出す

install.packages("ggmap")
library(ggmap)

2. 地名で緯度経度をグーグルから取得する。

LonLatData <- geocode("ナゴヤドーム")

3. 緯度経度データから住所をグーグルから取得する。

revgeocode(c(lon = LonLatData[1, 1], lat = LonLatData[1, 2]))

4. グーグルマップ画像を取得する。

地図の種類を指定:maptypeオプション;"terrain", "satellite", "roadmap", "hybrid"が指定可能
拡大の指定:zoomオプション;1-21の整数で指定可能

GMapData <- get_googlemap(center = c(lon = LonLatData[1, 1], lat = LonLatData[1, 2]),
                          zoom = 15, size = c(640, 640), scale = 2, format = "png8",
                          maptype = "hybrid", language = "jpn", sensor = FALSE,
                          messaging = FALSE, urlonly = FALSE, filename = "ggmapTemp",
                          color = "color", force = FALSE)

5. プロットする

ggmap(GMapData)

f:id:ryosukefujii0320:20160218191734p:plain

6. fromからtoまでのルートを取得する

取得ルートのモード指定:modeオプション;"driving", "walking", "bicycling", "transit"

RouteData <- route(from = "大曽根駅", to = "ナゴヤドーム", mode = "walking",
                   structure = "route", output = "simple")
RouteData
##      m    km     miles seconds    minutes       hours leg      lon
## 1   62 0.062 0.0385268      45 0.75000000 0.012500000   1 136.9368
## 2   95 0.095 0.0590330      75 1.25000000 0.020833333   2 136.9370
## 3  518 0.518 0.3218852     375 6.25000000 0.104166667   3 136.9378
## 4   26 0.026 0.0161564      19 0.31666667 0.005277778   4 136.9433
## 5   82 0.082 0.0509548      59 0.98333333 0.016388889   5 136.9435
## 6   83 0.083 0.0515762      59 0.98333333 0.016388889   6 136.9436
## 7    8 0.008 0.0049712       6 0.10000000 0.001666667   7 136.9445
## 8   37 0.037 0.0229918      27 0.45000000 0.007500000   8 136.9445
## 9   98 0.098 0.0608972      71 1.18333333 0.019722222   9 136.9445
## 10  86 0.086 0.0534404      62 1.03333333 0.017222222  10 136.9455
## 11 123 0.123 0.0764322      88 1.46666667 0.024444444  11 136.9463
## 12   7 0.007 0.0043498       5 0.08333333 0.001388889  12 136.9461
## 13 123 0.123 0.0764322      89 1.48333333 0.024722222  13 136.9461
## 14  41 0.041 0.0254774      30 0.50000000 0.008333333  14 136.9460
## 15 352 0.352 0.2187328     263 4.38333333 0.073055556  15 136.9465
## 16  NA    NA        NA      NA         NA          NA  NA 136.9483
##         lat
## 1  35.19149
## 2  35.19200
## 3  35.19209
## 4  35.19099
## 5  35.19078
## 6  35.19024
## 7  35.19006
## 8  35.18999
## 9  35.18970
## 10 35.18949
## 11 35.18924
## 12 35.18816
## 13 35.18809
## 14 35.18714
## 15 35.18704
## 16 35.18528

7. データをプロットする

qmap("大曽根駅, ナゴヤドーム", zoom = 15) +
geom_path(data = RouteData, aes(x = lon, y = lat),
          colour = "red", size = 1.5, lineend = "round")


f:id:ryosukefujii0320:20160218191715p:plain

詳細はパッケージの説明をご参照ください。

20160218
RF