日本地図を描く用事があったので昔やったようにもう一度やろうと思った。
だがしかし、Global Administrative Areaから地図データをダウンロードしてきたら、ファイル名の仕様が変わっていてちょっとてこずった。
library(spsurvey) jpn_GA <- vector("list", 3) for(i in seq(jpn_GA)){ jpn_GA[[i]] <- read.shape(paste("JPN_adm/JPN_adm", i-1, ".shp", sep="")) } par(mfcol=c(1, 3)) xlim<- c(128, 146) # 経度 ylim<- c(30, 44) # 緯度 for(i in seq(jpn_GA)){ par(mar=c(0, 0, 4, 0)) plot(jpn_GA[[i]], xlim=xlim, ylim=ylim, col=3) title(paste("JPN_adm", i-1, sep="")) }
都道府県区切りになっているのはJPN_adm1のようだ(図のタイトルを付け間違えた)。
いままで行ったことのある都道府県を描いてみる。
#行ったことのある都道府県 trip <- mapply(function(x) which(jpn_GA[[2]]$NAME_1 == x), c("Kyoto", "Chiba", "Tokyo", "Hokkaido", "Saitama", "Kanagawa", "Aichi", "Shiga", "Nara", "Osaka", "Wakayama", "Hyogo", "Okayama", "Hiroshima", "Kagawa", "Fukuoka", "Miyazaki", "Kumamoto", "Oita", "Okinawa")) cols <- replace(rep(0, length(jpn_GA[[2]]$NAME_1)), trip, 3) #色付け par(mar=c(0, 0, 0, 0)) plot(jpn_GA[[2]], xlim=xlim, ylim=ylim, col=cols)