Create a gist now

Instantly share code, notes, and snippets.

What would you like to do?
# パッケージの読み込み
library(rgl)
library(maptools)
library(spsurvey)
# ダウンロードしたzipファイルを解答してできたフォルダがあるディレクトリに行きましょう。(下のコードはデスクトップに置いてある場合です)
setwd('~/Desktop/ne_10m_admin_1_states_provinces_scale_rank')
# shpファイルの読み込み
shape.world <- read.shape(filename='ne_10m_admin_1_states_provinces_scale_rank.shp')
# 日本の座標が保存されている要素番号を抜き出す
jpn.num <- which(shape.world$sr_sov_a3=='JPN')
# 座標データを抜き出す
tmp.data <- shape.world@polygons
# 日本の座標データを格納するリストを作成
map.data <- list()
# 要素数
n <- length(jpn.num)
# 座標だけ抜き出してmap.dataに格納していく
for(i in 1:n){
num <- jpn.num[i]
tmp1 <- tmp.data[[num]]
tmp2 <- tmp1@Polygons
tmp3 <- tmp2[[1]]
tmp4 <- tmp3@coords
map.data[[i]] <- tmp4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment