砂防業務関連のデジタル化ノート

iPhone LiDAR（Scaniverse）とCloudCompareで作る！手軽で高精度な3D断面➡2D断面図の作成手順 iPhoneアプリ「Scaniverse」で取得した点群データから、CADで使える2次元の横断図（距離と標高のプロット）をサクッと作成する一連のワークフローをご紹介します！ 現地での調査からCAD作業に入るまでの時間を劇的に短縮できるので、ぜひ試してみてください。 STEP 1：Scaniverseでデータ取得 まずは現場でiPhoneを使い、対象となる地形や水路などをScaniverseでスキャンします。取得した3Dモデルは、点群データとして「LAS形式」でエクスポートし、PCに取り込みます。 STEP 2：CloudCompareへの取り込みとノイズ除去 PCで無料の点群処理ソフト「CloudCompare」を立ち上げ、LASデータを読み込みます。 自然地形をスキャンすると、雑草や落ち葉、スキャン時のブレによる空中のノイズが含まれるため、以下の方法で事前に綺麗にしておきます。 SORフィルターで自動除去： Tools > Clean > SOR filter を使い、空中に浮いた細かい孤立点を一掃します。 ハサミツールで手動カット： 明らかに邪魔な草などは、画面上部のハサミアイコン（Segmentツール）で囲って切り取ります。これが一番確実です。 STEP 3：3D断面（DXF）の抽出 綺麗になった点群から断面を切り出します。 メニューの「Cross Section（断面ツール）」を起動。 枠（クリッピングボックス）のサイズを調整します。ここで重要なのが 「スライスの厚み（Z軸など奥行き方向）を限界まで薄くする（デォでok）」 こと！分厚いと手前と奥のノイズを拾ってしまい、断面がギザギザになってしまいます。 パネルの「点群の抽出ボタン（緑の下矢印）」を押し、切り出された点群をDXF形式で保存します。 (注意：この時、出力された3DのDXFをLibreCADなどの2D専用CADで開いて上書き保存しないでください。高さデータが消えてしまいます) STEP 4：Pythonで2D DXFへ一括変換 CloudCompareから出力し...

岐阜県の土砂災害・水害情報をまとめたWeb地図 岐阜県の土砂災害・水害情報をまとめたWeb地図 岐阜県の土砂災害警戒区域・特別警戒区域を中心に、岐阜県CS立体図、洪水・内水などの水害情報、 避難場所、昔の地形図、旧市町村境界、地価、バス停などを一つの画面で切り替えて確認できるWeb地図です。 土砂災害区域は箇所番号・箇所名・住所から検索でき、地図上をクリックすると住所・緯度経度・ JGD2011平面直角座標もその場で確認できます。 岐阜県 土砂災害・防災Webマップ 地図を大きく表示 別タブで開く スマートフォンでの現在地表示： 地図左上の現在地ボタンを押すと、端末の位置情報を利用して現在地を表示できます。 初回はブラウザから位置情報の利用許可を求められます。 クリックした場所の住所・座標をその場で確認 地図上の任意の場所をクリックまたはタップすると、その地点の 住所・緯度・経度・JGD2011平面直角座標 を画面左下に表示します。 山間部や斜面でも「この場所はどの住所付近か」「座標はいくつか」を一度に確認できます。 平面直角座標はクリックした緯度経度から自動で変換します。 通常であれば、緯度経度を確認したあと国土地理院などの座標変換ツールへ入力して 平面直角座標を求める作業が必要ですが、このWeb地図ではその部分を自動化しています。 岐阜県ではJGD2011の第VII系（EPSG:6675） を使用します。 資料作成にも便利： 表示された住所・緯度経度・平面直角座標は文字として選択できます。 PCではドラッグ、スマートフォンでは長押ししてコピーできるので、 Excel、QGIS、メール、調査メモなどへそのまま貼り付けられます。 なお、Web地図上でクリックした地点を変換した参考値であり、測量成果そのものではありません。 岐阜県のGIS・オープンデータを活用 岐阜県には、行政情報や地域情報を地図上で公開する 「岐阜県域統合型GIS」 があり、さらに岐阜県オープンデータカタログでは 土砂災害警戒区域・特別警戒区域のGISデータなども公開されています。 土砂災害警戒区域のGISデータは、岐阜・西濃・中濃・東濃・飛騨...

1. 国土地理院ベクタータイルの「拡大すると消える」問題を解決する QGISで「国土地理院ベクター（最適化ベクトルタイル）」を表示した際、標準スタイルでは ズームレベル16.1以上に拡大すると建物や道路が消えてしまう 仕様になっています。これでは詳細な図面作成ができません。 ベクタタイルレイヤはプロパティから一括でズームレベルを変更できませんが、Pythonコンソールを使えば一瞬で全項目の制限を解除できます。 【解決手順】 レイヤパネルで「国土地理院ベクター」レイヤを選択。 メニューの [プラグイン] ＞ [Pythonコンソール] を開く。 下部の入力欄に以下の1行を貼り付けてEnter。 Python layer = iface.activeLayer(); styles = layer.renderer().styles(); [s.setMaxZoomLevel( 22 ) for s in styles]; layer.renderer().setStyles(styles); layer.triggerRepaint() これでレベル22まで拡大しても地物が表示され続けます。 2. 巨大点群から傾斜区分図を一括作成する（岐阜県仕様） 数千万点のデータから傾斜区分図を作る際、QGISのGUI操作ではフリーズのリスクが高まります。ここでは、用途に応じて使い分けられる「超高速版」 と 「高精度フルバージョン」の2つのスクリプトを紹介します。 傾斜区分の定義（5段階） 値 1：0° - 2° （緑色） 値 2：2° - 15° （黄緑色） 値 3：15° - 30° （薄い黄色） 値 4：30° - 45° （山吹色） 値 5：45° - 90° （オレンジ色） A. 【超高速版】とりあえず全体を素早く確認したい時 1000万点超えのデータでも 数分 で完了します。微地形は多少簡略化されますが、スピード重視の現場確認に最適です。 Python import os import numpy as np import tkinter as tk from tkinter import filedialog try : from osgeo import gdal, ogr, osr except ImportError: ...