JP2023122016A - 情報処理装置、情報処理プログラム及び情報処理方法 - Google Patents

Abstract

【課題】使用するリモートセンシングデータの画質及び撮像条件に関わらず、高精度に浸水域を推定可能な装置を提供する。【解決手段】情報処理装置１０は、第１画像＃１と、第１画像＃１よりも過去の日時に撮影された第２画像＃２との比較から、初期浸水域１１Ａを推定する初期浸水域推定部１１と、地形の標高を示す標高情報３０について地形勾配を解析することにより得られた窪地領域１２Ａに初期浸水域１１Ａの少なくとも一部が含まれる場合に、窪地領域１２Ａの浸水の程度を示す指標値を計算し、指標値が閾値以上の場合に、窪地領域１２Ａを浸水と判断し、浸水を示す判断結果に基づいて、窪地領域１２Ａを浸水域で充填する。【選択図】図１

Description

本開示は、情報処理装置、情報処理プログラム及び情報処理方法に関する。本開示は特に、水害発生時に家屋又は山林の浸水状況を推定する情報処理に関する。

洪水、大雨、台風といった水害発生時に、浸水被害を把握する手法として、従来は地上センサが使用されていた。地上センサを用いた浸水把握の手法は、エレベータ又は発電所といった特定のスポットの浸水状況把握には有用である。しかし、地上センサを用いた浸水把握の手法は、被害が広域に発生する状況において、どの地域がより被害が大きいかといった面方向の広がりを持った情報を取得することには向いていない。

広域の浸水被害を把握する手法として、リモートセンシングデータを活用した浸水深推定手法が提案されている（非特許文献１参照）。非特許文献１では、災害前後で取得された二時期のＳＡＲ衛星画像を比較して変化領域を検出し、浸水の範囲である浸水域を推定している。

"衛星ＳＡＲを用いた豪雨災害時の迅速な浸水範囲及び浸水深抽出の試み．"本田謙一ら， Ｊｏｕｒｎａｌ ｏｆ Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ Ｓｏｃｉｅｔｙ ｏｆ Ｊａｐａｎ，Ｖｏｌ．４０ Ｎｏ．３（２０２０），ｐｐ．１５８－１６２．

非特許文献１では、浸水域は有る時刻の第１の画像と、第１の画像に対して過去の時刻の第２の画像との２つのＳＡＲ衛星画像との比較から浸水域を算出する。一般的には２つの画像間で差が顕著に発生している箇所を浸水域とみなす。この場合、画像中のノイズ及び季節の変動の影響から、差が疎らに発生するケースが存在する。このように画像に依存した浸水域検出方法の場合、画質や撮像条件が浸水域の判定に影響を与える。このため、実際に浸水している領域であっても、その領域の画像間に差が発生していないことで、その領域が非浸水域と見做される可能性がある。これは、推定精度の面で課題である。

本開示は、使用するリモートセンシングデータの画質及び撮像条件に関わらず、高精度に浸水域を推定することを目的とする。

本開示の情報処理装置は、
第１画像と、前記第１画像よりも過去の日時に撮影された第２画像との比較から、前記第１画像と前記第２画像との両方に存在する領域であって浸水している領域を示す初期浸水域を推定する推定部と、
地形の標高を示す標高情報について地形勾配を解析することにより得られた窪地領域に前記初期浸水域の少なくとも一部が含まれる場合に、前記窪地領域の浸水の程度を示す指標値を計算し、前記指標値が閾値以上の場合に、前記窪地領域を浸水と判断し、浸水を示す判断結果に基づいて、前記窪地領域を浸水域で充填する浸水域充填部と、
を備える。

本開示の情報処理装置によれば、使用するリモートセンシングデータの画質あるいは撮像条件に関わらず、高精度に浸水域を推定可能である。

実施の形態１の図で、情報処理装置１０１のブロック図。 実施の形態１の図で、情報処理装置１０１のハードウェア構成を示す図。 実施の形態１の図で、情報処理方法の動作を示すフローチャート。 実施の形態１の図で、浸水窪地候補１３１Ａと初期浸水域１１Ａとを示す図。 実施の形態１の図で、変形例における情報処理装置２０のブロック図。 実施の形態１の図で、情報処理装置１０，２０のハードウェア構成を補足する図。

