UAVと水域可視化処理による河川地形の新しい計測手法の開発

Size: px

Start display at page:

Download "UAVと水域可視化処理による河川地形の新しい計測手法の開発"

まなたかはし
5 years ago
Views:

平成 28 年 6 月 16 日インフライノベーション研究会第 31 回講演会 @ 東京大学ダイワユビキタス学術研究館 3F ホール UAV による河川地形河床材料モニタリング手法の検討岐阜大学流域圏科学研究センター原田守啓本研究は,

1 平成 28 年 6 月 16 日インフライノベーション研究会第 31 東京大学ダイワユビキタス学術研究館 3F ホール UAV による河川地形河床材料モニタリング手法の検討岐阜大学流域圏科学研究センター原田守啓本研究は, 国土交通省河川砂防技術研究開発助成地域課題 ( 平成 27 年度 UAVと水域可視化処理による河川地形の新しい計測手法の開発に関する研究 ) により実施した. 本研究の実施にあたり, 国土交通省木曽川上流河川事務所調査課の多大なるご協力を得た. 1

2 従来のモニタリング手法は多様化するニーズに応えられるか? 河川地形計測定期縦横断測量数年に 1 度,200m 間隔の横断測量が基本航空レーザー測量 (LP) 水中に弱点 ALB(Airborne Laser Bathymetry) セグメント 1 における砂州の動態, セグメント 2 における高水敷掘削後の地形変化 ( 土砂堆積 ) 等を把握するには時間解像度空間解像度ともに不十分. 植物の動態河川水辺の国政調査 5 年に 1 度, 空中写真判読から河川植生図樹林化に至るプロセス ( 安定した草本群落 or 樹木の実生の定着樹林化 ) を監視するには, やはり時間解像度空間解像度ともに不十分. 河床材料地整事務所によって異なる ( 不定期?) 縦断的分布の把握に重点. 地形を構成する粒径の幅が広く, 空間的なばらつきが大きい中上流域における出水前後の河床環境の変化を捉えるには甚だ不十分.( ニーズ形成も不十分 ) 2

3 セグメント 2 における高水敷掘削後の変化の例木曽川水系揖斐川の自然堤防帯 32~39kp (I b 1/3,300) 洪水流下能力の向上を目的とした河道掘削施工年度掘削高さの設定により 14 の工区が存在低く掘った箇所は赤系高めの箇所は青系 Hori et al

4 高水敷掘削後の景観変化の例 (D 地区 ) H13 施工掘削後 6 年 4 ワンドたまり掘削後 11 年

5 高水敷掘削後の景観変化の例 (O 地区 ) H17 施工掘削後 2 年ワンドたまり 5 掘削後 7 年

6 変動量 (m) 地形変化 ( 掘削後の堆積速度 ) H14,H17,H20 の定期縦横断測量結果から掘削部分の変動量を算出グループ 1 グループ 2 グループ 3 A B C D E L G M F O H I J N H12 H13 H13 H13 H14 H15 H16 H16 H17 H17 H18 H18 H19 H19 H14-H17 の変動量堆積速度 ( 参考値 )= 上記の変動量 / 経過年数 (3 年 ) 堆積速度 ( 参考値 )5~12cm/ 年やや右肩上がりか? 6 堤内地となっている旧河道では自然堤防 1.5cm/ 年, 後背湿地 0.1~0.2cm/ 年 Hori et al.(2011)

7 H13 H19 H24 H13 H19 H24 H13 H19 H24 H13 H19 H24 面積割合掘削後の植生遷移 H13,H19,H24 の河川植生図に基づく群落変化 100% 90% 80% 70% 60% 50% 40% 30% 20% 10% 0% 水域が保たれている草地ヤナギ群落裸地ヤナギ群落草地ヤナギ群落ヤナギ類多年生草本 1 年生草本裸地管理地開放水面渇水位低水位平水位豊水位 7 掘削高さ - 地被面積割合の変化

