123
17． ＜参考①＞施設勢力圏域ごとの高齢者数等を集計
医療機関（診療所）の分布から施設勢力圏（ボロノイ分割圏域）を作成し、それに含まれ る町丁字の高齢者数等を集計することにより、施設勢力圏域ごとの指標算出を行います。 医療機関（診療所）データは本マニュアル１５．１で作成したものを使用します。17.1 ボロノイ分割圏域の作成
ボロノイ分割とは、複数個の点（母点）に対して、同一距離空間上の、他の点がどの母点 に近いかによって、領域分けされた分割のことです。利用者が最も近くの施設に行くとした 場合の施設勢力圏域を示す図となり、施設最適配置(マクドナルドの Mac GIS)や商圏分析(コ ンビニ)などに使われています。 (1) ベクタ(t)→ジオメトリツール(G)→ ボロノイポリゴン をクリック (2) ボロノイポリゴン ウィンドウが開く ① 「ポイントベクタレイヤの入力」を「診療所」に選択 ② 「ブラウズ」から、「施設勢力圏.shp」として「C:\GISjinzai\」に保存する ① ② ③ 施設勢力圏を作成 する元レイヤ124 (3) 「ジオメトリ」ウィンドウの質問には「Yes」を選択 (4) エラーが発生し、レイヤが追加されなかった場合は、 レイヤ→ベクタレイヤの追加 から、作成された「施設勢力圏.shp」を読み込む ↓エラー画面の例 無事に施設勢力圏.shp は作成されますが、自動的に読み込んでレイヤを追加する動作に失 敗しています。手動で改めて作成したレイヤを追加する必要があります。 次回以降は、(3)での質問に「No」を選択し、手動で作成されたレイヤを読み込むように するとエラーが発生しなくなります。 (5) 施設勢力圏のボロノイポリゴンが作成されます 施設勢力圏の shp が 作成される。
125
17.2 施設勢力圏の表示の変更
前図では、ポリゴンがすべて塗りつぶされて町丁・字などの情報と重ね合わせてみること ができないため、施設勢力圏ポリゴンの表示方法を変更し、重ね合わせて表示させるように します。 (1) レイヤ施設勢力圏をダブルクリック もしくは右クリック→プロパティ (2) レイヤプロパティが開く ① 「スタイル」タブを開く ② 「変更」をクリック ③ 「塗りつぶしスタイル」を「ブラシ無し」に変更 ④ その他、「ボーダーカラー」や「ボーダースタイル」を適宜変更して確定する。 ① ② ③ ④ ⑤ ⑥126
(3) 施設勢力圏が境界線のみで中抜きで表示されます。
施設勢力圏域が可視 化される。
127
17.3 町丁字と施設勢力圏の空間結合
施設勢力圏域ごとに含まれる町丁字を GIS で空間的に結合させます。町丁字の属性テー ブルに、どの施設の勢力圏域に属しているか情報を付与し、それをもとに高齢者数等の集計 を行います。 (1) ベクタ(t)→データマネジメントツール(D)→場所で属性を結合する をクリック (2) 場所で属性を結合するウィンドウが開く ① 対象ベクタレイヤを「H22 稲城市」に選択 ② ベクタレイヤを結合するを「施設勢力圏」に設定 ③ 「ブラウズ」をクリックし、「temp.shp」として「C:\GISjinzai\」に保存する。 (3) 「空間結合」ウィンドウの質問には「Yes」を選択 ① ② ③ ④ 属する施設の情報が付与 された町丁字の shp が作成 される。128 (4) エラーが発生し、レイヤが追加されなかった場合は、 レイヤ→ベクタレイヤの追加 から、作成された「temp.shp」を読み込む (5) 作成したレイヤ temp を右クリック→属性テーブルを開く 町丁字ごとにそれぞれの町丁字から最も近い施設(P14_007)が選ばれている。 どの施設勢力圏に属 するか情報が付与さ れている。
129
17.4 施設圏域ごとに町丁字の色分け
町丁目ごとに付与されたどの施設に属しているかの情報をもとに色を塗りわけ、施設勢力 圏を可視化します。 (1) レイヤ temp をダブルクリック もしくは右クリックしてプロパティ(P)をクリックする。 (2) レイヤプロパティが開く ① 「スタイル」タブを開く ② 「分類された」をクリック ③ 「カラム」を「P14_002」に変更 ※施設名が格納されたフィールドです。 ④ 色諧調を選択 ⑤ 「分類」をクリックし、 ⑥ 「OK」ボタンで確定させる。 ① ② ③ ④ ⑤ ⑥130 (3) それぞれ選択された勢力圏の施設ごとに色分けがなされます。 ※上図ではレイヤの順序を並び替えてあります。 属する施設ごとに町 丁字が塗り分けられ る。
131
17.5 施設勢力圏域ごとの高齢者数の集計
QGIS で勢力圏の施設情報を空間的に結合させた町丁字データを Excel に読み込ませるこ とによって、施設勢力圏域ごとの高齢者数の集計を行います。17.5.1 QGIS からのテーブルの出力
(1) レイヤ temp を右クリック→属性テーブルを開く (2) セクションを反転する ボタンを押して、すべてのレコードを選択する (3) 選択している行をクリップボードにコピーする ボタンを押す132
17.5.2 Excel 上での集計
(1) Excel を起動して、新しいブックに貼り付ける ※以降、Excel 上での作業について、本マニュアルでは Excel2010 を使用した画面となって います。 (2) 一旦、名前を付けて保存する。 マイドキュメントフォルダに 稲城市集計 2.xls として保存する。 (3) ユーザー設定の並べ替え 最も近い施設(P14_002) （リストの下のほうにあります）で並べ替える。 ③ ① ②133 (4) 神奈川県で違う列に格納されてしまっているデータは、選択して移動させます。 (5) データ-小計をクリック ※Excel2003 では「集計」 (6) 集計の設定を行う ① グループの基準(A) を P04_002 列（最寄りの施設名称） に設定 ② 集計の方法(U) を 合計 に設定する。 ③ 集計するフィールド(D) を T000573019（65 歳以上人口）にする ① ② ③ ④
134 (7) 画面左端にできる[2]のボタンを押す 施設勢力圏域ごとの高齢者数が集計できます。 ※図では、一部列を省略して表示させています。 (8) CA 列 T000573019 (高齢者数)を選択して、右クリック→切り取り(T) 拡大図 診療所勢力圏域ごとの高齢者数 が集計される。 診療所勢力圏域ごとの高齢者数 が集計される。
135 (9) DU 列で、右クリック→貼り付け(P) (10) 体裁を整える 表頭を「診療所、65 歳以上高齢者数」に修正する、罫線を入れるなど報告書向けに体裁 を整えます。 その後、上書き保存します。 作成した表は本マニュアル１６．の手順により、報告書形式で取りまとめることができま す。
136
18． ＜参考②＞OpenStreetMap から道路を取り込む
OpenStreetMap（http://osm.jp/）は商用・非商用を問わずに誰でも自由に利用、編集できる 地図データベースです。一般のボランティアによって編集された世界中の地図が無料で利用 できるほか、QGIS のプラグイン（追加機能）を用いてダウンロードし、QGIS 上に表示す ることができます。 一般の人が編集したデータであるため、公的機関発行のデータに比べ正確性・信頼性に劣 るという面もありますが、現在無料で利用することのできない日本の「道路」のシェープデ ータを得ることができます。 本章では、OpenStreetMap からデータをダウンロードし、道路ラインのデータを抽出して 表示させます。18.1 OpenStreetMap プラグインの有効化
QGIS 上で OpenStreetMap を扱うためのプラグインを有効化します。 (1) メニューより プラグイン(P)→プラグインの管理 をクリックする。(2) OpenStreetMap plugin のチェックをつけて有効化し、OK ボタンを押してウィンドウ を閉じる。
137
18.2 OpenStreetMap データのダウンロード
プラグインを使用し、OpenStreetMap からデータをダウンロードします (1) 稲城市全体が表示されるように表示を調整する。 (2) メニューより プラグイン(P)→OpenStreetMap→OSM データのダウンロード をクリ ックする。 (3) OSM データのダウンロード画面の「…」を押し、「C:\GISjinzai\」に「inagishi.osm」 として保存する。 ※範囲は概ね現在表示している領域の座標が自動的に入力されます。 縮小 地図移動 拡大138 (4) ダウンロードされ、表示されます 道路のデータは「inagishi lines」中に格納されています。 ※エラーが出た場合は、ネットワーク設定を確認（メニューより設定→オプションから、ネ ットワークタブ）する。もしくは、地図の表示範囲を小さくしてダウンロードする範囲を小 さくする。などの対応を行います。
18.3 道路データの抽出
ダウンロードした inagi lines は、道路データのみではなく、河川や鉄道路線のラインデー タなども含まれてしまっているため、道路データのみ抽出します。 (1) レイヤ inagi lines を右クリック→属性テーブルを開く 「highway」フィールドに値が入っているものが道路関係の地物です。 OSM からダウンロードされ たポイント、ライン、ポリゴ ンのデータが表示される。 “highway”に値が入っている ものが道路データである。139 (2) アドバンストサーチをクリックする。
(3) SQL where 句に、「highway != ''」と入力して、OK ボタンを押します。
!= は「と等しくない」という意の演算子で、highway 欄が’ ’(空)でないものを選択してい ます。 (4) レイヤ inagi lines を右クリック→選択を名前をつけて保存をクリックする。 ※フィールド欄 や演算子ボタン を選択して計算 式を作成するこ ともできます。
140 (5) ベクタレイヤに名前をつけて保存するウィンドウが開く。 ① 名前をつけて保存の「ブラウズ」を押し、「C:\GISjinzai\」に「稲城市道路.shp」 として保存する。 ② エンコーディングが「SHIFT-JIS」であることを確認する。 ③ CRS の「ブラウズ」を押す。
(6) 参照座標系の選択から、「JGD2000 / Japan Plane Rectangular CS IX EPSG:2451」を選 択して OK ボタンを押す。
① ②
141 (7) OK ボタンを押して保存する。 (8) 保存したファイルを読み込みます。 メニューからレイヤ(L)－ベクタレイヤの追加をクリックする。 ブラウズ から C:\GISjinzai\稲城市メッシュ人口.shp を選択し、「Open」ボタンを押す。 道路データが読み込まれ、表示されます。
OSM からダウンロードした「inagisi points」「inagisi lines」「inagisi polygons」は非表示にし てかまいません。また、「OSM 地物」パネルも閉じてかまいません。 非表示にする 閉じる 抽出した道路 shp が 作成される。 道路が表示される。
142
19． ＜参考③＞アドレスマッチング
アドレスマッチングとは、東京大学空間情報科学研究センターが提供する無料のサービス で、一覧データの住所情報もとに街区レベルの位置座標を付与することができます。 これにより、施設の一覧などを簡単に QGIS 上にプロットして表示させることができます。19.1 作成するデータの準備
施設の一覧を作成し、csv ファイルとして保存します。 住所を１セルに記述するようにします。19.2 アドレスマッチングでの座標データの結合
アドレスマッチングサービスを使って、作成したデータに緯度経度の情報を与えます。 (1) インターネットで「東京大学 CSV Geocoding Service」と入力し、検索する。東大のア ドレスマッチングサービスのサイトを開く。 （http://newspat.csis.u-tokyo.ac.jp/geocode-cgi/geocode.cgi?action=start）。 作成した一覧の例 （ここではコンビニの 住所一覧）143 (2) 下記のような画面が開かれるので以下のように設定する。 ① 対象範囲⇒データの範囲と測地系を選択する。ここでは「街区レベル（緯度経 度・世界測地系）」を選択 ② 住所を含むカラム番号⇒投入する csv ファイルの何列目にアドレスマッチング に使う住所情報があるかを入力する。例としたデータでは 2 列目なので 2 とす る。 ③ 入力ファイルの漢字コード⇒シフト JIS（SJIS）にする。 ④ 出力ファイルの漢字コード⇒入力ファイルと同じにする。 ⑤ マッチングオプション⇒特に操作不要。 ⑥ 変換したいファイル名⇒作成した一覧ファイルを選択する。 ⑦ 「送信」をクリックしてデータを送ると、しばらくたってからファイルのダウン ロード画面が現れるので、名前を付けて保存する。 ※文字化けしたファイルがダウンロードされた場合は、名前を適宜変更します。 (3) 保存したファイルを Excel で開く fX と fY がそれぞれ経度、緯度である。 iLvl がマッチングレベルで、1:都道府県、2:郡・支庁、3:市町村・23 区、4:政令市の区、 5:大字、6:丁目・小字、7:街区・地番、8:号・枝番となっている。 (4) Excel 上で「上書き保存」を行う。 「この形式でブックを保存しますか？」には「はい」を押す。 ダウンロードされたファイルそのままでは仕様上 QGIS で読み込むことができないため、 Excel で改めて保存することで読み込める形になります。 ① ② ③ ④ ⑥ ⑦ 緯度経度座標が 付与される
144
19.3 QGIS でアドレスマッチングしたデータの表示
アドレスマッチングで付与した緯度経度をもとに、地図上にデータを表示します。 (1) メニューより レイヤ→デリミテッドテキストレイヤの追加 をクリックする。 (2) デリミティッドテキストファイルからレイヤを作成 ウィンドウが開く。 ① ファイル名の「参照」を押し、アドレスマッチングした csv ファイルを開く ② レイヤ名を適宜設定する。 ③ 選択されたデリミタを選択し、「カンマ」のみにチェックをつける。 ④ X フィールドが「fX」、Y フィールドが「fY」になっていることを確認して、⑤OK ボタンを押す。 ① ② ③ ④ ⑤ 緯度経度が入力されてい る フィールドを選択する csv ファイルの 区切り文字を選択する （今回はカンマ）145
(3) 座標参照系の選択ウィンドウでは、「JGD2000 EPSG:4612」を選択する。
※(2)で「緯経度・世界測地系」を選択していることによります。
(4) 地図上にプロットされて表示されるので、これを shp ファイルとして保存します。 作成されたレイヤを右クリック→名前をつけて保存をクリックする。
146 (5) ベクタレイヤに名前をつけて保存するウィンドウが開く。 ① 名前をつけて保存の「ブラウズ」を押し、「C:\GISjinzai\」に名前をつけて保 存する。 ② エンコーディングが「SHIFT-JIS」であることを確認する。 ③ CRS の「ブラウズ」を押す。
(6) 参照座標系の選択から、「JGD2000 / Japan Plane Rectangular CS IX EPSG:2451」を選 択して OK ボタンを押す。
① ②
147 (7) 保存したファイルを読み込みます。 メニューからレイヤ(L)－ベクタレイヤの追加をクリックする。 ブラウズ から 先ほど保存したファイルを選択し、「Open」ボタンを押す。 施設のポイントが shp ファイルとして表示されます。 これをもとに、本マニュアル１４．や１５．などの空間的集計や圏域の分析が行えます。 作成したコンビニの データを地図上に 表示できた
