JP2007011582A - 洪水予測システム - Google Patents

Abstract

【課題】 従来の河川の流域情報は、国土地理院が発行する数値地図や現地測量を行って作成しているが、これらの情報を作成するためには膨大な人手と時間及び費用を要するのが現状である。また、時間の経過に伴い、気象、気候又は自然災害により、流域を形成する地形、地質又は地表面被覆植生が変化し、実態を把握することが困難となる。
【解決手段】 地球観測衛星の観測データであるリモートセンシングデータを高度解析処理する手段により、現状の問題点と困難を解決することが可能となる。その結果、短期間に最新の流域情報を安価に得ることが可能となる。また、定期的にかつ広域に情報を更新することができると共に、災害時における緊急情報の取得も可能となる。
【選択図】 図２

Description

本発明は、地球観測衛星からのリモートセンシングデータを利用し、豪雨時の河川、山岳部及び平野部における水の移動を計算することにより、洪水災害のシミュレーションを行う洪水予測システムに関するものである。

地球温暖化による異常気象に起因する集中豪雨は、大規模な洪水災害をもたらし、人的及び経済的損失も増大している。洪水被害の軽減には、信頼性の高い洪水予測技術の確立が不可欠である。

広域に渡る面的及び線的情報の収集方法として、地球観測衛星からのリモートセンシングデータの利用が考えられ、洪水予測技術に必要とされる調査作業を地球観測衛星からのリモートセンシングデータによって代替することにより、コストや期間が大幅に軽減されうる。

特許文献１に記載されているように、異なる位置にある衛星のセンサから所定の地域を地球観測衛星からのリモートセンシングすることにより得られたデジタル画像データを入力し、所定の地域の標高分布を示すデジタル標高データ、水害の起点位置データ及び流量データを算出して、水害が起こった場合の浸水状況をシミュレーションする水害シミュレーション方法及びシステムという発明も公開されている。
特開２００４−３４０７４３号公報

しかしながら、洪水予測に必要となる情報の中には、収集が困難な情報も多く、専門家による長期間の調査には膨大なコストも必要なことから、洪水予測技術の確立が遅れ、有効な洪水予防ができない地域が存在する。

そこで、本発明は、地球観測衛星からのリモートセンシングデータを利用し、豪雨時の河川、山岳部及び平野部における水の移動を計算することにより、洪水災害のシミュレーションを行う洪水予測システムを提供することを目的とするものである。

本発明は、上記の課題を解決するために、衛星画像データ及びＧＰＳ・ＧＩＳデータを使用して、海領域図、河川図、斜面方位図、水系図及び流域図のデータを作成する流域の抽出と、前記衛星画像データを使用して、土地被覆分布図のデータを作成する分布図の作成・編集と、前記流域の抽出及び分布図の作成・編集で作成したデータ並びに気象データを使用して、河川の上流域についてコンピュータで流出状況を計算する流域流出モデルによるシミュレーションと、前記ＧＰＳ・ＧＩＳデータを使用して、河川の分断地点であるノード、及び分断された河川であるセグメントごとにデータを作成するノード・セグメントの作成と、前記ＧＰＳ・ＧＩＳデータを使用して、単位領域であるセルごとにデータを作成するセルの作成と、前記流域流出モデルによるシミュレーション、ノード・セグメントの作成及びセルの作成で作成したデータを使用して、河川の下流域についてコンピュータで氾濫状況を計算する河道水理モデルによるシミュレーションと、前記河道水理モデルによるシミュレーションで作成したデータを編集して、コンピュータの画面に表示するシミュレーション結果の表示とからなることを特徴とする洪水予測システムの構成とした。

本発明は、以上の構成であるから以下の効果が得られる。第１に、リモートセンシングデータを利用して洪水災害のシミュレーションを行うことにより、リアルタイムに洪水を予測できる上、コストや期間も大幅に軽減することができる。

第２に、上流域を対象とする流域流出モデルと、下流域を対象とする河道水理モデルに分かれており、一度に全体を動かしたり、個別に動かしたりと、利用目的に応じたシミュレーションが可能となる。

本発明は、地球観測衛星からのリモートセンシングデータを利用し、豪雨時の河川、山岳部及び平野部における水の移動を計算することにより、洪水災害のシミュレーションを行うという目的を、衛星画像データ及びＧＰＳ・ＧＩＳデータを使用して、海領域図、河川図、斜面方位図、水系図及び流域図のデータを作成する流域の抽出と、前記衛星画像データを使用して、土地被覆分布図のデータを作成する分布図の作成・編集と、前記流域の抽出及び分布図の作成・編集で作成したデータ並びに気象データを使用して、河川の上流域についてコンピュータで流出状況を計算する流域流出モデルによるシミュレーションと、前記ＧＰＳ・ＧＩＳデータを使用して、河川の分断地点であるノード、及び分断された河川であるセグメントごとにデータを作成するノード・セグメントの作成と、前記ＧＰＳ・ＧＩＳデータを使用して、単位領域であるセルごとにデータを作成するセルの作成と、前記流域流出モデルによるシミュレーション、ノード・セグメントの作成及びセルの作成で作成したデータを使用して、河川の下流域についてコンピュータで氾濫状況を計算する河道水理モデルによるシミュレーションと、前記河道水理モデルによるシミュレーションで作成したデータを編集して、コンピュータの画面に表示するシミュレーション結果の表示とからなる構成とすることで実現した。

以下に、添付図面に基づいて、本発明である洪水予測システムについて詳細に説明する。

洪水予測システム１は、衛星画像データ２及びＧＰＳ・ＧＩＳデータ３を使用して、海領域図、河川図、斜面方位図、水系図及び流域図のデータを作成する流域の抽出６と、前記衛星画像データ２を使用して、土地被覆分布図のデータを作成する分布図の作成・編集７と、前記流域の抽出６及び分布図の作成・編集７で作成したデータ並びに気象データ４を使用して、河川の上流域についてコンピュータで流出状況を計算する流域流出モデルによるシミュレーション８と、前記ＧＰＳ・ＧＩＳデータ３を使用して、河川の分断地点であるノード、及び分断された河川であるセグメントごとにデータを作成するノード・セグメントの作成９と、前記ＧＰＳ・ＧＩＳデータ３を使用して、単位領域であるセルごとにデータを作成するセルの作成１０と、前記流域流出モデルによるシミュレーション８、ノード・セグメントの作成９及びセルの作成１０で作成したデータを使用して、河川の下流域についてコンピュータで氾濫状況を計算する河道水理モデルによるシミュレーション１１と、前記河道水理モデルによるシミュレーション１１で作成したデータを編集して、コンピュータの画面に表示するシミュレーション結果の表示１２とからなることを特徴とする。

