(57)【要約】【課題】車両に搭載されたカメラの映像から道路の路面上の水浮を検知し得る、水浮検知装置、水浮検知方法、及びプログラムを提供する。【解決手段】水浮検知装置は、道路の画像から、処理対象となる領域におけるヒストグラムを生成する、ヒストグラム生成部１１と、ヒストグラムを、各部分の面積が同一となるように、設定数で分割し、各部分の中心の輝度をそれぞれの初期値として算出する、初期値算出部１２と、部分毎に、初期値を用いて正規分布を推定する、正規分布推定部１３と、推定された正規分布のうち、設定条件に合致する正規分布を選択し、選択した正規分布を用いて新たな正規分布を生成し、新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する、特徴量算出部１４と、特徴量に基づいて、道路の路面上に前記水浮が発生しているかどうかを判定する、路面状況推定部１５とを備えている。【選択図】図１
2015088106.tif 000002
【特許請求の範囲】
【請求項１】
道路の画像から路面上の水浮を検知するための装置であって、
前記画像から、処理対象となる領域におけるヒストグラムを生成する、ヒストグラム生成部と、
生成された前記ヒストグラムを、各部分の面積が同一となるように、設定数で分割し、更に、各部分の中心の輝度をそれぞれの初期値として算出する、初期値算出部と、
前記部分毎に、当該部分について算出された前記初期値を用いて、正規分布を推定する、正規分布推定部と、
推定された前記正規分布のうち、設定条件に合致する正規分布を選択し、更に、選択した前記正規分布を用いて新たな正規分布を生成し、生成した前記新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する、特徴量算出部と、
算出された前記特徴量に基づいて、前記道路の路面上に前記水浮が発生しているかどうかを判定する、路面状況推定部と、
を備えていることを特徴とする水浮検知装置。
【請求項２】
前記ヒストグラム生成部が、生成した前記ヒストグラムに対して平滑化処理を実行し、
前記初期値算出部が、前記平滑化処理が実行された前記ヒストグラムを分割する、
請求項１に記載の水浮検知装置。
【請求項３】
前記画像が、設定されたフレームレートでカメラから出力されてきた画像であり、
前記ヒストグラム生成部、前記初期値算出部、前記正規分布推定部、及び前記特量算出部は、前記カメラから前記画像が出力される度に、処理を実行し、
前記路面状況推定部は、前記特徴量算出部によって前記特徴量が算出される度に、そのとき算出された特徴量と、過去に算出された特徴量とを平均し、更に、平均によって得られた値を設定数だけ累積し、累積によって得られた累積値に基づいて、前記道路に前記水浮が発生しているかどうかを判定する、
請求項１または２に記載の水浮検知装置。
【請求項４】
前記特徴量算出部が、設定された昼間の時間帯においては、前記新たな正規分布の分散のみを、前記特徴量として算出し、設定された夜間の時間においては、前記新たな正規分布の平均値及び分散値を、前記特徴量として算出する、
請求項１〜３のいずれかに記載の水浮検知装置。
【請求項５】
道路の画像から路面上の水浮を検知するための方法であって、
（ａ）前記画像から、処理対象となる領域におけるヒストグラムを生成する、ステップと、
（ｂ）前記（ａ）のステップで生成された前記ヒストグラムを、各部分の面積が同一となるように、設定数で分割し、更に、各部分の中心の輝度をそれぞれの初期値として算出する、ステップと、
（ｃ）前記部分毎に、当該部分について算出された前記初期値を用いて、正規分布を推定する、ステップと、
（ｄ）前記（ｃ）のステップで推定された前記正規分布のうち、設定条件に合致する正規分布を選択し、更に、選択した前記正規分布を用いて新たな正規分布を生成し、生成した前記新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する、ステップと、
（ｅ）前記（ｄ）のステップで算出された前記特徴量に基づいて、前記道路の路面上に前記水浮が発生しているかどうかを判定する、ステップと、
を有することを特徴とする水浮検知方法。
【請求項６】
前記（ａ）のステップにおいて、生成した前記ヒストグラムに対して、更に、平滑化処理を実行し、
前記（ｂ）のステップにおいて、前記平滑化処理が実行された前記ヒストグラムを分割する、
請求項５に記載の水浮検知方法。
【請求項７】
前記画像が、設定されたフレームレートでカメラから出力されてきた画像であり、
前記（ａ）のステップ、前記（ｂ）のステップ、前記（ｃ）のステップ、及び前記（ｄ）のステップは、前記カメラから前記画像が出力される度に、実行され、
前記（ｅ）のステップにおいて、前記（ｄ）のステップによって前記特徴量が算出される度に、そのとき算出された特徴量と、過去に算出された特徴量とを平均し、更に、平均によって得られた値を設定数だけ累積し、累積によって得られた累積値に基づいて、前記道路に前記水浮が発生しているかどうかを判定する、
請求項５または６に記載の水浮検知方法。
【請求項８】
前記（ｄ）のステップにおいて、設定された昼間の時間帯においては、前記新たな正規分布の分散のみを、前記特徴量として算出し、設定された夜間の時間においては、前記新たな正規分布の平均値及び分散値を、前記特徴量として算出する、
請求項５〜７のいずれかに記載の水浮検知方法。
【請求項９】
コンピュータによって、道路の画像から路面上の水浮を検知するためのプログラムであって、
前記コンピュータに、
（ａ）前記画像から、処理対象となる領域におけるヒストグラムを生成する、ステップと、
（ｂ）前記（ａ）のステップで生成された前記ヒストグラムを、各部分の面積が同一となるように、設定数で分割し、更に、各部分の中心の輝度をそれぞれの初期値として算出する、ステップと、
（ｃ）前記部分毎に、当該部分について算出された前記初期値を用いて、正規分布を推定する、ステップと、
（ｄ）前記（ｃ）のステップで推定された前記正規分布のうち、設定条件に合致する正規分布を選択し、更に、選択した前記正規分布を用いて新たな正規分布を生成し、生成した前記新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する、ステップと、
（ｅ）前記（ｄ）のステップで算出された前記特徴量に基づいて、前記道路の路面上に前記水浮が発生しているかどうかを判定する、ステップと、
を実行させるプログラム。
【請求項１０】
前記（ａ）のステップにおいて、生成した前記ヒストグラムに対して、更に、平滑化処理を実行し、
前記（ｂ）のステップにおいて、前記平滑化処理が実行された前記ヒストグラムを分割する、
請求項９に記載のプログラム。
【請求項１１】
前記画像が、設定されたフレームレートでカメラから出力されてきた画像であり、
前記コンピュータが、前記（ａ）のステップ、前記（ｂ）のステップ、前記（ｃ）のステップ、及び前記（ｄ）のステップを、前記カメラから前記画像が出力される度に、実行し、
前記（ｅ）のステップにおいて、前記（ｄ）のステップによって前記特徴量が算出される度に、そのとき算出された特徴量と、過去に算出された特徴量とを平均し、更に、平均によって得られた値を設定数だけ累積し、累積によって得られた累積値に基づいて、前記道路に前記水浮が発生しているかどうかを判定する、
請求項９または１０に記載のプログラム。
【請求項１２】
前記（ｄ）のステップにおいて、設定された昼間の時間帯においては、前記新たな正規分布の分散のみを、前記特徴量として算出し、設定された夜間の時間においては、前記新たな正規分布の平均値及び分散値を、前記特徴量として算出する、
請求項９〜１１のいずれかに記載のプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、道路の路面上の水浮を検出するための、検知装置及び検知方法に関し、更には、これらを実現するためのプログラムに関する。
【背景技術】
【０００２】
一般に、自動車にとって、雨天等によって道路の路面上に水浮が生じている状態は、事故発生率が高い状態といえるので、ドライバーは慎重に運転を行う必要がある。但し、時間帯、天候によっては、道路の路面上に水浮が生じているか否かを判断することが、ドライバーにとって難しい場合がある。このため、路面上の水浮の検知を行い、水浮が生じている場合に、そのことをドライバーに通知するシステムが、車両に搭載されれば、交通事故の発生が抑制されると考えられる。
【０００３】
また、近年、自動車においては、トラクションコントロール装置、アンチロックブレーキシステム、横滑り防止装置、といった各種電子デバイスが搭載されている。これらの電子デバイスは、路面上に水浮が存在する場合において、車両の安定性の確保に有効となるが、このとき、車両において路面上の水浮の検知が行われていれば、タイヤに発生するスリップを予測できるので、これらの電子デバイスによる効果をより高めることができる。
【０００４】
具体的には、トラクションコントロール装置は、スリップの予測結果に基づき、スリップが発生する前に、原動機の出力を制御することができる。また、同様に、アンチロックブレーキシステムは、スリップの予測結果に基づき、確実にタイヤのロックを回避できる。更に、同様に、横滑り防止装置は、スリップの予測結果に基づき、横滑りの発生前に各輪のブレーキを制御することができる。これらの結果、車両の安定性の確保がより確実なものとなる。
【０００５】
また、特許文献１は、道路の路面上の水浮を検知する技術の一例を開示している。具体的には、特許文献に開示された技術では、まず、道路に設置された監視カメラの路面映像から得られた特徴量（路面の色相の平均・分散、路面の彩度の平均・分散、路面の明度の平均・分散）の混合正規分布を用いて、種々の路面状況が定義される。次に、現在の路面映像から特徴量を求め、求めた特徴量を、路面状況を定義する混合正規分布に照合することによって、現在の路面状況が判定される。
【０００６】
従って、特許文献１に開示された技術を、車両に適用すれば、ドライバーに対して水浮の発生を通知したり、スリップの予測による車両安定性の更なる向上を図ったり、といったことが可能になると考えられる。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特許第４２２５８０７号
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１に開示された技術では、カメラからの映像に基づいて、路面状況を定義する必要があるため、この技術を、カメラからの映像が刻一刻と変化する車両に適用することは困難である。なお、特許文献１には、カメラの設置条件の変化に応じて路面状況の定義を更新する処理も開示されているが、この処理は、同じ場所での季節の変化などに対応するための処理であり、この処理を用いたとしても、特許文献１に開示された技術を車両に適用することは困難である。
【０００９】
［発明の目的］
本発明の目的の一例は、上記問題を解消し、車両に搭載されたカメラの映像から道路の路面上の水浮を検知し得る、水浮検知装置、水浮検知方法、及びプログラムを提供することにある。
【課題を解決するための手段】
【００１０】
上記目的を達成するため、本発明の一側面における水浮検知装置は、道路の画像から路面上の水浮を検知するための装置であって、
前記画像から、処理対象となる領域におけるヒストグラムを生成する、ヒストグラム生成部と、
生成された前記ヒストグラムを、各部分の面積が同一となるように、設定数で分割し、更に、各部分の中心の輝度をそれぞれの初期値として算出する、初期値算出部と、
前記部分毎に、当該部分について算出された前記初期値を用いて、正規分布を推定する、正規分布推定部と、
推定された前記正規分布のうち、設定条件に合致する正規分布を選択し、更に、選択した前記正規分布を用いて新たな正規分布を生成し、生成した前記新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する、特徴量算出部と、
算出された前記特徴量に基づいて、前記道路の路面上に前記水浮が発生しているかどうかを判定する、路面状況推定部と、
を備えていることを特徴とする。
【００１１】
また、上記目的を達成するため、本発明の一側面における水浮検知方法は、道路の画像から路面上の水浮を検知するための方法であって、
（ａ）前記画像から、処理対象となる領域におけるヒストグラムを生成する、ステップと、
（ｂ）前記（ａ）のステップで生成された前記ヒストグラムを、各部分の面積が同一となるように、設定数で分割し、更に、各部分の中心の輝度をそれぞれの初期値として算出する、ステップと、
（ｃ）前記部分毎に、当該部分について算出された前記初期値を用いて、正規分布を推定する、ステップと、
（ｄ）前記（ｃ）のステップで推定された前記正規分布のうち、設定条件に合致する正規分布を選択し、更に、選択した前記正規分布を用いて新たな正規分布を生成し、生成した前記新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する、ステップと、
（ｅ）前記（ｄ）のステップで算出された前記特徴量に基づいて、前記道路の路面上に前記水浮が発生しているかどうかを判定する、ステップと、
を有することを特徴とする。
【００１２】
更に、上記目的を達成するため、本発明の一側面におけるプログラムは、コンピュータによって、道路の画像から路面上の水浮を検知するためのプログラムであって、
前記コンピュータに、
（ａ）前記画像から、処理対象となる領域におけるヒストグラムを生成する、ステップと、
（ｂ）前記（ａ）のステップで生成された前記ヒストグラムを、各部分の面積が同一となるように、設定数で分割し、更に、各部分の中心の輝度をそれぞれの初期値として算出する、ステップと、
（ｃ）前記部分毎に、当該部分について算出された前記初期値を用いて、正規分布を推定する、ステップと、
（ｄ）前記（ｃ）のステップで推定された前記正規分布のうち、設定条件に合致する正規分布を選択し、更に、選択した前記正規分布を用いて新たな正規分布を生成し、生成した前記新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する、ステップと、
（ｅ）前記（ｄ）のステップで算出された前記特徴量に基づいて、前記道路の路面上に前記水浮が発生しているかどうかを判定する、ステップと、
を実行させることを特徴とする。
【発明の効果】
【００１３】
以上のように、本発明によれば、車両に搭載されたカメラの映像から道路の路面上の水浮を検知することができる。
【図面の簡単な説明】
【００１４】
【図１】図１は、本発明の実施の形態における水浮検知装置の構成を示すブロック図である。
【図２】図２は、本発明の実施の形態における水浮検知装置の動作を示すフロー図である。
【図３】図３は、車両に取り付けられたカメラから出力されてくる画像の一例を示す図である。
【図４】図４は、本発明の実施の形態で生成されるヒストグラムの一例を示す図である。
【図５】図５は、本発明の実施の形態におけるヒストグラムの分割と初期値の算出とを説明する図である。
【図６】図６は、本発明の実施の形態において推定された正規分布の一例を示す図である。
【図７】図７は、本発明の実施の形態において、選択された正規分布から生成された新たな正規分布の一例を示す図である。
【発明を実施するための形態】
【００１５】
（実施の形態）
以下、本発明の実施の形態における、水浮検知装置、水浮検知方法、及びプログラムについて、図１〜図７を参照しながら説明する。
【００１６】
［装置構成］
最初に、本実施の形態における水浮検知装置の構成について説明する。図１は、本発明の実施の形態における水浮検知装置の構成を示すブロック図である。
【００１７】
図１に示す本実施の形態における水浮検知装置１０は、道路の画像から路面上の水浮を検知する装置である。図１に示すように、水浮検知装置１０は、ヒストグラム生成部１１と、初期値算出部１２と、正規分布推定部１３と、特徴量算出部１４と、路面状況推定部１５とを備えている。
【００１８】
ヒストグラム生成部１１は、まず、道路の画像を取得し、取得した道路の画像から、処理対象となる領域におけるヒストグラムを生成する。初期値算出部１２は、ヒストグラム生成部１１によって生成されたヒストグラムを設定数で分割する。このとき、初期値算出部１２による分割は、分割によって得られる各部分の面積が同一となるように行われる。そして、初期値算出部１２は、各部分の中心の輝度をそれぞれの初期値として算出する。
【００１９】
正規分布推定部１３は、初期値算出部１２による分割で得られた部分毎に、各部分について算出された初期値を用いて、正規分布を推定する。特徴量算出部１４は、まず、正規分布推定部１３によって推定された正規分布のうち、設定条件に合致する正規分布を選択し、更に、選択した正規分布を用いて新たな正規分布を生成する。
【００２０】
また、特徴量算出部１４は、生成した新たな正規分布の平均値及び分散のいずれか又は両方を、特徴量として算出する。路面状況推定部１５は、特徴量算出部１４によって算出された特徴量に基づいて、道路の路面に水浮が発生しているかどうかを判定する。
【００２１】
このように、本実施の形態では、画像から生成したヒストグラムから正規分布を生成することによって、処理対象となる領域の特徴量が算出され、更に、この特徴量に基づいて道路の路面上に水浮が発生しているかどうかが判定される。つまり、本実施の形態では、特許文献１に開示された技術と異なり、路面状況を定義しなくても、水浮を検知できる。そして、この一連の処理は、設定されたフレームレートで画像が出力されてくる場合に、フレーム毎に実行可能な処理である。このため、本実施の形態によれば、車両に搭載されたカメラの映像から道路の路面上の水浮を検知することができる。
【００２２】
ここで、本実施の形態における水浮検知装置１０の構成について更に具体的に説明する。図１に示すように、本実施の形態では、水浮検知装置１０は、車両に搭載されるカメラ２０に接続されている。また、カメラ２０は、設定されたフレームレートで画像を出力する機能を備えている。例えば、フレームレートが３０ｆｐｓに設定されている場合であれば、カメラ２０は、１秒間に３０枚の画像を水浮検知装置１０に出力することになる。
【００２３】
また、本実施の形態では、カメラ２０が、設定されたフレームレートで画像を出力するため、ヒストグラム生成部１１、初期値算出部１２、正規分布推定部１３、及び特量算出部１４は、カメラ２０から画像が出力される度に、処理を実行する。
【００２４】
ヒストグラム生成１１は、上述したように、画像から、処理対象となる領域におけるヒストグラムを生成するが、このときの処理対象となる領域は、１つであっても良いし、複数であっても良い。処理対象となる領域の数及び大きさは、道路の路面を表している領域が処理対象となるように、カメラ２０の搭載位置、画角、解像度等に応じて、予め適宜設定される。
【００２５】
また、本実施の形態では、ヒストグラム生成部１１は、生成したヒストグラムに対して平滑化処理を実行でき、更に、初期値算出部１２は、平滑化処理が実行されたヒストグラムに対して分割処理を実行し、初期値を算出することができる。この場合、ヒストグラムにおけるノイズが除去され、初期値の値が適切なものとなる。
【００２６】
特徴量算出部１４は、本実施の形態では、設定された昼間の時間帯においては、新たな正規分布の分散のみを、特徴量として算出することができる。また、特徴量算出部１４は、設定された夜間の時間においては、新たな正規分布の平均値及び分散値を、特徴量として算出することができる。これは、昼間の時間帯に得られた画像は明るいため、分散のみでも水浮を検知することが可能であるためである。また、このような態様とした場合は、昼間の時間帯において、処理速度の向上が図られることになる。なお、各時間帯は、季節の変化に伴って適宜設定されれば良い。
【００２７】
更に、本実施の形態では、路面状況推定部１５は、特徴量算出部１４によって特徴量が算出される度に、そのとき算出された特徴量と、過去に算出された特徴量とを平均し、平均によって得られた値を設定数だけ累積することができる。そして、路面状況推定部１５は、累積によって得られた累積値に基づいて、道路の路面上に水浮が発生しているかどうかを判定することができる。
【００２８】
［装置動作］
次に、本実施の形態における水浮検知装置１０の動作について図２を用いて説明する。図２は、本発明の実施の形態における水浮検知装置の動作を示すフロー図である。以下の説明においては、適宜図１を参酌する。また、本実施の形態では、水浮検知装置１０を動作させることによって、水浮検知方法が実施される。よって、本実施の形態における水浮検知方法の説明は、以下の水浮検知装置１０の動作説明に代える。
【００２９】
図２に示すように、まず、ヒストグラム生成部１１は、カメラ２０から出力されてきた道路の画像を取得し、取得した道路の画像から、処理対象となる領域におけるヒストグラムを生成する（ステップＡ１）。
【００３０】
具体的には、例えば、カメラ２０から図３に示す画像が出力されているとする。この場合、ヒストグラム生成部１１は、図３に示す画像から、処理対象となる破線で囲まれた２つの領域を抽出し、抽出した２つの領域を用いて、図４に示すヒストグラムを生成する。図３は、車両に取り付けられたカメラから出力されてくる画像の一例を示す図である。図４は、本発明の実施の形態で生成されるヒストグラムの一例を示す図である。また、図３において斜線は影を表している。図４において、ヒストグラムの縦軸は画素数（頻度）、横軸は輝度である。
【００３１】
次に、ヒストグラム生成部１１は、ステップＡ１で生成されたヒストグラムに対して平滑化処理を実行する（ステップＡ２）。具体的には、ヒストグラム生成部１１は、ステップＡ１で生成されたヒストグラムにおいて、局所的に変化している部分を除去し、滑らかな曲線とする。なお、ヒストグラムの平滑化は、既存の手法を用いることによって行うことができる。また、図４において、破線は平滑化されたヒストグラムを表している。
【００３２】
次に、初期値算出部１２は、ステップＡ２による平滑化が行われたヒストグラムを設定数で分割し、分割によって得られた部分（以下「分割部分」と表記する。）の中心の輝度を、それぞれの初期値として算出する（ステップＡ３）。
【００３３】
具体的には、図５に示すように、本実施の形態では、初期値算出部１２は、ヒストグラムを４つに分割する。なお、４つに分割しているのは、路面の日陰の部分と、路面の日向の部分と、路面のペイントの日陰の部分と、路面のペイントの日向の部分とに対応する正規分布を推定するためである。図５は、本発明の実施の形態におけるヒストグラムの分割と初期値の算出とを説明する図である。
【００３４】
次に、正規分布推定部１３は、初期値算出部１２によって得られた分割部分毎に、各分割部分について算出された初期値を用いて、正規分布を推定する（ステップＡ４）。
【００３５】
具体的には、正規分布推定部１３は、ステップＡ４においては、既存のＥＭアルゴリズムを利用して、図６に示すように正規部分布を推定する。図６は、本発明の実施の形態において推定された正規分布の一例を示す図である。図６の例では、分割部分毎に、４つの正規分布が推定されている。また、本実施の形態では、推定された４つの正規分布をひとまとめにする場合は、「混合正規分布」と表記する。
【００３６】
但し、ＥＭアルゴリズムを利用する場合においては、ヒストグラムの分布形状によっては、ある箇所の分布について、２つの正規分布の間で競合が生じることがある。このため、正規分布推定部１３は、下記の条件（ａ１）及び（ａ２）のいずれかに合致する場合は、該当する２つの正規分布を統合する。なお、以下においてｎ、ｍは隣接している正規分布のいずれかに対応していることを示している。また、以下に示す、閾値ａ、閾値ｂ、閾値ｃは、事前の実験に基づいて設定される。
【００３７】
条件（ａ１）：隣接する正規分布それぞれの平均値μの差が閾値ａ以下である。
｜｜μｎ−μｍ｜｜＜閾値ａ
条件（ａ２）：隣接する正規分布それぞれの、確率分布ｐｎ（μｎ）及びｐｍ（μｍ）が閾値ｂ未満であり、且つ、それぞれの分散σ２ｎ及びσ２ｍも閾値ｃ未満である。
【００３８】
次に、特徴量算出部１４は、ステップＡ４で推定された正規分布のうち、設定条件に合致する正規分布を選択し、選択した正規分布を用いて新たな正規分布を生成する（ステップＡ５）。
【００３９】
具体的には、カメラ２０からの画像においては、路面以外の物体（路面のペイント等）が映り込んでいる場合があるが、正確な水浮検知の観点からは、この物体を除去する必要がある。また、路面においては、日陰になっている部分と日向になっている部分とが生じている場合があるが、正確な水浮検知の観点からは、どちらかの部分を除去する必要がある。このため、特徴量算出部１４は、混合正規分布の中から、以下の条件（ｂ１）〜（ｂ３）のいずれかに合致する正規分布を選択する。また、以下に示す、閾値ｄ及び閾値ｅも、事前の実験に基づいて設定される。図６の例では、最も右側に位置している正規分布以外の正規分布が選択される。
【００４０】
条件（ｂ１）：分散σ２が閾値ｄより大きい。
条件（ｂ２）：全正規分布のうち頂点の値が最大（最大値）となる正規分布と裾が重なり、且つ、その頂点の値と最大値との差が閾値ｅより大きい（頂点の値が小さすぎない）。
【００４１】
なお、上記条件(ｂ２)において、裾が重なっているかどうかは、まず、各正規分布について、標準偏差σと定数Ｎ（Ｎａ、Ｎｂ）との積算値を算出し、算出した積算値が、以下の条件（ｃ１）に合致する場合に重なりがあると判定する。また、下記の状件（ｃ１）における定数Ｎａ、Ｎｂの値は、画像が出力された時間帯（昼間及び夜間）、及び正規分布の頂点の値の大きさ等に応じて適宜設定される。
【００４２】
条件（ｃ１）：｜｜μｎ−μｍ｜｜＞｜｜Ｎａ×σｎ−Ｎｂ×σｍ｜｜
【００４３】
次に、特徴量算出部１４は、選択した正規分布を用いて新たな正規分布Ｎ（μ、σ２）を生成し、生成した新たな正規分布の平均値及び分散を、特徴量として算出する（ステップＡ６）。
【００４４】
具体的には、ステップＡ６では、図７に示すように、特徴量算出部１４は、既存の手法を用いて、選択した正規分布を統合して、新たな正規分布Ｎ（μ、σ２）を生成する。図７は、本発明の実施の形態において、選択された正規分布から生成された新たな正規分布の一例を示す図である。なお、図７では、新たな正規分布Ｎ（μ、σ２）は「統合された正規分布」と表記されている。
【００４５】
また、新たな正規分布Ｎ（μ、σ２）は、下記の数１〜数３によって規定される。なお、下記の数１〜数３においては、統合対象となる正規分布の数をｋとし、各正規分布に１〜ｋまでの番号が付与されているとする。数１〜数３において、ｗ（ｋ）は、統合対象となる各正規分布の分布重みであり、μ（ｋ）は統合対象となる各正規分布の平均値である。また、σ２（ｋ）は統合対象となる各正規分布の分散であり、ｒ（ｋ）は統合対象となる各正規分布の混合率である。更に、δ（ｋ）は、統合対象となる各正規分布を表すデルタ関数である。
【００４６】
【数１】
2015088106.tif 000003
【００４７】
【数２】
2015088106.tif 000004
【００４８】
【数３】
2015088106.tif 000005
【００４９】
また、ステップＡ６では、特徴量算出部１４は、時間帯が夜間である場合は、選択した正規分布のうち、平均値μの値が設定値以上となる正規分布を除外して、新たな正規分布Ｎ（μ、σ２）を生成することができる。これは、夜間において平均値μが高すぎる正規分布は、ライトの反射光の影響を強く受けている部分を示している可能性があるからである。
【００５０】
また、特徴量算出部１４は、本実施の形態では、設定された昼間の時間帯においては、新たな正規分布Ｎ（μ、σ２）の分散σ２のみを、特徴量として算出することができる。更に、特徴量算出部１４は、設定された夜間の時間帯においては、新たな正規分布Ｎ（μ、σ２）の平均値μ及び分散値σ２を、特徴量として算出することができる。
【００５１】
次に、路面状況推定部１５は、ステップＡ６で算出された特徴量に基づいて、道路の路面に水浮が発生しているかどうかを判定し、判定結果を外部に出力する（ステップＡ７）。
【００５２】
具体的には、路面状況推定部１５は、ステップＡ６で算出された特徴量と、それ以前に算出されたフレームＸ枚分（例えば２９９枚分）の特徴量とを用い、時間帯に応じて、分散の移動平均のみを、又はこれと平均値の分布平均との両方を算出し、得られた値を蓄積する。
【００５３】
続いて、路面状況推定部１５は、蓄積している値の数がＹ個（例えば９９個）ある場合は、これの累積値を求め、求めた累積値と閾値ｆとを比較する。そして、累積値が閾値ｆを超えている場合は、路面状況推定部１５は、道路の路面上に水浮が発生していると判断する。なお、閾値ｆは、水浮が発生している路面の画像を用いた実験によって取得することができる。また、蓄積している値の数がＹ個に達していない場合は、達するまで、ステップＡ１〜Ａ６と、分布平均及び移動平均の算出処理とが繰り返し実行される。
【００５４】
また、ステップＡ７の終了後に、新たな画像がカメラ２０から出力されてくると、再度、ステップＡ１が実行される。ステップＡ１〜Ａ７は、カメラ２０のフレームレートに合せて繰り返し実行される。
【００５５】
以上のように、ステップＡ１〜Ａ７が実行されると、路面の水浮が生じている部分を最も良く表す正規分布（統合された正規分布）が得られ、この正規部分布の最新のものと、過去のものとを用いて、水浮が発生しているかどうかが判断される。このため、本実施の形態によれば、水浮の判断結果の信頼性は高いものとなる。また、事前に閾値を設定する必要はあるが、場面毎に状況を定義する必要がないため、カメラ２０が車両に搭載されている場合であっても、道路の路面上の水浮を確実に検知することができる。
【００５６】
本実施の形態におけるプログラムは、コンピュータに、図２に示すステップＡ１〜Ａ７を実行させるプログラムであれば良い。このプログラムをコンピュータにインストールし、実行することによって、本実施の形態における水浮検知装置１０と水浮検知方法とを実現することができる。この場合、コンピュータのＣＰＵ（Central Processing Unit）は、ヒストグラム生成部１１、初期値算出部１２、正規分布推定部１３、特徴量算出部１４、及び路面状況推定部１５として機能し、処理を行なう。
【００５７】
また、本実施の形態におけるプログラムを実行可能なコンピュータは、特に限定されず、例えば、車載用のコンピュータが挙げられる。また、コンピュータは、スマートフォン、タブレット型端末等の携帯型情報端末であっても良い。このようなコンピュータを水浮検知装置として利用することも可能であるからである。
【産業上の利用可能性】
【００５８】
以上のように、本発明によれば、カメラの映像から道路の路面上の水浮を検知することができ、カメラは車両に搭載されていても良い。本発明は、道路の路面上の水浮を検知する必要性がある車両のシステム、道路施設等に有用である。
【符号の説明】
【００５９】
１０  水浮検知装置
１１  ヒストグラム生成装置
１２  初期値算出部
１３  正規分布推定部
１４  特徴量算出部
１５  路面状況推定部
２０  カメラ
【図１】
2015088106.tif 000006
【図２】
2015088106.tif 000007
【図３】
2015088106.tif 000008
【図４】
2015088106.tif 000009
【図５】
2015088106.tif 000010
【図６】
2015088106.tif 000011
【図７】
2015088106.tif 000012