令和になってからQGISを使っている..ではじまるQGISの記事を別のブログからこちらに転載します。元の記事は2020年5月に書いたものです。

河川の縦断面図をつくりたくて試行錯誤した。縦断面図という表現で合っているかわからないけど、河口からの距離と標高のグラフのやつ。よく、ヨーロッパのセーヌ川なんかに比べると日本の川は滝だ、という表現で示されるあれ。

今回はQGISと国土数値情報を使って作成した。ここに方法をメモしておく。
おおまかには、国土数値情報のシェープファイルをラスタ化して、QGISのプラグインである「Plofile tool」や「qProf」で河川等の任意の断面図を作成する。

ちなみに、国土地理院の基盤地図情報をダウンロードして、フリーのツールでDEMにして…という方法もあるみたいだけど、国土数値情報だけでできました。まぁ、元データは同じだろうから変換の過程が違うだけだろう（詳しいことは知らない）。

手順

まずはデータの準備。国土数値情報ダウンロードから標高・傾斜度3次メッシュデータをダウンロードする。ダウンロードしたZIPファイルを展開する。そして、QGISで表示する。
ベクタレイヤのマージを行う。（しなくてもよいがしておくと便利）
「プロセッシングツールボックス」－「ベクタテーブル」－「文字列を浮動小数点に変換」で、標高のフィールド（この例ではG04a_002）を文字列から数値に変換する。これをせず文字列のままだと、ラスタ化できない。
「ラスタ」－「変換」－「ベクタのラスタ化」を実行する。オプションは次の通り；
出力ラスタの単位：地理単位
水平方向の解像度：0.001
鉛直方向の解像度：0.001

ちなみに、解像度を適切な値にしないとうまくいかない（ここでハマりました）。地理単位が度なのかkm（キロメートル）なのかによっても違うと思うので、留意されたい。

QGISの画面右下のアイコンをクリックして、プロジェクトの座標参照系を「JGD200/Japan Plane Rectangular CSⅧ」あたりにする。これをすることで、距離の単位を度からkmにかえる（最初、距離が度でしか表示されず、ハマりました）。
プラグイン「Profile tool」を起動して、作成した標高のレイヤ（ラスタ）を追加する。あとはProfile toolを操作して必要な断面を指定、表示する。

豆知識