実施の形態の説明及び図面において、同じ要素及び対応する要素には同じ符号を付している。同じ符号が付された要素の説明は、適宜に省略又は簡略化する。以下の実施の形態では、「部」を、「回路」、「工程」、「手順」、「処理」又は「サーキットリ」に適宜読み替えてもよい。

実施の形態１．
図１は、情報処理装置１０の機能構成図である。
図２は、情報処理装置１０１のハードウェア構成を示す。

＊＊＊構成の説明＊＊＊
図１に示すように、情報処理装置１０１は、初期浸水域推定部１１、窪地検出部１２、浸水域充填部１３を備えている。浸水域充填部１３は、候補検出部１３１、浸水判断部１３２、浸水窪地充填部１３３を備えている。
初期浸水域推定部１１は推定部である。

情報処理装置１０は、コンピュータである。情報処理装置１０は、プロセッサ５１０を備える。情報処理装置１０は、プロセッサ５１０の他に複数のハードウェアを備える。複数のハードウェアは、主記憶装置５２０、補助記憶装置５３０、入力インタフェース５４０、出力インタフェース５５０及び通信インタフェース５６０である。プロセッサ５１０は、信号線５７０を介して、他のハードウェアと接続され、他のハードウェアを制御する。

情報処理装置１０は、機能要素として、初期浸水域推定部１１、窪地検出部１２、浸水域充填部１３を備える。
初期浸水域推定部１１、窪地検出部１２、浸水域充填部１３の機能は、情報処理プログラム５３１により実現される。

プロセッサ５１０は、情報処理プログラム５３１を実行する装置である。プロセッサ５１０が情報処理プログラム５３１を実行することで、初期浸水域推定部１１、窪地検出部１２、浸水域充填部１３の機能が実現される。プロセッサ５１０は、演算処理を行うＩＣ（Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）である。プロセッサ５１０の具体例は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）である。

主記憶装置５２０は記憶装置である。主記憶装置５２０の具体例は、ＳＲＡＭ（Ｓｔａｔｉｃ Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、ＤＲＡＭ（Ｄｙｎａｍｉｃ Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）である。主記憶装置５２０は、プロセッサ５１０の演算結果を保持する。

補助記憶装置５３０は、データを不揮発的に保管する記憶装置である。補助記憶装置５３０の具体例は、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）である。また、補助記憶装置５３０は、可搬記録媒体であってもよい。可搬記録媒体として、ＳＤ（登録商標）（Ｓｅｃｕｒｅ Ｄｉｇｉｔａｌ）メモリカード、ＮＡＮＤフラッシュ、フレキシブルディスク、光ディスク、コンパクトディスク、ブルーレイ（登録商標）ディスク、ＤＶＤ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｋ）がある。補助記憶装置５３０は、情報処理プログラム５３１を記憶している。

入力インタフェース５４０は、各装置からデータが入力されるポートである。出力インタフェース５５０は、各種機器が接続される。出力インタフェース５５０は、各種機器にプロセッサ５１０によってデータが出力されるポートである。通信インタフェース５６０は、プロセッサが他の装置と通信するための通信ポートである。

プロセッサ５１０は補助記憶装置５３０から情報処理プログラム５３１を主記憶装置５２０にロードする。プロセッサ５１０は、ロードされた情報処理プログラム５３１を主記憶装置５２０から読み込んで実行する。主記憶装置５２０には、情報処理プログラム５３１の他に、ＯＳ（Ｏｐｅｒａｔｉｎｇ Ｓｙｓｔｅｍ）も記憶されている。プロセッサ５１０は、ＯＳを実行しながら、情報処理プログラム５３１を実行する。情報処理装置１０は、プロセッサ５１０を代替する複数のプロセッサを備えてもよい。これら複数のプロセッサは、情報処理プログラム５３１の実行を分担する。それぞれのプロセッサは、プロセッサ５１０と同じように、情報処理プログラム５３１を実行する装置である。情報処理プログラム５３１によって利用、処理または出力されるデータ、情報、信号値及び変数値は、主記憶装置５２０、補助記憶装置５３０、または、プロセッサ５１０内のレジスタあるいはキャッシュメモリに記憶される。

