#はじめに
この記事は、FOSS4G 2015年アドベントカレンダーの前準備用の記事です。
面積按分をするために、e-Statから3次メッシュの人口統計情報をダウンロードし、PostgreSQLにデータを登録するまでの手順について説明します。

##必要なソフトウェア
今回用いるソフトは以下の2つです。

QGIS：最新版はこちらからダウンロードできます。
PostgreSQL：PostGISがSQL1文で導入できる比較的新しいバージョン9.1以上がおすすめです。ダウンロードはこちら

WindowsOSなら、どちらのソフトもGUIのインストーラーがあるので導入はそれほど難しくないと思います。また、PostGISはPostgreSQLのインストール後に別途インストール作業が必要です。こちらで手順をご確認下さい。

##データの入手
e-StatのWebサイトから「統計情報」を「メッシュ形状」をダウンロードします。
###ダウンロード

ダウンロードサイト「e-Stat 政府統計の総合窓口」
e-Stat → 地図で見る統計（統計GIS）→ データダウンロード
Step1: 統計調査（集計）を選択データ
調査名：国勢調査
データ：平成２２年国勢調査（国勢調査ー世界測地系１ｋｍメッシュ）
Step2：統計表を選択（複数選択可能）
「男女別人口総数及び世帯総数」を選択
「統計表各種データダウンロードへ」ボタンをクリック
Step3：地域選択
市区町村：「097：M5339」たくさんのメッシュがありますが、都心の5339を選択します。
「検索」ボタンを押下
Step4：データダウンロード
統計データのダウンロード（「M5339」と書かれたリンクをクリックします）
境界データのダウンロード（しばらくすると境界データのリストが表示されます。「世界測地系緯度経度・Shape形式をダウンロード」）

QGISで読込

ダウンロードしたデータをQGISで読み込みます。
データは統計データと境界データの２つに分かれているので、QGISで１つに結合し、結合したデータをPostgreSQLに投入します。

ファイルの集類	ファイル名称
統計データ	tblT000608S53390.zip
境界データ	S002005112009DDSWS05339.zip

データの加工

統計データはtblT000608S53390.zipを解凍するとテキストファイルが１つ得られます。
先頭行に英語の列名がが、2行目に日本語のが入っています。

QGISのCSVデータの読込機能では、複数行のカラム名に対応していないので、1行目のKEY_CODE以降を削除して上書き保存します。

境界データで読込

境界データを解凍して、shapeデータ「MESH05339.shp」をダブルクリックします。

QGISで表示されます

統計データの読込

デリミテッドテキストレイヤの読込
QGISでレイヤ → レイヤの追加 → デリミテッドテキストレイヤの追加　を選択します。

加工済みの統計データの選択
先頭行を編集したtblT000608S53390.txtを読み込みます。

- ファイル形式：CSV（コマンドで区切られた値） - レコードオプション:最初のレコードはフィールド名を保持している［チェック］ - ジオメトリ定義：ジオメトリなし（属性のみのテーブル）［チェック］上手く行くと、QGISに属性のみのテーブルが読み込まれます

境界データと統計データの結合

境界データのプロバティを表示
レイヤウインドウの境界データ「MESH05339」をダブルクリックして、プロバティを表示

データの結合
結合タブをクリック → 右下の「＋」ボタンをクリック
統計データと境界データの2つのデータを結合する列名を指定

- 結合するレイヤ：統計データの名称 - 結合フィールド：KEY_CODE - ターゲットフィールド：KEY_CODE 上記の設定で、統計情報のメッシュコードと境界データのメッシュコードが同一の行が結合し、１つのテーブルのように扱えます

shape形式で保存

レイヤを名前をつけて保存
レイヤウインドウの境界データ「MESH05339」を右クリックして、「名前をつけて保存」を選択します。
名称は特に拘る必要はありませんが、今回は「mesh5339.shp」という名称にしました。

形式を指定する
ファイルの形式を「ESRI Shapefile」と指定して下さい。
ファイル名を指定
「ブラウズ」ボタンをクリックし、ファイルの保存場所とファイル名を指定して下さい。
保存
「OK」ボタンをクリックして保存します。

#データベースへ投入
QGISの拡張機能である「SPIT」を使って、shapeデータをPostGISに投入します。
SPITがインストールされていれば、象さんのアイコンが表示されています。

作成されたテーブル
SPITが自動的に作成したテーブルを観察してみます。
テーブル構造

CREATETABLE

CREATE TABLE mesh5339
(
  gid serial NOT NULL, --ID
  "MESH1_ID" character varying(4), --1次メッシュコード
  "MESH2_ID" character varying(2), --2次メッシュコード
  "MESH3_ID" character varying(2), --3次メッシュコード
  "OBJ_ID" integer, --メッシュ形状のID
  "KEY_CODE" character varying(8), --メッシュコード
  "tblT000608" integer, --人口
  "tblT0006_1" integer, --男性人口
  "tblT0006_2" integer, --女性人口
  "tblT0006_3" character varying(254), --世帯数
  the_geom geometry(Polygon,4612),--メッシュ形状：測地系（新日本測地系）
  CONSTRAINT mesh5339_pkey PRIMARY KEY (gid)
)
WITH (
  OIDS=FALSE
);

データ例

お疲れ様でした。ここまで来れば準備OKです！