日本地図をプロットしてみよう

日本地図のデータは、こちらのウェブサイト の指示に従って、ダウンロードし、所定のディレクトリに保存しておきます。

都道府県の行政界をプロットするには、次のようなスクリプトを実行します。

library(maptools)
library(ggplot2)

japan <- readShapePoly("./map_tutor/japan_pref/jpn_pref.shp")

map1 <- ggplot(japan, aes(long, lat, group = group))
map1 <- map1 + geom_polygon(colour = "gray20", size = 0.2, fill = "white")
map1 <- map1 + xlim(c(120, 150))
map1 <- map1 + coord_equal()

画面に地図を表示するには次のようにします。

print(map1)

次に千葉県の部分を赤色で塗ってみましょう。

flag <- rep(FALSE, length(japan))
flag[japan$PREF == "Chiba"] <- TRUE

df <- fortify(spCbind(japan, flag), region = "flag")

map2 <- ggplot(df, aes(long, lat, group = group))
map2 <- map2 + geom_polygon(aes(fill = id))
map2 <- map2 + scale_fill_manual(values = c("white", "red"))
map2 <- map2 + geom_path(data = japan, aes(long, lat, group = group), colour = "gray20", 
    size = I(0.1))
map2 <- map2 + coord_equal() + labs(x = "", y = "") + theme(axis.ticks = element_blank(), 
    axis.text = element_blank()) + theme(legend.position = "none") + xlim(c(120, 
    150))

画面に地図を表示するには次のようにします。

print(map2)