情報処理プログラム５３１は、初期浸水域推定部１１、窪地検出部１２、浸水域充填部１３、候補検出部１３１、浸水判断部１３２、浸水窪地充填部１３３の「部」を「処理」、「手順」あるいは「工程」に読み替えた各処理、各手順あるいは各工程を、コンピュータに実行させるプログラムである。

また、方法は、コンピュータである情報処理装置１０が情報処理プログラム５３１を実行することにより行われる情報処理方法である。情報処理プログラム５３１は、コンピュータが読み取り可能な記録媒体に格納されて提供されてもよいし、プログラムプロダクトとして提供されてもよい。

＊＊＊動作の説明＊＊＊
図３を参照して、情報処理装置１０の動作を説明する。
図３は、情報処理装置１０の動作を説明するフローチャートである。情報処理装置１０の動作は、情報処理方法に相当する。また情報処理装置１０の動作は、情報処理プログラムの処理に相当する。

＜ステップＳ１＞
ステップＳ１では、第１画像＃１と第２画像＃２とから、初期浸水域推定部１１が初期浸水域１１Ａを推定する。
初期浸水域推定部１１と、第１画像＃１と、第１画像＃１よりも過去の日時に撮影された第２画像＃２との比較から、第１画像＃１と第２画像＃２との両方に存在する領域であって浸水している領域を示す初期浸水域１１Ａを推定する。以下に具体的に説明する。
第１画像＃１と第２画像＃２とは、飛翔体に搭載された撮影装置によって撮影されたリモートセンシングデータである。飛翔体とは、航空機、ＵＡＶの他、人工衛星である。初期浸水域推定部１１は、入力データとして、第１画像＃１と第２画像＃２とが入力される。第１画像＃１は、ある日時を示す時期Ｔ１に取得された、リモートセンシングデータである。第２画像＃２は、ある時期Ｔ１よりも古いある日時を示す時期Ｔ２に取得されたリモートセンシングデータである。第２画像＃２と第１画像＃１とは、同一の領域を含む。第２画像＃２と第１画像＃１とは、同一領域の差分が比較される。初期浸水域推定部１１は、第１画像＃１と第２画像＃２とを比較し、比較結果から、初期浸水域１１Ａを推定する。
ここで「初期浸水域１１Ａ」における「初期」とは、浸水域充填部１３へ入力される前の「浸水域」を意味する。浸水域充填部１３による補正処理が実施された「最終浸水域１３３Ａ」との対比において、補正処理が実施前を意味する「初期」を付している。初期浸水域推定部１１は、推定結果である初期浸水域１１Ａを、浸水域充填部１３の候補検出部１３１へ出力する。初期浸水域推定部１１による浸水の有無の推定は、例えば、以下のように（１）光学画像、または（２）ＳＡＲ画像、から抽出可能である。

（１）光学画像
光学画像から初期浸水域１１Ａを推定する場合、初期浸水域推定部１１は、近赤外帯域の画像を利用した水指標であるＮＤＷＩ（Ｎｏｒｍａｌｉｚｅｄ Ｄｉｆｆｅｒｅｎｃｅ ｗａｔｅｒ ｉｎｄｅｘ）を、用いることができる。ρ_ＮＩＲ、ρ_ＳＷＩＲを、それぞれある地点のＮＩＲ（Ｎｅａｒ Ｉｎｒｅａｒｅｄ）、ＳＷＩＲ（Ｓｈｏｒｔ―Ｗａｖｅ Ｉｎｆｒａｒｅｄ）帯域の計測値とする。その場合、ＮＤＷＩは以下に示す式（１）で表される。初期浸水域推定部１は、第１画像＃１と第２画像＃２との間でのＮＤＷＩの変化を検出することで、初期浸水域１１Ａを推定可能である。
ＮＤＷＩ＝［ρ_ＮＩＲ－ρ_ＳＷＩＲ］／［ρ_ＮＩＲ＋ρ_ＳＷＩＲ］ 式（１）
近年、赤外帯域の代わりに可視帯域を利用した浸水域の推定手法も提案されている。
初期浸水域推定部１１は、可視帯域を利用した手法で初期浸水域１１Ａを推定することもできる。
（２）ＳＡＲ画像
ＳＡＲ画像から初期浸水域１１Ａを推定する場合、非特許文献１の手法は、水害発生前後のＳＡＲ画像を比較し、後方散乱強度が低下した領域を浸水域として抽出する。初期浸水域推定部１１は、非特許文献１の手法を用いても良い。

