JP6084810B2 - 追尾処理装置、及び追尾処理方法 - Google Patents

Description

本発明は、物標の追尾処理に用いられる、追尾処理装置、及び追尾処理方法に関する。

物標から反射されたレーダ信号等を用いて、物標を追尾（ＴＴ：Target Tracking）する追尾装置が知られている（例えば、特許文献１〜３参照）。レーダ信号は、レーダアンテナを用いたスキャン動作によって得られる。レーダ信号は、複数の物標について、位置を特定するデータを含んでいる。得られたレーダ信号は、追尾装置へ出力される。追尾装置は、上記のデータを用いて、追尾中の物標の位置を特定する。追尾装置は、特定した位置等の情報を用いて、追尾中の物標の運動を推定する。

特許文献１に記載の装置は、複数の物標について、位置等を観測することで、探知データを取得する。この装置は、探知データ等を用いて、追尾物標の航跡を算出し、その後、当該航跡のデータを保持する。また、この装置は、探知データと航跡との関連付けを行う。

特許第３０４４２９６号公報（［請求項１］） 特許第４４８８１８５号公報 実開平５−２３１７２号公報

ところで、追尾装置は、追尾中の物標がどのような運動状態にあるかを、定期的に参照する。そして、追尾装置は、過去の算出結果と、最新時点のレーダ信号に含まれる情報と、を用いて、最新時点における追尾物標の運動状態を推定する。具体的には、例えば、追尾装置は、追尾中の物標について、過去の算出結果を基に、最新時点での位置を、推定位置として算出する。次いで、追尾装置は、最新時点のレーダ信号で特定される複数の物標のうち、上記推定位置に最も近い物標を、追尾物標であると判定する。追尾装置は、当該追尾対象について、観測されたデータを用いて、最新時点における追尾物標の運動を推定する。

しかしながら、追尾中の船舶の航行によって生じる波又はクラッタが、上記推定位置の近傍に発生する場合がある。この場合、追尾装置は、当該波又はクラッタを、追尾物標であると誤って判定するおそれがある。このような誤判定が生じた場合、追尾装置は、本来の物標とは異なる波又はクラッタについて、追尾処理を行うこととなり、正確な追尾処理を行うことができない。

また、追尾物標の位置については、予測誤差と、観測誤差とが生じる。予測誤差とは、追尾処理において、上記推定位置と、追尾物標の実際の位置とのずれ量をいう。また、観測誤差とは、追尾物標について、レーダ信号を基に特定される位置と、実際の位置とのずれ量をいう。そこで、上記予測誤差を、予め設定したモデルに従って導出し、且つ、上記観測誤差を、予め設定したモデルに従って導出することが考えられる。しかしながら、このような構成では、予測誤差と、観測誤差のそれぞれについて、モデル化が必要であり、実用に耐えるモデルの作成は、非常に困難となる。

本発明は、上記実情に鑑みることにより、複数の物標のなかから、追尾対象を、確実且つ簡易に検出することができる、追尾処理装置、及び追尾処理方法を提供することを目的とする。

（１）上記課題を解決するために、この発明のある局面に係わる追尾処理装置は、追尾対象として選択された物標を追尾するための追尾処理装置であって、検出部と、尤度判定部と、運動推定部と、面積判定部と、を備えている。前記検出部は、複数の物標について、代表点の情報を含む特徴情報を検出するように構成されている。前記尤度判定部は、複数の前記代表点の少なくとも一部について、それぞれ、前記追尾対象の代表点に該当する尤度を算出するように構成されている。前記運動推定部は、前記追尾対象として特定された前記物標の前記特徴情報を用いて、前記追尾対象の運動状態を推定するように構成されている。前記面積判定部は、面積を判定するように構成されている。前記尤度判定部は、算出された前記尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定する。

さらに、前記検出部は、複数の前記物標の少なくとも一部について、前記代表点の情報以外の特徴情報を、追加の特徴情報として検出する。前記尤度判定部は、追加の前記特徴情報について、前記追尾対象の特徴情報に該当する尤度を算出し、且つ、当該尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定する。前記検出部は、送信信号に対する前記追尾対象を含む複数の前記物標での反射波であるエコー信号を検出し、且つ、当該エコー信号によって特定されるエコー像を検出し、且つ、当該エコー像の面積を追加の前記特徴情報として検出する。前記面積判定部は、前記追尾対象の前記エコー像の前記面積が所定のしきい値を超えているか否かを判定する。前記運動推定部は、前記面積が所定のしきい値以下である場合、前記運動状態を推定する時点での前記追尾対象の前記特徴情報を用いて、前記運動状態を算出し、且つ、前記面積が所定のしきい値を超えている場合、前記運動状態を算出する時点での前記追尾対象の前記特徴情報を用いることなく、前記運動状態を算出する。

（２）より好ましくは、前記検出部は、送信信号に対する各前記物標での反射波であるエコー信号を用いて、各前記特徴情報を検出する。追加の前記特徴情報は、前記エコー信号によって特定されるエコー像の形状情報、前記エコー信号のレベル、各前記物標の周囲の状態を特定する情報、及び、前記エコー信号に関するドップラーシフト量の少なくとも一つを含む。

（３）好ましくは、各前記特徴情報は、数値によって表される情報であり、前記尤度判定部は、各前記特徴情報毎に、確率密度関数を用いて、各前記尤度を算出する。

（４）より好ましくは、前記確率密度関数は、各前記特徴情報毎に設定される。各前記確率密度関数において平均の値として用いられる値は、対応する特徴情報について、尤度が算出される前に予め推定された値である。前記尤度判定部は、各確率密度関数において、対応する前記特徴情報を独立変数として用いる。また、前記尤度判定部は、前記独立変数に従属する従属変数を、前記尤度として算出する。

（５）より好ましくは、前記尤度判定部は、同一物標についての前記尤度を互いに乗算することで、総合尤度を算出する。また、前記尤度判定部は、各前記物標のうち、前記総合尤度が最も高い前記物標を、前記追尾対象として特定する。

（６）より好ましくは、前記面積判定部は、確率密度関数を用いて前記面積を判定するように構成されている。前記確率密度変数の平均の値として用いられる値は、前記面積について、前記面積判定部で判定される前に予め推定された値である。前記所定のしきい値は、前記平均に基づいて設定される。

（７）上記課題を解決するために、この発明のある局面に係わる追尾処理方法は、追尾対象として選択された物標を追尾するための追尾処理方法であって、検出ステップと、尤度判定ステップと、運動推定ステップと、面積判定ステップと、を含む。前記検出ステップでは、複数の物標について、代表点の情報を含む特徴情報を検出する。前記尤度判定ステップでは、複数の前記代表点の少なくとも一部について、それぞれ、前記追尾対象の代表点に該当する尤度を算出する。前記運動推定ステップでは、前記追尾対象として特定された前記物標の前記特徴情報を用いて、前記追尾対象の運動状態を推定する。前記面積判定ステップでは、面積を判定する。前記尤度判定ステップでは、算出された前記尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定する。さらに、前記検出ステップでは、複数の前記物標の少なくとも一部について、前記代表点の情報以外の特徴情報を、追加の特徴情報として検出する。前記尤度判定ステップでは、追加の前記特徴情報について、前記追尾対象の特徴情報に該当する尤度を算出し、且つ、当該尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定する。前記検出ステップでは、送信信号に対する前記追尾対象を含む複数の前記物標での反射波であるエコー信号を検出し、且つ、当該エコー信号によって特定されるエコー像を検出し、且つ、当該エコー像の面積を追加の前記特徴情報として検出する。前記面積判定ステップでは、前記追尾対象の前記エコー像の前記面積が所定のしきい値を超えているか否かを判定する。前記運動推定ステップでは、前記面積が所定のしきい値以下である場合、前記運動状態を推定する時点での前記追尾対象の前記特徴情報を用いて、前記運動状態を算出し、且つ、前記面積が所定のしきい値を超えている場合、前記運動状態を算出する時点での前記追尾対象の前記特徴情報を用いることなく、前記運動状態を算出する。

本発明によると、複数の物標のなかから、追尾物標を、確実、且つ簡易に検出することができる、追尾処理装置、及び追尾処理方法を提供することができる。

本発明の実施形態に係る追尾処理装置を含む、レーダ装置のブロック図である。 自船と、物標エコー像との関係を説明するための模式的な平面図である。 物標エコー像（物標）に関して抽出されるデータを説明するためのデータ一覧表である。 エコー検出部で検出された、物標エコー像を示す模式的な平面図である。 関数更新部における、確率密度関数更新についての処理の流れの一例を説明するためのフローチャートである。 代表点について、関数更新部で更新される関数の一例を示すグラフ図である。 追尾物標の面積について、関数更新部で更新される関数の一例を示すグラフ図である。 尤度算出部における処理の流れの一例を説明するためのフローチャートである。 面積判定部での処理の流れの一例を説明するためのフローチャートである。 物標エコー像と、他の物標エコー像とについて、融合が生じていない場合を説明するための模式図である。 物標エコー像と、他の物標エコー像とについて、融合が生じている場合を説明するための模式図である。 物標エコー像に過大なぼやけが生じている状態を示す模式図である。 運動推定部による処理の流れの一例を説明するためのフローチャートである。 運動推定部での処理を説明するための模式図である。 運動推定部での処理を説明するための模式図である。

以下、本発明を実施するための形態について、図面を参照しつつ説明する。本発明は、追尾対象として選択された物標を追尾する、追尾処理装置として広く適用することができる。以下では、追尾対象として設定された物標を「追尾物標」という。また、以下では、図中同一または相当部分には、同一符号を付してその説明は繰り返さない。

図１は、本発明の実施形態に係る追尾処理装置３を含む、レーダ装置１のブロック図である。本実施形態のレーダ装置１は、例えば、漁船等の船舶に備えられる舶用レーダである。レーダ装置１は、主に他船等の物標の探知に用いられる。また、レーダ装置１は、追尾物標として選択された物標を追尾することが可能に構成されている。レーダ装置１は、複数の追尾物標を、同時に追尾可能に構成されている。レーダ装置１は、追尾物標の運動状態を推定するように構成されている。本実施形態では、レーダ装置１は、上記運動状態として、追尾物標の推定速度ベクトルを算出する。推定速度ベクトルとは、追尾物標について推定される、進行方向及び進行速度を示すベクトルである。レーダ装置１は、追尾物標の推定速度ベクトルを、画面に表示する。尚、以下では、レーダ装置１が備えられている船舶を「自船」という。

