Rで解析：GoogleMapやStamenMapを利用「ggmap」パッケージ

地図をプロットするパッケージは多くありますが、シェープファイルを用意したり、使用方法が複雑だったりすることがあります。APIの利用回数の制限がありますが、GoogleMapやStamenMapのデータを利用して地図をプロットする「ggmap」パッケージを紹介します。なお、通常の利用であればAPIの利用回数の制限は問題ないかと思います。

本パッケージでは、Google APIを利用したルートを地図上に表示することもできます。

例に、札幌駅から北海道大学へ行くルートを示します。札幌駅から北海道大学までは以外と近く、1.127kmの距離でした。

Stamen Map
http://maps.stamen.com

パッケージのバージョンは2.6.1。R version 3.2.1でコマンドを確認しています。

パッケージのインストール

下記コマンドを実行してください。

#パッケージのインストール
install.packages("ggmap")

実行コマンド

詳細はコマンド、パッケージヘルプを確認してください。

#パッケージの読み込み
library("ggmap")

#地名で緯度経度をグーグルから取得:geocodeコマンド
LonLatData <- geocode("sapporo odori", source = "google") lon lat 1 141.3477 43.0599 #緯度経度データから住所をグーグルから取得:revgeocodeコマンド revgeocode(c(lon = LonLatData[1, 1], lat = LonLatData[1, 2])) [1] "Nish 7chome Dori, Ōdōrinishi, Chūō-ku, Sapporo-shi, Hokkaidō 060-0042, Japan" #グーグルマップ画像の取得:get_googlemapコマンド #地図の種類を指定:maptypeオプション;"terrain", "satellite", "roadmap", "hybrid"が指定可能 #拡大の指定:zoomオプション;1-21の整数で指定可能 GMapData <- get_googlemap(center = c(lon = LonLatData[1, 1], lat = LonLatData[1, 2]), zoom = 16, size = c(640, 640), scale = 2, format = "png8", maptype = "hybrid", language = "jpn", sensor = FALSE, messaging = FALSE, urlonly = FALSE, filename = "ggmapTemp", color = "color", force = FALSE) #プロット ggmap(GMapData) #fromからtoまでのルートを取得:routeコマンド #取得ルートのモード指定:modeオプション;"driving", "walking", "bicycling", "transit" RouteData <- route(from = "sapporo station", to = "hokkaido university", mode = "walking", structure = "route", output = "simple") RouteData m km miles seconds minutes hours leg lon lat 1 104 0.104 0.0646256 74 1.2333333 0.020555556 1 141.3508 43.06866 2 12 0.012 0.0074568 11 0.1833333 0.003055556 2 141.3499 43.06888 3 108 0.108 0.0671112 79 1.3166667 0.021944444 3 141.3498 43.06891 4 71 0.071 0.0441194 51 0.8500000 0.014166667 4 141.3485 43.06879 5 58 0.058 0.0360412 42 0.7000000 0.011666667 5 141.3484 43.06942 6 68 0.068 0.0422552 48 0.8000000 0.013333333 6 141.3482 43.06993 7 57 0.057 0.0354198 40 0.6666667 0.011111111 7 141.3480 43.07052 8 139 0.139 0.0863746 100 1.6666667 0.027777778 8 141.3479 43.07103 9 1127 1.127 0.7003178 811 13.5166667 0.225277778 9 141.3464 43.07100 10 NA NA NA NA NA NA NA 141.3401 43.07805 #データをプロット qmap("sapporo station, hokkaido university", zoom = 15) + geom_path(data = RouteData, aes(x = lon, y = lat), colour = "red", size = 1.5, lineend = "round") #Stamen Mapからデータを取得:Stamen Mapコマンド #取得エリアの設定 #大体の札幌 bbox <- c(left = 141.2, bottom = 43.0, right = 141.5, top = 43.15) #取得画像枚数が多くなるzoom設定だとプロットうまくいかない場合があります #maptype = "watercolor" ggmap(get_stamenmap(bbox, maptype = "watercolor", zoom = 13)) #maptype = "toner-lite" ggmap(get_stamenmap(bbox, maptype = "toner-lite", zoom = 11)) [/code]