Created
October 31, 2016 15:48
-
-
Save masayukeeeee/025531edb17547aa2803766ad18887a4 to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
# パッケージの読み込み | |
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