図１は、本発明である洪水予測システムの構成を示すブロック図である。洪水予測システム１は、コンピュータ５において実行されるプログラムであり、衛星画像データ２、ＧＰＳ・ＧＩＳデータ３及び気象データ４の外部データを使用する。

コンピュータ５は、少なくとも入力装置５ａ、記憶装置５ｂ、出力装置５ｃ及び中央処理装置５ｄからなり、必要に応じて、描画機能をサポートする装置などを追加することができる。

入力装置５ａは、コンピュータ５にデータを入力したり、プログラムの実行指令を与える装置であり、マウス、キーボード又はタッチパネル等が該当する。尚、入力されたデータ等は、一時的に記憶装置５ｂに記憶される。

記憶装置５ｂは、作業領域となる主記憶装置（メモリ）と、保存領域となる補助記憶装置とからなる。補助記憶装置には、ハードディスクやデータベース等があるが、フレキシブルディスクやコンパクトディスク等の記録媒体も含めるものとする。

洪水予測システム１のプログラムは、補助記憶装置に保存されており、入力装置５ａから実行指令が与えられると、主記憶装置に読み込まれ実行される。尚、プログラムは、処理又は機能ごとに部品化することもできる。

出力装置５ｃは、記憶装置５ｂに記憶されているデータ等を出力する装置であり、ディスプレイ（モニタ）又はプリンタ等が該当する。尚、データは、グラフ又は画像などに編集して出力することもできる。

中央処理装置（ＣＰＵ）５ｄは、コンピュータ５の中枢となる装置であり、プログラムの実行を行う。実行に際しては、プログラムを解析し、データに対して計算を行う演算命令や、入力装置５ａ、記憶装置５ｂ及び出力装置５ｃに対する制御命令を出す。

衛星画像データ２は、人工衛星のリモートセンシング技術を利用して取得したデータであり、地質調査や資源探索に加えて、温暖化、砂漠化、オゾン層破壊、酸性雨又は鉱物研究等の地球規模の環境変化を把握することができる。

ＧＰＳ・ＧＩＳデータ３は、ＧＰＳ（全地球無線測位システム）とＧＩＳ（地理情報システム）を組み合わせたデータであり、デジタル化された地図データに統計データや位置の持つ属性情報が付加される。

尚、ＧＰＳとは、人工衛星から発射した電波の到達時間などから地球上の電波受信者の位置を３次元測位するシステムであり、ＧＩＳとは、地理的な様々な情報に関連付けなどの処理を行い，データ化された地図として視覚的に表示するシステムである。

気象データ４は、洪水予測システム１の対象となる地域の気象観測所又は気象衛星などから取得したデータである。過去に観測したデータであるが、予測したデータを用いることもできる。

衛星画像データ２、ＧＰＳ・ＧＩＳデータ３及び気象データ４の外部データは、コンピュータ５の記憶装置５ｂに取り込まれる。外部データの取込みには、記録媒体を利用したり、インターネット又は専用通信回線などを利用する。尚、有線だけでなく無線を利用することもできる。

図２は、本発明である洪水予測システムの流れを示すフローチャートである。洪水予測システム１は、流域の抽出６、分布図の作成・編集７、流域流出モデルによるシミュレーション８、ノード・セグメントの作成９、セルの作成１０、河道水理モデルによるシミュレーション１１、及びシミュレーション結果の表示１２の処理からなる。

流域の抽出６は、衛星画像データ２及びＧＰＳ・ＧＩＳデータ３から抽出した情報を基に、海領域図、河川図、斜面方位図、水系図及び流域図についてのデータを作成する。分布図の作成・編集７は、衛星画像データ２から抽出した情報を基に、土地被覆分布図についてのデータを作成する。

流域流出モデルによるシミュレーション８は、気象データ４から抽出した情報、流域の抽出６で作成したデータ、及び分布図の作成・編集７で作成したデータを基に、流域流出モデルを作成し、コンピュータ５を用いて模擬実験を行うことで洪水を予測する。

ノード・セグメントの作成９は、ＧＰＳ・ＧＩＳデータ３から抽出した情報、及び流域の抽出６で作成したデータを基に、河道区間図及び河道の地理についてのデータを作成する。セルの作成１０は、衛星画像データ２及びＧＰＳ・ＧＩＳデータ３から抽出した情報を基に、セル図及びセルの地理についてのデータを作成する。

河道水理モデルによるシミュレーション１１は、気象データ４から抽出した情報、並びに流域流出モデルによるシミュレーション８、ノード・セグメントの作成９及びセルの作成１０で作成したデータを基に、河道水理モデルを作成し、コンピュータ５を用いて模擬実験を行うことで洪水を予測する。

シミュレーション結果の表示１２は、衛星画像データ２及び気象データ４から抽出した情報、並びにセルの作成１０、流域流出モデルによるシミュレーション８及び河道水理モデルによるシミュレーション１１で作成したデータを基に、洪水予測の結果をコンピュータ５で編集して画面に表示する。

図３は、本発明である洪水予測システムの処理の流れを示す図である。洪水予測システム１の処理うち、流域の抽出６、分布図の作成・編集７、ノード・セグメントの作成９及びセルの作成１０はデータ抽出・モデル化処理１３としてまとめることができ、流域流出モデルによるシミュレーション８及び河道水理モデルによるシミュレーション１１は洪水予測シミュレーション１４としてまとめることができる。

データ抽出・モデル化処理１３は、洪水予測シミュレーション１４を行うために必要となるデータを取得及び変換する手段であり、洪水予測シミュレーション１４は、洪水予測を行う手段であり、シミュレーション結果の表示１２は、洪水予測シミュレーション１４の結果を編集して可視化する手段である。

洪水予測システム１の各処理は、更に機能ごとに部品化することができる。各機能は、外部とやり取りをするインタフェースとして入力パラメータと出力パラメータを設けることができ、外部からデータを受け取り、演算等した後、外部に結果のデータを渡すことができる。

尚、各機能の形態はそれぞれ異なり、例えば、データベースからデータを入力して演算した結果のデータを別のデータベースに出力したり、入力されたデータを編集して出力装置５ｃに表示したり、コンピュータ５内の装置を制御する等である。

