【QGIS】PLATEAUデータを表示する

## はじめに
この記事は#30DayMapChallenge2022 27日目の記事です。
テーマはMusicです。
QGISを使って[音楽のまち](https://www.ongakunomachi.jp/)、川崎市のPLATEAUデータを表示してみます。

![image](https://raw.githubusercontent.com/tjukanovt/30DayMapChallenge/main/images/flyers/30dmc-2022.png)

https://30daymapchallenge.com/

## QGISとは
フリーでオープンソースの地理情報システム
ほとんどのUnixプラットフォーム、Windows、macOS上で動作する

## PLATEAUとは
国土交通省が主導する3D都市モデル整備・活用・オープンデータ化プロジェクトです。
3D都市モデルには、家屋やビルなどの「建築物」、都市計画区域などの「都市計画決定情報」、土地の用途を示した「土地利用」、都市のインフラである「道路」や「橋梁」など、都市に関するデータが格納されています。

## PLATEAUデータについて
### データ形式
PLATEAUデータには、様々なデータ形式があります。
例えば、CityGML、3D Tiles、MVT、ファイルジオデータベース、OBJ、FBXなどなど。。
今回はCityGMLを採用します。

### CityGML
OGC（Open Geospatial Consotium）で提唱されている、3D都市モデルを扱うための標準的なフォーマットです。
XML形式のデータです。
CityGMLは、CityGML 1.0、CityGML 2.0、CityGML 3.0の3種類があり、現在3D都市モデルはCityGML 2.0に準拠しています。


### データの取得
データは[G空間情報センター](https://www.geospatial.jp/)からダウンロード可能です。
また、PLATEAU公式が[PLATEAU配信サービス（試験運用）](https://github.com/Project-PLATEAU/plateau-streaming-tutorial)で3D Tilesを配信しています。
今回は、G空間情報センターからCityGML形式のデータをダウンロードします。

### LOD
LOD（Level of Detail）とは、3D都市モデルに含まれている地物の詳細度です。
建築物はLOD0からLOD4までが定義されています。

![image](https://www.mlit.go.jp/plateau/uploads/2022/11/about_citygml_02.png)
画像の出典：[国土交通省 Project PLATEAU](https://www.mlit.go.jp/plateau/learning/tpc01-2/#p1_3_2)

## 事前準備
既にインストール済みのQGIS（バージョン3.22.3）を使います。
QGISを開きます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2667439/13c801a9-4ba2-7996-cb45-16f32c999cd1.png)

## データダウンロード
G空間情報センターから[川崎市のPLATEAUデータ](https://www.geospatial.jp/ckan/dataset/plateau-14130-kawasaki-shi-2020)（CityGML形式）をダウンロードします。
ダウンロードしたZIPファイルを展開します。
ファイルサイズが大きので、ダウンロードと解凍に少し時間がかかります。。

## データ表示
今回は、建物データを表示してみます。
建物データは以下のディレクトリにありました。
```
14130_kawasaki-shi_2020_citygml_5_op\udx\bldg
```

gmlファイルを選択してドラッグ&ドロップします。
または、レイヤタブからレイヤを追加>ベクタレイヤを追加をクリックして、gmlファイルを選択します。
![image](https://user-images.githubusercontent.com/83101503/204134998-b96a6c03-1c6e-4fc7-bc09-0d7dacec29b9.png)

ソースタイプはファイルを選択し、ソースからgmlファイルを選択します。
入力できたら「追加」をクリックします。
![image](https://user-images.githubusercontent.com/83101503/204135039-7caa13a0-6aec-4ce2-a1f7-8e400c45508e.png)

![image](https://user-images.githubusercontent.com/83101503/204135295-e700af6e-d6de-4ee0-b3c2-ece77851a54e.png)

表示できました！
地物が少ないのですぐに表示できましたが、地物が多いと表示だけでかなり時間がかかることがあります。。
LOD1の地物は表示できますが、LOD2の地物は表示できないようです。

## 3D表示
レイヤを右クリックして、プロパティをクリックします。
単一定義（single）を選択し、OKをクリックします。
![image](https://user-images.githubusercontent.com/83101503/204135470-1e0bdf96-e480-479e-9df9-c93d4abc8145.png)

プロジェクトの座標参照系を変更します。
左下の「EPSG:6668」と書かれている欄をクリックします。
6677を適用して、OKをクリックします。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2667439/afbfdea5-f8f7-261b-af5b-5a7cd1df1ebc.png)

ビュータブから、新規3Dマップビューをクリックします。
![image.png](https://user-images.githubusercontent.com/83101503/204135803-15091e9d-4893-4bcb-8d53-48a72ebd95e6.png)

3Dで表示できました！
QGISのプラグイン「Qgis2threejs」を使うと同様に3D表示が可能です。（私はこちらの方をよく使います）

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2667439/908c178b-adcb-0e9c-a263-281cd63547f7.png)


## 参考文献
https://www.mlit.go.jp/plateau/learning/tpc03-1/

https://www.mlit.go.jp/plateau/learning/tpc03-2/

https://www.mlit.go.jp/plateau/learning/tpc05-1

https://www.mlit.go.jp/plateau/learning/tpc01-1/

https://www.mlit.go.jp/plateau/learning/tpc01-2