(57)【特許請求の範囲】
【請求項１】
計算機システムを制御するプログラムであって、
前記計算機システムは、前記プログラムが格納される記憶装置と、前記記憶装置に格納されたプログラムを実行するプロセッサと、を備え、
前記プログラムは、
浸水域の輪郭線上の点の座標を入力として受け付ける第１手順と、
前記浸水域の輪郭線を含む複数の区画の浸水位を、前記各区間の地盤の標高に所定の浸水深を加算することによって算出する第２手順と、
前記浸水域内の一つ以上の第１地点における浸水深を入力として受け付ける第３手順と、
前記第１地点を含む前記区画の浸水位を、前記第１地点を含む前記区画の地盤の標高に前記第１地点の浸水深を加算することによって算出する第４手順と、
前記第２手順及び前記第４手順によって算出された浸水位を空間補間することによって、前記浸水域に含まれる前記区画のうち、前記第１地点を含まない区画の浸水位を算出する第５手順と、
前記第５手順によって算出された前記区画の浸水位から、前記区画の地盤の標高を減算することによって、前記区画の浸水深を算出する第６手順と、
前記算出された前記区画の浸水深を出力する第７手順と、を前記プロセッサに実行させることを特徴とするプログラム。
【請求項２】
前記第５手順において、前記第１地点を含まない区画の浸水位は、前記第２手順及び前記第４手順によって算出された浸水位、及び、前記第１地点を含まない区画の浸水位が、ラプラス方程式を満たすように算出されることを特徴とする請求項１に記載のプログラム。
【請求項３】
前記第５手順において、前記第１地点を含まない区画である第１区画の浸水位は、前記第１区画に隣接する複数の区画の浸水位に所定の計数を乗じ、前記所定の計数を乗じた複数の区画の浸水位を加算することによって算出され、
前記第５手順は、
前記第１区画の浸水位を算出する手順と、
前記算出された値と、前回算出された前記第１区画の浸水位との差を算出する手順と、
前記算出された差の絶対値が所定の閾値より大きい場合、前記第１区画及び前記第１区画に隣接する前記複数の区画の浸水位を算出する手順と、を含むことを特徴とする請求項１に記載のプログラム。
【請求項４】
前記区画は、前記浸水域を含む地域を所定の間隔で分割することによって生成されることを特徴とする請求項１に記載のプログラム。
【請求項５】
前記計算機システムは、さらに表示装置を備え、
前記第７手順は、前記入力された第１地点の浸水深と、前記算出された各区画における浸水深と、の少なくとも一方を前記表示装置に表示させる第８手順を含むことを特徴とする請求項１に記載のプログラム。
【請求項６】
前記計算機システムは、前記浸水域を含む地域の地図情報を保持し、
前記第７手順は、前記地図情報に基づいて、前記浸水域を含む地域の地図を前記表示装置に表示させる第９手順をさらに含み、
前記第９手順によって表示される地図には、前記浸水域の輪郭線が表示され、
前記入力された第１地点の浸水深と、前記算出された各区画における浸水深と、の少なくとも一方は、前記第９手順によって表示された地図の上に表示されることを特徴とする請求項５に記載のプログラム。
【請求項７】
前記第７手順において、前記算出された各区画における浸水深の表示は、前記地図上の前記各区画に対応する領域のうち、前記浸水域内の領域に、前記浸水深と対応付けられた色彩を表示することによって実行されることを特徴とする請求項６に記載のプログラム。
【請求項８】
前記計算機システムは、
さらに表示装置及び入力装置を備え、
一つ以上の災害の識別情報と、前記各災害に対応して既に入力されている前記第１地点の浸水深と、前記各災害が発生した地域の地図情報と、を保持し、
前記プログラムは、さらに、
前記災害の識別情報を前記表示装置に表示させる第１０手順と、
前記表示された識別情報のうち一つの選択を前記入力装置によって受け付ける第１１手順と、
前記選択された識別情報によって識別される災害における浸水域を含む地域の地図、及び、その災害に対応して既に入力されている前記第１地点の浸水深を、前記保持された地図情報に基づいて前記表示装置に表示させる第１２手順と、を前記プロセッサに実行させ、
前記第１手順は、前記第１２手順によって表示された地図の上の点の指定を受け付けることによって実行され、
前記第３手順は、前記第１２手順によって表示された地図の上の点の指定を受け付け、さらに、前記指定された点に対応する地点の浸水深の入力を受け付けることによって実行されることを特徴とする請求項１に記載のプログラム。
【請求項９】
サーバ計算機と、前記サーバ計算機にネットワークを介して接続されるクライアント計算機と、を備える計算機システムであって、
前記サーバ計算機は、第１プロセッサと、前記第１プロセッサに接続される第１記憶装置と、を備え、
前記クライアント計算機は、第２プロセッサと、前記第２プロセッサに接続される第２記憶装置と、前記第２プロセッサに接続される表示装置と、前記第２プロセッサに接続される入力装置と、を備え、
前記第１プロセッサは、
浸水域の輪郭線上の点の座標を、前記ネットワークを介して、入力として受け付ける第１手順と、
前記浸水域の輪郭線を含む複数の区画の浸水位を、前記各区間の地盤の標高に所定の浸水深を加算することによって算出する第２手順と、
前記浸水域内の一つ以上の第１地点における浸水深を、前記ネットワークを介して、入力として受け付ける第３手順と、
前記第１地点を含む前記区画の浸水位を、前記第１地点を含む前記区画の地盤の標高に前記第１地点の浸水深を加算することによって算出する第４手順と、
前記第２手順及び前記第４手順によって算出された浸水位を空間補間することによって、前記浸水域に含まれる前記区画のうち、前記第１地点を含まない区画の浸水位を算出する第５手順と、
前記第５手順によって算出された前記区画の浸水位から、前記区画の地盤の標高を減算することによって、前記区画の浸水深を算出する第６手順と、を実行することを特徴とする計算機システム。
【請求項１０】
前記第５手順において、前記第１地点を含まない区画の浸水位は、前記第２手順及び前記第４手順によって算出された浸水位、及び、前記第１地点を含まない区画の浸水位が、ラプラス方程式を満たすように算出されることを特徴とする請求項９に記載の計算機システム。
【請求項１１】
前記第１プロセッサは、前記入力された第１地点の浸水深と、前記算出された各区画における浸水深と、の少なくとも一方を、前記ネットワークを介して前記クライアント計算機に送信することを指示し、
前記第２プロセッサは、前記入力された第１地点の浸水深と、前記算出された各区画における浸水深と、の少なくとも一方を前記表示装置に表示させる第７手順を実行することを特徴とする請求項９に記載の計算機システム。
【請求項１２】
前記計算機システムは、一つ以上の災害の識別情報と、前記各災害に対応して既に入力されている前記第１地点の浸水深と、前記各災害が発生した地域の地図情報と、を保持し、
前記第２プロセッサは、さらに、
前記災害の識別情報を前記表示装置に表示させる第９手順と、
前記表示された識別情報のうち一つの選択を前記入力装置によって受け付ける第１０手順と、
前記選択された識別情報によって識別される災害における浸水域を含む地域の地図、及び、その災害に対応して既に入力されている前記第１地点の浸水深を、前記保持された地図情報に基づいて前記表示装置に表示させる第１１手順と、
前記第１１手順によって表示された地図の上の複数の点の指定、及び、前記指定された複数の点のうち少なくとも一つに対応する地点の浸水深の入力を受け付ける第１２手順と、
前記第１２手順において指定された複数の点の少なくとも一つの座標、及び、前記第１２手順において入力された浸水深を、前記サーバ計算機に前記ネットワークを介して送信する第１３手順と、を実行し、
前記第１手順は、前記第１３手順によって送信された座標を受け付けることによって実行され、
前記第３手順は、前記第１３手順によって送信された浸水深を受け付けることによって実行されることを特徴とする請求項９に記載の計算機システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本願明細書で開示される技術は、水害の現地調査を支援するプログラムに関する。本願明細書で開示される技術は特に、水害に遭った地域の現地調査によって得られた調査地点ごとの浸水深を空間補間するプログラムに関する。
【背景技術】
【０００２】
近年、集中豪雨や台風などによる水害が頻発し、効果的な被害軽減策が求められている。被害軽減のためには被害の状況を調べ、原因を推定することが有効である。この様な背景から、実際に起きた水害の調査が行われている。調査においては、主に浸水深と浸水域が計測される。ここで浸水深とは、ある水害においてある地点で発生した最大の水深である。また、浸水域とは、ある水害において浸水した地域である。浸水深が０となる等値線が浸水域の境界であると定義される。
【０００３】
浸水深の調査方法を記載した文献として、非特許文献１が挙げられる。調査の方法は以下の通りである；非特許文献１の筆者らは、ある水害で浸水した地域内の１４６地点において浸水深を計測した。計測には、レーザ距離計（ｌａｓｅｒ  ｒａｎｇｅｆｉｎｄｅｒ）と呼ばれる測量用の装置が利用された。この装置によって、ある点と、その点から見える別の点との標高差を計測することができる。彼らは、標高点と、その点から見える水面の痕跡との標高差を計測した。ここで、標高点とは、行政機関によって緯度、経度及び標高が公表されている地点であり、道路など行政機関が管理している領域に存在する。日本の都市では、おおむね１００ｍから１０００ｍ程度の間隔で標高点が設置されている。標高点の標高の精度は０．１ｍ程度と考えられている。また、水面の痕跡とは、建物、塀、街路樹、電柱などの都市構造物に残された泥などであり、この水面の痕跡は、その水害においてその地点で発生した最も高い水面に一致すると考えられる。非特許文献１の筆者らは、このように計測した標高差を、浸水深とした。さらに、彼らは、浸水深に標高点の標高を加えた値を浸水位（すなわち、浸水した時の水面の標高）とし、浸水位の分布図を作成した。
【０００４】
また、浸水域の調査結果が非特許文献２に記載されている。非特許文献２には、浸水域は航空写真から推定したと記載されている。非特許文献２によれば、浸水した領域を線分で囲むことによって、浸水域の形状が地図上で表現されている。このように、従来、浸水深の現地調査と浸水域の現地調査とは、それぞれ独立な方法によってなされていた。
【０００５】
ところで近年、航空レーザ計測と呼ばれる技術が開発され、地表面の標高を０．１ｍ程度の精度、５ｍ程度の空間分解能で計測できるようになった。この航空レーザ計測技術を用いれば、５ｍ間隔で、標高点と同程度の精度の標高を得ることができるようになると考えられる。
【０００６】
この高精度の地形データを用いた水害に関連する技術が開発されている。その例として、特許文献１に記載された「標高別地物占有率計測方法及びこれを用いた浸水深補正方法」が挙げられる。特許文献１は、都市部のような密集市街地において正確に氾濫流の挙動を把握するため、高密度標高データから地物占有率を算出する方法、及び、その地物占有率を用いてシミュレーション精度を向上させる方法を開示している。この方法ではまず、格子内の標高データを「地盤」と「地物」のどちらかに分離する。次に、ある標高ｈにおける地物占有率ｒをｒ（ｈ）＝ｎ（ｈ）／Ｎとして算出する。ここで、Ｎ：格子内の標高データ点総数、ｎ（ｈ）：ある標高ｈにおける「地物」の標高データ点数である。さらに、このｒ（ｈ）を用いて、氾濫シミュレーションで求めた水位ｈを補正する。この方法によって、地物が氾濫流に大きな影響を与える都市部において、シミュレーション結果の精度を向上できると記載されている。
【特許文献１】特開２００５−１７２６３４号公報
【非特許文献１】山本博文，福井市街足羽川左岸および鯖江市河和田地区における浸水被害について，平成１６年７月新潟・福島、福井豪雨災害に関する調査研究，平成１６年度科学研究費補助金（特別研究促進費（１））研究報告書，課題番号１６８００００１，ｐ１２１，２００７．
【非特許文献２】Ｒｏｂｅｒｔ  Ｊ．Ｃｏｎｎｅｌｌ，  Ｄａｖｉｄ  Ｊ．Ｐａｉｎｔｅｒ，  ａｎｄ  Ｃｏｒｎｅｌ  Ｂｅｆｆａ：Ｔｗｏ−ｄｉｍｅｎｓｉｏｎａｌ  ｆｌｏｏｄ  ｐｌａｎｅ  ｆｌｏｗ．ＩＩ：ｍｏｄｅｌ  ｖａｌｉｄａｔｉｏｎ，Ｊｏｕｒｎａｌ  ｏｆ  ｈｙｄｒｏｌｏｇｉｃ  ｅｎｇｉｎｅｅｒｉｎｇ，Ａｍｅｒｉｃａｎ  ｓｏｃｉｅｔｙ  ｏｆ  ｃｉｖｉｌ  ｅｎｇｉｎｅｅｒｓ，Ｖｏｌｕｍｅ  ６，Ｉｓｓｕｅ  ５，ｐｐ．４０６−４１５，２００１．
【発明の開示】
【発明が解決しようとする課題】
【０００７】
浸水被害は、浸水位ではなく、浸水深に強く依存する。この浸水深は小さい空間スケールで変化するため、浸水被害の面的な分布を知るためには、調査した地点における浸水深だけでなく浸水深の面的な分布が必要である。ところが従来の技術では、浸水深の面的な分布を求めることができなかった。例えば、非特許文献１では、複数地点で計測した浸水位から浸水位の面的な分布を求める一方、浸水深については調査した地点以外の値を求めていない。
【０００８】
一般に、複数地点における情報からその情報の分布を推定するためには空間補間が用いられる。そこで発明者らは、まず浸水深を空間補間することを試みた。ところがこの方法には、以下２つの問題があることがわかった。
【０００９】
第１の問題は、得られる浸水深分布が不自然なことである。これは、浸水深の空間スケールより粗い空間スケールで計測した値を空間補間したためであると考えられる。
【００１０】
第２の問題は、得られる浸水深分布が、現地調査で得られる浸水域と矛盾するおそれがあることである。上述したように、浸水深と浸水域とはそれぞれ独立な現地調査によって得られる。ところが、これらの情報は独立ではない。浸水域の定義から、浸水域の境界の外側にある地点では浸水深が０、浸水域の境界の内側にある地点では浸水深が０より大とならなければならない。浸水域が浸水深分布に与えるこの制約を、浸水深のみを用いた補間によって得られた浸水深分布が満たすとは限らない。
【００１１】
第１の問題を解決する方法として、非特許文献１にあるように浸水位を計測し、浸水深分布を求める方法が考えられる。浸水位から標高を減算すると浸水深が得られるのは定義から明らかであるため、浸水位分布から標高の分布を減算すると浸水深分布が得られることも明らかであるといえる。ところが、上述したように浸水位の計測には測量用の装置が必要であるため、通常の現地調査では浸水位が計測されない。通常の現地調査では、測量用の装置が不要な浸水深のみが計測される。したがって、この方法は通常の現地調査には適用できないといえる。
【００１２】
第１の問題を解決する別の方法として、浸水深の計測間隔を空間的に細かくする方法が考えられる。さらに、浸水深が０となる近傍で浸水深の計測間隔を空間的にほぼ無限小にすれば、浸水域を浸水深分布によって表現できるため、第１の問題だけでなく、第２の問題も解決できる。ところが、この方法は莫大な手間がかかる。したがって、この方法は現実的でないといえる。
【００１３】
以上のように、産業的な観点からは、従来の現地調査によって得られる浸水深及び浸水域を入力として、浸水域内の浸水深分布を適切に推定する方法が望ましいと言える。
【００１４】
本発明はこれらの問題に鑑みてなされたものであり、複数地点の浸水深から、浸水深の面的な分布を推定することを目的とする。
【課題を解決するための手段】
【００１５】
本願で開示する代表的な発明は、計算機システムを制御するプログラムであって、前記計算機システムは、前記プログラムが格納される記憶装置と、前記記憶装置に格納されたプログラムを実行するプロセッサと、を備え、前記プログラムは、浸水域の輪郭線上の点の座標を入力として受け付ける第１手順と、前記浸水域の輪郭線を含む複数の区画の浸水位を、前記各区間の地盤の標高に所定の浸水深を加算することによって算出する第２手順と、前記浸水域内の一つ以上の第１地点における浸水深を入力として受け付ける第３手順と、前記第１地点を含む前記区画の浸水位を、前記第１地点を含む前記区画の地盤の標高に前記第１地点の浸水深を加算することによって算出する第４手順と、前記第２手順及び前記第４手順によって算出された浸水位を空間補間することによって、前記浸水域に含まれる前記区画のうち、前記第１地点を含まない区画の浸水位を算出する第５手順と、前記第５手順によって算出された前記区画の浸水位から、前記区画の地盤の標高を減算することによって、前記区画の浸水深を算出する第６手順と、前記算出された前記区画の浸水深を出力する第７手順と、を前記プロセッサに実行させることを特徴とする。
【発明の効果】
【００１６】
本発明の一実施形態によれば、調査地点以外の地点における浸水深を算出できる。したがって、いくつかの地点において浸水深を調査した結果に基づいて、ある地域全体の浸水深分布を算出できる。これによって、水害の全体像を把握しやすくなる。
【００１７】
また、入力データとして、特殊な装置で計測する必要がある浸水位ではなく、一般的な装置で計測できる浸水深を利用できる。これによって現地調査の手間が軽減される。
【発明を実施するための最良の形態】
【００１８】
最初に、本発明の概要を説明する。
【００１９】
本発明では、浸水深をより空間スケールの大きい量である浸水位に変換し、さらに浸水位が従う自然法則を用いて浸水位を空間補間する。本発明は発明者らの３つの発見の上に成り立っている。本発明は、これらの発見に基づく仮定が妥当と見なせる精度で有効であるといえる。
【００２０】
第１の発見は、浸水深の空間スケールに比べて浸水位の空間スケールが大きいと仮定できることである。この発見は、水理学における比水頭の理論（ｓｐｅｃｉｆｉｃ−ｈｅａｄ）に基づく。この理論では、水害のように流れが遅い場合、水深に比べ水位は変化しにくいことがいえる。この仮定が成り立つとき、ある間隔で計測した浸水深を補間した結果が不自然であっても、同じ間隔で計測した浸水位を補間した結果が自然である可能性がある。
【００２１】
第２の発見は、近年得られるようになった標高データを用いれば、任意の場所の標高を推定できると仮定できることである。例えば航空レーザ計測によって得られるデータの計測間隔は、標高の空間スケールよりも小さいと考えられる。したがって、このデータを利用すれば、任意の場所の標高を精度よく推定できると考えられる。任意の場所の標高が得られれば、任意の場所で計測した浸水深を、その値に標高を加えた値である浸水位に変換できる。同様に、浸水位分布を、その分布から標高分布を減じた分布である浸水深分布に変換できる。
【００２２】
第３の発見は、浸水位ｗｌｅｖがラプラス方程式（Ｌａｐｌａｃｅ  ｅｑｕａｔｉｏｎ）
ｄｉｖ（ｇｒａｄ（ｗｌｅｖ））＝０  ・・・式（１）
を満たすと仮定できることである。ただし、ｄｉｖ（）は、引数の発散（ｄｉｖｅｒｇｅｎｃｅ）を返す関数、ｇｒａｄ（）は、引数の勾配（ｇｒａｄｉｅｎｔ）を返す関数である。式（１）が成り立つとき、浸水域内の浸水位分布を求める問題は、ディリクレ問題（Ｄｉｒｉｃｈｌｅｔ  ｐｒｏｂｌｅｍ）として知られる数学的問題に帰着する。ディリクレ問題とは、ラプラス方程式を満たす関数に境界条件として領域の境界における値を与え、領域の内部の分布を求める問題である。
【００２３】
式（１）の仮定は、水害において静水圧平衡（ｈｙｄｒｏｓｔａｔｉｃ  ｂａｌａｎｃｅ）が成り立つと見なせるため、妥当であると考えられる。静水圧平衡が成り立つとき、圧力ｐの空間変化ｇｒａｄ（ｐ）は水位ｗｌｅｖの空間変化ｇｒａｄ（ｗｌｅｖ）と比例関係にある。また、圧力ｐはラプラス方程式ｄｉｖ（ｇｒａｄ（ｐ））＝０を満たすことが知られている。したがって、静水圧平衡が成り立つとき、水位ｗｌｅｖはラプラス方程式を満たすと見なせる。
【００２４】
式（１）に与える境界条件は、浸水域の境界における浸水位である。浸水域の定義から、浸水域の境界では浸水深が０である。したがって浸水域の境界における浸水位は標高と同じ値になる。
【００２５】
これらの発見に基づき発明者らが考案した、浸水深の補間方法は以下のステップからなる；  まず、空間を任意の形状の区画に分割し、区画ごとに標高を与える。この際、精度及び空間分解能が十分な標高データを利用する。次に、浸水深を計測した地点を含む区画における浸水位を算出する。同様に、浸水域の境界線を含む区画における浸水位を算出する。さらに、浸水域に含まれる区画の浸水位を、浸水位がラプラス方程式を満たすことを仮定して補間する。区画ごとの浸水位から、標高を減算して浸水深を算出する。以上のステップによって、区画ごとの浸水深が得られる。
【００２６】
以下、本発明の実施形態を添付図面に基づいて説明する。想定されるユーザが特定の一人であるシステムを第１の実施の形態として、想定されるユーザが不特定多数であるシステムを第２の実施の形態として示す。
【００２７】
最初に、本発明の第１の実施の形態について説明する。
【００２８】
図１は、本発明の第１の実施の形態のコンピュータシステムの構成を示す概略ブロック図である。
【００２９】
本実施の形態のコンピュータシステムは、ＰＣ（Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ）１００によって構成される。ＰＣ１００は、入出力部１１０を備える。ユーザは、入出力部１１０からシステムを操作する。入出力部１１０は、例えばキーボード１１１及びマウス１１２等からなる入力部と、例えばディスプレイ１１３等からなる出力部とを備える。ユーザは、入力部から命令を入力し、出力部から結果を確認することができる。
【００３０】
ＰＣ１００は、さらに、メモリ１２０及びＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）１３０を備える。メモリ１２０には入力支援プログラム１４１及び浸水深補間プログラム１４２が展開される。これらの入力支援プログラム１４１と浸水深補間プログラム１４２とは、互いにデータ及び命令を送受信することができる。これらの入力支援プログラム１４１及び浸水深補間プログラム１４２に基づいて、ＣＰＵ１３０が命令を実行する。したがって、以下の説明において入力支援プログラム１４１及び浸水深補間プログラム１４２が実行する処理は、実際には、入力支援プログラム１４１及び浸水深補間プログラム１４２に基づいて、ＣＰＵ１３０によって実行される。
【００３１】
なお、図１には、入力支援プログラム１４１及び浸水深補間プログラム１４２を独立したプログラムとして記載しているが、これらのプログラムは、一つのプログラムを構成するプログラムモジュールであってもよい。
【００３２】
この入力支援プログラム１４１及び浸水深補間プログラム１４２は、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）１９０に格納された地図データ１９１、浸水域データ１９２、浸水地点データ１９３、浸水深分布データ１９４、標高データ１９５、及び災害データ１９６を必要に応じて読み出す。ＨＤＤ１９０は、ＰＣ１００に内蔵されてもよいし、ＰＣ１００の外部に接続されてもよい。
【００３３】
浸水域データ１９２、浸水地点データ１９３、浸水深分布データ１９４及び災害データ１９６については、後で詳細に説明する（それぞれ、図７、図９、図１２及び図３参照）。地図データ１９１は、少なくとも、災害データ１９６に登録されている災害における浸水域を含む地域の地図情報を含む。標高データ１９５は、少なくとも、浸水域を含む地域の地盤の標高値を含む。本実施の形態は、標高データ１９５が、０．１ｍ程度の精度、５ｍ程度の空間分解能を有することを前提としている。このような標高データ１９５は、例えば航空レーザ計測のような技術によって得られる。
【００３４】
図１では、入力支援プログラム１４１と浸水深補間プログラム１４２とが同一のメモリ１２０に展開される構成を示した。ただし、入力支援プログラム１４１と浸水深補間プログラム１４２とが互いにデータ及び命令を送受信できる限り、それぞれのプログラムが異なるメモリに展開されてもよい。いずれの場合であっても、入力支援プログラム１４１と浸水深補間プログラム１４２とはプロセス間通信によってデータ及び命令を送受信すればよい。
【００３５】
図２は、本発明の第１の実施の形態のＰＣ１００によって実行される第１の処理を示すシーケンス図である。
【００３６】
以下、図２にしたがってシステムの処理シーケンスを説明する。
【００３７】
ＰＣ１００で入力支援プログラム１４１が起動すると、入力支援プログラム１４１は、ＨＤＤ１９０から災害データ１９６を読み出す（ステップ２０１）。データを読み出した入力支援プログラム１４１は、ウィンドウ３００からなる第１のＧＵＩ（Ｇｒａｐｈｉｃａｌ  Ｕｓｅｒ  Ｉｎｔｅｒｆａｃｅ）をディスプレイ１１３に出力する（ステップ２０２）。
【００３８】
図３は、本発明の第１の実施の形態の入力支援プログラム１４１が提供する第１のＧＵＩを示す説明図である。
【００３９】
図３のＧＵＩは、ディスプレイ１１３の所定の領域に表示されるウィンドウ３００からなり、キーボード１１１及びマウス１１２による入力を受け付ける。このＧＵＩを用いてユーザは、興味のある水害のレコードを選択する（ステップ２０３）。本実施の形態では、ユーザが選択ボタン３０１を選択状態にし、さらに「次へ」ボタン３０２を押下することによって水害のレコードが選択される。
【００４０】
ウィンドウ３００には、１以上のレコードを含むテーブル３１０が表示される。テーブル３１０は、ＨＤＤ１９０に格納された災害データ１９６に基づいて生成される。このテーブル３１０の各レコードには、ユーザが災害を一意に特定するために用いられる情報が記載される。図３の例では、災害の名称３１２及び災害にあった代表的な地点（すなわち代表点）の住所３１３が記載される。さらに、災害データ１９６の各レコードに、所定の方法によって生成した一意な番号（すなわち災害ＩＤ）３１１が与えされる。さらに、図３の例では、空白レコード３０３が用意される。災害データ１９６に存在しない災害についてユーザが興味を持った場合、この空白レコードに、ユーザ自身がキーボード１１１、マウス１１２等を用いてレコードを入力することができる。
【００４１】
入力支援プログラム１４１は、ステップ２０３で選択されたレコードの代表点３１３周辺の地図データ１９１を読み出す。さらに入力支援プログラム１４１は、ステップ２０３で選択されたレコードの災害ＩＤ３１１と一致する災害ＩＤを有するレコード（以下、該当レコードと記載する）を、浸水域データ１９２、浸水地点データ１９３、及び浸水深分布データ１９４のそれぞれから読み出す（ステップ２０４）。なおステップ２０４の結果は、
Ｃａｓｅ  Ａ：  該当するレコードが浸水域データに存在しない場合
Ｃａｓｅ  Ｂ：  該当するレコードが浸水域データに存在する場合
の２通りに分類できる。Ｃａｓｅ  Ａの場合、浸水域データだけでなく、浸水地点データ、浸水深分布データのどちらにも該当するレコードが存在しない。Ｃａｓｅ  Ｂの場合、浸水地点データ、浸水深分布データにそれぞれ該当するレコードが存在する可能性がある。
【００４２】
入力支援プログラム１４１は、ステップ２０４で読み出した地図データ１９１に基づいて、ウィンドウ４００からなる第２のＧＵＩをディスプレイ１１３に出力する（ステップ２０５）。
【００４３】
図４は、本発明の第１の実施の形態の入力支援プログラム１４１が提供する第２のＧＵＩの初期状態を示す説明図である。
【００４４】
図４のＧＵＩは、ディスプレイ１１３の所定の領域に表示されるウィンドウ４００からなり、キーボード１１１及びマウス１１２による入力を受け付ける。ウィンドウ４００には、少なくとも地図４１０が表示される。さらに、ウィンドウ４００は少なくとも、ユーザに多角形をした浸水域を指定すること、及び、１以上の地点の浸水深を指定することを可能にする機能を有する。後述するように、浸水域及び浸水深を指定するために、マウスポインタ４０１が使用される。
【００４５】
地図４１０は、選択された災害における浸水域を含む地域の地図である。地図４１０の縮尺は、浸水分布の空間スケールと同程度であることが望ましい。典型的な都市であれば、１／２５，０００から１／５００が適当と考えられる。地図４１０には水路４１１及び道路４１２等を表示することができる。
【００４６】
図５は、本発明の第１の実施の形態のＰＣ１００によって実行される第２の処理を示すシーケンス図である。
【００４７】
まずユーザは、ウィンドウ４００で浸水域の形状を入力する（ステップ５０１）。このステップは、前述したステップ２０４でＣａｓｅ  Ａの場合に必要である。Ｃａｓｅ  Ｂの場合、このステップで作成対象となるデータがすでにステップ２０４で取得されているため、ステップ５０１は省略できる。
【００４８】
図６は、本発明の第１の実施の形態において、浸水域を入力中の状態にある第２のＧＵＩを示す説明図である。
【００４９】
このＧＵＩを用いてユーザは、ステップ２３０で選択した水害の浸水域の形状を入力する。本実施の形態では、ユーザは、浸水域の形状をポリゴン（３以上の頂点からなる２次元閉多角形）として入力する。ユーザは、マウスポインタ４０１で地図４１０の１点を指示し、マウス１１２をクリックすることによって、ポリゴンの頂点を指定する。このとき指示される地図４１０の１点は、浸水域の輪郭線上の任意の点である。浸水域の輪郭線とは、浸水深が０である領域（すなわち浸水していない領域）と、浸水深が０より大きい領域（すなわち浸水している領域）との境界線である。
【００５０】
このように、浸水域の輪郭線上の点を入力する操作を、浸水域を周回する方向に３回以上繰り返すことによって、３以上の頂点が指定される。ユーザが確定ボタン４０２を押下すると、入力支援プログラム１４１は、入力された順に頂点をたどることによって、ポリゴン６２５を作成する。そして、入力支援プログラム１４１は、このポリゴン６２５に相当する浸水域データ１９２を作成する。
【００５１】
図７は、本発明の第１の実施の形態の浸水域データ１９２のデータ構造を示す説明図である。
【００５２】
浸水域データ１９２は、少なくとも、災害を一意に識別するＩＤ（識別子）、浸水域を一意に識別するＩＤ、浸水域を示すポリゴンの頂点の位置及び頂点の順番を含む。本実施の形態では、浸水域データ１９２は、災害ＩＤ７０１、浸水域ＩＤ７０２、頂点ＩＤ７０３、次の頂点ＩＤ７０４、頂点の経度７０５及び緯度７０６からなる。
【００５３】
災害ＩＤ７０１は、図３の災害ＩＤ３１１と同様、災害を一意に識別する識別子である。
【００５４】
浸水域ＩＤ７０２は、浸水域の識別子である。浸水域、すなわちポリゴンごとに一意の番号が浸水域ＩＤ４０２として割り当てられる。
【００５５】
頂点ＩＤ７０３は、入力されたポリゴンの頂点の識別子である。
【００５６】
次の頂点ＩＤ７０４は、頂点ＩＤ７０３によって識別される頂点の次に入力された頂点の識別子である。
【００５７】
経度７０５及び緯度７０６は、頂点ＩＤ７０３によって識別される頂点の経度及び緯度である。
【００５８】
このようにして浸水域データが入力された後、浸水深補間プログラム１４１は、ユーザに浸水地点データの入力を促す（ステップ５０２）。このステップ５０２は省略できる。前述したステップ２０４でＣａｓｅ  Ａの場合、このステップ５０２を省略すると浸水地点レコードが存在しないことになる。前述したステップ２０４でＣａｓｅ  Ｂの場合、このステップ５０２を省略すると浸水地点レコードが存在しない可能性があることになる。浸水地点レコードが存在しない場合、後述するステップ５０５において、浸水域６２５の境界線における浸水深のみに基づく空間補間を実行する。
【００５９】
図８は、本発明の第１の実施の形態において、浸水地点を入力中の状態にある第２のＧＵＩを示す説明図である。
【００６０】
ユーザは、浸水域６２５に含まれる１以上の地点の浸水深を入力することができる。本実施の形態では、入力支援プログラム１４１は、マウスポインタ４０１で指定された地点にアイコン８２４を表示させる。ユーザは、キーボード１１１によってその地点の浸水深をｃｍ単位で入力することができる。
【００６１】
なお、ユーザが入力位置を把握しやすいよう、浸水深の入力待ち状態にあるアイコン８２４は、浸水深の入力が済んだアイコン８２１、８２２、８２３とは異なる態様で表示されることが望ましい。図８の例では、アイコン８２４は破線の円形として、アイコン８２１、８２２及び８２３は実線の円形として表示される。
【００６２】
最後にユーザが確定ボタン４０２を押下することで、アイコン８２１、８２２、８２３及び８２４に相当する浸水地点データ１９３が生成される。なお、ユーザには浸水域の内側の地点を浸水地点として指定させる。これを実現するために、入力支援プログラム１４１は、ユーザが浸水域の内側の地点をクリックした場合のみに、浸水深の入力待ち状態にあるアイコン８２４を出現させればよい。
【００６３】
図９は、本発明の第１の実施の形態の浸水地点データ１９３のデータ構造を示す説明図である。
【００６４】
浸水地点データ１９３は、少なくとも、災害を一意に識別するＩＤ、浸水域を一意に識別するＩＤ、調査地点の位置、及び、その地点の浸水深の情報を含む。例えば本実施の形態では、浸水地点データ１９３は１以上のレコードから構成される。ただし、浸水域データに存在しない浸水域ＩＤが、浸水地点データの浸水域ＩＤに存在しないようにする。つまり、ユーザが浸水深をまだ入力していない場合、浸水地点データ１９３はレコードを含まない。１レコードは、災害ＩＤ９０１、浸水域ＩＤ９０２、調査地点の経度９０３、調査地点の緯度９０４、及び調査地点の浸水深９０５から構成される。
【００６５】
災害ＩＤ９０１は、図３の災害ＩＤ３１１と同様、災害を一意に識別する識別子である。
【００６６】
浸水域ＩＤ９０２は、図７の浸水域ＩＤ７０２と同様、浸水域を識別する識別子である。
【００６７】
経度９０３、緯度９０４及び浸水深９０５は、ユーザによって指定された調査地点の座標及びその地点の浸水深として入力された値である。
【００６８】
このようにして浸水地点データが入力された後、入力支援プログラム１４１はユーザに確定ボタン４０２の押下を促す。ユーザが確定ボタン４０２を押下すると、ステップ５０１で作成された浸水域データ１９２及びステップ５０２で作成された浸水地点データ１９３が確定する（ステップ５０３）。
【００６９】
次に入力支援プログラム１４１は、浸水深補間プログラム１４２に浸水深補間の要求を送信する（ステップ５０４）。この際、入力支援プログラム１４１は、ステップ５０３で確定した浸水域データ１９２及び浸水地点データ１９３を浸水深補間プログラム１４２に送信する。
【００７０】
入力支援プログラム１４１から浸水深補間補間の要求を受信した浸水深補間プログラム１４２は浸水深補間処理を実行する（ステップ５０５）。
【００７１】
図１０は、本発明の第１の実施の形態の浸水深補間処理（ステップ５０５）の詳細を示すＰＡＤ（Ｐｒｏｂｌｅｍ  Ａｎａｌｙｓｉｓ  Ｄｉａｇｒａｍ）図である。
【００７２】
この処理は、ステップ１００１から１０１３よりなる。この処理に先立ち、あらかじめ所定の方法により空間に隙間なく格子を作成しておく。ここで、格子とは空間を複数の小領域に分割した場合の１つの小領域（すなわち区画）である。
【００７３】
図１１は、本発明の第１の実施の形態における浸水地点、浸水域及び格子のそれぞれの位置関係を示す説明図である。
【００７４】
本発明は格子分割方法に依存しないため、本発明において各格子は任意の形状及び任意のサイズであってよい。ただし、全ての格子が均一の形状である場合、そうでない場合と比べ浸水深補間処理の計算負荷が軽くなる。また、格子間隔が粗い（すなわち各格子のサイズが大きい）場合、格子間隔が細かい場合と比べ計算負荷が軽くなる一方、得られる浸水深分布も粗くなる。したがって、計算負荷等と、得たい浸水深分布の空間分解能とを考慮して、格子の形状及びサイズを決めればよい。
【００７５】
例えば、浸水域を含む地域を東西方向及び南北方向に所定の間隔で分割することによって得られる区画を格子としてもよい。あるいは、住宅の敷地を格子としてもよい。
【００７６】
本実施の形態では、ＵＴＭ（Ｕｎｉｖｅｒｓａｌ  Ｔｒａｎｓｖｅｒｓｅ  Ｍｅｒｃａｔｏｒ）座標系において東西方向に５ｍ、南北方向に５ｍごとの格子に分割する方法を採った。これは、細かい空間スケールの分布に適したメッシュデータ表現である。また、原点からその格子までの東西方向に沿った距離、及び南北方向に沿った距離を並べた２数を格子ＩＤとした。この２数（ｉ，ｊ）は以下の式で算出すればよい。
【００７７】
ｉ＝ｆｌｏｏｒ（ｘ／ｄｘ）；  ｊ＝ｆｌｏｏｒ（ｙ／ｄｙ）；
ただし、ｘ及びｙは、それぞれ、ＵＴＭ座標におけるｘ座標値及びｙ座標値である。ｄｘ及びｄｙは、それぞれ、東西方向の格子のサイズ及び南北方向の格子のサイズである。ｉ及びｊは、それぞれ、原点から東西方向に沿って数えた格子の数、及び、原点から南北方向に沿って数えた格子の数である。ｆｌｏｏｒ（）は、引数の小数部分を切り捨てて整数を返す関数である。
【００７８】
なお、格子ＩＤは、各格子を一意に識別するものである限り、上記以外の方法によって定められてもよい。
【００７９】
ステップ１００１において、浸水深補間プログラム１４２は、浸水地点を含む格子についてステップ１００２を繰り返す。凡例１１０１にある通り、浸水地点は例えばアイコン８２４が表示された位置である。したがって、ステップ１００２における処理対象である格子は、例えば格子１１２１である。
【００８０】
ステップ１００２において、浸水深補間プログラム１４２は、処理対象の格子（例えば格子１１２１）における浸水位を算出する。浸水位の値は、浸水地点データ１９３の浸水深９０５に標高を加算することによって算出される。格子における標高は、ＨＤＤ１９０に格納されている標高データ１９５を参照することによって取得される。
【００８１】
ステップ１００３において、浸水深補間プログラム１４２は、浸水域の境界線（すなわち浸水域の輪郭線）を含む格子についてステップ１００４を繰り返す。図１１において浸水域の境界線は、境界線１１１１である。したがって、ステップ１００４で処理対象とする格子は、例えば格子１１２２である。なお、浸水域の境界線を含む格子について、浸水地点データ１９３が入力されている場合、その格子はすでにステップ１００２で処理されているため、ステップ１００４では処理されない。
【００８２】
ステップ１００４において、浸水深補間プログラム１４２は、格子における浸水位を算出する。浸水位の算出方法は、ステップ１００２と同様である。ただしステップ１００２とは異なり、浸水域の境界線では浸水深として０を与える。なお、浸水域の調査方法によっては、浸水域の境界線で浸水深が０とならず、若干の誤差を持つ場合も考えられる。この場合は、その誤差を標高に加えた値を浸水位とすればよい。
【００８３】
ステップ１００５において、浸水深補間プログラム１４２は、補間対象の格子を選択する。補間対象の格子は、浸水域の境界線より内側にあり、かつ浸水地点データを含まず、かつ浸水域の境界線を含まない格子である。図１１において補間対象の格子は例えば格子１１２３である。
【００８４】
以降、ステップ１００６から１０１１において、浸水深補間プログラム１４２は、ラプラス方程式ｄｉｖ（ｇｒａｄ（ｗｌｅｖ））＝０を解く。本実施例では均一の格子を用いているため、式（１）で示したラプラス方程式は、例えば以下の式（２）のように離散化される。
ｗｌｅｖ［ｉ，ｊ］＝１／６＊（（ｗｌｅｖ［ｉ−１，ｊ］＋ｗｌｅｖ［ｉ＋１，ｊ］＋ｗｌｅｖ［ｉ，ｊ−１］＋ｗｌｅｖ［ｉ，ｊ＋１］）＋１／２＊（ｗｌｅｖ［ｉ−１，ｊ−１］＋ｗｌｅｖ［ｉ＋１，ｊ−１］＋ｗｌｅｖ［ｉ−１，ｊ＋１］＋ｗｌｅｖ［ｉ＋１，ｊ＋１］））  ・・・式（２）
ただし、ｗｌｅｖ［ｉ，ｊ］は格子（ｉ，ｊ）における水位ｗｌｅｖを表す。式（２）を解くためには、平均化処理、あるいは、ＳＯＲ法（ｓｕｃｃｅｓｓｉｖｅ  ｏｖｅｒ−ｒｅｌａｘａｔｉｏｎ  ｍｅｔｈｏｄ）などの行列解法が利用できる。本実施の形態では、理解のしやすい平均化処理を使う方法をステップ１００６から１０１１に示す。
【００８５】
なお、式（１）の離散化方法は式（２）に限られない。式（２）では８近傍の格子を用いて離散化したが、例えば４近傍の格子を用いて離散化した以下の式（３）を用いてもよい。
ｗｌｅｖ［ｉ，ｊ］＝１／４＊（（ｗｌｅｖ［ｉ−１，ｊ］＋ｗｌｅｖ［ｉ＋１，ｊ］＋ｗｌｅｖ［ｉ，ｊ−１］＋ｗｌｅｖ［ｉ，ｊ＋１］））  ・・・式（３）
このように、本実施の形態において、式（１）で示したラプラス方程式を解くことは、処理対象の格子に隣接する複数の格子の浸水位に所定の計数（上記の例では、１／６、１／２、１／４等）を乗じた値を算出し、それら算出された複数の値を加算することによって実現される。
【００８６】
なお、８近傍の格子とは、処理対象の格子に隣接する８個の格子である。具体的には、処理対象の格子の北、北西、西、南西、南、南東、東及び北東に隣接する８個の格子である。一方、４近傍の格子とは、処理対象の格子に隣接する４個の格子である。具体的には、処理対象の格子の北、西、南及び東に隣接する４個の格子である。
【００８７】
ステップ１００６において、浸水深補間プログラム１４２は、選択された格子が存在する間ステップ１００７以降を繰り返す。
【００８８】
ステップ１００７において、浸水深補間プログラム１４２は、その時点で選択された格子についてステップ１００８以降を繰り返す。
【００８９】
ステップ１００８において、浸水深補間プログラム１４２は、処理対象の格子における浸水位を、式（２）を用いて算出する。あるいは、ステップ１００８において、上記の式（３）が使用されてもよい。
【００９０】
ステップ１００９において、浸水深補間プログラム１４２は、処理対象の格子の浸水位の変化量（すなわち、現時点で算出されている浸水位と、前回のステップ１００８が実行される前の浸水位との差の絶対値）を算出する。浸水深補間プログラム１４２は、この変化量が閾値以下である場合はステップ１０１０を、そうでない場合はステップ１０１１を実行する。ここで閾値は、求めたい補間精度よりも小さい値とする。閾値の具体的な値は経験的に決めればよい。発明者らは０．１ｍ程度の補間精度が妥当であると考え、閾値を０．０５ｍにした。
【００９１】
ステップ１０１０において、浸水深補間プログラム１４２は、処理対象の格子の選択を解除する。解除された格子は、再び処理対象として選択されない限り、以後ステップ１００８以降の処理がなされない。
【００９２】
ステップ１０１１において、浸水深補間プログラム１４２は、処理対象の格子の８近傍にある補間対象の格子を選択状態にする。ステップ１００８において式（２）が使用された場合、ステップ１０１１において、浸水深補間プログラム１４２は、処理対象の格子の４近傍にある補間対象の格子を選択状態にする。なお、処理対象の格子はすでに選択状態であるため、その後ステップ１０１０において選択を解除されるまで、繰り返しステップ１００８が実行される。
【００９３】
ステップ１０１２において、浸水深補間プログラム１４２は、このステップ以前のステップで浸水位が与えられた格子についてステップ１０１３を実行する。
【００９４】
ステップ１０１３において、浸水深補間プログラム１４２は、処理対象の格子における浸水位から浸水深を算出する。具体的には、処理対象の格子の浸水位からその格子の標高を減算することによって浸水深が算出される。
【００９５】
以上、図１０に示したステップによって浸水深補間処理（ステップ５０５）が完了する。これによって、ユーザが入力した浸水域形状及び浸水地点データを反映した浸水深分布データが生成された。
【００９６】
図１２は、本発明の第１の実施の形態の浸水深分布データ１９４のデータ構造を示す説明図である。
【００９７】
浸水深分布データ１９４は、少なくとも、災害を一意に識別する災害ＩＤ１２０１、浸水域を一意に識別する浸水域ＩＤ１２０２、格子を一意に識別する格子ＩＤ１２０３、及びその格子の浸水深１２０４の情報を含む。
【００９８】
浸水深補間プログラム１４２は、ステップ５０５で作成した浸水深分布データ１９４を、ＨＤＤ１９０に格納する。さらに、浸水深補間プログラム１４２は、ステップ５０４で受信した浸水域データ１９２及び浸水地点データ１９３をＨＤＤ１９０に格納する（ステップ５０６）。
【００９９】
さらに、浸水深補間プログラム１４２は、作成した浸水深分布データ１９４を入力支援プログラム１４１に送信する（ステップ５０７）。入力支援プログラム１４１は、受信した浸水深分布データを地図４１０に出力する（ステップ５０８）。
【０１００】
図１３は、本発明の第１の実施の形態の入力支援プログラム１４１が提供する第２のＧＵＩが浸水深分布を表示した状態を示す説明図である。
【０１０１】
図１３を図８と比較すると、ユーザがステップ５０２で入力した浸水地点にあるアイコンが、入力待ち状態のアイコン８２４から確定状態のアイコン１３２４に変化した。さらに、浸水深分布１３２５が表示された。なお、浸水域の境界線を含む格子（たとえば格子１３２２）については、格子全体をそのまま表示するのではなく、格子のうち浸水域に含まれる部分のみが表示される。その結果、ユーザが入力した浸水域形状と、最終的に表示される浸水域の形状とが等しくなるため、ユーザの直感に沿う。
【０１０２】
浸水深分布１３２５は、いかなる方法で表示されてもよい。例えば、浸水深分布１３２５は、各格子を浸水深１２０４に応じた色に塗り分けることによって浸水深１２０４を表示してもよい。具体的には、浸水深１２０４の値（又はその値の範囲）を所定の色と対応付け、各格子の浸水深１２０４に対応する色を浸水深分布１３２５に表示してもよい。あるいは、浸水深１２０４に対応する所定の模様又は図形等が表示されてもよい。凡例１３４１には、色（又は模様又は図形等）と浸水深１２０４との対応を記載することが望ましい。
【０１０３】
図１３の例では、０ｃｍから５０ｃｍまでの浸水深が右下がりの細い斜線によって表示され、５０ｃｍを超える浸水深が右上がりの太い斜線によって表示される。
【０１０４】
ステップ５０９において、入力支援プログラム１４１は、ユーザに浸水地点データの入力、あるいはプログラムの終了命令の入力を促し、ユーザからの入力を受け付ける。このステップは前述したステップ２０４でＣａｓｅ  Ａ、Ｃａｓｅ  Ｂのいずれの場合であっても省略できない。本実施の形態では、ユーザがステップ５０２と同じ操作、すなわち地図４１０の１以上の地点の浸水深を入力すると、入力支援プログラム１４１はステップ５０２以降の処理を再び繰り返す。一方、ユーザが終了命令ボタン４０３を押下すると、入力支援プログラム１４１は終了する。
【０１０５】
以上の本発明の第１の実施の形態では、図１に示すように、ＰＣ１００で入力支援プログラム１４１と浸水深補間プログラム１４２とを動作させるコンピュータシステムを示した。しかし、上記のコンピュータシステムと同様の処理は、クライアントコンピュータとサーバコンピュータとがインターネットで接続されているコンピュータシステムで実行されてもよい。この場合、サーバコンピュータは多数のクライアントコンピュータから浸水地点データを受け付けることが可能になる。
【０１０６】
図１４は、本発明の第２の実施の形態のコンピュータシステムの構成を示す概略ブロック図である。
【０１０７】
図１４の計算機システムは、クライアントコンピュータ１０１及びサーバコンピュータ１０２によって構成される。クライアントコンピュータ１０１およびサーバコンピュータ１０２は、インターネット１５０を介して接続されている。
【０１０８】
クライアントコンピュータ１０１は、入出力部１１０を備える。入出力部１１０は、第１の実施の形態において説明したものと同じ構成である。クライアントコンピュータ１０１は、さらに、メモリ１２１及びＣＰＵ１３１を備える。メモリ１２１に展開された入力支援プログラム１４１に基づいて、ＣＰＵ１３１が命令を実行する。したがって、入力支援プログラム１４１が実行する処理は、実際には、入力支援プログラム１４１に基づいて、ＣＰＵ１３１によって実行される。
【０１０９】
サーバコンピュータ１０２は、メモリ１２２及びＣＰＵ１３２を備える。メモリ１２２に展開された浸水深補間プログラム１４２に基づいて、ＣＰＵ１３２が命令を実行する。したがって、浸水深補間プログラム１４２が実行する処理は、実際には、浸水深補間プログラム１４２に基づいて、ＣＰＵ１３２によって実行される。この浸水深補間プログラム１４２は、ＨＤＤ１９０に格納された地図データ１９１、浸水域データ１９２、浸水地点データ１９３、浸水深分布データ１９４、標高データ１９５、および災害データ１９６を必要に応じて読み出す。ＨＤＤ１９０は、サーバコンピュータ１０２に内蔵されてもよいし、サーバコンピュータ１０２の外部に接続されてもよい。
【０１１０】
図１４は、クライアントコンピュータ１０１とサーバコンピュータ１０２とがインターネット１５０を介して接続される構成を示すが、インターネット１５０は、他のいかなる種類のネットワークによって置き換えられてもよい。
【０１１１】
第２の実施の形態の入力支援プログラム１４１及び浸水深補間プログラム１４２は、第１の実施の形態の入力支援プログラム１４１及び浸水深補間プログラム１４２と同様の処理を実行する。このため、第２の実施の形態においてこれらのプログラムが実行する処理についての説明は省略する。ただし、第２の実施の形態において、入力支援プログラム１４１と浸水深補間プログラム１４２との間で送受信されるデータは、インターネット１５０を経由する。
【０１１２】
本発明の第１及び第２の実施の形態によれば、粗い間隔で計測された浸水深に基づいて、細かい浸水深分布を推定することができる。さらに、本発明の第２の実施の形態によれば、インターネットを通じて浸水地点データを収集し、そのデータに基づいて浸水深分布を推定することができる。このため、一般市民が計測した水害に関するデータを有効利用することができる。このため、従来行われている浸水実績の現地調査にかかる時間と手間を軽減することができる。
【図面の簡単な説明】
【０１１３】
【図１】本発明の第１の実施の形態のコンピュータシステムの構成を示す概略ブロック図である。
【図２】本発明の第１の実施の形態のＰＣによって実行される第１の処理を示すシーケンス図である。
【図３】本発明の第１の実施の形態の入力支援プログラムが提供する第１のＧＵＩを示す説明図である。
【図４】本発明の第１の実施の形態の入力支援プログラムが提供する第２のＧＵＩの初期状態を示す説明図である。
【図５】本発明の第１の実施の形態のＰＣによって実行される第２の処理を示すシーケンス図である。
【図６】本発明の第１の実施の形態において、浸水域を入力中の状態にある第２のＧＵＩを示す説明図である。
【図７】本発明の第１の実施の形態の浸水域データのデータ構造を示す説明図である。
【図８】本発明の第１の実施の形態において、浸水地点を入力中の状態にある第２のＧＵＩを示す説明図である。
【図９】本発明の第１の実施の形態の浸水地点データのデータ構造を示す説明図である。
【図１０】本発明の第１の実施の形態の浸水深補間処理の詳細を示すＰＡＤ図である。
【図１１】本発明の第１の実施の形態における浸水地点、浸水域及び格子のそれぞれの位置関係を示す説明図である。
【図１２】本発明の第１の実施の形態の浸水深分布データのデータ構造を示す説明図である。
【図１３】本発明の第１の実施の形態の入力支援プログラムが提供する第２のＧＵＩが浸水深分布を表示した状態を示す説明図である。
【図１４】本発明の第２の実施の形態のコンピュータシステムの構成を示す概略ブロック図である。
【符号の説明】
【０１１４】
１００  ＰＣ
１４１  入力支援プログラム
１４２  浸水深補間プログラム
１９２  浸水域データ
１９３  浸水地点データ
１９４  浸水深分布データ
１９５  標高データ
【図１】
0005403726.tif 000002
【図２】
0005403726.tif 000003
【図９】
0005403726.tif 000004
【図１０】
0005403726.tif 000005
【図３】
0005403726.tif 000006
【図４】
0005403726.tif 000007
【図５】
0005403726.tif 000008
【図６】
0005403726.tif 000009
【図７】
0005403726.tif 000010
【図８】
0005403726.tif 000011
【図１１】
0005403726.tif 000012
【図１２】
0005403726.tif 000013
【図１３】
0005403726.tif 000014
【図１４】
0005403726.tif 000015