SASによる逆ジオコーディング(緯度経度から都道府県・市を求める)

SASによるジオコーディングはアメリカやカナダのデータは標準装備されていて非常にやりやすいのですが、日本のデータでやろうとするとshapfileから用意しなければならず、非常にやりづらいです。一応、動くものが作れたので備忘録としてメモ。

shapefile

シェープファイル (Shapefile) は、他の地理情報システム(GIS)間でのデータの相互運用におけるオープン標準として用いられるファイル形式である。例えば、井戸、川、湖などの空間要素がベクタ画像である点 (数学)、線分、多角形で示され、各要素に固有名称や温度などの任意の属性を付与できる。(Wikipediaより)

今回は、国土数値情報　ダウンロードサービスの「2014年版行政区域データ全国のファイル」を使います。ダウンロード及びzipの解凍は手動でやってもいいのですが、今回はRで緯度経度から都道府県・市区町村を求める - Qiitaを参考にしてRを使いました。

library(magrittr)
library(dplyr)
library(kokudosuuchi)

## このサービスは、「国土交通省　国土数値情報（カテゴリ名）」をもとに加工者が作成
## 以下の国土数値情報ダウンロードサービスの利用約款をご確認の上ご利用ください：
## 
## http://nlftp.mlit.go.jp/ksj/other/yakkan.html

# 2014年版 行政区域データ 全国のファイルをダウンロードして解凍
dl.url <- kokudosuuchi::getKSJURL(identifier = "N03") %>%
  dplyr::filter(year == 2014, areaCode == 0) %>% use_series(zipFileUrl)

download.file(dl.url,
              destfile = "N03-150101_GML.zip")
unzip("N03-150101_GML.zip")