以上に説明した（１）光学画像、（２）ＳＡＲ画像の手法は、リモートセンシングデータの一例である。本明細書の情報処理装置に対しては、センサと計測装置との組み合わせにより、様々なリモートセンシングデータを入力できる。

＜ステップＳ２＞
ステップＳ２では、窪地検出部１２が標高情報３０を取得し、取得した標高情報３０から窪地領域１２Ａを検出する。窪地検出部１２には、後述の標高情報３０が入力される。窪地検出部１２は、標高情報３０を解析し、解析範囲内から窪地領域１２Ａを検出する。

＜標高情報３０＞
標高情報３０としては、リモートセンシング技術または測量によって作成されたＤＥＭ（Ｄｉｇｉｔａｌ Ｅｌｅｖａｔｉｏｎ Ｍｏｄｅｌ）情報が利用可能である。ＤＥＭは、レーザ測量または写真測量などの手法で作成される。ＤＥＭは一般ユーザが取得可能なオープンデータである。一般ユーザが取得可能なオープンデータの例として、「国土地理院：基盤地図情報 ｈｔｔｐｓ：／／ｆｇｄ．ｇｓｉ．ｇｏ．ｊｐ／ｄｏｗｎｌｏａｄ／ｍｅｎｕ．ｐｈｐ」がある。国土地理院が５ｍメッシュＤＥＭ５Ａを公開している。全世界ではＳＲＴＭ－３（Ｓｈｕｔｔｌｅ Ｒａｄｏｒ Ｔｏｐｏｇｒａｐｈｙ Ｍｉｓｓｉｏｎ）が公開されているほか、垂直分解能の高いＤＥＭデータも入手可能である。

＜窪地領域１２Ａ＞
窪地検出部１２は、窪地を勾配のない空間的に接続された領域と見做す。これにより、窪地検出部１２は、窪地領域１２Ａを、標高情報３０から検出する。窪地検出部１２は、標高情報３０の２次元数値勾配を算出し、勾配の絶対値が予め設定した閾値以下となる座標を算出する。窪地検出部１２は、算出した座標で構成される連結要素のうち、面積が第１閾値以上、第２閾値未満となる領域を、窪地領域１２Ａとして検出する。

もしくは、別の実装例として、窪地検出部１２は、以下の処理で窪地領域１２Ａを検出する。窪地検出部１２は、標高情報３０の２次元数値勾配から各画素の流向を算出し、流向がいずれの方向にも当てはまらない座標を算出する。窪地検出部１２は、算出した座標で構成される連結要素のうち、面積が第１閾値以上、第２閾値未満となる領域を、窪地領域１２Ａとして検出する。

＜ステップＳ３＞
ステップＳ３では、浸水域充填部１３は、地形の標高を示す標高情報３０について地形勾配を解析することにより得られた窪地領域に初期浸水域１１Ａの少なくとも一部が含まれる場合に、窪地領域の浸水の程度を示す指標値を計算する。この窪地領域は、後述の浸水窪地候補である。また指標値の例は、後述のように、割合Ｒまたは面積Ａである。浸水域充填部１３は、指標値が閾値以上の場合に、窪地領域を浸水と判断し、浸水を示す判断結果に基づいて、窪地領域を浸水域で充填する。以下に詳しく説明する。
浸水域充填部１３は、初期浸水域１１Ａと、窪地領域１２Ａとから、浸水している窪地領域を判断する。浸水域充填部１３は、浸水していると判断した窪地領域で初期浸水域１１Ａを充填することで、初期浸水域１１Ａを高精度化した最終浸水域１３３Ａを算出する。浸水域充填部１３は、窪地検出部１２によって検出された窪地領域１２Ａに、ある割合以上で初期浸水域１１Ａが存在する場合には、該当の窪地領域全体を初期浸水域に充填して最終浸水域１３３Ａとして出力する。
以下に、浸水域充填部１３の備える候補検出部１３１、浸水判断部１３２、浸水窪地充填部１３３の動作を説明する。

＜候補検出部１３１＞
候補検出部１３１は、初期浸水域１１Ａが少しでも含まれる窪地領域を、浸水窪地候補１３１Ａとして検出する。

