Skip to content

Instantly share code, notes, and snippets.

@masayukeeeee
Created October 31, 2016 15:48
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save masayukeeeee/025531edb17547aa2803766ad18887a4 to your computer and use it in GitHub Desktop.
Save masayukeeeee/025531edb17547aa2803766ad18887a4 to your computer and use it in GitHub Desktop.
# パッケージの読み込み
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