Rで解析：お手軽！インタラクティブに地図でデータを確認できる「leafletR」パッケージ

地図にデータをプロットし、インタラクティブに数値を確認できる資料を作成できる「leafletR」パッケージをご紹介します。面白いのでぜひ使用してみてください。

パッケージバージョンは0.4-1。windows11のR version 4.2.2で確認しています。

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

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

if(!require("leafletR", quietly = TRUE)){
  devtools::install_github("chgrl/leafletR")
}

実行コマンド

「フィジーの地震」とオリジナル「東京駅の位置」をプロットするコードを紹介します。ベースとなる地図はbase.mapで”osm”, “tls”, “mqosm”, “mqsat”, “water”, “toner”, “tonerbg” , “tonerlite”, “positron”, “darkmatter”が選べます。

「フィジーの地震」プロットコード

#&#12521;&#12452;&#12502;&#12521;&#12522;&#12398;&#35501;&#12415;&#36796;&#12415;
library("leafletR")

#&#12487;&#12540;&#12479;&#12398;&#35501;&#12415;&#36796;&#12415;
data(quakes)

#&#12487;&#12540;&#12479;&#12398;&#20316;&#25104;
q.dat <- toGeoJSON(data = quakes[1:99,], dest = tempdir(), name = "quakes")

#&#12503;&#12525;&#12483;&#12488;&#12473;&#12479;&#12452;&#12523;&#12398;&#35373;&#23450;
q.style <- styleGrad(prop = "mag", breaks = seq(4, 6.5, by=0.5), 
                        style.val = rev(heat.colors(5)), leg = "Richter Magnitude", 
                        fill.alpha = 0.7, rad = 8)

#&#22320;&#22259;&#12398;&#20316;&#35069;
q.map <- leaflet(data = q.dat, dest = tempdir(), title = "Fiji Earthquakes", 
                    base.map = "osm", style = q.style, popup = "mag", incl.data = TRUE)

#&#12502;&#12521;&#12454;&#12470;&#12391;&#34920;&#31034;
browseURL(q.map)

「東京駅の位置」プロットコード

#&#12521;&#12452;&#12502;&#12521;&#12522;&#12398;&#35501;&#12415;&#36796;&#12415;
library("leafletR")

#&#26481;&#20140;&#39365;&#12398;&#20301;&#32622;&#12487;&#12540;&#12479;&#12434;&#20316;&#25104;
LocateData <- as.data.frame(cbind(lat = 35.681283, long = 139.766092, landname = "Tokyo Station"))
q.dat <- toGeoJSON(data = LocateData, dest = tempdir(), name = "Tokyo Station")

#&#12503;&#12525;&#12483;&#12488;&#12473;&#12479;&#12452;&#12523;&#12398;&#35373;&#23450;
sty <- styleSingle(col = "red", fill = "red")

#&#22320;&#22259;&#12398;&#20316;&#35069;
q.map <- leaflet(data = q.dat, dest = tempdir(), title = "Tokyo Station", 
                    base.map = "osm", style = sty, popup = "landname", incl.data=TRUE)

#&#12502;&#12521;&#12454;&#12470;&#12391;&#34920;&#31034;
browseURL(q.map)

出力例

・フィジーの地震

・東京駅の位置

少しでも、あなたの解析が楽になりますように！！