流域の抽出６は、海領域図の作成６ａ、河川図の作成６ｂ、斜面方位図の作成６ｃ、水系図の作成６ｄ、及び流域図の作成６ｅの機能からなる。

海領域図の作成６ａは、衛星画像データ２のうち反射率データ２ｂを読み込み、反射率の違いから開水面を検出し、画面上に描画及び画像を海領域図データ６ｆとして保存する。尚、開水面とは、衛星の観測している上空に対し、遮蔽のない水域のことである。また、水域とは、湖沼、海又は河川などの多量の水が定常的に存在する領域である。

河川図の作成６ｂは、手動操作により洪水予測システム１に認識させる河川を地図上に定義し、河川図データ６ｇを作成する。ＧＰＳ・ＧＩＳデータ３の河川情報３ｄ又は海領域図の作成６ａで作成した海領域図データ６ｆを参考にして地図上に手動で河川を定義していくが、河川の終点は、必ず海領域図データ６ｆの海又は画面端に至るようにする。

斜面方位図の作成６ｃは、衛星画像データ２のうち数値標高モデルデータ２ａから斜面方位図データ６ｈを自動作成する。数値標高モデルデータ２ａに対し、河川図の作成６ｂで作成された河川図データ６ｇを用いて高度補正を施すことにより、水理解析に適した斜面方位図データ６ｈを得ることができる。

水系図の作成６ｄは、斜面方位図の作成６ｃで作成された斜面方位図データ６ｈに基づいて地図を構成する各点からの落水線を計算し、落水線が高頻度に通過する点を水系として自動抽出することにより水系図データ６ｉを作成する。尚、落水線とは、傾斜方向を結んでいくことによって出来る線であり、降水の流れていく経路とみなされる。また、水系とは、降水時に多量の水が通過する経路のことである。

流域図の作成６ｅは、水系図の作成６ｄで作成された水系図データ６ｉに示される水系上のある点について、その点に至る全ての落水線を斜面方位図の作成６ｃで作成された斜面方位図データ６ｈに基づいて計算し、結合することにより、その点を排水先とする流域を自動抽出し、流域図データ６ｊを作成する。尚、流域とは、降水がある点へと集まる範囲のことである。

分布図の作成・編集７は、土地被覆分布図の作成７ａ、土壌分布図の編集７ｂ、パラメータテーブルの編集７ｃ、及びパラメータテーブルの編集７ｄの機能からなる。

土地被覆分布図の作成７ａは、衛星画像データ２のうち反射率データ２ｂを読み込み、反射率の違いからトレーニングデータに基づいて土地被覆の分類を行い、土地被覆分布図データ７ｅを作成する。

尚、トレーニングデータとは、土地被覆の分類の指標となる反射率であり、対象の土地の反射率が指標のうちどれに最も近いかによって分類する。また、土地被覆とは、植生、舗装面又は水域など土地の表面の状態である。

土壌分布図の編集７ｂは、ＧＰＳ・ＧＩＳデータ３のうち土壌分布図データ３ａの内容を編集する。パラメータテーブルの編集７ｃは、土地被覆分布図データ７ｅに対し、地表面粗度や貯留損失量などの被覆パラメータ３ｂを与え、内容を編集する。パラメータテーブルの編集７ｄは、土壌分布図データ３ａに対し、毛管吸上げ、透水係数、土湿不足などの土壌パラメータ３ｃを与え、内容を編集する。

流域流出モデルによるシミュレーション８は、計算制御ファイルの編集８ａ、流域特性計算８ｂ、及び流域流出シミュレーション８ｃの機能からなる。計算制御ファイルの編集８ａは、シミュレーション開始日時、シミュレーション期間、計算単位時間などの流域流出モデルによるシミュレーション８に関する設定を変更する。

流域特性計算８ｂは、以下の３つの計算を行い、上流域の地理データ８ｄを作成する。第１に、流域図の作成６ｅで作成された流域図データ６ｊの各流域に対して、衛星画像データ２の数値標高モデルデータ２ａから面積、平均勾配及び幅を計算する。

第２に、ＧＰＳ・ＧＩＳデータ３の被覆パラメータ３ｂと土地被覆分布図の作成７ａで作成した土地被覆分布図７ｅから、土地被覆分布の面積で被覆パラメータの平均を取ることにより初期損失雨量及び粗度を計算する。

第３に、ＧＰＳ・ＧＩＳデータ３の土壌分布図データ３ａと土壌パラメータ３ｃから、土壌分布の面積で土壌パラメータの平均を取ることにより初期土湿不足、毛管吸上げ及び水理伝導度を自動計算する。

流域流出シミュレーション８ｃは、気象データ４の時系列降水量データ４ｄ及び月平均蒸発量データ４ｅ並びに流域特性計算８ｂで作成された上流域の地理データ８ｄを用いて、上流域の地理データ８ｄの各流域の流出量である時系列流出量データ８ｅを自動算出する。尚、流出量とは、流域の排水点から流域外へ排出される水の量である。

ノード・セグメントの作成９は、河道分割９ａ、及び河道特性入力９ｂの機能からなる。河道分割９ａは、流域図の作成６ｅで作成された河川図データ６ｇの河川をＧＰＳ・ＧＩＳデータ３の河川情報３ｄの計測地点位置に基づいて分断し、河道区間図データ９ｃを作成する。

河道特性入力９ｂは、河道分割９ａで作成した河道区間図データ９ｃの各ノードと各セグメントに対し、ＧＰＳ・ＧＩＳデータ３の河川情報３ｄを基にパラメータを手動入力し、河道の地理データ９ｄを作成する。尚、ノードとは、河川の分断地点のことであり、セグメントとは、分断された河川、即ち２つのノードに挟まれた河道の区間のことである。

セルの作成１０は、氾濫原分割１０ａ、及びセル特性入力１０ｂの機能からなる。尚、セルとは、浸水状況をシミュレートする領域の最小単位のことであり、セル単位にデータを作成する。

氾濫原分割１０ａは、衛星画像データ２の反射率データ２ｂ及びＧＰＳ・ＧＩＳデータ３の氾濫原情報３ｅを参考にして、浸水しうる領域を、セル単位に手動で入力し、セル図１１ｃを作成する。

セル特性入力１０ｂは、氾濫原分割１０ａで作成したセル図データ１０ｃの各セルに対し、ＧＰＳ・ＧＩＳデータ３の氾濫原情報３ｅを基にパラメータを手動入力し、セルの地理データ１０ｄを作成する。

