【R】国勢調査 小地域 Shapefile を使って緯度経度から住所を求める(市までわかっているときver.)
はじめに
@u_riboさんによるこんなツイートがある。
一部で話題の(?)e-Statが提供する国勢調査小地域のShapefileをRからダウンロードする関数を書きました (gifが高速過ぎる orz) コードは https://t.co/Dd2R58ov94 pic.twitter.com/SRo2H1Fk4W
— Uryu Shinya (@u_ribo) 2017年7月14日
要するに、shapefileの細かい版ですね。細かいので、市レベルでファイルが分かれているようです。
今回は、これを使って逆ジオコーディングをやってみます。ただし、「やってみた」レベルなので今回は知りたい緯度経度が何市にあるかまでは絞れているとします。 市を見つけるのには、RであればRで緯度経度から都道府県・市区町村を求める、SASであればSASによる逆ジオコーディング(緯度経度から都道府県・市を求める)が参考になります。
続きを読むSASで緯度経度から地域メッシュコードを計算する
はじめに
検索してもでてこなかったので、緯度経度からその点が含まれる3次メッシュ(基準地域メッシュ)コードを算出する関数をSASで作成しました。(僕が知らないだけでもしかしたら既にあるのかもしれませんが。)
地域メッシュコードの概要については統計局ホームページ/地域メッシュ統計についてに、その算出方法については
緯度経度からメッシュコードを作成する方法に解説が載っていますが、標準地域メッシュコーどは約1km四方で、大体の粒度は地図上で標準地域メッシュを確認するページで確認できます。ちなみに、JR田町駅周辺だとこのくらいの粒度です。
SASによる逆ジオコーディング(緯度経度から都道府県・市を求める)
SASによるジオコーディングはアメリカやカナダのデータは標準装備されていて非常にやりやすいのですが、日本のデータでやろうとするとshapfileから用意しなければならず、非常にやりづらいです。一応、動くものが作れたので備忘録としてメモ。
shapefile
シェープファイル (Shapefile) は、 他の地理情報システム(GIS)間でのデータの相互運用におけるオープン標準として用いられるファイル形式である。 例えば、井戸、川、湖などの空間要素がベクタ画像である点 (数学)、線分、多角形で示され、各要素に固有名称や温度などの任意の属性を付与できる。(Wikipediaより)
今回は、国土数値情報 ダウンロードサービスの「2014年版 行政区域データ 全国のファイル」を使います。 ダウンロード及びzipの解凍は手動でやってもいいのですが、今回はRで緯度経度から都道府県・市区町村を求める - Qiitaを参考にしてRを使いました。
続きを読む