＜浸水判断部１３２＞
浸水判断部１３２は、浸水窪地候補１３１Ａが浸水か非浸水かを判断し、判断結果である浸水判断結果１３２Ａを、浸水窪地充填部１３３へ出力する。
浸水判断結果１３２Ａは、「浸水」と「非浸水」とのどちらかを示す。
以下に具体的に説明する。
図４は、浸水窪地候補１３１Ａと初期浸水域１１Ａとを示す。図０４では、浸水窪地候補１３１Ａに初期浸水域１１Ａの少なくとも一部が含まれている。
浸水判断部１３２は、浸水窪地候補１３１Ａの面積の中に、どれくらいの割合で、もしくはどれくらいの面積の、初期浸水域１１Ａが含まれるかを算出する。浸水判断部１３２は以下の（１）（２）（３）のいずれかによって、「浸水窪地候補１３１Ａが浸水している」と判断することができる。
（１）浸水判断部１３２は、算出した割合が、割合の閾値ＴＨｒ以上の場合、「浸水窪地候補１３１Ａが浸水している」と判断する。ここで割合を割合Ｒとすれば、
割合Ｒ＝［初期浸水域１１Ａの面積］／［浸水窪地候補１３１Ａの面積］
において、
Ｒ≧ＴＨｒ
で、浸水判断部１３２は、「浸水」と判断する。
（２）浸水判断部１３２は、算出した面積が、面積の閾値ＴＨａ以上の場合、「浸水窪地候補１３１Ａが浸水している」と判断する。
算出した面積は、浸水窪地候補１３１Ａの面積の中における、初期浸水域１１Ａの面積である。
初期浸水域１１Ａの面積を面積Ａすれば、
Ａ≧ＴＨａ
で、浸水判断部１３２は「浸水」と判断する。
（３）浸水判断部１３２は、算出した割合Ｒが閾値ＴＨｒ以上、かつ、算出した面積Ａが閾値ＴＨａ以上の場合に、浸水窪地候補１３１Ａが浸水していると判断する。
Ｒ≧ＴＨｒ、かつ、Ａ≧ＴＨａ
で、浸水判断部１３２は、「浸水」と判断する。
なお、候補検出部１３１は、（１）（２）（３）のいずれかで判断した場合に、該当しなかった浸水窪地候補１３１Ａは、「非浸水」と判断する。

また、浸水窪地候補１３１Ａの面積中、一定未満の割合で、もしくは一定未満の面積で、初期浸水域１１Ａが含まれる場合、「浸水窪地候補１３１Ａが浸水していない」と判断してもよい。すなわち、浸水判断部１３２は、以下の（４）（５）（６）のいずれかによって、「浸水窪地候補１３１Ａが浸水していない」と判断してもよい。逆に言えば、浸水判断部１３２が以下の（４）（５）（６）のいずれかで判断した場合に、該当しなかった浸水窪地候補１３１Ａは、「浸水」と判断する。
（４）浸水判断部１３２は、算出した割合Ｒが、閾値ＴＨｒ未満の場合、「浸水窪地候補１３１Ａが浸水していない」と判断する。
Ｒ＜ＴＨｒ
で、浸水判断部１３２は「非浸水」と判断する。
（５）浸水判断部１３２は、算出した面積Ａが、面積の閾値ＴＨａ未満の場合、「浸水窪地候補１３１Ａが浸水していない」と判断する。
Ａ＜ＴＨａ
で、浸水判断部１３２は、「非浸水」と判断する。
（６）浸水判断部１３２は、算出した割合Ｒが閾値ＴＨｒ未満、かつ、算出した面積Ａが閾値ＴＨａ未満の場合に、「浸水窪地候補１３１Ａが浸水していない」と判断する。
Ｒ＜ＴＨｒ、かつ、Ａ＜ＴＨａ
で、浸水判断部１３２は、「非浸水」と判断する。

