(富士山の)断面図をQGISで作成する

![](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/948ab61e-8611-8faa-2d5b-517630544b06.gif)

この記事は、[MIERUNE Advent Calendar 2021](https://qiita.com/advent-calendar/2021/mierune)　6日目の記事です。

## はじめに

みなさまご存知のように、あまたの芸人の顔を持つわたくしですが、QGIS講習会芸人仮面をつけてレクチャーを行う際、けっこう聞かれるのが
>断面図書けますか

というご質問です(年に2回くらい)。

そこでそっけなく

---

はい。
まず、DEM（数値標高モデル）データ準備して
そのあと[ProfileTool](https://plugins.qgis.org/plugins/profiletool/)プラグインインストールして使ってください。


　　終
制作・著作
━━━━━
　ⓃⒽⓀ

---

…と答えるのは野暮ってものです。

ここで出てくる[ProfileTool](https://plugins.qgis.org/plugins/profiletool/)は2012年にローンチされ、QGISプラグインの中でも2021/12/6現在[1548個あるQGISプラグインの中でも8位のダウンロード数(68万回)もされている](https://plugins.qgis.org/plugins/?sort=-downloads)花形プラグインなのですが、日本語解説はほぼ無さげだったりします。
なので、ここに記しておきたいとおもいます。
（英語読め (ﾟДﾟ)ｺﾞﾙｧ!! という🪓な方は初級な人類に酷すぎるのでナシの方向で🤫 ）

流れとしては、

1. 準備
    - DEMを準備する
    - プラグインを入れる
1. つかう
    - 断面引く
    - よしなに

の順番になります。

## 準備
### DEMの準備

DEMはお好きなデータをご準備ください。きれいな断面図を書きたい方はDEMが精密であればあるほど、細かい断面が得られます。
数年かけてお手製の[ミマツダイヤグラム](https://ja.wikipedia.org/wiki/%E3%83%9F%E3%83%9E%E3%83%84%E3%83%80%E3%82%A4%E3%83%A4%E3%82%B0%E3%83%A9%E3%83%A0)などを駆使して計測したデータでもよいですし、[ますますロックな静岡県さん](https://digital-construction.jp/administration/18)の[点群オープンデータ](https://pointcloud.pref.shizuoka.jp/)などからもってきてもOKです。

基盤地図情報から精密なDEMを得たい方は、ヤングなMIERUNEメンバーががんばって開発公開したプラグイン[QuickDEM4JP](https://plugins.qgis.org/plugins/QuickDEM4JP/)を使ってもよいですね。

今回は精度はともかくさっくり作りたいので、その姉妹プラグインである　[ElevationTile4JP](https://plugins.qgis.org/plugins/ElevationTile4JP/)を用いて富士山周辺のズーム14のDEMをもってきます。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/e28d43ef-574d-750d-af70-ff0ea56650a9.png)
(青富士ということで陰影透過＋青系段彩)

※DEMを容易にもって来れるすばらしいプラグインに関する解説は下記からどうぞ。

**QuickDEM4JP**
[国土地理院の標高データ（DEM）をQGIS上でサクッとGeoTIFFを作って可視化するプラグインを公開しました！（Terrain RGBもあるよ）](https://qiita.com/nokonoko_1203/items/b99aa733cb215305f8aa)

**ElevationTile4JP**
[国土地理院の標高タイルを2クリックでダウンロード→可視化するQGISプラグインを公開！](https://qiita.com/nokonoko_1203/items/51cc0990cea4c96c4565)


### プラグインのインストール
よく使う機能とは思いますが、プラグインメニュー⇒プラグインの管理とインストール
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/4f61028e-2d54-16bf-4a51-c776065fe3ac.png)

すべて⇒検索窓に`profile`と打つとprofile toolが出てきますので、選択後インストール

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/1ae5cee4-57da-b263-4fc5-db6e92ef067c.png)

メニューバーに ~~できそこないの王蟲~~ みどりのお山アイコンが出てきたらインストール成功です。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/de41f192-c45e-2968-61d2-2bfd8061fc63.png)

## つかう
### 断面を引く

まず測定されるDEMレイヤをビューで選択してください。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/d76891dd-9027-8def-3098-2c9c87d84d5c.png)

profile toolプラグインアイコンを押すと下部にウインドウが出てきますので、測定レイヤをつくるために`Add Layer`を押す
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/3023e736-2ab6-b731-d1bb-0438350f344e.png)

とりあえず富士山火口の北西付近でクリックし、南東方向でもう一度左のダブルクリックして、斜めに赤線を引いてみましょう。
もうすでに下のグラフにあの形が出てきますね。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/8555f549-3d11-6dea-eda2-428ec24e12a2.png)


`Save as`ボタンの左のセレクトから断面図の保存形式が色々と選べます。CAD使いな方はDXFなどでも良いかもしれませんね。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/0af52586-478c-df86-5a50-7c144da2dea5.png)

こちらはPNGで保存したサンプルです
![fuji.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/3b0beab1-da63-b179-5873-882274cb43db.png)

### よしなに

画像で保存だけだとちょっと物足りないという方には、数値のデータでエクスポートも可能です。`Table`タブをクリックすると、
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/190c7ee3-5eee-4921-51e3-9c3b0214c37b.png)

