qgisで地理院地図の空中写真を使うサンプル

#空中写真とは
　航空測量や地図作成のため飛行機などで撮影した写真のこと。
　詳しくは国土地理院の[解説ページ](https://www.gsi.go.jp/CHIRIKYOUIKU/syasin.html)を参照のこと。

#どこで入手するのか?
##地図・空中写真閲覧サービス
　普通の人が自由に見られる空中写真は、まず[国土地理院の空中写真閲覧サービス](http://mapps.gsi.go.jp/maplibSearch.do#1)があげられる。
 閲覧する前に利用条件の確認が行われるので、確認の上、同意する。
　そうすると、日本地図が表示されるので、見たいところを選ぶ。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/c58cafb7-f4b8-130c-89eb-fe3577612102.png)
　画面の左側に次のような選択画面があるので、見たい年代等を設定して検索する。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/94a18d8f-cad1-07cc-89cf-bf5a3def8bbc.png)
　そうすると、次のように候補が表示される(下図は秋田駅周辺)。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/68885051-2f0e-bed5-06f7-2739a7a39b3a.png)
　みたいところを選択すると、単写真で見られる。
　必要な写真は[(一財)日本地図センターのネット通販](http://www.jmc.or.jp/phototsuhan.html)で購入できるとのこと。

##gisソフトで表示する
　昔の写真を見ていて、現在と比較したい場合が良くある。
　この場合、地理院地図を利用する。
　秋田市の場合、現時点では1960年代と1975年前後の空中写真が閲覧できる。

#QGISの背景地図として利用する
##qgisの対応状況
　現在（2019/8/18現在）、公開されているバージョンのWindows用では利用できる。
　他のOSで試したところ、次のとおり。
　FreeBSD 11.3のパッケージでは不可
　Raspbian Busterではqgis2.18がインストールされる。このパージョンでもWindowsなら利用できるはずだが不可。
　現時点ではWindows版を利用するのが間違いないみたい（Macは持っていないので不明）。
##qgisでの設定
　qgisにレイヤを追加する「ブラウザ」に必要な情報を設定する。
　ブラウザの【XYZ Tiles】に地理院地図の情報を設定する。
　完成図は次のとおり。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/c0e1672a-451d-ffc0-d9f5-1ca288b98c52.png)
　まず【XYZ Tiles】を選択し、マウスの右ボタンを押すと次のダイアログが開く。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/dc79ec0d-882d-5a0b-7f01-7ea7f7db0bb1.png)
　名前の項目は自分がわかりやすい情報を入れる。
　URLに入力する情報は[国土地理院のページ](https://maps.gsi.go.jp/development/ichiran.html)で調べる。
　空中写真は次のあたりにまとまっている。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/85c35d27-58af-b2c0-ce36-7c7c51d81259.png)
　例えば「空中写真（1961～1969）」のリンクを開く。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/8a8bbd46-e0bf-95dd-8a4b-5b7496a9d81a.png)
　ここのURL、ズームレベルを入力する。
　![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/cab79457-f4ca-ef46-b96e-cced3fb4cbc8.png)
　【OK】を押すと【XYZ Tiles】に追加されているので、これをダブルクリックすると、【レイヤ】に追加される。
　この際注意点として、どこを描画しようとしているのかわからないので、先に標準地図を設定しておき、見たいところを表示しながら作業した方が良い。
『ポイント』
　まず標準地図を表示し、見たいところを指定する。
　標準地図を指定しただけだと、このように世界地図になってるので、提供範囲の限られている空中写真は表示されませんね。当然（反省）。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/ecbf0d0d-81c2-c523-08f3-8e55a3bed715.png)