＜浸水窪地充填部１３３＞
浸水窪地充填部１３３は、浸水判断部１３２の判断結果に基づいて、窪地領域である浸水窪津候補を浸水域で充填する。以下に説明する。
（１）浸水域充填部１３の浸水窪地充填部１３３は、「浸水」を示す浸水判断結果１３２Ａに基づいて、窪地領域である浸水窪地候補１３１Ａに含まれる初期浸水域１１Ａの一部を浸水域として用いることにより、浸水窪地候補１３１Ａを浸水域で充填する。浸水窪地充填部１３３は、この充填結果を最終浸水域１３３Ａとして出力する。
（２）あるいは、浸水域充填部１３は、「浸水」を示す浸水判断結果１３２Ａに基づいて、
浸水窪地候補１３１Ａの全体を浸水域として用いることにより、浸水窪地候補１３１Ａを浸水域で充填する。浸水窪地充填部１３３は、この充填結果を最終浸水域１３３Ａとして出力する。
（３）また、浸水域充填部１３は、指標値が閾値未満の場合に、浸水窪地候補１３１Ａを「非浸水」と判断し、「非浸水」を示す浸水判断結果１３２Ａに基づいて、浸水窪地候補１３１Ａに含まれる初期浸水域１１Ａの一部を初期浸水域１１Ａから除外する。浸水窪地充填部１３３は、この充填結果を反映した初期浸水域１１Ａを、最終浸水域１３３Ａとして出力する。

＊＊＊実施の形態１の効果の説明＊＊＊
以上のように、情報処理装置１０によれば、入力画像である第１画像＃１、第２画像＃２の画質と撮像条件依存との少なくもいずれかで発生する、疎らな初期浸水域１１Ａを補正し、高精度に浸水域を推定することができる。即ち、使用するリモートセンシングデータの画質及び撮像条件に関わらず、高精度に浸水域を推定できる。

なお、実施の形態１の情報処理装置１０では、窪地検出部１２で検出された窪地領域１２Ａを用いて、初期浸水域１１Ａの改善処理を行っているが、このような例に限定されない。例えば、情報処理装置１０は、情報処理装置１０とは異なる他のシステムあるいは装置のような外部システムから、窪地領域１２Ａを取得し、取得した窪地領域１２Ａを浸水域充填部１３が用いることにより、最終浸水域１３３Ａを算出してもよい。この構成を変形例として以下に説明する。

＜変形例＞
図５は、情報処理装置１０の変形例である情報処理装置２０の機能構成図を示す。図５は図１に対応する。情報処理装置２０のハードウェア構成は、図２の情報処理装置１０のハードウェア構成と同様である。情報処理装置２０は、外部システムから、窪地領域を取得し、浸水域充填部１３に入力される点以外は情報処理装置１０と同じ構成である。

変形例の窪地領域１２Ａは、例えば、自治体が公開している窪地情報を利用可能である。自治体が公開している窪地情報は、地形情報から、周囲と比較して地盤面が低い領域を手作業で検出したものがある。情報処理装置１０のように窪地検出部１２が、標高情報３０から自動で窪地領域１２Ａを算出する場合と比較して、手作業で検する場合、高精度に窪地を検出できている場合がある。

＊＊＊変形例の効果の説明＊＊＊
情報処理装置２０によれば、標高情報３０を利用する必要が無く、高精度に浸水域を推定できる点で、解析の迅速性が高くなる。また、情報処理装置２０によれば、利用する窪地情報の精度によっては、更に高精度に浸水域を推定可能である。

＜ハードウェア構成の補足＞
情報処理装置１０、２０では情報処理装置１０の機能がソフトウェアで実現される場合を図２で述べた。しかし、情報処理装置１０、２０の機能がハードウェアで実現されてもよい。
図６は、情報処理装置１０，２０の機能がハードウェアで実現される構成を示す。以下では情報処理装置１０を例に説明する。図６の電子回路９０は、情報処理装置１０の、初期浸水域推定部１１、窪地検出部１２、浸水域充填部１３の機能を実現する専用の電子回路である。電子回路９０は、信号線９１に接続している。電子回路９０は、具体的には、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ロジックＩＣ、ＧＡ、ＡＳＩＣ、または、ＦＰＧＡである。ＧＡは、Ｇａｔｅ Ａｒｒａｙの略語である。ＡＳＩＣは、Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔの略語である。ＦＰＧＡは、Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙの略語である。情報処理装置１０の構成要素の機能は、１つの電子回路で実現されてもよいし、複数の電子回路に分散して実現されてもよい。また、情報処理装置１０の構成要素の一部の機能が電子回路で実現され、残りの機能がソフトウェアで実現されてもよい。