このような表が見えるようになり、下部のボタンからクリップボードなどに吐き出すことが可能です。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/cff543ce-b551-f042-ef26-34068bc31b6f.png)

試しに真ん中の「座標付きで」コピーを選んで、表計算ソフトにペーストするとこのようになります。縦横が変わったり、列名が入らなかったりするのでご注意を。ここから先はいかようにでもご活用ください。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/f2d26e39-b62c-2606-0d92-84101f814f38.png)

また、`Create Temporary layer`を押すと、各測定ポイントのスクラッチレイヤが作成されますので、データ確認にも使えます。

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/0a3a6f31-74bf-3501-83e2-c5b7a22c8788.png)


### このほかにも

- 今回はあえて断面ラインを手書きとして`Temporary polyline`で書いてみましたが、例えば林道ラインや河川キロポストのような、既存の選択ラインフィーチャー`Selected polyline`やレイヤ全体`Selected Layer`もProfileタブで断面として指定できます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/8c27eb08-5935-71df-9bca-4dabb7658a5f.png)

- Settingタブで、`Full resolution enabled` にチェックをつけると、最大1000ポイント/セグメントだったポイント数が、解像度の最大まで検出するようになるので、より詳細な断面がとれます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/35e2263e-42a5-9972-bfd2-17c5e72c2c41.png)


- `Show cursor`オプションをオンにしてグラフをなぞると、地図でもどの位置の断面かもわかる表示になるので、データ確認に便利ですね。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/e21289e3-e0dd-1491-a8cf-723b4df7b9c8.png)

![](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/948ab61e-8611-8faa-2d5b-517630544b06.gif)


- 同じようなエリアでラスタDEMレイヤを増やすと比較もできます。たとえば、土砂災害の起きる前と起きた後の表示なども一目瞭然になるし、数値としても残せますね。下記では仮に富士山全体をラスタ演算で250ｍ隆起させた場合のレイヤ（緑）との比較をしめしています。

 ![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/ebd09fd4-7c9f-d3a6-63b1-96826c2adf2b.png)

- 挙動がバグって適当に書いたラインのゴミが残る場合がありますが、プロジェクトを保存後QGISを再起動すると消えるのでご安心を

![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/100567/e792c1ff-c805-0b4a-48d0-870b40160c79.png)

- 英語ですが、[この動画](https://www.youtube.com/watch?v=UD0Oumv5y1w)とか説明聞かなくてもわかりやすいかとおもいます

### 最後に

けっこうProfile tool使えますね。まとめながら「そうだ！兵庫のDSMとかPLATEUでやると夜景マステとか作れて、クリスマスに向けて大富豪になれる！」と思ったのですが、最低ロット450巻の納期が一ヶ月と聞いてやめました。でもそのうちやってみるかもしれませんので、そのときまでおたのしみに。よいクリスマスを。