(57)【要約】【課題】本発明は、正解航跡候補を含めた航跡候補をその状況に適切に生成保持することができ、探知状況によらず高い追尾精度を維持することができる目標追尾装置を提供する。【解決手段】多次元相関行列生成部１２１は、センサ２００の観測値のうち最新の観測値に対応する最新の観測フレームから新しい順に所定の指定フレーム数の観測フレームを選択する。多次元相関行列生成部１２１は、その選択した観測フレームに対応する航跡候補と、その航跡候補についてのコストと含む多次元相関行列を生成する。多次元相関行列生成部１２１は、多次元相関行列の生成に際して、全航跡候補と相関する観測値が存在しない観測フレームである探知抜けフレームを多次元相関行列の観測フレームから除外し、探知抜けフレームよりも前の観測フレームを選択する。【選択図】図１
2010286275.tif 000002
【特許請求の範囲】
【請求項１】
センサの観測値を用いて目標物の航跡を推定する目標追尾装置であって、
前記センサの観測値を取得し、その取得した観測値を既存の航跡候補に割り当てて、既存の航跡候補を更新し、その更新した航跡候補についての尤度に関する指標であるコストを算出する航跡候補生成部と、
前記センサの観測値のうち最新の観測値に対応する最新の観測フレームから新しい順に所定の指定フレーム数の観測フレームを選択し、その選択した観測フレームに対応する航跡候補と、その航跡候補についてのコストと含む多次元相関行列を生成する多次元相関行列生成部と、
前記多次元相関行列生成部によって生成された前記多次元相関行列を用いて、多次元の相関決定を行い、当該観測フレームで航跡候補の組み合わせを抽出する多次元相関決定部と、
抽出された航跡候補の組み合わせを次の観測フレームでの既存の航跡候補として設定するとともに、抽出された航跡候補について航跡確立する航跡決定部と、
次の観測フレームでの既存の航跡候補について、航跡表示用の平滑値と、次の観測フレームで取得する観測値を当該既存の航跡候補に割り当てるための予測値とを算出する航跡候補追尾処理部と
を備え、
前記多次元相関行列生成部は、前記多次元相関行列を生成する際に、
選択した観測フレームが、航跡候補に割り当てられる観測値が存在しない探知抜けフレームであれば、その探知抜けフレームを除外して、
航跡候補に割り当てられる観測値が存在する前記指定フレーム数の探知フレームを選択するために、前記探知抜けフレームよりも前の観測フレームを選択する
ことを特徴とする目標追尾装置。
【請求項２】
前記多次元相関行列生成部は、前記多次元相関行列を生成する際に、最新の観測フレームからの前記探知抜けフレームを含めた観測フレームの個数が、所定の上限フレーム数に達した場合には、前記探知抜けフレームをも選択し、選択した探知フレーム及び探知抜けフレームを用いて前記多次元相関行列を生成する
ことを特徴とする請求項１記載の目標追尾装置。
【請求項３】
前記航跡候補生成部による航跡候補生成結果に基づいて、複数の航跡が互いに近接するか否かを判断し、互いに近接する複数の航跡については同一のクラスタで処理し、互いに離れた航跡については異なる複数のクラスタで処理するように、クラスタの生成・統合・分離を行って、クラスタに航跡候補を振り分けるクラスタ処理部
をさらに備え、
前記多次元相関行列生成部は、
前記クラスタ処理部によるクラスタ生成・分離の際には、クラスタに振り分けられた航跡候補と、その航跡候補についてのコストとを含む多次元相関行列を生成し、
前記クラスタ処理部によるクラスタ統合の際には、統合される複数のクラスタのうち最も長い観測フレーム数を用いているクラスタの観測フレーム数に対応させて、前記多次元相関行列を生成する
ことを特徴とする請求項１又は請求項２に記載の目標追尾装置。
【請求項４】
前記航跡候補生成部による航跡候補生成結果に基づいて、複数の航跡が互いに近接するか否かを判断し、互いに近接する複数の航跡については同一のクラスタで処理し、互いに離れた航跡については異なる複数のクラスタで処理するように、クラスタの生成・統合・分離を行って、クラスタに航跡候補を振り分けるクラスタ処理部
をさらに備え、
前記多次元相関行列生成部は、
前記クラスタ処理部によるクラスタ生成・分離の際には、クラスタに振り分けられた航跡候補と、その航跡候補についてのコストとを含む多次元相関行列を生成し、
前記クラスタ処理部によるクラスタ統合の際には、統合後の航跡候補の観測フレーム数である統合時フレーム数が所定値となるように、統合される複数のクラスタのそれぞれで多次元相関行列を生成し、前記クラスタ処理部によるクラスタ統合とともに当該複数のクラスタのそれぞれの前記多次元相関行列を結合し、指定フレーム数の前記多次元相関行列を生成する
ことを特徴とする請求項１又は請求項２に記載の目標追尾装置。
【請求項５】
前記航跡候補生成部による航跡候補生成結果に基づいて、複数の航跡が互いに近接するか否かを判断し、互いに近接する複数の航跡については同一のクラスタで処理し、互いに離れた航跡については異なる複数のクラスタで処理するように、クラスタの生成・統合・分離を行って、クラスタに航跡候補を振り分けるクラスタ処理部と、
統合後の航跡候補の観測フレーム数である統合時フレーム数を、可変及び固定のいずれか一方に設定する統合時フレーム数ルール設定部と
をさらに備え、
前記クラスタ処理部によるクラスタ生成・分離の際には、クラスタに振り分けられた航跡候補と、その航跡候補についてのコストとを含む多次元相関行列を生成し、
前記クラスタ処理部によるクラスタ統合の際には、
前記統合時フレーム数ルール設定部により統合時フレーム数が可変に設定されている場合に、統合される複数のクラスタのうち最も長い観測フレーム数を用いているクラスタの観測フレーム数に対応させて、前記多次元相関行列を生成し、
前記統合時フレーム数ルール設定部により統合時フレーム数が固定に設定されている場合に、前記統合時フレーム数が所定値となるように、統合される複数のクラスタのそれぞれで多次元相関行列を生成し、前記クラスタ処理部によるクラスタ統合とともに当該複数のクラスタのそれぞれの前記多次元相関行列を結合し、前記所定値の前記多次元相関行列を生成する
ことを特徴とする請求項１又は請求項２に記載の目標追尾装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、例えばレーダ等のセンサシステムに用いられ、センサからの目標位置についての観測値に基づいて目標の航跡を生成する目標追尾装置に関する。
【背景技術】
【０００２】
目標追尾装置では、一般的に、センサからの目標位置の観測値を用いて目標の航跡を生成する際の相関を多次元的に解く技術が使用される。このような技術では、探知状況が比較的良好な場合には、所定の指定フレーム数で多次元の相関決定が行われる。他方、探知状況が比較的悪く探知抜けが多く発生している場合には、指定フレーム数内の観測フレーム数で航跡候補の区別がつかずに正解航跡が破棄されるような状況において、指定フレーム数を超えた数の観測フレームを用いて、多次元の相関決定を行う。即ち、各探知状況に応じて、航跡候補の推定に用いる観測フレーム数を自動的に変化させる。このように、各探知状況に応じた適切な観測フレーム数で相関決定を行うことによって、比較的精度のよい目標追尾を実現する。
【０００３】
この多次元的に相関を解く技術としては、複数サンプリング時刻に跨る観測値の相関をまとめて計算して航跡を生成するアルゴリズムである複数フレーム割当法がある（例えば、特許文献１及び非特許文献１参照）。このような複数フレーム割当法では、既存航跡とＮフレーム時刻に跨る観測値との相関決定をＮ次元の割当問題に変換して解く。このＮ次元割当問題は、以下の（１）のように表すことができる。
【０００４】
【数１】
2010286275.tif 000003
【０００５】
この式（１）のＣｉ2・・・ＣｉＮ＋１は、既存航跡及び観測値ｉ１・・・ｉＮ＋１から航跡を生成する際に用いられる航跡尤度に関するコストであり、通常航跡尤度の逆数の対数値である。変数ｚｉ１・・・ｚｉＮ＋１は、次の式（２）の通りである。また、（１）の割当問題には、次の式（３）のような制約条件がある。
【０００６】
【数２】
2010286275.tif 000004
【０００７】
ここで、ｎフレーム目における航跡候補は、前回（ｎ−１フレーム目）の航跡候補ｘｎ−１(＋)を基に算出された予測値ｘｎ(−)から、相関関係にある当該フレームの観測値を割り当てることで更新される。ここでは、一例として、線形カルマンフィルタを想定した場合の航跡候補の更新、及び残差共分散行列を用いた相関判定について、簡単に説明する。
【０００８】
航跡候補のｘｎ−１（＋）のｎフレームでの予測値は、ｎ−１からｎフレームにおける推移行列をΦ（ｎ−１）、ｎフレームでの観測行列をＨ（ｎ）とすると、次の式（４）の通りである。
【数３】
2010286275.tif 000005
【０００９】
また、ｎ−１フレームでの平滑誤差共分散行列Ｐｎ−１（＋）、駆動雑音行列Ｑ（ｎ−１）から予測誤差共分散行列Ｐｎ（−）は、次の式（５）の通りである。
【数４】
2010286275.tif 000006
【００１０】
残差共分散行列Ｓｎは、観測誤差共分散行列Ｂｎを用いて以下の式（６）のように求まる。
【数５】
2010286275.tif 000007
【００１１】
そして、複数フレーム割当法では、残差行分散行列Ｓｎを用いて、以下の式（７）が成立する場合に、航跡候補と観測値ｚ（ｎ）は相関があると判定する。ここで、ｄはゲートサイズであり観測値の持つ自由度Ｍのカイ２乗分布表から定められる事前に設定されたパラメータである。
【数６】
2010286275.tif 000008
【００１２】
ｎフレームにおける航跡候補ｉ（尤度γｉ）がｎ−１フレームにおける航跡候補ｌ（尤度γｌ）から生成されたとすると、探知確率ＰＤ、ゲート捕捉率ＰＧ及び航跡候補生成における観測値との相関度ｇｉ，ｌを用いた際に、探知抜けによる航跡候補の生成における尤度γｉ，０と、観測値による航跡候補の生成における尤度γｉ,ｌとは、以下の式（８）の通りである。そして、更新された各航跡候補の尤度は、式（８）の演算処理により更新される。
【数７】
2010286275.tif 000009
【００１３】
なお、既存航跡候補と観測値の割り当て結果とを用いた相関決定の処理を効率化させるための技術としては、互いに観測値を共有する航跡候補からなるクラスタを生成し、クラスタ毎に処理を実行する方式がある（例えば、特許文献２参照）。
【先行技術文献】
【特許文献】
【００１４】
【特許文献１】特開２００７−２１２２４４号公報
【特許文献２】特開平８−２７１６１７号公報
【非特許文献】
【００１５】
【非特許文献１】A.B.Poore and A.J.RobertsonIII "A New Lagrangian Relaxation Based Algorithm for a Class of Multidimensional Assignment Problems" Computational Optimization and Applications, 8, 129-150, 1997.
【発明の概要】
【発明が解決しようとする課題】
【００１６】
特許文献１及び非特許文献１のような複数フレーム割当法を用いた従来の目標追尾装置では、探知抜けフレームも多次元の相関決定に利用されていた。このため、指定フレーム数内で多次元の相関決定の観測フレームに探知抜けフレームが比較的多く含まれている場合に、その指定フレーム数内で航跡候補の区別がつかず指定フレーム数内で起点となる観測値が同一である航跡候補が複数存在しても、１フレーム前に選択された航跡候補のみを抽出し、他の航跡候補を削除せざるを得なかった。
【００１７】
この結果、一時的な尤度の低下により選択されなかった正解航跡候補（目標信号からなる観測値を多く有する航跡）が次フレームの観測値を割り当てる前に削除されてしまい、各フレームで割り当てられた観測値の組み合わせが目標らしいと判断されるまでの時間である航跡確立の時間が長くなったり、追尾維持（確立した航跡と目標からなる観測値とを各フレームで割り当てる）に失敗したりするといった現象が発生していた。
【００１８】
ここで、図９，１０に、従来装置における追尾確立が遅延する様子を示す。図９の表は、各フレームでの観測値数、及びその観測値が目標からか不要信号かを示す。また、図９の行列は、各航跡候補の各フレームでの観測値の割当を示し、指定フレーム数（この例では５フレーム）では航跡候補の区別がつかず正解航跡が削除される様子を示す。図１０の航跡グラフにおいて、図９の探知状況及び行列のような情報を用いた場合には、航跡確立が遅延していることがわかる。
【００１９】
また、従来装置において、正解航跡候補の破棄を回避するために単純に指定フレーム数を多くした場合には、探知状況がよく探知抜けが発生せずに少ない指定フレーム数で十分航跡候補が区別できるにもかかわらず、長い指定フレーム数となることで不要に航跡候補が多くなる。この結果、過去の観測フレームを多く使うことにより過去の情報の影響を強く受けることから、例えば旋回等の高機動運動への追従性が損なわれたり、演算負荷が増大したりするという問題が生じる。
【００２０】
具体的に、図１１の矢示Ａのように、目標物の旋回動作等の高機動運動が生じた場合に、１０フレーム分の観測値を利用するときには、矢示Ｂのように、過去の観測値の影響を受けて目標の運動に対する追従性が低下し、目標物が直進しているものとして追尾確立してしまうこととなり、追尾精度が低下してしまう。なお、５フレーム分の観測値を利用するときには、矢示Ｃのように、比較的良好な追従性を示すこととなる。
【００２１】
この発明は、上記のような課題を解決するためになされたものであり、正解航跡候補を含めた航跡候補を適切に生成保持することができ、探知状況によらず高い追尾精度を維持することができる目標追尾装置を得ることを目的とする。
【課題を解決するための手段】
【００２２】
この発明に係る目標追尾装置は、センサの観測値を用いて目標物の航跡を推定するものであって、前記センサの観測値を取得し、その取得した観測値を既存の航跡候補に割り当てて、既存の航跡候補を更新し、その更新した航跡候補についての尤度に関する指標であるコストを算出する航跡候補生成部と、前記センサの観測値のうち最新の観測値に対応する最新の観測フレームから新しい順に所定の指定フレーム数の観測フレームを選択し、その選択した観測フレームに対応する航跡候補と、その航跡候補についてのコストと含む多次元相関行列を生成する多次元相関行列生成部と、前記多次元相関行列生成部によって生成された前記多次元相関行列を用いて、多次元の相関決定を行い、当該観測フレームで航跡候補の組み合わせを抽出する多次元相関決定部と、抽出された航跡候補の組み合わせを次の観測フレームでの既存の航跡候補として設定するとともに、抽出された航跡候補について航跡確立する航跡決定部と、次の観測フレームでの既存の航跡候補について、航跡表示用の平滑値と、次の観測フレームで取得する観測値を当該既存の航跡候補に割り当てるための予測値とを算出する航跡候補追尾処理部とを備え、前記多次元相関行列生成部は、前記多次元相関行列を生成する際に、選択した観測フレームが、航跡候補に割り当てられる観測値が存在しない探知抜けフレームであれば、その探知抜けフレームを除外して、航跡候補に割り当てられる観測値が存在する前記指定フレーム数の探知フレームを選択するために、前記探知抜けフレームよりも前の観測フレームを選択するものである。
【発明の効果】
【００２３】
この発明に係る目標追尾装置によれば、多次元相関行列生成部が、多次元相関行列を生成する際に、選択した観測フレームが探知抜けフレームであれば、その探知抜けフレームを除外して、探知抜けフレームよりも前の観測フレームを選択するので、探知抜けフレームを多次元の相関決定の観測フレームに含めないことで探知状況に合わせて自動的に多次元の相関決定に用いる絶対的なフレーム数を変化させ、航跡候補の区別ができるようにして多くの航跡候補を残すことが可能となることにより、正解航跡候補を含めた航跡候補を適切に生成保持することができ、探知状況によらず高い追尾精度を維持することができる。
【図面の簡単な説明】
【００２４】
【図１】この発明の実施の形態１による目標追尾装置を示すブロック図である。
【図２】図１の目標追尾装置の動作を示すフローチャートである。
【図３】図１の目標追尾装置の動作を示すフローチャートである。
【図４】実施の形態１の目標追尾装置及び従来装置のそれぞれの航跡候補選択動作を対比して説明するための説明図である。
【図５】この発明の実施の形態２による目標追尾装置を示すブロック図である。
【図６】図５の目標追尾装置の動作を示すフローチャートである。
【図７】目標番号が付与された２つの航跡候補群で最新フレームの観測値の取り合いが発生した場合の処理を説明するための説明図である。
【図８】目標番号が付与された２つの航跡候補群で最新フレームの観測値の取り合いが発生した場合の処理を説明するための説明図である。
【図９】従来装置における追尾確立が遅延する様子を説明するための説明図である。
【図１０】従来装置における追尾確立が遅延する様子を説明するための説明図である。
【図１１】指定フレームのフレーム長を長くした場合の追尾確立の様子を説明するための説明図である。
【発明を実施するための形態】
【００２５】
以下、この発明を実施するための形態について、図面を参照して説明する。
実施の形態１．
図１は、この発明の実施の形態１による目標追尾装置を示すブロック図である。
図１において、目標追尾装置１００は、例えばレーダ等のセンサ２００及び例えばモニタ等の表示装置３００に接続されている。また、目標追尾装置１００は、航跡候補生成処理系（航跡候補生成処理部）１１０、多次元相関行列生成部１２１、多次元相関決定部１３１、航跡決定部１４１及び航跡候補追尾処理系（航跡候補追尾処理部）１５０を有している。
【００２６】
航跡候補生成処理系１１０は、観測値割当部１１１及びコスト算出部１１２を含んでいる。観測値割当部１１１は、航跡候補追尾処理系１５０から受けた予測値を用いて、センサ２００から受けた観測値を既存の航跡候補に割り当てて、既存の航跡候補を更新する。また、コスト算出部１１２は、観測値割当部１１１によって更新された各航跡候補について、航跡尤度の逆数の対数値であるコスト（航跡尤度に関するコスト）を算出する。このコストは、予測値に対して観測値が尤もらしいかどうか（信頼度が高いか否か）を判別するための指標である。
【００２７】
多次元相関行列生成部１２１は、センサ２００の観測値のうち最新の観測値に対応する最新の観測フレームから新しい順に所定の指定フレーム数の観測フレームを選択する。そして、多次元相関行列生成部１２１は、その選択した観測フレームに対応する航跡候補と、その航跡候補についてのコストと含む多次元相関行列を生成する。
【００２８】
ここで、多次元相関行列生成部１２１は、多次元相関行列の生成に際して、全航跡候補と相関する観測値が存在しない観測フレームである探知抜けフレームを多次元相関行列の観測フレームから除外し、探知抜けフレームよりも前（過去）の観測フレームを選択する。即ち、多次元相関行列生成部１２１は、探知抜けフレームを多次元相関行列の観測フレームから除外することによって、全航跡候補と相関する観測値が存在する探知フレームを指定フレーム数分選択する。
【００２９】
多次元相関決定部１３１は、多次元相関行列生成部１２１によって生成された多次元相関行列を用いて多次元の相関決定を行い、当該観測フレームでの尤もらしい航跡候補の組み合わせを抽出する。航跡決定部１４１は、多次元相関決定部１３１によって抽出された航跡候補の組み合わせと、指定フレーム数（探知抜けフレームを除く）の前までの割り当てられた観測値が同一である航跡候補以外を削除し、残った航跡候補を次の観測フレームでの既存の航跡候補（選択航跡候補）として保持する。また、航跡決定部１４１は、抽出された航跡候補に、目標番号が付与されておらず、かつその航跡候補に対応する尤度が一定閾値（事前に設定されたパラメータ）を超えた場合には、航跡を確立し、その航跡候補に目標番号を付与する。
【００３０】
航跡候補追尾処理系１５０は、平滑値算出部１５１及び予測値算出部（航跡選択モデル参照部）１５２を含んでいる。平滑値算出部１５１は、航跡決定部１４１によって保持されている航跡候補の（航跡表示用の）平滑値を算出する。また、平滑値算出部１５１は、航跡決定部１４１によって目標番号が付与され、多次元相関決定部１３１によって抽出された航跡候補についての平滑値を表示装置３００に送り、その航跡候補についての平滑値を表示装置３００に表示させる。予測値算出部１５２は、先の式（４）を用いた演算処理を行い、当該フレームの各既存航跡候補の予測値を算出する。
【００３１】
ここで、目標追尾装置１００は、演算処理部（ＣＰＵ）、記憶部（ＲＯＭ、ＲＡＭ等）及び信号入出力部を持ったハードウェア（図示せず）により構成することができる。目標追尾装置１００のハードウェアの記憶部には、航跡候補生成処理系１１０、多次元相関行列生成部１２１、多次元相関決定部１３１、航跡決定部１４１及び航跡候補追尾処理系１５０の機能を実現するためのプログラムが格納されている。
【００３２】
次に、動作について説明する。図２，３は、図１の目標追尾装置１００の動作を示すフローチャートである。ここでは、センサ２００が新たに生成した観測値に対応する観測フレームをｋフレームとし、このｋフレーム目の処理について説明する。図２，３において、目標追尾装置１００の観測値割当部１１１は、センサ２００からｋフレーム目の観測値を取得する（ステップＳ１０１）。また、予測値算出部１５２は、先の式（４）を用いた演算処理を行うことによって、ｋフレーム目の各既存航跡候補の予測値を算出し、その算出した予測値を観測値割当部１１１へ送る（ステップＳ１０２）。
【００３３】
そして、観測値割当部１１１は、ｋフレーム目の観測値及び予測値から、航跡と観測値との割当を行う（ステップＳ１０３）。この後、コスト算出部１１２は、その割り当てられた観測値を用いて航跡候補を生成し、既存の航跡候補と観測値の相関度から算出した航跡候補の尤度から、航跡候補のコストを算出する（ステップＳ１０４）。なお、予測値算出部１５２、観測値割当部１１１及びコスト算出部１１２は、ステップＳ１０２〜Ｓ１０４の処理を、複数の航跡候補のそれぞれについて実行する（航跡候補分ループする）。
【００３４】
多次元相関行列生成部１２１は、当該観測フレームが探知抜けフレームであるか否かを判定する（ステップＳ１０５）。このときに、多次元相関行列生成部１２１は、当該観測フレームが探知抜けフレームでないことを確認した場合には、図３に示す多次元相関行列の設定処理を行う（ステップＳ２００）。他方、多次元相関行列生成部１２１は、当該観測フレームが探知抜けフレームであることを確認した場合には、前回の処理で割り当てられた観測値を流用し（ステップＳ１０６）、ステップＳ１０７へ移行する。
【００３５】
ここで、図３のステップＳ２００では、多次元相関行列生成部１２１は、利用する最大フレーム数である探索フレーム数（探索フレーム数ＮＳ）、及び多次元相関行列の設定フレーム数（Ｎｉ）を、次の式（９），（１０）のように設定する（ステップＳ２０１）。
【００３６】
【数８】
2010286275.tif 000010
【００３７】
そして、多次元相関行列生成部１２１は、当該観測フレームが探知抜けフレームであるか否かを判定する（ステップＳ２０２）。このときに、多次元相関行列生成部１２１は、当該観測フレームが探知抜けフレームでないことを確認した場合には、全航跡候補の当該観測フレームにおける観測値を多次元相関行列に設定する（ステップＳ２０３）。そして、多次元相関行列生成部１２１は、設定フレーム数を次の式（１１）のように更新する（ステップＳ２０４）。
【００３８】
【数９】
2010286275.tif 000011
【００３９】
この後、多次元相関行列生成部１２１は、設定フレーム数が多次元の相関決定の指定フレーム数に達したか否かを判定する（ステップＳ２０５）。このときに、多次元相関行列生成部１２１は、設定フレーム数が多次元の相関決定の指定フレーム数に達したことを確認した場合には、探査ループから抜ける。
【００４０】
これに対して、多次元相関行列生成部１２１は、設定フレーム数が多次元の相関決定の指定フレーム数に達していないことを確認した場合（ステップＳ２０５のＮＯ方向）、又は当該フレームが探知抜けフレームであることを確認した場合（ステップＳ２０２のＹＥＳ方向）には、探索フレームを１つ過去のフレームに更新し（ステップＳＳ２０６）、探索ループにより、ステップＳ２０２からの処理を繰り返す。
【００４１】
そして、多次元相関行列生成部１２１は、設定フレーム数が指定フレーム数に達したことを確認した場合、又は設定フレーム数が探索フレーム数（上限フレーム数）に達したことを確認した場合には、設定フレーム数を多次元相関行列のフレーム数とし（ステップＳ２０７）、図２のステップＳ１０７へ移行する。
【００４２】
ステップＳ１０７では、多次元相関行列生成部１２１は、多次元相関行列に設定されている各航跡候補に対応するコストとして、ステップＳ１０４で算出したコストを設定する。この後、多次元相関決定部１３１は、多次元相関行列生成部１２１によって生成された多次元相関行列を用いて多次元の相関決定を行い、尤もらしい航跡候補の組み合わせを抽出する（ステップＳ１０８）。
【００４３】
そして、航跡決定部１４１は、図２，３の一連の処理で使用した全ての航跡候補のうち、多次元相関決定部１３１によって抽出された航跡候補の組み合わせに含まれる航跡候補と、指定フレーム数（探知抜けフレームを除く）分割り当てられた観測値が同一である航跡候補以外を削除し、残った航跡候補を次の観測フレーム（ｋ＋１フレーム目）での既存の航跡候補として保持する。また、航跡決定部１４１は、抽出された航跡候補に目標番号が付与されておらず、かつ尤度が一定閾値以上である場合には、その航跡候補を確立航跡とし、目標番号を付与する（ステップＳ１０９）。
【００４４】
平滑値算出部１５１は、次フレームの既存の航跡候補とした航跡候補について、現在のフレーム（ｋフレーム目）での平滑値を算出する（ステップＳ１１０）。また、平滑値算出部１５１は、平滑値の算出処理を、次フレームでの既存の航跡候補のそれぞれについて実行する（航跡候補分ループする）。そして、その算出した平滑値を表示装置３００へ送り、目標追尾装置１００の一連の処理が終了する。
【００４５】
ここで、平滑値算出部１５１による平滑値の算出処理において、追尾フィルタに線形カルマンフィルタを用いた場合に、ｎフレーム目における平滑値ｘｎ（＋）、平滑誤差共分散行列Ｐｎ（＋）、ゲインＫ（ｎ）は、次の式（１２）〜（１４）に示すように記述できる。
【００４６】
【数１０】
2010286275.tif 000012
なお、平滑値を求める方法としては、線形カルマンフィルタを用いる他に、例えば、指定フレーム数分の観測値からの線形最小自乗法等を用いることもできる。
【００４７】
次に、実施の形態１の目標追尾装置１００と従来装置とを対比して説明する。図４は、実施の形態１の目標追尾装置１００及び従来装置のそれぞれの航跡候補選択動作を対比して説明するための説明図である。図４（ａ）は、従来装置における航跡候補抽出動作を示し、図４（ｂ）は、実施の形態１の目標追尾装置１００の航跡候補抽出動作を示す。
【００４８】
従来装置では、図４（ａ）に示すように、判断対象の破線枠の５フレームの指定フレームのうち探知フレームが大多数を占めており、判断対象の航跡候補は２種類である。このため、既存の航跡候補として保持される航跡候補が、同じ起点を有する２つの航跡候補となり、正解航跡候補が保持されずに破棄されてしまう。
【００４９】
これに対して、実施の形態１の目標追尾装置１００では、図４（ｂ）に示すように、探知抜けフレームについては除外して、航跡候補群から航跡候補を選択し、同じ起点を有する３つ以上の航跡候補が既存の航跡候補として保持される。この結果、実施の形態１の目標追尾装置１００では、正解航跡候補が保持される。
【００５０】
従って、実施の形態１の目標追尾装置によれば、多次元相関行列生成部１２１が、多次元相関行列を生成する際に、選択した観測フレームが探知抜けフレームであれば、その探知抜けフレームを除外して、探知抜けフレームよりも前の観測フレームを選択する。この構成により、探知抜けフレームを多次元の相関決定の観測フレームに含めないことで探知状況に合わせて自動的に多次元の相関決定に用いる絶対的なフレーム数を変化させ、航跡候補の区別ができるようにして多くの航跡候補を残すことが可能となる。この結果、正解航跡候補を含めた航跡候補を適切に生成保持することができる。よって、新たに観測値が得られた場合に多次元の相関決定において正解航跡が抽出されるようになり、相関性能を向上させることができる。
【００５１】
これに加えて、センサ２００の探知状況が良好な場合においては、指定フレーム数の観測フレームを利用して、航跡候補を生成するので、過去の観測値の影響を受けにくいことにより、目標物の高機動運動に対する追従性の低下を抑えることができ、探知状況によらず高い追尾精度を維持することができる。
【００５２】
また、最新の観測フレームからの探知抜けフレームを含めた観測フレームの個数が、探索フレーム数（上限フレーム数）に達した場合には、探知抜けフレームをも選択し、選択した探知フレーム及び探知抜けフレームを用いて多次元相関行列を生成する。この構成により、探知抜けフレームが比較的長く連続するような状況下においても、過去の観測値の影響を受けにくくすることができる。
【００５３】
なお、実施の形態１の目標追尾装置１００の動作においては、図３の探索フレーム数分のループ処理を省略することもできる。
【００５４】
実施の形態２．
先の実施の形態１では、全航跡候補を１つの多次元相関行列として多次元の相関決定処理を実行していた。しかしながら、航跡候補が観測値を共有しないような場合には、多次元相関行列を個別に生成し、多次元の相関決定処理を実行することによって、多次元相関行列の要素数を減少させ、多次元の相関決定処理に要する演算負荷を抑えることが可能である。この発明の実施の形態２では、観測値を共有する航跡候補の集合であるクラスタを用いて、探知抜けフレームを除いたフレームでの多次元の相関決定を実行するための形態について説明する。
【００５５】
図５は、この発明の実施の形態２による目標追尾装置１００を示すブロック図である。図５において、実施の形態２の目標追尾装置１００は、実施の形態１の目標追尾装置１００の構成に加えて、クラスタ生成統合分離部（クラスタ処理部）１６１及び統合時フレーム数ルール設定部１７１をさらに有している。なお、目標追尾装置１００のハードウェアの記憶部には、クラスタ生成統合分離部１６１及び統合時フレーム数ルール設定部１７１の機能を実現するためのプログラムがさらに格納されている。
【００５６】
クラスタ生成統合分離部１６１は、複数のクラスタに存在する航跡が同一観測値と相関がとれた場合にクラスタを統合する。また、クラスタ生成統合分離部１６１は、どのクラスタの航跡候補とも相関がとれない観測値について新たにクラスタを生成する。さらに、クラスタ生成統合分離部１６１は、クラスタ内で観測値を取り合わない航跡候補が生成された場合には各航跡を異なるクラスタに分離する。実施の形態２の多次元相関行列生成部１２１は、クラスタ生成統合分離部１６１によるクラスタの操作（統合・分離・生成）がされた際に、多次元の相関決定のためのフレーム数を設定し、新たに多次元相関行列を生成する。
【００５７】
ここで、クラスタの統合時には、各クラスタでの探知抜けフレームが異なるため、多次元相関行列で利用する観測フレーム数が統合前と大きく異なる場合がある。この場合、指定フレーム数で観測値の割り当てでは区別がつかない航跡候補が発生する。そこで、多次元相関行列生成時に、指定フレーム内で区別がつかない航跡候補を削除するか、多次元の相関決定では一時的にそのフレーム数を長くして、多次元の相関決定後に指定フレーム内での航跡決定を実行するかのどちらかの処理を実行する。このどちらの処理を実行するかについては、統合時フレーム数ルール設定部１７１が設定し、その設定に応じて、多次元相関行列生成部１２１がいずれかの処理を行い、多次元相関行列を生成する。
【００５８】
次に、動作について説明する。ここで、実施の形態２の目標追尾装置１００の動作は、クラスタの操作に関する処理であるステップＳ３０１〜Ｓ３０５が追加されている点が実施の形態１の動作とは異なっている。ここでは、実施の形態１との違いについて説明する。図６は、図５の目標追尾装置１００の動作を示すフローチャートである。なお、目標追尾装置１００の各機能は、図６のステップＳ１０２〜１１０，Ｓ２００，Ｓ３０２〜Ｓ３０６の処理をクラスタ毎に行う。
【００５９】
図６において、観測値割当部１１１がセンサ２００から新たに観測値を取得すると（ステップＳ１０１）、予測値算出部１５２、観測値割当部１１１及びコスト算出部１１２がステップＳ１０２〜Ｓ１０４のそれぞれの処理を行う。そして、クラスタ生成統合分離部１６１は、クラスタ操作を行う（ステップＳ３０１）。
【００６０】
具体的に、クラスタ生成統合分離部１６１は、複数のクラスタに存在する航跡候補が同一観測値と相関がとれた場合には、観測値を共有する航跡を有するクラスタを統合し一つのクラスタとする。また、クラスタ生成統合分離部１６１は、どのクラスタの航跡とも相関がとれない観測値について、新たにクラスタを生成する。さらに、クラスタ生成統合分離部１６１は、クラスタ内で観測値を取り合わない航跡が生成された場合には、各航跡を別のクラスタに分離する。
【００６１】
この後、多次元相関行列生成部１２１は、当該観測フレームにおけるクラスタ操作を判別する（ステップＳ３０２）。このときに、クラスタ操作が「操作なし」、「クラスタ生成」又は「クラスタ分離」である場合には、目標追尾装置１００の各機能は、それぞれのクラスタ内で、実施の形態１と同様にステップＳ１０５〜Ｓ１１０・Ｓ２００の処理を行う。
【００６２】
他方、クラスタ操作が「クラスタ統合」である場合には、多次元相関行列生成部１２１は、統合時フレーム数ルール設定部１７１のルール設定に応じて、当該観測フレームにおける指定フレーム数を決定する（ステップＳ３０３）。そして、多次元相関行列生成部１２１は、多次元相関行列の設定処理を行って、指定フレーム数に応じた多次元相関行列を、先の図３に示す観測値割当処理を実行して生成する（ステップＳ３０４）。そして、多次元相関行列生成部１２１は、統合時フレーム数ルール設定部１７１のルールが「可変」であるか「固定」であるかを確認する（ステップＳ３０５）。
【００６３】
ここで、統合時フレーム数ルール設定部１７１は、多次元の相関決定を実行するための目標追尾装置１００（計算機）の演算負荷の許容量が大きい場合には「可変」を設定し、その許容量が少ない場合には「固定」を設定する。具体的に、統合時フレーム数ルール設定部１７１は、多次元の相関決定を実行するための目標追尾装置１００の演算負荷の許容量と、事前に設定されたルール設定用パラメータとを用いた演算により、ルールを「固定」及び「可変」のいずれかに設定する。なお、観測値が存在していることがクラスタ統合前提条件であることから、クラスタ統合が発生した場合に、当該観測フレームが探知抜けフレームとなることはない。
【００６４】
そして、ルールが「固定」の場合、多次元相関行列生成部１２１は、多次元相関行列において同一の観測値が割り当てられている複数の航跡候補が存在する場合に、その中から最高尤度の航跡候補を残して、他の航跡を削除する（ステップＳ３０６）。そして、目標追尾装置１００の各機能は、ステップＳ１０７〜ステップＳ１０９のコスト決定から航跡決定までの処理を実行し、ステップＳ１１０へ移行する。
【００６５】
他方、ルールが「可変」の場合、多次元相関行列生成部１２１は、多次元相関行列において同一の観測値が複数の航跡候補に割り当てられることはないことから、目標追尾装置１００の各機能は、ステップＳ１０７〜ステップＳ１０９のコスト決定から航跡決定までの処理を実行し、ステップＳ１１０へ移行する。
【００６６】
従って、図７，８に示すように、目標番号が付与された２つの航跡候補群で最新フレームの観測値の取り合いが発生した際には、ルールが「固定」の場合と、ルールが「可変」とで、異なる処理がされる。そして、ステップＳ１１０では、平滑値算出部１５１は、最終的に残った航跡候補の当該観測フレームにおける平滑値を算出する。
【００６７】
上記のような実施の形態２の目標追尾装置によれば、多次元相関行列生成部１２１は、クラスタ生成統合分離部１６１によるクラスタ生成・分離の際には、クラスタに振り分けられた航跡候補と、その航跡候補についてのコストとを含む多次元相関行列を生成する。また、多次元相関行列生成部１２１は、クラスタ生成統合分離部１６１によるクラスタ統合の際には、統合される複数のクラスタのうち最も長い観測フレーム数を用いているクラスタの観測フレーム数に対応させて、多次元相関行列を生成する。この構成により、各クラスタの探知抜けフレームの発生状況の違いからクラスタ統合において多次元相関行列を生成する場合や、前回からの航跡候補を全て残す場合等の多次元相関行列のフレーム数を一時的に長くする必要が生じた場合に、その状況に対応した多次元相関行列を生成することができる。
【００６８】
また、多次元相関行列生成部１２１は、クラスタ生成統合分離部１６１によるクラスタ統合の際には、統合時フレーム数が所定値となるように、統合される複数のクラスタのそれぞれで多次元相関行列を生成し、クラスタ生成統合分離部１６１によるクラスタ統合とともに当該複数のクラスタのそれぞれの多次元相関行列を結合し、統合時フレーム数の多次元相関行列を生成する。この構成により、目標追尾装置１００の許容負荷の制限によって長いフレーム数の多次元の相関決定を実行できない場合には、多次元相関行列の生成時に前回からの航跡候補が削除されることから、目標追尾装置１００の許容負荷に対応する多次元相関行列を生成することができる。
【００６９】
さらに、統合時フレーム数ルール設定部１７１は、多次元の相関決定を実行するための目標追尾装置１００の演算負荷の許容量と、ルール設定用パラメータとを用いた演算により、ルールを「固定」及び「可変」のいずれかに設定する。そして、多次元相関行列生成部１２１は、統合時フレーム数ルール設定部１７１により設定されたルールに対応する多次元相関行列を生成する。この構成により、多次元相関行列のフレーム数を一時的に長くする必要が生じた場合や、目標追尾装置１００の許容負荷の制限によって長いフレーム数の多次元の相関決定を実行できない場合の異なる状況に対応した多次元相関行列を生成することができ、目標追尾装置１００の演算負荷の許容量を超えることなく、その許容量の範囲内で最も良い精度の多次元の相関決定の結果を出力することができる。
【００７０】
なお、実施の形態２では、多次元相関行列生成部１２１が図７，８に示すような２つの方式の多次元相関行列の生成処理を実行可能となっていた。しかしながら、多次元相関行列生成部１２１は、いずれか一方の方式の多次元相関行列の生成処理のみを実行可能であってもよい。この場合の目標追尾装置１００の構成は、統合時フレーム数ルール設定部１７１を省略した構成となる。
【符号の説明】
【００７１】
１００  目標追尾装置、１１０  航跡候補生成処理系（航跡候補生成処理部）、１１１  観測値割当部、１１２  コスト算出部、１２１  多次元相関行列生成部、１３１  多次元相関決定部、１４１  航跡決定部、１５０  航跡候補追尾処理系（航跡候補追尾処理部）、１５１  平滑値算出部、１５２  予測値算出部、１６１  クラスタ生成統合分離部（クラスタ処理部）、１７１  統合時フレーム数ルール設定部、２００  センサ、３００  表示装置。
【図１】
2010286275.tif 000013
【図２】
2010286275.tif 000014
【図３】
2010286275.tif 000015
【図４】
2010286275.tif 000016
【図５】
2010286275.tif 000017
【図６】
2010286275.tif 000018
【図７】
2010286275.tif 000019
【図８】
2010286275.tif 000020
【図９】
2010286275.tif 000021
【図１０】
2010286275.tif 000022
【図１１】
2010286275.tif 000023