河道水理モデルによるシミュレーション１１は、計算制御ファイルの編集１１ａ、河川・セル結合１１ｂ、及び河道水理シミュレーション１１ｃの機能からなる。計算制御ファイルの編集１１ａは、シミュレーション開始日時、シミュレーション期間、計算単位時間などの河道水理モデルによるシミュレーション１１に関する設定を変更する。

河川・セル結合１１ｂは、河道特性入力９ｂで河川情報３ｄを基に作成した河道の地理データ９ｄのノード及びセグメント群と、セル特性入力１０ｂで氾濫原情報３ｅを基に作成したセルの地理データ１０ｄのセル群との間に、水を相互交換する関係を定義し、水理的に結合し、下流域・河川の地理データ１１ｄを作成する。

河道水理シミュレーション１１ｃは、気象データ４の時系列降水量データ４ｄ、時系列海面高度データ４ａ、初期水位データ４ｂ及び初期流量データ４ｃ、並びに河川・セル結合１１ｂで作成した下流域・河川の地理データ１１ｄ及び流域流出シミュレーション８ｃで作成した時系列流出量データ８ｅを用いて、下流域・河川の地理データ１１ｄの各セルの水位である時系列水位データ１１ｅを自動算出する。

シミュレーション結果の表示１２は、２次元洪水描画１２ａ、３次元洪水描画１２ｂ、及びハイエト−ハイドログラフ描画１２ｃの機能からなる。

２次元洪水描画１２ａは、氾濫原分割１０ａで作成したセル図データ１０ｃのセルの範囲に、河道水理シミュレーション１１ｃで作成した時系列水位データ１１ｅの水位を色分けして２次元洪水マップデータ１２ｄを作成し表示する。尚、時間変化を追って表示させることも可能である。

３次元洪水描画１２ｂは、２次元洪水描画１２ａで作成した平面の地図画像に衛星画像データ２のうち数値標高モデルデータ２ａの標高データを組み合わせて３次元洪水マップデータ１２ｅを作成し、３次元に立体化して表示する。

ハイエト−ハイドログラフ描画１２ｃは、気象データ４のうち時系列降水量データ４ｄと流域流出シミュレーション８ｃで作成した時系列流出量データ８ｅを用いて、ハイエト−ハイドログラフデータ１２ｆを作成し、ハイエト−ハイドログラフを描画する。尚、ハイエト−ハイドログラフとは、降水量と流出量の時間変化を表したグラフである。

図４は、本発明である洪水予測システムのデータの流れを示す図である。衛星画像データ２、ＧＰＳ・ＧＩＳデータ３及び気象データ４は、洪水予測システム１を実行するに際して予め用意される外部データであり、その他に、洪水予測システム１の機能を実行することにより作成される内部データがある。外部データ及び内部データは、主としてデータベースである。

データベースは、データベース管理システムによって管理され、洪水予測システム１の継続的な運用におけるデータ管理の煩雑さを軽減し、毎回のシミュレーションの実施手順を簡便にすることができる。

衛星画像データ２は、数値標高モデルデータ２ａ、及び反射率データ２ｂからなり、いずれも地球観測衛星に搭載されたセンサにより、洪水予測システム１のシミュレーション対象地域について観測されたデータである。

数値標高モデルデータ２ａは、ピクセルごとに高度情報が定められた経緯度付きの地図データであり、斜面方位図の作成６ｃ、流域特性計算８ｂ及び３次元洪水描画１２ｂの機能の入力データとなる。

反射率データ２ｂは、ピクセルごとに放射計で記録された地表面の電波反射率が定められた経緯度付きの地図データであり、海領域図の作成６ａ、土地被覆分布図の作成７ａ及び氾濫原分割１０ａの機能の入力データとなる。

ＧＰＳ・ＧＩＳデータ３は、土壌分布図データ３ａ、被覆パラメータ３ｂ、土壌パラメータ３ｃ、河川情報３ｄ、及び氾濫原情報３ｅからなる。土壌分布図データ３ａは、洪水予測システム１のシミュレーション対象地域の土壌分布の情報を示す地図データであり、流域特性計算８ｂの機能の入力データとなる。

被覆パラメータ３ｂは、初期損失雨量と粗度について、洪水予測システム１のシミュレーション対象地域に分布する植生の種類ごとの数値をまとめたデータであり、流域特性計算８ｂの機能の入力データとなる。尚、初期損失雨量とは、降水のうち植生に遮られ地表面に到達しない雨量のことであり、粗度とは、植生による地表面流の流れにくさを示す指標である。

土壌パラメータ３ｃは、初期土湿不足、毛管吸上げ、水理伝導度について、洪水予測システム１のシミュレーション対象地域に分布する土壌の種類ごとの数値をまとめたデータであり、流域特性計算８ｂの機能の入力データとなる。尚、土壌パラメータ３ｃにおける土壌分布は、土壌分布図データ３ａと同一のものを使用する。

土壌パラメータ３ｃにおいて、初期土湿不足とは、平時における土の空気含有率、即ち飽和時の水分含有量に対する平時の空気含有量の割合であり、毛管吸上げとは、土が毛管作用により水を吸い上げる強さであり、水理伝導度とは、土の湿度による浸透能変化の度合いの係数である。

河川情報３ｄは、河川の流れる位置と形状、河川の水位ごとの垂直断面積と河床粗度及びその計測地点を示すデータであり、河川図の作成６ｂ、河道分割９ａ、河道特性入力９ｂ及び河川・セル結合１１ｂの機能の入力データとなる。

氾濫原情報３ｅは、過去の洪水実績により定められる最大浸水範囲、浸水状況及びその範囲の詳細な地形情報であり、氾濫原分割１０ａ、セル特性入力１０ｂ及び河川・セル結合１１ｂの機能の入力データとなる。

気象データ４は、時系列海面高度データ４ａ、初期水位データ４ｂ、初期流量データ４ｃ、時系列降水量データ４ｄ、及び月平均蒸発量データ４ｅからなる。

時系列海面高度データ４ａは、洪水予測システム１のシミュレーション対象となる期間の海面の高さの時系列データであり、河道水理シミュレーション１１ｃの機能の入力データとなる。

初期水位データ４ｂは、洪水予測システム１のシミュレーション開始時における河川及び氾濫原の水位のデータであり、河道水理シミュレーション１１ｃの機能の入力データとなる。