　表示するところを、こんな感じで指定する。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/39b73390-a4e0-2f5e-eb31-f0bf916c30ca.png)
　次に空中写真をダブルクリックする。するとめでたく表示される。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/31f49623-4e50-2e9c-615b-e22d64b5a5b8.png)
　ズームインすると、こんな感じで見られます。これは秋田駅の周辺。
　機関車トーマスでおなじみの扇形車庫と転車台があったことがわかります（なお、[小樽市市立総合博物館](https://www.city.otaru.lg.jp/simin/sisetu/museum/iron_horse/)では、軽便鉄道で利用されていた蒸気機関車と転車台などが動態保存されていて、乗せてもらうこともできるのでお薦めです）。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/ebf3a1f6-b7e9-077b-3b38-2285bebb91db.png)

##見るだけでも楽しいけど
　昔の写真を見るだけでも楽しいけど、現在との比較とかできれば、もっといろいろな情報が見えてくる。
　まず背景地図にしている標準地図を見えるようにしてみる。

###標準地図を透かしてみる。
　空中写真のレイヤの透明度を調整すると、背景の地図が見えてくる。
###自分の情報を載せてみる。
　昔の空中写真を見る場合、現在との比較を行いたい場合が多いと思う。
　例えば、家を建てるときの土地選びでは、学校や駅が近くスーパーもあるとか便利の良さに目が行きがちだけど、建物の寿命を考えたらできるだけ地盤が良く、水はけが良い土地を選んだ方が良い。
　昔の地形図を探して調べるとか、もっと手軽には地名で推測するという方法もある。例えば[『地名に隠された由来、旧地名の調べ方』](https://allabout.co.jp/gm/gc/25963/)のような記事を参照のこと。
　しかし、昔の空中写真を見れば、数十年前までという限界はあるけど、その土地がどんなところだったか一目瞭然になる。さかのぼれる時代が数十年前といっても、地方都市の場合、宅地化が進展したのは1960年代以降のところが多いわけだし、実用的には問題ないと思う。
　さて、具体的にどうするか。
　実際に調べたいポイントを表示するレイヤを作成し、空中写真を背景地図にすれば良い。
　まずqgisのメニューバーから次のアイコン![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/5cf8f069-0c03-28ef-b32d-2a5af06e8d1c.png)をクリックする。
　次のようなダイアログが開く。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/a440a474-e342-53e6-8343-2b512de3d8ce.png)
これはShapefileを作成する。
　ファイル名には作成したデータを保存するファイル名、ファイルエンコーディングは文字のエンコーデイング、ジオメトリタイプはデータ型（ポイント、ポリゴンなど）を指定する。
　ここでEPSGを適切に指定しないと、作成したデータが地球上のあらぬところに飛んでいってしまうので注意。
　私の場合、公共座標系がなじみ深いので指定している。なお公共座標系のEPSGは住んでいる（おおよそ）都道府県で異なるので適切なものを調べること。
　ShapefileはdBaseのデータファイルを構成要素にしており、制限付きながら地物の属性を保存できる。フィールドに保存したいデータ型に応じて追加する。このあたりはAccessのテーブルと考え方は同じ。
　次に編集したいレイヤを指定した上で、qgisを編集モードにする。
　マウスで編集したいレイヤをクリックすると､次のようになる（例では今作成したレイヤ）。
　![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/2b5f9029-18c9-47b2-be3c-ba6dcd4f3dd5.png)
　次にメニューバーの鉛筆![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/f83cb36d-1706-fd6b-be22-dfe9c76682d9.png)をクリックする。
　メニューバーで編集に関連するアイコンが有効になる。
　![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/af2ae4d0-3ed9-ff08-33f2-ac7699877ddc.png)
　![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/bcc41b3c-66c6-a535-ced7-644519aa7a5a.png)のアイコンが地物を追加するアイコン。
　これでポイントを追加すると、追加した地物（ポイント）の属性を入力するダイアログが開く。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/89b3f064-8df7-c7a2-36d4-6b38640e535f.png)
　【OK】を押し、FDのアイコン（保存）を押す。
　サンプルとして、秋田駅東口のアルヴェにポイントを打ってみる。
　作成した地物が見えない場合、レイヤの上下やシンポロジを確認する。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/791be5cd-cb19-0281-0b33-f204fcc086c7.png)
　ここまでできたら、背景地図を過去の見てみたい空中写真に変える。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/a25fe834-f444-f44d-39ed-f696913d39c2.png)
　クリップの位置がずれてしまった・・・。
　現在、アルヴェがある位置には1960年代、転車台と扇形車庫があったことがわかる。
　基盤地図情報の建築物を重ねてみると、こんな感じ。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/347462/f1f62d6d-2046-3c8b-13bb-c8f4033cb793.png)
　今と昔を比較するだけでも、結構楽しいと思います。