新鮮とれたて最新国土地理院空中写真をQGISで読み込みたい

![Jan-25-2023 23-45-26.gif](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/3efeaf68-bca6-071c-33a5-64075a671fe7.gif)




:::note info
意外と知られていない(n=5, [へーシャ](https://mierune.co.jp))ので、インターネッツの海に放流してみます。
:::


## はじめに
国土地理院さんは長年様々な[地図タイル](https://maps.gsi.go.jp/development/ichiran.html)を無料で公開してくださっています。
ジオ界の国民栄誉賞ものです。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/af047ee9-3f4f-f20a-1807-46fcb19cf6c5.png)


たとえば、災害に関するタイルなどは現場の様々なかたのご尽力により発災後すぐアップされます(たとえば[こことか](https://maps.gsi.go.jp/development/ichiran.html#t20190618yamagata_tsuruokamurakami_0620do))。311の時に全国有志でﾋｲﾋｲ言いながら現地の凄惨なオルソ写真をインターネットバケツリレーしていた老兵としては隔世の感があります。皆様ありがとうございます。

あと、QGISのビギナーな方はここから空中写真や標準地図(地形図)などのURLをコピーして、多くの業務や可視化に役立てていることでしょう。やり方は[ここにありますし](https://gis-oer.github.io/gitbook/book/materials/QGIS/QGIS.html#%E5%9C%B0%E5%9B%B3%E3%82%BF%E3%82%A4%E3%83%AB%E3%81%AE%E8%AA%AD%E3%81%BF%E8%BE%BC%E3%81%BF)、一覧のXMLを読ませれば一括登録できたりします(すみませんそのうちupdateします)。こんな地図タイルがない時代なんてジオジサンたちはですね、その、CD-ROMという円盤状の七色に輝くメディアから…(約2000文字略)

さてこのタイルリスト、長年多くの優良データが蓄積されてきた分、よーく読まないと見落としてしまう情報が時々あります。
今回は最新の空中写真を見るためのTipsについてご紹介したいとおもいます。

## 空中写真タイルはたくさんあります

上記の地理院さんが配信してくださっているタイルのうち、よく使われる空中写真タイルとしては

- [ベースマップ(シームレス)写真](https://maps.gsi.go.jp/development/ichiran.html#seamlessphoto)
`https://cyberjapandata.gsi.go.jp/xyz/seamlessphoto/{z}/{x}/{y}.jpg`

- [簡易空中写真（2004年～）](https://maps.gsi.go.jp/development/ichiran.html#airphoto)
`https://cyberjapandata.gsi.go.jp/xyz/airphoto/{z}/{x}/{y}.png`

- [電子国土基本図（オルソ画像）（2007年～）](https://maps.gsi.go.jp/development/ichiran.html#ort)
`https://cyberjapandata.gsi.go.jp/xyz/ort/{z}/{x}/{y}.jpg`

などがあります（このへんのネーミングはちょっと初学者にはわかりにくいので将来的な整理が望まれますね）
あとほかにも、[米軍撮影による終戦直後の写真](https://maps.gsi.go.jp/development/ichiran.html#ort_USA10)なども提供されており、QGISのハンズオンなどをではこれをみると結構もりあがります。

## 最新空中写真を使うには

しかし、上記のベースマップの写真などを使っている方から、これ最新の撮影データでないんじゃないのーという声をちらほら聞きます。
たしかにそれぞれのタイルの撮影年度はすぐわからないのですがけっこう古い状況のママの地区もあったりして、目的に合致しない場合があるようです。

::: note warn
だからといってG●●mapsのタイルがわーキレイだー!あたらしー!と読み込んでトレースしたり成果物に使う~~暴挙~~冒険はおすすめできません。
節子、それあかん。読めちゃう/使って良いは別の問題だよ！　ジオジサンとの約束だよ！
参考：[行政サイトでウェブ地図を使う際の注意点などまとめ](https://qiita.com/nyampire/items/5fd06107f25bc12a526f#google-maps)
:::

それはさておき、最新の空中写真がほしいというあなたに今回紹介したいのは国土地理院さんが令和2年(2020年)に提供開始された[年度別空中写真](https://maps.gsi.go.jp/development/ichiran.html#nendo)というタイルについてです。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/031d8418-611e-48a1-9bf6-44a9ef6682f5.png)

タイルURLはこのように定義されています
`https://cyberjapandata.gsi.go.jp/xyz/nendophoto{year}/{z}/{x}/{y}.png`
この説明だけ読むと、ちょっとすぐ意味がわかりにくいかもしれません。
1. `{year}`は`{z}/{x}/{y}.png`みたいにそのままでよいの？
2. 何年にどこで撮られているの？
3. ベースマップ写真タイルとかに反映されてるんじゃないの？


### 読み方

まず、HPにも書いてありますこの定義を
`https://cyberjapandata.gsi.go.jp/xyz/nendophoto{year}/{z}/{x}/{y}.png`
をQGISで読むには、`{year}`部分に読みたい年代を記述するとOKです、{}も不要です。

たとえば2021年度であれば
`https://cyberjapandata.gsi.go.jp/xyz/nendophoto2021/{z}/{x}/{y}.png`

をタイルURLとしてQGISなどに読み込めば閲覧可能です、でも西暦3250年とかはまだ提供されてないようです、ごめんなさい。

そして、QGISで見るには、経験者さんがよくやる[このやり方](https://gis-oer.github.io/gitbook/book/materials/QGIS/QGIS.html#%E5%9C%B0%E5%9B%B3%E3%82%BF%E3%82%A4%E3%83%AB%E3%81%AE%E8%AA%AD%E3%81%BF%E8%BE%BC%E3%81%BF)で
- [地理院淡色地図](https://maps.gsi.go.jp/development/ichiran.html#pale)
`https://cyberjapandata.gsi.go.jp/xyz/pale/{z}/{x}/{y}.png`
- 2021年度のタイル
`https://cyberjapandata.gsi.go.jp/xyz/nendophoto2021/{z}/{x}/{y}.png`
を読み込んでみましょう

…何も写真が出ない…
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/f62b7046-d316-4865-738e-fa0730db806a.png)

これは仕様どおりで、ズームレベル14以上(QGISだと1/70000くらい？)に拡大しないと表示がされないようになっています。しかし、ズーム14をしらみつぶしに探すのはなかなか大変です。
そこで、撮影年度がわかるサイトがありまして、こちらで見ることもできます

[国土地理院撮影・航空レーザ範囲](https://psgsv2.gsi.go.jp/koukyou/kihonsatsuei/index_photo_area2.html)

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/341416df-5f48-c12b-d61c-a87eff2343c4.png)

…しかし、ちょっと使いづらい、しかもレイヤは和暦でタイルURLは西暦で…しかも、KMLて…[とても、つらい](https://www.google.com/search?q=%E3%81%A8%E3%81%A6%E3%82%82%E3%81%A4%E3%82%89%E3%81%84+%E4%B8%89%E5%9B%BD%E5%BF%97&tbm=isch&ved=2ahUKEwjsyJDB_OL8AhWDAKYKHeKpAkEQ2-cCegQIABAA&oq=%E3%81%A8%E3%81%A6%E3%82%82%E3%81%A4%E3%82%89%E3%81%84&gs_lcp=CgNpbWcQARgCMgQIIxAnMgUIABCABDIFCAAQgAQyBwgAEAQQgAQyBQgAEIAEMgUIABCABDIFCAAQgAQyBQgAEIAEMgcIABAEEIAEUABYAGDaDWgAcAB4AIABWogBWpIBATGYAQCqAQtnd3Mtd2l6LWltZ8ABAQ&sclient=img&ei=L0LRY6z1CYOBmAXi04qIBA&bih=1216&biw=1499)

どこぞにこの索引GISデータとかタイルがないかと探したのですが、よくわからなかったです（あったらごめんなさい）
なので、カッとなってついうっかり現状のKMLぜんぶ引っこ抜いて1つのGeoJSONにしました。[こちらからどうぞ](https://github.com/Yfuruchin/GSI_photo_tile/blob/main/H19_R4_GSIPhoto.geojson)（あくまで2023/1現在のものですのでご了承ください）

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/ffc7a6e5-8f2e-8c00-4e32-eafc3393ae2e.png)


このGeoJSONをDLしてきてQGISで読んで色付けなどするとこんなかんじになります
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/7e5f27fa-b2c1-0676-d4b0-1e5990fcb954.png)

…おお、けっこう撮られていますね。これがあるとほしい年度での撮影位置もわかるので助かるかと。あと、R4-R5年に撮影される（された）レーザー計測範囲なども入っています。

### 確認

手前味噌ですが、北海道小樽らへんも2021年に撮影されているので、先程の2021年度タイルを上に載せた状態でズームしていきますと、2018年に開通した[小樽ジャンクション](https://ja.wikipedia.org/wiki/%E5%B0%8F%E6%A8%BD%E3%82%B8%E3%83%A3%E3%83%B3%E3%82%AF%E3%82%B7%E3%83%A7%E3%83%B3)がはっきりと写っています。

![2021.gif](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/85b2406d-5cf4-3682-74d8-e00157b6e722.gif)


その一方ベースマップ(シームレス)写真では反映されていないことがわかります。おそらくちょっと前の空中写真が使われている理由のようです。

![シームレス.gif](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/9336bef9-aeda-6bcd-9e75-0315c583dca6.gif)


もちろんすべてのタイルを比較をしたわけでもなく、このあたりの撮影ログは別途どこかにあるとは思いますが、「あー地理院の写真新しくないよねー」とお嘆きの方がいましたら、G●●gl…などをQGISで読み取るなどといったスリル満点のアウツオブガバナンスな危ない橋を渡るのではなく、一度拙GeoJSONとともに地理院さんの年度別レイヤを見ていただけるとよいのではないでしょうか？

※ちょっと見た感じですと、場所によってはシームレスと2021年度がマージされてたりしたので、まずは比較されてみることをおすすめしますし、年度ごと撮影範囲とあわせることで撮影年の推定には使えるかもしれません（しかしなんで地域によってマージされてるされてないがあるのか謎なのですが）。

**2023/1現在、2022年度のデータはまだどちらのタイルにも未反映です

## まとめ

ということで
1. `{year}`は`{z}/{x}/{y}.png`みたいにそのままでよいの？
→
`https://cyberjapandata.gsi.go.jp/xyz/nendophoto2021/{z}/{x}/{y}.png`
とかでQGISに放り込みましょう

2. 何年にどこで撮られているの？
→
[このサイト](https://psgsv2.gsi.go.jp/koukyou/kihonsatsuei/index_photo_area2.html)で探すか
[このGeoJSON](https://github.com/Yfuruchin/GSI_photo_tile/blob/main/H19_R4_GSIPhoto.geojson)をQGISで読み込んで確認しましょう
 
3. ベースマップ写真タイルとかに反映されてるんじゃないの？
→概ね反映されているぽいではありますが、一部反映されてない場所もあるので、確認するのはアリかとおもいます。あと、マージされると撮影年度が確認できなくなるので、逆算で辿れるのは便利かもしれないですね

+α.　[ベースマップ(シームレス)写真](https://maps.gsi.go.jp/development/ichiran.html#seamlessphoto)もそこそこマージされてるのでおすすめです、G●●leを読んでくるよりよほどよい(しつこい)
`https://cyberjapandata.gsi.go.jp/xyz/seamlessphoto/{z}/{x}/{y}.jpg`

----

出典:[地理院タイル](https://maps.gsi.go.jp/development/ichiran.html) 
いつもありがとうございます