初期流量データ４ｃは、洪水予測システム１のシミュレーション開始時における河川及び氾濫原の流量のデータであり、河道水理シミュレーション１１ｃの機能の入力データとなる。

時系列降水量データ４ｄは、洪水予測システム１のシミュレーション対象となる期間の降水量の時系列データであり、流域流出シミュレーション８ｃ、河道水理シミュレーション１１ｃ及びハイエト−ハイドログラフ描画１２ｃの機能の入力データとなる。

月平均蒸発量データ４ｅは、洪水予測システム１のシミュレーション対象となる月の水の蒸発量の平均値のデータであり、流域流出シミュレーション８ｃの機能の入力データとなる。

データ抽出・モデル化処理１３において作成される内部データは、海領域図データ６ｆ、河川図データ６ｇ、斜面方位図データ６ｈ、水系図データ６ｉ、流域図データ６ｊ、土地被覆分布図データ７ｅ、河道区間図データ９ｃ、及びセル図データ１０ｃである。

海領域図データ６ｆは、海領域図の作成６ａによって出力される水域画像ファイルのうち、手動操作によって海の範囲のみを残した画像であり、河川図の作成６ｂの入力データとなる。

河川図データ６ｇは、洪水予測システム１に認識させる河川を定義したデータであり、斜面方位図の作成６ｃ及び河道分割９ａの入力データとなる。尚、河川図データ６ｇを画像として編集し画面に描画すると、図６に示す通りになる。

斜面方位図データ６ｈは、地図を構成する各点において、点が周囲８方向のうちどの方向に最も大きく傾斜しているかを色分けしたデータであり、水系図の作成６ｄ及び流域図の作成６ｅの入力データとなる。尚、斜面方位図データ６ｈを画像として編集し画面に描画すると、図７に示す通りになる。

水系図データ６ｉは、洪水予測システム１の対象地域中の全ての水系を示すデータであり、流域図の作成６ｅの入力データとなる。尚、水系図データ６ｉを画像として編集し画面に描画すると、図８に示す通りになる。

流域図データ６ｊは、洪水予測システム１に認識させる流域を定義したデータであり、流域特性計算８ｂの入力データとなる。尚、流域図データ６ｊを画像として編集し画面に描画すると、図９に示す通りになる。

土地被覆分布図データ７ｅは、洪水予測システム１の対象地域中の全ての点について、土地被覆を分類したデータであり、流域特性計算８ｂの入力データとなる。尚、土地被覆分布図データ７ｅを画像として編集し画面に描画すると、図１０に示す通りになる。

河道区間図データ９ｃは、洪水予測システム１に認識させる河川と、河川の分断地点を定義するデータであり、河道特性入力９ｂの入力データとなる。河道の地理データ９ｄは、各ノード及び各セグメントのパラメータを定義したデータであり、河川・セル結合１１ｂの入力データとなる。

セル図データ１０ｃは、洪水予測システム１に認識させるセルを定義したデータであり、セル特性入力１０ｂ及び２次元洪水描画１２ａの入力データとなる。尚、セル図データ１０ｃを画像として編集し画面に描画すると、図１１に示す通りになる。セルの地理データ１０ｄは、各セルのパラメータを定義したデータであり、河川・セル結合１１ｂの入力データとなる。

洪水予測シミュレーション１４において作成される内部データは、上流域の地理データ８ｄ、時系列流出量データ８ｅ、下流域・河川の地理データ１１ｄ及び時系列水位データ１１ｅである。

上流域の地理データ８ｄは、各流域ごとに、面積、勾配、幅、不浸透域の割合、不浸透域の初期損失雨量、不浸透域の粗度、浸透域の初期損失雨量、浸透域の粗度、初期土湿不足、毛管吸上げ、及び水理伝導度を定めたデータであり、流域流出シミュレーション８ｃの入力データとなる。尚、不浸透域とは、舗装などにより水の浸透が起きない領域のことである。

時系列流出量データ８ｅは、上流域の地理データ８ｄの各流域における気象条件が気象データ４の時系列降水量データ４ｄ及び月平均蒸発量データ４ｅの場合の流出量の時系列データであり、河道水理シミュレーション１１ｃ及びハイエト−ハイドログラフ描画１２ｃの入力データとなる。

下流域・河川の地理データ１１ｄは、各ノード、各セグメント及び各セルごとのパラメータと、各ノード、各セグメント及び各セル間の水理的関係性を定義したデータであり、河道水理シミュレーション１１ｃの入力データとなる。

時系列水位データ１１ｅは、下流域・河川の地理データ１１ｄの各セルにおける気象条件が気象データ４の時系列降水量データ４ｄ、時系列海面高度データ４ａ、初期水位データ４ｂ及び初期流量データ４ｃ、並びに時系列流出量データ８ｅの場合の水位の時系列データであり、２次元洪水描画１２ａの入力データとなる。

シミュレーション結果の表示１２において作成される内部データは、２次元洪水マップデータ１２ｄ、３次元洪水マップデータ１２ｅ及びハイエト−ハイドログラフデータ１２ｆである。

２次元洪水マップデータ１２ｄは、時系列水位データ１１ｅの内容を地図上に色分けしたデータであり、３次元洪水描画１２ｂの入力データとなる。尚、２次元洪水マップデータ１２ｄを画像として編集し画面に描画すると、図１３に示す通りになる。

３次元洪水マップデータ１２ｅは、２次元洪水マップデータ１２ｄの内容を衛星画像データ２の数値標高モデルデータ２ａを組み合わせて立体的にしたデータである。尚、３次元洪水マップデータ１２ｅを画像として編集し画面に描画すると、図１４に示す通りになる。

ハイエト−ハイドログラフデータ１２ｆは、気象データ４の時系列降水量データ４ｄと時系列流出量データ８ｅをハイエト−ハイドログラフ方式にしたデータである。尚、ハイエト−ハイドログラフデータ１２ｆを画像として編集し画面に描画すると、図１２に示す通りになる。

図５は、本発明である洪水予測システムのモデルの役割を示す図である。洪水予測システム１には、山岳部から河川へと水が流れ出ていく現象を計算する流域流出モデル１５と、河川内及び平野部における水の移動を計算する河道水理モデル１５ａとがある。

流域流出モデル１５は、河川の上流域についてのシミュレーションである。山岳部の降雨１５ｂが蒸発１５ｃ又は浸透１５ｄの影響を受けながら集水１５ｅし、河川に流出１５ｆする状況を計算する。