河川地形計測手法を比較した例ヘリ LP デジタル航空カメラ +SfM 1400m 300m 空間スケール別のコスト比較 ( リーチスケールの UAV-SfM のコストを 1 とした場合 ) 対象スケール航測技術水系セグメントリーチデジ航 50 10 2 ヘリLP 100 15 3 UAV 100 13 1 UAV-SfM 60m

8 河川地形計測手法を比較した例ヘリ LP デジタル航空カメラ +SfM 1400m 300m 空間スケール別のコスト比較 ( リーチスケールの UAV-SfM のコストを 1 とした場合 ) 対象スケール航測技術水系セグメントリーチデジ航ヘリLP UAV UAV-SfM 60m 従来の定期縦横断測量と比較して, 情報量は圧倒的に多い. 樹木, 草本等の植物が繁茂する区間, 水深が浅い場面では,UAV-SfM が有利. 要求される計測精度, 規模, 頻度に応じた使い分けを提案. 佐貫方城, 渡辺敏, 宮田真考, 草加大輝 : 3 種の航空測量技術を使用した河道地形の効率的測量の実装展開に向けた比較検討, 河川技術論文集,21, pp ,

9 SfM / MVS 技術 Structure from Motion / Multi-View Stereo 異なるアングルから撮影された複数の画像から, 被写体とカメラの相対的な位置関係を推定し, 被写体の三次元形状を再構成する技術対象は地形に限らない. ありとあらゆるものの三次元形状を復元可能. UAV 空撮画像でなくとも, 同じ被写体に対してラップを確保した複数の写真であれば, 普通のデジカメを手持ち撮影した画像でも良い. 使いやすいユーザーインターフェースを備えた市販ソフトから, フリーウェアまで. 個人的な印象では,Agisoft Photoscan が最もシェアが高い? 9

: The development of an automated correction procedure for digital photogrammetry for the study of wide,

10 UAV-SfM の河川地形計測への適用例 Westaway, R. M., Lane, S. N., Hicks, D. M.: The development of an automated correction procedure for digital photogrammetry for the study of wide, shallow, gravel-bed rivers, Earth Surface Processes and Landforms, 25, 2, pp , Javernick L, Brasington J, Caruso B.: Modeling the topography of shallow braided rivers using Structure-from-Motion photogrammetry. Geomorphology. 213, pp ,

11 本報のあらまし UAV-SfM による写真地形計測の適用範囲を水中部に拡大することにより, 河川管理に資する付加価値の高い調査手法を提供する. セグメント 1( 扇状地 ) 区間における現地調査を実施し, 画像処理による水中可視化による水中地形計測精度の向上について検討. UAV 空撮により得られた高精細な画像に基づき, 河床表層材料の粒度分布を推定する手順について検討. 電動小型 UAV これらの検討にあたって運用のノウハウを整理することに重点をおいて検討を進めました. 11

12 方法長良川扇状地区間における UAV 空撮と河床形状計測の実施 1 回目 : 増水時 2 回目 : 平水時水中可視化! アクアスコープ ( 朝日航洋 ) による処理流量の異なる条件で 2 度の UAV 空撮を実施水中河床形状の検証用データを別途 ADCP により取得 12

13 13

14 14

15 15

16 16

17 17

18 18

19 19

20 20

) 撮影写真ラップ率 60% 以上にて撮影未処理画像データセット処理済画像データセット 3 次元モデル ( 画像処理なし ) 画像処理の有無による水中河床形状計測精度の検証

21 UAV 空撮写真 * 画像処理なし * 水中可視化処理後アクアスコープ ( 朝日航洋 ) による処理 αuav (amuse-oneself 社製 ) 飛行高度は対地高度 150m ( 部分的に 50m,100m も.) 撮影写真ラップ率 60% 以上にて撮影未処理画像データセット処理済画像データセット 3 次元モデル ( 画像処理なし ) 画像処理の有無による水中河床形状計測精度の検証水深補正方法の検討 SfM ソフトウェア (Agisoft Photoscan) 河床高の推定値 3 次元モデル ( 画像処理あり ) 河床高の実測値 ADCP 超音波流速計 M9(Sontek 社製 ) + エンジン付きボート澪筋に沿って往復しながら河床形状を計測. ADCP 計測結果に基づく標高 ( 水深 ) 分布 21

