"動画画像について日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出装置であって、時系列で連続的に与えられるフレーム単位のカラー原画像を、三原色の各画素値を有する画素の集合体データとして入力する画像入力部と、入力された原画像を逐次格納する原画像格納部と、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した平均画像とに基づいて、新たな平均画像を作成する平均画像作成部と、作成された平均画像を逐次格納する平均画像格納部と、新たな原画像が入力されるたびに、直前に作成した参照画像に基づいて新たな参照画像を作成する参照画像作成部と、作成された参照画像を逐次格納する参照画像格納部と、前記原画像格納部に格納された原画像と、前記平均画像格納部に格納された平均画像と、前記参照画像格納部に格納された参照画像と、を用いた比較処理により、日照変化領域と非日照変化領域とを区別する日照変化画像を作成する画像比較部と、前記日照変化画像を格納する日照変化画像格納部と、前記日照変化画像を出力する画像出力部と、を備え、前記画像比較部は、前記原画像格納部に格納された第i番目の原画像、前記平均画像格納部に格納された第i-1番目の平均画像、前記参照画像格納部に格納された第i-1番目の参照画像から、それぞれ対応する所定位置にある画素の三原色の各画素値を読み出す画素値読出部と、前記画素値読出部によって読み出された画素値に基づいて、前記第i番目の原画像の前記所定位置にある画素が前景画素か背景画素かを判定するとともに、日照変化画素か非日照変化画素かを判定する条件判定部と、第i番目の日照変化画像を構成する前記所定位置の画素の画素値として、前記条件判定部が日照変化画素と判定した場合には日照変化領域を示す画素値を、前記条件判定部が非日照変化画素と判定した場合には非日照変化領域を示す画素値を、それぞれ前記日照変化画像格納部に書き込む画素値書込部と、を有し、前記条件判定部は、第i番目の原画像を構成する画素の画素値piと第i-1番目の平均画像を構成する画素の画素値ai-1とを比較して類似するか否かを判定する第1判定部と、第i番目の原画像を構成する画素の画素値piと第i-1番目の参照画像を構成する画素の画素値ui-1とを比較して類似するか否かを判定するとともに、非類似との判定がなされた場合には、更に、類似境界の所定近傍範囲内にあるか否かを判定する第2判定部と、前記第1判定部により類似との判定結果が得られ、かつ、前記第2判定部により非類似かつ所定近傍範囲内との判定結果が得られた場合には、日照変化画素、その他の場合には、非日照変化画素との判定を行うとともに、前記第1判定部による判定結果および前記第2判定部による判定結果の少なくとも一方が類似を示す場合には背景画素、双方とも非類似を示す場合には前景画素との判定を行う総合判定部と、を有し、前記平均画像作成部は、第1番目の原画像をそのまま第1番目の平均画像とし、第i番目の原画像と第i-1番目の平均画像とについて、互いに対応する位置にある画素の三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を各画素値とする画素の集合からなる画像を第i番目i≧2の平均画像とし、前記参照画像作成部は、第1番目の原画像をそのまま第1番目の参照画像とし、第i-1番目の参照画像を構成する個々の画素のうち、前記条件判定部によって第i番目の原画像の対応画素が前景画素と判定された画素については、更新せずに、もとの画素値をそのまま維持し、前記条件判定部によって第i番目の原画像の対応画素が背景画素と判定された画素については、当該画素と第i-1番目の平均画像の対応する位置にある画素とについて、三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を新画素値として更新することにより、第i番目i≧2の参照画像を作成することを特徴とする日照変化領域検出装置。請求項1に記載の日照変化領域検出装置において、平均画像作成部が、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の平均画像A1として平均画像格納部へ格納する処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の平均画像Aiを、ai=1-wai-1wpi但し、aiは、平均画像Aiの所定位置の画素の所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1なる演算式を用いて作成することを特徴とする日照変化領域検出装置。請求項2に記載の日照変化領域検出装置において、参照画像作成部が、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の参照画像U1として参照画像格納部へ格納する処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の参照画像Uiの各画素値を、背景画素と判定された画素については、ui=1-vui-1vai-1但し、uiは、平均画像Uiの所定位置の画素の所定色の画素値、ui-1は、平均画像Ui-1の前記所定位置の画素の前記所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、vは、所定の重みを示すパラメータv<1なる演算式を用いて求められた画素値uiに更新し、前景画素と判定された画素については、ui=ui-1として、更新せずに、もとの画素値をそのまま維持することにより決定することを特徴とする日照変化領域検出装置。請求項3に記載の日照変化領域検出装置において、重みを示すパラメータw,vを、w>vとなるように設定することを特徴とする日照変化領域検出装置。請求項2~4のいずれかに記載の日照変化領域検出装置において、パラメータwもしくはvまたはこれら双方を、ユーザの操作入力によって任意の値に設定するパラメータ設定部を更に備えることを特徴とする日照変化領域検出装置。請求項1~5のいずれかに記載の日照変化領域検出装置において、第1判定部および第2判定部が、比較対象となる一方の画素の画素値を基準画素値、他方の画素の画素値を比較画素値として、三原色の各画素値を各座標軸にとった三次元座標系において、前記基準画素値に対応する座標に位置する基準点Qと前記比較画素値に対応する座標に位置する比較点qをとり、前記基準点Qを中心とする所定サイズの基準立体と前記比較点qとの位置関係を調べ、前記比較点qが前記基準立体の外部に位置すると判定できる場合に非類似との判定を行い、内部に位置すると判定できる場合には類似との判定を行い、第2判定部が、前記基準点Qを中心として配置された前記基準立体よりも大きなサイズをもつ拡張基準立体と前記比較点qとの位置関係を調べ、前記比較点qが前記基準立体の外部、かつ、前記拡張基準立体の内部に位置すると判定できる場合に、類似境界の所定近傍範囲内にあるとの判定を行うことを特徴とする日照変化領域検出装置。請求項6に記載の日照変化領域検出装置において、第1判定部が、原画像の画素の画素値を基準画素値、平均画像の画素の画素値を比較画素値として判定を行い、第2判定部が、原画像の画素の画素値を基準画素値、参照画像の画素の画素値を比較画素値として判定を行うことを特徴とする日照変化領域検出装置。請求項6または7に記載の日照変化領域検出装置において、第1判定部および第2判定部が、三次元座標系の原点Oと基準点Qとの距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この楕円を前記基準軸Zを中心に回転させることにより得られる同一の回転楕円体Eもしくはその近似体を基準立体として用いることを特徴とする日照変化領域検出装置。請求項8に記載の日照変化領域検出装置において、第1判定部および第2判定部が、原点Oを中心とし半径D-αの球C1の内側に比較点qが位置するもしくは原点Oを中心とし半径Dαの球C2の外側に比較点qが位置するという第1の条件を満足するか否かを判定する第1の条件判定部と、前記第1の条件判定部が否定的判定を行った場合に原点Oを頂点、基準軸Zを中心軸とし、回転楕円体Eと円で接する円錐ξもしくはその近似体の外側に比較点qが位置するという第2の条件を満足するか否かを判定する第2の条件判定部と、前記第2の条件判定部が否定的判定を行った場合に回転楕円体Eの外側に比較点qが位置するという第3の条件を満足するか否かを判定する第3の条件判定部と、前記第1の条件判定部、前記第2の条件判定部、前記第3の条件判定部のいずれかが肯定的判定を行った場合に、非類似との判定結果を示す判定値を出力し、前記第1の条件判定部、前記第2の条件判定部、前記第3の条件判定部のすべてが否定的判定を行った場合に、類似との判定結果を示す判定値を出力する判定値出力部と、を有することを特徴とする日照変化領域検出装置。請求項9に記載の日照変化領域検出装置において、第2判定部の判定値出力部が、第1の条件判定部および第2の条件判定部が否定的判定を行い、第3の条件判定部が肯定的判定を行った場合に、類似境界の所定近傍範囲内との判定結果を示す判定値を出力し、その他の場合には、所定近傍範囲外との判定結果を示す判定値を出力することを特徴とする日照変化領域検出装置。請求項8に記載の日照変化領域検出装置において、第1判定部および第2判定部が、パラメータhおよびkの値を保持するパラメータ保持部と、画素値読出部が読み出した基準画素値R,G,Bおよび比較画素値r,g,bに基づいて、D2=R2G2B2およびd2=r2g2b2なる演算を行い、値D2および値d2を算出する基本演算部と、前記基本演算部が算出した値D2および値d2と前記パラメータ保持部に保持されているパラメータhの値とを用いて、1-h2D2>d2もしくは1h2D2<d2という条件を満足するか否かを判定するための演算を行う第1の条件判定部と、前記第1の条件判定部が否定的判定を行った場合に、画素値読出部が読み出した基準画素値R,G,Bおよび比較画素値r,g,bに基づいて、γ=RrGgBbなる演算を行い、値γを算出し、前記基本演算部が算出した値D2および値d2と前記パラメータ保持部に保持されているパラメータkの値と前記値γとを用いて、γ2/D2d2<1-k2/22という条件を満足するか否かを判定するための演算を行う第2の条件判定部と、前記第2の条件判定部が否定的判定を行った場合に、前記基本演算部が算出した値D2および値d2と前記パラメータ保持部に保持されているパラメータhおよびkの値と前記第2の条件判定部が算出した前記値γとを用いて、α2=h2D2β2=k2D2x2=D2γ2/D2-2γy2=d2-γ2/D2なる演算を行い、当該演算結果を用いて、β2x2α2y2-α2β2>0という条件を満足するか否かを判定するための演算を行う第3の条件判定部と、前記第1の条件判定部、前記第2の条件判定部、前記第3の条件判定部のいずれかが肯定的判定を行った場合に、非類似との判定結果を示す判定値を出力し、前記第1の条件判定部、前記第2の条件判定部、前記第3の条件判定部のすべてが否定的判定を行った場合に、類似との判定結果を示す判定値を出力する判定値出力部と、を有することを特徴とする日照変化領域検出装置。請求項11に記載の日照変化領域検出装置において、第2判定部の判定値出力部が、第1の条件判定部および第2の条件判定部が否定的判定を行い、第3の条件判定部が肯定的判定を行った場合に、類似境界の所定近傍範囲内との判定結果を示す判定値を出力し、その他の場合には、所定近傍範囲外との判定結果を示す判定値を出力することを特徴とする日照変化領域検出装置。請求項8~12のいずれかに記載の日照変化領域検出装置において、パラメータhおよびkを、ユーザの操作入力によって任意の値に設定するパラメータ設定部を更に備えることを特徴とする日照変化領域検出装置。請求項1~13のいずれかに記載の日照変化領域検出装置としてコンピュータを動作させるプログラム。請求項1~13のいずれかに記載の日照変化領域検出装置として機能する半導体集積回路。動画画像について移動物体を検出する移動物体検出装置であって、時系列で連続的に与えられるフレーム単位の原画像によって構成される共通動画画像について、それぞれ異なる判定基準を用いて移動物体を検出する複数n台の移動物体検出ユニットと、前記共通動画画像について、日照変化に起因した変動が生じている日照変化領域を随時検出する日照変化領域検出ユニットと、前記日照変化領域検出ユニットによって検出された日照変化領域の全領域に対する割合を示す日照変化率を随時算出する日照変化率算出ユニットと、前記日照変化率に基づいて、前記複数n台の移動物体検出ユニットのうちのいずれか1台を随時選択し、選択された移動物体検出ユニットが検出した移動物体を最終的な検出結果として出力する選択ユニットと、を備えることを特徴とする移動物体検出装置。請求項16に記載の移動物体検出装置において、個々の移動物体検出ユニットが、入力された原画像を、移動物体を含まないと推定される背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出し、n台の移動物体検出ユニットには、それぞれ互いに異なる類似判定基準が設定されており、選択ユニットが、日照変化率が高いほど、より緩い基準が設定されたユニットを選択することを特徴とする移動物体検出装置。請求項17に記載の移動物体検出装置において、各移動物体検出ユニットが、原画像と背景画像とを個々の画素ごとに比較し、比較対象となる一方の画素の画素値を基準画素値、他方の画素の画素値を比較画素値として、三原色の各画素値を各座標軸にとった三次元座標系において、前記基準画素値に対応する座標に位置する基準点Qと前記比較画素値に対応する座標に位置する比較点qをとり、前記基準点Qを中心とする所定サイズの基準立体と前記比較点qとの位置関係を調べ、前記比較点qが前記基準立体の外部に位置すると判定できる場合に非類似との判定を行い、内部に位置すると判定できる場合には類似との判定を行い、第j番目2≦j≦nの移動物体検出ユニットが採用する基準立体が、第j-1番目の移動物体検出ユニットが採用する基準立体を包含する、より大きなサイズとなっており、選択ユニットが、日照変化率が高いほど、よりサイズの大きな基準立体を採用するユニットを選択することを特徴とする移動物体検出装置。請求項18に記載の移動物体検出装置において、各移動物体検出ユニットが、原画像の画素の画素値を基準画素値、背景画像の画素の画素値を比較画素値として判定を行うことを特徴とする移動物体検出装置。請求項18または19に記載の移動物体検出装置において、各移動物体検出ユニットが、三次元座標系の原点Oと基準点Qとの距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この楕円を前記基準軸Zを中心に回転させることにより得られる回転楕円体Eもしくはその近似体を基準立体として用い、第j番目1≦j≦nの移動物体検出ユニットが採用する基準立体のパラメータh,kをそれぞれhj,kjとしたときに、hj>hj-1,kj>kj-1に設定されている但し、2≦j≦nことを特徴とする移動物体検出装置。請求項17~20のいずれかに記載の移動物体検出装置において、n=2に設定することにより、2台の移動物体検出ユニットを設け、第1の移動物体検出ユニットは第1の類似判定基準に基づく判定を行い、第2の移動物体検出ユニットは前記第1の類似判定基準よりも緩い第2の類似判定基準に基づく判定を行い、選択ユニットは、日照変化率算出ユニットによって算出された日照変化率が、所定のしきい値以下の場合には前記第1の移動物体検出ユニットを選択し、前記しきい値を超える場合には前記第2の移動物体検出ユニットを選択することを特徴とする移動物体検出装置。請求項16~21のいずれかに記載の移動物体検出装置において、日照変化領域検出ユニットとして、請求項1~13のいずれかに記載の日照変化領域検出装置を用いたことを特徴とする移動物体検出装置。請求項16~22のいずれかに記載の移動物体検出装置としてコンピュータを動作させるプログラム。請求項16~22のいずれかに記載の移動物体検出装置として機能する半導体集積回路。動画画像について日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出方法であって、コンピュータが、時系列で連続的に与えられるフレーム単位のカラー原画像を、三原色の各画素値を有する画素の集合体データとして入力する画像入力段階と、コンピュータが、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した平均画像とに基づいて、新たな平均画像を作成する平均画像作成段階と、コンピュータが、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した参照画像とに基づいて、新たな参照画像を作成する参照画像作成段階と、コンピュータが、原画像と平均画像と参照画像とを用いた比較処理により、日照変化領域と非日照変化領域とを区別する日照変化画像を作成する画像比較段階と、コンピュータが、前記日照変化画像を出力する画像出力段階と、を有し、前記画像比較段階は、第i番目の原画像、第i-1番目の平均画像、第i-1番目の参照画像から、それぞれ対応する所定位置にある画素の三原色の各画素値を読み出す画素値読出ステップと、第i番目の原画像を構成する画素の画素値piと第i-1番目の平均画像を構成する画素の画素値ai-1とが類似するか否かを判定する第1判定処理と、第i番目の原画像を構成する画素の画素値piと第i-1番目の参照画像を構成する画素の画素値ui-1とが類似するか否かを判定するとともに、非類似との判定がなされた場合には、更に、類似境界の所定近傍範囲内にあるか否かを判定する第2判定処理と、を行い、第i番目の原画像の前記所定位置にある画素について、前記第1判定処理により類似との判定結果が得られ、かつ、前記第2判定処理により非類似かつ所定近傍範囲内との判定結果が得られた場合には、日照変化画素、その他の場合には、非日照変化画素との判定を行うとともに、前記第1判定処理による判定結果と前記第2判定処理による判定結果との少なくとも一方が類似を示す場合には背景画素、双方とも非類似を示す場合には前景画素との判定を行う条件判定ステップと、第i番目の日照変化画像を構成する前記所定位置の画素の画素値として、前記条件判定ステップにより日照変化画素と判定された場合には日照変化領域を示す画素値を、前記条件判定ステップにより非日照変化画素と判定された場合には非日照変化領域を示す画素値を、それぞれ与えることにより、第i番目の日照変化画像を作成する日照変化画像作成ステップと、を有し、前記平均画像作成段階では、第1番目の原画像をそのまま第1番目の平均画像とし、第i番目i≧2の平均画像を作成する際には、第i番目の原画像と第i-1番目の平均画像とについて、互いに対応する位置にある画素の三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を各画素値とする画素の集合からなる画像を第i番目の平均画像とし、前記参照画像作成段階では、第1番目の原画像をそのまま第1番目の参照画像とし、第i番目i≧2の参照画像を作成する際には、第i-1番目の参照画像を構成する個々の画素のうち、前記条件判定ステップにおいて第i番目の原画像の対応画素が前景画素と判定された画素については、更新せずに、もとの画素値をそのまま維持し、前記条件判定ステップにおいて第i番目の原画像の対応画素が背景画素と判定された画素については、第i-1番目の平均画像の対応する位置にある画素との間で、三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を新画素値として更新することにより、第i番目の参照画像を作成することを特徴とする日照変化領域検出方法。請求項25に記載の日照変化領域検出方法において、平均画像作成段階では、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の平均画像A1とする処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の平均画像Aiを、ai=1-wai-1wpi但し、aiは、平均画像Aiの所定位置の画素の所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1なる演算式を用いて作成し、参照画像作成段階では、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の参照画像U1とする処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の参照画像Uiの各画素値を、背景画素と判定された画素については、ui=1-vui-1vai-1但し、uiは、平均画像Uiの所定位置の画素の所定色の画素値、ui-1は、平均画像Ui-1の前記所定位置の画素の前記所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、vは、所定の重みを示すパラメータv<1、かつw>vなる演算式を用いて求められた画素値uiに更新し、前景画素と判定された画素については、ui=ui-1として、更新せずに、もとの画素値をそのまま維持することにより、参照画像Uiを作成することを特徴とする日照変化領域検出方法。動画画像について移動物体を検出する移動物体検出方法であって、コンピュータが、時系列で連続的に与えられるフレーム単位の原画像によって構成される共通動画画像について、複数n通りの異なる判定基準を用いて複数n通りの移動物体検出結果を得る個別検出段階と、コンピュータが、前記共通動画画像について、日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出段階と、コンピュータが、前記日照変化領域検出段階によって検出された日照変化領域の全領域に対する割合を示す日照変化率を算出する日照変化率算出段階と、コンピュータが、前記日照変化率に基づいて、前記複数n通りの移動物体検出結果のいずれか1つを選択し、選択された検出結果を最終的な検出結果として出力する検出結果出力段階と、を有することを特徴とする移動物体検出方法。請求項27に記載の移動物体検出方法において、個別検出段階では、入力された原画像を、移動物体を含まないと推定される背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出する処理を、互いに異なる複数n通りの類似判定基準を用いて実行し、検出結果出力段階では、日照変化率が高いほど、より緩い類似判定基準を用いて得られた検出結果を選択することを特徴とする移動物体検出方法。コンピュータが、時系列で連続的に与えられるフレーム単位の原画像によって構成される動画画像について、入力された原画像を背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出する移動物体検出方法において、前記動画画像について、日照変化に起因した変動が生じている日照変化領域を検出し、検出した日照変化領域の全領域に対する割合を示す日照変化率を算出し、この日照変化率が高いほど、より緩い類似判定基準を用いた検出を行うことを特徴とする移動物体検出方法。上述したように、移動物体を検出する基本原理は、新たな入力画像を背景画像と比較し、時間的に変化している領域を認識することにある。しかしながら、画像の時間変化は、必ずしも移動物体に起因するものではなく、雲の動きによる明るさの変化、風の作用による樹木の揺れ、波の動きによる海面の変化など、様々な自然現象にも起因する。このため、特に、屋外の監視システムなどから取り込まれた動画画像に基づいて、移動物体の検出を行う場合、このような自然現象による画像変化を移動物体による画像変化と混同した誤検出が生じやすい。前掲の特許文献1,2には、このような問題を解決するための一手法が開示されているが、入力画像に対して、小領域もしくはブロック単位で画素値のヒストグラムをとるなどの統計処理を施し、当該統計処理によって得られた特徴量を複雑なモデルに当てはめて移動物体の領域を推定する必要があるため、演算処理が極めて複雑になるという問題がある。このため、一般的なプロセッサでは処理時間が長くかかり、リアルタイムでの処理を行うことが困難である。また、より正確な検出を行うためには、多数のフレームにわたって入力画像を格納しておく必要があり、大容量の画像バッファメモリが必要になる。このため、リアルタイムでの監視を目的とした安価な監視システムに利用するには不適当である。一方、上述した非特許文献1に開示されている手法では、画素単位での比較を行えばよいので、複雑な統計処理を行ったり、複雑なモデルを構築したりする必要はない。しかしながら、領域単位での統計処理を省き、画素単位での比較結果のみによって検出を行うため、誤検出が生じやすい。このような誤検出を生む最も大きな要因は、上述したとおり、雲の動き、樹木の揺れ、波の動きなどの自然現象である。一般に、移動物体が存在しない状況であっても、雲の動き、樹木の揺れ、波の動きなどの自然現象によって画像に変動が生じ、その結果、移動物体か否かの判定に誤りが生じることになる。これらの自然現象による画像の変動は、結局、被写体に対する太陽光による照明状態が変動することに起因する。そこで、本願では、このような太陽光による被写体に対する照明状態の短期的な変動を日照変化と呼ぶことにする。地下室や遮光された部屋内で撮影された動画画像であれば、このような日照変化の影響を受けることはないが、屋外で撮影された動画画像や、窓のある室内で撮影された動画画像では、日照変化の影響を避けることができない。このように、動画画像について日照変化に起因した変動を解析することは、移動物体の検出のみならず、当該動画画像の特徴を把握する上で非常に重要である。しかしながら、従来、日照変化に起因した変動に対する解析は十分に行われておらず、当然ながら、移動物体を検出する上で、日照変化に起因した誤検出に対する対策も十分ではない。本発明の第1の目的は、与えられた動画画像について、日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出装置を提供することにあり、第2の目的は、この日照変化領域の検出結果を利用して、より正確な検出を行うことが可能な移動物体検出装置を提供することにある。本発明の第1の態様は、動画画像について日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出装置において、時系列で連続的に与えられるフレーム単位のカラー原画像を、三原色の各画素値を有する画素の集合体データとして入力する画像入力部と、入力された原画像を逐次格納する原画像格納部と、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した平均画像とに基づいて、新たな平均画像を作成する平均画像作成部と、作成された平均画像を逐次格納する平均画像格納部と、新たな原画像が入力されるたびに、直前に作成した参照画像に基づいて新たな参照画像を作成する参照画像作成部と、作成された参照画像を逐次格納する参照画像格納部と、原画像格納部に格納された原画像と、平均画像格納部に格納された平均画像と、参照画像格納部に格納された参照画像と、を用いた比較処理により、日照変化領域と非日照変化領域とを区別する日照変化画像を作成する画像比較部と、日照変化画像を格納する日照変化画像格納部と、日照変化画像を出力する画像出力部と、を設け、画像比較部は、原画像格納部に格納された第i番目の原画像、平均画像格納部に格納された第i-1番目の平均画像、参照画像格納部に格納された第i-1番目の参照画像から、それぞれ対応する所定位置にある画素の三原色の各画素値を読み出す画素値読出部と、画素値読出部によって読み出された画素値に基づいて、第i番目の原画像の所定位置にある画素が前景画素か背景画素かを判定するとともに、日照変化画素か非日照変化画素かを判定する条件判定部と、第i番目の日照変化画像を構成する所定位置の画素の画素値として、条件判定部が日照変化画素と判定した場合には日照変化領域を示す画素値を、条件判定部が非日照変化画素と判定した場合には非日照変化領域を示す画素値を、それぞれ日照変化画像格納部に書き込む画素値書込部と、を有するようにし、条件判定部は、第i番目の原画像を構成する画素の画素値piと第i-1番目の平均画像を構成する画素の画素値ai-1とを比較して類似するか否かを判定する第1判定部と、第i番目の原画像を構成する画素の画素値piと第i-1番目の参照画像を構成する画素の画素値ui-1とを比較して類似するか否かを判定するとともに、非類似との判定がなされた場合には、更に、類似境界の所定近傍範囲内にあるか否かを判定する第2判定部と、第1判定部により類似との判定結果が得られ、かつ、第2判定部により非類似かつ所定近傍範囲内との判定結果が得られた場合には、日照変化画素、その他の場合には、非日照変化画素との判定を行うとともに、第1判定部による判定結果および第2判定部による判定結果の少なくとも一方が類似を示す場合には背景画素、双方とも非類似を示す場合には前景画素との判定を行う総合判定部と、を有するようにし、平均画像作成部は、第1番目の原画像をそのまま第1番目の平均画像とし、第i番目の原画像と第i-1番目の平均画像とについて、互いに対応する位置にある画素の三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を各画素値とする画素の集合からなる画像を第i番目i≧2の平均画像とし、参照画像作成部は、第1番目の原画像をそのまま第1番目の参照画像とし、第i-1番目の参照画像を構成する個々の画素のうち、条件判定部によって第i番目の原画像の対応画素が前景画素と判定された画素については、更新せずに、もとの画素値をそのまま維持し、条件判定部によって第i番目の原画像の対応画素が背景画素と判定された画素については、当該画素と第i-1番目の平均画像の対応する位置にある画素とについて、三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を新画素値として更新することにより、第i番目i≧2の参照画像を作成するようにしたものである。本発明の第2の態様は、上述した第1の態様に係る日照変化領域検出装置において、平均画像作成部が、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の平均画像A1として平均画像格納部へ格納する処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の平均画像Aiを、ai=1-wai-1wpi但し、aiは、平均画像Aiの所定位置の画素の所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1なる演算式を用いて作成するようにしたものである。本発明の第3の態様は、上述した第2の態様に係る日照変化領域検出装置において、参照画像作成部が、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の参照画像U1として参照画像格納部へ格納する処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の参照画像Uiの各画素値を、背景画素と判定された画素については、ui=1-vui-1vai-1但し、uiは、平均画像Uiの所定位置の画素の所定色の画素値、ui-1は、平均画像Ui-1の前記所定位置の画素の前記所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、vは、所定の重みを示すパラメータv<1なる演算式を用いて求められた画素値uiに更新し、前景画素と判定された画素については、ui=ui-1として、更新せずに、もとの画素値をそのまま維持することにより決定するようにしたものである。本発明の第4の態様は、上述した第3の態様に係る日照変化領域検出装置において、重みを示すパラメータw,vを、w>vとなるように設定したものである。本発明の第5の態様は、上述した第2~第4の態様に係る日照変化領域検出装置において、パラメータwもしくはvまたはこれら双方を、ユーザの操作入力によって任意の値に設定するパラメータ設定部を更に備えるようにしたものである。本発明の第6の態様は、上述した第1~第5の態様に係る日照変化領域検出装置において、第1判定部および第2判定部が、比較対象となる一方の画素の画素値を基準画素値、他方の画素の画素値を比較画素値として、三原色の各画素値を各座標軸にとった三次元座標系において、基準画素値に対応する座標に位置する基準点Qと比較画素値に対応する座標に位置する比較点qをとり、基準点Qを中心とする所定サイズの基準立体と比較点qとの位置関係を調べ、比較点qが基準立体の外部に位置すると判定できる場合に非類似との判定を行い、内部に位置すると判定できる場合には類似との判定を行い、第2判定部が、基準点Qを中心として配置された基準立体よりも大きなサイズをもつ拡張基準立体と比較点qとの位置関係を調べ、比較点qが基準立体の外部、かつ、拡張基準立体の内部に位置すると判定できる場合に、類似境界の所定近傍範囲内にあるとの判定を行うようにしたものである。本発明の第7の態様は、上述した第6の態様に係る日照変化領域検出装置において、第1判定部が、原画像の画素の画素値を基準画素値、平均画像の画素の画素値を比較画素値として判定を行い、第2判定部が、原画像の画素の画素値を基準画素値、参照画像の画素の画素値を比較画素値として判定を行うようにしたものである。本発明の第8の態様は、上述した第6または第7の態様に係る日照変化領域検出装置において、第1判定部および第2判定部が、三次元座標系の原点Oと基準点Qとの距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この楕円を基準軸Zを中心に回転させることにより得られる同一の回転楕円体Eもしくはその近似体を基準立体として用いるようにしたものである。本発明の第9の態様は、上述した第8の態様に係る日照変化領域検出装置において、第1判定部および第2判定部が、原点Oを中心とし半径D-αの球C1の内側に比較点qが位置するもしくは原点Oを中心とし半径Dαの球C2の外側に比較点qが位置するという第1の条件を満足するか否かを判定する第1の条件判定部と、第1の条件判定部が否定的判定を行った場合に原点Oを頂点、基準軸Zを中心軸とし、回転楕円体Eと円で接する円錐ξもしくはその近似体の外側に比較点qが位置するという第2の条件を満足するか否かを判定する第2の条件判定部と、第2の条件判定部が否定的判定を行った場合に回転楕円体Eの外側に比較点qが位置するという第3の条件を満足するか否かを判定する第3の条件判定部と、第1の条件判定部、第2の条件判定部、第3の条件判定部のいずれかが肯定的判定を行った場合に、非類似との判定結果を示す判定値を出力し、第1の条件判定部、第2の条件判定部、第3の条件判定部のすべてが否定的判定を行った場合に、類似との判定結果を示す判定値を出力する判定値出力部と、を有するようにしたものである。本発明の第10の態様は、上述した第9の態様に係る日照変化領域検出装置において、第2判定部の判定値出力部が、第1の条件判定部および第2の条件判定部が否定的判定を行い、第3の条件判定部が肯定的判定を行った場合に、類似境界の所定近傍範囲内との判定結果を示す判定値を出力し、その他の場合には、所定近傍範囲外との判定結果を示す判定値を出力するようにしたものである。本発明の第11の態様は、上述した第8の態様に係る日照変化領域検出装置において、第1判定部および第2判定部が、パラメータhおよびkの値を保持するパラメータ保持部と、画素値読出部が読み出した基準画素値R,G,Bおよび比較画素値r,g,bに基づいて、D2=R2G2B2およびd2=r2g2b2なる演算を行い、値D2および値d2を算出する基本演算部と、基本演算部が算出した値D2および値d2とパラメータ保持部に保持されているパラメータhの値とを用いて、1-h2D2>d2もしくは1h2D2<d2という条件を満足するか否かを判定するための演算を行う第1の条件判定部と、第1の条件判定部が否定的判定を行った場合に、画素値読出部が読み出した基準画素値R,G,Bおよび比較画素値r,g,bに基づいて、γ=RrGgBbなる演算を行い、値γを算出し、基本演算部が算出した値D2および値d2とパラメータ保持部に保持されているパラメータkの値と値γとを用いて、γ2/D2d2<1-k2/22という条件を満足するか否かを判定するための演算を行う第2の条件判定部と、第2の条件判定部が否定的判定を行った場合に、基本演算部が算出した値D2および値d2とパラメータ保持部に保持されているパラメータhおよびkの値と第2の条件判定部が算出した値γとを用いて、α2=h2D2β2=k2D2x2=D2γ2/D2-2γy2=d2-γ2/D2なる演算を行い、当該演算結果を用いて、β2x2α2y2-α2β2>0という条件を満足するか否かを判定するための演算を行う第3の条件判定部と、第1の条件判定部、第2の条件判定部、第3の条件判定部のいずれかが肯定的判定を行った場合に、非類似との判定結果を示す判定値を出力し、第1の条件判定部、第2の条件判定部、第3の条件判定部のすべてが否定的判定を行った場合に、類似との判定結果を示す判定値を出力する判定値出力部と、を有するようにしたものである。本発明の第12の態様は、上述した第11の態様に係る日照変化領域検出装置において、第2判定部の判定値出力部が、第1の条件判定部および第2の条件判定部が否定的判定を行い、第3の条件判定部が肯定的判定を行った場合に、類似境界の所定近傍範囲内との判定結果を示す判定値を出力し、その他の場合には、所定近傍範囲外との判定結果を示す判定値を出力するようにしたものである。本発明の第13の態様は、上述した第8~12の態様に係る日照変化領域検出装置において、パラメータhおよびkを、ユーザの操作入力によって任意の値に設定するパラメータ設定部を更に備えるようにしたものである。本発明の第14の態様は、上述した第1~13の態様に係る日照変化領域検出装置を、コンピュータに専用のプログラムを組み込むことにより構成したものである。本発明の第15の態様は、上述した第1~13の態様に係る日照変化領域検出装置を、半導体集積回路によって構成したものである。本発明の第16の態様は、動画画像について移動物体を検出する移動物体検出装置において、時系列で連続的に与えられるフレーム単位の原画像によって構成される共通動画画像について、それぞれ異なる判定基準を用いて移動物体を検出する複数n台の移動物体検出ユニットと、共通動画画像について、日照変化に起因した変動が生じている日照変化領域を随時検出する日照変化領域検出ユニットと、日照変化領域検出ユニットによって検出された日照変化領域の全領域に対する割合を示す日照変化率を随時算出する日照変化率算出ユニットと、日照変化率に基づいて、複数n台の移動物体検出ユニットのうちのいずれか1台を随時選択し、選択された移動物体検出ユニットが検出した移動物体を最終的な検出結果として出力する選択ユニットと、を設けたものである。本発明の第17の態様は、上述した第16の態様に係る移動物体検出装置において、個々の移動物体検出ユニットが、入力された原画像を、移動物体を含まないと推定される背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出し、n台の移動物体検出ユニットには、それぞれ互いに異なる類似判定基準が設定されており、選択ユニットが、日照変化率が高いほど、より緩い基準が設定されたユニットを選択するようにしたものである。本発明の第18の態様は、上述した第17の態様に係る移動物体検出装置において、各移動物体検出ユニットが、原画像と背景画像とを個々の画素ごとに比較し、比較対象となる一方の画素の画素値を基準画素値、他方の画素の画素値を比較画素値として、三原色の各画素値を各座標軸にとった三次元座標系において、基準画素値に対応する座標に位置する基準点Qと比較画素値に対応する座標に位置する比較点qをとり、基準点Qを中心とする所定サイズの基準立体と比較点qとの位置関係を調べ、比較点qが基準立体の外部に位置すると判定できる場合に非類似との判定を行い、内部に位置すると判定できる場合には類似との判定を行い、第j番目2≦j≦nの移動物体検出ユニットが採用する基準立体が、第j-1番目の移動物体検出ユニットが採用する基準立体を包含する、より大きなサイズとなっており、選択ユニットが、日照変化率が高いほど、よりサイズの大きな基準立体を採用するユニットを選択するようにしたものである。本発明の第19の態様は、上述した第18の態様に係る移動物体検出装置において、各移動物体検出ユニットが、原画像の画素の画素値を基準画素値、背景画像の画素の画素値を比較画素値として判定を行うようにしたものである。本発明の第20の態様は、上述した第18または第19の態様に係る移動物体検出装置において、各移動物体検出ユニットが、三次元座標系の原点Oと基準点Qとの距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この楕円を基準軸Zを中心に回転させることにより得られる回転楕円体Eもしくはその近似体を基準立体として用い、第j番目1≦j≦nの移動物体検出ユニットが採用する基準立体のパラメータh,kをそれぞれhj,kjとしたときに、hj>hj-1,kj>kj-1に設定されている但し、2≦j≦nようにしたものである。本発明の第21の態様は、上述した第17~第20の態様に係る移動物体検出装置において、n=2に設定することにより、2台の移動物体検出ユニットを設け、第1の移動物体検出ユニットは第1の類似判定基準に基づく判定を行い、第2の移動物体検出ユニットは第1の類似判定基準よりも緩い第2の類似判定基準に基づく判定を行い、選択ユニットは、日照変化率算出ユニットによって算出された日照変化率が、所定のしきい値以下の場合には第1の移動物体検出ユニットを選択し、しきい値を超える場合には第2の移動物体検出ユニットを選択するようにしたものである。本発明の第22の態様は、上述した第16~第21の態様に係る移動物体検出装置において、日照変化領域検出ユニットとして、上述した第1~第13の対応に係る日照変化領域検出装置を用いるようにしたものである。本発明の第23の態様は、上述した第16~第22の態様に係る移動物体検出装置を、コンピュータに専用のプログラムを組み込むことにより構成したものである。本発明の第24の態様は、上述した第16~22の態様に係る移動物体検出装置を、半導体集積回路によって構成したものである。本発明の第25の態様は、動画画像について日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出方法において、コンピュータが、時系列で連続的に与えられるフレーム単位のカラー原画像を、三原色の各画素値を有する画素の集合体データとして入力する画像入力段階と、コンピュータが、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した平均画像とに基づいて、新たな平均画像を作成する平均画像作成段階と、コンピュータが、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した参照画像とに基づいて、新たな参照画像を作成する参照画像作成段階と、コンピュータが、原画像と平均画像と参照画像とを用いた比較処理により、日照変化領域と非日照変化領域とを区別する日照変化画像を作成する画像比較段階と、コンピュータが、日照変化画像を出力する画像出力段階と、を行い、画像比較段階を、第i番目の原画像、第i-1番目の平均画像、第i-1番目の参照画像から、それぞれ対応する所定位置にある画素の三原色の各画素値を読み出す画素値読出ステップと、第i番目の原画像を構成する画素の画素値piと第i-1番目の平均画像を構成する画素の画素値ai-1とが類似するか否かを判定する第1判定処理と、第i番目の原画像を構成する画素の画素値piと第i-1番目の参照画像を構成する画素の画素値ui-1とが類似するか否かを判定するとともに、非類似との判定がなされた場合には、更に、類似境界の所定近傍範囲内にあるか否かを判定する第2判定処理と、を行い、第i番目の原画像の所定位置にある画素について、第1判定処理により類似との判定結果が得られ、かつ、第2判定処理により非類似かつ所定近傍範囲内との判定結果が得られた場合には、日照変化画素、その他の場合には、非日照変化画素との判定を行うとともに、第1判定処理による判定結果と第2判定処理による判定結果との少なくとも一方が類似を示す場合には背景画素、双方とも非類似を示す場合には前景画素との判定を行う条件判定ステップと、第i番目の日照変化画像を構成する所定位置の画素の画素値として、条件判定ステップにより日照変化画素と判定された場合には日照変化領域を示す画素値を、条件判定ステップにより非日照変化画素と判定された場合には非日照変化領域を示す画素値を、それぞれ与えることにより、第i番目の日照変化画像を作成する日照変化画像作成ステップと、によって構成し、平均画像作成段階では、第1番目の原画像をそのまま第1番目の平均画像とし、第i番目i≧2の平均画像を作成する際には、第i番目の原画像と第i-1番目の平均画像とについて、互いに対応する位置にある画素の三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を各画素値とする画素の集合からなる画像を第i番目の平均画像とし、参照画像作成段階では、第1番目の原画像をそのまま第1番目の参照画像とし、第i番目i≧2の参照画像を作成する際には、第i-1番目の参照画像を構成する個々の画素のうち、条件判定ステップにおいて第i番目の原画像の対応画素が前景画素と判定された画素については、更新せずに、もとの画素値をそのまま維持し、条件判定ステップにおいて第i番目の原画像の対応画素が背景画素と判定された画素については、第i-1番目の平均画像の対応する位置にある画素との間で、三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を新画素値として更新することにより、第i番目の参照画像を作成するようにしたものである。本発明の第26の態様は、上述した第25の態様に係る日照変化領域検出方法において、平均画像作成段階では、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の平均画像A1とする処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の平均画像Aiを、ai=1-wai-1wpi但し、aiは、平均画像Aiの所定位置の画素の所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1なる演算式を用いて作成し、参照画像作成段階では、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の参照画像U1とする処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の参照画像Uiの各画素値を、背景画素と判定された画素については、ui=1-vui-1vai-1但し、uiは、平均画像Uiの所定位置の画素の所定色の画素値、ui-1は、平均画像Ui-1の前記所定位置の画素の前記所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、vは、所定の重みを示すパラメータv<1、かつw>vなる演算式を用いて求められた画素値uiに更新し、前景画素と判定された画素については、ui=ui-1として、更新せずに、もとの画素値をそのまま維持することにより、参照画像Uiを作成するようにしたものである。本発明の第27の態様は、動画画像について移動物体を検出する移動物体検出方法において、コンピュータが、時系列で連続的に与えられるフレーム単位の原画像によって構成される共通動画画像について、複数n通りの異なる判定基準を用いて複数n通りの移動物体検出結果を得る個別検出段階と、コンピュータが、共通動画画像について、日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出段階と、コンピュータが、日照変化領域検出段階によって検出された日照変化領域の全領域に対する割合を示す日照変化率を算出する日照変化率算出段階と、コンピュータが、日照変化率に基づいて、複数n通りの移動物体検出結果のいずれか1つを選択し、選択された検出結果を最終的な検出結果として出力する検出結果出力段階と、を行うようにしたものである。本発明の第28の態様は、上述した第27の態様に係る移動物体検出方法において、個別検出段階では、入力された原画像を、移動物体を含まないと推定される背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出する処理を、互いに異なる複数n通りの類似判定基準を用いて実行し、検出結果出力段階では、日照変化率が高いほど、より緩い類似判定基準を用いて得られた検出結果を選択するようにしたものである。本発明の第29の態様は、コンピュータが、時系列で連続的に与えられるフレーム単位の原画像によって構成される動画画像について、入力された原画像を背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出する移動物体検出方法において、動画画像について、日照変化に起因した変動が生じている日照変化領域を検出し、検出した日照変化領域の全領域に対する割合を示す日照変化率を算出し、この日照変化率が高いほど、より緩い類似判定基準を用いた検出を行うようにしたものである。本発明に係る日照変化領域検出装置では、過去に入力された複数の原画像の平均的な特徴を有する平均画像が逐次作成され、更に、この平均画像に基づいて参照画像が逐次作成される。そして、新たに入力された原画像は、平均画像と比較されるとともに参照画像とも比較され、両方の比較結果がともに非類似とされた領域が前景領域移動物体の領域、それ以外の領域が背景領域と判定される。ここで、平均画像は、新たな原画像が入力されるたびに更新されるため、原画像の変化に追随して変化を遂げる動的画像になる。これに対して、参照画像は、背景領域と判定された部分に対しては平均画像を利用した更新が行われるが、前景領域と判定された部分に対しては更新が行われない。このため、原画像の変化に追随した参照画像の変化は緩慢になる。このように、平均画像と参照画像とは、原画像の変化に対する追随態様が互いに異なっている。しかも、本願発明者が行った実験によると、新たに入力された原画像のうち、平均画像には類似しており、参照画像には非類似であるが類似境界の所定近傍範囲内にあるという条件を満たす領域は、日照変化に起因した変動が生じている領域である可能性が高い。そこで、上記条件を満たす領域を日照変化領域として出力することにより、日照変化領域の検出を行うことができるようになる。特に、回転楕円体を用いて画素の類否判定を行う実施形態の場合、三次元色空間上において、一方の画素値をプロットした座標点を中心とする回転楕円体内に、他方の画素値をプロットした座標点が入るか否かに基づいて、比較対象となるカラー画素の類否判定を行うようにしたため、演算処理の負担を大幅に軽減することが可能になる。一方、本発明に係る移動物体検出装置では、それぞれ異なる判定基準を用いて移動物体を検出する複数n台の移動物体検出ユニットが並列動作する。しかも、日照変化領域検出ユニットにより、その時点における日照変化領域が検出され、日照変化率算出ユニットにより、その時点における日照変化率が算出される。この日照変化率は、その時点における日照変化に起因した画像変動の程度を示すパラメータとなる。そこで、当該日照変化率に基づいて、その時点で最適と推定されるユニットの検出結果を選択し、これを最終的な検出結果として出力するようにすれば、日照変化が生じていた場合でも、当該日照変化の程度に応じた最適な判定基準に基づく検出結果を得ることができるようになり、より正確な移動物体の検出が可能になる。特に、個々の移動物体検出ユニットに、原画像と背景画像との類似度を判定させ、類似度が低い領域を移動物体の領域として検出する実施形態をとる場合は、個々の移動物体検出ユニットごとに、それぞれ異なる類似判定基準を設定しておき、日照変化率が高い時点ほど、より緩い基準が設定されたユニットを選択するようにすればよい。そうすれば、日照変化が激しい時点では、類似判定基準が緩くなり、類似度が低い部分でも類似と判定され、背景領域として認識されることになる。その結果、日照変化に基づいて変動が生じた背景領域が、誤って移動物体の領域前景領域として誤検出されるのを防ぐことができる。逆に、日照変化が穏やかな時点では、類似判定基準が厳しくなり、類似度が低い部分は非類似と判定され、移動物体の領域前景領域として認識されることになる。その結果、移動物体の領域が誤って背景領域として誤検出されるのを防ぐことができる。以下、本発明を図示する実施形態に基づいて説明する。本発明は、日照変化領域検出装置および移動物体検出装置に係るものであるが、ここでは、まず、§1~§3において、従来の一般的な移動物体検出装置について述べ、§4~§9において、改良型の移動物体検出装置について述べる。この改良型の移動物体検出装置は、特願2011-020228号明細書に開示されている発明であるが、本発明に係る日照変化領域検出装置に密接に関連した構成を有しているため、説明の便宜上、まず、この改良型の移動物体検出装置についての説明を行うことにする。そして、続く§10において、本発明に係る日照変化領域検出装置の基本原理を説明し、§11において、本発明に係る日照変化領域検出装置の基本構成を述べ、§12において、本発明に係る日照変化領域検出装置における具体的処理を述べる。また、§13では、本発明に係る移動物体検出装置の基本構成を述べ、最後の§14では、本発明を方法として捉えた場合の説明を行う。<<<§1.一般的な移動物体検出の基本原理>>>はじめに、一般的な移動物体検出の基本原理を説明する。この基本原理は、動画画像について移動物体を検出するためのものであり、三原色の各画素値を有する画素の集合体として、時系列で連続的にフレーム単位のカラー画像が与えられることを前提としている。たとえば、図1の左側に示すように、フレーム単位の画像P10,P20,P30,P40がこの順番で与えられたものとしよう。図示の例は、街頭に設置された定点ビデオカメラによる撮影画像を示すものであり、右から左へと1台の車両移動物体が通過する状態が示されている。なお、一般的なビデオカメラの場合、1秒間に30フレーム程度の周期で連続して撮影画像が出力されるので、実際には、より細かな時間間隔で多数のフレーム画像が得られることになるが、ここでは、説明の便宜上、図示のとおり、4枚のフレーム画像P10,P20,P30,P40が順番に与えられたものとしよう。これら4枚の画像を見ると、画像P10は街の背景のみであるが、画像P20では、右側に車両の一部が侵入したため、背景上に車両の一部が重なった画像になっている。同様に、画像P30,P40も、背景上に車両が重なった画像になっている。背景は静止しているため、画像P10~P40について共通になるが、車両は移動物体であるため、画像P10~P40では異なる。ここでは、各画像上において、背景を構成する領域を背景領域B、車両移動物体を構成する領域を前景領域Fと呼ぶことにする。ここで、移動物体検出の目的は、与えられた各画像P10~P40について、それぞれ背景領域Bと前景領域Fとを区別することにある。図1の右側に示す各画像M10~M40は、それぞれ左側に示す各画像P10~P40について、背景領域Bと前景領域Fとの区別を示す画像である。図にハッチングを施した領域が背景領域Bであり、白地の領域が前景領域Fである。いずれの領域かを識別できればよいので、これらの画像M10~M40は、二値画像として与えられる。たとえば、背景領域Bに所属する画素には画素値1を与え、前景領域Fに所属する画素には画素値0を与えることにすれば、画像M10~M40は、1ビットの画素値をもった画素の集合体からなる画像データとして用意することができる。一般的な移動物体検出装置によれば、与えられた個々のフレーム画像のそれぞれについて、対応する二値画像を作成することができる。結局、移動物体検出処理の本質は、時系列で連続的に与えられた個々のカラー画像について、それぞれ背景領域Bと前景領域Fとの区別を示す二値画像を作成する処理ということができる。そこで、以下、便宜上、時系列で連続的に与えられるカラー画像を原画像と呼び、この原画像に基づいて作成される二値画像をマスク画像と呼ぶことにする。図1に示す例では、各原画像P10~P40のそれぞれについて、マスク画像M10~M40が作成された例が示されている。これらマスク画像を利用すれば、原画像上のどの領域が移動物体であるかを認識できるので、移動物体の形状やサイズを把握したり、移動物体の部分をズームアップして、車両のナンバーを確認したり、人物を特定したりする処理を行うことが可能になる。また、画面上で移動物体をズームアップしたまま追跡したり、移動経路を把握したりすることも可能になる。もちろん、このようなマスク画像の作成を、1枚の原画像のみから行うことは困難である。たとえば、図示の原画像P20のみが静止画像として与えられた場合、人間の脳は様々な情報から車両の存在を認識することができるが、この1枚の静止画像のみから、コンピュータに車両の存在を認識させるには、複雑なアルゴリズムに基づく処理が必要となる。そこで、一般的な移動物体検出アルゴリズムでは、既に述べたとおり、入力画像と背景画像との差を求めることにより、当該差の領域を移動物体の領域前景領域Fと判断する手法が採られる。たとえば、予め原画像P10を背景画像として指定しておけば、原画像P20が入力された時点で両者を比較することにより、差の領域異なる領域を前景領域Fと認識し、図示のようなマスク画像M20を作成することが可能になる。差の領域の認識手法としては、前述したとおり、小領域単位で画素値の統計処理を施す方法も知られているが、入力画像と背景画像とを画素単位で比較して、両画素の画素値が類似している場合、当該画素は背景領域Bを構成する画素であるとし、類似していない場合、当該画素は移動物体を構成する前景領域Fを構成する画素であるとする方法が最も簡便な方法である。背景画像としては、予め移動物体が存在しない状態たとえば、画像P10に示す状態で撮影した静止画像を利用することも可能であるが、天候や時刻などの要因により、照明環境は時々刻々と変化してゆくものであり、ある特定の時点で撮影した静止画像をそのまま背景画像として継続して利用するのは不適切である。そこで、通常は、過去に入力された複数の原画像に基づいて、これら原画像の平均的な特徴を有する平均画像を作成し、この平均画像を背景画像として利用する手法が採られる。平均画像は新たな画像が入力されるたびに逐次更新されるので、天候や時刻などの要因により、照明環境が変化したとしても、当該変化に追随した適切な背景画像として機能する。図2は、一般的な移動物体検出の基本原理を示す図である。ここで、図の上段に示す画像P1,P2,,Pi-2,Pi-1は、時系列で連続的に与えられるフレーム単位の原画像であり、これまでに第1番目の原画像P1から第i-1番目の原画像Pi-1まで、i-1フレーム分の画像が与えられた状態が示されている。一方、図の中段左に示す画像Piは、新たに与えられた第i番目の原画像であり、図の中段右に示す画像Ai-1は、図の上段に示す過去に与えられたi-1フレーム分の原画像P1~Pi-1の平均的な特徴を有する平均画像である。新たに入力された原画像Piについてのマスク画像Miは、当該原画像Piと平均画像Ai-1とを比較することによって作成される。図2の下段は、このようにして作成されたマスク画像Miを示すものである。具体的には、原画像Pi上の各画素の画素値を、平均画像Ai-1上の対応する画素の画素値と比較し、類似範囲内と判定された画素については背景領域Bを示す画素値たとえば、1を与え、類似範囲外と判定された画素については前景領域Fを示す画素値たとえば、0を与えることにより、二値画像からなるマスク画像Miを作成すればよい。続いて、第i1番目の原画像Pi1が入力されたときには、過去に入力された第i番目の原画像Piまでの平均的な特徴を有する平均画像Aiが新たに作成され、原画像Pi1と平均画像Aiとを比較することにより、マスク画像Mi1が作成される。このように、平均画像は逐次更新されてゆくため、天候や時刻などの要因により、照明環境が変化した場合でも、各時点に適した背景画像としての機能を果たすことができる。なお、平均画像を作成する際には、時系列に沿った重みづけを行うようにし、最新の原画像の情報に重みをおいた加重平均を求めるようにするのが好ましい。たとえば、過去10分間に入力された原画像について重み1、それ以前に入力された原画像について重み0を与えて平均を求めれば、常に、最新10分間に得られた原画像についての平均画像を得ることができ、照明環境が変化した場合にも対応することができる。もちろん、常に最新の原画像ほど重みづけが大きくなるように、重みづけをきめ細かく設定することもできる。結局、この図2に示す例の場合、平均画像を作成する際には、第i番目の原画像Piが入力されたときにi=1,2,、当該原画像Piを含めた過去の原画像について、それぞれ対応する位置の画素の各色別の画素値の重みづけ平均値を算出し、当該平均値をもつ画素の集合体からなる第i番目の平均画像Aiを作成すればよい。そして、第i1番目の原画像Pi1が入力されたときには、第i1番目のマスク画像Mi1を作成するために、原画像Pi1と、第i番目の平均画像Aiとを比較すればよい。ただ、1秒間に30フレームという一般的なフレームレートで動画画像が入力された場合、過去10分間の原画像であっても、その画像データをすべて蓄積しておくには、画像データのバッファにかなりの記憶容量が必要になる。そこで、実用上は、第i番目の原画像Piが与えられたときに、第i-1番目の平均画像Ai-1と、当該第i番目の原画像Piとの2枚の画像に基づいて、第i番目の平均画像Aiを作成する手法が採られる。図3は、このような手法に基づく平均画像の作成方法の一例を示す図である。まず、最初の原画像P1が与えられたときには、当該原画像P1をそのまま最初の平均画像A1とする処理を行う。そして、以後、図3の上段右に示す第i番目の原画像Piが入力されるたびにi=2,3,、図3の上段左に示す第i-1番目の平均画像Ai-1を利用して、図3の中段に示す第i番目の平均画像Aiを、図3の下段に示すai=1-wai-1wpi式1なる演算式を用いて作成すればよい。ここで、aiは、平均画像Aiの所定位置の画素の所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1である。要するに、直前に与えられた原画像Pi-1までの原画像について算出された第i-1番目の平均画像Ai-1の画素値ai-1と、新たに与えられた第i番目の原画像Piの画素値piとについて、重みwを考慮した平均値が算出され、当該算出値を第i番目の平均画像Aiの画素値aiとする処理が行われることになる。このように、新たな原画像Piが入力されるたびに、平均画像Aiを更新してゆき、しかも更新には直前の平均画像Ai-1のみを用いるようにすれば、過去の原画像データを何フレーム分にもわたって蓄積保持しておく必要がないので、処理を行うために必要なバッファの容量を大幅に節約することができる。図3に示すとおり、重みwの値を大きくすればするほど、最新の原画像Piの画素値が平均画像Aiの画素値に大きな影響を与えることになる。図4は、図3に示す方法に基づいて作成される平均画像上の1画素の画素値変遷プロセスを示す図である。この例では、重みw=0.01に設定した具体例が示されている。たとえば、第1番目i=1の原画像P1の特定位置の画素の画素値がp1=100であったとすると、第1番目の平均画像A1の当該特定位置の画素の画素値a1は、そのままa1=100となる。続いて、第2番目の原画像P2の当該特定位置の画素の画素値がp2=100であったとすると、第2番目の平均画像A2の当該特定位置の画素の画素値a2は、a1とp2との重みつき平均0.99×1000.01×100を計算することにより、a2=100になる。図示の例では、原画像の当該特定位置の画素値は、p1~p4まで100のまま変化しないため、平均画像の当該特定位置の画素値も、a1~a4まで100を維持する。ところが、第5番目の原画像P5では、画素値p5=200に急変したため、第5番目の平均画像A5の画素値a5は、図4の下段の式に示すとおり、a4とp5との重みつき平均0.99×1000.01×200を計算することにより、a5=101になる。更に、第6番目の原画像P6および第7番目の原画像P7では、画素値p6=p7=200が維持されているため、図4の下段の式に示すとおり、第6番目の平均画像A6および第7番目の平均画像A7では、画素値a6=102、画素値a7=103、と徐々に増加している。この例のように、重みwを0.01程度の値に設定すると、原画像の画素値が急激に変化しても、平均画像の画素値は直ちに変化することはなく、原画像の画素値にゆっくりと近づいてゆくことになる。したがって、図1に示す例のように、車両が通過する動画が与えられたとしても、一過性の画素値変化は、平均画像の画素値に大きな変化をもたらすことはなく、平均画像は背景画像としての機能を果たすことができる。一方、時刻が日中から夕暮れに変わった場合など、照明環境に変化が生じた場合、当該変化は多数のフレーム画像にわたって持続するため、平均画像の画素値も当該変化に追従して変化することになる。したがって、日中には日中の背景画像に適した平均画像が得られ、夕暮れには夕暮れの背景画像に適した平均画像が得られることになる。移動物体の領域を示すマスク画像Mは、このような方法で作成された平均画像を背景画像として、個々の原画像を比較することによって作成することができる。図5は、このようにして作成されるマスク画像Mの画素値決定プロセスを示す図である。第i番目の原画像Piについてのマスク画像Miは、次のような方法で決定される画素値をもつ画素の集合体からなる二値画像である。すなわち、マスク画像Miの特定位置の画素の画素値miは、原画像Piの当該特定位置の画素の画素値piと、平均画像Ai-1の当該特定位置の画素の画素値ai-1との比較によって決定され、両画素値が類似範囲内であれば、背景領域B内の画素であることを示す画素値たとえば、1が与えられ、両画素値が類似範囲外であれば、前景領域F内の画素であることを示す画素値たとえば、0が与えられる。なお、実際には、原画像は三原色の各画素値たとえば、RGB表色系の場合、R,G,Bの3色の画素値をもったカラー画像であるので、原画像Piを構成する各画素の画素値piや、平均画像Aiを構成する各画素の画素値aiは、いずれも色ごとに独立した3つの値から構成されていることになる。そして、式1の演算は、個々の色ごとにそれぞれ独立して実行されることになる。たとえば、平均画像Ai上の特定位置の画素の画素値aiは、R色の画素値、G色の画素値、B色の画素値という3組の画素値によって構成され、R色の画素値は、過去の原画像のR色の画素値の重みつき平均として得られた値になる。したがって、図5における画素値piと画素値ai-1との比較は、単なる2つの値の比較ではなく、3組の画素値と3組の画素値との比較ということになる。当然ながら、画素の類否判定は、この3組の画素値の組み合わせからなる色同士が類似するか否かを判定する処理ということになる。以下、この類否判定を行うための一般的な方法を、§2において詳述する。<<<§2.画素の一般的な類否判定方法>>>さて、§1で述べた移動物体検出方法では、原画像入力画像と平均画像背景画像とを画素単位で比較して、両画素が類似している場合は背景領域Bの画素、非類似の場合は前景領域Fの画素という判定を行う。このように画素単位での類否判定を行う方法は、領域単位で特徴量を求める方法に比べて処理が単純になり、演算負担は比較的軽くなる。しかしながら、画素の類否判定が移動物体の検出精度を左右するファクターになるため、どのような方法で、どの程度の範囲までを類似範囲と判定するかが非常に重要である。たとえば、§1では、図1の画像P10を平均画像背景画像とし、画像P20を原画像入力画像とした場合、両者の比較により、画像M20のようなマスク画像が得られることを述べた。ここで、マスク画像M20の前景領域Fは、両画像P10,P20を比較した結果、画素が類似範囲外と判定された領域であり、背景画像Bは、画素が類似範囲内と判定された領域である。この場合、類似範囲を適切に設定しないと、正しいマスク画像を得ることはできない。カラー画像の場合、画素の類否は、両画素の色が近似しているか否かに基づいて判定する必要がある。そのためには、一方の画素の3組の画素値と他方の画素の3組の画素値との総合的な比較を行う必要がある。たとえば、RGB表色系の場合、カラー画像を構成する個々の画素には、R,G,Bの3種類の画素値が与えられるが、たとえ、画素値Rが同一であっても、それだけでは両画素の色が類似していることにならない。通常、このようなカラー画像の画素の比較は、三次元の色空間上で、双方の画素値に対応する座標点をプロットし、両座標点の空間的な位置関係を調べることによって行われる。以下、各画素がRGB表色系の画素値をもっている場合について、色空間上での類否判定を行う代表的な方法をいくつか説明しよう。図6は、画素の類否判定の最も単純な方法を示す三次元色空間図である。図示のとおり、三原色R,G,Bの画素値を各座標軸にとったRGB三次元直交座標系に、座標点AおよびPがプロットされている。ここで、座標点ARa,Ga,Baは、背景を示す平均画像A上の画素の画素値Ra,Ga,Baをプロットした点であり、座標点PRp,Gp,Bpは、判定対象となる原画像P上の画素の画素値Rp,Gp,Bpをプロットした点である。ここで、この三次元色空間上での2点A,Pの距離をδとすれば、当該距離δは、両画素の類似度を示す1つの指標になる。δ=0であれば、両画素は全く同一の画素値をもった画素になり、δが大きくなればなるほど、3つの画素値の総合的な類似度は低下する。したがって、予め所定の閾値δthを設定しておき、δ<δthであれば類似範囲内、δ≧δthであれば類似範囲外と判定するようにすれば、一応の判定基準を設定することができる。あるいは、原点Oを起点として、座標点Aへ向かうベクトルVaと座標点Pへ向かうベクトルVpとを定義し、両ベクトルVa,Vpのなす角θを、両画素の類似度を示す指標として用いることもできる。θ=0の場合、2点A,Pは必ずしも同一の点にはならないが、原点Oを通る同一直線上の点になるため、少なくとも三原色R,G,Bの配合割合で決まる両画素の色相は類似したものになる。結局、θが小さければ、色相の類似性は高まることになるので、指標δを用いた判定基準と同様に、予め所定の閾値θthを設定しておき、θ<θthであれば類似範囲内、θ≧θthであれば類似範囲外と判定するようにすれば、一応の判定基準を設定することができる。一方、より正確な判定を行う方法として、三原色の各画素値を各座標軸にとった三次元座標系において、一方の画素値に対応する座標に位置する基準点を中心として所定サイズの基準立体を配置し、他方の画素値に対応する座標点との位置関係を調べ、当該座標点が基準立体の外部に位置すると判定できる場合に非類似との判定を行い、内部に位置すると判定できる場合には類似との判定を行う方法が知られている。たとえば、前掲の非特許文献1には、適切な判定基準を設定するために、基準立体として円柱を用いる類否判定方法が提案されている。図7は、この円柱モデルを用いる方法の基本原理を示す三次元色空間図である。この図7においても、図6と同様に、背景を示す平均画像A上の画素の画素値Ra,Ga,Baが座標点Aとしてプロットされ、判定対象となる原画像P上の画素の画素値Rp,Gp,Bpが座標点Pとしてプロットされている。ここでも、ベクトルVa,Vpは、原点Oを起点として、各点A,Pへ向かうベクトルである。このモデルでは、図示のとおり、座標点Aを中心として所定半径をもった長さLの円柱Cが定義されている。別言すれば、円柱Cは、ベクトルVa上に定義された基準軸Zを中心軸とし、その中心点が座標点Aにくるように配置された円柱ということになる。そして、この円柱Cの内部の空間領域を、背景を示す平均画像Aを構成する画素の類似範囲と定める。したがって、図示の例のように、座標点Pに対応する画素値をもった画素は、円柱Cの外側に位置するため類似範囲外と判定され、前掲領域Fを構成する画素と認定されることになる。逆に、座標点Pが円柱Cの内部に位置していれば、類似範囲内と判定され、背景領域Bを構成する画素と認定される。ここに示す円柱モデルの特徴は、ベクトルVaの方向を向いた基準軸Zを長手方向とする円柱領域が設定される点である。この円柱Cの半径は、色相の類似範囲を決めるパラメータとして機能し、長さLは明度の類似範囲を決めるパラメータとして機能する。実際の景色では、太陽光による照明環境の変化、風による樹木等の明暗の変化などによって、背景領域Bの画素値に時間的変動が生じることになるが、このような変動は、色相よりも明度に大きな影響を及ぼすと考えられる。したがって、このモデルのように、基準軸Zを長手方向とする円柱Cに基づいて類似範囲を設定することは理にかなっている。<<<§3.従来提案されている移動物体検出装置>>>ここでは、§1および§2で述べた原理に基づいて移動物体の検出を行う装置の基本構成を、図8のブロック図を参照しながら説明する。この装置は、特願2009-226146号明細書に開示された装置であり、与えられた動画画像について移動物体を検出する機能をもった移動物体検出装置である。図示のとおり、この装置は、画像入力部110、原画像格納部120、平均画像作成部130、平均画像格納部140、画像比較部200、マスク画像格納部310、画像出力部320、パラメータ設定部330によって構成される。画像入力部110は、時系列で連続的に与えられるフレーム単位のカラー原画像Pを、三原色の各画素値を有する画素の集合体データとして入力する構成要素であり、原画像格納部120は、こうして入力された原画像を逐次格納する構成要素である。一方、平均画像作成部130は、過去に入力された複数の原画像Pに基づいて、これら原画像の平均的な特徴を有する平均画像Aを逐次作成する構成要素であり、平均画像格納部140は、作成された平均画像Aを逐次格納する構成要素である。この平均画像Aを作成する具体的な処理機能については、図3を参照して§1で述べたとおりである。すなわち、平均画像作成部130は、最初の原画像P1が与えられたときには、当該原画像P1をそのまま最初の平均画像A1として平均画像格納部140に格納する処理を行う。そして、以後、第i番目の原画像Piが入力されるたびにi=2,3,、この原画像Piと平均画像格納部140に格納されていた第i-1番目の平均画像Ai-1を利用して、第i番目の平均画像Aiを作成し、これを平均画像格納部140に格納する。このとき、平均画像Aiの所定位置の画素の所定色の画素値をaiとすれば、画素値aiは、既に述べたとおり、ai=1-wai-1wpi式1なる演算式を用いて算出される。ここで、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1である。なお、各画素値は、実際には、§2で述べたとおり、たとえば、RGB表色系の場合、R,G,Bの3種類の画素値から構成され、上記演算は、個々の色ごとの画素値についてそれぞれ別個に行われる。図8に示す原画像格納部120や平均画像格納部140は、バッファメモリによって構成することができる。平均画像作成部130が、図3に示す方法で重みwをパラメータとして平均画像を逐次作成し、これを平均画像格納部140に格納してゆくようにすれば、原画像格納部120には、常に処理に必要な最新の原画像のみが格納されるようにし、平均画像格納部140には、常に処理に必要な最新の平均画像Aのみが格納されるようにすればよいので、バッファメモリの容量を節約することができる。画像比較部200は、原画像格納部120に格納された原画像Pと、平均画像格納部140に格納された平均画像Aとを比較し、前景領域Fと背景領域Bとを区別するマスク画像Mを作成する機能を果たす構成要素であり、図示のとおり、画素値読出部210、類否判定部220、画素値書込部230によって構成されている。画素値読出部210は、比較対象となる一対の画像すなわち、原画像格納部120に格納されている原画像Pと、平均画像格納部140に格納されている平均画像Aのうち、いずれか一方の画像の所定位置の画素の画素値を基準画素値R,G,Bとして読み出し、他方の画像の対応位置の画素の画素値を比較画素値r,g,bとして読み出す処理を行う。類否判定部220は、比較画素値r,g,bが、基準画素値R,G,Bについて設定された所定の類似範囲に入っているか否かを判定する処理を行う。具体的には、図7に示す円柱モデルを利用する場合、三原色の各画素値を各座標軸にとった三次元座標系において、基準画素値R,G,Bに対応する座標に位置する基準点AR,G,Bと比較画素値r,g,bに対応する座標に位置する比較点Pr,g,bをとり、基準点Aを中心とする所定サイズの円柱Cと比較点Pとの位置関係を調べ、比較点Pが円柱Cの外部に位置する場合には非類似との判定を行い、内部に位置する場合には類似との判定を行う。このとき、円柱Cのサイズは、所定の判定パラメータによって決定される。そして、画素値書込部230は、マスク画像Mを構成する所定位置画素値読出部210の読出対象となった位置の画素の画素値として、類否判定部220の判定結果に応じた値を定め、マスク画像格納部310に書き込む処理を行う。すなわち、類否判定部220が類似と判定した場合には背景領域Bを示す画素値たとえば、1を、類否判定部220が非類似と判定した場合には前景領域Fを示す画素値たとえば、0を、それぞれマスク画像格納部310に書き込む処理を行う。マスク画像格納部310は、こうして作成されたマスク画像Mを格納する構成要素であり、画像出力部320は、このマスク画像Mを出力する構成要素である。一方、パラメータ設定部330は、平均画像作成部130による平均画像作成処理に利用されるパラメータw図3に示す重みwと、類否判定部220の類否判定処理に利用される判定パラメータ図7に示す円柱Cの半径と長さLとを、ユーザの操作入力によって任意の値に設定する機能を有する。ユーザは、必要に応じて、これらパラメータの値を調整することにより、より精度の高い検出結果を得ることができる。もっとも、これらのパラメータの一部もしくは全部は、固定値にしておくことも可能である。全部のパラメータ値を固定値にする実施例では、パラメータ設定部330を設ける必要はない。なお、実用上、この図8に示す移動物体検出装置は、コンピュータに専用のプログラムを組み込むことにより構築することができ、図8に各ブロックで示す構成要素は、コンピュータのハードウエアとソフトウエアとの組み合わせによって実現することができる。もちろん、この移動物体検出装置を半導体集積回路として構成することも可能である。この移動物体検出装置による検出精度を向上させるには、各パラメータの値を利用環境に合致した最適値に設定することが重要である。また、類否判定部220において適切な判定アルゴリズムを用いることも重要である。ここでは、図7に示す円柱モデルを利用した例を述べたが、たとえば、§7で述べる回転楕円体モデルを利用すると、一般的に検出精度の向上を図ることができる。このように、図8に示す移動物体検出装置の検出精度を向上させるために、様々な工夫を施すことが可能であるが、本願発明者が行った実験によると、時刻や天候などによる照明変動の影響を完全に排除することは困難であり、依然として、誤検出が生じる問題がある。図1に示す各原画像P10~P40には、背景部分に全く同じ画像が描かれているが、実際の風景では、雲の動きで太陽光による照明環境が変化し、風により樹木の葉の位置や向きが微妙に変化する。したがって、類否判定部220における画素の類否判定基準を厳しくするとたとえば、図7に示す円柱モデルの場合、円柱Cのサイズを小さくすると、実際には背景部分であるのに、非類似との誤判定が生じ、前景領域F内の画素との認定がなされてしまう。逆に、画素の類否判定基準を緩めると、実際には前景部分であるのに、類似との誤判定が生じ、背景領域B内の画素との認定がなされてしまう。もちろん、過去に入力された多数の原画像に対して、より複雑な統計処理を施して背景画像を決定するようにすれば、検出精度の向上を図ることが可能であるが、そのためには演算処理が極めて複雑になり、一般的なプロセッサによるリアルタイム処理を行うことが困難になる。また、多数のフレームにわたって入力画像を格納しておく必要があるため、大容量の画像バッファメモリが必要になる。<<<§4.改良型の移動物体検出装置>>>ここで述べる改良型の移動物体検出装置は、図8に示した従来提案されている移動物体検出装置を改良したものである。前述したとおり、この改良型の移動物体検出装置は、特願2011-020228号明細書に開示されている装置であるが、本発明に係る日照変化領域検出装置に密接に関連した構成を有しているため、説明の便宜上、本発明に係る日照変化領域検出装置の説明を行う前に、この改良型の移動物体検出装置の構成および動作についての説明を行うことにする。図9は、この改良型の移動物体検出装置の基本構成を示すブロック図である。図示のとおり、この装置は、画像入力部110、原画像格納部120、平均画像作成部130、参照画像作成部135、平均画像格納部140、参照画像格納部145、画像比較部400、マスク画像格納部310、画像出力部320、パラメータ設定部340によって構成され、動画画像について移動物体を検出する移動物体検出装置として機能する。この図9に示す各構成要素のうち、画像入力部110、原画像格納部120、平均画像作成部130、平均画像格納部140の機能は、図8に同符号で示す各構成要素の機能と全く同じである。すなわち、画像入力部110は、時系列で連続的に与えられるフレーム単位のカラー原画像を、三原色の各画素値を有する画素の集合体データとして入力する機能を果たし、原画像格納部120は、入力された原画像を逐次格納する機能を果たす。また、平均画像作成部130は、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した平均画像とに基づいて、新たな平均画像を作成する機能を果たし、平均画像格納部140は、作成された平均画像を逐次格納する機能を果たす。たとえば、画像入力部110によって第i番目の原画像Piが入力され、原画像格納部120に格納された場合、平均画像作成部130は、平均画像格納部140から直前に作成した平均画像Ai-1を読み出し、原画像格納部120から新たな原画像Piを読み出し、これらの画像に基づいて、新たな平均画像Aiを作成し、これを平均画像格納部140に格納する。より具体的には、平均画像作成部130は、第1番目の原画像P1をそのまま第1番目の平均画像A1として平均画像格納部140に格納し、第i番目i≧2の原画像Piが入力された場合、この第i番目の原画像Piと第i-1番目の平均画像Ai-1とについて、互いに対応する位置にある画素の三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を各画素値とする画素の集合からなる画像を第i番目i≧2の平均画像Aiとして平均画像格納部140に格納する。一方、参照画像作成部135および参照画像格納部145は、改良型の移動物体検出装置において新たに付加された構成要素である。参照画像作成部135は、画像入力部110によって新たな原画像が入力されるたびに、直前に作成した参照画像に基づいて新たな参照画像を作成する機能を果たし、参照画像格納部145は、作成された参照画像を逐次格納する機能を果たす。これらの具体的な機能については、後に詳述する。また、画像比較部400は、図8に示す画像比較部200と同様に、画像を個々の画素単位で比較し、類否判定を行う機能をもった構成要素であるが、その比較対象には、新たに参照画像が加わることになる。すなわち、画像比較部400は、原画像格納部120に格納された原画像Pと、平均画像格納部140に格納された平均画像Aと、参照画像格納部145に格納された参照画像Uと、を用いた比較処理により、前景領域Fと背景領域Bとを区別するマスク画像Mを作成する機能を果たす。図示のとおり、画像比較部400は、画素値読出部410、前背判定部450、画素値書込部460によって構成されているが、これら各構成要素の機能については後に詳述する。図9に示すマスク画像格納部310および画像出力部320の機能も、図8に同符号で示す各構成要素の機能と全く同じである。すなわち、マスク画像格納部310は、画像比較部400によって作成されたマスク画像Mを格納する機能を果たし、画像出力部320は、当該マスク画像Mを出力する機能を果たす。パラメータ設定部340は、この装置で利用されるパラメータを設定する構成要素である。この装置の場合、平均画像作成部130で用いられる重みw、参照画像作成部135で用いられる重みv、前背判定部450で用いられる判定パラメータたとえば、図7に示す円柱モデルを利用した類否判定を行う場合は、円柱Cの半径と長さLを、パラメータ設定部340によって設定することができる。ユーザは、必要に応じて、パラメータ設定部340に対して操作入力を行い、これらパラメータの値を調整し、より精度の高い検出結果を得ることができる。もちろん、これらのパラメータの一部もしくは全部は、固定値にしておくことも可能であり、全部のパラメータ値を固定値にする実施例では、パラメータ設定部340を設ける必要はない。図9に示す改良型の移動物体検出装置が、図8に示す従来の移動物体検出装置と大きく異なる点は、参照画像作成部135および参照画像格納部145を新たに付加した点と、画像比較部200の代わりに画像比較部400を設けた点である。そこで、まず、画像比較部400の処理機能について説明する。図示のとおり、画像比較部400は、画素値読出部410、前背判定部450、画素値書込部460によって構成されており、前背判定部450は、第1判定部420、第2判定部430、総合判定部440によって構成されている。画素値読出部410は、画像入力部110が第i番目の原画像Piを入力し、これを原画像格納部120に格納するたびに、この原画像格納部120に格納されている第i番目の原画像Pi、平均画像格納部140に格納されている第i-1番目の平均画像Ai-1、参照画像格納部145に格納されている第i-1番目の参照画像Ui-1から、それぞれ対応する所定位置にある画素の画素値を読み出す処理を行う。ここでは、対応する所定位置にある各画素の画素値をそれぞれpi,ai-1,ui-1で表すことにする。原画像P,平均画像A,参照画像Uは、いずれも同じ二次元画素配列をもった画像であり、対応する所定位置にある各画素とは、この二次元画素配列において同一座標x,yに位置する画素ということになる。こうして画素値読出部410が読み出した、同一座標x,yに位置する画素の画素値pi,ai-1,ui-1は、前背判定部450に与えられる。なお、前述したとおり、実際には、各画素は、三原色R,G,Bの成分を有しているので、画素値pi,ai-1,ui-1は、それぞれ三原色R,G,Bの成分ごとに独立した値をもつことになる。たとえば、画素値piは、実際には、R色成分pir,G色成分pig,B色成分pibという3通りの画素値の集合体であるが、ここでは便宜上、単に画素値piと呼ぶことにする。前背判定部450は、画素値読出部410によって読み出された画素値pi,ai-1,ui-1に基づいて、第i番目の原画像Piの所定位置座標x,yの位置にある画素が前景画素か背景画素かを判定する処理を行う。ここでは、判定結果をJなる論理値で表す。具体的には、判定結果がJ=0であれば前景画素、J=1であれば背景画素であることを示す。前背判定部450は、第1判定部420、第2判定部430、総合判定部440によって構成されており、座標x,yに位置する画素についての判定を次のようなプロセスで実行する。まず、第1判定部420は、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の平均画像Ai-1を構成する画素の画素値ai-1とを比較して類似するか否かを判定する。一方、第2判定部430は、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の参照画像Ui-1を構成する画素の画素値ui-1とを比較して類似するか否かを判定する。そして、総合判定部440は、第1判定部420による判定結果および第2判定部430による判定結果の少なくとも一方が類似を示す場合には背景画素との判定を行い判定結果J=1を出力し、双方とも非類似を示す場合には前景画素との判定を行い判定結果J=0を出力する。第1判定部420および第2判定部430は、2つの画素間の類否判定を行うことになるが、その具体的な方法は§2で例示したとおりである。たとえば、図7に示す円柱モデルを利用する場合、第1判定部420は、次のような方法で、画素値piと画素値ai-1との類否判定を行うことができる。まず、画素値ai-1を構成する三原色の各画素値に基づいて、RGB三次元座標系に基準点Aをプロットし、この基準点Aを中心とする所定サイズの基準立体として円柱Cを配置する。次に、画素値piを構成する三原色の各画素値に基づいて、RGB三次元座標系に比較点Pをプロットする。そして、円柱Cと比較点Pとの位置関係を調べ、比較点Pが円柱Cの外部に位置する場合には非類似との判定を行い、内部に位置する場合には類似との判定を行えばよい。第2判定部430も、同様の方法で、画素値piと画素値ui-1との類否判定を行うことができる。総合判定部440は、これら2通りの判定結果に基づき、少なくとも一方が類似を示す場合には判定結果J=1背景画素との判定結果を出力し、双方とも非類似を示す場合には判定結果J=0前景画素との判定結果を出力する。このような判定結果の出力は、第1判定部420および第2判定部430の判定結果の論理和として得ることができる。すなわち、第1判定部420および第2判定部430が、類似との判定結果を論理値1として出力し、非類似との判定結果を論理値0として出力するようにしておけば、総合判定部440は、これらの各判定結果の論理和として、最終的な判定結果Jを出力することができる。画素値書込部460は、第i番目のマスク画像Miを構成する所定位置座標x,yの位置の画素の画素値として、前背判定部450が背景画素と判定した場合には背景領域を示す画素値この例では、判定結果Jと同じ論理値1を、前背判定部450が前景画素と判定した場合には前景領域を示す画素値この例では、判定結果Jと同じ論理値0を、それぞれマスク画像格納部310に書き込む処理を行う。以上、画像比較部400が、特定の座標x,yに位置する1画素について実行する具体的な処理を述べたが、このような処理が、1枚の画像を構成する全画素について繰り返し実行される。こうして、原画像格納部120に格納された第i番目の原画像Piを構成する全画素について、それぞれ背景画素か前景画素かの判定が行われることになり、その判定結果に基づいて、マスク画像格納部310内に第i番目i≧2のマスク画像Miが得られることになる。こうして得られたマスク画像Miにおいて、前景画素からなる領域上例の場合、画素値0をもった領域は、前景領域、すなわち、移動物体が存在すると判断できる領域ということになる。続いて、参照画像作成部135および参照画像格納部145の処理機能について説明する。参照画像作成部135および参照画像格納部145の役割は、平均画像作成部130および平均画像格納部140の役割と対比して考えることができる。後者は、平均画像Aを作成しこれを逐次格納する役割を果たすのに対して、前者は、参照画像Uを作成しこれを逐次格納する役割を果たす。そして、平均画像Aも参照画像Uも、画像比較部400において、原画像Pと比較されることになる。この比較処理では、平均画像Aも参照画像Uも、原画像Pに対する背景としての意味合いをもつ。ただ、平均画像Aと参照画像Uとは、その作成プロセスが異なるため、原画像に動的変化が生じた場合、その追従性に差が生じることになる。平均画像作成部130による平均画像Aの作成処理は、既に§3で述べたとおりである。すなわち、平均画像作成部130は、まず、初期設定として、第1番目の原画像P1をそのまま第1番目の平均画像A1として、平均画像格納部140に格納する処理を行う。続いて、第2番目以降の原画像が入力されると、第i番目の原画像Piと第i-1番目の平均画像Ai-1とについて、互いに対応する位置にある画素の三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を各画素値とする画素の集合からなる画像を第i番目i≧2の平均画像Aiとして、平均画像格納部140に格納する処理を行う。図10は、図9に示す移動物体検出装置の各構成要素で実行される演算内容を示す図である。図の上段には、平均画像作成部130で実行される演算内容が式1Aとして示されている。この式1Aは、図3に示す式1と同じものである。要するに、平均画像作成部130は、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の平均画像A1として平均画像格納部140へ格納する処理を行い図10上段の但し書、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の平均画像Aiを、ai=1-wai-1wpi式1A但し、aiは、平均画像Aiの所定位置の画素の所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1なる演算式を用いて作成することになる。これに対して、参照画像作成部135による参照画像Uの作成処理は、次のようなプロセスによって行われる。すなわち、参照画像作成部135は、まず、初期設定として、第1番目の原画像P1をそのまま第1番目の参照画像U1として、参照画像格納部145に格納する処理を行う。続いて、第2番目以降の原画像が入力されると、直前の参照画像に基づいて、次のような方法で新たな参照画像を作成する処理を行う。すなわち、第i番目i≧2の原画像Piが入力され、この原画像Piを構成する個々の画素について、前背判定部450において前景画素か背景画素かの判定が行われたら、参照画像作成部135は、参照画像格納部145に格納されている第i-1番目の参照画像Ui-1を構成する個々の画素のうち、前背判定部450によって第i番目の原画像Piの対応画素が前景画素と判定された画素については、更新せずに、もとの画素値ui-1をそのまま維持し、前背判定部450によって第i番目の原画像Piの対応画素が背景画素と判定された画素については、当該画素と第i-1番目の平均画像Ai-1の対応する位置にある画素とについて、三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を新画素値として更新することにより、第i番目i≧2の参照画像Uiを作成する処理を行う。図10の中段には、参照画像作成部135で実行される演算内容が示されている。すなわち、参照画像作成部135は、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の参照画像U1として参照画像格納部へ格納する処理を行い図10中段の最下行、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の参照画像Uiの各画素値を、次のような方法で決定する処理を行う。この決定処理では、前背判定部450による判定結果が利用される。すなわち、前背判定部450によって背景画素と判定された画素判定結果J=1が得られた画素については、ui=1-vui-1vai-1式1B但し、uiは、平均画像Uiの所定位置の画素の所定色の画素値、ui-1は、平均画像Ui-1の前記所定位置の画素の前記所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、vは、所定の重みを示すパラメータv<1なる演算式を用いて求められた画素値uiに更新し、前景画素と判定された画素判定結果J=0が得られた画素については、ui=ui-1式1Cとして、更新せずに、もとの画素値をそのまま維持する。図10の下段には、前背判定部450で実行される判定処理の概念が表として示されている。前述したとおり、前背判定部450における判定処理は、第1判定部420で行われる画素値piとai-1との類否判定処理表の横欄と、第2判定部430で行われる画素値piとui-1との類否判定処理表の縦欄と、これらの判定結果に基づいて総合判定部440で行われる総合判定処理と、によって構成される。図10下段の表の横欄は、第1判定部420による判定結果の別類似○/非類似×を示し、縦欄は、第2判定部430による判定結果の別類似○/非類似×を示している。総合判定部440は、少なくとも一方が類似○を示す場合には背景画素J=1、双方とも非類似×を示す場合には前景画素J=0との判定を行うことになる。参照画像作成部135は、この判定結果Jを利用して、図10の中段に示す処理を実行する。なお、実用上、この図9に示す改良型の移動物体検出装置も、コンピュータに専用のプログラムを組み込むことにより構築することができ、図9に各ブロックで示す構成要素は、コンピュータのハードウエアとソフトウエアとの組み合わせによって実現することができる。もちろん、この移動物体検出装置を半導体集積回路として構成することも可能である。<<<§5.画像変遷の実例>>>図11および図12は、図9に示す改良型の移動物体検出装置における各画像の変遷の一例を示す図である。前述したとおり、この装置を監視システムなどに利用する場合、実際には、原画像Pは1秒間に30フレーム程度の周期で入力されることになるが、ここでは、便宜上、平均画像Aと参照画像Uとの変遷態様の相違を説明するのに都合が良い非現実的なモデルを用いた変遷例を示すことにする。図11および図12に股がって示されている縦方向軸は、時間軸に対応しており、t1~t15は時間軸上の各時点を示す。また、横方向に並べた3つの画像は、それぞれ各時点における原画像P,平均画像A,参照画像Uを示している。<第1の時点t1>まず、第1の時点t1において、図示のような原画像P1が与えられたものとしよう。ここでは、便宜上、原画像P1を真っ白な画像として示しているが、この真っ白な画像は、街頭に設置された監視カメラから見た背景を示す画像であるものとする。このような第1番目の原画像P1が与えられると、平均画像作成部130は、当該原画像P1をそのまま第1番目の平均画像A1として平均画像格納部140に格納する。同様に、参照画像作成部135は、当該原画像P1をそのまま第1番目の参照画像U1として参照画像格納部145に格納する。したがって、図11の最上段に示すとおり、第1の時点t1では、図示のような画像P1,A1,U1が得られることになる。なお、原画像P1,平均画像A1,参照画像U1内に示された小さな正方形は、それぞれ同一位置x,yの画素p1,a1,u1を示している。以下、個々の画像について、この特定位置x,yに配置された特定の画素を代表画素として、この代表画素についてのみ、その画素値の変遷を追ってゆくことにする。ここで、各画素p1,a1,u1の画素値も画素と同じ符号で示すことにすれば、時点t1では、当然ながら、p1=a1=u1である。<第2の時点t2>続いて、第2の時点t2において、図示のような原画像P2が与えられたものとしよう。ここで、原画像P2は原画像P1と全く同じ画像であり、監視カメラから見た背景に何ら変わりがなかった例を示している。したがって、画素値p2=p1である。次に、特定位置x,yに配置された画素についての前背判定部450による判定処理を考えると、第1判定部420では、画素p2と、1つ前の時点の画素a1との類否判定が行われ、第2判定部430では、画素p2と、1つ前の時点の画素u1との類否判定が行われる。図に実線で示す矢印は、この類否判定の相手先を示しており、矢印上の記号○/×は判定結果○は類似/×は非類似を示している。図示のとおり、p2=a1であり、p2=u1であるから、いずれも類似○との判定結果が得られ、前背判定部450は、J=1なる判定結果を出す。これは、画素p2が背景画素と判定されたことを意味する。さて、平均画像作成部130は、このような判定結果に依存することなしに、常に、図10に示す式1Aに従って、平均画像Aの画素値aを更新する処理を行う。すなわち、時点t2における平均画像A2の特定位置x,yに配置されたの画素の画素値a2は、式1Aに従って、a2=1-wa1wp2なる演算によって求められ、画素値a1はa2に更新される。図の画素a1から画素a2に向かう破線および画素p2から画素a2に向かう破線は、式1Aに従って、画素値a1および画素値p2に基づいて画素値a2が決定されることを示している。もっとも、実際には、a1=p2であるから、a2=a1となり、画素値自身に変動は生じない。これに対して、参照画像作成部135は、図10に示すとおり、判定結果Jに依存して異なる処理を行うことになる。ここに示す例の場合、時点t2についての判定結果は、上述したとおり、J=1になるので、参照画像作成部135は、図10に示す式1Bに従って、参照画像Uの画素値uを更新する処理を行う。すなわち、時点t2における参照画像U2の特定位置x,yに配置されたの画素の画素値u2は、式1Bに従って、u2=1-vu1va1なる演算によって求められ、画素値u1はu2に更新される。図の画素u1から画素u2に向かう破線および画素a1から画素u2に向かう破線は、式1Bに従って、画素値a1および画素値u1に基づいて画素値u2が決定されることを示している。もっとも、実際には、a1=u1であるから、u2=u1となり、画素値自身に変動は生じない。<第3の時点t3>続いて、第3の時点t3において、図示のような原画像P3が与えられたものとしよう。ここで、原画像P3は、視界に車両Vが侵入してきたため、背景上に車両V移動物体が重なった状態の画像になっている。したがって、画素値p3≠p2である。図では、画素p3を含む車両Vの部分を二重斜線ハッチング領域として示してある。この時点t3における特定位置x,yに配置された画素についての前背判定部450による判定処理を考えると、第1判定部420では、画素p3と、1つ前の時点の画素a2との類否判定が行われ、第2判定部430では、画素p3と、1つ前の時点の画素u2との類否判定が行われる。ここでは、説明の便宜上、車両Vの部分図の二重斜線ハッチング部分の色が背景部分図の白地部分と大きく異なっているものとしよう。そうすると、p3≠a2であり、p3≠u2であるから、いずれも非類似×との判定結果が得られ、前背判定部450は、J=0なる判定結果を出す。これは、画素p3が前景画素と判定されたことを意味する。したがって、第3のマスク画像M3において、特定位置x,yの画素は前景領域移動物体の領域に含まれる。さて、平均画像作成部130は、このような判定結果に依存することなしに、常に、図10に示す式1Aに従って、平均画像Aの画素値aを更新する処理を行う。すなわち、時点t3における平均画像A3の特定位置x,yに配置されたの画素の画素値a3は、式1Aに従って、a3=1-wa2wp3なる演算によって求められ、画素値a2はa3に更新される。ここで、画素値a3は、画素値p3の影響を受けるため、画素値a2から変動することになる。図では、このようにして画素値が変動した画素a3に中濃ハッチングを施して示すことにする。一方、参照画像作成部135は、図10に示すとおり、判定結果J=0の場合は、値を更新しない。したがって、1つ前の時点の画素値u2をそのまま画素値u3として維持することになる。図の画素u2から画素u3に向かう破線は、更新が行われずに、元の画素値がそのまま維持されることを示している。<第4の時点t4>続いて、第4の時点t4において、図示のような原画像P4が与えられたものとしよう。この原画像P4では、車両Vが視界から消え、元の原画像P1,P2と同じ状態に戻っている。したがって、画素値p4≠p3である。この時点t4における前背判定部450による判定処理を考えると、第1判定部420では、画素p4と、1つ前の時点の画素a3との類否判定が行われ、第2判定部430では、画素p4と、1つ前の時点の画素u3との類否判定が行われる。ここでは、説明の便宜上、図の中濃ハッチング部分で示す色が、白地部分で示す色の類似範囲内にギリギリ入っていたものとしよう。実際、重みwの値を0.01程度に設定すると、車両Vが視界を通り過ぎた程度では、平均画像Aの画素値変動はそれほど大きくなく、画素値p4とa3とが類似範囲に入るケースは現実的である。また、画素値u3は更新されずに画素値u2を維持しているので、当然、画素値p4=u3である。よって、いずれの判定においても類似○との判定結果が得られ、前背判定部450は、J=1なる判定結果を出す。これは、画素p4が背景画素と判定されたことを意味する。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a4=1-wa3wp4なる演算によって求められ、画素値a3はa4に更新される。ここで、画素値a4は、画素値p4の影響を受けるため、画素値a3から変動する。図では、このようにして画素値が変動した画素a4を白地の画素として示してある。これに対して、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u4=1-vu3va3なる演算によって求められ、画素値u3はu4に更新される。図では、このようにして画素値が変動した画素u4に薄いハッチングを施して示すことにする。<第5の時点t5>続いて、第5の時点t5において、図示のような原画像P5が与えられたものとしよう。この原画像P5では、再び視界に車両Vが侵入してきたため、背景上に車両V移動物体が重なった状態の画像になっている。したがって、画素値p5≠p4である。図では、画素p5を含む車両Vの部分を二重斜線ハッチング領域として示してある。この時点t5における前背判定部450による判定処理を考えると、第1判定部420では、画素p5と、1つ前の時点の画素a4との類否判定が行われ、第2判定部430では、画素p5と、1つ前の時点の画素u4との類否判定が行われる。図示の例では、いずれも非類似×との判定結果が得られ、前背判定部450は、J=0なる判定結果を出している。これは、画素p5が前景画素と判定されたことを意味する。したがって、第5のマスク画像M5において、特定位置x,yの画素は前景領域移動物体の領域に含まれる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a5=1-wa4wp5なる演算を行い、画素値a4はa5に更新される。図では、このようにして画素値が変動した画素a5に中濃ハッチングを施して示してある。一方、参照画像作成部135は、判定結果J=0を受けて、値を更新しない。したがって、1つ前の時点の画素値u4をそのまま画素値u5として維持することになる。<第6の時点t6>続いて、第6の時点t6において、図示のような原画像P6が与えられたものとしよう。ここでは、時点t5で侵入してきた車両Vが、そのままの位置で駐車してしまった場合を考えてみよう。この場合、原画像P6は原画像P5と全く同じになる。この時点t6における前背判定部450による判定処理を考えると、第1判定部420では、画素p6と、1つ前の時点の画素a5との類否判定が行われ、第2判定部430では、画素p6と、1つ前の時点の画素u5との類否判定が行われる。図示の例では、いずれも非類似×との判定結果が得られ、前背判定部450は、J=0なる判定結果を出している。これは、画素p6が前景画素と判定されたことを意味する。したがって、第6のマスク画像M6において、特定位置x,yの画素は前景領域移動物体の領域に含まれる。一般に、重みw,vを小さく設定しておけば、原画像Pに対する平均画像Aや参照画像Uの追随速度は小さくなるので、車両Vが駐車していた場合でも、当該車両Vが背景に溶け込み、背景の一部として認識されるまでには時間がかかることになる。したがって、図示の例では、車両Vは時点t6ではまだ移動物体として認識されることになる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a6=1-wa5wp6なる演算を行い、画素値a5はa6に更新される。図では、このようにして画素値が変動した画素a5に斜線ハッチングを施して示してある。一方、参照画像作成部135は、判定結果J=0を受けて、値を更新しない。したがって、1つ前の時点の画素値u5をそのまま画素値u6として維持することになる。結果的に、時刻t4の画素値u4が時刻t6においても画素値u6として維持されることになる。<第7の時点t7>続いて、第7の時点t7において、図示のような原画像P7が与えられたものとしよう。この例では、時点t5で侵入してきた車両Vが駐車状態にあり、原画像P7は原画像P6と全く同じ状態である。この時点t7における前背判定部450による判定処理を考えると、第1判定部420では、画素p7と、1つ前の時点の画素a6との類否判定が行われ、第2判定部430では、画素p7と、1つ前の時点の画素u6との類否判定が行われる。ここでは、第1判定部420では類似○との判定結果が得られ、第2判定部430では非類似×との判定結果が得られたものとしよう。平均画像Aは新たな原画像Pにより常に更新されるのに対して、参照画像Uは条件付き更新が行われるので、原画像Pの変化に対する追随速度は両者で異なる。このため、2つの判定部による判定結果に差が生じることになる。図示の例では、画素a6の色は画素p7の色にかなり近くなってきているため類似○との判定結果が得られたが、画素u6の色は、まだ画素p7の色に対して差があるため非類似×との判定結果が得られている。図10の前背判定部450の判定表に示すとおり、一方が類似○、他方が非類似×との判定結果が得られた場合、前背判定部450は、当該画素を背景画素と判定し、J=1なる判定結果を出力する。したがって、第7のマスク画像M7において、特定位置x,yの画素は背景領域に含まれる。すなわち、時点t7において、駐車中の車両Vは背景に溶け込んだものと認識されることになる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a7=1-wa6wp7なる演算を行い、画素値a6はa7に更新される。図では、このようにして画素値が変動した画素a6に斜線ハッチングを施して示してある実際には、画素値a6とa7とは若干異なるが、図では便宜上、同一ハッチングを施して示す。一方、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u7=1-vu6va6なる演算を行い、画素値u6をu7に更新する。図では、このようにして画素値が変動した画素u7に中濃ハッチングを施して示してある。<第8の時点t8>続いて、第8の時点t8において、図示のような原画像P8が与えられたものとしよう。この例では、時点t5で侵入してきた車両Vが駐車状態にあり、原画像P8は原画像P7と全く同じ状態である。この時点t8における前背判定部450による判定処理を考えると、第1判定部420では、画素p8と、1つ前の時点の画素a7との類否判定が行われ、第2判定部430では、画素p8と、1つ前の時点の画素u7との類否判定が行われる。ここでは、第1判定部420では類似○との判定結果が得られ、第2判定部430では非類似×との判定結果が得られたものとしよう。図示の例では、画素a7の色は画素p8の色にかなり近くなってきているため類似○との判定結果が得られたが、画素u7の色は、まだ画素p8の色に対して差があるため非類似×との判定結果が得られている。その結果、前背判定部450は、当該画素を背景画素と判定し、J=1なる判定結果を出力する。したがって、第8のマスク画像M8において、特定位置x,yの画素は背景領域に含まれる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a8=1-wa7wp8なる演算を行い、画素値a7はa8に更新される。図では、このようにして画素値が変動した画素a8に斜線ハッチングを施して示してある実際には、画素値a7とa8とは若干異なるが、図では便宜上、同一ハッチングを施して示す。一方、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u8=1-vu7va7なる演算を行い、画素値u7をu8に更新する。図では、このようにして画素値が変動した画素u8に斜線ハッチングを施して示してある。<第9の時点t9>続いて、図12に示す第9の時点t9について説明する図12の最上段に示す第8の時点t8の各画像は、図11の最下段に示す第8の時点t8の各画像と同じものである。ここでは、第9の時点t9において、図示のような原画像P9が与えられたものとしよう。この例では、時点t5で侵入してきた車両Vが依然として駐車状態にあり、原画像P9は原画像P8と全く同じ状態である。この時点t9における前背判定部450による判定処理を考えると、第1判定部420では、画素p9と、1つ前の時点の画素a8との類否判定が行われ、第2判定部430では、画素p9と、1つ前の時点の画素u8との類否判定が行われる。ここでは、いずれの判定結果も類似○との判定結果が得られたものとしよう。すなわち、車両Vが長時間駐車していたため、平均画像Aにおいても参照画像Uにおいても、この車両Vの領域が背景に溶け込んだ状態になり、画素p9は、いずれの判定部においても背景画素と判定されるに至ったことになる。したがって、第9のマスク画像M9において、特定位置x,yの画素は背景領域に含まれる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a9=1-wa8wp9なる演算を行い、画素値a8はa9に更新される。図では、このようにして画素値が変動した画素a9に斜線ハッチングを施して示してある実際には、画素値a8とa9とは若干異なるが、図では便宜上、同一ハッチングを施して示す。一方、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u9=1-vu8va8なる演算を行い、画素値u8をu9に更新する。図では、このようにして画素値が変動した画素u9に斜線ハッチングを施して示してある実際には、画素値u8とu9とは若干異なるが、図では便宜上、同一ハッチングを施して示す。<第10の時点t10>続いて、第10の時点t10において、図示のような原画像P10が与えられたものとしよう。この原画像P10では、これまで駐車していた車両Vが視界から消え、第1番目の原画像P1と同じ状態に戻っている。したがって、画素値p10≠p9である。この時点t10における前背判定部450による判定処理を考えると、第1判定部420では、画素p10と、1つ前の時点の画素a9との類否判定が行われ、第2判定部430では、画素p10と、1つ前の時点の画素u9との類否判定が行われる。画素a9および画素u9の色は、いずれも車両Vの色にかなり近づいていたため、図示の例では、いずれも非類似×との判定結果が得られ、前背判定部450は、J=0なる判定結果を出している。これは、画素p10が前景画素と判定されたことを意味する。したがって、第10のマスク画像M10において、特定位置x,yの画素は前景領域移動物体の領域に含まれる。結局、この例の場合、第10の原画像P10上には車両Vは存在していないが、第10のマスク画像M10上では、車両Vがかつて存在していた領域が前景領域Fとして検出されることになる。もちろん、現実的な監視カメラの撮影画像の場合は、時点t9から時点t10に至るときに車両Vが忽然と消えることはなく、車両Vは画面上を徐々に移動して去ってゆくことになるので、実際には、徐々に移動中の車両Vの占有領域がマスク画像M上に前景領域Fとして現れることになる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a10=1-wa9wp10なる演算を行い、画素値a9はa10に更新される。図では、このようにして画素値が変動した画素a10に中濃ハッチングを施して示してある。一方、参照画像作成部135は、判定結果J=0を受けて、値を更新しない。したがって、1つ前の時点の画素値u9をそのまま画素値u10として維持することになる。<第11の時点t11>続いて、第11の時点t11において、図示のような原画像P11が与えられたものとしよう。この例では、原画像P11は、視界に車両Vが再び侵入してきたため、背景上に車両V移動物体が重なった状態の画像になっている。したがって、画素値p11≠p10である。この時点t11における特定位置x,yに配置された画素についての前背判定部450による判定処理を考えると、第1判定部420では、画素p11と、1つ前の時点の画素a10との類否判定が行われ、第2判定部430では、画素p11と、1つ前の時点の画素u10との類否判定が行われる。ここでは、第1判定部420では非類似×との判定結果が得られ、第2判定部430では類似○との判定結果が得られたものとしよう。図示の例では、画素a10の色は画素p10の影響を受けて背景の色に近づいていたため、車両の一部を構成する画素p11に対しては非類似×との判定結果が得られたが、画素u10は更新されなかったため、画素p11に対して類似○との判定結果が得られたことになる。図10の前背判定部450の判定表に示すとおり、一方が類似○、他方が非類似×との判定結果が得られた場合、前背判定部450は、当該画素を背景画素と判定し、J=1なる判定結果を出力する。したがって、第11のマスク画像M11において、特定位置x,yの画素は背景領域に含まれる。すなわち、ここに示す非現実的なモデルの場合、時点t11に出現した車両Vは、マスク画像M11上では移動物体として検出されない。ただ、後述するように、マスク画像M12上で検出されることになる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a11=1-wa10wp11なる演算を行い、画素値a10はa11に更新される。図では、このようにして画素値が変動した画素a11に斜線ハッチングを施して示してある。一方、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u11=1-vu10va10なる演算を行い、画素値u10をu11に更新する。図では、このようにして画素値が変動した画素u11に斜線ハッチングを施して示してある実際には、画素値u10とu11とは若干異なるが、図では便宜上、同一ハッチングを施して示す。<第12の時点t12>続いて、第12の時点t12について説明する。ここでは、第12の時点t12において、図示のような原画像P12が与えられたものとしよう。この例では、車両Vが視界から消え、第1番目の原画像P1と同じ状態に戻っている。したがって、画素値p12≠p11である。この時点t12における前背判定部450による判定処理を考えると、第1判定部420では、画素p12と、1つ前の時点の画素a11との類否判定が行われ、第2判定部430では、画素p12と、1つ前の時点の画素u11との類否判定が行われる。画素a11および画素u11の色は、いずれも車両Vの色にかなり近いため、図示の例では、いずれも非類似×との判定結果が得られ、前背判定部450は、J=0なる判定結果を出している。これは、画素p12が前景画素と判定されたことを意味する。したがって、第12のマスク画像M12において、特定位置x,yの画素は前景領域移動物体の領域に含まれる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a12=1-wa11wp12なる演算を行い、画素値a11はa12に更新される。図では、このようにして画素値が変動した画素a12に中濃ハッチングを施して示してある。一方、参照画像作成部135は、判定結果J=0を受けて、値を更新しない。したがって、1つ前の時点の画素値u11をそのまま画素値u12として維持することになる。<第13の時点t13>続いて、第13の時点t13において、図示のような原画像P13が与えられたものとしよう。この例では、画像内に新たな車両の侵入はなく、原画像P13は第1番目の原画像P1と全く同じ状態である。この時点t13における前背判定部450による判定処理を考えると、第1判定部420では、画素p13と、1つ前の時点の画素a12との類否判定が行われ、第2判定部430では、画素p13と、1つ前の時点の画素u12との類否判定が行われる。ここでは、第1判定部420では類似○との判定結果が得られ、第2判定部430では非類似×との判定結果が得られたものとしよう。前背判定部450は、この判定結果に基づいて当該画素を背景画素と判定し、J=1なる判定結果を出力する。したがって、第13のマスク画像M13において、特定位置x,yの画素は背景領域に含まれる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a13=1-wa12wp13なる演算を行い、画素値a12はa13に更新される。図では、このようにして画素値が変動した画素a13に薄いハッチングを施して示してある。一方、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u13=1-vu12va12なる演算を行い、画素値u12をu13に更新する。図では、このようにして画素値が変動した画素u13に中濃ハッチングを施して示してある。<第14の時点t14>続いて、第14の時点t14について説明する。ここでは、第14の時点t14において、図示のような原画像P14が与えられたものとしよう。この例では、依然として画像内に車両の侵入はなく、原画像P14は第1番目の原画像P1と全く同じ状態である。この時点t14における前背判定部450による判定処理を考えると、第1判定部420では、画素p14と、1つ前の時点の画素a13との類否判定が行われ、第2判定部430では、画素p14と、1つ前の時点の画素u13との類否判定が行われる。ここでは、いずれも類似○との判定結果が得られたものとしよう。前背判定部450は、この判定結果に基づいて当該画素を背景画素と判定し、J=1なる判定結果を出力する。したがって、第14のマスク画像M14において、特定位置x,yの画素は背景領域に含まれる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a14=1-wa13wp14なる演算を行い、画素値a13はa14に更新される。図では、このようにして画素値が変動した画素a14を白地で示してある。一方、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u14=1-vu13va13なる演算を行い、画素値u13をu14に更新する。図では、このようにして画素値が変動した画素u9に薄いハッチングを施して示してある。<第15の時点t15>最後に、第15の時点t15について説明する。ここでは、第15の時点t15において、図示のような原画像P15が与えられたものとしよう。この例では、依然として画像内に車両の侵入はなく、原画像P15は第1番目の原画像P1と全く同じ状態である。この時点t15における前背判定部450による判定処理を考えると、第1判定部420では、画素p15と、1つ前の時点の画素a14との類否判定が行われ、第2判定部430では、画素p15と、1つ前の時点の画素u14との類否判定が行われる。ここでは、いずれも類似○との判定結果が得られたものとしよう。前背判定部450は、この判定結果に基づいて当該画素を背景画素と判定し、J=1なる判定結果を出力する。したがって、第15のマスク画像M15において、特定位置x,yの画素は背景領域に含まれる。平均画像作成部130は、このような判定結果に依存することなしに、式1Aに従って、a15=1-wa14wp15なる演算を行い、画素値a14はa15に更新される。一方、参照画像作成部135は、判定結果がJ=1になるので、式1Bに従って、u15=1-vu14va14なる演算を行い、画素値u14をu15に更新する。<<<§6.平均画像と参照画像との双方を用いる利点>>>これまで、§3において、従来提案されている移動物体検出装置の構成図8を説明し、§4において、改良型の移動物体検出装置の構成図9を説明した。いずれの装置を用いた場合も、新たに入力された原画像Pについて、背景領域Bと前景領域Fとを区別するマスク画像Mを得ることができ、前景領域Fは、移動物体が占めると推定される領域ということになる。ただ、図8に示す従来装置の場合、マスク画像Mは、画像比較部200における第i番目の原画像Piと第i-1番目の平均画像Ai-1との比較結果に基づいて得られるのに対して、図9に示す改良型の移動物体検出装置の場合、マスク画像Mは、画像比較部400における第i番目の原画像Piと第i-1番目の平均画像Ai-1との比較結果第1判定部420の判定結果と、第i番目の原画像Piと第i-1番目の参照画像Ui-1との比較結果第2判定部430の判定結果と、の双方に基づいて得られることになる。別言すれば、図8に示す従来装置の場合、原画像Pと平均画像Aとの比較により移動物体の検出を行っていたのに対し、図9に示す改良型の移動物体検出装置の場合、原画像Pを平均画像Aと参照画像Uとの双方と比較することにより移動物体の検出を行うことになる。そこで、ここでは、平均画像Aと参照画像Uとの双方を用いた移動物体検出の利点を説明する。平均画像Aと参照画像Uとの大きな相違点は、§5に画像変遷の実例を示したように、平均画像Aの各画素は、新たな原画像Pが入力されるたびに毎回必ず更新されるのに対して、参照画像Uの各画素は、必ずしも毎回更新されるわけではない、という点である。なお、ここで言う更新とは、必ずしも画素値の変動を指すものではなく、前回の画素値をそのまま維持するのではなく、新たな画素値を与えることを意味する。したがって、たとえば、図10上段に示す式1Aにおいて、pi=ai-1であった場合、結果的にai=ai-1になるが、第i番目の平均画像Aiの画素値として新たな画素値aiが与えられているので、値が更新されていることになる。同様に、図10中段に示す式1Bにおいて、ai-1=ui-1であった場合、結果的にui=ui-1になるが、第i番目の参照画像Uiの画素値として新たな画素値uiが与えられているので、値が更新されていることになる。これに対して、図10中段の式1Cui=ui-1は、前回の画素値をそのまま維持することを意味しており、値は更新されていないことになる。上述したとおり、平均画像Aの各画素は、新たな原画像Pが入力されるたびに必ず更新されるため、原画像Pの変化に追随して変化を遂げる動的画像になる。もちろん、平均画像Aを構成する各画素の画素値は、重みwを考慮した重みつき平均として算出されるため、原画像Pの画素値が急激に変動しても、平均画像Aの画素値がこれに追随して急激に変動することはない。これは、図4に示す画素値piとaiとの変動を比べれば容易に理解できよう。しかしながら、平均画像Aは、新たな原画像Pの影響を受けて、必ず変化する画像であるので、たとえば、図1に示すような原画像P10~P40が与えられた場合、通過した車両の影響が平均画像Aにも及ぶことになる。より具体的には、赤色の車両が通過した場合、平均画像Aの一部も若干赤みを帯びるように変化する。これに対して、参照画像Uの各画素は、新たな原画像Pが入力されるたびに必ずしも更新されるものではない。すなわち、前背判定部450において、背景領域J=1と判定された部分に対しては平均画像Aを利用した更新が行われるが図10の式1B、前景領域J=0と判定された部分に対しては更新は行われない図10の式1C。このため、原画像Pの変化に追随した参照画像Uの変化は、平均画像Aに比べて緩慢になる。実際、図11および図12に示す実例を見ると、時点t3において車両Vが通過すると、平均画像A3の画素a3は、これに追随して若干変化しているが、参照画像U3の画素u3は、更新されないため、画素u2のまま変化していない。同様に、時点t10,t12を見ても、平均画像Aは原画像Pの変動に追随して変化しているのに対して、参照画像Uは変化していない。要するに、図11および図12に示す時点t1~t15の画像変遷を見ればわかるとおり、平均画像Aは、原画像Pの変動に追随して逐次変化するのに対して、参照画像Uの変化は緩慢であり、原画像Pが急激に変動しても、直ちに追随することはない。改良型の移動物体検出装置では、新たに入力された原画像Pは、上述したように追随態様が互いに異なっている平均画像Aおよび参照画像Uの双方と比較され、少なくとも一方に類似している領域は背景領域と判定されることになる。したがって、たとえば、雲の動きなどによって原画像Pの背景部分に短期的な日照変化が生じた場合でも、平均画像Aおよび参照画像Uの少なくとも一方に類似している領域は、背景領域Bと判定されることになる。その結果、日照変化などの自然現象に起因して変化した背景領域が前景領域F移動物体の領域として誤検出される可能性を低減することができ、より正確な検出を行うことができるようになる。これが、平均画像Aと参照画像Uとの双方を用いて移動物体検出を行う利点である。なお、図9に示す改良型の移動物体検出装置には、パラメータ設定部340が設けられており、ユーザは、パラメータ設定部340に対する操作入力によって、各構成要素で用いられるパラメータを任意の値に設定することができる。既に述べたとおり、この装置で設定可能なパラメータは、第1判定部420および第2判定部430の類否判定の基準となる判定パラメータたとえば、図7の円柱モデルを利用した判定を行う場合は、円柱Cの半径および長さLと、平均画像作成部130で用いる重みwと、参照画像作成部135で用いる重みvである。最適なパラメータ値は、装置の実際の使用環境に基づいて変わるので、実用上は、個々の使用環境ごとに試行錯誤によって最適なパラメータ値を設定するようにすればよい。ところで、本願発明者が行った実験によると、平均画像作成部130で用いる重みwと、参照画像作成部135で用いる重みvとについて、w>vなる関係が得られるような設定を行うと、より良好な検出結果が得られることが判明した。上述したとおり、参照画像Uは平均画像Aに比べて、原画像Pの変動に対して緩慢な追従を行うが、参照画像Uを作成する際に用いる重みvを、平均画像Aを作成する際に用いる重みwに比べて小さく設定することにより、参照画像Uの追従の緩慢度を更に上げることができる。その結果、より良好な検出結果が得られるようになるものと考えられる。なお、パラメータ設定部340を設けた場合には、w>vなる関係が得られるように、ユーザの操作入力に制限を課する機能を付加しておくのが好ましい。たとえば、w>vなる条件に反するようなパラメータ設定を行う操作入力が行われた場合には、エラー表示を行い、当該操作入力を受け付けないような機能を設けておけばよい。あるいは、たとえば、v=0.1×wというような式により、vとwとの関係を固定しておき、ユーザの操作入力によって、パラメータwのみが設定できるようにしておいてもよい。実際、本願発明者は、屋外に設置した監視カメラから毎秒30フレームの原画像を取り込んで、重みw=0.01、重みv=0.001なる設定で実験を行ったところ、日照変化による誤検出を実用的なレベルまで低減させた良好な検出結果を得ることができた。<<<§7.回転楕円体モデルを用いた類否判定>>>図9に示す改良型の移動物体検出装置における第1判定部420および第2判定部430では、2つの画素についての類否判定が行われる。このような類否判定の一般的な方法は、既に§2においていくつかの例を紹介したとおりである。ただ、本願発明者は、改良型の移動物体検出装置を設計するにあたって特に有効な類否判定方法を着想したので、以下、この方法を述べておくこの方法は、特願2009-226146号明細書に開示されている方法である。ここに示す類否判定方法は、基本的には、比較対象となる一方の画素の画素値を基準画素値、他方の画素の画素値を比較画素値として、三原色の各画素値を各座標軸にとった三次元座標系において、基準画素値に対応する座標に位置する基準点Qと比較画素値に対応する座標に位置する比較点qをとり、基準点Qを中心とする所定サイズの基準立体と比較点qとの位置関係を調べ、比較点qが基準立体の外部に位置すると判定できる場合に非類似との判定を行い、内部に位置すると判定できる場合には類似との判定を行うものである。§2で述べた図7に示す例は、基準立体として円柱Cを用いる例であるが、ここで述べる例は、基準立体として回転楕円体を用いるようにした点に特徴がある。ここでは、第1判定部420による類否判定に、この回転楕円体モデルを用いた場合を説明するが、この回転楕円体モデルは、第2判定部430による類否判定にも同様に利用することができる。図13は、図7に示す円柱Cを回転楕円体Eaに置き換えた回転楕円体モデルを示す三次元色空間図である。ここで、類似範囲を示す回転楕円体Eaは、平均画像Aの画素値を示す座標点Aを中心点とし、ベクトルVa上に定義された基準軸Zaを長軸方向とする楕円を、基準軸Zaを回転軸として回転させることにより得られる回転体である。この回転楕円体Eaのサイズは、長軸方向の長さLa1と短軸方向の長さLa2によって規定される。この回転楕円体モデルを用いた場合でも、類否判定の方法は、前述の円柱モデルを用いた場合の判定方法と同様であり、回転楕円体Eaの内部領域にプロットされる画素は、座標点Aに対応する画素値をもった画素に対して類似する画素と判定されることになる。具体的には、図示の例の場合、座標点Pに対応する画素値をもった画素は、回転楕円体Eaの外側に位置するため類似範囲外と判定され、両画素は非類似と判定される。逆に、座標点Pが回転楕円体Eaの内部に位置していれば、両画素は類似と判定される。一方、図14は、図13に示す回転楕円体モデルのバリエーションを示す三次元色空間図である。図示のとおり、このモデルにおける類似範囲を示す回転楕円体Epは、判定対象となる原画像Pの画素値を示す座標点Pを中心点とし、ベクトルVp上に定義された基準軸Zpを長軸方向とする楕円を、基準軸Zpを回転軸として回転させることにより得られる回転体である。この回転楕円体Epのサイズは、長軸方向の長さLp1と短軸方向の長さLp2によって規定される。図13に示すモデルの場合、背景を示す平均画像Aの画素値Ra,Ga,Baをプロットした座標点Aを基準として、類似範囲を示す回転楕円体Eaが定義され、比較対象となる原画像Pの画素値Rp,Gp,Bpをプロットした座標点Pが、回転楕円体Eaの内側にあるか外側にあるかを判定していたのに対して、図14に示すモデルの場合、比較対象となる原画像Pの画素値Rp,Gp,Bpをプロットした座標点Pをプロットした座標点Pを基準として、類似範囲を示す回転楕円体Epが定義され、背景を示す平均画像Aの画素値Ra,Ga,Baをプロットした座標点Aが、回転楕円体Epの内側にあるか外側にあるかを判定することになる。要するに、前者では、新たに入力された原画像Pの画素が、背景を示す過去の平均画像Aの画素の類似範囲に入っているか否かを判定する手法をとるのに対して、後者では、背景を示す過去の平均画像Aの画素が、新たに入力された原画像Pの画素の類似範囲に入っているか否かを判定する手法をとることになる。結局、両者の相違は、比較対象となる2つの画素のどちらを基準にして回転楕円体を定義するかという点だけである。そこで、以下の説明では、比較対象となる2つの画像のうち、一方の画像の画素値を基準画素値R,G,Bと表記し、他方の画像の画素値を比較画素値r,g,bと表記することにする。また、三次元色空間上で基準画素値R,G,Bをプロットした座標点QR,G,Bを基準点Qと呼び、比較画素値r,g,bをプロットした座標点qr,g,bを比較点qと呼ぶことにする。類否判定は、基準点Qを中心として定義された回転楕円体と比較点qとの位置関係に基づいてなされることになる。図15は、このような一般的な表記にしたがって、回転楕円体モデルを図示した三次元色空間図である。原点Oを起点として基準点Qへ向かうベクトルをベクトルVとすれば、基準軸Zは、このベクトルV上の軸になる。そして、回転楕円体Eは、基準点Qを中心点とし、ベクトルV上に定義された基準軸Zを長軸方向とする楕円を、基準軸Zを回転軸として回転させることにより得られる回転体である。ここで、比較点qr,g,bが回転楕円体Eの内部に位置すれば、両画素は類似と判定され、回転楕円体Eの外部に位置すれば、両画素は非類似と判定される。図13に示すモデルは、類否判定を行う際に、平均画像Ai-1上の画素ai-1の画素値を基準画素値R,G,Bとして基準点QR,G,Bをプロットし、原画像Pi上の画素piの画素値を比較画素値r,g,bとして比較点qr,g,bをプロットしたものである。これに対して、図14に示すモデルは、類否判定を行う際に、原画像Pi上の画素piの画素値を基準画素値R,G,Bとして基準点QR,G,Bをプロットし、平均画像Ai-1上の画素ai-1の画素値を比較画素値r,g,bとして比較点qr,g,bをプロットしたものである。図7に示す円柱モデルに比べて、図15に示す回転楕円体モデルは、次の2つの利点を有している。第1の利点は、時刻や天候などによる照明変動の影響を排除した、より正確な移動物体の検出が可能になる点である。本願発明者は、ビデオカメラで撮影した様々な実写画像について、図7に示す円柱モデルに基づく類否判定処理を採り入れた移動物体の検出と、図15に示す回転楕円体モデルに基づく類否判定処理を採り入れた移動物体の検出とを試みたところ、後者の検出の方がより正確である実験結果が得られた。これは、円柱モデルの場合、円柱の長手方向の端部周縁領域において、誤検出が行われる可能性が高くなるためと考えられる。図16は、円柱モデルに対する回転楕円体モデルの利点を示す断面図である。図示の例において、円柱Cの断面と回転楕円体Eの断面とを比較すると、いずれも基準点Qを中心として、基準軸Zの方向を長手方向とする図形であるが、図にハッチングを施す領域が、円柱Cでは内部領域と判定されるのに対して、回転楕円体Eでは外部領域と判定されることがわかる。したがって、比較点qが、このハッチング領域に位置する場合、一方では類似範囲内と判定され、他方では類似範囲外と判定される。もちろん、いずれの判定結果が正しいものであるかは、一概には決定できないが、少なくとも実写画像を用いた移動物体の検出処理に利用する限りにおいて、後者の判定結果が正しい結果となる可能性が高いことになる。すなわち、円柱モデルの場合、比較点qが図のハッチング領域内に位置する場合、判定対象となる画素は背景画像に類似しているため背景領域B内の画素とされることになるが、実際には、背景画像に対して非類似であり、前景領域F内の画素とすべきであった、というケースが比較的高い確率で生じてしまうものと考えられる。図16に示すとおり、断面形状で比較すると、回転楕円体Eの断面は曲面からなる楕円であるのに対して、円柱Cの断面は直線からなる矩形であり、実写画像を対象物とした類否判定の境界線としては、断面が曲面からなる楕円の方が適していることは、直感的にも理解できよう。このように、少なくとも、実写画像に基づく移動物体の検出処理に利用する場合、円柱モデルよりも回転楕円体モデルの方が検出精度が向上するという第1の利点が得られる。一方、回転楕円体モデルを採用することにより得られる第2の利点は、演算負担の軽減を図ることができる点である。この第2の利点の詳細については、§8で述べることにするが、ここでは、演算負担の軽減に寄与する工夫のひとつを簡単に述べておく。この工夫は、回転楕円体Eのサイズを、三次元座標系の原点Oと基準点Qとの距離Dに比例した値に設定することである。具体的には、図15に示すとおり、距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、基準軸Zを中心に回転させることにより回転楕円体Eを定義すればよい。このような条件で回転楕円体Eのサイズを決定すると、基準点Qが原点Oから離れれば離れるほど、回転楕円体Eのサイズは大きくなり、それだけ類似範囲が広くなる。これは、より明度の高い部分、すなわち、画像のハイライト部ほど、類似範囲が広くなることを意味するが、実用上、そのような類似範囲設定がなされても支障は生じない。一方、類似範囲内か否か、すなわち、回転楕円体Eの内部か外部かを判定するための演算負担は、長軸半径をα=hD、短軸半径をβ=kDと設定することによりすなわち、楕円のサイズを、距離Dに比例した値に設定することにより、大幅に軽減することができる。その詳細は§8で述べることにする。ところで、本願発明者が行った実験によると、回転楕円体Eのサイズを、距離Dに比例した値に設定する場合は、図13のモデルよりも、図14のモデルを用いた方が、検出精度が向上することが判明した。図14のモデルを採用した場合、図15において、基準点QR,G,Bは、判定対象となる原画像Pの画素の画素値を基準画素値としてプロットした点になり、比較点qr,g,bは、背景となる平均画像A第2判定部430に適用する場合は、参照画像Uの画素の画素値を比較画素値としてプロットした点になる。上述したとおり、図13のモデルと図14のモデルは、原画像Pの画素と平均画像A参照画像Uの画素とのいずれを基準にして回転楕円体を定義するか、という点において異なる。しかも、回転楕円体のサイズを、距離Dに比例した値に設定すれば、定義される回転楕円体のサイズが、両者では異なることになる。すなわち、図13に示すモデルにおける回転楕円体Eaの長軸方向の長さLa1および短軸方向の長さLa2は、ベクトルVaの長さに、それぞれパラメータ2hおよび2kを乗じることにより得られるのに対して、図14に示すモデルにおける回転楕円体Epの長軸方向の長さLp1および短軸方向の長さLp2は、ベクトルVpの長さに、それぞれパラメータ2hおよび2kを乗じることにより得られる。したがって、図示の例の場合、回転楕円体Eaのサイズよりも、回転楕円体Epのサイズの方が大きくなっており、より広い類似範囲が設定されている。もちろん、常識的には、図14に示すモデルよりも図13に示すモデルを採用するのが一般的であろう。すなわち、予め背景とみなすことができる類似範囲を回転楕円体Eaとして定義しておき、新たに入力された画像の画素を示す点Pが、この類似範囲内に入るか否かを判定することにより、当該画素が移動物体上の点であるか否かを認識する、という方針を採るのが、一般的な考え方である。これに対して、図14に示すモデルは、新たに入力された画像の画素の類似範囲を回転楕円体Ebとして定め、その中に、背景画素を示す点Aが入るか否かを判定する手法を採るものであり、一般常識に基づく手法とは逆の手法と言うことができる。しかしながら、本願発明者が行った実験によると、少なくとも実写画像を用いた移動物体の検出処理に利用する限りにおいて、図13に示すモデルを採用するよりも、図14に示すモデルを採用した方が、より正しい判定結果が得られた。これは、これまで述べた方法で、平均画像Aの画素値ai第2判定部430に適用する場合は、参照画像Uの画素値uiを決定すると、様々な要因により、新たに入力されてくる原画像Pの画素値piが大きく変動しても、画素値ai,uiの値が当該変動に追従するまでにはある程度の時間がかかるため、類似範囲の広狭を、原画像Pの画素値piに応じて定めた方が、正しい判定が可能になるためと考えられる。結局、改良型の移動物体検出装置における画素の類否判定を、基準立体モデル円柱モデルや回転楕円体モデルを利用して行う場合、第1判定部420は、原画像Pの画素の画素値を基準画素値、平均画像Aの画素の画素値を比較画素値として判定を行い、第2判定部430は、原画像Pの画素の画素値を基準画素値、参照画像Uの画素の画素値を比較画素値として判定を行うようにするのが好ましい。以上、移動物体の検出に利用する画素の類否判定処理方法において、円柱モデルを採用する方法に比べて、回転楕円体モデルを採用する利点として、より正確な移動物体の検出が可能になる、という第1の利点を述べるとともに、演算負担を軽減することができるという第2の利点についても簡単に説明した。続く§8では、この第2の利点についての詳細を述べることにする。<<<§8.回転楕円体モデルを利用した実用的な類否判定手順>>>一般に、三次元空間上において、ある1点が、円柱の内側にあるか外側にあるかを判定する場合も、回転楕円体の内側にあるか外側にあるかを判定する場合も、座標値を用いた幾何学演算が必要になる。しかしながら、回転楕円体の内側にあるか外側にあるかを判定する場合は、演算式を工夫し、また、一部に近似演算を取り込むことにより、演算負担を大幅に軽減することが可能である。ここでは、図15に示す回転楕円体モデルにおいて、比較点qr,g,bが、基準点QR,G,Bを中心として定義された回転楕円体Eの内側にあるか外側にあるかを判定する実用的な手順を述べる。ここで述べる手順は、3つの条件判定段階によって構成される。図17は、そのうちの第1の条件判定段階の原理を示す断面図である。図において、下方に示す原点Oは、図15の三次元色空間図に示されているRGB三次元直交座標系の原点Oであり、当該座標系上に、基準点QR,G,Bと比較点qr,g,bとがプロットされている座標軸は図示省略。基準軸Zは、原点Oと基準点Qとを結ぶ線であり、回転楕円体Eは、この基準軸Zを長軸方向とし、基準点Qを中心とする楕円を、基準軸Zを回転軸として回転させることによって得られる回転体である。図に断面を示す楕円は、長軸半径をα、短軸半径をβとする楕円であるが、§7で述べたとおり、原点Oと基準点Qとの距離をDとしたときに、所定のパラメータh,kについて、α=hD但し、h<1、β=kD但し、k<hとなるように設定されている。第1の条件判定段階では、比較点qが、図のハッチング領域内にあるかどうかの判定が行われる。すなわち、図の下方のハッチング領域は、原点Oを中心とし半径D-αの球C1の内側領域であるから、原点Oと比較点qとの距離をdとすれば、D-α>dなる条件が満足されれば、比較点qは、球C1の内側に位置することになる。球C1は、回転楕円体Eの原点Oに対する最近接点E1において、回転楕円体Eに接する球であり、球C1の内側領域は、当然、回転楕円体Eの外側領域になる。一方、図の上方のハッチング領域は、原点Oを中心とし半径Dαの球C2の外側領域であるから、原点Oと比較点qとの距離をdとすれば、Dα<dなる条件が満足されれば、比較点qは、球C2の外側に位置することになる。球C2は、回転楕円体Eの原点Oに対する最遠接点E2において、回転楕円体Eに接する球であり、球C2の外側領域は、当然、回転楕円体Eの外側領域になる。図18は、図17に示す第1の条件判定原理に基づく判定式を示す図である。図18は、比較点qが球C1の内側領域にあるための条件式D-α>dの変形プロセスを示す。すなわち、α=hDを代入すると、条件式は1-hD>dとなり、両辺を2乗すると、1-h2D2>d2式2が得られる。一方、図18は、比較点qが球C2の外側領域にあるための条件式Dα<dの変形プロセスを示す。すなわち、α=hDを代入すると、条件式は1hD<dとなり、両辺を2乗すると、1h2D2<d2式3が得られる。ここで、基準点Qの座標値がR,G,B、比較点qの座標値がr,g,bであるから、図18に示すように、D=√R2G2B2d=√r2g2b2であり、両辺を2乗すると、D2=R2G2B2d2=r2g2b2となる。したがって、式2,式3を書き直せば、図18に示すように、1-h2R2G2B2>r2g2b2式2′1h2R2G2B2<r2g2b2式3′なる条件式が得られる。ここで、hは予め設定された所定のパラメータであり、R,G,Bおよびr,g,bは、基準点Qおよび比較点qの座標値すなわち、比較対象となる2つの画素の三原色の画素値である。式2′および式3′は、2乗演算と和差演算のみで構成される演算式であるから、演算負担は非常に軽いものになる。こうして、式2′もしくは式3′のいずれかの条件を満足する結果が得られた場合、比較点qは、回転楕円体Eの外側に位置すると判断できるので、基準点Qの類似範囲外との判定を行うことができる。一方、式2′および式3′のいずれの条件も満足しない場合は、図17において、比較点qは、球C1の外側かつ球C2の内側の領域に位置することになる。したがって、その場合には、以下に述べる第2の条件判定段階が引き続き行われる。図19は、この第2の条件判定段階の原理を示す断面図である。第2の条件判定段階では、比較点qが、図19のハッチング領域内にあるかどうかの判定が行われる。このハッチング領域は、原点Oを頂点、基準軸Zを中心軸とし、回転楕円体Eと円で接する円錐ξの外側領域として定義される領域である。円錐ξの半頂角は、図示のとおり角φになる。図示の断面図では、円錐ξが接点Tにおいて楕円Eに接している状態が示されているが、実際には、円錐ξと回転楕円体Eとは、接線円Tで線接触することになる。三次元空間上での円錐ξと回転楕円体Eとの関係は、コーンに卵を挿入した状態を考えれば容易に理解できよう。結局、この第2の条件判定段階では、比較点qが、円錐ξの内側にあるか、外側にあるかを判定することになる。円錐ξと回転楕円体Eとが接する接線円Tの位置は、幾何学的な演算によって求めることができ、そのような接線円Tを通る円錐ξを示す幾何学的な方程式を求めることも可能である。しかしながら、そのような方程式は比較的複雑な式になる。そこで、ここでは、円錐ξの代わりに、図19に破線で示す円錐ξ′を利用した近似的な判定を行うことにより、演算負担の軽減を図ることにする。図19に実線で描いた円錐ξは、接線円Tで回転楕円体Eに接する図形であるが、破線で描いた円錐ξ′は、交線円Sで回転楕円体Eと交差する図形である。ここで、交線円Sは、原点Oを中心とする半径Dの球面C3と回転楕円体Eとの交線を構成する円であり、基準点Qと交線円S上の任意の点との距離は図示のとおりfである。結局、この近似的な判定では、原点Oを頂点、基準軸Zを中心軸とし、交線円Sを通る円錐ξ′の外側に比較点qが位置するか否かという判定が行われることになる。図示のとおり、円錐ξ′の半頂角φ′は、円錐ξの半頂角φよりも若干小さくなり、円錐ξ′は、円錐ξよりも若干小さな円錐になる。このため、円錐ξ′の外側かつ円錐ξの内側の点の場合、実際は回転楕円体Eの内部の点であるにもかかわらず、円錐ξ′を代用する近似的な判定方法では、回転楕円体Eの外部の点と判定されてしまうことになる。したがって、回転楕円体Eの内外を厳密に判定するという観点では、円錐ξ′を代用する判定方法は不正確である。ただ、そもそも回転楕円体E自体が、移動物体を厳密に検出する機能をもつ立体であるわけではなく、あくまでも1つの指標を提供するものであるから、円錐ξ′を代用する近似的な判定方法を用いたとしても、大きな支障は生じない。そこで、比較点qが、この円錐ξ′の内側にあるのか外側にあるのかを幾何学的に判定するための具体的な方法を考えてみる。いま、図20に示すように、三辺の長さがそれぞれa,b,cである一般的な三角形を考える。ここで、2辺b,cの挟角をφとすると、余弦定理により、cosφ=b2c2-a2/2bc式4が成り立つ。続いて、図20に示す三角形OQSを考える。この図20は、図19に示す三角形OQSをそのまま抜き出して示したものである。交線円Sは、前述したとおり、原点Oを中心とする半径Dの球面C3と回転楕円体Eとの交線であるから、2点OSの距離は、2点OQの距離に等しくDになる。したがって、三角形OQSは二等辺三角形になり、2つの等辺の挟角はφ′、第3の辺の長さはfである。この三角形OQSに、式4に示す余弦定理を適用すると、cosφ′=D2D2-f2/2D2式5が成り立つ。ここで、図19を見ると、長さfは、回転楕円体Eの短軸半径βに近い値になることがわかる。そこで、近似的に、f≒βとすると、β=kDなので、f≒kDになり、図20の下段右に示すとおり、cosφ′≒1-k2/2式6なる式が得られる。次に、図21に示す三角形OQqを考える。この図21は、図19に示す4点O,Q,S,qをそのまま抜き出してプロットしたものである。2点OQ間の距離はDであり、2点Oq間の距離はdである。ここで、この長さD,dをもった2辺の挟角をθとし、第3の辺の長さをtとする。また、線分OQと線分OSとのなす角φ′は、前述した円錐ξ′の半頂角である。さて、図19において、円錐ξの代わりに円錐ξ′を用いて、比較点qがハッチング領域にあるための条件を求めると、θ>φ′であればよいことがわかる。ここで、座標値画素値R,G,B,r,g,bが正の値をとるものとすれば、角度θおよびφ′はいずれも鋭角になる。図21の右側にも記載したとおり、一般に、θ,φ′を正の鋭角とすれば、θ>φ′の場合、cosθ<cosφ′であり、cos2θ<cos2φ′である。式6より、cos2φ′≒1-k2/22であるから、結局、cos2θ<1-k2/22式7が成り立てば、比較点qは、円錐ξ′の外側領域にあると近似的に判断できる。続いて、図21に示す三角形OQqについて、式4の余弦定理を適用すれば、図21の下段に示すように、cosθ=D2d2-t2/2Dd式8が成り立つ。このとき、図22にも記載したとおり、D2=R2G2B2d2=r2g2b2t2=R-r2G-g2B-b2だから、これを式8に代入して整理すると、cosθ=γ/Dd式9但し、γ=RrGgBbが得られる。式9の両辺を2乗すれば、cos2θ=γ2/D2d2式10である。よって、式7より、γ2/D2d2<1-k2/22式11が成り立てば、比較点qは、円錐ξ′の外側領域にあると近似的に判断できる。ここで、γ,D,dを画素値を用いて書き直せば、図22の下段に示すように、RrGgBb2/R2G2B2r2g2b2<1-k2/22式11′が成り立てば、比較点qは、円錐ξ′の外側領域にあると近似的に判断できる。既に述べたとおり、図19において、円錐ξ′は、円錐ξに近似的に代用することができるので、上記式11′が成り立てば、比較点qは、回転楕円体Eの外側領域にある、との近似的判定を行うことができる。もちろん、この近似的な判定結果は、円錐ξ′を円錐ξに代用し、かつ、f≒βとの近似を行うことによって得られる判定結果であるので、厳密な判定結果ではない。しかしながら、そもそも回転楕円体E自体が、移動物体を厳密に検出する機能をもつ立体であるわけではないので、このような近似的な判定を行っても、大きな支障は生じない。別言すれば、上述した近似的な判定手法は、比較点qが回転楕円体Eの外部にあるかどうかを判定する代わりに、その近似体の外部に位置するかどうかを判定する手法ということになる。このようにして、式11′を満足する結果が得られた場合、比較点qは、回転楕円体Eの近似体の外側に位置すると判断できるので、基準点Qの類似範囲外との判定を行うことができる。一方、式11′を満足しない場合は、比較点qは、図19においてハッチングが施されていない円錐ξの内側に位置すると近似的に判断されたことになるので、最終的に、第3の条件判定段階が行われる。図23は、この第3の条件判定段階の原理を示す断面図である。第3の条件判定段階では、比較点qが、図23のハッチング領域内にあるかどうかの判定が行われる。このハッチング領域は、回転楕円体Eの外側領域であり、いわば第3の条件判定段階は、§7で述べた回転楕円体モデルを用いた類否判定として行うべき本来の判定段階ということになる。このため、条件判定に用いられる式は、楕円の方程式を用いた式になり、これまでの条件判定式に比べると、若干複雑にならざるを得ない。いま、図23に示すように、基準点Qの位置に原点をもち、紙面上に定義されたローカルなXY二次元直交座標系をとり、この二次元座標系上において、楕円Eと比較点qとの位置関係を考えることにしよう。ここで、X軸は回転楕円体Eの長軸方向の座標軸であり、Y軸は回転楕円体Eの短軸方向の座標軸である。もちろん、基準点QのXY座標系での座標は0,0になる。したがって、比較点qは、RGB三次元直交座標系ではqr,g,bなる三次元の座標値をもち、XYローカル座標系ではqx,yなる二次元の座標値をもつ。ここで、比較点qの座標値x,yを、2点OQ間の距離D、2点Oq間の距離d、そして挟角θを用いて表すと、図23の右下にも記載したとおり、x=D-dcosθ、y=dsinθ式12が成り立つ。一方、図24に示すとおり、楕円の公式からは、x2/α2y2/β2=1式13が得られるので、fx,y=β2x2α2y2-α2β2式14なる判定式を定義すれば、比較点qは、fx,y<0なら楕円内部の点fx,y=0なら楕円上の点fx,y>0なら楕円外部の点ということになる。したがって、β2x2α2y2-α2β2>0式15であれば、比較点qは、楕円外部の点ということになる。ここで、XYローカル座標系を、基準軸Zを中心軸として回転しても、上式が成り立つことに変わりはないので、式15は、三次元立体にも拡張することができる。すなわち、式15の条件が満足されていれば、比較点qは、回転楕円体Eの外側領域の点ということになる。なお、式15において、各変数の値は、図24の下半分に示すように、次のようにして求めることができる。まず、α=hD、β=kDであるから、α2=h2D2、β2=k2D2式16である。一方、式12より、x=D-dcosθだから、両辺を2乗すると、x2=D2d2cos2θ-2Ddcosθとなり、これに式9および式10を代入すれば、x2=D2γ2/D2-2γ式17が得られる但し、γ=RrGgBb。また、式12より、y=dsinθだから、両辺を2乗すると、y2=d2sin2θ=d21-cos2θとなり、これに式10を代入すれば、y2=d2-γ2/D2式18が得られる但し、γ=RrGgBb。かくして、式16~式18によって、α2,β2,x2,y2の各値を求めた後に、式15の条件判定式を演算すれば、比較点qが、回転楕円体Eの内側領域の点か、外側領域の点かを判定することができる。図25は、これまで述べてきた類否判定方法の手順を示す流れ図である。この類否判定方法によれば、三原色の各画素値R,G,Bを有する第1のカラー画素と、三原色の各画素値r,g,bを有する第2のカラー画素と、の類否判定を行うことができる。なお、実際には、この図25に示す各手順は、コンピュータによって実行されることになる。まず、ステップS11では、回転楕円体Eのサイズを決めるパラメータh,kを設定する。前述したとおり、パラメータhは、回転楕円体Eの長軸半径αを決定する値h<1であり、長軸半径はα=hDなる式で定められる。また、パラメータkは、回転楕円体Eの短軸半径βを決定する値k<hであり、短軸半径はβ=kDなる式で定められる。これらパラメータh,kの値は、背景画像の特徴たとえば、屋外か屋内か、街中か郊外か、といった種別や、移動物体として侵入してくる前景画像の特徴たとえば、車両か人間か、といった種別を考慮して、できるだけ高い検出精度が得られるような値に適宜設定することになる。本願発明者が行った実験によれば、一般的な用途の場合、h=0.45~0.20、k=0.15~0.05程度の値に設定すれば、比較的高い検出精度が得られる。図9に示す移動物体検出装置の場合、パラメータ設定部340に対する操作入力により、パラメータh,kの値を任意に設定することができる。続く、ステップS12では、基準点QR,G,Bと比較点qr,g,bの座標値R,G,Bおよびr,g,bが取り込まれる。これらの座標値は、それぞれ2つのカラー画像の三原色の画素値として与えられた値である。これら座標値が取り込まれたら、準備段階は完了である。すなわち、コンピュータには、パラメータhh<1およびkk<h、画素値R,G,Bおよびr,g,bの値が取り込まれた状態となり、これらの値を用いて、以下の演算処理を実行する準備が整ったことになる。以下、ステップS13~S18に示す手順は、図15に示すように、三原色の各画素値を各座標軸にとった三次元座標系において、画素値R,G,Bに対応する座標をもつ基準点Qを中心とする所定サイズの回転楕円体Eと、画素値r,g,bに対応する座標をもつ比較点qとについて、後者が前者の内部にあるか外部にあるかを調べる位置関係判定段階の手順である。具体的には、既に述べたとおり、三次元座標系の原点Oと基準点Qとの距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この基準軸Zを中心に回転させることにより得られる回転楕円体Eを定義する。そして、比較点qがこの回転楕円体Eの内部にあるか外部にあるかを判定することになる。なお、実際には、前述したような近似を行うため、論理的には、回転楕円体Eの代わりに、その近似体について、内部にあるか外部にあるかの判定が行われることになる。この位置関係判定段階は、3つの条件判定段階によって構成される。まず、ステップS13に示す第1の条件判定段階では、図17に示すように、原点Oを中心とし半径D-αの球C1の内側に比較点qが位置するもしくは原点Oを中心とし半径Dαの球C2の外側に比較点qが位置するという条件を満足するか否かの判定が行われる。そのためには、図18で説明したとおり、1-h2D2>d2式2もしくは1h2D2<d2式3という条件を満足するか否かを判定するための演算を行えばよい。ここで、D2=R2G2B2d2=r2g2b2であるから、結局、ステップS13のブロック内に記載したように、1-h2R2G2B2>r2g2b2式2′もしくは1h2R2G2B2<r2g2b2式3′という条件を満足するか否かを判定する演算を行えばよい。このステップS13に示す第1の条件判定段階において、肯定的な判定がなされた場合式2′もしくは式3′の条件を満足する場合は、ステップS14からステップS19へと進み、位置関係判定段階は終結する。一方、否定的判定がなされた場合には、ステップS14からステップS15へと進む。ステップS15に示す第2の条件判定段階では、図19に示すように、原点Oを頂点、基準軸Zを中心軸とし、回転楕円体Eと円で接する円錐ξの外側に比較点qが位置するという第2の条件を満足するか否かの判定が行われる。ただ、前述したとおり、実用上は、円錐ξの代わりに円錐ξ′を用い、しかもf≒βと近似する式が用いられる。このような近似的な判定手法を採る場合は、第2の条件は回転楕円体Eの近似体の外側に比較点qが位置するという条件になる。この近似的な判定手法を採るのであれば、図22に示すとおり、γ2/D2d2<1-k2/22式11という条件を満足するか否かを判定するための演算を行えばよい。ここで、D2=R2G2B2d2=r2g2b2、γ=RrGgBbであるから、結局、ステップS15のブロック内に記載したように、RrGgBb2/R2G2B2r2g2b2<1-k2/22式11′という条件を満足するか否かを判定する演算を行えばよい。このステップS15に示す第2の条件判定段階において、肯定的な判定がなされた場合式11′の条件を満足する場合は、ステップS16からステップS19へと進み、位置関係判定段階は終結する。一方、否定的判定がなされた場合には、ステップS16からステップS17へと進む。ステップS17に示す第3の条件判定段階では、図23に示すように、回転楕円体Eの外側に比較点qが位置するという第3の条件を満足するか否かの判定が行われる。そのためには、ステップS17のブロック内に記載したように、β2x2α2y2-α2β2>0式15という条件を満足するか否かを判定するための演算を行えばよい。ここで、各変数は、図24で説明したとおり、α2=h2D2、β2=k2D2式16x2=D2γ2/D2-2γ式17y2=d2-γ2/D2式18であり、また、既に述べたとおり、D2=R2G2B2d2=r2g2b2γ=RrGgBbである。このステップS17に示す第3の条件判定段階において、肯定的な判定がなされた場合式15の条件を満足する場合は、ステップS18からステップS19へと進み、位置関係判定段階は終結する。一方、否定的判定がなされた場合には、ステップS18からステップS20へと進み、やはり位置関係判定段階は終結する。ここで注目すべき点は、以上の第1~第3の条件判定段階で行われる条件判定式には、加減乗除演算のみしか含まれていない点である。すなわち、一般的な幾何学演算に含まれている開平演算や三角関数演算は一切含まれていない。したがって、演算装置の負担は比較的軽いものになり、安価なプロセッサを用いた場合でも、リアルタイム処理が可能になる。最後に、ステップS19もしくはステップS20に示す判定結果出力段階が行われる。すなわち、上述したとおり、第1の条件判定段階ステップS13、第2の条件判定段階ステップS15、第3の条件判定段階ステップS17のいずれかにおいて肯定的判定がなされた場合には、ステップS19へと進み、第1のカラー画素と第2のカラー画素とが非類似との判定結果が出力される。一方、第1の条件判定段階ステップS13、前記第2の条件判定段階ステップS15、前記第3の条件判定段階ステップS17のすべてにおいて否定的判定がなされた場合には、ステップS20へ進み、第1のカラー画素と第2のカラー画素とが類似するとの判定結果が出力される。要するに、この図25の流れ図に示す手順によれば、回転楕円体Eもしくはその近似体の外部に比較点qがあると判断された場合には、第1のカラー画素と第2のカラー画素とが非類似との判定結果が出力され、回転楕円体Eもしくはその近似体の内部に比較点qがあると判断された場合には、第1のカラー画素と第2のカラー画素とが類似するとの判定結果が出力されることになる。ステップS15の第2の条件判定段階は、ステップS13の第1の条件判定段階において否定的判定がなされた場合にのみ実行され、同様に、ステップS17の第3の条件判定段階は、ステップS15の第2の条件判定段階において否定的判定がなされた場合にのみ実行されることになるので、無駄な演算を省くことが可能になり、演算負担の軽減を図ることができる。特に、パラメータh,kの値を固定値にした運用を行う場合、ステップS13の演算式2′,3′に用いる1-h2および1h2の値、ステップS15の演算式11′に用いる1-k2/22の値は、予め決められた定数として用意しておくことができるので、演算負担を更に軽減することができる。実際には、まず、R2G2B2の値およびr2g2b2の値を計算してしまえば、ステップS13の演算式2′,3′に係る演算は、上記定数を利用した極めて単純な演算になる。同様に、ステップS15の演算式11′に係る演算も、RrGgBbの値を計算してしまえば、既に計算済みのR2G2B2の値およびr2g2b2の値と、上記定数を利用することにより、極めて単純な演算になる。<<<§9.前背判定部の実用的な構成>>>ここでは、図9に示す前背判定部450において、画素の類否判定に、§8で述べた実用的な類否判定手順を採用する場合の構成を説明する。図26は、この前背判定部450の構成を示すブロック図である。既に述べたとおり、前背判定部450は、第1判定部420、第2判定部430、総合判定部440によって構成されている。ここで、第1判定部420は、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の平均画像Ai-1を構成する画素の画素値ai-1とを比較して類似するか否かを判定する処理を行い、第2判定部430は、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の参照画像Ui-1を構成する画素の画素値ui-1とを比較して類似するか否かを判定する処理を行う。一方、総合判定部440は、図示の表に示すとおり、第1判定部420による判定結果および第2判定部430による判定結果の少なくとも一方が類似を示す場合には背景画素との判定を行い判定結果J=1を出力し、双方とも非類似を示す場合には前景画素との判定を行い判定結果J=0を出力する。第1判定部420および第2判定部430が、類似との判定結果を論理1、非類似との判定結果を論理0として出力するようにしておけば、総合判定部440は、これらの論理値の論理和を判定結果Jとして出力することができるので、総合判定部440を論理和回路によって構成することができる。図27は、図26に示す第1判定部420の詳細構成を示すブロック図である。この第1判定部420は、§8で述べた回転楕円体モデルを利用した実用的な類否判定手順を採用し、図示のとおり、基本演算部421、第1の条件判定部422、第2の条件判定部423、第3の条件判定部424、判定値出力部425、パラメータ保持部426によって構成されている。パラメータ保持部426は、パラメータhおよびkの値を保持する構成要素であり、図9に示すパラメータ設定部340を設けた場合には、このパラメータ設定部340によって設定されたパラメータhおよびkの値がパラメータ保持部426に書き込まれることになる。パラメータhおよびkの値を固定にする場合には、当該固定値がパラメータ保持部426に保持される。基本演算部421は、図9に示す画素値読出部410が読み出した基準画素値R,G,B画素piの画素値と比較画素値r,g,b画素ai-1の画素値とに基づいて、D2=R2G2B2およびd2=r2g2b2なる演算を行い、値D2および値d2を算出する構成要素である。ここで算出された値D2および値d2は、第1の条件判定部422、第2の条件判定部423、第3の条件判定部424における演算で利用される。第1の条件判定部422は、図17に示すように、原点Oを中心とし半径D-αの球C1の内側に比較点qが位置するもしくは原点Oを中心とし半径Dαの球C2の外側に比較点qが位置するという第1の条件を満足するか否かを判定する構成要素である。そのために、基本演算部421が算出した値D2および値d2とパラメータ保持部426に保持されているパラメータhの値とを用いて、1-h2D2>d2式2もしくは1h2D2<d2式3という条件を満足するか否かを判定するための演算を行う。そして、条件満足との肯定的判定が行われた場合、当該判定結果は、判定値出力部425に報知される。一方、条件満足せずとの否定的判定が行われた場合、当該判定結果は、第2の条件判定部423に報知される。第2の条件判定部423は、第1の条件判定部422から、条件満足せずとの否定的判定結果が報知された場合に限り、条件判定処理を実行する構成要素である。ここで行われる処理の目的は、図19に示すように、原点Oを頂点、基準軸Zを中心軸とし、回転楕円体Eと円で接する円錐ξの外側に比較点qが位置するという第2の条件を満足するか否かを判定することにある。ただ、実用上は、演算負担を軽減するために、円錐ξの代わりに円錐ξ′を用いた判定を行うようにするのが好ましい。すなわち、原点Oを中心とする半径Dの球面C3と回転楕円体Eとの交線を構成する円を交線円Sとしたときに、原点Oを頂点、基準軸Zを中心軸とし、交線円Sを通る円錐ξ′の外側に比較点qが位置するという第2の条件を満足するか否かを判定すればよい。このとき、演算負担を更に軽減するために、交線円S上の任意の点と基準点Qとの距離fが短軸半径βに等しいとおいた近似式を用いて近似的な判定を行うようにするのが好ましい。このような近似的な判定を行う場合、第2の条件判定部423は、画素値読出部410が読み出した基準画素値R,G,Bおよび比較画素値r,g,bに基づいて、γ=RrGgBbなる演算を行い、値γを算出し、更に、基本演算部421が算出した値D2および値d2とパラメータ保持部426に保持されているパラメータkの値と、算出した値γとを用いて、γ2/D2d2<1-k2/22式11という条件を満足するか否かを判定するための演算を行えばよい。結局、この式11の演算式で示される条件は、原点Oを頂点、基準軸Zを中心軸とし、回転楕円体Eと円で接する円錐ξの近似体の外側に比較点qが位置するという条件になる。そして、この第2の条件判定部423において、条件満足との肯定的判定が行われた場合、当該判定結果は、判定値出力部425に報知される。一方、条件満足せずとの否定的判定が行われた場合、当該判定結果は、γの値とともに、第3の条件判定部424に報知される。第3の条件判定部424は、第2の条件判定部423から、条件満足せずとの否定的判定結果が報知された場合に限り、条件判定処理を実行する構成要素である。ここで行われる処理の目的は、図23に示すように、回転楕円体Eの外側に比較点qが位置するという第3の条件を満足するか否かを判定することにある。そのため、第3の条件判定部424は、基本演算部421が算出した値D2および値d2と、パラメータ保持部426に保持されているパラメータhおよびkの値と、第2の条件判定部423が算出した値γとを用いて、α2=h2D2、β2=k2D2式16x2=D2γ2/D2-2γ式17y2=d2-γ2/D2式18なる演算を行い、当該演算結果を用いて、β2x2α2y2-α2β2>0式15という条件を満足するか否かを判定するための演算を行う。そして、この第3の条件判定部424における条件満足との肯定的判定結果もしくは条件満足せずとの否定的判定結果は、判定値出力部425に報知される。判定値出力部425は、第1の条件判定部422、第2の条件判定部423、第3の条件判定部424のいずれかが条件満足との肯定的判定を行った場合に、非類似との判定結果論理値0を出力し、第1の条件判定部422、第2の条件判定部423、第3の条件判定部424のすべてが条件満足せずとの否定的判定を行った場合に、類似との判定結果論理値1を出力する。以上、図27を参照して、第1判定部420の詳細構成を述べたが、第2判定部430の構成もこれと全く同様である。なお、第1判定部420による類否判定に用いられる回転楕円体と第2判定部430による類否判定に用いられる回転楕円体とは、両判定部における判定基準を統一するために、同一サイズの回転楕円体とし、三次元座標系上の同一位置に配置するのが好ましい。既に述べたとおり、上記条件判定式には、加減乗除演算のみしか含まれていないため、一般的な幾何学演算に必要とされる開平演算や三角関数演算を行う必要はない。したがって、演算装置の負担は比較的軽いものになり、前背判定部450を安価なプロセッサを用いて構成した場合でも、リアルタイム処理が可能になる。なお、これまで述べた条件判定式には、不等号>や不等号<を含む演算式が含まれているが、これらの演算式の一部もしくは全部について、不等号>の代わりに不等号≧を用いた演算式、もしくは不等号<の代わりに不等号≦を用いた演算式を用いてもかまわない。すなわち、各条件判定式は、比較点qが何らかの図形の内側領域にあるか外側領域にあるかを示すものであるが、当該図形の境界面上の点は、当該図形の内側領域に含ませても、外側領域に含ませてもかまわない。また、これまで述べた実施形態では、三原色を示す符号として、便宜上、R,G,Bもしくはr,g,bという符号を用いているが、これらの符号は、必ずしも赤,緑,青なる三原色を意味するものではない。一般的には、ビデオカメラで撮影した画像は、赤,緑,青なる三原色の各画素値をもった画素の集合からなるデータとして与えられるが、本発明は、赤,緑,青なる三原色を用いた表色系に限定されるものではなく、他の色からなる三原色を用いた別な表色系にも適用可能である。その場合、各演算式に用いられているR,G,Bもしくはr,g,bという符号は、当該別な表色系の三原色の画素値を示すことになる。<<<§10.本発明に係る日照変化領域検出の基本原理>>>さて、これまで§1~§3において、従来の一般的な移動物体検出装置について述べ、§4~§9において、その改良型の移動物体検出装置について述べた。本発明に係る日照変化領域検出装置は、この改良型の移動物体検出装置を開発するプロセスにおける副産物として案出されたものであり、実際、その構成要素の大部分は、改良型の移動物体検出装置の構成要素と共通する。そこで、ここでは、本発明に係る日照変化領域検出の基本原理を説明する。まず、図26に示す前背判定部450において実行される前景画素/背景画素の判定処理と、本発明に係る日照変化領域検出装置において実行される日照変化画素/非日照変化画素の判定処理との関係を説明しよう。既に述べたとおり、図26において、第1判定部420は、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の平均画像Ai-1を構成する画素の画素値ai-1とを比較して類似するか否かを判定する処理を行い、第2判定部430は、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の参照画像Ui-1を構成する画素の画素値ui-1とを比較して類似するか否かを判定する処理を行う。そして、総合判定部440は、図示の表に示すとおり、第1判定部420による判定結果および第2判定部430による判定結果の少なくとも一方が類似を示す場合には背景画素との判定を行い判定結果J=1を出力し、双方とも非類似を示す場合には前景画素との判定を行い判定結果J=0を出力する。図28は、この総合判定部440による判定結果のバリエーションを回転楕円体モデルで示す図であり、CASE1~4までの4通りの判定結果を示すものである。ここで、点piは第i番目の原画像の画素値、点ai-1は第i-1番目の平均画像の画素値、点ui-1は第i-1番目の参照画像の画素値を、それぞれ色空間三原色を各座標軸にとった三次元座標系上にプロットした点である。また、符号Eで示されている楕円は、基準立体として用いられる回転楕円体Eを示しており、この回転楕円体Eは点piが中心点となるように配置されている。第1判定部420は、点ai-1が点piの類似範囲にあるか否か、すなわち、点ai-1が回転楕円体Eの内部にあるか否かを判定する処理を行い、第2判定部430は、点ui-1が点piの類似範囲にあるか否か、すなわち、点ui-1が回転楕円体Eの内部にあるか否かを判定する処理を行う。具体的には、CASE1のように、点ai-1,点ui-1がともに回転楕円体Eの内部双方とも類似○にあれば、総合判定部440はJ=1背景画素との判定を行う。また、CASE2もしくは3のように、一方が内部類似○、他方が外部非類似×にある場合も、総合判定部440はJ=1背景画素との判定を行う。これに対して、CASE4のように、点ai-1,点ui-1がともに回転楕円体Eの外部双方とも非類似×にあれば、総合判定部440はJ=0前景画素との判定を行う。結局、この図28に示されている4通りの判定結果CASE1~4は、図26の総合判定部440の表に示されている4通りの判定結果に対応する。このように、平均画像Aとの比較を行う第1判定部420および参照画像Uとの比較を行う第2判定部430のうち、少なくとも一方で類似○との判定がなされた場合に背景画素と判定する利点は、既に§6で説明したとおり、原画像Pの背景部分に短期的な日照変化等が生じた場合でも、当該変化領域を構成する画素を前景画素として誤検出する可能性を低減することができ、より正確な移動物体の検出を行うことができるようになる点である。これは、裏を返せば、一方では類似○との判定がなされたのに、他方では非類似×との判定がなされたケースには、日照変化の影響を受けたケースが含まれていると言える。本願発明者は、このような点に着眼して、屋外の様々な環境で撮影された多数の動画画像について解析を行ったところ、新たに入力された原画像のうち、平均画像には類似、参照画像には非類似であるが類似境界の所定近傍範囲内にあるという条件を満たす領域は、日照変化に起因した変動が生じている領域である可能性が高い、という事実を確認することができた。本発明に係る日照変化領域検出装置は、上記条件を満たす領域を日照変化領域として出力することにより、日照変化領域の検出を行うものである。続いて、本発明に係る日照変化領域検出の基本原理を、図29を参照しながら、もう少し具体的に説明してみよう。図29に示すCASE3A,CASE3Bは、いずれも図28に示すCASE3の範疇に入るものである。すなわち、点ai-1は点piの類似範囲内回転楕円体Eの内部にあるが、点ui-1は類似範囲外回転楕円体Eの外部にある例である。この場合、総合判定部440は、J=1背景画素なる判定結果を出力することになる。CASE3AとCASE3Bとの相違は、点ui-1の位置である。すなわち、いずれの場合も、点ui-1が類似範囲外回転楕円体Eの外部に位置する点に変わりはないが、CASE3Aでは、点ui-1が類似境界回転楕円体Eの表面の近傍に位置するのに対して、CASE3Bでは、点ui-1が遠方に位置している。したがって、図にハッチングを施して示すように、類似境界の外側回転楕円体Eの外側に所定の近傍範囲ηを設定しておけば、点ui-1がこの近傍範囲η内にある場合をCASE3A、近傍範囲η外にある場合をCASE3Bと分けることができるなお、図29にハッチングを施して示す領域は、近傍範囲ηの断面部分であり、実際の近傍範囲ηは、回転楕円体Eの表面を取り巻く立体になる。本願発明者が屋外の様々な環境で撮影された多数の動画画像について解析を行ったところ、CASE3Aとの判定結果が得られた画素は、日照変化画素、すなわち、日照変化に起因した変動が生じている背景領域の画素である可能性が高いことが確認できた。ここで、日照変化とは、太陽光による被写体に対する照明状態の短期的な変動を意味し、雲の動き、樹木の揺れ、波の動き、雨滴や雪片などの瞬間的に消えるノイズに基づく照明状態の変動を指すものである。なお、昼から夜にかけての長期的な明るさ変動は、本願にいう日照変化ではない。結局、ここに示すモデルでは、CASE3Aの場合は日照変化画素、それ以外の場合CASE3B、およびCASE1,2,4の場合は非日照変化画素との判定を行えばよい。これが本発明に係る日照変化領域検出の基本原理である。このように、CASE3Aの場合に日照変化画素である可能性が高い理由は、次のように考えられる。まず、CASE3は、もともと第1判定部420による判定結果では類似○とされたが、第2判定部430による判定結果では非類似×とされたケースである。別言すれば、毎回更新が行われる平均画像と比較すると類似であるが、変化が緩慢な参照画像と比較すると非類似、と判定されたケースということになる。一方、日照変化は、太陽光による被写体に対する照明状態の短期的な変動であり、たとえば、雲の短期的な動きによって生じる変化である。このような、日照変化の影響を受けた画素である日照変化画素の動的ふるまいは、CASE3Aの特徴に合致する。すなわち、毎回更新が行われる平均画像は、短期的な日照変化に追従可能なため類似との判定結果が得られるが、変化が緩慢な参照画像は、これに追従できずに非類似となり、しかも、非類似であっても、類似に近い非類似となる、という特徴に合致する。もちろん、CASE3Aとの判定結果が得られた画素が日照変化画素である信憑性は、近傍範囲ηの設定内容に依存する。近傍範囲ηを狭く設定しすぎると、本来は日照変化画素であるのにCASE3Bと判定され、非日照変化画素との誤判定がなされることになる。逆に、近傍範囲ηを広く設定しすぎると、本来は非日照変化画素であるのにCASE3Aと判定され、日照変化画素との誤判定がなされることになる。そこで、実用上、日照変化領域の検出を高精度で行うためには、実際に用いる監視カメラについて得られる動画画像に基づいて、最適な近傍範囲ηを試行錯誤で決定するようにすればよい。なお、日照変化は、主として明るさに関する変化であり、画素に関して、色相よりも明度に重点をおいた変化を奏するものである。したがって、図29に示す回転楕円体Eを用いたモデルの場合、近傍範囲ηは、楕円の長軸方向図15に示すモデルにおけるZ軸方向への広がりが、楕円の短軸方向への広がりよりも大きくなるように設定するのが理にかなっている。<<<§11.本発明に係る日照変化領域検出装置>>>ここでは、本発明に係る日照変化領域検出装置の基本構成とその動作を説明する。図30は、本発明に係る日照変化領域検出装置の基本構成を示すブロック図である。前述したとおり、この日照変化領域検出装置は、§4で述べた改良型の移動物体検出装置を開発するプロセスにおける副産物として案出されたものであり、その構成要素の大部分は、図9に示す改良型の移動物体検出装置の構成要素と共通する。そこで、図30において、図9に示す構成要素と同一の構成要素については同符号を付して示す。図30に示すとおり、この装置は、画像入力部110、原画像格納部120、平均画像作成部130、参照画像作成部135、平均画像格納部140、参照画像格納部145、画像比較部400A、日照変化画像格納部510、画像出力部520、パラメータ設定部350によって構成され、動画画像について日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出装置として機能する。この図30に示す各構成要素のうち、画像入力部110、原画像格納部120、平均画像作成部130、平均画像格納部140、参照画像作成部135、参照画像格納部145の機能は、図9に同符号で示す各構成要素の機能と全く同じである。すなわち、画像入力部110は、時系列で連続的に与えられるフレーム単位のカラー原画像を、三原色の各画素値を有する画素の集合体データとして入力する機能を果たし、原画像格納部120は、入力された原画像を逐次格納する機能を果たす。また、平均画像作成部130は、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した平均画像とに基づいて、新たな平均画像を作成する機能を果たし、平均画像格納部140は、作成された平均画像を逐次格納する機能を果たす。一方、参照画像作成部135は、画像入力部110によって新たな原画像が入力されるたびに、直前に作成した参照画像に基づいて新たな参照画像を作成する機能を果たし、参照画像格納部145は、作成された参照画像を逐次格納する機能を果たす。なお、後述するように、条件判定部450Aは、図9に示す前背判定部450の機能をそっくり含んでおり、前景画素/背景画素のいずれかを示す判定結果Jを出力する機能を有する。したがって、参照画像作成部135は、図9に示す例と同様に、この判定結果Jに基づいて、参照画像の画素値を更新するか否かを決定することができる。また、画像比較部400Aは、図9に示す画像比較部400と同様に、原画像格納部120に格納された原画像Pと、平均画像格納部140に格納された平均画像Aと、参照画像格納部145に格納された参照画像Uと、を用いた比較処理を行う構成要素である。ただ、その主目的は、前景領域Fと背景領域Bとを区別するマスク画像Mを作成することではなく実際には、原画像Pを構成する個々の画素が前景画素/背景画素のいずれであるかを判定結果Jとして出力する機能を有しているため、マスク画像Mを作成することもできる、日照変化領域と非日照変化領域とを区別する日照変化画像Sを作成することである。図示のとおり、この画像比較部400Aは、画素値読出部410A、条件判定部450A、画素値書込部460Aによって構成されているが、これら各構成要素の機能については後に詳述する。図30に示す日照変化画像格納部510および画像出力部520の機能は、図9に示すマスク画像格納部310および画像出力部320の機能に対応するものである。ただ、格納および出力される画像が、マスク画像Mではなく、日照変化画像Sである点が異なる。日照変化画像Sは、原画像格納部120に格納されている原画像Pについての日照変化領域と非日照変化領域とを示す二値画像である図では、日照変化領域をハッチング領域として示すことにする。パラメータ設定部350は、この装置で利用されるパラメータを設定する構成要素である。具体的には、平均画像作成部130で用いられる重みw、参照画像作成部135で用いられる重みv、条件判定部450Aで用いられる判定パラメータたとえば、回転楕円体モデルを利用した類否判定を行う場合は、パラメータh,k、そして図29に示す近傍範囲ηを決定するためのパラメータを、パラメータ設定部350によって設定することができる。ユーザは、必要に応じて、パラメータ設定部350に対して操作入力を行い、これらパラメータの値を調整し、より精度の高い検出結果を得ることができる。もちろん、これらのパラメータの一部もしくは全部は、固定値にしておくことも可能であり、全部のパラメータ値を固定値にする実施例では、パラメータ設定部350を設ける必要はない。図30に示す日照変化領域検出装置が、図9に示す改良型の移動物体検出装置と異なる主要点は、画像比較部400の代わりに画像比較部400Aを設けた点である。図示のとおり、画像比較部400Aは、画素値読出部410A、条件判定部450A、画素値書込部460Aによって構成されており、条件判定部450Aは、第1判定部420A、第2判定部430A、総合判定部440Aによって構成されている。画素値読出部410Aの機能は、図9に示す画素値読出部410の機能と同じである。すなわち、画像入力部110が第i番目の原画像Piを入力し、これを原画像格納部120に格納するたびに、この原画像格納部120に格納されている第i番目の原画像Pi、平均画像格納部140に格納されている第i-1番目の平均画像Ai-1、参照画像格納部145に格納されている第i-1番目の参照画像Ui-1から、それぞれ対応する所定位置にある画素の画素値を読み出す処理を行う。ここでも、対応する所定座標x,yにある各画素の画素値をそれぞれpi,ai-1,ui-1で表すことにする。こうして画素値読出部410Aが読み出した、同一座標x,yに位置する画素の画素値pi,ai-1,ui-1は、条件判定部450Aに与えられる。なお、前述したとおり、実際には、各画素は、三原色R,G,Bの成分を有しているので、画素値pi,ai-1,ui-1は、それぞれ三原色R,G,Bの成分ごとに独立した値をもつことになる。上述したとおり、条件判定部450Aは、図9に示す前背判定部450の機能をそっくり含んでおり、前景画素/背景画素のいずれかを示す判定結果Jを出力する機能を有する。すなわち、条件判定部450Aは、画素値読出部410Aによって読み出された画素値pi,ai-1,ui-1に基づいて、第i番目の原画像Piの所定位置座標x,yの位置にある画素が前景画素か背景画素かを判定する処理を行い、判定結果Jを出力する。具体的には、判定結果がJ=0であれば前景画素、J=1であれば背景画素であることを示す。この判定結果Jは、参照画像作成部135による画素値の更新判断に利用される。もちろん、この判定結果Jを利用すれば、マスク画像Mを作成することも可能であるが、この装置の目的は、マスク画像Mを作成することではなく、日照変化画像Sを作成することであるため、この装置では、マスク画像Mの作成は行われない。条件判定部450Aには、図9に示す前背判定部450にはなかった付加機能が備わっている。すなわち、条件判定部450Aは、画素値読出部410Aによって読み出された画素値に基づいて、第i番目の原画像Piの所定位置にある画素piが前景画素か背景画素かを判定する機能だけでなく、日照変化画素か非日照変化画素かを判定する機能も備わっている。以下、この点を説明する。図示のとおり、条件判定部450Aは、第1判定部420A、第2判定部430A、総合判定部440Aによって構成されており、座標x,yに位置する画素についての判定を次のようなプロセスで実行する。まず、第1判定部420Aの機能は、図9に示す第1判定部420の機能と全く同じである。すなわち、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の平均画像Ai-1を構成する画素の画素値ai-1とを比較して類似するか否かを判定する。一方、第2判定部430Aは、図9に示す第2判定部430と同様に、第i番目の原画像Piを構成する画素の画素値piと第i-1番目の参照画像Ui-1を構成する画素の画素値ui-1とを比較して類似するか否かを判定する機能を有している。その上、非類似との判定がなされた場合には、更に、類似境界の所定近傍範囲内にあるか否かを判定する機能を有している。具体的には、図29にCASE3Aとして示すように、画素値ui-1に対応する点が予め設定された近傍範囲η内の点であるか否かを判定すればよい。また、総合判定部440Aは、第1判定部420Aによる判定結果と第2判定部430Aによる判定結果とに基づいて、判定結果JおよびNを出力する機能を有している。ここで、判定結果Jを出力する機能は、図9に示す総合判定部440の機能と全く同じである。すなわち、第1判定部420Aによる判定結果および第2判定部430Aによる判定結果の少なくとも一方が類似を示す場合には背景画素との判定を行い判定結果J=1を出力し、双方とも非類似を示す場合には前景画素との判定を行い判定結果J=0を出力する。一方、総合判定部440Aには、判定結果Nを出力する機能も備わっている。この機能は、第1判定部420Aにより類似との判定結果が得られ、かつ、第2判定部430Aにより非類似かつ所定近傍範囲内との判定結果が得られた場合には、日照変化画素との判定を行い判定結果N=1を出力し、その他の場合には、非日照変化画素との判定を行い判定結果N=0を出力する機能である。画素値書込部460Aは、第i番目の日照変化画像Siを構成する所定位置座標x,yの位置の画素の画素値として、条件判定部450Aが日照変化画素と判定した場合には日照変化領域を示す画素値この例では、判定結果Nと同じ論理値1を、条件判定部450Aが非日照変化画素と判定した場合には非日照変化領域を示す画素値この例では、判定結果Nと同じ論理値0を、それぞれ日照変化画像格納部510に書き込む処理を行う。以上、画像比較部400Aが、特定の座標x,yに位置する1画素について実行する具体的な処理を述べたが、このような処理が、1枚の画像を構成する全画素について繰り返し実行される。こうして、原画像格納部120に格納された第i番目の原画像Piを構成する全画素について、それぞれ日照変化画素か非日照変化画素かの判定が行われることになり、その判定結果に基づいて、日照変化画像格納部510内に第i番目i≧2の日照変化画像Siが得られることになる。上例の場合、こうして得られた日照変化画像Siにおいて、画素値1をもつ画素からなる領域が、日照変化領域、すなわち、日照変化に起因した変動が生じていると推定される領域ということになる。図31は、図30に示す日照変化領域検出装置における条件判定部450Aの構成を示すブロック図である。上述したとおり、第1判定部420Aからは、原画像Pと平均画像Aとを画素同士で比較した類否結果○/×が出力される。また、第2判定部430Aからは、原画像Pと参照画像Uとを画素同士で比較した類否結果○/×とともに、非類似×との結果が得られた場合には、類似境界の所定近傍範囲内/外を示す結果が出力される。総合判定部440Aは、これらの結果を受けて、判定結果NおよびJを出力することになる。図31の総合判定部440Aのブロック内には、このような判定処理の概念が表として示されている。ここで、表の横欄は、第1判定部420Aによる判定結果の別類似○/非類似×を示し、縦欄は、第2判定部430Aによる判定結果の別類似○/非類似×を示している。総合判定部440Aは、少なくとも一方が類似○を示す場合には背景画素J=1、双方とも非類似×を示す場合には前景画素J=0との判定を行うことになる。この点は、図10の下段に示す前背判定部450による判定処理と全く同様である。ただ、図31に示す総合判定部440Aの表は、図28に示すCASE3に対応する左下欄が近傍範囲内CASE3Aと近傍範囲外CASE3Bとに分かれている。すなわち、第1判定部420Aが類似○を示し、第2判定部430Aが非類似×を示している場合は、更に、第2判定部430Aから出力される近傍範囲内/外を示す情報に基づいて、CASE3AとCASE3Bとに分けて取り扱い、CASE3Aの場合にのみ日照変化画素N=1との判定結果を出力し、その他の場合は非日照変化画素N=0との判定結果を出力する。なお、近傍範囲内/外を示す情報は、第1判定部420Aが非類似×を示している場合には必要ないので、演算負担を軽減するためには、まず、第1判定部420Aによる判定を行い、その結果を第2判定部430Aに通知し、第2判定部430Aは、画素piとui-1とは非類似×との判定結果が得られ、かつ、第1判定部420Aから画素piとai-1とは類似○との判定結果が通知された場合に限り、近傍範囲内/外の判定を行うようにしてもよい。図32は、図30に示す日照変化領域検出装置によって作成される日照変化画像の一例を示す図である。上述したとおり、画像入力部110が、第i番目の原画像Piを入力するたびに、画像比較部400Aの比較処理によって、この原画像Piに対応した第i番目の日照変化画像Siが出力されることになる但し、第1番目の原画像P1に対応する日照変化画像S1は作成できないので、日照変化画像Siについては、i≧2となる。こうして得られた日照変化画像Siは、画素値0判定結果N=0が得られた画素もしくは画素値1判定結果N=1が得られた画素をもつ画素の集合からなる二値画像であり、画素値1をもつ画素から構成される領域は日照変化領域、画素値0をもつ画素から構成される領域は非日照変化領域を示す。すなわち、日照変化画像Siは、原画像Piについて、どの部分が日照変化に起因した変動が生じている領域かを示す画像ということになる。このように、動画画像について日照変化に起因した変動が生じている領域を把握することは、当該動画画像の特徴を把握する上で有用である。たとえば、図32に示すような原画像Piに日照変化画像Siを重ねてディスプレイ画面上に表示すれば、路面や家屋の屋根の部分が日照変化の影響を大きく受けている領域であるといった事項を視覚的に認識することができるようになる。なお、実用上、この図30に示す日照変化領域検出装置も、コンピュータに専用のプログラムを組み込むことにより構築することができ、図30に各ブロックで示す構成要素は、コンピュータのハードウエアとソフトウエアとの組み合わせによって実現することができる。もちろん、この移動物体検出装置を半導体集積回路として構成することも可能である。<<<§12.本発明に係る日照変化領域検出装置の実用的な実施形態>>>§11で述べたとおり、本発明に係る日照変化領域検出装置は、§4で述べた改良型の移動物体検出装置の構成の一部を改変することにより実現することができる。したがって、この日照変化領域検出装置においても、§7~§9で述べた回転楕円体モデルを利用した類否判定を行い、演算負担の軽減を図ることが可能である。ただ、日照変化領域検出装置を構成する第2判定部430Aには、付加的な判定処理が必要になる。すなわち、第2判定部430Aは、比較対象となる一対の画素について、類似するか否かの判定を行うとともに、非類似との判定がなされた場合には、更に、類似境界の所定近傍範囲内にあるか否かの判定を行う必要がある。前者の判定には、§8もしくは§9で述べた演算方法を適用することにより、演算負担の軽減が可能になる。ここでは後者の判定について、演算負担を軽減する方法を説明する。後者の判定は、図29に示すように、点ui-1が近傍範囲η内の点であるか否かの判定であり、当該判定を行うためには、当然ながら、予め近傍範囲ηを幾何学的に定義しておく必要がある。図33は、RGB三次元座標系上における近傍範囲ηの設定例を示す図である。図示の例では、図15に示すモデルと同様に、基準点Qを中心として配置された回転楕円体が基準立体Eとして定義されている。ここで、画素ui-1に対応する比較点qが、この回転楕円体Eの内部の点であれば類似、外部の点であれば非類似と判定されることになる。そして、非類似と判定された場合には、更に、比較点qが近傍範囲η内にあるか否かの判定が行われる。近傍範囲ηとしては、たとえば、基準立体Eと同様に基準点Qを中心として配置された基準立体Eよりも大きなサイズをもつ拡張基準立体E′を定義し、基準立体Eの外側かつ拡張基準立体E′の内側の領域として定義することができる。この場合、第2判定部430Aは、基準立体Eおよび拡張基準立体E′と比較点qとの位置関係を調べ、比較点qが基準立体Eの外部、かつ、拡張基準立体E′の内部に位置すると判定できる場合に、類似境界の所定近傍範囲内にあるとの判定を行うことになる。図示の例のように、基準立体Eおよび拡張基準立体E′の双方を回転楕円体とするのであれば、たとえば、回転楕円体Eと回転楕円体E′とを相似形の立体として定義すればよい。前述したとおり、回転楕円体Eは、三次元座標系の原点Oと基準点Qとの距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この楕円を基準軸Zを中心に回転させることにより得られる回転体として定義することができる。したがって、回転楕円体EのサイズをK倍K>1した立体として回転楕円体E′を定義するのであれば、回転楕円体E′は、値Kαを長軸半径とし、値Kβを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この楕円を基準軸Zを中心に回転させることにより得られる回転体として定義することができる。図33に示す例は、このような方法で、回転楕円体EおよびE′を定義した例であり、図にハッチングを施した部分が近傍範囲ηということになる。しかしながら、図33に例示する近傍範囲ηは、必ずしも効率的な設定例ではない。この例のような近傍範囲ηを設定した場合、比較点qが近傍範囲η内の点であるか否かを判定するためには、比較点qが回転楕円体Eの外側の点であり、かつ、回転楕円体E′の内側の点であることを確認する必要があり、回転楕円体Eに対する幾何学演算と、回転楕円体E′に対する幾何学演算との双方が必要になる。そこで、ここでは、基準立体として回転楕円体Eを採用し、§8もしくは§9で述べた演算方法を適用して類否判定を行うことを前提とした場合に、極めて効率的な近傍範囲ηの設定例を、図34を参照して説明する。この図34は、図19に対応する図であり、同一符号で示す構成要素は同一の構成要素である。ここに示す例の場合、ハッチングを施して示す領域が、近傍範囲ηとして設定される。もっとも、図34は断面図であるため、ハッチングを施して示す領域は、この近傍範囲ηの断面を示しているにすぎない。実際の近傍範囲ηは、このハッチングが施された平面領域を、基準軸Zまわりに回転することによって得られる回転体からなる領域である。この図34に示す近傍範囲ηは、回転楕円体Eの外側近傍に位置する領域であり、日照変化画素か否かの判定に利用する近傍範囲に相応しい領域である。特に、回転楕円体Eの内部の点を基準軸Z方向に移動させると、この近傍範囲η内に入る関係になっている。これは、日照変化に基づく明度変化が生じた場合に、回転楕円体Eの内部の点が近傍範囲η内に移動してくる可能性が高いことを意味しており、このような観点からも、日照変化画素か否かの判定に利用するのに適している。実際、本願発明者が行った実験では、図34に示すような近傍範囲ηを設定した場合、様々な環境で撮影された動画画像に対して、良好な日照変化領域の検出ができた。しかも、図34に示すような近傍範囲ηを設定した場合、比較点qが近傍範囲η内の点であるか否かを判定するための特別な演算が一切不要になるという利点が得られる。その理由は、§8で述べた類否判定手順を参照すれば容易に理解できる。すなわち、図25の流れ図では、ステップS13,S15,S17において、それぞれ所定の演算が行われ、ステップS14,S16,S18において、演算結果が所定の条件を満足しているか否かの判定が行われている。ここで、ステップS14の判定は、比較点qが図17のハッチング領域内にあるか否かの判定であり、ステップS16の判定は、比較点qが図19のハッチング領域内にあるか否かの判定であり、ステップS18の判定は、比較点qが図23のハッチング領域内にあるか否かの判定である。ここで、比較点qが図34のハッチング領域内すなわち、近傍範囲η内にあった場合に、図25の流れ図でどのような判定がなされるかを考えてみる。すると、ステップS14では否定的判定がなされ、ステップS16では否定的判定がなされ、ステップS18では肯定的判定がなされることがわかる。結局、図25の流れ図を実行した結果、ステップS14で否定的判定、ステップS16で否定的判定、ステップS18で肯定的判定がなされたとすれば、その時点で、比較点qが近傍範囲η内にあると判定することができることになる。これは、第2判定部430Aが、図25の流れ図に従った類否判定手順のみを実行すれば足り、その結果、非類似との判定結果が得られたとしても、比較点qが近傍範囲η内にあるか否かの付加的な判定手順を行う必要がないことを意味する。これは、演算負担を軽減する上で大きなメリットである。図27には、図9に示す移動物体検出装置に用いるための第1判定部420第2判定部430も同様の具体的な構成例を示した。この図27に示す構成例の場合、判定値出力部425は、類似○論理値1か非類似×論理値0かを示す判定結果しか出力することができない。これに対して、図31に示す第2判定部430Aは、非類似×の場合には、更に、近傍範囲内か近傍範囲外かを示す判定結果を出力する機能が必要になる。図34に示すような近傍範囲ηを設定すれば、図27に示す判定値出力部425に、非常に単純な機能を付加するだけで、近傍範囲内か近傍範囲外かを示す判定結果を出力させることが可能になり、図31に示す第2判定部430Aとして利用可能になる。具体的には、図27に示す構成において、判定値出力部425が、第1の条件判定部422および第2の条件判定部423が否定的判定を行い、第3の条件判定部424が肯定的判定を行った場合に、類似境界の所定近傍範囲η内との判定結果を示す判定値を出力し、その他の場合には、所定近傍範囲η外との判定結果を示す判定値を出力するような機能を付加すればよい。そうすれば、何ら余分な演算要素を付加することなく、図27に示す第1判定部420すなわち、第2判定部430を図31に示す第2判定部430Aとして利用することができる。<<<§13.本発明に係る移動物体検出装置>>>ここでは、本発明に係る移動物体検出装置の基本構成およびその動作を説明する。ここで述べる移動物体検出装置の特徴は、それぞれ異なる判定基準を用いて移動物体を検出する複数n台の移動物体検出ユニットを用意し、日照変化の状態に応じて、その時点で最適と推定されるユニットの検出結果を選択して出力する点にある。したがって、この移動物体検出装置では、日照変化の状態を検出する装置が必須の構成要素になる。図35は、n=2に設定し、2台の移動物体検出ユニットを用いて構成した移動物体検出装置の基本構成を示すブロック図である。図示のとおり、この移動物体検出装置は、第1の移動物体検出ユニット610、第2の移動物体検出ユニット620、選択ユニット710、日照変化率算出ユニット720、日照変化領域検出ユニット730によって構成されており、動画画像について移動物体を検出する機能を有する。第1の移動物体検出ユニット610、第2の移動物体検出ユニット620および日照変化領域検出ユニット730には、移動物体の検出対象となる共通の動画画像Hが与えられる。この共通動画画像Hは、これまで述べてきた動画画像と同様に、時系列で連続的に与えられるフレーム単位の原画像によって構成される。ここで重要な点は、2台の移動物体検出ユニット610,620が、それぞれ異なる判定基準を用いて移動物体の検出を行う点である。個々の移動物体検出ユニット610,620としては、時系列で連続的に与えられるフレーム単位の原画像に基づいて、何らかの判定基準を参照して移動物体の検出を行う機能をもった装置であれば、どのような装置を用いてもかまわない。一般的には、入力された原画像を、移動物体を含まないと推定される何らかの背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出する装置を用いればよい。具体的には、§1~§3に従来装置として述べた移動物体検出装置を用いることもできるし、§4~§9に改良型の移動物体検出装置として述べた装置を用いることもできる。実用上は、2台の移動物体検出ユニット610,620として、同じアルゴリズムで移動物体の検出を行う装置を用い、当該アルゴリズムで利用される判定基準のみを両者で異ならせるようにすればよい。たとえば、2台の移動物体検出ユニット610,620として、§3で述べた従来型の移動物体検出装置図8に示す装置を用いた場合、類似判定基準は、類否判定部220内の判定パラメータによって設定されるので、予めパラメータ設定部330を用いて、第1の移動物体検出ユニット610として利用される装置については第1の類似判定基準を設定し、第2の移動物体検出ユニット620として利用される装置については、第1の類似判定基準とは異なる第2の類似判定基準を設定しておけばよい。もちろん、2台の移動物体検出ユニット610,620として、§4で述べた改良型の移動物体検出装置図9に示す装置を用いることもできる。この場合も、予めパラメータ設定部340を用いて、それぞれ異なる類似判定基準を設定しておけばよい。あるいは、本願で述べていない別な移動物体検出装置を移動物体検出ユニットとして利用してもかまわない。図8に示す従来型装置の場合も、図9に示す改良型装置の場合も、具体的な類似判定基準は基準立体回転楕円体の位置およびサイズである。ここで、基準立体のサイズは、類似判定基準のレベルを調整するパラメータになる。すなわち、原画像と背景画像図8や図9に示す例の場合、平均画像や参照画像との類否判断を行う場合、両画像を個々の画素ごとに比較し、比較対象となる一方の画素の画素値を基準画素値、他方の画素の画素値を比較画素値として、三原色の各画素値を各座標軸にとった三次元座標系において、基準画素値に対応する座標に位置する基準点Qと比較画素値に対応する座標に位置する比較点qをとり、基準点Qを中心とする所定サイズの基準立体と比較点qとの位置関係を調べ、比較点qが基準立体の外部に位置すると判定できる場合に非類似との判定を行い、内部に位置すると判定できる場合には類似との判定を行うことになる。結局、基準立体のサイズが大きくなれば、類似判定基準は緩くなり、類似と判定される確率が高くなるが、基準立体のサイズが小さくなれば、類似判定基準は厳しくなり、類似と判定される確率は低くなる。したがって、第1の類似判定基準として設定される基準立体のサイズと、第2の類似判定基準として設定される基準立体のサイズとが異なるようにしておけば、2台の移動物体検出ユニット610,620には、互いに異なる類似判定基準を設定することができる。ここで、2台の移動物体検出ユニット610,620は、同時に並列して動作するので、共通動画画像Hに基づく移動物体の検出結果として、第1の移動物体検出ユニット610からはマスク画像M1が出力され、第2の移動物体検出ユニット620からはマスク画像M2が出力される。ただ、両者では、類似判定基準が異なっているので、マスク画像M1およびM2は必ずしも同一にはならない。一方、日照変化領域検出ユニット730は、与えられた共通動画画像Hについて、日照変化に起因した変動が生じている日照変化領域を随時検出する機能を果たす。この日照変化領域検出ユニット730としては、たとえば、§11で述べた日照変化領域検出装置図30に示す装置をそのまま利用することができる。この場合、共通動画画像Hを構成する1フレーム分の画像を、図30に示す原画像Pとして与えることにより、個々のフレーム単位で日照変化画像Sが得られることになる。もちろん、日照変化領域検出ユニット730としては、何らかの方法で日照変化領域を検出する機能をもった装置であれば、どのような装置を利用してもかまわない。日照変化率算出ユニット720は、こうして日照変化領域検出ユニット730によって検出された日照変化領域の全領域に対する割合を示す日照変化率ρを随時算出する構成要素である。すなわち、日照変化領域検出ユニット720からは、図32の右に示されているような日照変化画像Siが逐次与えられるので、日照変化率算出ユニット720は、この画像Siの全領域に対する日照変化領域図では、ハッチングを施して示す領域の割合を算出して、これを日照変化率ρとして出力する処理を行うことになる。具体的には、図示の例の場合、画素値N=1をもった画素数をカウントすれば、日照変化率ρ=カウント値/全画素数として算出できる。選択ユニット710は、この日照変化率ρに基づいて、2台の移動物体検出ユニット610,620のうちのいずれか1台を随時選択し、選択された移動物体検出ユニットが検出した移動物体を最終的な検出結果として出力する。各移動物体検出ユニット610,620には、それぞれ互いに異なる類似判定基準が設定されているが、選択ユニット710は、日照変化率ρが高いほど、より緩い基準が設定されたユニットを選択する。たとえば、図35に示す実施形態において、第1の移動物体検出ユニット610が第1の類似判定基準に基づく判定を行い、第2の移動物体検出ユニット620が、この第1の類似判定基準よりも緩い第2の類似判定基準に基づく判定を行うような設定がなされていた場合を考えてみよう。具体的には、第1の移動物体検出ユニット610が類否判定に用いる基準立体のサイズよりも、第2の移動物体検出ユニット620が類否判定に用いる基準立体のサイズの方が大きく設定されていた場合を考えればよい。この場合、選択ユニット710は、日照変化率算出ユニット720によって算出された日照変化率ρが、所定のしきい値以下の場合には第1の移動物体検出ユニット610を選択し、当該しきい値を超える場合には第2の移動物体検出ユニット620を選択する。図示の例の場合、選択ユニット710は、第1の移動物体検出ユニット610が選択されたときには、第1の移動物体検出ユニット610が作成したマスク画像M1を最終的なマスク画像M移動物体の検出結果として出力し、第2の移動物体検出ユニット620が選択されたときには、第2の移動物体検出ユニット620が作成したマスク画像M2を最終的なマスク画像M移動物体の検出結果として出力することになる。もちろん、日照変化領域検出ユニット730が出力する日照変化画像Sの内容は時間とともに共通動画画像Hとして与えられるフレーム単位の原画像ごとに変化してゆくため、日照変化率算出ユニット720で算出される日照変化率ρも時間とともに変化してゆく。このため、選択ユニット710による選択結果も時間とともに変化し、選択ユニット710からは、ある期間は第1の移動物体検出ユニット610が作成したマスク画像M1が出力され、別なある期間は第2の移動物体検出ユニット620が作成したマスク画像M2が出力されることになる。結局、選択ユニット710が行う選択処理は、日照変化率ρが高いほど、より緩い類似判定基準を用いた検出結果を出力する処理ということになり、理にかなった処理になる。すなわち、日照変化率ρは、その時点における日照変化に起因した画像変動の程度を示すパラメータであり、この日照変化率ρに基づいて、その時点で最適と推定されるユニットの検出結果を選択し、これを最終的な検出結果として出力するようにすれば、日照変化が生じていた場合でも、当該日照変化の程度に応じた最適な判定基準に基づく検出結果を得ることができるようになり、より正確な移動物体の検出が可能になる。より具体的に説明すれば、上例の場合、第1の移動物体検出ユニット610には厳しい類似判定基準が設定され、第2の移動物体検出ユニット620には緩い類似判定基準が設定されている。ここで、日照変化率ρが低い期間たとえば、濃霧により、全体的にコントラストの低い画像が得られている期間は、第1の移動物体検出ユニット610の検出結果の方がより適切な結果になると推定される。なぜなら、日照変化の影響で色が大きく変化する可能性は低いので、類似判定基準を厳しくしても、日照変化の影響で色が変化した背景領域が移動物体の領域前景領域として誤検出される可能性も低くなるためである。このような環境下で、類似判定基準を緩くすると、移動物体の領域前景領域が背景領域として誤検出される可能性が高まり、好ましくない。したがって、日照変化率ρが低い期間は、厳しい類似判定基準を採用した第1の移動物体検出ユニット610の検出結果を選択して出力した方が、より正確な検出が期待できる。一方、日照変化率ρが高い期間たとえば、強い日差しがあり、全体的にコントラストの高い画像が得られている期間は、第2の移動物体検出ユニット620の検出結果の方がより適切な結果になると推定される。なぜなら、日照変化の影響で色が大きく変化する可能性が高くなるので、類似判定基準を厳しくすると、日照変化の影響で色が変化した背景領域が移動物体の領域前景領域として誤検出される可能性が高まり、好ましくない。類似判定基準を緩くすれば、このような誤検出の可能性を低減させることができる。したがって、日照変化率ρが高い期間は、緩い類似判定基準を採用した第2の移動物体検出ユニット620の検出結果を選択して出力した方が、より正確な検出が期待できる。具体的には、選択ユニット710は、日照変化率ρのしきい値を、たとえば3と定めておき、ρ≦3の場合には第1の移動物体検出ユニット610を選択し、ρ>3の場合には第2の移動物体検出ユニット620を選択する処理を行えばよい。しきい値の設定は、監視カメラの実際の設置環境に応じて、より好ましい結果が得られるように適宜調整するようにすればよい。以上、2台の移動物体検出ユニット610,620を用いる実施形態を示したが、本発明で用いる移動物体検出ユニットの台数は必ずしも2台に限定されるものではなく、それぞれ異なる判定基準を用いて移動物体を検出する機能を有する複数n台の移動物体検出ユニットを用いてもかまわない。この場合は、厳しい/緩いという2段階の類似判定基準ではなく、よりきめの細かな合計n段階の類似判定基準が設定されることになり、選択ユニット710は、日照変化率ρに基づいて、この複数n台の移動物体検出ユニットのうちのいずれか1台を随時選択し、選択された移動物体検出ユニットが検出した移動物体を最終的な検出結果として出力することになる。具体的には、色空間上に配置した基準立体を類否判定に用いる場合、複数n台の移動物体検出ユニットが採用する基準立体のサイズがそれぞれ異なるような設定を行えばよい。すなわち、第j番目2≦j≦nの移動物体検出ユニットが採用する基準立体が、第j-1番目の移動物体検出ユニットが採用する基準立体を包含する、より大きなサイズとなっているようにすればよい。この場合、選択ユニット710は、日照変化率ρが高いほど、よりサイズの大きな基準立体を採用するユニットを選択することになる。基準立体として、回転楕円体を用いる場合、各移動物体検出ユニットは、三次元座標系の原点Oと基準点Qとの距離Dに、所定のパラメータh但し、h<1を乗じて求まる値α=hDを長軸半径とし、所定のパラメータk但し、k<hを乗じて求まる値β=kDを短軸半径とする楕円を、原点Oと基準点Qとを結ぶ基準軸Z上に長軸が重なるように、かつ、基準点Qが中心点となるように配置し、この楕円を前記基準軸Zを中心に回転させることにより得られる回転楕円体Eもしくはその近似体を基準立体として用いることになる。この場合、第j番目1≦j≦nの移動物体検出ユニットが採用する基準立体のパラメータh,kをそれぞれhj,kjとしたときに、hj>hj-1,kj>kj-1に設定すればよい但し、2≦j≦n。この図35に示す移動物体検出装置も、実用上は、コンピュータに専用のプログラムを組み込むことにより構築することができ、図35に各ブロックで示す構成要素は、コンピュータのハードウエアとソフトウエアとの組み合わせによって実現することができる。もちろん、この移動物体検出装置を半導体集積回路として構成することも可能である。<<<§14.本発明に係る日照変化領域検出および移動物体検出方法>>>最後に、本発明を方法の発明として捉えた場合に、当該方法を構成する個々の手順を説明する。まず、本発明に係る日照変化領域検出方法は、動画画像について日照変化に起因した変動が生じている日照変化領域を検出する方法であり、実用上、この方法を構成する各手順はコンピュータによって実行される。すなわち、この方法は、コンピュータが、時系列で連続的に与えられるフレーム単位のカラー原画像を、三原色の各画素値を有する画素の集合体データとして入力する画像入力段階と、コンピュータが、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した平均画像とに基づいて、新たな平均画像を作成する平均画像作成段階と、コンピュータが、新たな原画像が入力されるたびに、当該新たな原画像と直前に作成した参照画像とに基づいて、新たな参照画像を作成する参照画像作成段階と、コンピュータが、原画像と平均画像と参照画像とを用いた比較処理により、日照変化領域と非日照変化領域とを区別する日照変化画像を作成する画像比較段階と、コンピュータが、日照変化画像を出力する画像出力段階と、によって構成される。ここで、画像比較段階は、第i番目の原画像、第i-1番目の平均画像、第i-1番目の参照画像から、それぞれ対応する所定位置にある画素の三原色の各画素値を読み出す画素値読出ステップを有している。また、画像比較段階は、第i番目の原画像を構成する画素の画素値piと第i-1番目の平均画像を構成する画素の画素値ai-1とが類似するか否かを判定する第1判定処理と、第i番目の原画像を構成する画素の画素値piと第i-1番目の参照画像を構成する画素の画素値ui-1とが類似するか否かを判定するとともに、非類似との判定がなされた場合には、更に、類似境界の所定近傍範囲内にあるか否かを判定する第2判定処理と、を行い、第i番目の原画像の前記所定位置にある画素について、第1判定処理により類似との判定結果が得られ、かつ、第2判定処理により非類似かつ所定近傍範囲内との判定結果が得られた場合には、日照変化画素、その他の場合には、非日照変化画素との判定を行うとともに、第1判定処理による判定結果と第2判定処理による判定結果との少なくとも一方が類似を示す場合には背景画素、双方とも非類似を示す場合には前景画素との判定を行う条件判定ステップを有している。更に、画像比較段階は、第i番目の日照変化画像を構成する前記所定位置の画素の画素値として、条件判定ステップにより日照変化画素と判定された場合には日照変化領域を示す画素値を、条件判定ステップにより非日照変化画素と判定された場合には非日照変化領域を示す画素値を、それぞれ与えることにより、第i番目の日照変化画像を作成する日照変化画像作成ステップを有している。また、平均画像作成段階では、第1番目の原画像をそのまま第1番目の平均画像とし、第i番目i≧2の平均画像を作成する際には、第i番目の原画像と第i-1番目の平均画像とについて、互いに対応する位置にある画素の三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を各画素値とする画素の集合からなる画像を第i番目の平均画像とする処理が行われる。一方、参照画像作成段階では、第1番目の原画像をそのまま第1番目の参照画像とし、第i番目i≧2の参照画像を作成する際には、第i-1番目の参照画像を構成する個々の画素のうち、条件判定ステップにおいて第i番目の原画像の対応画素が前景画素と判定された画素については、更新せずに、もとの画素値をそのまま維持し、条件判定ステップにおいて第i番目の原画像の対応画素が背景画素と判定された画素については、第i-1番目の平均画像の対応する位置にある画素との間で、三原色の各画素値の重みつき平均値をそれぞれ求め、得られた各平均値を新画素値として更新することにより、第i番目の参照画像を作成する処理が行われる。なお、実用上は、平均画像作成段階において、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の平均画像A1とする処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の平均画像Aiを、ai=1-wai-1wpi但し、aiは、平均画像Aiの所定位置の画素の所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、piは、原画像Piの前記所定位置の画素の前記所定色の画素値、wは、所定の重みを示すパラメータw<1なる演算式を用いて作成するようにする。また、参照画像作成段階では、最初の原画像P1が入力されたときに、当該原画像P1をそのまま最初の参照画像U1とする処理を行い、以後、第i番目の原画像Piが入力されるたびにi=2,3,、第i番目の参照画像Uiの各画素値を、背景画素と判定された画素については、ui=1-vui-1vai-1但し、uiは、平均画像Uiの所定位置の画素の所定色の画素値、ui-1は、平均画像Ui-1の前記所定位置の画素の前記所定色の画素値、ai-1は、平均画像Ai-1の前記所定位置の画素の前記所定色の画素値、vは、所定の重みを示すパラメータv<1、かつw>vなる演算式を用いて求められた画素値uiに更新し、前景画素と判定された画素については、ui=ui-1として、更新せずに、もとの画素値をそのまま維持することにより、参照画像Uiを作成するようにする。一方、本発明に係る移動物体検出方法は、動画画像について移動物体を検出する方法であり、実用上、この方法を構成する各手順はコンピュータによって実行される。すなわち、この方法は、コンピュータが、時系列で連続的に与えられるフレーム単位の原画像によって構成される動画画像について、入力された原画像を背景画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出する移動物体検出方法において、動画画像について、日照変化に起因した変動が生じている日照変化領域を検出し、検出した日照変化領域の全領域に対する割合を示す日照変化率を算出し、この日照変化率が高いほど、より緩い類似判定基準を用いた検出を行うことを特徴とする。より具体的には、本発明に係る移動物体検出方法は、コンピュータが、時系列で連続的に与えられるフレーム単位の原画像によって構成される共通動画画像について、複数n通りの異なる判定基準を用いて複数n通りの移動物体検出結果を得る個別検出段階と、コンピュータが、共通動画画像について、日照変化に起因した変動が生じている日照変化領域を検出する日照変化領域検出段階と、コンピュータが、日照変化領域検出段階によって検出された日照変化領域の全領域に対する割合を示す日照変化率を算出する日照変化率算出段階と、コンピュータが、日照変化率に基づいて、複数n通りの移動物体検出結果のいずれか1つを選択し、選択された検出結果を最終的な検出結果として出力する検出結果出力段階と、を有する。ここで、個別検出段階では、入力された原画像を、移動物体を含まないと推定される背景画像§4で述べた改良型の移動物体検出装置を利用して個別検出段階を行う場合は、平均画像および参照画像と比較し、両者の類似度が所定の類似判定基準に満たないと判定された非類似領域を移動物体の領域として検出する処理を、互いに異なる複数n通りの類似判定基準を用いて実行し、検出結果出力段階では、日照変化率が高いほど、より緩い類似判定基準を用いて得られた検出結果を選択すればよい。"