河道水理モデル１５ａは、河川の下流域についてのシミュレーションである。河川への流出１５ｆを増水量として受け取って、河川の水移動１５ｇ及び平野部への氾濫１５ｈを計算する

図６は、本発明である洪水予測システムの河川図データを画面表示した図である。地図の画像上に、河川図データ６ｇが河川１６ａとして描画される。尚、海領域図データ６ｆも海１６として描画される。

図７は、本発明である洪水予測システムの斜面方位図データを画面表示した図である。地図の画像上の各地点１６ｂそれぞれについて、斜面方位図データ６ｈからどの方向に傾斜しているか判断し、８色で色分けして表示する。尚、地点１６ｂは、コンピュータ５の画面に表示した場合は、ピクセル等に相当する。

図８は、本発明である洪水予測システムの水系図データを画面表示した図である。地図の画像上に、河川図データ６ｇの河川１６ａに加えて、水系図データ６ｉが水系１６ｃとして描画される。

図９は、本発明である洪水予測システムの流域図データを画面表示した図である。地図の画像上に、流域図データ６ｊが流域１６ｄとして描画される。

図１０は、本発明である洪水予測システムの土地被覆分布図データを画面表示した図である。地図の画像上の各地点１６ｅそれぞれについて、土地被覆分布図データ７ｅから土地の表面の状態を判断し、色分けして表示する。

図１１は、本発明である洪水予測システムのセル図データを画面表示した図である。地図の画像上に、セル図データ１０ｃのセル１６ｆを重ねて表示する。

図１２は、本発明である洪水予測システムのハイエト−ハイドログラフを画面表示した図である。ハイエト−ハイドログラフデータ１２ｆは、降水量１６ｇと流出量１６ｈのグラフを同時に示したものである。

降水量１６ｇのグラフは、左側を縦軸、下側を横軸とする折れ線グラフであり、流出量１６ｈのグラフは、右側を縦軸、上側を横軸とする逆向きの棒グラフである。降水量１６ｇのグラフも流出量１６ｈのグラフも時間変化を表しており、降水量１６ｇが多くなると、流出量１６ｈも多くなることが分かる。

図１３は、本発明である洪水予測システムの２次元洪水マップデータを画面表示した図である。地図の画像上の各セル１６ｉそれぞれについて、２次元洪水マップデータ１２ｄから水位を判断し、色分けして表示する。

図１４は、本発明である洪水予測システムの３次元洪水マップデータを画面表示した図である。立体的に表示した地図の画像上の各セル１６ｊそれぞれについて、３次元洪水マップデータ１２ｅから水位を判断し、色分けして表示する。

図１５は、本発明である洪水予測システムの勾配方向算出を示す図であり、図１６は、本発明である洪水予測システムの流域決定を示す図である。

本発明である洪水予測システム１は、粗い精度の数値標高モデルデータ２ａから、正確に流域１６ｄを決定できるという特長を持ち、斜面方位図の作成６ｃにおいて河川１６ａの位置及び形状の情報を加味し、流域図の作成６ｅを行うことにより実現する。

流域図の作成６ｅにおいて流域図データ６ｊを作成する場合、数値標高モデルデータ２ａから斜面方位図データ６ｈを作成して勾配方向１６ｋを算出し、斜面方位図データ６ｈ上で落水線追跡１７ｄを行う。

斜面方位図データ６ｈ上のある地点１６ｂにおいて、地点１６ｂに水があれば、地点１６ｂに定められた勾配方向１６ｋは、水が流れていく方向と解することができる。一の地点１６ｂから次々に勾配方向１６ｋの地点１６ｂを辿った軌跡、即ち、落水線１６ｌは、降水時の水の流路を示すこととなる。

一の地点１６ｂを指定したとき、他の全ての地点１６ｂから落水線追跡１７ｄを実行する。指定した地点１６ｂに到達する落水線１６ｌを全てまとめた領域が、指定した地点１６ｂに降水が流れ込む範囲、即ち、流域１６ｄとなる。

尚、流域決定１７ｅを行うために斜面方位図データ６ｈを作成する場合、海１６又は広大な湖沼等のような流路の終端となる場所を除き、全ての地点１６ｂの勾配方向１６ｋが定義されなければならない。

勾配方向算出１７では、対象の地点１６ｂの高度と、その周囲、即ち、東、東南、南、南西、西、西北、北、北東の８方向にある地点１６ｂの高度を比較し、最も低い地点１６ｂの方向を勾配方向として定義するが、対象の地点１６ｂ自身が最も低い場合は、勾配方向を定義することができない。

前記のように、対象の地点１６ｂが窪んでいる場合には、数値標高モデルデータ２ａに対し、窪みを埋める処理を施し、周囲の地点１６ｂのうちいずれかが最も低くなるようにした上で、勾配方向を定義する。勾配方向算出１７ａでは、中心の地点１６ｂが８であった場合、１０に補正するなどして処理する。

尚、地点１６ｂの高度は、精度の高いデータで扱われるので、周囲の地点１６ｂが同じ高度になることはほとんどない。仮に、周囲に最も低い地点１６ｂが２箇所以上存在した場合、いずれか一箇所が選択されるが、地図全体に対して一地点１６ｂは極微小であるため、誤差として大きな影響はない。

しかし、流域決定１７ｅにおいては、正確な斜面方位図データ６ｈを前提としているため、正確な数値標高モデルデータ２ａが必要となる。そこで、誤差を出来るだけ少なくするために、以下の方法を採る。

数値標高モデルデータ２ａから斜面方位図データ６ｈを作成する際、河川図データ６ｇを読み込み、高度と勾配方向１６ｋの自動書換えを行う。まず、河川図データ６ｇで河川１６ａとして定義されている地点１６ｂの高度を最低値、例えば、−９や−９９９９といった値に置き換える。

次に、河川１６ａ以外の地点１６ｂについての勾配方向算出１７ｂは、窪み対策をした上で勾配方向を算出し、河川１６ａについての勾配方向算出１７ｃは、河川図データ６ｇに定義されている流れの方向を勾配方向として算出する。

河川１６ａの位置が斜面方位図データ６ｈ上に刻み付けられるため、水の集中先が誤った場所になるのを防ぐことができ、結果として、流域決定１７ｅの正確さを飛躍的に向上させることができる。

また、水系図の作成６ｄを行い、流域算出位置指定１７ｆに反映させる。斜面方位図データ６ｈ上の全ての地点１６ｂで落水線追跡１７ｄを実行し、落水線１６ｌの通過頻度が高い地点１６ｂを水系１６ｃとして抽出する。

