基盤地図情報1mメッシュDEMから遺跡立体図を作成する

[基盤地図情報サイト](https://www.gsi.go.jp/kiban/)で1mメッシュのDEMデータの公開が始まっています。1mメッシュのDEMをダウンロードして遺跡立体図（[永惠2024]( https://sitereports.nabunken.go.jp/ja/article/121416)）を作成する手順です。
ここで紹介する手順は永惠裕和さんの「遺跡の地形判読・記録のための、遺跡立体図・縄張図の作成」の手順に沿った再現手順です。気前よく手順を公開してくださった永惠さんに感謝いたします。

# 基盤地図情報のダウンロード

基盤地図情報サイトから基盤地図情報ダウンロードページへ進みます。

![01.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/904944f5-9de7-80b8-2f23-7a175ee84af2.png)

「ファイル選択」へ進みます。
![02.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/708fad80-57b2-a8d5-7925-385ed9ed7a21.png)

1. 「検索条件指定」→1mメッシュ
1. 「選択方法指定」→地図上で選択
白抜きのメッシュが1mメッシュが公開されているところで、グレーアウトしているところは未公開図幅

![03.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/52544543-d046-0605-25fc-7b9d273390b7.png)

- 必要な図幅を選択して、「ダウンロードファイル確認へ」をクリック
- 図幅は複数選択できる

![04.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/ea14156f-73b6-ce7f-da26-1325eff32fa1.png)

PackDLMap.zipというファイルがダウンロードされるので解凍します。
![06.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/a360ab54-67ee-a01e-14ff-7e2963231f05.png)

解凍するとさらにzipファイルとxmlファイルがいくつかできます。使うのはzipファイルです。
![07.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/3ce05a4b-10b8-6ab6-30d3-4a3d2d4d80b8.png)

# QGISでGeoTIFFに変換する

「Quick DEM for JP」というプラグインを使用します。インストールされていない方は「プラグインの管理とインストール」からインストールしてください。
「プラグイン」→「Quick_DEM_for_JP」→「Quick_DEM_for_JP」

![08.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/f46ec7e5-47ad-1f82-936e-ef67a8d2f315.png)

3点リーダーをクリック

![09.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/9b8f3948-3b20-a652-308d-277a8a0d5671.png)


右下が「*.zip」になっていることを確認して解凍したディレクトリに入っているzipファイルをクリックします。複数の図幅をまとめて読み込むこともできます。

![10.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/7b9f9829-d1ab-7c33-5ced-c697ea335710.png)

ファイル名を入力して「保存」

![11.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/9a343ef3-8953-ffd5-39fd-cd73fb635fef.png)

CRSを選択して「OK」クリック

![12.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/99c68e2e-3ad4-6bdc-8d5b-eb5e5c946292.png)

処理が始まります。1図幅あたり数分かかると思います。

![13.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/63c83eb1-8a01-2ee5-f077-852b96adc568.png)

処理が終わると、DEMが表示されます。

![14.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/bb0c9991-5fe5-9a5f-0242-6a135b9c6df2.png)

# 新たな地形指標を作成する。

## 傾斜と曲率

まずは「プロセッシングツールボックス」を開きます。

![15.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/a914079d-3ca7-3b5a-dc7a-db4daf8e3dea.png)

「GRASS」→「Raster(*.r)」と進みます。

![19.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/46b39c4e-bad2-55e1-8ab3-40eb65e89514.png)

下の方にある「r.slope.aspect」を選択します。
なお、永惠さんはSAGA GISを使用されていますが、2025年2月10日現在、私の環境ではSAGA GISをQGISから利用することができませんでした。「Processing Saga NextGen Provider」を使う方法もうまくいっていないので、代替方法を紹介します。ここではLongitudinal Curvetureの代わりにProfile Curvetureを使用します。2つの曲率はいずれも傾斜方向の断面の曲率で、細かい違いはありますがおよそ似たものと理解できます。）

![20.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/29f4c38d-18b9-1055-0ff2-f0f201c552ee.png)


1. 「Elevation」→GeoTIFFファイルのレイヤーを指定します
1. 「Format for reputing the slope」→「degrees」

![21.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/f2dd6a04-e1b4-cd67-ecf5-6aaf21238295.png)


下の方にスクロールします。
1. 「Slope」と「Profile curvature」の「アルゴリズムの終了後に出力ファイルを開く」にチェック
1. 「実行」をクリック


![22.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/e9f65884-c798-41d5-7271-fc7fd8eacc18.png)

このようなものが出力されれば成功です。

![23.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/4b7d6873-2404-0ae1-1894-f7ce66b67a3d.png)

なお、このままではQGISを終了するとデータが消えてしまうので、レイヤーを右クリック→「エクスポート」→「名前をつけて保存」からデータを保存してください。

![25.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/1113b222-02af-c5d9-fe0c-3f83e753bf5b.png)

## 地形起伏指標

「ラスタ」→「解析」→「TRI(Terrain Ruggendness Index)」と進みます。

![26.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/d77ddc36-7cea-eb79-03f9-2d6cfa88681c.png)

![27.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/a166bfbb-9dbd-deb0-e96b-ee64e04055b2.png)

以上の手順で、DEM、slope、profile curvature、TRIの4つのレイヤーが生成されれば完成です。

![28.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/0dc663b0-6f09-2148-37cf-a193f5bd7b2f.png)

# 地図の見た目を整える

## 標高データ

標高データレイヤをダブルクリックしてレイヤプロパティを開きます。

1. 「レンダリングタイプ」→「単バンド疑似カラー」
1. 「カラーランプ」→Spectral
1. 「カラーランプを反転」をクリック
1. 「モード」→「等間隔分類」
1. 「分類数」→4とか5にする
1. 「OK」

![31.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/bf14a9d9-a99d-a10e-6107-9d56daa33cbb.png)


![30.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/e8b9f62e-7921-d146-cc09-6cf01d0ce68a.png)

## 傾斜量
シンポロジタブを選択

1. 「レンダリングタイプ」→「単バンドグレー」
1. 「グラデーション」→「白から黒」
1. 「混合モード」→「乗算（Multiply）」

![31.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/f4449216-54cb-c70c-034c-2685c61233e0.png)

透明度タブを選択

1. 「不透明度」→75%
1. 「OK」

![32.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/36d57be6-511e-341f-8134-bd76a8057a4e.png)

![33.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/870c4c28-cbf7-ff5e-9a3c-120af20f7799.png)

## Profile Curvature

シンポロジタブを選択

1. 「レンダリングタイプ」→「単バンド疑似カラー」
1. 最小値＝-0.6　最大値＝0.6
1. 「カラーランプ」→「Blues」
1. 「混合モード」→「乗算（Multiply）」

![34.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/6219d86a-92b6-345e-f1df-a66018999eaf.png)


透明度タブを選択

1. 「不透明度」→75%
1. 「OK」

![35.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/1f4d4a3c-3930-cdd8-9dfb-f435bd72282e.png)

![36.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/813418e7-7c36-6ff0-a90b-6e9cc4c0e27a.png)

## TRI

1. 「レンダリングタイプ」→「単バンド疑似カラー」
1. 「カラーランプ」→「Reds」
1. 「混合モード」→「乗算（Multiply）」

![37.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/5db6e734-2ce3-cfd2-ff3b-80e4d9422271.png)

透明度タブを選択

1. 「不透明度」→75%
1. 「OK」

![38.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/6f470430-552d-45ac-81bc-963653a3cc2f.png)

## 完成！！

以上で完成です。

![39.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/8e02d2a5-7266-bfe2-e8e5-001b724caea4.png)

ここは箱館戦争遺跡として知られる二股台場周辺ですが、鶉山道を遮断する2つの大きな塹壕は十分に目視確認できます。

![40.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/422ac2f0-12ad-0bc6-9a3e-50b86729147e.png)

こちらは二股沢川に面した鉱山の火薬庫跡がよく見えます。

![41.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/8f9e0568-afb0-5348-33da-628e520305f9.png)

ここは斜面崩壊の痕跡でしょうか？

![42.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/60438/1aadfd81-e45a-e723-d2bd-a040a6acc5b2.png)

# まとめ

眺めていると興味がつきない遺跡立体図です。紹介した手順は長いようですが、慣れればどなたでも簡単にできますので、ぜひ、「お家フィールドワーク」を楽しんでください。