図１に示すように、レーダ装置１は、アンテナユニット２と、追尾処理装置３と、表示器４と、を備えている。

アンテナユニット２は、アンテナ５と、受信部６と、Ａ／Ｄ変換部７と、を含んでいる。

アンテナ５は、指向性の強いパルス状電波を送信可能なレーダアンテナである。また、アンテナ５は、物標からの反射波であるエコー信号を受信するように構成されている。即ち、物標のエコー信号は、アンテナ５からの送信信号に対する、物標での反射波である。レーダ装置１は、パルス状電波を送信してからエコー信号を受信するまでの時間を測定する。これにより、レーダ装置１は、物標までの距離ｒを検出することができる。アンテナ５は、水平面上で３６０°回転可能に構成されている。アンテナ５は、パルス状電波の送信方向を変えながら（アンテナ角度を変えながら）、電波の送受信を繰り返し行うように構成されている。以上の構成で、レーダ装置１は、自船周囲の平面上の物標を、３６０°にわたり探知することができる。

なお、以下の説明では、パルス状電波を送信してから次のパルス状電波を送信するまでの動作を「スイープ」という。また、電波の送受信を行いながらアンテナを３６０°回転させる動作を「スキャン」と呼ぶ。以下では、あるスキャンのことを、「ｎスキャン」といい、ｎスキャンから１個前のスキャンのことを、「ｎ−１スキャン」という。同様に、ｎスキャンからｍ個前のスキャンのことを、「ｎ−ｍスキャン」という。尚、ｎ、ｍは、何れも自然数である。また、以下では、動作について特に説明なき場合、ｎスキャン時点におけるレーダ装置１の動作を説明する。

受信部６は、アンテナ５で受信したエコー信号を検波して増幅する。受信部６は、増幅したエコー信号を、Ａ／Ｄ変換部７へ出力する。Ａ／Ｄ変換部７は、アナログ形式のエコー信号をサンプリングし、複数ビットからなるデジタルデータ（エコーデータ）に変換する。ここで、上記エコーデータは、アンテナ５が受信したエコー信号の強度（信号レベル）を特定するデータを含んでいる。Ａ／Ｄ変換部７は、エコーデータを、追尾処理装置３へ出力する。

追尾処理装置３は、複数の物標の中から選択された物標を、追尾物標として特定し、且つ、当該追尾物標の追尾処理を行うように構成されている。より具体的には、追尾処理装置３は、追尾物標の推定速度ベクトル、及び追尾物標の推定位置等を算出するように構成されている。

追尾処理装置３は、ＣＰＵ、ＲＡＭ及びＲＯＭ（図示せず）等を含むハードウェアを用いて構成されている。また、追尾処理装置３は、ＲＯＭに記憶された追尾処理プログラムを含む、ソフトウェアを用いて構成されている。

上記追尾処理プログラムは、本発明に係る追尾処理方法を、追尾処理装置３に実行させるためのプログラムである。上記ハードウェアとソフトウェアとは、協働して動作するように構成されている。これにより、追尾処理装置３を、スイープメモリ８、信号処理部９、エコー検出部１０、及び追尾処理部１１等として機能させることができる。追尾処理装置３は、１スキャン毎に、以下に説明する処理を行うように構成されている。

追尾処理装置３は、信号処理部９と、エコー検出部（検出部）１０と、追尾処理部１１と、を有している。

信号処理部９は、フィルタ処理等を施すことにより、エコーデータに含まれる干渉成分と、不要な波形データと、を除去する。また、信号処理部９は、物標エコー像に関するエコーデータの特徴情報の検出を、行うように構成されている。信号処理部９は、処理したエコーデータを、エコー検出部１０へ出力する。

エコー検出部１０は、物標エコー像の検出と、物標エコー像に関するエコーデータの特徴情報の検出とを、行うように構成されている。即ち、エコー検出部１０は、物標エコー像検出部と、特徴情報抽出部と、を含んでいる。信号処理部９及びエコー検出部１０によって、物標の特徴情報を検出するための検出部が、構成されている。

エコー検出部１０は、信号処理部９からエコーデータを読み出すときの読出しアドレスに基づいて、当該エコーデータに対応する位置までの距離ｒを求める。また、アンテナ５からエコー検出部１０へは、当該アンテナ５が現在どの方向を向いているか（アンテナ角度θ）を示すデータが出力されている。以上の構成で、エコー検出部１０は、エコーデータを読み出す際には、当該エコーデータに対応する位置を、距離ｒとアンテナ角度θとの極座標で取得することができる。

エコー検出部１０は、エコーデータに対応する位置に物標が存在するか否かを検出するように構成されている。エコー検出部１０は、例えば、エコーデータに対応する位置の信号レベル、即ち、信号強度を判別する。エコー検出部１０は、信号レベルが所定のしきいレベル値以上である位置には、物標が存在していると判別する。

次いで、エコー検出部１０は、物標が存在している範囲を検出する。エコー検出部１０は、例えば、物標が存在している一まとまりの領域を、物標エコー像が存在している領域として検出する。このようにして、エコー検出部１０は、エコーデータを基に、物標エコー像を検出する。この物標エコー像の外郭形状は、物標の外郭形状と略合致する。但し、エコーデータに含まれるノイズ等に起因して、この物標エコー像の外郭形状と、物標の外郭形状とは、わずかに異なる。次に、エコー検出部１０は、エコーデータを用いて、物標エコー像に関連する特徴情報を抽出する。

図２は、自船１００と、物標エコー像１２０との関係を説明するための模式的な平面図である。図２では、物標エコー像１２０は、矩形の像として例示されている。また、図２では、物標エコー像１２０によって特定される物標１３０が示されている。図２では、物標１３０の外郭形状は、物標エコー像１２０と合致した状態で表示されている。

図１及び図２に示すように、極座標系では、自船１００の位置としての自船位置Ｍ１を基準として、自船位置Ｍ１からの直線距離が、距離ｒとして示され、自船位置Ｍ１周りの角度が、角度θとして示される。本実施形態では、自船位置Ｍ１は、アンテナ５の位置に相当する。エコー検出部１０は、物標エコー像１２０の代表点Ｐの抽出に際しては、自船位置Ｍ１を中心とする、リング状部分の一部形状の像１１０を用いる。この像１１０は、第１直線１１１、第２直線１１２、第１円弧１１３、及び第２円弧１１４によって囲まれた領域の像である。

第１直線１１１は、物標エコー像１２０の後縁１２０ａのうち自船位置Ｍ１に最も近い点と、自船位置Ｍ１と、を通る直線である。第２直線１１２は、物標エコー像１２０の前縁１２０ｂのうち自船位置Ｍ１に最も近い点と、自船位置Ｍ１と、を通る直線である。第１円弧１１３は、物標エコー像１２０のうち自船位置Ｍ１から最も近い部分１２０ｃを通る円弧である。第１円弧１１３の曲率中心点は、自船位置Ｍ１である。第２円弧１１４は、物標エコー像１２０のうち自船位置Ｍ１から最も遠い部分１２０ｄを通る円弧である。第２円弧１１４は、第１円弧１１３と同心である。

図３は、物標エコー像１２０（物標１３０）に関して抽出されるデータを説明するためのデータ一覧表である。図２及び図３に示すように、本実施形態では、信号処理部９及びエコー検出部１０は、協働して、物標エコー像１２０（物標１３０）について、下記１２個のデータを、特徴情報データとして抽出する。即ち、エコー検出部１０は、エコーデータ、及び物標エコー像１２０の画像データを基に、下記１２個のテキストデータを抽出する。本実施形態では、１２個のテキストデータとは、フラグのデータ２０１と、距離ｒｐのデータ２０２と、終了角度θｅのデータ２０３と、角度幅θｗのデータ２０４と、最前縁距離ｒｎのデータ２０５と、最後縁距離ｒｆのデータ２０６と、面積ａｒ（形状情報）のデータ２０７と、代表点Ｐの座標データ２０８と、エコーレベルｅｃのデータ２０９と、隣接距離（物標の周囲の状態を特定する情報）ａｄのデータ２１０と、ドップラーシフト量ｄｓのデータ２１１と、時刻ｔｍのデータ２１２と、である。

上記のフラグは、例えば、物標エコー像１２０について、所定の状態であるか否か、等を特定するフラグである。このフラグは、「１」又は「０」等に設定されるように構成されている。

距離ｒｐは、自船位置Ｍ１から物標エコー像１２０の代表点Ｐまでの直線距離である。本実施形態では、代表点Ｐは、像１１０の図心点である。終了角度θｅは、物標エコー像１２０の検出が終わった時点における、前述のアンテナ角度θである。角度幅θｗは、物標エコー像１２０について、自船位置Ｍ１回りの角度方向の幅である。角度幅θｗは、第１直線１１１と第２直線１１２とがなす角度でもある。最前縁距離ｒｎは、物標エコー像１２０の部分１２０ｃと、自船位置Ｍ１との距離である。最後縁距離ｒｆは、物標エコー像１２０の部分１２０ｄと、自船位置Ｍ１との距離である。面積ａｒは、リング状部分の一部形状の像１１０における面積であり、本実施形態では、物標エコー像１２０の面積として扱われる。