水系１６ｃと河川図データ６ｇに定義された河川１６ａとの交点に当たる場所が、集まった降水が河川１６ａに流れ出す出口であり、流域算出位置指定１７ｆの際に、流域算出位置１６ｍの候補となる。

流域決定１７ｅでは、可視化した流域算出位置１６ｍを指定し、流域算出位置１６ｍに至る全ての落水線１６ｌを結合することにより、流域１６ｄを抽出する。流域決定１７ｅにより、水の排出先を範囲分けすることができる。

図１７は、本発明である洪水予測システムの流域の幅の算出を示す図であり、図１８は、本発明である洪水予測システムの流域の勾配の算出を示す図である。

本発明である洪水予測システム１は、決定した流域１６ｄについて、地形的なパラメータと地質及び植生によるパラメータを算出することができ、流域特性計算８ｂにおいて、流域図データ６ｊに、数値標高モデルデータ２ａ、土地被覆分布図データ７ｅ、土壌分布図３ａ、被覆パラメータ３ｂ及び土壌パラメータ３ｃを組み合わせることで行う。

流域１６ｄについて算出されるパラメータは、面積、幅、勾配、不浸透域の占める割合、不浸透域の地表面粗度、浸透域の地表面粗度、不浸透域の初期損失雨量、浸透域の初期損失雨量、初期土湿不足、毛管吸上げ及び水理伝導度の１１個である。

面積は、一地点１６ｂ当たりの単位面積を予め定義し、流域１６ｄに含まれる地点１６ｂ数に単位面積を乗ずることにより計算する。以下、面積に関する計算は、同様にして行うものとする。

幅については、まず、流域１６ｄに含まれる水系１６ｃの中から主たるものを１つ選び、水系１６ｃの両端を延長することにより流域１６ｄを２つの領域に分断する。尚、水系１６ｃの選択と延長は、落水線１６ｌの通過頻度が最も高い場所を計算することにより自動的に行われる。

次に、水系１６ｃで分断された２つの領域のうち、大きい方の面積をＡ、小さい方の面積をＢ、分断した線の長さをＬとしたとき、（１＋２×Ｂ／（Ａ＋Ｂ））×Ｌの計算式により、流域１６ｄの幅算出１８を行う。

勾配については、まず、流域図データ６ｊと数値標高モデルデータ２ａを組み合わせて、流域１８ｂの正面１８ｃに示すように、最も高い場所である最高点１８ｄと最も低い場所である最低点１８ｅの高度を求め、高さの差を予め決めておいた自然数ｎで等分したものをΔｈとする。

最も低い場所での流域１６ｄの水平断面を取り、土に埋もれる部分の断面積１８ｇを求めて、断面積１８ｇに等しい面積を持つ円の半径Ｒ０を算出する。次に、最も低い場所からΔｈだけ上昇した場所の水平断面を取り、断面積１８ｇが等しい円の半径Ｒ１を算出し、同様にして、最も高い場所における半径Ｒｎまで、高さΔｈの間隔で、合計ｎ＋１個の半径を計算する。

また、流域１６ｄの境界線の長さ、即ち、流域１６ｄの周囲長を計算する。一地点１６ｂの縦及び横の単位長を予め定めておき、境界線上の各地点１６ｂについて、境界１８ｈとして繋がる両隣りの地点１６ｂとの位置関係に従って長さを計算する。

例えば、境界１８ｈで両隣りが西と東南という繋がり方をしている地点１６ｂでは、縦の単位長の半分の２乗と横の単位長の２乗との和の平方根を取ることで求める。各地点１６ｂにおける境界１８ｈの長さを全て合計することにより、流域１６ｄの境界線の長さを求めることができる。

流域１８ｂの平面１８ｆにおいて、流域１６ｄの境界線の長さをＬ、流域１６ｄの面積をＡとしたとき、（Δｈ×Ｌ×ΣＲｉ）／（Ａ×Ｒ０）の計算式により、流域１６ｄの勾配算出１８ａを行う。尚、ΣＲｉ＝Ｒ０＋Ｒ１＋…＋Ｒｎである。

不浸透域の占める割合は、流域図データ６ｊと数値標高モデルデータ２ａと土地被覆分布図７ｅを組み合わせて計算する。土地被覆分布図７ｅから舗装面に当たる場所を識別し、対象の流域１６ｄにおける舗装面の面積を数値標高モデルデータ２ａから計算し、当該流域１６ｄ全体の面積との比を取ることで求める。

地表面粗度及び初期損失雨量は、被覆パラメータ３ｂに定義されており、不浸透域の地表面粗度及び不浸透域の初期損失雨量は、土地被覆分布図データ７ｅで不浸透域として定められている分類について、流域１６ｄに占める面積によって重みを付け、平均を計算することにより値を求める。

同様に、浸透域の地表面粗度及び浸透域の初期損失雨量は、土地被覆分布図データ７ｅで浸透域として定められている分類について、流域１６ｄに占める面積によって重みを付け、平均を計算することにより値を求める。

初期土湿不足、毛管吸上げ及び水理伝導度は、土壌パラメータ３ｃに定義されており、土壌分布図データ３ａの各分類について、流域１６ｄに占める面積によって重みを付け、平均を計算することにより値を求める。ただし、面積のうち、土地被覆分布図データ７ｅで不浸透領域として定められている領域を除くものとする。

以上のように、本発明である洪水予測システム１は、地球観測衛星からのリモートセンシングデータを利用して洪水災害のシミュレーションを行うことにより、リアルタイムに洪水を予測できる上、コストや期間も大幅に軽減することができる。

また、上流域を対象とする流域流出モデルと、下流域を対象とする河道水理モデルに分かれており、一度に全体を動かしたり、個別に動かしたりと、利用目的に応じたシミュレーションが可能となる。

