(57)【要約】【課題】水中にいかなる物体をも設置することなく、河川等を被写体とする映像信号のみから流水領域を検出できるようにする。【解決手段】固定カメラＣＭを用いて撮影された流水領域及び非流水領域が含まれる映像を複数フレームにわたって加算するフレーム加算部２と、加算された映像のエッジ成分を強調させる空間フィルタ処理部３と、空間フィルタ処理部３により生成された映像における各画素をエッジ画素か否かに分類する二値化処理部４と、得られる二値画像をもとに流水領域及び非流水領域を判別する領域分割部５Ａとを備え、水中にいかなる物体をも設置することなく、河川等を被写体とする映像信号のみから流水領域を検出できるようにする。【選択図】図１
2007256254.tif 000002
【特許請求の範囲】
【請求項１】
固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像が供給され、供給された映像を時間方向に積分する積分手段と、
上記積分手段により得られた映像のエッジ成分を強調させる空間フィルタ処理手段と、
上記空間フィルタ処理手段で処理された映像における各画素に対してエッジ画素であるか否かを画素毎に判定し、当該判定結果に基づいて上記流水領域及び上記非流水領域を判別する領域判別手段とを備えることを特徴とする流水領域検出システム。
【請求項２】
上記領域判別手段は、
上記空間フィルタ処理手段で処理された映像に対し、画素毎にエッジ画素であるか否かの判定を行ってエッジ画素又は非エッジ画素に分類する二値化処理手段と、
上記二値化処理手段により分類された上記エッジ画素又は上記非エッジ画素の分布に基づいて、映像内領域を上記流水領域及び上記非流水領域に二分割する領域分割手段とを備えることを特徴とする請求項１記載の流水領域検出システム。
【請求項３】
上記領域分割手段は、上記エッジ画素又は上記非エッジ画素を対象画素として、当該対象画素数を上記流水領域と非流水領域との境界が伸びる方向に関して積算し、さらに得られた対象画素数の分布から境界画素数を決定し、当該境界画素数をしきい値として上記得られた対象画素数に基づき上記流水領域及び上記非流水領域に分割することを特徴とする請求項２記載の流水領域検出システム。
【請求項４】
固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像が供給され、供給された映像を時間方向に積分する積分手段と、
上記積分手段により得られた映像のエッジ成分を強調させる空間フィルタ処理手段と、
上記空間フィルタ処理手段で処理された映像における各ライン毎に求めた特徴量に基づいて上記流水領域及び上記非流水領域を判別する領域判別手段とを備えることを特徴とする流水領域検出システム。
【請求項５】
上記領域判別手段は、
上記空間フィルタ処理手段で処理された映像の各ライン毎に、上記特徴量として分散値又は絶対値和を算出する特徴量計算手段と、
上記特徴量計算手段により算出された特徴量に基づいて、映像内領域を上記流水領域及び上記非流水領域に二分割する領域分割手段とを備えることを特徴とする請求項４記載の流水領域検出システム。
【請求項６】
上記積分手段は、供給される映像を複数のフレームにわたって加算することを特徴とする請求項１〜５の何れか１項に記載の流水領域検出システム。
【請求項７】
上記積分手段により複数のフレームにわたって映像を加算する前に、当該映像に非線形処理を施すことを特徴とする請求項６記載の流水領域検出システム。
【請求項８】
上記固定カメラを用いて撮影された映像から上記流水領域及び非流水領域が含まれる領域映像を抽出して上記積分手段に供給する領域抽出手段をさらに備えることを特徴とする請求項１〜７の何れか１項に記載の流水領域検出システム。
【請求項９】
上記固定カメラを用いて撮影された映像から互いに領域の異なる複数の上記領域映像を抽出し、抽出した各領域映像毎に処理して得られる各領域映像での上記流水領域と上記非流水領域の判別結果に基づいて、上記撮影された映像における上記流水領域及び上記非流水領域を判別することを特徴とする請求項８記載の流水領域検出システム。
【請求項１０】
上記空間フィルタ処理手段で用いる空間フィルタは、上記流水領域に対するフィルタ処理により得られる結果の分散値が最小となるようにフィルタ係数値が計算され随時更新されることを特徴とする請求項１〜９の何れか１項に記載の流水領域検出システム。
【請求項１１】
固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分工程と、
上記積分工程で得られた映像のエッジ成分を強調させる空間フィルタ処理工程と、
上記空間フィルタ処理工程で処理された映像における各画素に対してエッジ画素であるか否かを画素毎に判定する判定工程と、
上記判定工程での判定結果に基づいて、上記流水領域及び上記非流水領域を判別する領域判別工程とを有することを特徴とする流水領域検出方法。
【請求項１２】
固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分工程と、
上記積分工程で得られた映像のエッジ成分を強調させる空間フィルタ処理工程と、
上記空間フィルタ処理工程で処理された映像の各ライン毎に特徴量を算出する特徴量計算工程と、
上記特徴量計算工程で算出された特徴量に基づいて、上記流水領域及び上記非流水領域を判別する領域判別工程とを有することを特徴とする流水領域検出方法。
【請求項１３】
固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分ステップと、
上記積分ステップで得られた映像のエッジ成分を強調させる空間フィルタ処理ステップと、
上記空間フィルタ処理ステップで処理された映像における各画素に対してエッジ画素であるか否かを画素毎に判定する判定ステップと、
上記判定ステップでの判定結果に基づいて、上記流水領域及び上記非流水領域を判別する領域判別ステップとをコンピュータに実行させるためのプログラム。
【請求項１４】
固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分ステップと、
上記積分ステップで得られた映像のエッジ成分を強調させる空間フィルタ処理ステップと、
上記空間フィルタ処理ステップで処理された映像の各ライン毎に特徴量を算出する特徴量計算ステップと、
上記特徴量計算ステップで算出された特徴量に基づいて、上記流水領域及び上記非流水領域を判別する領域判別ステップとをコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、流水領域検出システム、流水領域検出方法、及びプログラムに関し、詳しくは、撮影して得られる映像信号に基づく河川等の流水領域検出技術に関する。
【背景技術】
【０００２】
河川の氾濫による水害を事前に察知して水防活動に役立てるために、河川敷の各所に水位計や監視カメラを設置し、それらにより得られるデータを一箇所に集約するシステムが、国土交通省が設置するテレメータをはじめとして多数実現されている。また、傾斜模様が描かれた量水板を河川中に設置し、その量水板を含む画像を取り込んで画像処理を行うことで、水位計を設置せずに河川映像から河川の水位を検出する方法が提案されている（例えば、特許文献１参照。）。
【０００３】
【特許文献１】特開平９−１６１０７６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、流水を妨げる物体を河川の中に設置することは、洪水等の災害の要因となり得るため、河川法により厳しく制限管理されており、河川において量水板が設置できる場所は限られている。その一方で、できるだけ多くの箇所から水位データを収集することが、より正確で信頼できる防災情報を生成するためには欠かせない。
【０００５】
本発明は、水中にいかなる物体をも設置することなく、河川等を被写体とする映像信号のみから、その河川等の流水領域を検出できるようにすることを目的とする。
【課題を解決するための手段】
【０００６】
本発明の流水領域検出システムは、固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像が供給され、供給された映像を時間方向に積分する積分手段と、上記積分手段により得られた映像のエッジ成分を強調させる空間フィルタ処理手段と、上記空間フィルタ処理手段で処理された映像における各画素に対してエッジ画素であるか否かを画素毎に判定し、当該判定結果に基づいて上記流水領域及び上記非流水領域を判別する領域判別手段とを備えることを特徴とする。
また、本発明の流水領域検出システムは、固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像が供給され、供給された映像を時間方向に積分する積分手段と、上記積分手段により得られた映像のエッジ成分を強調させる空間フィルタ処理手段と、上記空間フィルタ処理手段で処理された映像における各ライン毎に求めた特徴量に基づいて上記流水領域及び上記非流水領域を判別する領域判別手段とを備えることを特徴とする。
本発明の流水領域検出方法は、固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分工程と、上記積分工程で得られた映像のエッジ成分を強調させる空間フィルタ処理工程と、上記空間フィルタ処理工程で処理された映像における各画素に対してエッジ画素であるか否かを画素毎に判定する判定工程と、上記判定工程での判定結果に基づいて、上記流水領域及び上記非流水領域を判別する領域判別工程とを有することを特徴とする。
また、本発明の流水領域検出方法は、固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分工程と、上記積分工程で得られた映像のエッジ成分を強調させる空間フィルタ処理工程と、上記空間フィルタ処理工程で処理された映像の各ライン毎に特徴量を算出する特徴量計算工程と、上記特徴量計算工程で算出された特徴量に基づいて、上記流水領域及び上記非流水領域を判別する領域判別工程とを有することを特徴とする。
本発明のプログラムは、固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分ステップと、上記積分ステップで得られた映像のエッジ成分を強調させる空間フィルタ処理ステップと、上記空間フィルタ処理ステップで処理された映像における各画素に対してエッジ画素であるか否かを画素毎に判定する判定ステップと、上記判定ステップでの判定結果に基づいて、上記流水領域及び上記非流水領域を判別する領域判別ステップとをコンピュータに実行させることを特徴とする。
また、本発明のプログラムは、固定カメラを用いて撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分する積分ステップと、上記積分ステップで得られた映像のエッジ成分を強調させる空間フィルタ処理ステップと、上記空間フィルタ処理ステップで処理された映像の各ライン毎に特徴量を算出する特徴量計算ステップと、上記特徴量計算ステップで算出された特徴量に基づいて、上記流水領域及び上記非流水領域を判別する領域判別ステップとをコンピュータに実行させることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、河川等を被写体として固定カメラで撮影された流水領域及び非流水領域が含まれる映像を時間方向に積分して流水領域におけるエッジ成分を抑圧し、その映像における非流水領域におけるエッジ成分を強調させる空間フィルタを画像全体に施す。そして、空間フィルタの施された映像の各画素に対してエッジ成分であるか否かを判定し、その結果に基づいて流水領域及び非流水領域を判別する。あるいは、空間フィルタの施された映像における各ライン毎に求めた特徴量に基づいて流水領域及び非流水領域を判別する。これにより、水中にいかなる物体をも設置することなく、河川等を被写体とする映像信号のみから、その流水領域を検出することができる。
【発明を実施するための最良の形態】
【０００８】
以下、本発明の実施形態を図面に基づいて説明する。
【０００９】
＜第１の実施形態＞
図１は、本発明の第１の実施形態による流水領域検出システムの構成例を示すブロック図である。第１の実施形態による流水領域検出システムは、固定カメラＣＭを用いた撮影で得られる河川等を被写体とする映像信号から、時間方向の積分（例えば本実施形態ではフレーム加算）及び空間フィルタ処理に基づき、その河川等における流水領域を検出する。
【００１０】
第１の実施形態による流水領域検出システムは、図１に示すように、領域抽出部１、フレーム加算部２、空間フィルタ処理部３、二値化処理部４、及び領域分割部５Ａを有する。以下、図１に示した流水領域検出システムの各機能部について詳細に説明する。なお、以下においては、必要に応じて図２に示すような実際の河川画像について流水領域及び非流水領域を検出、判別することで、河川の水位検出を行う場合を適宜例示して説明する。
【００１１】
（領域抽出部１）
領域抽出部１は、固定カメラＣＭを用いて撮影された撮像映像から、半分に陸部が、もう半分には流水部が写っているような領域を数フレームにわたり抽出する。ここで、流水部とは、撮像映像における水面等の水が写りこんだ領域であり、陸部とは、流水部以外の領域（例えば、背景画像等の領域）である。
【００１２】
例えば、領域抽出部１は、図２（ａ）に示すように上半分が陸部、下半分が流水部であるような領域画像を撮像映像から抽出する。なお、領域画像を抽出する際、陸部が下で流水部が上、又は陸部及び流水部がそれぞれ左と右、あるいは右と左に写っていても良く、流水領域の検出は可能である。また、撮影映像全体の半分が陸部、もう半分が流水部となるように固定カメラＣＭのアングル等を調整して撮影を行い、得られた撮影映像そのものを領域画像としても良い。
【００１３】
ここで、領域画像における陸部には模様のある静止した物体が写っており，流水部においては水面が常に動いていることが望ましい。また、領域画像において水面は水平であることが望ましく、そうなるように撮影画像（領域画像）を９０度回転させるか拡大や縮小をするなどのアフィン変換を施すことも効果的である。また、夜間に撮影した映像については、例えば赤外線照射した結果として、陸部に何らかの模様が写ってさえいれば、流水部には何も写っていなくても良く、流水領域を検出可能である。
【００１４】
領域画像の抽出は、ある時点でユーザーが抽出位置を指定した後は、流水領域の検出結果を基に得られた水位に応じて水面が常に領域画像の中心付近に位置するようにして、ある時間間隔で自動的に実施するようにしても良い。
【００１５】
（フレーム加算部２）
フレーム加算部２は、領域抽出部１により抽出された領域画像について、複数のフレームを同期加算する。すなわち、フレーム加算部２は、領域抽出部１により抽出された領域画像を時間方向に積分する。これにより、動いている流水部がぼやけた、あるいは滑らかとなった画像が得られる。つまり、流水部のエッジ成分（高周波成分）が抑圧される。その一方で、静止している陸部は特に変化しない。このことは、例えば図２（ａ）及び図２（ｂ）を比較することで確認できる。図２（ｂ）は、フレーム加算部２により複数フレームを同期加算して得られたものである。
【００１６】
ただし、複数のフレームの領域画像を加算する際、単純に加算を行うと画素値も単純に増大してしまい要求されるリソース等も増大するので、例えば加算したフレーム数で画素値を割るなどして画素値の値域拡大を防止することが望ましい。また、例えば、画像を加算する際に、画素値に重み係数を乗じたり、あるいは画素値に関しビットシフトを行ったりしても良い。
【００１７】
例えば、フレーム加算部２は、４枚のフレームを加算する場合には，各フレームにおける画素値に０．２５を乗じるか、あるいはそれと等価になるように画素値を右に２ビットシフトさせた後に、既にフレームメモリに格納されている画素値と加算し、加算結果をフレームメモリに新たに格納する。これにより、加算される枚数分のフレームメモリを用意することなく、少ないフレームメモリでフレーム加算を実現することができる。
【００１８】
また、フレームメモリに格納されている画素値と、領域抽出部１により抽出された領域画像の画素値とを加算して２で割った後、その加算結果をフレームメモリに新たに格納する。この処理を繰り返すことで，複数のフレームの領域画像が重み加算された画像が生成される。あるいは、固定カメラＣＭで撮像する際に、シャッターを比較的長い時間に亘り開けておき、露光時間を比較的長くすることでも、上述したように複数のフレームを加算した場合と同等の効果を得ることができる。
【００１９】
また、複数のフレームにわたって画像を加算する前に、加算する画像に対して、例えば時間方向のメディアンフィルタ処理といった非線形処理を施すようにしても良い。このように、加算する領域画像に非線形処理を施すことで、固定カメラＣＭの前を遮る鳥や虫、降雨雪粒子、あるいは車のライトや水滴といった一時的に発生した外乱による影響を取り除くことができ、安定した水位検出を行うことが可能となる。
【００２０】
ここで、詳細は後述するが本実施形態では、映像においてエッジ画素を多く含む領域を陸部分、エッジ画素が少ない領域を流水部分として、陸部及び流水部の判別を行う。そして、陸部と流水部の境界の位置を水面として認識することで、図２（ａ）に示したような水位を算出することが可能となる。
【００２１】
（空間フィルタ処理部３）
空間フィルタ処理部３は、フレーム加算部２により生成された画像に対して、流水部には少ないが陸部には多く存在するエッジ成分、つまり流水部と陸部の差異を強調させるための空間フィルタ処理を施す。例えば、図２（ｂ）に示した画像に、空間フィルタ処理部３によるフィルタ処理を施すことで、図２（ｃ）に示す画像が得られる。
【００２２】
本実施形態では、空間フィルタ処理部３は、画像の画素値と空間フィルタ係数値との畳み込み処理により、画像のエッジ成分を強調させる。その空間フィルタの係数値としては、微分処理に相当するものが効果的である。
ここで、空間フィルタの伝達関数を、
【００２３】
【数１】
2007256254.tif 000003
【００２４】
と表現すると、空間フィルタ処理部３での処理に用いる空間フィルタの係数値としては、上記式（１）においてＴ１＝Ｔ２＝１として、
【００２５】
【数２】
2007256254.tif 000004
【００２６】
と表されるＰｒｅｗｉｔｔオペレータや、又は、
【００２７】
【数３】
2007256254.tif 000005
【００２８】
と表されるＳｏｂｅｌオペレータが効果的である。これらは、どちらも画像における横方向に微分処理を施し、かつ縦方向に平均化処理を施したことに等しい。
また、空間フィルタの伝達関数を、
【００２９】
【数４】
2007256254.tif 000006
【００３０】
と表した場合の、Ｈ（ｚ）＝１−ｚ-1なるフィルタ及びＨ（ｚ）＝１＋ｚ-1なるフィルタを画像における横方向にそれぞれ１回づつ適用し、かつＨ（ｚ）＝１＋ｚ-1なるフィルタを縦方向に２回適用することで、上述したＳｏｂｅｌオペレータと等価な処理となる。これらのフィルタの組み合わせや適用回数を適宜調整することで、より効果的な処理が可能となることもある。
【００３１】
また、例えば、上記式（１）においてＴ１やＴ２の値を１よりも大きな整数値とし、ｈ（ｋ1，ｋ2）を任意に与えることで、自由度が高くなるため、より効果的な処理が可能となる場合もある。
特に、流水部に対して空間フィルタを適用した結果の分散値が最小となるように、最小自乗法により空間フィルタのフィルタ係数値を計算するようにしても良い。
あるいは、空間フィルタ処理後の分散値について、流水部と陸部の比が最大となるように空間フィルタのフィルタ係数値やタップ数等のパラメータを決定するようにしても良い。例えば、空間フィルタ処理後の流水部の分散値σW2が一定という制約条件の下で、陸部の分散値σL2を最大化するようにする。すなわち、下記式（５）に示す評価関数Ｌを最大化するように空間フィルタに係るパラメータを決定しても良い。
【００３２】
【数５】
2007256254.tif 000007
【００３３】
この問題は、ラグランジュの未定乗数法により固有値問題として解けることが知られている。なお、逆に、陸部の分散値が一定という制約条件の下で流水部の分散値を最小化するように空間フィルタに係るパラメータを決定しても良い。
また、画像によっては上述のような画像における横方向ではなく、縦方向に微分処理を施した方が効果的な場合もある。また、平均化処理を施す必要がない場合もある。
【００３４】
上述したような空間フィルタの方向、係数値、タップ数等のパラメータは、時々刻々と変化する撮影環境（状況）に応じて、リアルタイムで最適なものに適宜更新することが望ましい。それにより、日照の変化や水質・濁度・流速などの変化に伴う撮影画像の変動に対応することができ、精度が高くかつ安定した流水領域検出を行うことができる。
【００３５】
（二値化処理部４）
二値化処理部４は、空間フィルタ部３により生成された画像における個々の画素に対して、エッジ画素であるか否かの判定を画素毎に行い、判定結果に基づいて各画素をエッジ画素か否か（非エッジ画素）に分類する。さらに、二値化処理部４は、その分類に応じて画素値を例えば白（エッジ画素）又は黒（非エッジ画素）の何れかに二値化し、図２（ｄ）に一例を示すようような二値画像を得る。
【００３６】
ここで、二値化処理部４にて二値化処理を行う際のしきい値については，例えば空間フィルタ部３により生成された画像の画素値に応じて計算すれば良く、例えばその平均値やそのヒストグラムの重心値などを適用できる。
【００３７】
（領域分割部５Ａ）
領域分割部５Ａは、二値化処理部４により生成された画像に対して、エッジ画素を多く含む領域を陸部（非流水領域）、それ以外を流水部（流水領域）と判別して領域を二分割する。この判別処理は、二値化処理部４により生成された画像におけるエッジ画素の分布又は非エッジ画素の分布に基づいてなされる。
【００３８】
例えば、このようにして得られた２つの領域の境界を出力することで、図２（ａ）に示したように陸部と流水部の境界として水位が算出される。
ここで、二値化処理部４及び領域分割部５Ａにより、本発明の領域判別手段が構成される。
【００３９】
二値化処理部４により生成された画像において、領域を陸部（非流水領域）及び流水部（流水領域）に二分割する領域分割部５Ａにおける処理については、様々な形態が考えられる。以下、図３〜図６を参照して、領域分割部５Ａにおける処理の一例について説明する。
【００４０】
図３は、領域分割部５Ａの構成例を示すブロック図であり、図４は、図３に示す領域分割部５Ａによる領域分割の原理を説明するための図である。
【００４１】
（横方向ヒストグラム計算部１１）
横方向ヒストグラム計算部１１は、図１に示した二値化処理部４により生成された二値画像について、各行毎に幾つのエッジ画素が存在するかを数える。つまり、横方向ヒストグラム計算部１１は、二値画像における各行毎にエッジ画素数を積算し、エッジ画素の横方向のヒストグラムを生成する。
【００４２】
例えば、図４（ａ）に示す二値画像において、上から一行目には６個のエッジ画素、二行目には７個のエッジ画素が存在している。したがって、横方向ヒストグラム計算部１１により生成される横方向のヒストグラム（図４（ｂ））には、それぞれの行、すなわち縦位置について、エッジ画素数が６及び７としてプロットされる。
【００４３】
（縦方向ヒストグラム計算部１２）
縦方向ヒストグラム計算部１２は、横方向ヒストグラム計算部１１により生成された横方向ヒストグラムを参照し、各エッジ画素数の値について、対応する縦位置が幾つ存在するかを数える。つまり，縦方向ヒストグラム計算部１２は、横方向ヒストグラムにおける各エッジ画素数に関する縦方向のヒストグラムを生成する。
【００４４】
例えば，図４（ｂ）に示した横方向ヒストグラムにおいて、左から一列目には２個のエッジ画像数、二列目には１個のエッジ画像数がある。縦方向ヒストグラム計算部１２により生成される横方向のヒストグラム（図４（ｃ））には、それぞれの列、すなわちエッジ画素数について、エッジ画素数の縦方向積算値が２及び１としてプロットされる。
【００４５】
（境界値決定部１３）
縦方向ヒストグラム計算部１２により生成された縦方向ヒストグラムには、陸部に関連するヒストグラム部分と、流水部に関連するヒストグラム部分の２つのクラスタ（分布）が存在する。したがって、これら２つのクラスタを分離する境界部分を決定することができる。境界値決定部１３は、その境界部分（境界画素数）を決定する。この境界部分の決定方法としては種々の方法を用いることができる。
【００４６】
例えば、境界値決定部１３は、クラス内の分散が小さく、かつクラス間の距離が大きくなるように境界を定めても良いし、又はクラス内分散・クラス間分散比が最大となるように境界を定めても良い。あるいは、境界値決定部１３は、縦方向ヒストグラムを双峰性のガウス分布に帰着させて２つのクラスの境界を決定しても良い。
【００４７】
（水位計算部１４）
水位計算部１４は、境界値決定部１３により決定された縦方向ヒストグラムにおける２つのクラスの境界に基づいて、それに対応する横方向ヒストグラムにおける境界を定め、更にそれに対応する二値画像の境界を定める。これにより、画像が陸部と流水部の領域に二分割され、その境界を水位として出力する。
【００４８】
具体的には、境界値決定部１３により決定された境界をしきい値として、それに対応する図４（ｃ）に示した縦方向ヒストグラムにおける境界Ｐ１が、図４（ｂ）に示した横方向ヒストグラムにおける対応点Ｐ２を経て、二値画像における陸部と流水部との境界である水位Ｐ３が求められる。実際の河川画像に適用した場合を図５に示している。
【００４９】
なお、上述した説明では、画像において水面が水平である（陸部と流水部の境界が横方向に伸びている）場合を一例として示したが、画像において水面が垂直である（陸部と流水部の境界が縦方向に伸びている）場合には、各列毎にエッジ画素数を積算してエッジ画素の縦方向のヒストグラムを生成し、それを参照して各エッジ画素数に関する横方向のヒストグラムを生成するようにすれば良い。
【００５０】
以上、説明したように第１の実施形態によれば、河川等を被写体として固定カメラＣＭを用いて撮影された撮像映像から、領域抽出部１により陸部と流水部の境界、すなわち陸部と流水部を含む領域画像を抽出し、フレーム加算部２により各領域について数フレーム分を加算して流水部が滑らかになった画像を生成する。さらに、その画像におけるエッジ成分を強調させる空間フィルタを空間フィルタ処理部３にて画像全体に施し、得られた画像の各画素に対してエッジ成分か否かに分類する二値化処理を二値化処理部４で施し、領域分割部５Ａによりエッジ成分を多く含む領域を陸部、それ以外を流水部と判別して領域分割を行う。これにより、陸部及び流水部をそれぞれ検出でき、その二領域の境界を水位と判定することで、水中にいかなる物体をも設置することなく、河川等を被写体とする映像信号のみから、その流水領域を検出することができる。
【００５１】
＜第２の実施形態＞
次に、本発明の第２の実施形態について説明する。
第２の実施形態による流水領域検出システムは、固定カメラを用いた撮影で得られる河川等を被写体とする映像信号から、時間方向の積分（フレーム加算等）及び空間フィルタ処理に基づき、その河川等における流水領域を検出する。
【００５２】
図８は、第２の実施形態による流水領域検出システムの構成例を示すブロック図である。なお、この図８において、図１に示したブロック等と同一の機能を有するブロック等には同一の符号を付し、重複する説明は省略する。
図８に示すように、第２の実施形態による流水領域検出システムは、領域抽出部１、フレーム加算部２、空間フィルタ処理部３、ライン特徴量計算部６、及び領域分割部５Ｂを有する。
【００５３】
（空間フィルタ処理部３）
空間フィルタ処理部３は、パラメータ計算部７を有し、フレーム加算部２により生成された画像に対して、流水部には少ないが陸部には多く存在するエッジ成分、つまり流水部と陸部の差異を強調させるための空間フィルタ処理を施す。例えば、領域抽出部１により撮像映像から抽出した図９（ａ）に示す領域画像を、フレーム加算部２により複数フレームを同期加算して得られる図９（ｂ）に示す画像に、空間フィルタ処理部３によるフィルタ処理を施すことで、図９（ｃ）に示す画像が得られる。
【００５４】
本実施形態では、空間フィルタ処理部３は、画像の画素値と空間フィルタ係数値との畳み込み処理により、画像のエッジ成分を強調させる。空間フィルタの係数値としては、微分処理に相当するものが効果的である。
【００５５】
この空間フィルタの係数値は、空間フィルタ処理部３内のパラメータ計算部７により算出される。パラメータ計算部７は、空間フィルタの方向、係数値、タップ数等のパラメータを算出し、時々刻々と変化する撮影環境（状況）に応じて、最適なものに適宜更新する処理をリアルタイムで行う。なお、パラメータ計算部７による空間フィルタの係数値やタップ数等のパラメータの決定方法については、第１の実施形態での空間フィルタに係るパラメータの決定方法と同様であるので説明は省略する。
【００５６】
このように、パラメータ計算部７が、時々刻々と変化する撮影環境に応じて、リアルタイムで最適な空間フィルタに係るパラメータを算出して適宜更新することにより、日照の変化や水質・濁度・流速などの変化に伴う撮影画像の変動に対応することができ、精度が高くかつ安定した流水領域検出を行うことができる。また、例えば固定カメラＣＭの設置位置を変えても、空間フィルタに係るパラメータを自動的に最適化することができる。
【００５７】
（ライン特徴量計算部６）
ここで、図９（ａ）に示した画像をフレーム加算部２により複数フレームにわたって同期加算することにより、図９（ｂ）に示したように、動いている流水部については高周波成分が低減された滑らかな信号となり、静止している陸部については高周波成分が維持された明確な輪郭を持つ信号となる。こうした違いを顕在化させるために、図９（ｂ）に示した画像に対して空間フィルタ処理部３により空間フィルタ処理を施すことで、流水部に対して陸部の輪郭（高周波）成分が強調された図９（ｃ）に示した画像が得られる。
【００５８】
ライン特徴量計算部６は、このように空間フィルタ処理部３で処理された画像の各ラインの特徴量として、画像の各ライン毎に分散値を算出する。図９（ｃ）に示した画像に基づいて、ライン特徴量計算部６により算出した特徴量を図９（ｄ）に示す。図９（ｄ）に示されるように、算出された特徴量の値は、概して陸部では大きく、流水部では小さくなっている。
【００５９】
なお、ライン特徴量計算部６は、画像の各ラインの特徴量として分散値を算出するようにしているが、画像の各ライン毎に絶対値和を算出するようにしても良い。絶対値和を算出するようにした場合には、分散値を算出する場合に比べて、乗算処理が含まれないために演算量を低減させることができる。
【００６０】
（領域分割部５Ｂ）
領域分割部５Ｂは、ライン特徴量計算部６により算出された各ラインの特徴量に基づいて、画像における陸部（非流水領域）と流水部（流水領域）とを判別して領域を二分割する。
【００６１】
ライン特徴量計算部６により算出された特徴量に基づいて、領域を陸部（非流水領域）及び流水部（流水領域）に二分割する領域分割部５Ｂにおける処理については、様々な形態が考えられる。例えば、図９（ｄ）に示される特徴量を鍵型図形Ｐ１１に近似し、この鍵のくびれ部分を陸部と流水部の境界として水位とすることもできる。また、例えば、クラス内分散・クラス間分散比最大基準により２クラスに分離しするといった線形判別法を適用することもできる。このようにして得られる２つの領域の境界を水位として出力する。
ここで、ライン特徴量計算部６及び領域分割部５Ｂにより、本発明の領域判別手段が構成される。
【００６２】
以上、説明したように第２の実施形態によれば、河川等を被写体として固定カメラＣＭを用いて撮影された撮像映像から、領域抽出部１により陸部と流水部の境界、すなわち陸部と流水部を含む領域画像を抽出し、フレーム加算部２により各領域について数フレーム分を加算して流水部が滑らかになった画像を生成する。さらに、パラメータ計算部６により撮影状況等に応じて空間フィルタのパラメータを適宜変更しつつ、その画像におけるエッジ成分を強調させる空間フィルタを空間フィルタ処理部３にて画像全体に施し、得られた画像の各ライン毎に特徴量をライン特徴量計算部６で算出し、領域分割部５Ｂにより各ラインの特徴量に基づいて画像を流水部と陸部の２クラスに分類して領域分割を行う。これにより、陸部及び流水部をそれぞれ検出でき、その二領域の境界を水位と判定することで、水中にいかなる物体をも設置することなく、河川等を被写体とする映像信号のみから、その流水領域を検出することができる。
【００６３】
また、第２の実施形態によれば、ラインの特徴量に基づいて領域分割を行って水位を検出することで、第１の実施形態で行っていた二値化処理や縦方向及び横方向のヒストグラム計算が不要となるとともに、画素単位の処理ではなくライン単位での処理となる。これにより、システム全体での演算処理量が低減され、高精度な水位計算処理システム（処理回路等）の省電力化や小型化が可能となる。また、二値化処理においては、設定する閾値の調整が容易ではなく、その決定方法に依存して水位計算の精度が異なり、ヒストグラム計算においては、二値化処理の結果を累積するために二値化誤差をも累積してしまうが、第２の実施形態では、このような処理を用いないことで、より高精度な水位検出を行うことができる。
【００６４】
なお、上述した第１及び第２の実施形態では、領域抽出部１は、固定カメラＣＭを用いて撮影された撮像映像から１つの領域画像を抽出するようにしているが、図６に示すように、１枚のフレームに対して複数の領域画像（図６に示す例では、領域１、領域２、及び領域３）を抽出するようにしても良い。その場合、複数の領域画像において、その一部領域が重複していても構わない。これらの各領域画像のそれぞれについて、上述した処理により陸部及び流水部を判別し、得られた複数の結果から総合的に判断して最終的な陸部及び流水部の領域分割を行うようにしても良い。例えば、複数の結果の平均値を最終的な陸部及び流水部の境界としたり、多数決法により最終的な陸部及び流水部の領域を決定することで、固定カメラＣＭのレンズに雨水が付着したり、降雨や降雪あるいは霧による画像障害が発生して流水部が検出しにくい場合でも、流水部を高精度かつ安定して検出することが可能となる。
【００６５】
（本発明の他の実施形態）
なお、上述した実施形態は、コンピュータのＣＰＵ又はＭＰＵ、ＲＡＭ、ＲＯＭ等で構成可能なものであり、ＲＡＭやＲＯＭに記憶されたプログラムが動作することによって実現でき、上記プログラムは本発明の実施形態に含まれる。また、コンピュータが上述の実施形態の機能を果たすように動作させるプログラムを、例えばＣＤ−ＲＯＭのような記録媒体に記録し、コンピュータに読み込ませることによって実現できるものであり、上記プログラムを記録した記録媒体は本発明の実施形態に含まれる。上記プログラムを記録する記録媒体としては、ＣＤ−ＲＯＭ以外に、フレキシブルディスク、ハードディスク、磁気テープ、光磁気ディスク、不揮発性メモリカード等を用いることができる。
また、コンピュータがプログラムを実行し処理を行うことにより、上述の実施形態の機能が実現されるプログラムプロダクトは、本発明の実施形態に含まれる。上記プログラムプロダクトとしては、上述の実施形態の機能を実現するプログラム自体、上記プログラムが読み込まれたコンピュータ、ネットワークを介して通信可能に接続されたコンピュータに上記プログラムを提供可能な送信装置、及び当該送信装置を備えるネットワークシステム等がある。
また、コンピュータが供給されたプログラムを実行することにより上述の実施形態の機能が実現されるだけでなく、そのプログラムがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）又は他のアプリケーションソフト等と共同して上述の実施形態の機能が実現される場合や、供給されたプログラムの処理のすべて又は一部がコンピュータの機能拡張ボードや機能拡張ユニットにより行われて上述の実施形態の機能が実現される場合も、かかるプログラムは本発明の実施形態に含まれる。
【００６６】
例えば、図７に示すようなコンピュータ機能７００を有し、そのＣＰＵ７０１により上述した実施形態での動作が実施される。
コンピュータ機能７００は、図７に示すように、ＣＰＵ７０１と、ＲＯＭ７０２と、ＲＡＭ７０３と、キーボード（ＫＢ）７０９のキーボードコントローラ（ＫＢＣ）７０５と、表示部としてのＣＲＴディスプレイ（ＣＲＴ）７１０のＣＲＴコントローラ（ＣＲＴＣ）７０６と、ハードディスク（ＨＤ）７１１及びフレキシブルディスク（ＦＤ）７１２のディスクコントローラ（ＤＫＣ）７０７と、ネットワークインタフェースカード（ＮＩＣ）７０８とが、システムバス７０４を介して互いに通信可能に接続された構成としている。
ＣＰＵ７０１は、ＲＯＭ７０２又はＨＤ７１１に記憶されたソフトウェア、又はＦＤ７１２より供給されるソフトウェアを実行することで、システムバス７０４に接続された各構成部を総括的に制御する。
すなわち、ＣＰＵ７０１は、上述したような動作を行うための処理プログラムを、ＲＯＭ７０２、ＨＤ７１１、又はＦＤ７１２から読み出して実行することで、上述した実施形態での動作を実現するための制御を行う。
ＲＡＭ７０３は、ＣＰＵ７０１の主メモリ又はワークエリア等として機能する。
ＫＢＣ７０５は、ＫＢ７０９や図示していないポインティングデバイス等からの指示入力を制御する。
ＣＲＴＣ７０６は、ＣＲＴ７１０の表示を制御する。
ＤＫＣ７０７は、ブートプログラム、種々のアプリケーション、ユーザファイル、ネットワーク管理プログラム、及び上記処理プログラム等を記憶するＨＤ７１１及びＦＤ７１２とのアクセスを制御する。
ＮＩＣ７０８はネットワーク７１３上の他の装置と双方向にデータをやりとりする。
【００６７】
なお、上記実施形態は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、またはその主要な特徴から逸脱することなく、様々な形で実施することができる。
【図面の簡単な説明】
【００６８】
【図１】本発明の第１の実施形態による流水領域検出システムの構成例を示す図である。
【図２】図１に示す流水領域検出システムで処理された画像の一例を示す図である。
【図３】領域分割部の構成例を示す図である。
【図４】図３に示す領域分割部による領域分割の原理を説明するための図である。
【図５】図３に示す領域分割部による領域分割を実際の河川画像に適用した場合を示す図である。
【図６】領域抽出部による領域画像抽出の他の例を示す図である。
【図７】本実施形態における流水領域検出システムを実現可能なコンピュータ機能を示すブロック図である。
【図８】本発明の第２の実施形態による流水領域検出システムの構成例を示す図である。
【図９】第２の実施形態による流水領域検出システムでの処理を説明するための図である。
【符号の説明】
【００６９】
１  領域抽出部
２  フレーム加算部
３  空間フィルタ処理部
４  二値化処理部
５Ａ、５Ｂ  領域分割部
６  ライン特徴量計算部
７  パラメータ計算部
１１  横方向ヒストグラム計算部
１２  縦方向ヒストグラム計算部
１３  境界値決定部
１４  水位計算部
【図１】
2007256254.tif 000008
【図２】
2007256254.tif 000009
【図３】
2007256254.tif 000010
【図４】
2007256254.tif 000011
【図５】
2007256254.tif 000012
【図６】
2007256254.tif 000013
【図７】
2007256254.tif 000014
【図８】
2007256254.tif 000015
【図９】
2007256254.tif 000016