QGISでTINを用いた標高修正ラスタの作成手順

QGISで標高ラスタを部分修正する：ラインからTIN生成までの全手順

この記事では、以下の流れで標高修正ラスタを作成する方法を解説します：

  1. ラインを作成
  2. 等間隔で点を生成
  3. 点に標高（例：40m）を一括入力
  4. TINメッシュを作成
  5. TINをラスタに変換
  6. 元のDEMと合成して新しいラスタを出力

1. ラインを作成

QGISで「ライン」レイヤを新規作成し、修正したい区間に線を引きます。保存しておきます。

2. 等間隔に点を作成

メニューの ベクタ ▶ ジオメトリツール ▶ 線上に点を配置 を使用し、例えば 10m 間隔で点を生成します。

3. 点に標高を一括入力

生成されたポイントレイヤを右クリック ▶ 属性テーブル ▶ フィールド演算機 で、以下のように数値を一括で入力します：

新規フィールド: z （実数）
式: 40

4. TINメッシュを作成

加工 ▶ TINメッシュの作成 を使って、ポイントレイヤからTIN（三角網）を作成します。

5. TINメッシュをラスタに変換

「TINメッシュからラスタを作成」ツールを使用し、値フィールドに z を選択し、出力解像度を DEM と合わせて 1m または 5m に指定します。

補足：「内挿する領域」に元のDEMを指定すると、グリッドが揃いやすくなります。

6. TINとDEMを合成（Python）

以下のPythonスクリプトで、TINの値がある部分だけをDEMに上書きし、new付きの新しいラスタを作成します。

# merge_raster_gui.py を実行
# → TINとDEMを選ぶと、"DEMファイル名 + _new.tif" で保存されます

📥 ダウンロード：

merge_raster_gui.py をダウンロード

7. 結果を確認

出力された *_new.tif をQGISで開いて、変更箇所が反映されているか確認します。

まとめ

  • QGISのGUIとPythonを併用することで、柔軟なラスタ編集が可能になる
  • 「ライン → 点 → TIN → ラスタ」という流れで、局所的な標高修正が視覚的に行える
  • 同様のスクリプトは他のファイルにも使い回し可能

この記事が、あなたの地形解析や災害対応の一助となれば幸いです。

コメント

  1. eon2025年7月15日 22:20

    import rasterio
    from rasterio.warp import reproject, Resampling
    import numpy as np
    from pathlib import Path
    from tkinter import Tk, filedialog

    # --- ファイル選択 GUI ---
    root = Tk()
    root.withdraw()
    print("🟦 部分修正ラスタ（TIN）を選んでください")
    tin_file = Path(filedialog.askopenfilename(title="部分修正ラスタを選択"))

    print("🟩 全体ラスタ（DEM）を選んでください")
    dem_file = Path(filedialog.askopenfilename(title="全体ラスタを選択"))

    # 出力ファイル名生成
    out_file = dem_file.with_name(dem_file.stem + "_new.tif")

    # --- ラスタの読み込みとリプロジェクト ---
    with rasterio.open(dem_file) as dem_src, rasterio.open(tin_file) as tin_src:
    dem_meta = dem_src.meta.copy()
    dem_array = dem_src.read(1)
    dem_nodata = dem_meta.get("nodata", -9999)

    aligned = np.full_like(dem_array, fill_value=dem_nodata, dtype=np.float32)

    reproject(
    source=rasterio.band(tin_src, 1),
    destination=aligned,
    src_transform=tin_src.transform,
    src_crs=tin_src.crs,
    dst_transform=dem_src.transform,
    dst_crs=dem_src.crs,
    resampling=Resampling.nearest,
    dst_nodata=dem_nodata,
    )

    # --- 合成処理（TINがある場所を優先） ---
    merged = np.where(np.isnan(aligned) | (aligned == dem_nodata), dem_array, aligned)

    # --- 出力 ---
    dem_meta.update(dtype="float32", nodata=dem_nodata)
    with rasterio.open(out_file, "w", **dem_meta) as dst:
    dst.write(merged, 1)

    print(f"✅ 統合完了: {out_file}")

    返信削除

コメントを投稿

このブログの人気の投稿

石川県：土砂災害（特別）警戒区域＋ＣＳ立体図