22 標高 [m] 標高 [m] 標高 [m] 標高 [m] 標高 [m] 標高 [m] 結果 1 水中河床形状推定精度の検証河床形状が複雑な水制周辺に 6 本の検査断面を設定し断面図を作図 3 次元モデル 4 ケースと ADCP 実測値を比較陸上部はどのケースも概ね一致水中部にはケースによりばらつきがある Transect Transect Transect 水制左岸堤防天端からの距離 [m] Transect 左岸堤防天端からの距離 [m] 水制 Transect 左岸堤防天端からの距離増水時 [m] ( 画像処理なし ) 増水時 ( 画像処理あり ) 平水時 ( 画像処理なし ) 平水時 ( 画像処理あり ) ADCP Transect2 Transect4 Transect6 9 Transect Transect6 6 6 水制左岸堤防天端からの距離 [m] 左岸堤防天端からの距離 [m] 左岸堤防天端からの距離 [m] 22

23 3 次元モデルから求めた見かけの水深結果 2 水中河床形状推定精度の検証 ( 平水時 ) 画像処理なし画像処理あり空撮写真から推定された河床地形は見かけの水深を光の屈折率である 1.33 倍するだけの処理で河床高を表現しうる処理の定式化が可能. ADCP 計測結果の水深コンディションがよければ水深 2m 程度までは計測可能水中可視化処理によってより水深が大きい領域まで計測が可能となる 23

24 3 次元モデルから求めた見かけの水深結果 2 水中河床形状推定精度の検証 ( 増水時 ) 画像処理なし画像処理あり ADCP 計測結果の水深水深の影響よりも増水による濁り水面の擾乱による反射といった撮影コンディションのほうが支配的な影響を及ぼしている現場での運用では, 空撮の実施適否の判断が重要になる. 27

25 UAV-SfM による河川地形計測の一般化水域可視化処理 UAV-SfMの実施オルソ画像 DEMモデルマスク抽出水部マスクの作成水部 DEM 水位縦断形の作成見かけの水深分布 Apparent Water Depth 左右岸の水際線と DEM から水位縦断形を推定見かけの水深 AWD= 水位 -DEM 1.33 倍するだけ. + ノイズ処理? 補正された水深分布 Corrected Water Depth 補正 DEM= 水位 - 補正された水深 CWD 陸部 DEM 補正された水部 DEM 合成あるいは上書き水深補正された DEM モデル 28

0 本ケースにおいて検出限界と予想した粒径 (40mm) 面積格子法画像解析法 1 10 100 1000 diameter[mm] 撮影高度 50m の UAV 空撮画像

26 結果 3 河床材料調査への活用の試行 UAV 空撮画像を用いて砂州表層河床材料の粒度分布の推定を試みた本検討では現場技術者が河床材料調査を目的としたUAV 空撮を行うための運用のノウハウを整理することを目的に撮影及び画像解析を試行した本ケースにおいて検出限界と予想した粒径 (40mm) 面積格子法画像解析法 diameter[mm] 撮影高度 50m の UAV 空撮画像輪郭抽出法による画像解析例 (BASEGRAIN, スイス工科大 ) 試行の結果撮影に先立って画像解析で検出したい粒子径と撮影機材の性能に基づいて飛行高度等の計画を立案することが可能 29

27 空撮による表層河床材料調査あらかじめ現地踏査を行い, 現地に見られる河床材料径の幅を大まかに把握して, 目標精度を決定. 撮像素子解像度 =4592[pixel] 撮像素子サイズ =17.3[mm] 焦点距離 =14[mm] 画像解析により判別可能な粒子径の目安を, 撮影機材と飛行高度から計算. 対地高度 =50[m] 目標とする粒子の検出精度に合わせて,UAV の撮影高度を逆算することにより, 効率的なフライト計画を立案できる. 撮影領域幅 = 撮像素子サイズ x ( 対地高度 / 焦点距離 )=61.8[m] 空間分解能 = 撮影領域幅 / 撮像素子解像度 = 13.5[mm] 粒子の輪郭抽出のため, 最低 3pixel 必要とすれば本ケースでは,40mm が画像計測の下限と想定された. 30