本発明である洪水予測システムの構成を示すブロック図である。 本発明である洪水予測システムの流れを示すフローチャートである。 本発明である洪水予測システムの処理の流れを示す図である。 本発明である洪水予測システムのデータの流れを示す図である。 本発明である洪水予測システムのモデルの役割を示す図である。 本発明である洪水予測システムの河川図データを画面表示した図である。 本発明である洪水予測システムの斜面方位図データを画面表示した図である。 本発明である洪水予測システムの水系図データを画面表示した図である。 本発明である洪水予測システムの流域図データを画面表示した図である。 本発明である洪水予測システムの土地被覆分布図データを画面表示した図である。 本発明である洪水予測システムのセル図データを画面表示した図である。 本発明である洪水予測システムのハイエト−ハイドログラフを画面表示した図である。 本発明である洪水予測システムの２次元洪水マップデータを画面表示した図である。 本発明である洪水予測システムの３次元洪水マップデータを画面表示した図である。 本発明である洪水予測システムの勾配方向算出を示す図である。 本発明である洪水予測システムの流域決定を示す図である。 本発明である洪水予測システムの流域の幅の算出を示す図である。 本発明である洪水予測システムの流域の勾配の算出を示す図である。

符号の説明

１ 洪水予測システム
２ 衛星画像データ
２ａ 数値標高モデルデータ
２ｂ 反射率データ
３ ＧＰＳ・ＧＩＳデータ
３ａ 土壌分布図データ
３ｂ 被覆パラメータ
３ｃ 土壌パラメータ
３ｄ 河川情報
３ｅ 氾濫原情報
４ 気象データ
４ａ 時系列海面高度データ
４ｂ 初期水位データ
４ｃ 初期流量データ
４ｄ 時系列降水量データ
４ｅ 月平均蒸発量データ
５ コンピュータ
５ａ 入力装置
５ｂ 記憶装置
５ｃ 出力装置
５ｄ 中央処理装置
６ 流域の抽出
６ａ 海領域図の作成
６ｂ 河川図の作成
６ｃ 斜面方位図の作成
６ｄ 水系図の作成
６ｅ 流域図の作成
６ｆ 海領域図データ
６ｇ 河川図データ
６ｈ 斜面方位図データ
６ｉ 水系図データ
６ｊ 流域図データ
７ 分布図の作成・編集
７ａ 土地被覆分布図の作成
７ｂ 土壌分布図の編集
７ｃ パラメータテーブルの編集
７ｄ パラメータテーブルの編集
７ｅ 土地被覆分布図データ
８ 流域流出モデルによるシミュレーション
８ａ 計算制御ファイルの編集
８ｂ 流域特性計算
８ｃ 流域流出シミュレーション
８ｄ 上流域の地理データ
８ｅ 時系列流出量データ
９ ノード・セグメントの作成
９ａ 河道分割
９ｂ 河道特性入力
９ｃ 河道区間図データ
９ｄ 河道の地理データ
１０ セルの作成
１０ａ 氾濫原分割
１０ｂ セル特性入力
１０ｃ セル図データ
１０ｄ セルの地理データ
１１ 河道水理モデルによるシミュレーション
１１ａ 計算制御ファイルの編集
１１ｂ 河川・セル結合
１１ｃ 河道水理シミュレーション
１１ｄ 下流域・河川の地理データ
１１ｅ 時系列水位データ
１２ シミュレーション結果の表示
１２ａ ２次元洪水描画
１２ｂ ３次元洪水描画
１２ｃ ハイエト−ハイドログラフ描画
１２ｄ ２次元洪水マップデータ
１２ｅ ３次元洪水マップデータ
１２ｆ ハイエト−ハイドログラフデータ
１３ データ抽出・モデル化処理
１４ 洪水予測シミュレーション
１５ 流域流出モデル
１５ａ 河道水理モデル
１５ｂ 降雨
１５ｃ 蒸発
１５ｄ 浸透
１５ｅ 集水
１５ｆ 流出
１５ｇ 水移動
１５ｈ 氾濫
１６ 海
１６ａ 河川
１６ｂ 地点
１６ｃ 水系
１６ｄ 流域
１６ｅ 地点
１６ｆ セル
１６ｇ 降水量
１６ｈ 流出量
１６ｉ セル
１６ｊ セル
１６ｋ 勾配方向
１６ｌ 落水線
１６ｍ 流域算出位置
１７ 勾配方向算出
１７ａ 勾配方向算出
１７ｂ 勾配方向算出
１７ｃ 勾配方向算出
１７ｄ 落水線追跡
１７ｅ 流域決定
１７ｆ 流域算出位置指定
１８ 幅算出
１８ａ 勾配算出
１８ｂ 流域
１８ｃ 正面
１８ｄ 最高点
１８ｅ 最低点
１８ｆ 平面
１８ｇ 断面積
１８ｈ 境界

Claims (5)

1. 衛星画像データ及びＧＰＳ・ＧＩＳデータを使用して、海領域図、河川図、斜面方位図、水系図及び流域図のデータを作成する流域の抽出と、前記衛星画像データを使用して、土地被覆分布図のデータを作成する分布図の作成・編集と、前記流域の抽出及び分布図の作成・編集で作成したデータ並びに気象データを使用して、河川の上流域についてコンピュータで流出状況を計算する流域流出モデルによるシミュレーションと、前記ＧＰＳ・ＧＩＳデータを使用して、河川の分断地点であるノード、及び分断された河川であるセグメントごとにデータを作成するノード・セグメントの作成と、前記ＧＰＳ・ＧＩＳデータを使用して、単位領域であるセルごとにデータを作成するセルの作成と、前記流域流出モデルによるシミュレーション、ノード・セグメントの作成及びセルの作成で作成したデータを使用して、河川の下流域についてコンピュータで氾濫状況を計算する河道水理モデルによるシミュレーションと、前記河道水理モデルによるシミュレーションで作成したデータを編集して、コンピュータの画面に表示するシミュレーション結果の表示とからなることを特徴とする洪水予測システム。
2. 斜面方位図データを作成する際、衛星画像データに対し、河川図データにより河川を定義した上で勾配方向を算出することを特徴とする請求項１に記載の洪水予測システム。
3. 水系図データを作成する際、斜面方位図データが示す勾配方向に従って落水線追跡を行い、前記落水線の通過頻度の高い地点を抽出することを特徴とする請求項１又は請求項２に記載の洪水予測システム。
4. 流域図データを作成する際、河川図データの河川と水系図データの水系との交点を流域算出位置として指定し、前記流域算出位置に至る全ての落水線を結合することにより、流域を抽出することを特徴とする請求項１、請求項２又は請求項３に記載の洪水予測システム。
5. 流域流出モデルによるシミュレーションにおいて、流域図データに、衛星画像データ、土地被覆分布図データ及びＧＰＳ・ＧＩＳデータを組み合わせて、地形的なパラメータと地質及び植生によるパラメータを算出することができることを特徴とする請求項１、請求項２、請求項３又は請求項４に記載の洪水予測システム。
   

Images