エコーレベルｅｃは、物標エコー像１２０を特定するエコー信号の、強度を示している。この強度は、物標エコー像１２０を特定するエコー信号のピーク強度であってもよいし、当該エコー信号の強度の平均値であってもよい。本実施形態では、隣接距離ａｄは、例えば、隣接する２つの物標エコー像１２０間の距離である。ドップラーシフト量ｄｓは、例えば、アンテナ５から放射されたパルス信号の周波数と、物標エコー像１２０によって特定される物標１３０で反射したエコー信号の周波数と、の差である。ドップラーシフト量ｄｓを基に、物標エコー像１２０によって特定される物標１３０と、自船１００との相対速度を求めることが可能である。時刻ｔｍは、物標エコー像１２０を検出した時点の時刻である。尚、エコー像１２０のデータは、予備のデータ領域を含んでいてもよい。本実施形態では、信号処理部９及びエコー検出部１０は、各物標エコー像１２０について、上記１２個の特徴情報を抽出する。これら１２個の特徴情報のうち、代表点Ｐの位置以外の特徴情報は、本発明の「追加の特徴情報」の一例である。また、これら１２個の特徴情報は、何れも、数値によって表される情報である。

エコー検出部１０で検出された複数の物標エコー像１２０の一例を、図４に示している。図４は、エコー検出部１０で検出された、物標エコー像１２０を示す模式的な平面図である。図４では、一例として、ｎスキャン時点における４つの物標エコー像１２０（１２１，１２２，１２３，１２４）を示している。図４では、物標エコー像１２０（１２１，１２２，１２３，１２４）の形状は、それぞれ、物標１３０（１３１，１３２，１３３，１３４）の形状と合致している。

物標エコー像１２１によって特定される物標１３１，物標エコー像１２２によって特定される物標１３２，及び物標エコー像１２３によって特定される物標１３３は、例えば、小型船舶である。物標エコー像１２４によって特定される物標１３４は、例えば、大型船舶である。エコー検出部１０は、ｎスキャン時点において、代表点Ｐとして、物標エコー像１２１の代表点Ｐ１（ｎ）と、物標エコー像１２２の代表点Ｐ２（ｎ）と、物標エコー像１２３の代表点Ｐ３（ｎ）と、物標エコー像１２４の代表点Ｐ４（ｎ）と、を検出している。以下では、物標１３１が、追尾物標１４０である場合を例に説明する。

図１〜図４に示すように、エコー検出部１０は、各物標エコー像１２０についての特徴情報データ２０１〜２１２を、追尾処理部１１へ出力する。

追尾処理部１１は、複数の物標１２０のなかから、追尾物標１４０を特定し、且つ、当該追尾物標１４０の追尾処理を行うように構成されている。追尾物標１４０は、例えば、表示器４に表示された、複数の物標１３０を示すシンボル等を基に、オペレータによって選択される。オペレータによる、追尾物標１４０の選択指令は、例えば、オペレータが操作装置（図示せず）を操作することにより、発せられる。本実施形態では、追尾処理部１１は、特に説明無き限り、Ｘ−Ｙ座標系を基準に、処理を行うように構成されている。

追尾処理部１１は、ｎスキャン時点（最新のスキャン時点）における、追尾物標１４０の推定速度ベクトルＶ１（ｎ）を算出するように構成されている。また、追尾処理部１１は、推定速度ベクトルＶ１（ｎ）を表示器４に表示させるように、構成されている。

追尾処理部１１は、特徴情報メモリ１２と、選別部１３と、関連付け部（尤度判定部）１４と、面積判定部１５と、運動推定部１６と、を有している。

特徴情報メモリ１２は、信号処理部９及びエコー検出部１０から出力されたデータを蓄積するように構成されている。また、特徴情報メモリ１２は、後述する関数ｆ１〜ｆ５に関連するデータを蓄積するように構成されている。特徴情報メモリ１２は、（ｎ−Ｔ）スキャン時点から、ｎスキャン時点までの各時点における、全ての物標エコー像１２０について、特徴情報データ２０１〜２１２を蓄積している。尚、定数Ｔは、予め設定されている値であり、例えば、数十程度である。

選別部１３は、選別領域Ｓ１（ｎ）を設定するように構成されている。選別領域Ｓ１（ｎ）は、ｎスキャン時点における、追尾物標１４０の代表点（追尾代表点）ＰＲ（ｎ）を選別するための領域である。選別部１３は、エコー検出部１０、関連付け部１４、及び運動推定部１６に接続されている。

選別部１３は、エコー検出部１０から、ｎスキャン時点における、各物標エコー像１２０の代表点Ｐ（Ｐ１（ｎ）〜Ｐ４（ｎ））の座標データを受け付ける。また、選別部１３は、（ｎ−１）スキャン時点を算出基準時点とする、追尾物標１４０の推定速度ベクトルＶ１（ｎ−１）のデータを、受け付ける。本実施形態では、選別部１３は、運動推定部１６から、上記推定速度ベクトルＶ１（ｎ−１）のデータを受け付ける。

次に、選別部１３は、追尾物標１４０の推定位置ＰＲｅ（ｎ−１）を検出する。この推定位置ＰＲｅ（ｎ−１）は、推定速度ベクトルＶ１（ｎ−１）によって特定される位置である。推定位置ＰＲｅ（ｎ−１）は、ｎスキャン時点において、追尾物標１４０の代表点ＰＲ（ｎ）が到達すると推定される位置である。選別部１３は、この推定位置ＰＲｅ（ｎ−１）を中心にして、所定の半径を有する選別領域Ｓ１（ｎ）を設定する。図４は、選別領域Ｓ１内に、物標エコー像１２１，１２２の代表点Ｐ１（ｎ），Ｐ２（ｎ）が含まれた状態を例示している。選別部１３は、代表点Ｐ１（ｎ），Ｐ２（ｎ）のそれぞれについて、座標データ２０８を、関連付け部１４へ出力する。

関連付け部１４は、推定位置ＰＲｅ（ｎ−１）と、ｎスキャン時点における追尾物標１４０の代表点ＰＲ（ｎ）と、を関連付ける処理を行う。本実施形態では、関連付け部１４は、選別領域Ｓ１（ｎ）内の代表点Ｐ１（ｎ），Ｐ２（ｎ）のそれぞれについて、追尾物標１４０の代表点ＰＲ（ｎ）に該当する尤度Ｌｈ１００，Ｌｈ２００を算出する。関連付け部１４は、この尤度Ｌｈ１００，Ｌｈ２００に基づいて、ｎスキャン時点における追尾物標１４０代表点ＰＲ（ｎ）を特定する。

関連付け部１４は、関数更新部２１と、尤度算出部２２と、尤度比較部２３と、を有している。

関数更新部２１は、関連付け処理に用いる関数を更新するように構成されている。この関数として、本実施形態では、確率密度関数（ガウス関数、正規分布関数）が用いられる。本実施形態において、確率密度関数は、平均、及び分散を有する関数である。本実施形態では、確率密度関数は、過去の追尾結果を基に更新される。

本実施形態では、確率密度関数に用いられる特徴情報として、代表点ＰＲ（ｎ）の位置情報と、追尾物標１４０についての面積ａｒの情報と、追尾物標１４０についてのエコーレベルｅｃの情報と、追尾物標１４０についての隣接距離ａｄの情報と、追尾物標１４０につてのドップラーシフト量ｄｓの情報と、が挙げられる。本実施形態では、１つの確率密度関数ｐｄｆが設定されている。

確率密度関数ｐｄｆは、代表点ＰＲの位置の次元の関数ｆ１と、面積ａｒの次元の関数ｆ２と、エコーレベルｅｃの次元の関数ｆ３と、隣接距離ａｄの次元の関数ｆ４と、ドップラーシフト量ｄｓの次元の関数ｆ５と、によって特定される関数である。関数ｆ１〜ｆ５は、それぞれ、確率密度関数であり、確率密度関数ｐｄｆは、多次元の確率密度関数である。関数更新部２１は、関数ｆ１〜ｆ５を更新することで、確率密度関数ｐｄｆを更新する。尚、以下では、各上記関数ｆ１〜ｆ５を総称していう場合に、単に関数ｆという場合がある。

以下では、関数更新部２１における、確率密度関数ｐｄｆを更新する処理の流れの一例を説明する。図５は、関数更新部２１における、確率密度関数ｐｄｆの更新についての処理の流れの一例を説明するためのフローチャートである。追尾処理装置３は、以下に示すフローチャートの各ステップを、図示しないメモリから読み出して実行する。このプログラムは、外部からインストールできる。このインストールされるプログラムは、例えば記録媒体に格納された状態で流通する。追尾処理装置３での処理は、各スキャン時点毎に行われる。また、本実施形態では、フローチャートについて説明する場合には、フローチャート以外の図面についても、適宜参照する。

本実施形態では、関数更新部２１は、最尤推定処理を行う。具体的には、図５に示すように、関数更新部２１は、確率密度関数ｐｄｆを更新するための特徴情報を選択する（ステップＳ１０１）。関数更新部２１は、例えば、まず、代表点ＰＲの位置情報を選択する。即ち、関数更新部２１は、代表点ＰＲの位置に関する関数ｆ１を更新する。

図６は、代表点ＰＲについて、関数更新部２１で更新される関数ｆ１の一例を示すグラフ図である。関数ｆ１は、平均μ１１，μ１２、及び分散σ１１^２，σ１２^２を有している。本実施形態において、平均は、各関数ｆにおける推定値に相当する。即ち、本実施形態の各関数ｆにおいて、平均の値として用いられる値は、対応する特徴情報について、尤度Ｌｈが算出される前に予め推定された値である。図６のグラフのｘ軸は、ｘ座標を示しており、ｙ軸は、ｙ座標を示している。また、図６のグラフの縦軸（ｚ軸）は、尤度Ｌｈを示している。

関数更新部２１は、関数ｆ１における平均μ１１，μ１２を算出する（ステップＳ１０２）。具体的には、関数更新部２１は、追尾物標１４０の代表点ＰＲ（ｎ−１）について、ｘ座標の平均μ１１、及びｙ座標の平均μ１２を算出する。

本実施形態では、平均μ１１は、下記式（１）によって算出される。

尚、Δｄｘ_ｉは、任意のｉスキャン時点における、代表点ＰＲ（ｉ）のｘ座標の予測誤差である。より具体的には、Δｄｘ_ｉは、ｉスキャン時点における、代表点ＰＲ（ｉ）のｘ座標の（観測値−予測値）である。この場合の観測値とは、ｉスキャン時点でエコー検出部１０によって検出された、代表点ＰＲ（ｉ）のｘ座標である。また、予測値とは、（ｉ−１）スキャン時点において運動推定部１６で予測された、代表点ＰＲ（ｉ）のｘ座標である。したがって、平均μ１１は、（ｎ−Ｔ）スキャン時点から（ｎ−１）スキャン時点までの各スキャン時点における、代表点ＰＲのｘ座標の予測誤差を平均した値である。