28 3. 結果と考察 Results and discussion / 粒度分布 GSD 砂州における河床材料分布の実態下流上流 1つの砂州の中でも, 面的な分級が生じており, 場所によって粒径が異なる. さらに, 表層と準表層は粒度分布が相当傾向が異なる.( 表面は粗粒化している ) 出水により, 地形変化を生じると同時に粒度分布の変化も伴う. 31

29 UAV によるモニタリング手法は多様化するニーズに応えられるか? 河川地形計測定期縦横断測量数年に 1 度,200m 間隔の横断測量が基本航空レーザー測量 (LP) 水中に弱点 ALB(Airborne Laser Bathymetry) UAV-SfM は, 河川地形の変化を高精度に計測しうる. 水中部は撮影条件の影響を大きく受けるものの, 陸部の地形変化のモニタリングには十分供しうる. 植物の動態河川水辺の国政調査 5 年に 1 度, 空中写真判読から河川植生図高解像度な画像を, 植物の生活環に合わせた頻度で撮影できれば, モニタリング手法として機能しうる ( 未検証 ). 高水敷掘削の場合, 掘削後 1,2 年が勝負河床材料地整事務所によって異なる ( 不定期?) 縦断的分布の把握に重点. 表層粒度分布のうち, ある程度の大きさの粒子は画像解析で捉えうる. しかし, 準表層, 水面下の粒度分布を捉えずして, 土砂動態は把握しきれない. 32

30 一般的なニーズシーズの関係土木分野では, 伝統的に, 行政ニーズと技術シーズが一体となって, 研究開発がすすめられてきた ( ように思われる ). ニーズの発信データの提供試験施工等河川管理者のニーズノウハウの蓄積技術者研究者の技術シーズ研究成果の発表, 新技術の提案提供 33

31 計測技術をめぐるニーズシーズ考他分野において発展した汎用的計測技術 ( 例 :UAV-SfM) は, 河川管理者のニーズだけでなく, 研究者のニーズも満たしうる. ニーズはやや異なるものの, 河川管理者研究者全てがユーザーである. フィールドでの協働によるノウハウの開発と蓄積, 共有へこれに応える民間ニーズも含む河川管理者のニーズ現場での運用のノウハウ ( 計測精度, コスト, 調査計画手法, 実施可否の判断 ) 汎用的計測技術例 :UAV-SfM フィールドでの協働ノウハウの共有研究者のニーズ既存技術を土台とした新しい計測手法の開発痒いところに手が届くモニタリング技術 [ 計測技術を研究対象としない多くの研究者 ] [ とくに現場技術者 ] 多様化した例 : 河道掘削後の地形変化, 植物の侵入 [ 計測技術が研究対象の研究者 ] 研究対象とする現象の観測データ ( 河川地形, 河床材料, 植物等 ) 数値計算の初期値, 検証データ現地調査の効率化( 測量の代替手段等 ) 例 : 石礫床河川の河床環境の数値計算モデル 34

資料 -5 第 5 回岩木川魚がすみやすい川づくり検討委員会現地説明資料平成 28 年 12 月 2 日東北地方整備局青森河川国道事務所

資料 -5 第 5 回岩木川魚がすみやすい川づくり検討委員会現地説明資料平成 28 年月 2 日東北地方整備局青森河川国道事務所現地説明資料富士見橋経年変化富士見橋は 51.8k 付近に H7~H22 の河川水辺の国勢調査で早瀬が確認しており H5~ で近傍で最深河床高の低下したことで平水流量時の水深が 0.2~0.4m の浅場 ( 瀬 ) が減少したと推定されるがその後も早瀬が確認されている