石川県の土砂災害（特別）警戒区域まとめ 2025 年 5 月 16 日 時点 石川県砂防課が公表する 土砂災害（特別）警戒区域指定状況 （PDF）を整理し、市町別に 土石流・地すべり・急傾斜地崩壊 の警戒区域（イエロー）と 特別警戒区域（レッド）を集計しました。 市町 警戒区域 (イエロー) 特別警戒区域 (レッド) 合計 土石流 地すべり 急傾斜 土石流 地すべり 急傾斜 輪島市 374 217 375 313 0 372 966 金沢市 230 82 247 170 0 235 559 能登町 240 68 218 210 0 209 526 七尾市 192 45 201 165 0 201 438 珠洲市 91 76 164 83 0 163 331 小松市 181 13 66 162 0 64 260 白山市 139 20 91 73 0 84 250 志賀町 130 15 100 106 0 99 245 穴水町 93 12 115 79 0 113 220 加賀市 78 21 101 65 0 96 200 津幡町 78 48 57 56 0 56 183 中能登町 130 4 36 99 0 35 170 能美市 64 1 41 59 0 40 106 宝達志水町 33 30 53 26 0 52 116 かほく市 12 2 30 5 0 30 44 川北町 0 0 0 0 0 0 0 野々市市 0 0 0 0 0 0 0 内灘町 0 0 5 0 0 5 5 合計 2065 654 1900 1671 0 1854 4619 ＣＳ立体図＋土砂災害警戒区域 MAP を開く 関連資料・データ ▸ 土砂災害（特別）警戒区域指定状況 PDF（R7.5.16） ▸ 石川県土砂災害（特別）警戒区域図 オープンデータ（SHP・GeoJSON 等） ▸ CS 立体図を使った 地形判読マニュアル（林野庁） ハザードマップの入手方法 お住まいの 市役所・町役場の防災担当窓口 や自...
続きを読む

日本でよく使う EPSG コード 一覧

1. 経緯度（地理）座標系 区分 名称 EPSG 主な用途・備考 世界標準（地理） WGS 84（緯度・経度） 4326 GNSS/GPS の測位やデータ交換で使用。単位は「度」。 Web地図（投影） WGS 84 / Pseudo-Mercator（Web Mercator） 3857 Google/OSM/国土地理院タイルなどの 表示 に使用。単位は「メートル」。別名：Popular Visualization / Google Mercator。 日本最新 JGD2011（緯度・経度） 6668 測地成果2011。公共測量の基準。平面直角座標系は下表参照。 過去の全国基準 JGD2000（緯度・経度） 4612 2002–2011 年に使用。JGD2011 と混在時は変換推奨。 旧日本測地系 Tokyo Datum 4301 1940年代〜2002年まで使用。WGS84とは数百mのズレあり。 2. 平面直角座標系（Gauss–Krüger 投影） 系 (ローマ数字) 系 (数字) 中央経線 JGD2011 EPSG JGD2000 EPSG 適用都道府県（概要） I 1 129°30′E 6669 2443 長崎、鹿児島西部 II 2 131°00′E 6670 2444 福岡・佐賀・熊本・大分・宮崎、鹿児島...
続きを読む

設計定数を求めるための代表N値について

代表 N 値（設計用採用値）の決め方 ― 要点まとめ 1．実務でよく尋ねられる「出典」の悩み（改訂） 発注者から「この整理手法の出典はどこか」と問われることが多い一方で、道路土工・道路橋示方書・地盤調査法など主要文献には 「平均値±1/2σで採用する」と明記した箇所はほとんど見当たりません。 しかし、 『報告書作成マニュアル［土質編］』（社団法人 全国地質調査業協会連合会） には、 ばらつきデータから設計用値を定める 三手法 （ a： 平均、 b： 最小／最大、 c： 平均と標準偏差による補正）が明記されています。 特に c) として次式の例が示され、 設計値 ＝ 平均値 − (1/2) × 標準偏差 （平均値を標準偏差で補正して設計値とする） と記載されています（第2章 2.7.3、p.74）。 したがって、本稿で用いる 「平均 − 1/2σ」 という整理は単なる経験則ではなく、 公的マニュアルに裏づけのある実務手法 です。 なお同マニュアルは、標準偏差の重み係数（1/2 など）の選び方には議論があること、すなわち k の選定はデータ量・ばらつき・重要度等に応じて調整する余地があることも併記しています。 出典： 『報告書作成マニュアル［土質編］』第2章 2.7.3「設計用土質定数の設定方法」p.74（a〜c の三手法と「平均 − (1/2) 標準偏差」の例示） 2．建築基礎構造設計規準・同解説（1974）の記述 地盤の場合、土質試験などによって得られる諸数値にかなりのばらつきがあるのがむしろ普通であり、このようにばらつきのある測定値の取り扱いについては、ばらつきの原因とその程度、基礎底面と各試料採取位置との相対的な関係など、複雑な要素が絡むため一概には言えない。しかし通常、1地点で同一土質層について得られた測定値は、その平均値を用いて差し支えないとする。そのうえで、とくにばらつきが著しい場合などは、平均値よりさらに安全側を見込んで 設計用採用値 ＝ （平均値） ± （標準偏差）／2 とする。ここで複号 ± は、算定結果が安全側となる符号をとる...
続きを読む