平均μ１２は、平均μ１１と同様にして算出される。即ち、平均μ１２は、（ｎ−Ｔ）スキャン時点から（ｎ−１）スキャン時点までの各スキャン時点における、代表点ＰＲのｙ座標の予測誤差を平均した値である。

次に、関数更新部２１は、関数ｆ１の分散σ１１^２，σ１２^２を算出する（ステップＳ１０３）。即ち、関数更新部２１は、関数ｆ１のｘ軸方向についての分散σ１１^２と、ｙ軸方向についての分散σ１２^２と、を算出する。

本実施形態では、分散σ１１^２は、下記式（２）によって算出される。

分散σ１２^２は、平均μ１１と同様にして算出される。即ち、分散σ１２^２は、式（２）のΔｄｘ_ｉをΔｄｙ_ｉに置き換え、且つ、μ１１をμ１２に置き換えた式によって算出される。このように、関数更新部２１は、平均μ１１，μ１２、及び分散σ１１^２，σ１２^２を算出する。関数更新部２１は、平均μ１１，μ１２、及び分散σ１１^２，σ１２^２を算出することで、関数ｆ１を更新する。

次に、関数更新部２１は、関数ｆ１〜ｆ５の全てが更新されたか否かを判定する（ステップＳ１０４）。この場合、関数ｆ２，ｆ３，ｆ４ｆ，ｆ５は、未だ更新されていない（ステップＳ１０４でＮＯ）。よって、関数更新部２１は、ステップＳ１０１〜Ｓ１０４の処理を繰り返す。

具体的には、関数更新部２１は、確率密度関数ｐｄｆを更新するための特徴情報を選択する（ステップＳ１０１）。この場合、関数更新部２１は、例えば、特徴情報として、面積ａｒの情報を選択する。即ち、関数更新部２１は、関数ｆ２を更新する。

図７は、ｎスキャン時点における追尾物標１４０の面積ａｒについて、関数更新部２１で更新される関数ｆ２の一例を示すグラフ図である。図７に示すように、関数ｆ２は、平均μ２、分散σ２^２を有している。図７のグラフの横軸は、面積ａｒを示している。図７のグラフの縦軸は、尤度Ｌｈを示している。

関数更新部２１は、関数ｆ２における平均μ２を算出する（ステップＳ１０２）。平均μ２は、関数更新部２１によって、例えば、以下のようにして算出される。即ち、関数更新部２１は、まず、特徴情報メモリ１２を参照する。これにより、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの各スキャン時点について、追尾物標１４０として設定された物標１３１における面積ａｒのデータ２０７を参照する。そして、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの間における、上記面積ａｒの平均値を算出する。関数更新部２１は、この平均値を、平均μ２として設定する。尚、関数更新部２１は、（ｎ−１）スキャン時点での物標１３１についての面積ａｒを、平均μ２として設定してもよい。

次に、関数更新部２１は、関数ｆ２の分散σ２^２を算出する（ステップＳ１０３）。関数更新部２１は、例えば、最尤推定処理を行う。これにより、関数更新部２１は、関数ｆ２の分散σ２^２を算出する。このように、関数更新部２１は、平均μ２及び分散σ２^２を算出することで、関数ｆ２を更新する。

次に、関数更新部２１は、関数ｆ１〜ｆ５の全てが更新されたか否かを判定する（ステップＳ１０４）。この場合、関数ｆ３，ｆ４，ｆ５は、未だ更新されていない（ステップＳ１０４でＮＯ）。よって、関数更新部２１は、ステップＳ１０１〜Ｓ１０４の処理を繰り返す。

具体的には、関数更新部２１は、確率密度関数ｐｄｆを更新するための特徴情報を選択する（ステップＳ１０１）。この場合、関数更新部２１は、例えば、特徴情報として、エコーレベルｅｃを選択する。即ち、関数更新部２１は、関数ｆ３を更新する。尚、関数ｆ３は、図示していないけれども、関数ｆ２と同様の形状となる。

関数更新部２１は、関数ｆ３における平均μ３を算出する（ステップＳ１０２）。この場合のエコーレベルｅｃの平均μ３は、関数更新部２１によって、例えば、以下のようにして算出される。即ち、関数更新部２１は、まず、特徴情報メモリ１２を参照する。これにより、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの各スキャン時点について、追尾物標１４０として特定された部物標１３１のエコーレベルｅｃのデータ２０９を参照する。そして、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの間における、上記エコーレベルｅｃの平均値を算出する。関数更新部２１は、この平均値を、平均μ３として設定する。尚、関数更新部２１は、（ｎ−１）スキャン時点での物標１３１についてのエコーレベルｅｃを、平均μ３として設定してもよい。

次に、関数更新部２１は、関数ｆ３の分散σ３^２を算出する（ステップＳ１０３）。関数更新部２１は、例えば、最尤推定処理を行う。これにより、関数更新部２１は、関数ｆ３の分散σ３^２を算出する。このように、関数更新部２１は、平均μ３及び分散σ３^２を算出することで、関数ｆ３を更新する。

次に、関数更新部２１は、関数ｆ１〜ｆ５の全てが更新されたか否かを判定する（ステップＳ１０４）。この場合、関数ｆ４，ｆ５は、未だ設定されていない（ステップＳ１０４でＮＯ）。よって、関数更新部２１は、ステップＳ１０１〜Ｓ１０４の処理を繰り返す。

具体的には、関数更新部２１は、確率密度関数ｐｄｆを更新するための特徴情報を選択する（ステップＳ１０１）。この場合、関数更新部２１は、例えば、特徴情報として、隣接距離ａｄを選択する。即ち、関数更新部２１は、関数ｆ４を設定する。尚、関数ｆ４は、図示していないけれども、関数ｆ２と同様の形状となる。

関数更新部２１は、関数ｆ４における平均μ４を算出する（ステップＳ１０２）。平均μ４は、関数更新部２１によって、例えば、以下のようにして算出される。即ち、関数更新部２１は、まず、特徴情報メモリ１２を参照する。これにより、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの各スキャン時点について、追尾物標１４０として特定された物標１３１についての隣接距離ａｄのデータ２１０を参照する。そして、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの間における、上記隣接距離ａｄの平均値を算出する。関数更新部２１は、この平均値を、平均μ４として設定する。尚、関数更新部２１は、（ｎ−１）スキャン時点での物標１３１についての上記隣接距離ａｄを、平均μ４として設定してもよい。

次に、関数更新部２１は、関数ｆ４の分散σ４^２を算出する（ステップＳ１０３）。関数更新部２１は、例えば、最尤推定処理を行う。これにより、関数更新部２１は、関数ｆ４の分散σ４^２を算出する。このように、関数更新部２１は、平均μ４及び分散σ４^２を算出することで、関数ｆ４を更新する。

次に、関数更新部２１は、関数ｆ１〜ｆ５の全てが設定されたか否かを判定する（ステップＳ１０４）。この場合、関数ｆ５は、未だ更新されていない（ステップＳ１０４でＮＯ）。よって、関数更新部２１は、ステップＳ１０１〜Ｓ１０４の処理を繰り返す。

具体的には、関数更新部２１は、確率密度関数ｐｄｆを更新するための特徴情報を選択する（ステップＳ１０１）。この場合、関数更新部２１は、例えば、特徴情報として、ドップラーシフト量ｄｓを更新する。即ち、関数更新部２１は、関数ｆ５を更新する。尚、関数ｆ５については、図示していないけれども、関数ｆ２と同様となる。

関数更新部２１は、関数ｆ５における平均μ５を算出する（ステップＳ１０２）。この場合の隣接距離ａｄの平均μ５は、関数更新部２１によって、例えば、以下のようにして算出される。即ち、関数更新部２１は、まず、特徴情報メモリ１２を参照する。これにより、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの各スキャン時点について、追尾物標１４０として特定された物標１３１についてのドップラーシフト量ｄｓのデータ２１１を参照する。そして、関数更新部２１は、（ｎ−１）スキャン時点から、（ｎ−Ｔ）スキャン時点までの間における、上記ドップラーシフト量ｄｓの平均値を算出する。関数更新部２１は、この平均値を、平均μ５として設定する。尚、関数更新部２１は、（ｎ−１）スキャン時点での物標１３１についての上記ドップラーシフト量ｄｓを、平均μ５として設定してもよい。

次に、関数更新部２１は、関数ｆ５の分散σ５^２を算出する（ステップＳ１０３）。関数更新部２１は、例えば、最尤推定処理を行う。これにより、関数更新部２１は、関数ｆ５の分散σ５^２を算出する。このように、関数更新部２１は、平均μ５及び分散σ５^２を算出することで、関数ｆ５を更新する。

次に、関数更新部２１は、関数ｆ１〜ｆ５の全てが更新されたか否かを判定する（ステップＳ１０４）。この場合、関数ｆ１〜ｆ５の全てが更新されている（ステップＳ１０４でＹＥＳ）。よって、関数更新部２１は、関数ｆ１〜ｆ５によって特定される確率密度関数ｐｄｆのデータを、尤度算出部２２へ出力し（ステップＳ１０５）、処理を終了する。上記の構成により、関数ｆ１〜ｆ５は、それぞれ、独立して更新される。また、関数ｆ１〜ｆ５は、それぞれ、固定された関数ではなく、各スキャン時点毎に変化する関数である。

尤度算出部２２及び尤度比較部２３は、設定された確率密度関数ｐｄｆを基に、選別領域Ｓ１（ｎ）内の物標１３１，１３２の何れが追尾物標１４０であるかの判定処理を行う。以下では、尤度算出部２２及び尤度比較部２３における処理の流れの一例を説明する。図８は、尤度算出部２２及び尤度比較部２３における処理の流れの一例を説明するためのフローチャートである。

