maptoolsで日本の地図を描く
readShapePoly関数を使うため、maptoolsをlibraryで読み込みます。
library(maptools)
3種類の日本地図を描くため、GADM(https://gadm.org/index.html)から日本のデータ(shapefile)をダウンロード・展開し、使用しているdirectoryにデータを保存しておきます。
fortify (ggplot2の関数)でreadShapePolyで読み込んだshapefileを一般的なデータフレームに変換します。
Japan_shp0 <- readShapePoly("gadm36_JPN_shp/gadm36_JPN_0.shp") Japan_map0 <- fortify(Japan_shp0) ggplot(Japan_map0, aes(x = long, y =lat, group = group)) + geom_path()
参考までにJapan_map0の内容を見てみます。
> str(Japan_map0) 'data.frame': 562995 obs. of 7 variables: $ long : num 138 138 138 138 138 ... $ lat : num 34.6 34.6 34.6 34.6 34.6 ... $ order: int 1 2 3 4 5 6 7 8 9 10 ... $ hole : logi FALSE FALSE FALSE FALSE FALSE FALSE ... $ piece: Factor w/ 1398 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ... $ id : chr "0" "0" "0" "0" ... $ group: Factor w/ 1398 levels "0.1","0.2","0.3",..: 1 1 1 1 1 1 1 1 1 1 ... >
Japan_shp1 <- readShapePoly("gadm36_JPN_shp/gadm36_JPN_1.shp") Japan_map1 <- fortify(Japan_shp1) ggplot(Japan_map1, aes(x = long, y =lat, group = group)) + geom_path()
Japan_shp2 <- readShapePoly("gadm36_JPN_shp/gadm36_JPN_2.shp") Japan_map2 <- fortify(Japan_shp2) ggplot(Japan_map2, aes(x = long, y =lat, group = group)) + geom_path()