プロセッサ５１０と電子回路９０の各々は、プロセッシングサーキットリあるいはサーキットリとも呼ばれる。情報処理装置１０において、初期浸水域推定部１１、窪地検出部１２、浸水域充填部１３の機能がサーキットリにより実現されてもよい。

以上、実施の形態１について説明した。実施の形態１に含まれる複数の技術事項のうち、１つの技術事項を部分的に実施しても構わない。あるいは、実施の形態１に含まれる技術事項どうしを、部分的に組み合わせて実施しても構わない。

１０，２０ 情報処理装置、１１ 初期浸水域推定部、１１Ａ 初期浸水域、１２ 窪地検出部、１２Ａ 窪地領域、１３ 浸水域充填部、１３１ 候補検出部、１３１Ａ 浸水窪地候補、１３２ 浸水判断部、１３２Ａ 浸水判断結果、１３３ 浸水窪地充填部、１３３Ａ 最終浸水域、３０ 標高情報、５１０ プロセッサ、５２０ 主記憶装置、５３０ 補助記憶装置、５４０ 入力インタフェース、５５０ 出力インタフェース、５６０ 通信インタフェース、５７０ 信号線、９０ 電子回路、９１ 信号線。

Claims (8)

1. 第１画像と、前記第１画像よりも過去の日時に撮影された第２画像との比較から、前記第１画像と前記第２画像との両方に存在する領域であって浸水している領域を示す初期浸水域を推定する推定部と、
   地形の標高を示す標高情報について地形勾配を解析することにより得られた窪地領域に前記初期浸水域の少なくとも一部が含まれる場合に、前記窪地領域の浸水の程度を示す指標値を計算し、前記指標値が閾値以上の場合に、前記窪地領域を浸水と判断し、浸水を示す判断結果に基づいて、前記窪地領域を浸水域で充填する浸水域充填部と、
   を備える情報処理装置。
2. 前記浸水域充填部は、
   前記窪地領域に含まれる前記初期浸水域の一部を前記浸水域として用いることにより、前記窪地領域を充填する請求項１に記載の情報処理装置。
3. 前記浸水域充填部は、
   前記窪地領域の全体を前記浸水域として用いることにより、前記窪地領域を充填する請求項１に記載の情報処理装置。
4. 前記浸水域充填部は、
   前記指標値が閾値未満の場合に、前記窪地領域を非浸水と判断し、非浸水を示す判断結果に基づいて、前記窪地領域に含まれる前記初期浸水域の一部を前記初期浸水域から除外する請求項１から請求項３のいずれか１項に記載の情報処理装置。
5. 前記情報処理装置は、さらに、
   前記標高情報を取得し、取得した前記標高情報から窪地領域を検出する窪地検出部を備える請求項１から請求項４のいずれか１項に記載の情報処理装置。
6. 前記第１画像と前記第２画像とは、
   飛翔体に搭載された撮影装置によって撮影されたリモートセンシングデータである請求項１から請求項５のいずれか１項に記載の情報処理装置。
7. コンピュータに、
   第１画像と、前記第１画像よりも過去の日時に撮影された第２画像との比較から、前記第１画像と前記第２画像との両方に存在する領域であって浸水している領域を示す初期浸水域を推定する推定処理と、
   地形の標高を示す標高情報について地形勾配を解析することにより得られた窪地領域に前記初期浸水域の少なくとも一部が含まれる場合に、前記窪地領域の浸水の程度を示す指標値を計算し、前記指標値が閾値以上の場合に、前記窪地領域を浸水と判断し、浸水を示す判断結果に基づいて、前記窪地領域を浸水域で充填する浸水域充填処理と、
   を実行させる情報処理プログラム。
8. コンピュータが、
   第１画像と、前記第１画像よりも過去の日時に撮影された第２画像との比較から、前記第１画像と前記第２画像との両方に存在する領域であって浸水している領域を示す初期浸水域を推定し、
   地形の標高を示す標高情報について地形勾配を解析することにより得られた窪地領域に前記初期浸水域の少なくとも一部が含まれる場合に、前記窪地領域の浸水の程度を示す指標値を計算し、前記指標値が閾値以上の場合に、前記窪地領域を浸水と判断し、浸水を示す判断結果に基づいて、前記窪地領域を浸水域で充填する、
   情報処理方法。

Images