図８に示すフローチャートでは、尤度算出部２２は、選別領域Ｓ１（ｎ）の物標１３１，１３２毎に、追尾物標１４０に該当する尤度としての総合尤度Ｌｈ１００，２００を算出する。総合尤度Ｌｈ１００は、後述する尤度Ｌｈ１１，Ｌｈ１２，Ｌｈ１３，Ｌｈ１４，Ｌｈ１５の積である。総合尤度Ｌｈ２００は、後述する尤度Ｌｈ２１，Ｌｈ２２，Ｌｈ２３，Ｌｈ２４，Ｌｈ２５の積である。尤度比較部２３は、総合尤度Ｌｈ１００，２００に基づいて、物標１３１，１３２のなかから、追尾物標１４０を特定する。

具体的には、尤度算出部２２は、まず、尤度算出が行われる物標１３０を選択する（ステップＳ２０１）。本実施形態では、尤度算出部２２は、選別領域Ｓ１（ｎ）内の物標１３０（１３１，１３２）のうちの、物標１３１を選択する。次に、尤度算出部２２は、確率密度関数ｐｄｆの関数ｆ１〜ｆ５のなかから１つを選択する（ステップＳ２０２）。

この場合、尤度算出部２２は、まず、関数ｆ１を選択する。次に、尤度算出部２２は、物標１３１の代表点Ｐ１（ｎ）について、尤度Ｌｈ１１を算出する（ステップＳ２０３）。具体的には、尤度算出部２２は、関数ｆ１のデータと、物標１３１の代表点Ｐ１（ｎ）の座標データ２０８と、を参照する。そして、（ｎ−１）スキャン時点で推定された、代表点ＰＲ（ｎ）の予測値のｘ座標と、代表点Ｐ１（ｎ）のｘ座標との差を、独立変数として、関数ｆ１（図６参照）のｘ座標に代入する。また、（ｎ−１）スキャン時点で推定された、代表点ＰＲ（ｎ）の予測値のｙ座標と、代表点Ｐ１（ｎ）のｙ座標との差を、独立変数として、関数ｆ１（図６参照）のｙ座標に代入する。これにより、関数ｆ１における従属変数が、尤度Ｌｈ１１として算出される。即ち、物標１３１の代表点Ｐ１（ｎ）について、追尾物標１４０の代表点ＰＲ（ｎ）に該当する尤度Ｌｈ１１が、算出される。

次に、尤度算出部２２は、物標１３１に関する尤度Ｌｈ１１〜Ｌｈ１５の全てが算出されたか否かを判定する（ステップＳ２０４）。この場合、尤度Ｌｈ１２〜Ｌｈ１５は、未だ算出されていない（ステップＳ２０４でＮＯ）。よって、尤度算出部２２は、ステップＳ２０２〜Ｓ２０４の処理を繰り返す。

具体的には、尤度算出部２２は、確率密度関数ｐｄｆの関数ｆ１〜ｆ５のなかから１つを選択する（ステップＳ２０２）。

この場合、尤度算出部２２は、関数ｆ２を選択する。次に、尤度算出部２２は、物標１３１の面積ａｒについて、追尾物標１４０の面積ａｒに該当する尤度としての尤度Ｌｈ１２を算出する（ステップＳ２０３）。具体的には、尤度算出部２２は、関数ｆ２のデータと、物標１３１の面積ａｒのデータ２０７と、を参照する。そして、尤度算出部２２は、ｎスキャン時点における物標１３１の面積ａｒを、独立変数として、関数ｆ２に代入する。これにより、関数ｆ２における従属変数が、尤度Ｌｈ１２として算出される。即ち、物標１３１のエコーレベルｅｃについての尤度Ｌｈ１２が、算出される。

次に、尤度算出部２２は、物標１３１に関する尤度Ｌｈ１１〜Ｌｈ１５の全てが算出されたか否かを判定する（ステップＳ２０４）。この場合、尤度Ｌｈ１３〜Ｌｈ１５は、未だ算出されていない（ステップＳ２０４でＮＯ）。よって、尤度算出部２２は、ステップＳ２０２〜Ｓ２０４の処理を繰り返す。

具体的には、尤度算出部２２は、関数ｆ３〜ｆ５のなかから１つを選択する（ステップＳ２０２）。

この場合、尤度算出部２２は、関数ｆ３を選択する。次に、尤度算出部２２は、物標１３１のエコーレベルｅｃについて、追尾物標１４０のエコーレベルｅｃに該当する尤度としての尤度Ｌｈ１３を算出する（ステップＳ２０３）。具体的には、尤度算出部２２は、関数ｆ３のデータと、物標１３１のエコーレベルｅｃのデータ２０９と、を参照する。そして、尤度算出部２２は、ｎスキャン時点における物標１３１についてのエコーレベルｅｃを、独立変数として、関数ｆ３に代入する。これにより、関数ｆ３における従属変数が、尤度Ｌｈ１３として算出される。即ち、物標１３１についてのエコーレベルｅｃの尤度Ｌｈ１３が、算出される。

次に、尤度算出部２２は、物標１３１に関する尤度Ｌｈ１１〜Ｌｈ１５の全てが算出されたか否かを判定する（ステップＳ２０４）。この場合、尤度Ｌｈ１４〜Ｌｈ１５は、未だ算出されていない（ステップＳ２０４でＮＯ）。よって、尤度算出部２２は、ステップＳ２０２〜Ｓ２０４の処理を繰り返す。

具体的には、尤度算出部２２は、関数ｆ４，ｆ５のなかから１つを選択する（ステップＳ２０２）。

この場合、尤度算出部２２は、関数ｆ４を選択する。次に、尤度算出部２２は、隣接距離ａｄについて、追尾物標１４０の隣接距離ａｄに該当する尤度としての尤度Ｌｈ１４を算出する（ステップＳ２０３）。具体的には、尤度算出部２２は、関数ｆ４のデータと、物標１３１についての隣接距離ａｄのデータ２１０と、を参照する。そして、尤度算出部２２は、ｎスキャン時点における物標１３１についての隣接距離ａｄを、独立変数として、関数ｆ４に代入する。これにより、関数ｆ４における従属変数が、尤度Ｌｈ１４として算出される。即ち、物標１３１における隣接距離ａｄについて、尤度Ｌｈ１４が、算出される。

次に、尤度算出部２２は、物標１３１に関する尤度Ｌｈ１１〜Ｌｈ１５の全てが算出されたか否かを判定する（ステップＳ２０４）。この場合、尤度Ｌｈ１５は、算出されていない（ステップＳ２０４でＮＯ）。よって、尤度算出部２２は、ステップＳ２０２〜Ｓ２０４の処理を繰り返す。

具体的には、尤度算出部２２は、関数ｆ５を選択する（ステップＳ２０２）。

次に、尤度算出部２２は、ドップラーシフト量ｄｓについて、追尾物標１４０のドップラーシフト量ｄｓに該当する尤度としての尤度Ｌｈ１５を算出する（ステップＳ２０３）。具体的には、尤度算出部２２は、関数ｆ５のデータと、物標１３１のドップラーシフト量ｄｓのデータ２１１と、を参照する。そして、尤度算出部２２は、ドップラーシフト量ｄｓを、独立変数として、関数ｆ５に代入する。これにより、関数ｆ５における従属変数が、尤度Ｌｈ１５として算出される。即ち、物標１３１についてのドップラーシフト量ｄｓの尤度Ｌｈ１５が、算出される。

尚、図８に示すフローチャートでは、尤度算出部２２は、各上記尤度Ｌｈ（Ｌｈ１１〜Ｌｈ１５）を、個別に算出している。しかしながら、尤度算出部２２は、行列計算により、各上記尤度Ｌｈ（Ｌｈ１１〜Ｌｈ１５）を一括して算出してもよい。

次に、尤度算出部２２は、尤度Ｌｈ１１〜Ｌｈ１５の全てが算出されたか否かを判定する（ステップＳ２０４）。この場合、尤度Ｌｈ１１〜Ｌｈ１５が算出されている（ステップＳ２０４でＹＥＳ）。よって、尤度算出部２２は、物標１３１についての総合尤度Ｌｈ＝尤度Ｌｈ１１×Ｌｈ１２×尤度Ｌｈ１３×ｈ１４×Ｌｈ１５を、算出する（ステップＳ２０５）。

次に、尤度算出部２２は、選別領域Ｓ１（ｎ）内の全ての物標１３０（１３１，１３２）について、総合尤度が算出されているか否かを判定する（ステップＳ２０６）。この場合、物標１３２についての総合尤度Ｌｈ２００は、算出されていない（ステップＳ２０６でＮＯ）。

したがって、尤度算出部２２は、ステップＳ２０１〜Ｓ２０５の処理を繰り返す。即ち、尤度算出部２２は、物標１３１についての尤度Ｌｈ１１，Ｌｈ１２，Ｌｈ１３，Ｌｈ１４，Ｌｈ１５の算出処理と同様の処理を行う。これにより、尤度算出部２２は、物標１３２についての尤度Ｌｈ２１，Ｌｈ２２，Ｌｈ２３，Ｌｈ２４，Ｌｈ２５を算出する。尚、尤度Ｌｈ２１，Ｌｈ２２，Ｌｈ２３，Ｌｈ２４，Ｌｈ２５は、それぞれ、尤度Ｌｈ１１，Ｌｈ１２，Ｌｈ１３，Ｌｈ１４，Ｌｈ１５と同様にして算出される。次いで、尤度算出部２２は、尤度Ｌｈ２１，Ｌｈ２２，Ｌｈ２３，Ｌｈ２４，Ｌｈ２５を乗算することで、総合尤度Ｌｈ２００を算出する（ステップＳ２０５）。

尤度算出部２２が、選別領域Ｓ１（ｎ）内の全ての物標１３１，１３２について、総合尤度Ｌｈ１００，Ｌｈ２００を算出した場合（ステップＳ２０６でＹＥＳ）、尤度比較部２３は、ステップＳ２０７の処理を行う。即ち、尤度比較部２３は、総合尤度Ｌｈ１００，Ｌｈ２００を用いて、追尾物標１４０を特定する（ステップＳ２０７）。本実施形態では、尤度比較部２３は、総合尤度Ｌｈ１００，Ｌｈ２００のうち、最も高い値を有する総合尤度Ｌｈ１００を選択する。

