方位記号やスケールバーをplotコマンドで出力した図へ追記できるパッケージの紹介です。rasterパッケージとの組み合わせも面白いかもしれません。
方位記号とスケールバーを描写する「addnortharrow」と「addscalebar」、地名から緯度経度を取得する「searchbbox」コマンドを紹介します。
パッケージバージョンは0.2.3。windows 11のR version 4.1.3で確認しています。
パッケージのインストール
下記、コマンドを実行してください。
#パッケージのインストール install.packages("prettymapr")
実行コマンド
詳細はコメント、パッケージのヘルプを確認してください。
#パッケージの読み込み library("prettymapr") ###方位記号の描写:addnortharrowコマンド##### #描写の準備 plot(0:1, 0:1, type = "n", axes = FALSE, xlab = "", ylab = "") box() #描写位置の指定 Position <- c("bottomleft", "bottomright", "topleft", "topright") #色の指定 Cols <- c("#e1e6ea", "#505457", "#4b61ba", "#a87963") #矢印の描写 for(n in 1:4){ addnortharrow(pos = Position[n], padin = c(0.15, 0.15), scale = n, lwd = 1, border = "black", cols = c("white", Cols[n]), text.col = "black") } ######## ###スケールバーの描写:addscalebarコマンド##### #描写の準備 plot(0:1, 0:1, type = "n", axes = FALSE, xlab = "", ylab = "") box() #描写位置の指定 Position <- c("bottomleft", "bottomright", "topleft", "topright") #色の指定 Cols <- c("#e1e6ea", "#505457", "#4b61ba", "#a87963") #スタイルの指定 Styles <- c("bar","ticks", "bar","ticks") #緯度経度システムの設定 #lat/lon,Spherical Mercator otherwiseの順 LatLon <- c(4326, 3857, 4326, 3857) #矢印の描写 for(n in 1:4){ addscalebar(plotunit = "cm", plotepsg = LatLon[n], widthhint = 0.25, unitcategory = "metric", htin = n, padin = c(0.15, 0.15), style = "bar", bar.cols = c(Cols[n], "white"), lwd = 1, linecol = "black", tick.cex = 0.7, labelpadin = 0.08, label.cex = 0.8, label.col = "black", pos = Position[n]) } ######## #緯度経度の取得:searchbboxコマンド #PickPoint.io APIを利用 GetArea <- searchbbox("KYOTO") #確認 GetArea # min max #x 135.55901 135.87844 #y 34.87492 35.32122
出力例
・addnortharrowコマンド
・addscalebarコマンド
少しでも、あなたの解析が楽になりますように!!