地理学の無料ソフトウェアを使って、地域の動物病院の情勢を理解したつもりになる方法。近所の広告を考えたり、新規開業の先生にいかがでしょうか。12ステップ。

1 QGISをダウンロードhttps://qgis.org/ja/site/

2 地図データ(shapefile)を政府統計から入手。自分の関心ある地域の地図データを選択。https://www.e-stat.go.jp/gis/statmap-search?page=1&type=2&aggregateUnitForBoundary=A&toukeiCode=00200521&toukeiYear=2015&serveyId=A002005212015&coordsys=2&format=shape

3 世帯数データ(自分の都道府県)を政府統計から入手https://www.e-stat.go.jp/gis/statmap-search?page=1&type=1&toukeiCode=00200521

4 動物病院検索地域の動物病院を探す。例えば以下のサイト。https://pet.caloo.jp/https://itp.ne.jp/

住所と名前をエクセルに保存していく。手動で十分。ウェブスクレイピング技術があれば憂いなし。初めてスクレイピングやるならExcelVBAが最短ルートと思われる。Name, place〇〇動物病院, ○○県○○市123-456〇〇クリニック, ○○県○○市789-0125 ジオコーディング以下の2つの方法のどちらかで住所から地球上の座標(地球上のx,y座標)を入手。これをジオコーディングという。

5-1 google maps geocording APIを利用する方法。

googleはWGS84(球面におかれた座標)。例Name, place,ycoordinate, xcoordinate〇〇動物病院,○○県○○市123-456, 35.1211854, 138.4506843

5-2 ジオコーディング用のエクセル(国土地理院)https://renkei2.gsi.go.jp/renkei/130326mapsh_gijutu/mapsh_gijutu190731.zip4で得た住所をこのエクセルにコピペして、座標を入手。「住所→座標値」を押す→ geojson ファイルを出力 →「動物病院.geojson」として保存。QGISでの読み込みはこちらの方が楽。

6 地図データを読み込ませる2で入手したファイルの中から、○○.shp(h27~~.shp)のファイルをレイヤパネル(左下)にドラッグドロップしてQGISに読み込ませる。

7 世帯数データを読み込ませるデリミティッドテキストレイヤを選択(コンマの形をしたアイコン)。3で入手したテキストデータ(tblT~~~.txt)を選択。読み込む際の設定は、入力shift-JIS、ジオメトリなしを選択。

8 動物病院データを読み込ませる。5-1の場合デリミティッドテキストレイヤ(コンマの形をしたアイコン)を選択。「動物病院.csv」を読み込ませる。x座標とy座標は、ぞれぞれエクセルの列(xcoordinateとycoordinate)を選択する。座標系は、WGS84を選択。これだけで、地図上の動物病院の位置に〇が表示されるはず。表示されていない場合は、このデータをレイヤパネルで、一番上にドラッグして持ってくる。この操作は、当該地図を地図を重ねる順番として一番上に持ってきたという意味。5-2の場合「動物病院.geojson」をレイヤパネルにドラッグしてQGISに読み込ませる。

9 2で得た地図データに3の世帯数データを組み込むレイヤパネルで、地図のデータ(h27~~~)を右クリック、プロパティ>結合>+ボタン>レイヤの結合で、2の世帯数データを選択(tblT~~~)>結合フィールドをKEY_CODEに設定>ターゲットフィールドをKEY_CODEに設定>OKを押すこの操作の意味は、地図データと世帯数データが地域ごとに同じ数字の列(KEY_CODE)を持っているため、2つのデータを同じ数字の列を用いてデータを結合している。

10 9で得た世帯情報付きの地図を世帯数ごとに色を付ける。レイヤパネルで、地図のデータ(h27~~~)を右クリック、プロパティ>スタイル>段階に分けられた>カラムの横にあるシグマを押す>フィールドと値でSETAIを選択しダブルクリック>OKを押す>分類を押す>OKを押す。これで、世帯数ごとに色分けが完了。

11 動物病院からの距離半径1kmを付け足す)動物病院のデータの上で右クリック>名前をつけて保存>CRSを関東であればJGD2011/UTM zone 54N 、関西であればJGD2011/UTM zone 53N、九州であればJGD2011/UTM zone 52Nに設定。新しい名前をつけて保存。ベクタ>空間演算ツール>固定距離バッファ>入力レイヤは今作ったファイルを選択>距離をメートル単位で入力。つまり1000を入力>OK。動物病院から半径1㎞の円が完成。操作の意味は、座標系を球形(WGS84)から直交系に変換して、距離の単位を経度緯度からメートルに変換。

12 ベースマップの追加)国土交通省のQGISマニュアルを参照。これの全国最新写真を選択。http://nlftp.mlit.go.jp/ksj/other/manual.pdf

以下のようなものがつくれます。名前は実在の動物病院ではありません。

category