本実施形態では、図４に示されているように、物標１３１に関する物標エコー像１２１の推定位置ＰＲｅ（ｎ−１）と、物標エコー像１２１の代表点Ｐ１（ｎ）の位置とは、距離Ｄ１だけ離隔している。また、推定位置ＰＲｅ（ｎ−１）と、物標１３２に関する物標エコー像１２２の代表点Ｐ２（ｎ）の位置とは、距離Ｄ２だけ離隔している。ｎスキャン時点では、距離Ｄ１は、距離Ｄ２よりも、わずかに大きい。即ち、図６に示すように、尤度Ｌｈ１１＜Ｌｈ２１である。

一方で、図７に示すように、ｎスキャン時点において、物標１３１の面積ａｒと、面積ａｒの平均μ２とは、Δａｒ１だけ異なっている。また、ｎスキャン時点において、物標１３２の面積ａｒと、面積ａｒの平均μ２とは、Δａｒ２だけ異なっている。Δａｒ１は、Δａｒ２よりも明らかに小さい。即ち、Ｌｈ１２＞Ｌｈ２２である。

また、本実施形態では、ｎスキャン時点において、物標１３１についての他の特徴情報の尤度と、物標１３２についての対応する他の特徴情報の尤度とは、略同じである。即ち、物標１３１についての尤度Ｌｈ１３，Ｌｈ１４，Ｌｈ１５は、それぞれ、物標１３２についてのＬｈ２３，Ｌｈ２４，Ｌｈ２５と略同じである。

したがって、物標１３１についての総合尤度Ｌｈ１００は、物標１３２についての総合尤度Ｌｈ２００よりも大きい（Ｌｈ１００＞Ｌｈ２００）。この場合、尤度比較部２３は、物標１３１を、ｎスキャン時点における追尾物標１４０として特定する。

図１に示すように、尤度比較部２３が追尾物標１４０を特定した後、面積判定部１５は、追尾物標１４０についての物標エコー像１２１の面積ａｒを、判定する。これにより、面積判定部１５は、物標エコー像１２１について、他の物標エコー像１２０との融合、又は、過大なぼやけが生じているか否か、を判定する。

図９は、面積判定部１５での処理の流れの一例を説明するためのフローチャートである。面積判定部１５は、まず、関数ｆ２のデータを、関数更新部２１から読み出す（ステップＳ３０１）。面積判定部１５は、更に、追尾物標１４０の物標エコー像１２１について、ｎスキャン時点での面積ａｒのデータ２０７を、特徴情報メモリ１２から読み出す（ステップＳ３０２）。

次に、面積判定部１５は、ｎスキャン時点における、物標エコー像１２１の面積ａｒが、所定のしきい値より大きいか否かを判定する。この場合のしきい値ｔｈは、例えば、（μ２＋３×σ２）であり、平均μ２に基づいて設定されている。σ２は、関数ｆ２における標準偏差であり、分散σ２^２の正の平方根である。物標エコー像１２１の面積ａｒが、上記所定のしきい値ｔｈよりも大きい場合（ステップＳ３０３でＹＥＳ）、面積判定部１５は、ｎスキャン時点での物標エコー像１２１に、他のエコー像１２０との融合、又は過大なぼやけが生じていると判定する（ステップＳ３０４）。この場合、面積判定部１５は、関連付け部１４で代表点ＰＲ（ｎ）を特定できなかったものとして扱い、代表点ＰＲ（ｎ）のデータを、運動推定部１６へ出力しないで（ステップＳ３０５）、処理を完了する。

一方、物標エコー像１２１の面積ａｒが、所定のしきい値ｔｈ以下である場合（ステップＳ３０３でＮＯ）、面積判定部１５は、ｎスキャン時点での物標エコー像１２１に、他の物標エコー像１２０との融合、及び、過大なぼやけの何れも生じていないと判定する（ステップＳ３０６）。この場合、面積判定部１５は、代表点ＰＲ（ｎ）のデータを、運動推定部１６へ出力し（ステップＳ３０７）、処理を完了する。

上記した、（ａ）物標エコー像１２１と、他の物標エコー像１２０との融合、及び（ｂ）物標エコー像１２１のぼやけについて、以下に説明する。

まず、（ａ）物標エコー像１２１と、他の物標エコー像１２０（１２５）との融合が生じる理由を、図１０及び図１１等を用いて説明する。図１０は、物標エコー像１２１と、他の物標エコー像１２５とについて、融合が生じていない場合を説明するための模式図である。図１１は、物標エコー像１２１と、他の物標エコー像１２５とについて、融合が生じている場合を説明するための模式図である。尚、物標エコー像１２５によって特定される物標１３５として、他船、クラッタ等を例示することができる。

図１、図１０及び図１１を参照しながら説明すると、アンテナユニット２で検出された各物標エコー像１２０（１２１，１２５）の外郭形状は、アンテナユニット２における観測誤差に起因して、ぼんやりした形状となる。即ち、物標エコー像１２０（１２１，１２５）の外郭形状の大きさは、物標１３０（１３１，１３５）の外郭形状の大きさよりも大きくなる。しかしながら、図１０に示すように、物標１３１と、物標１３５とが十分に離隔している場合、物標エコー像１２１と、物標エコー像１２５とは、融合しない。即ち、追尾処理部１１は、物標エコー像１２１と、物標エコー像１２５とを、別個の物標であると判定する。

一方、物標１３１と、物標１３５との距離が小さい場合、図１１に示すように、物標エコー像１２１の外郭部と、物標エコー像１２５の外郭部とは、融合してしまう。即ち、物標エコー像１２１と、物標エコー像１２５とが融合してしまう。この場合、追尾処理部１１は、物標エコー像１２１と、物標エコー像１２５とを、１つの物標エコー像であると、誤って判定する。その結果、エコー検出部１０は、上記１つの物標エコー像の中心点ＰＦ１を、物標エコー像１２１の代表点として検出してしまう。即ち、エコー検出部１０は、中心点ＰＦ１を、追尾物標１４０の代表点ＰＲ（ｎ）として検出してしまう。この場合、運動推定部１６は、この代表点ＰＦ１を用いて追尾フィルタ処理を行うこととなる。その結果、運動推定部１６で算出される推定速度ベクトルＶ１（ｎ）は、実際の速度ベクトルに対する誤差が大きくなってしまう。

物標エコー像１２１，１２５の融合が、ｎスキャン時点で生じた場合、ｎスキャン時点での物標エコー像１２１の面積ａｒは、極端に大きくなる。即ち、追尾物標１４０についての面積ａｒが、所定のしきい値ｔｈを超える。よって、面積判定部１５は、この追尾物標の面積ａｒを判定することにより、物標エコー像１２１と、他の物標エコー像１２５とが融合しているか否かを判定できる。

次に、上記（ｂ）物標エコー像１２１に過大なぼやけが生じる理由を、図１２等を用いて説明する。図１２は、物標エコー像１２１に過大なぼやけが生じている状態を示す模式図である。図１及び図１２を参照しながら説明すると、アンテナユニット２における、アンテナ角度θ方向の分解能等に起因して、物標エコー像１２１は、特に、アンテナ角度θ方向に大きくぼやけた形状として検出される場合がある。その結果、ｎスキャン時点での物標エコー像１２１の面積ａｒが、過度に大きくなる場合がある。この場合、物標エコー像１２１の面積ａｒが、所定のしきい値ｔｈを超える。また、この場合、物標エコー像１２１における代表点ＰＦ２の位置は、物標エコー像１２１における真の代表点Ｐ１（ｎ）の位置から大きくずれてしまう。この場合に、運動推定部１６が、代表点ＰＦ２を用いて追尾フィルタ処理を行うと、運動推定部１６で算出される推定速度ベクトルＶ１（ｎ）は、実際の速度ベクトルに対する誤差が大きくなる。そこで、面積判定部１５は、追尾物標１４０（物標エコー像１２１）の面積ａｒを判定することにより、物標エコー像１２１の面積ａｒが過度に大きくなっているか否かを検出する。そして、面積ａｒがしきい値ｔｈを超えている場合、面積判定部１５は、エコー像１２１に過大なぼやけが生じていると判定する。

次に、運動推定部１６について説明する。運動推定部１６は、追尾フィルタ処理を行うように構成されている。これにより、運動推定部１６は、ｎスキャン時点における、追尾物標１４０の推定速度ベクトルＶ１（ｎ）を算出する。尚、追尾フィルタとして、α−βフィルタ、カルマン・フィルタ等を例示することができる。運動推定部１６は、推定速度ベクトルＶ１（ｎ）を特定するデータを、選別部１３、及び表示器４へ出力する。推定速度ベクトルＶ１（ｎ）のデータは、（ｎ＋１）スキャン時点での追尾フィルタ処理に用いられる。

次に、運動推定部１６での処理の一例を説明する。図１３は、運動推定部１６における処理の流れの一例を説明するためのフローチャートである。図１４は、運動推定部１６での処理を説明するための模式図である。図１４は、追尾物標１４０を特定する物標エコー像１２１について、他の物標エコー像１２５との融合、及び、過大なぼやけの何れも生じていない場合を示している。図１５は、運動推定部１６での処理を説明するための模式図である。図１５は、追尾物標１４０を特定する物標エコー像１２１について、他の物標エコー像１２５との融合、又は過大なぼやけが生じている場合を示している。

まず、運動推定部１６は、（ｎ−１）スキャン時点における、当該運動推定部１６での追尾フィルタ処理の結果を特定するデータを、読み出す（ステップＳ４０１）。この場合、上記追尾フィルタ処理の結果は、図１４及び図１５にそれぞれ示すように、（ｎ−１）スキャン時点における追尾物標１４０についての、平滑位置ＰＲｓ（ｎ−１）の座標と、推定位置ＰＲｅ（ｎ−１）の座標と、推定速度ベクトルＶ１（ｎ−１）と、を含んでいる。

