"広角レンズを有する撮像手段を設け、各実空間座標位置に予め設定された3次元寸法の仮想ブロックを想定し、この仮想ブロックが画像上に表示されるサイズを記憶する記憶手段と、この記憶手段に記憶された表示サイズから一つの画素の周辺画素を一群のデータとして画像処理を行うためのテンプレートサイズを設定するテンプレートサイズ設定手段と、このテンプレートサイズ設定手段により設定されたテンプレートサイズを用いて前記撮像手段により撮像された画像を処理する画像処理手段とを備える制御手段を設けたことを特徴とする車両周辺監視装置。前記記憶手段に記憶される表示サイズは、画像上の画素位置と実空間座標と対応させた計測マップを用いて算出され、計測マップ無効範囲では無効範囲と有効範囲との境界座標の表示サイズデータで補完することを特徴とする請求項1に記載の車両周辺監視装置。前記記憶手段に記憶される表示サイズは、前記仮想ブロックが画像上に表示されるサイズと同等サイズの近似円盤の半径とすることを特徴とする請求項1又は請求項2に記載の車両周辺監視装置。前記記憶手段に記憶される表示サイズは、前記仮想ブロックが画像上に表示されるサイズと同等サイズの近似円盤の直径とし、この直径が奇数の場合にはテンプレートを方形とし、前記直径が偶数の場合にはテンプレートを円形とすることを特徴とする請求項1又は請求項2に記載の車両周辺監視装置。前記画像処理手段は、前記テンプレートサイズ設定手段により設定されたテンプレートサイズに基づいて画像フィルタリングを実行することを特徴とする請求項1~請求項4のいずれか1項に記載の車両周辺監視装置。ところで、従来、駐車支援装置の撮像手段のような歪みをもつ広視野の撮像手段においては、実空間位置画像平面上の位置関係を構成する方法、及び校正結果をまとめた計測マップを使用するものがある。しかし、一般的に、画像処理の際には、輝度変化の影響の他、ノイズによる画質への影響、路面や対象物の面肌等の不要なテクスチャを除去するために、対象画素周辺のデータを平均化するデジタルフィルタ処理が使用される。また、近年、コンピュータ性能の向上により、対象物の検出に際し、同様に、周辺の画素からなるデータ群の分散等を基にした技術HARRISオペレータ特徴量検出法、KLT法等を使用する機会が増してきている。このような処理の場合、一般に、周辺画素を1群とするデータ以下テンプレートというの取り扱いは、画面に対して一様に2次元画素寸法で定義されるため、屋外で使用する撮像手段等と対象との間隔車間距離が多々異なる場合においては、必要な情報を取り逃すおそれがある。よって、最適なテンプレートサイズの決定に関しては、距離に基づく指標及び3次元の実寸に基づく定義が必要であると考えられる。しかしながら、特に、後方駐車支援の撮像手段等に広角レンズのレンズ収差による画像歪み等の影響が存在するため、距離に関する基準の取得が困難となり、テンプレートサイズの定義には、さらに検討の余地があった。そこで、この発明の目的は、実寸と相関性の高いテンプレートサイズを設定することができ、その結果、良好な画像処理結果を得ることができ、また、不必要に大きいテンプレートサイズを設定することがないため、計算量の増加を防ぐことができる車両周辺監視装置を提供することにある。この発明は、広角レンズを有する撮像手段を設け、各実空間座標位置に予め設定された3次元寸法の仮想ブロックを想定し、この仮想ブロックが画像上に表示されるサイズを記憶する記憶手段と、この記憶手段に記憶された表示サイズから一つの画素の周辺画素を一群のデータとして画像処理を行うためのテンプレートサイズを設定するテンプレートサイズ設定手段と、このテンプレートサイズ設定手段により設定されたテンプレートサイズを用いて前記撮像手段により撮像された画像を処理する画像処理手段とを備える制御手段を設けたことを特徴とする。この発明の車両周辺監視装置は、実寸と相関性の高いテンプレートサイズを設定して良好な画像処理結果を得ることができ、また、不必要に大きいテンプレートサイズを設定することがないため、計算量の増加を防ぐことができる。この発明は、良好な画像処理結果を得ること、また、計算量の増加を防ぐ目的を、実寸と相関性の高いテンプレートサイズを設定したり、不必要に大きいテンプレートサイズを設定することをなくするようにして実現するものである。図1~図26は、この発明の実施例を示すものである。図1において、1は車両に搭載される車両周辺監視装置である。この車両周辺監視装置1は、広角レンズ2を有して車両後方を撮像する撮像手段3と、入力側で撮像手段3に連絡した制御手段4と、この制御手段4の出力側に連絡した表示手段5とを備える。撮像手段3は、車両に取り付けたカメラ等の機器からなり、車両周辺を撮像する。表示手段5は、車両の車室内に配置されたカーナビゲーションシステムのモニタ等の表示機器からなり、撮像手段3により撮像された画像等を表示する。制御手段4は、各実空間座標位置に予め設定された3次元寸法の仮想ブロックを想定し、この仮想ブロックが画像上に表示されるサイズを記憶する記憶手段6と、この記憶手段6に記憶された表示サイズから一つの画素の周辺画素を一群のデータとして画像処理を行うためのテンプレートサイズを設定するテンプレートサイズ設定手段7と、このテンプレートサイズ設定手段7により設定されたテンプレートサイズを用いて撮像手段3により撮像された画像を処理する画像処理手段8とを備える。これにより、実寸と相関性の高いテンプレートサイズを設定することができ、その結果、良好な画像処理結果を得ることができ、また、不必要に大きいテンプレートサイズを設定することがないため、計算量の増加を防ぐことができる。記憶手段6に記憶される表示サイズは、画像上の画像位置と実空間座標と対応させた計測マップを用いて算出され、計測マップ無効範囲では無効範囲と有効範囲との境界座標の表示サイズデータで補完する。これにより、画像上の全範囲にわたって、実寸と相関性の高いテンプレートサイズを設定することができ、その結果、良好な画像処理結果を得ることができ、また、不必要に大きいテンプレートサイズを設定することがないため、計算量の増加を防ぐことができる。また、記憶手段6に記憶される表示サイズは、前記仮想ブロックが画像上に表示されるサイズと同等サイズの近似円盤の半径とする。これにより、表示サイズを記憶するのにデータ量を小さくできるため、記憶手段6のメモリ容量を削減することができ、計算効率を向上させることができる。更に、記憶手段6に記憶される表示サイズは、前記仮想ブロックが画像上に表示されるサイズと同等サイズの近似円盤の直径とし、この直径が奇数の場合にはテンプレートを方形とし、前記直径が偶数の場合にはテンプレートを円形とする。これにより、表示サイズを近似円盤の直径とし、直径が奇数の場合と偶数の場合とでテンプレート形状を変えるため、表示サイズを近似円盤の半径とする場合よりも画像処理の精度を上げることができる。画像処理手段8は、テンプレートサイズ設定手段7により設定されたテンプレートサイズに基づいて画像フィルタリングを実行する。これは、特に、画像フィルタリングに有効である。この実施例においては、以下に述べるように、実空間上に想定した仮想ブロックにより実際の3次元スケールとその画像平面上での表示サイズに対応したテンプレート処理を行う。先ず、この発明に係る経緯内容を述べるために、その背景について説明する。一般的に画像処理の際には、輝度変化の影響の他、ノイズによる画質への影響、路面や対象物の面肌等の不要なテクスチャを除去するために、対象の画素周辺のデータを使用し、平均化処理等を行うデジタルフィルタ処理が使用される。また、近年のコンピュータ等の演算性能の向上により、対象物の検出に際し、同様に周辺画素からなるデータ群の分散等を基にした技術HARRISオペレータ特徴量検出法、KLT法等を使用する場合にもテンプレートを単位とする処理が行われる。例として、デジタルフルィルタとして、最も一般的である平均値フィルタを取り上げる。この平均値フィルタにおいては、入力画像の注目画素について、その位置をU、V、輝度値に対してPU、V、周辺画素を1群とするデータ以下テンプレートというサイズを、N×Nサイズと代数で示すと、テンプレート単位での画質補正では、下記の[数1]で計算した値で注目画素PU、Vを置き換える操作を行う。図2に示すように、例えば、7×7テンプレートの構成では、中心のマスが着目画素であって、この着目画素の周りが周辺画素となる。例えば、図3に示すような入力画像を取り上げる。そして、この図3の入力画像へSOBEL等のエッジ検出フィルタ処理を加えた場合、路面テクスチャの影響で、図4に示すように、目的外の小エッジが出力される等で好ましくない影響が生じる。よって、事前に平均化フィルタによる画質改善等の作業を施しておくことが望まれる。3×3テンプレート及び7×7テンプレートにて画素を事前に平均値平均化フィルタにより改善した画像を、図5、図6に示す。図5では、3×3テンプレート平均化フィルタによる入力画像画質補正を示す。図6では、7×7テンプレート平均化フィルタによる入力画像画質補正を示す。但し、例に挙げた画像の場合、次のような課題が見受けられる。3×3テンプレートでは、画質改善効果が弱い。7×7テンプレートでは、画質改善効果は強くなるが、遠方に映る後続車等ではテクスチャ情報が失われるためピンボケ画像となり必要な情報形状特性を得られなくなり、また、計算量が増える。更に、必要でない画素部分背景や遠方に小さく映る対象等に対して7×7サイズで計算を行うと、計算量の増加と共に、上述のピンボケの両問題が発生するため、特に課題となる。このように通常のテンプレート単位のデジタルフィルタ処理では、画面に対し一様に2次元画像サイズで設定したテンプレートを計算単位に使用するため、屋外で広範囲を捉えた撮像手段を使用する場合等、監視対象と撮像手段の相対距離車間距離が多々異なる対象を扱うため、離れた位置にある対象の必要な情報を取り逃すおそれがある。よって、最適なテンプレートサイズの決定に間しては、3次元での相対距離や位置関係及び実際の3次元スケールに対応した最適なテンプレートサイズを設定にするための指標が必要であると考えられる。しかしながら、特に、後方駐車支援の撮像手段等で頻繁に使用される広角レンズを持つ仕様では、入力画像にはレンズ収差による画像歪み等が関与する等の課題が存在するため、距離に関する基準の取得が通常困難となる。よって、任意の画素位置において実空間と相関のとれた最適なテンプレート選択には、さらに検討の余地があった。以下に、この解決手段について述べる。上述したように、画像をデジタルフィルタ処理により最適に補正を行うためには、対象との実空間3次元上での相対位置距離を計測した上で、設定した実空間スケールと相間のとれた最適なサイズのテンプレートを選択し切り替える手法が望ましい。しかしながら、従来の手法では、広角レンズのレンズ収差歪を含む画像では、画像上の画素位置と実空間の対応が困難である。そのため、先ず、実空間3次元上での相対位置距離の計測作業で行き詰る。従来、広角レンズを使用した後方確認システムの障害認識方法では、計測マップ手法により、歪みを持つ広角レンズを含めて表示される画面上での画素位置と実空間座標の関係を円滑に相関づけることができる手法を得ている。よって、撮像手段と対象の位置距離の関係についての課題は、解決できた。次に、必要とする3次元の実空間スケールと相関の取れたテンプレートサイズ定義手段を取り上げる。後方車両接近監視システムにあっては、近距離にある車両全体のスケールを画像上で1つのテンプレートで取り扱うよう対応するのは非常に計算コストが高く、さらに、歪みの影響等により認識等に及ぼす変形の度合いが増加することが伺え、非効率的であると考えた。よって、今回対象とするテンプレートに相当する実空間スケールを、車両のある1部位程度のそれ程大きくない範囲へ着目して作業を進める。例えば、路面等の荒れを考慮する場合、同一部位を指す場合の実空間上のスケールを5~30cm程度と設定しても路面等の荒れたテクスチャヘの対応が取れるテンプレートを入手できると予測される。続いて、計測マップと同じ分解能を持つデータテーブルを用意し、撮像手段の校正情報や計測マップと照らし合わせながら、実空間上で任意に設定したスケールの物体部位が計測マップ上の任意な実空間位置[X、Yconst、Z]に存在するものと仮想定義し、物体が画面上に表示されるサイズを計算処理により求めた上で、その表示サイズを最適なテンプレートサイズとしてデータテーブル上へ記録する。実際にフィルタリング処理で使用する際には、注目画素位置[U、V]に対応するデータテーブル登録値を参照しならがら読み取った登録テンプレートサイズで処理を、円滑に切り替えできるよう対処する。次いで、実空間位置に相関するテンプレートサイズで設定する方法を説明する。1、カメラキャリブレーションの実施先ず、システム構築時に使用する撮像手段をキャリブレーションして校正情報を取得する。使用する撮像手段が広角レンズを有する場合、レンズ歪みによる画像の変形を考慮した歪みパラメータも計算する。レンズ歪みを考慮したキャリブレーションは、広角レンズを使用した後方確認システムの障害認識方法と略同様に実施されるため、ここでは、その説明を省略する。これは、例えば、下記の[数2]~[数7]で示される。上記の[数2]~[数7]の式において、各パラメータは、以下の通りである。[Ud、Vd]:入力画像上の画素位置歪あり[Xcd、Ycd]:イメージプレーン上のカメラ座標位置歪あり[Xcm、Ycm]:イメージプレーン上のカメラ座標位置歪なし[Xwm、Yconst、Zwm]:実空間上のワールド座標位置Yconstについては後述[rnn、tn]:実空間上⇔カメラ座標系変換パラメータ位置姿勢変換情報[nn]:ひずみパラメータ[pn]:対数ひずみパラメータ[Rcm]:イメージプレーン上、原点からの距離[Rck]:イメージプレーン上、原点からの距離。近似モデル式を対数近似と通常近似で切り替える。2、切断面高さYconstでの計測マップの作成撮像手段からの入力画像は2次元U、Vであるため、そのままでは実空間の3次元パラメータX、Y、Zを全て特定することはできない。よって、拘束条件として、高さ方向で実空間での規定高さY位置を定める。この位置を切断面高さYconstと称呼する。この切断面高さYconstと実空間のXZ位置の入力画面上での関係を、図7に示す。同手法にて、Yconst=0の位置で、画像上の座標[Ud、Vd]に対応する実空間座標[Xwm、Ywm=0、Zwm]の座標系に対するXZ位置データの分布を計測マップとして定義する。計測マップと入力画像の対応については、図8A、Bに例を示す。図8Aは、X計測マップを画像へ重ねたものである。図8Bは、Z計測マップを画像へ重ねたものである。また、図9A、Bには、基準高さ位置Yでスライスしたイメージを示す。また、図9Bでは、対象フローと対象外フローとを示す。3、仮想球体の定義次に、登録するテンプレートサイズの計算作業に取り掛かる。実空間上の[Xwm、Yconst、Zwm]座標と画面上の[Ud、Vd]との位置相関が取れるのは、設定したYconstに対する計測マップが有効である画像範囲である。よって、各計測マップ上の[Xwm、Yconst、Zwm]位置に仮想的に設定した球体図10の円形で示すが存在している状況を想定し、想定した仮想球体が画像上で表示される場合の仮想球体が画面上に表示されるサイズを計算する。球体は撮像手段に対して正投影されるため、図10の楕円の円周とその内部で示される円形として映る。よって、画面上の表示では、図11に示すように、略円盤型に表示される。この実施例では、上述したように、テンプレートサイズに対する実スケールとの相関を車両部位単位で着目した5~30cm程度と設定したので、半径5cmの球体を仮定し取り扱う。また、ここで、仮想的に設定する物体を球体と定義したのは、ローカライズされた際、形状と表示に関わる偏重性を最も低く抑えられるためである。また、球体の場合、表示サイズの計算を後述の仮想円盤による代理計算で推定する際、回転に関する姿勢や形状維持に都合が良い。4、仮想円盤の設定方法の実施例仮想球体の画面上での投影サイズを求める。実際に仮想球体を両面に描画し、その画素数から表示サイズを求めても構わないが、ここでは、上述したように、仮想球体が画面上で撮像手段に正対する仮想円盤として表示されることを利用し、直接表示サイズを計算するのではなく、仮想円盤の表示面積の計算を利用して表示サイズを求める。以下に、計算手順について説明する。手順1先ず、カメラ座標系原点レンズ焦点位置をワールド座標系へ変換する。カメラ座標系原点レンズ焦点位置は、[Xcm、Ycm、Zcm]=[Xcd、Ycd、Zcd]=[0、0、0]であるので、カメラ座標系原点のワールド系座標を[Xwc、Ywc、Zwc]とすると、その位置は、上式の[数2]~[数7]より、下記の[数8]、[数9]のように求めることができる。レンズ焦点位置は、イメージプレーン上の中心となるのでレンズ歪の有無に関わらず、上記の座標変換式が適用できる。手順2次に、その計算対象となる、任意のある仮想球体の中心位置を、実空間上で[Xp、Yp、Zp]と定める。[Xp、Yp、Zp]は、使用する計測マップの分解能に依存する実空間座標であり、計測マップ上で示される[X、Z]位置の有限個の組み合わせの内の一つを指す。[Xwc、Ywc、Zwc]から[Xp、Yp、Zp]へのベクトル[Vecwx、Vecwy、Vecwz]を求める。図12には、仮想円盤球体の定義モデル図を示す。手順3上記の手順2で求めた[Vecwx、Vecwy、Vecwz]を上記の[数2]~[数7]によりワールド座標系からカメラ座標系へ変換し、さらに、単位ベクトル[eVeccx、eVeccy、eVeccz]へ変換してベクトル方向を計算する。続いて、カメラ座標系にて、カメラ座標系原点レンズ焦点周りに、単位ベクトル[eVeccx、eVeccy、eVeccz]について座標系3次元軸まわりの各回転角度[θc、φc、ψc]を求める。本処理では、後述する手順4で定義する円盤をカメラ座標系で回転させる場合の前処理を行っている。後述する手順4で示されるような単位ベクトル[1、0、0]に垂直な円盤については、Xc軸周りの回転と面積変化とは完全に独立するどれだけθc回転しても、表示面積は変化しない。よって、θc=0.0固定と設定することができる。θc=0とした場合、回転変換より[eVecx、eVecy、eVecz]と[θc、φc、ψc]との関係は、下記の[数10]の形で、容易に求まる。よって、下記の[数11]のように、[φc、ψc]を求めることができる。この[数10]、[数11]において、角度[θc、φc、ψc]がカメラ座標系で原点、すなわち撮像手段から見た仮想球体を示す姿勢角度となる。手順4次に、カメラ座標原点位置に立ち戻り、半径rで初期姿勢角度を、角度[θc0、φc0、ψc0]=[0、0、0]とする仮想円盤を設定する。この実施例では、この半径rに対して、例えば、5cmを適応する。仮想円盤の円周上の点を360度/N分割して、各円周上の分割点[Xc1、Yc1、Zc1]~[Xcn、Ycn、Zcn]位置を、[数12]のように求める。この場合、分割点数Nは任意とする。Nの数が大きい程、精度は向上するが、計算量が増加する。但し、当該処理を他のマップやデータテーブルの作成等、開発段階のオフライン処理に適用するのであれば、計算量の増加は、大きな障害とはならない。手順5各円周上のN点の座標を、カメラ座標原点周りに、角度[θc、φc、ψc]回転変換する。[Xci、Yci、Zci]→回転変換→[Xciθ、Yciφ、Zciψ]手順6[Xθciθ、Yciφ、Zciψ]で求まった各頂点を、仮想円盤の各座標に対して、上式の[数2]~[数7]の計算によりカメラ座標系[Xciθ、Yciφ、Zciψ]からワールド座標系[Xwiθ、Ywiφ、Zwiψ]に変換する。さらに、上記の手順2で得た[Vecwx、Vecwy、Vecwz]で実空間上をシフト移動させ、実空間ワールド座標系での仮想円盤の頂点[Xwmiθ、Ywmiφ、Zwmiψ]へ変換する。[Xwmiθ、Ywmiφ、Zwmiψ]は、下記の[数13]のように求められる。図13には、上記の手順4~手順6における仮想円盤球体の変換操作のイメージを示す。手順7[Xwmiθ、Ywmiφ、Zwmiψ]の全N点を、上式の[数2]~[数7]の透視座標変換により3次元空間座標から広角レンズの2次元画像平面上[Ud1、Vd1]~[UdN、VdN]に投影変換する。手順8画面上の投影点[Ud1、Vd1]~[UdN、VdN]に囲まれた領域について、表示面積Aを求める。上述したように、表示面積Aは実際に画像に[Ud1、Vd1]~[UdN、VdN]で囲まれた領域を描画し、標準画素数から求めることもできるが、円盤内部に相当する位置を計算し、円盤全体を描画する作業、描画後の描画点数をカウントする作業を計測マップの有効範囲内で全ての画素に対して適用する処理は、負荷が大きく実用性が低い。そこで、以下のように仮想円盤をN角形と見立てて、ヘロンの公式等を利用しN角形を構成するN-1個の3角形の面積を求め総和をとることで、近似的に表示面積Aを取得する方法を採用する。先ず、[Ud1、Vd1]~[UdN、VdN]と同様に、仮想円盤中心の[Xp、Yp、Zp]に関して投影座標変換を行い、画像平面上の位置[Ud0、Vd0]を求める。[Ud0、Vd0]と[Udi、Vdi]、[Udi1、Vdi1]の隣り合う2頂点で囲まれる3角形に着目する。[Ud0、Vd0]、[Udi、Vdi]、[Udi1、Vdi1]の3点を結ぶそれぞれの辺[ai、bi、ci]を、下記の[数14]で求める。そして、[U1、V1]~[UN、VN]まで繰り返し実行してヘロンの面積公式で求めた各分割三角形の面積総和により仮想円盤の画面上での表示面積Aを、下記の[数15]で求める。この場合、[Ud0、Vd0]あるいは[Ud1、Vd1]~[UdN、VdN]の一部が表示画面外に変換されるときは、該当辺を画面の枠との交点でクリップし、画面の隅とクリップ点とで囲まれる領域を計算する。そして、全ての点が画面外に変換されるとき、三角形が一つも形成されないような位置は、表示面積Aを零0でクリップし無効扱いとする。求めた表示面積Aは、[pixe12]の次元を持つ値である。表示面積Aは、距離[Xp、Yconst、Zp]に対する規定スケールの仮想ブロック今回は5cmの仮想球体を良好に相関づけた有効な指標である。この登録サイズAを計測マップと同じ分解能を持つデータテーブルを準備し記録する。例えば、8bit計測マップならX、Zで、256×256サイズのデータテーブルを用意する。また、仮想円盤の表示サイズは、計測マップの有効範囲図14A、図14Bで示すでしか計算できない。よって、8bit計測マップの場合は、Xが±127、Zが0~256分解能のデータテーブルサイズとなる。図14には、X計測マップ、Z計測マップを示し、無効領域計測マップでは実空間位置を測定できない領域が存在している。図15には、8bit計測マップに対応したデータテーブルを示し、登録した仮想円盤表示サイズを輝度で表現し、色が明るい画素位置程、登録半径サイズが大きく車間距離=近、暗部輝度が0値の位置は、遠すぎて半径1以下になった領域、あるいは、近すぎて計測不可となった領域を示す。上記のデータテーブル上の表示サイズAについては、距離情報X、Zの位置は計測マップの分解能に沿ったものであり、画面表示の際には、図14で示す有効エリアにしか使用できない。そのままフィルタ処理を適用してしまうと、無効領域と有効領域とで全く画質の異なる画像に改質されてしまう。図16には、画像サイズに対応した拡張データテーブルを示し、登録した仮想円盤表示サイズを輝度で表現したものである。図17には、拡張データテーブルによる平均化フィルタを示し、図15で無効輝度が0値となった仮想円盤表示サイズ未登録位置で、フィルタリング効果が得られず、処理可能不能領域で処理ムラが発生していることが分かる。データテーブルを計測マップ分解能に沿ったサイズから画像サイズヘ拡張する。拡張の際は、画像の[U、V]座標に相当する計測マップ[X、Z]を取得し、[X、Z]により参照したデータテーブル上の仮想球体表示サイズ値を拡張データテーブル上ヘコピーする。無効領域については、有効領域との境界に位置する登録データを用いて穴埋め補完する。この一貫した手順を踏まえることで、画面と同じサイズの拡張データテーブルを構築できる。図18には、補完済み拡張データテーブルを示し、図16での拡張データテーブルの無効箇所を補完している。図19には、補完拡張データテーブルによる平均化フィルタを示し、全域でフィルタリング効果が発揮され、処理ムラの発生を回避可能となる。次いで、実際のフィルタリングに関する表示面積Aを登録したマップの運用方法を、図20のフローチャートで示す。この図20では、拡張データテーブルの作成までと、実際のフィルタリング処理を示した例を示す。図20に示すように、先ず、カメラキャリブレーションを実施し、カメラパラメータを取得しS01、X、Z半径マップを作成し図8A、図8B参照、マップ対応画素位置の仮想球体表示サイズを計算しS02図11参照、データテーブルを作成し、仮想球体表示サイズ登録計測マップ分解能サイズ例えば8bit256×256としS03図15参照、そして、拡張データテーブルを作成し、仮想球体表示サイズ登録画面サイズ例えば、VGA640×480としS04図18参照、デジタルフィルタリング作業をし、図21に示すように、中心の着目画素を基点に、矢印の順番に沿って周辺画素データで補完する。そして、使用した補完データの総数が拡張データテーブルに登録した球の表示面積Aと等価になった時点A=i×j個で、補完データの追加を終了し、着目画素のデータを確定するS05。5、半径マップの作成仮想円盤の大きさを示す場合、上記の表示面積Aではなく、表示面積Aと同等サイズの近似円盤として取り扱い、表示面積Aの代わりに近似円盤による近似半径を拡張データテーブル上に登録することで、さらに、次のようなメリットを得ることができる。登録した近似半径を利用したテンプレートサイズ決定を行うことで、実スケールとの対応のとれたテンプレートを取得できる。デジタルフィルタ型の各種画像処理をテンプレート単位の操作で実施できるため、計算効率が向上する。また、あらかじめ検出した半径rの全ての組み合わせに対してテンプレートを用意し、半径rに応じた参照のためのオフセットアドレステーブル図22参照を用意しておくことで、さらに、作業効率を向上できる。図22には、オフセットアドレステーブルの例画像サイズ幅640高さ480を示し、注目画素のX、Y座標に対するオフセットを加えることで、計算手順を効率化できる。例えば、図23に示すように、画素[j、i]位置に登録されている半径データが2ならば、5×5テンプレート上下左右2マスを適用したフィルタリングを行う。周囲の画素参照の際に、各画素のアドレス計算には、図22のUVオフセット値を参照すれば、効率化できる。拡張データテーブルヘの登録サイズをAの平方根レベルで削減できるため、メモリ削減低コスト化に有効である。近似半径Rad[Xp、Yconst、Zp]の計算方法を、下記の[数16]に示す。表示面積Aの代わりに、半径rを登録した拡張データテーブルを、以降半径マップと呼称する。デジタルフィルタリングで処理を行う場合は、半径マップより参照した半径を上下左右に持つテンプレートを利用する。半径マップと登録半径によるテンプレートサイズの調整で処理を実行する方が、上記の拡張周辺データと仮想球体表示サイズに対するスパイラルな移動を行う計算よりも運用が容易であるため、利便性を向上できる。高さ方向を規定した位置で対象物の位置情報を推定するシステムについては、対象物の測定高さが規定高さに近い場合は、半径マップに登録された値を基に表示サイズを推定し、フィルタリング処理に使用するカーネルのサイズを切り替えることで、車間距離に対応したより効果的なフィルタ効果を得ることができる。ところで、規定高さから外れた対象物の部位や背景等の物体に誤った感度のフィルタが設定されることとなり、作業への課題となるおそれも考えられる。例えば、図24に示すA点においては、基準高さ付近の部位は、仮想球体より求めた半径サイズに対応した適切なフィルタリング効果を得ることができるので、表面荒れ等の影響を緩和したマッチング探索処理を行うことができる。しかし、図24に示すB点においては、基準高さから外れた部位は、適切なフィルタリング効果を得ることはできない例えば、特微量が失われたぼやけた形状になる。よって、探索は基本高さ付近の部位を注視した形で実施されるので、図24のB点の位置はノイズとして処理される候補と考えられる。元々、何度か取り上げてきた計測マップのように規定高さに位置する対象に着目して検知抽出するシステムでは、規定高さから外れた部位の情報はそもそも位置の特定が困難な対象として処理対象外として取り扱う。このため、フィルタ操作した画像上でA→B点までの輝度変化が急激なものではなく一定の連続性が保てていれば、マッチング処理に間する影響は小さく、広角レンズを使用した後方確認システムの障害認識方法による車両検出手段等へ十分応用できる。変形例1ガウシアンフィルタ荷重平均フィルタへの展開第一定義した仮想球体スケールの単位で表面に存在するテクスチャをフィルタリングにより除去し平滑化する場合、除去対象となる表面のテクスチャのバラツキを画像上の輝度の分散として捉えれば、平均化フィルタだけではなく、ガウシアンフィルタ等による荷重平均フィルタへの応用も可能である。例えば、半径マップの半径登録値を2倍した直径を、±分散範囲として扱い、ガウシアンフィルタを作成する。半径rのテンプレートに対応した荷重平均フィルタに関する補正操作と荷重係数に関わる式を、下記の[数17]~[数19]に示す。上記の式において、σはガウス分布による分散値で対象となる仮想球面スケールの表面粗さを示し、実験的に取得した固定値を使用する。また、Kは、wj、i総和によるPU、Vへの補正項浮動小数表記の際は不要K=1.0である。更に、荷重平均係数wj、iについては、毎回計算を行わずに対応できるよう12r2個のテーブルを事前に作成し、テンプレート計算時の位置j、iに応じて参照する方式が効率的である。変形例2ガウシアンフィルタ加重平均フィルタへの展開第二上記の変形例1では、分散値σを固定し、テンプレートのサイズを半径マップからの取得値により変化させ、荷重平均処理を行う手法を提言した。この変形例2では、操作するためのカーネルサイズを固定し、分散値σを半径マップからの取得値により変更する手段を取り上げる。つまり、分散値の範囲が実空間スケールで表した仮想球体表面に存在するテクスチャであると仮定して操作を行うと言い換えることができる。車間距離が短くなり表示が大きくなるもの程、表面のテクスチャに対する解像度が向上するためばらつきが大きくなると解釈する。この場合のガウシアン荷重平均フィルタに間する補正処理と荷重係数に関わる式を、下記の[数20]~[数22]に示す。ここで、Kは、wj、i総和によるPU、Vへの補正項浮動小数表記の際は不要K=1.0である。また、τrU、Vは、分散値を示す。また、分散値rは半径マップから参照した半径rU,Vに依存する。荷重平均係数wj、iについては、毎回計算を行わずに対応できるようi×j×r個のテーブルを事前に作成し、半径マップrの値に応じて係数参照する方式が効率的である。応用例1SUZANオペレータヘの円形オペレータの応用画像上の特徴点抽出手法の中にSUZANオペレータと呼ばれる円形カーネルを使用した手法が一般的に知られている。よって、半径マップからテンプレートを作成する作業にならって、円形のマスクを作成し、SUZANオペレータに適用することで、車間距離や実スケールに沿った特徴点検出処理を行うことができる。応用例2HARRISオペレータLUCASKANADEオプティカルフロー画像上の特徴点抽出手法の中に、HARRISオペレータと呼ばれる手法が一般的に知られている。また、画像上の特徴点の移動成分の計測手法中に関して、LUCASKANADEオプティカルフローと呼ばれる手法が一般的に知られている。HARRISオペレータやLUCASKANADEオプティカルフロー処理操作には、精度向上のため重み係数として、しばしばガウシアンフィルタが採用される。このガウシアンフィルタ係数についても、上記の変形例1、2で定義した手法により距離に応じて切り替えたガウシアンフィルタを採用し、距離やスケールに応じたフィルタの切り分けを行うことで、検出性能を向上することができる。また、言うまでもなく、取り扱うデータ群周辺画素の数を半径マップ登録値で切り替えることで、精度を向上することができる。応用例3その他フィルタヘの展開LOGフィルタ他、DOG法、SHIFT法SURF法、CANNY法LOGLAPRACIANofGAUSSIAN法については、荷重平均と2次微分を同時に実施するデジタルフィルタである。LOGに距離と実スケールを考慮したフィルタの選択適応も有効であると判断される。その補正係数式のみ抜粋する。上記の変形例2に相当するLOGフィルタは、下記の[数23]で表わされる。ここで、Kは、wj、i総和によるPU、Vへの補正項浮動小数表記の際は不要K=1.0である。変形例3に相当するLOGフィルタは、下記の[数24]で表わされる。ここで、Kは、wj、i総和によるPU、Vへの補正項浮動小数表記の際は不要K=1.0である。LOGフィルタは、ゼロクロス判定等の極値判定法を併用しエッジ検出や特徴点検出を行うための画像改質手法であるが、計算負荷が高く、しばしば異なる分散値を使用したガウシアンフィルタによる荷重平均画像の差分値により代理演算を行うDOGDOG=DifferenceOfGaussian法が代用される。当然DOGでは、ガウシアンフィルタによる荷重平均演算を使用するため変形例1あるいは変形例2に沿って距離やスケールに応じた補正切り替え手法が有効に機能すると考えられる。つまり、変形例1ではσの値を、変形例2では計算基準となるテンプレートのサイズを調整しDOGで使用する画像を作成する。また、荷重平均操作を元画像と元画像を等倍縮小した画像に適用し、それぞれサイズの違う画像同士の倍率について考慮しながらDOG処理を行うピラミッドDOG操作に基づいたSIFT法、SURF法等の特徴点検出手法に関する用例がコンピュータ性能の向上に伴い近年増加しているが、これらの処理についても有効性が予測される。また、CANNYエッジアルゴリズムの中にもゼロクロス判定は応用されており、その他アルゴリズムを考慮しても、半径マップの適用範囲は広いと考えられる。変形例3直径マップと円形カーネルマスク設定による精度向上方法半径マップに登録された半径よりテンプレートを作成する。この場合、半径rを上下左右方向に持つテンプレートは、奇数倍で増加するサイズで表現されることとなり、感度調整について若干その精度が劣化する。拡張データテーブル内に登録するデータサイズは倍増するが、登録データを半径から直径データヘと置き換え、さらに、次の手順でテンプレートとテンプレート内部の計算カーネルを切り替えることで、よりきめ細かい作業を行うことができる。直径データからテンプレートを作成する場合に、奇数の直径のケースと偶数の直径のケースで再現するためのテンプレート情報を独立して確保する。半径rに基づくテンプレートは、通常奇数サイズとなる。よって、直径マップに登録された値が奇数の場合は、通常の方形テンプレートをデータロードして作業を行う。直径マップに登録された値が偶数の場合は、1つ上の奇数テンプレートサイズでテンプレートを確保し、内部の計算カーネルの形状を円形となる形式で登録しておき、利用する周辺のデータについて参照総数を調整する。円形の内部カーネルについては、上記の応用例1のSUZANオペレータを参考できる。直径マップ登録値が2の場合は、図25Aに示すように、3×3テンプレートで円形の内部計算カーネルを使用しフィルタリングを行う。直径マップ登録値が3の場合は、図25Bに示すように、3×3テンプレートで方形の内部計算カーネルを使用しフィルタリングを行う。直径マップ登録値が4の場合は、図26Aに示すように、5×5テンプレートで円形の内部計算カーネルを使用しフィルタリングを行う。直径マップ登録値が5の場合は、図26Bに示すように、5×5テンプレートで方形の内部計算カーネルを使用しフィルタリングを行う。一方、直径マップ登録値が1以下の場合は、フィルタリング操作をせず、元のデータをそのまま使用する。即ち、この実施例においては、歪みを持つ広角レンズについて、カメラキャリブレーションによりカメラパラメータを取得し、計測マップのX,Z位置を作成し、各画素位置に対応した3次元位置情報を取得し、画面に映り込む対象の画素位置からテンプレートを切り替えるための距離の指標を取得し、さらに、入力画像と同じサイズのデータテーブルを用意する=計測マップサイズも同様。そして、各画素位置に対し、仮想的に規定の3次元寸法に基づく仮想ブロックの設定方法を定義する。例えば、半径30cmの球体等である。そして、計測マップで定義した各画素位置XYZ座標について、この仮想ブロックを仮想設定し、カメラパラメータにより画面上へ描画する。その表示サイズを元にテンプレートのサイズを決定する。この例の場合、30cmの球体が、画面上に映し出される場所により1pixe1になったり、10pixe1になったりすることで、テンプレートのサイズを切り替えられる。そのサイズは、モデルの実寸サイズに応じて設定されたものであるので、理想的な最適テンプレートサイズ切り替え手段と高い相関を持つ選択手段をえることができる。また、上記のキャリブレーションを広角レンズに基づく手法で行うことで、レンズ収差を持つ撮像手段による画像へも適用できる。上記の用意したデータテーブル上の対応する画素位置に、上記のサイズを定義する。計測マップ有効範囲について、全ての画素についてデータテーブルヘ指標を記録する。この計測マップの無効範囲については、有効範囲との境界に沿ったデータで補完しテーブルを完成する。この場合、テーブルヘの登録指標を仮想ブロックの表示サイズ値から、サイズ領域形状に基づく近似円盤に見立てた近似半径として登録すること半径マップで、テーブルの実メモリサイズの削減ができ、また、実際の画像フィルタリング処理に際して、より円滑に行うことができる。球体を仮想モデルに利用する効果について補足する。半径マップに対して、例えば、仮想球体の半径を単位円1m等と設定しておくことにより、実際に、50cmの仮想球体サイズを基準とした処理を行う際に、ゲイン1/2を加えることで、代替処理が可能となる。この場合、偏重性の少ない球体をモデルに使用したことで、他の形状モデルを使うよりも発生する誤差を少なくする効果が期待される。なお、この発明の変形例としては、半径マップを応用した1次方向直線方向偏重性の探索によるノイズ耐性の強い高速白線検出アルゴリズム、半径マップを応用したオプティカルフロー探索が可能である。"