次に、運動推定部１６は、面積判定部１５から代表点ＰＲ（ｎ）のデータが出力されているか否かを判定する（ステップＳ４０２）。代表点ＰＲ（ｎ）のデータが出力されている場合（ステップＳ４０２でＹＥＳ）、即ち、物標エコー像１２１において、他の物標エコー像１２５との融合、及び、過大なぼやけの何れも生じていない場合、運動推定部１６は、ｎスキャン時点で観測された代表点ＰＲ（ｎ）のデータを用いて、追尾処理を行う（ステップＳ４０３〜Ｓ４０４）。

具体的には、運動推定部１６は、追尾物標１４０の代表点ＰＲ（ｎ）の座標データを、尤度算出部２２から読み出す（ステップＳ４０３）。その後、運動推定部１６は、ステップＳ４０１，Ｓ４０３で読み出したデータを用いて、追尾フィルタ処理を行う（ステップＳ４０４）。これにより、図１４に示すように、運動推定部１６は、ｎスキャン時点における、物標エコー像１２１についての平滑位置ＰＲｓ（ｎ）と、推定位置ＰＲｅ（ｎ）と、推定速度ベクトルＶ１（ｎ）と、を算出する。

運動推定部１６は、上記追尾フィルタ処理の結果を特定するデータを、特徴情報メモリ１２と、選別部１３と、表示器４と、に出力する（ステップＳ４０５）。

一方、代表点ＰＲ（ｎ）のデータが出力されていない場合（ステップＳ４０２でＮＯ）、即ち、物標エコー像１２１に、他の物標エコー像１２０との融合、又は過大なぼやけが生じている場合（ステップＳ４０２でＮＯ）、運動推定部１６は、ｎスキャン時点での代表点ＰＲ（ｎ）を用いることなく、追尾処理を行う（ステップＳ４０４）。

この場合、図１５に示すように、運動推定部１６で算出された推定速度ベクトルＶ１（ｎ）は、推定速度ベクトルＶ１（ｎ−１）と同じ向き及び大きさである。また、推定速度ベクトルＶ１（ｎ）の始点は、推定速度ベクトルＶ１（ｎ−１）の終点である。その後、運動推定部１６は、ステップＳ４０５の処理を行う。以上が、運動推定部１６における処理の流れの一例である。

表示器４は、例えばカラー表示可能な液晶ディスプレイである。表示器４は、各物標エコー像１２０の画像データを用いて、表示画面に、各物標エコー像１２０を表示する。また、表示器４は、推定速度ベクトルＶ１（ｎ）を、画像として表示する。これにより、表示器４の表示画面には、追尾物標１４０（物標エコー像１２１）について、推定速度ベクトルＶ１（ｎ）を示す画像が表示される。レーダ装置１のオペレータは、表示器４に表示されたレーダ映像を確認することにより、追尾物標１４０の運動状態を確認することができる。

以上説明したように、追尾処理装置３によると、関連付け部１４の尤度算出部２２は、選別領域Ｓ１（ｎ）の代表点Ｐ１（ｎ），Ｐ２（ｎ）について、追尾物標１４０の代表点ＰＲ（ｎ）に該当する尤度Ｌｈ１１，Ｌｈ２１を算出する。そして、関連付け部１４の尤度比較部２３は、算出された尤度Ｌｈ１１，Ｌｈ２１に基づいて、複数の物標１３１，１３２のなかから、追尾物標１４０を特定する。このように、関連付け部１４は、統計的手法を用いることで、追尾物標１４０の代表点ＰＲ（ｎ）、即ち代表点Ｐ１（ｎ）を、より正確に判定することができる。また、関連付け部１４は、尤度Ｌｈ１１，２１の算出に際して、各物標１３０についての予測誤差、及び観測誤差等を考慮したモデルを設定する必要が無い。即ち、尤度算出部２２は、複雑なモデルに従って演算を行う必要が無い。このため、関連付け部１４における演算処理は、簡易である。したがって、追尾処理装置３によると、複数の物標１３０（１３１，１３２）のなかから、追尾物標１４０を確実に、且つ、簡易に検出することができる。

また、追尾処理装置３によると、関連付け部１４は、選別領域Ｓ１（ｎ）内に位置する各物標１３１，１３２について、代表点Ｐ１（ｎ），Ｐ２（ｎ）の情報以外の特徴情報を参照する。そして、関連付け部１４の尤度算出部２２は、追尾物標１４０の特徴情報に該当する尤度Ｌｈ１２〜Ｌｈ１５；Ｌｈ２２〜Ｌｈ２５を算出する。また、関連付け部１４の尤度比較部２３は、算出した尤度Ｌｈ１１〜Ｌｈ１５；Ｌｈ２１〜Ｌｈ２５に基づいて、各物標１３１，１３２のなかから、追尾物標１４０を特定する。このように、関連付け部１４は、各物標１３１，１３２の代表点Ｐ１（ｎ），Ｐ２（ｎ）についての尤度Ｌｈ１１，２１に加え、尤度Ｌｈ１２〜Ｌｈ１５；Ｌｈ２２〜Ｌｈ２５を用いて、追尾物標１４０を特定できる。したがって、複数の物標１３１，１３２のなかから、追尾物標１４０を、より確実に判別することができる。また、前述したように、尤度Ｌｈを用いた処理は、複雑なモデルに従う複雑な演算処理と異なり、簡易な演算処理で済む。したがって、代表点Ｐ（ｎ）以外の追加の特徴情報について尤度判定を行う場合でも、関連付け部１４での演算負荷は、少なくて済む。

また、例えば、船舶の航行によって生じる波、又はクラッタ等を、エコー検出部１０が、物標として検出した場合を考える。この場合、上記の波又はクラッタ等の面積についての尤度は、物標１３１，１３２についての対応する尤度Ｌｈ１２，２２と比べて、格段に小さい。よって、尤度比較部２３は、上記波又はクラッタ等を、追尾物標１４０として誤って判定することを抑制できる。

また、例えば、前述の特許文献１，２には、追尾対象における１つの特徴情報としての航跡情報に基づく、多重仮説による追尾処理が開示されている。これに対し、追尾処理装置３は、複数の特徴情報を用いて、複数の物標１３０のなかから追尾物標１４０を特定している。これにより、複数の物標１３０のなかから、追尾物標１４０を、より確実に判別することができる。

また、追尾処理装置３によると、関連付け部１４は、各エコー像１２０に関して、面積ａｒの情報、エコー信号のレベルｅｃ、隣接距離ａｄの情報、及び、ドップラーシフト量ｄｓの少なくとも一つを用いる。関連付け部１４は、これらの特徴情報を用いることにより、複数の物標１３０のなかから、追尾物標１４０を、より正確に特定することができる。より具体的には、特徴情報として、面積ａｒの情報が用いられている。これにより、関連付け部１４は、追尾物標１４０（物標１３１）の形状と異なる物標１３２を、誤って追尾物標１４０であると判別することを抑制できる。例えば、追尾物標が大型船舶である場合には、小型船舶からなる物標、又はクラッタに起因して検出される物標を、関連付け部１４が追尾対象として検出することを、より確実に抑制できる。また、特徴情報として、エコーレベルｅｃが用いられている。これにより、関連付け部１４は、追尾物標１４０とはエコーレベルｅｃの異なる物標１３０を、誤って追尾物標１４０であると判別することを抑制できる。また、特徴情報として、隣接距離ａｄの情報が用いられている。これにより、関連付け部１４は、追尾物標１４０（物標１３１）と、クラッタ等の物標とを、より確実に判別することができる。また、特徴情報として、ドップラーシフト量ｄｓが用いられている。これにより、エコー検出部１０は、自船１００と、物標１３０とが向かい合う方向における、自船１００と物標１３０との相対速度を検出することができる。したがって、関連付け部１４は、この相対速度を考慮して、複数の物標１３０のなかから、追尾物標１４０を特定できる。その結果、例えば、追尾物標１４０（物標１３１）と、他の物標１３０とがすれ違う際において、尤度比較部２３は、当該他の物標１３０を、追尾物標１４０であると判別することを防止できる。

また、追尾処理装置３によると、関連付け部１４の尤度算出部２２は、確率密度関数としての関数ｆ（ｆ１〜ｆ５）を用いて、尤度Ｌｈ１１〜Ｌｈ１５；Ｌｈ２１〜Ｌｈ２５を算出する。このように、尤度算出部２２は、関数ｆを用いることにより、尤度Ｌｈ１１〜Ｌｈ１５；Ｌｈ２１〜Ｌｈ２５を、統計的に、且つ、簡易な演算で算出できる。

また、追尾処理装置３によると、関連付け部１４の尤度算出部２２は、各関数ｆ（ｆ１〜ｆ５）において、エコー検出部１０で検出された特徴情報（代表点Ｐ１（ｎ）及びＰ２（ｎ）の座標、面積ａｒ、エコーレベルｅｃ、隣接距離ａｄ、ドップラーシフト量ｄｓ）を、独立変数として用いる。また、尤度算出部２２は、上記独立変数に従属する従属変数を、尤度Ｌｈ１１〜Ｌｈ１５；Ｌｈ２１〜Ｌｈ２５として算出する。その結果、代表点Ｐ２（ｎ）の位置が、推定位置ＰＲｅ（ｎ−１）に近いという理由のみに基づいて、当該代表点Ｐ２（ｎ）を、追尾物標１４０の代表点ＰＲ（ｎ）として特定することは、防止される。即ち、関連付け部１４は、代表点Ｐ１（ｎ），Ｐ２（ｎ）の情報以外の、他の特徴情報における尤度Ｌｈ１２〜Ｌｈ１５；Ｌｈ２２〜Ｌｈ２５を考慮して、追尾物標１４０を特定する。これにより、関連付け部１４は、各物標１３１，１３２のなかから、追尾物標１４０を、より正確に特定できる。

また、追尾処理装置３によると、関連付け部１４の尤度算出部２２は、各物標１３１，１３２毎に、総合尤度Ｌｈ１００，Ｌｈ２００を算出する。そして、尤度比較部２３は、各物標１３１，１３２のうち、総合尤度Ｌｈ１００，Ｌｈ２００が最も高い物標１３１を、追尾物標１４０として特定する。即ち、尤度比較部２３は、各物標１３１，１３２について、複数の特徴情報に基づく尤度Ｌｈ１１〜Ｌｈ１５；Ｌｈ２１〜Ｌｈ２５を、総合的に判定する。これにより、関連付け部１４は、各物標１３１，１３２のなかから、追尾物標１４０を、より確実に特定できる。

また、追尾処理装置３によると、面積判定部１５は、追尾物標１４０を特定する物標エコー像１２１における面積ａｒが、所定のしきい値ｔｈを超えているか否かを判定する。また、運動推定部１６は、面積ａｒが、所定のしきい値ｔｈ以下である場合、ｎスキャン時点での代表点ＰＲ（ｎ）の情報を用いて、推定速度ベクトルＶ１（ｎ）を算出する。また、運動推定部１６は、面積ａｒが、所定のしきい値ｔｈを超えている場合には、ｎスキャン時点での代表点ＰＲ（ｎ）の情報を用いることなく、推定速度ベクトルＶ１（ｎ）を算出する。このような構成により、運動推定部１６は、代表点ＰＲ（ｎ）について、正確な運動推定に適した情報が得られている場合には、当該情報を用いて、追尾物標１４０の運動を、より正確に算出できる。一方、運動推定部１６は、代表点ＰＲ（ｎ）について、正確な運動推定に適した情報を有していない場合も、当該情報用いることなく、追尾物標１４０の運動を、より正確に算出できる。

また、面積判定部１５は、追尾物標１４０の物標エコー像１２１の面積ａｒについて、関数ｆ２を用いて、判定を行う。このため、アンテナユニット２での観測誤差等に起因して、追尾物標１４０の物標エコー像１２１の面積ａｒが経時的に変化する場合でも、面積判定部１５は、当該面積ａｒを、より正確に判定できる。

また、特許文献３には、追尾対象の観測位置に飛躍（飛び）が生じてしまっている状態への対応を行うための構成が開示されている。しかしながら、特許文献３に記載の構成では、上記した、エコー像同士の融合を検出する精度は、十分とはいえない。これに対し、追尾処理装置３によると、面積判定部１５は、関数ｆ２を用いて、面積ａｒを判定する。その結果、面積判定部１５は、追尾物標１４０を特定する物標エコー像１２１が、他のエコー像１２５と融合した場合に、当該融合を、より確実に検出できる。

また、追尾処理装置３によると、面積判定部１５は、関数ｆ２での平均μ２に基づいて、しきい値ｔｈを設定している。そして、面積判定部１５は、しきい値ｔｈを基準として、追尾物標１４０についての面積ａｒの大きさを判定している。このような構成によると、上記面積ａｒの平均μから、大きくずれている値が面積ａｒとして検出されている場合、面積判定部１５は、追尾物標１４０を特定する物標エコー像１２１について、他の物標エコー像１２５との融合、又は、過大なぼやけが生じていると判定する。このような構成により、面積判定部１５は、アンテナユニット２での観測誤差に左右され難い状態で、物標エコー像１２１の面積ａｒを、より正確に判定できる。

以上、本発明の実施形態について説明したけれども、本発明は上述の実施の形態に限られるものではなく、特許請求の範囲に記載した限りにおいて様々な変更が可能である。例えば、次のように変更して実施してもよい。

（１）上述の実施形態では、物標についての代表点情報以外の追加の特徴情報として、物標エコー像についての面積の情報等を用いる形態を例に説明した。しかしながら、この通りでなくてもよい。例えば、上記追加の特徴情報として、エコー検出部で検出された物標の航跡情報を用いてもよい。上記追加の特徴情報として、航跡情報を用いた場合、錯綜するように航行する複数の物標のなかから、追尾物標を、より精度よく特定できる。

（２）上述の実施形態では、形状情報として、物標エコー像の面積の情報を用いる形態を例に説明した。しかしながら、この通りでなくてもよい。例えば、形状情報として、物標エコー像の幅を用いてもよい。また、形状情報として、物標エコー像の長さを用いてもよい。

（３）また、上述の実施形態では、物標エコー像についての中心点を、代表点として用いる形態を例に説明した。しかしながら、この通りでなくてもよい。例えば、物標エコー像のうち、自船に対して最も近い点である最前縁点等を、代表点として用いてもよい。

（４）また、上述の実施形態では、追尾処理装置が、船舶用の追尾処理装置である形態を例に説明した。しかしながら、本発明は、船舶用の追尾処理装置に限らず、他の物標用の追尾処理装置として適用することができる。

本発明は、追尾処理装置、及び追尾処理方法として、広く適用することができる。

３ 追尾処理装置
１０ エコー検出部（検出部）
９ 信号処理部（検出部）
１９ 関連付け部（尤度判定部）
１３０ 物標
１４０ 追尾物標（追尾対象）
Ｌｈ１１，Ｌｈ２１ 追尾対象の代表点に該当する尤度
Ｐ 代表点

Claims (7)

1. 追尾対象として選択された物標を追尾するための追尾処理装置であって、
   複数の物標について、代表点の情報を含む特徴情報を検出するための、検出部と、
   複数の前記代表点の少なくとも一部について、それぞれ、前記追尾対象の代表点に該当する尤度を算出する、尤度判定部と、
   前記追尾対象として特定された前記物標の前記特徴情報を用いて、前記追尾対象の運動状態を推定する、運動推定部と、
   面積を判定するための面積判定部と、を備え、
   前記尤度判定部は、算出された前記尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定し、
   前記検出部は、複数の前記物標の少なくとも一部について、前記代表点の情報以外の特徴情報を、追加の特徴情報として検出し、
   前記尤度判定部は、追加の前記特徴情報について、前記追尾対象の特徴情報に該当する尤度を算出し、且つ、当該尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定し、
   前記検出部は、送信信号に対する前記追尾対象を含む複数の前記物標での反射波であるエコー信号を検出し、且つ、当該エコー信号によって特定されるエコー像を検出し、且つ、当該エコー像の面積を追加の前記特徴情報として検出し、
   前記面積判定部は、前記追尾対象の前記エコー像の前記面積が所定のしきい値を超えているか否かを判定し、
   前記運動推定部は、前記面積が所定のしきい値以下である場合、前記運動状態を推定する時点での前記追尾対象の前記特徴情報を用いて、前記運動状態を算出し、且つ、前記面積が所定のしきい値を超えている場合、前記運動状態を算出する時点での前記追尾対象の前記特徴情報を用いることなく、前記運動状態を算出することを特徴とする、追尾処理装置。
2. 請求項１に記載の追尾処理装置であって、
   前記検出部は、送信信号に対する各前記物標での反射波であるエコー信号を用いて、各前記特徴情報を検出し、
   追加の前記特徴情報は、前記エコー信号によって特定されるエコー像の形状情報、前記エコー信号のレベル、各前記物標の周囲の状態を特定する情報、及び、前記エコー信号に関するドップラーシフト量の少なくとも一つを含むことを特徴とする、追尾処理装置。
3. 請求項１又は請求項２に記載の追尾処理装置であって、
   各前記特徴情報は、数値によって表される情報であり、
   前記尤度判定部は、各前記特徴情報毎に、確率密度関数を用いて、各前記尤度を算出することを特徴とする、追尾処理装置。
4. 請求項３に記載の追尾処理装置であって、
   前記確率密度関数は、各前記特徴情報毎に設定され、
   各前記確率密度関数において平均の値として用いられる値は、対応する特徴情報について、尤度が算出される前に予め推定された値であり、
   前記尤度判定部は、各確率密度関数において、対応する前記特徴情報を独立変数として用い、前記独立変数に従属する従属変数を、前記尤度として算出することを特徴とする、追尾処理装置。
5. 請求項４に記載の追尾処理装置であって、
   前記尤度判定部は、
   同一物標についての前記尤度を互いに乗算することで、総合尤度を算出し、且つ、
   各前記物標のうち、前記総合尤度が最も高い前記物標を、前記追尾対象として特定することを特徴とする、追尾処理装置。
6. 請求項１乃至請求項５の何れか１項に記載の追尾処理装置であって、
   前記面積判定部は、確率密度関数を用いて前記面積を判定するように構成されており、
   前記確率密度変数の平均の値として用いられる値は、前記面積について、前記面積判定部で判定される前に予め推定された値であり、
   前記所定のしきい値は、前記平均に基づいて設定されることを特徴とする、追尾処理装置。
7. 追尾対象として選択された物標を追尾するための追尾処理方法であって、
   複数の物標について、代表点の情報を含む特徴情報を検出するための、検出ステップと、
   複数の前記代表点の少なくとも一部について、それぞれ、前記追尾対象の代表点に該当する尤度を算出する、尤度判定ステップと、
   前記追尾対象として特定された前記物標の前記特徴情報を用いて、前記追尾対象の運動状態を推定する、運動推定ステップと、
   面積を判定するための面積判定ステップと、を含み、
   前記尤度判定ステップでは、算出された前記尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定し、
   前記検出ステップでは、複数の前記物標の少なくとも一部について、前記代表点の情報以外の特徴情報を、追加の特徴情報として検出し、
   前記尤度判定ステップでは、追加の前記特徴情報について、前記追尾対象の特徴情報に該当する尤度を算出し、且つ、当該尤度に基づいて、複数の前記物標のなかから前記追尾対象を特定し、
   前記検出ステップでは、送信信号に対する前記追尾対象を含む複数の前記物標での反射波であるエコー信号を検出し、且つ、当該エコー信号によって特定されるエコー像を検出し、且つ、当該エコー像の面積を追加の前記特徴情報として検出し、
   前記面積判定ステップでは、前記追尾対象の前記エコー像の前記面積が所定のしきい値を超えているか否かを判定し、
   前記運動推定部ステップでは、前記面積が所定のしきい値以下である場合、前記運動状態を推定する時点での前記追尾対象の前記特徴情報を用いて、前記運動状態を算出し、且つ、前記面積が所定のしきい値を超えている場合、前記運動状態を算出する時点での前記追尾対象の前記特徴情報を用いることなく、前記運動状態を算出することを特徴とする、追尾処理方法。

Images