(57)【特許請求の範囲】
【請求項１】
時系列的に取得された目標信号から得られる複数の目標位置に関するデータに基づいて目標の追尾を行う目標追尾装置において、
あるタイミングで取得された複数の目標位置と、過去のタイミングで取得された複数の目標位置または過去のタイミングで取得された目標位置に基づき求められたあるタイミングにおける複数の目標予測位置である複数の被マッチング目標位置との対応付けを行うマッチング処理手段と、
マッチング処理手段で被マッチング目標位置と対応付けられた目標位置を、被マッチング目標位置に対応する目標識別データと関連付けて格納するマッチングデータテーブルと、
過去のタイミングで取得された目標位置または過去のタイミングで取得された目標位置及び過去の２つ以上のタイミングで取得された目標位置に基づき求めた目標速度から目標の予測位置を求めて、該予測位置を中心としたウィンドウを追尾対象の目標毎に計算するウィンドウ計算手段と、
ウィンドウ計算手段で計算されたウィンドウ内にあるデータから目標を抽出する目標抽出手段と、を備え、
前記目標抽出手段は、ウィンドウ内に目標とするべきデータが２つ以上存在するとき、またはウィンドウ内に目標とするべきデータが存在しないときに、前記マッチングデータテーブルを参照して、追尾対象の目標の目標識別データに関連付けられた目標位置を抽出することを特徴とする目標追尾装置。
【請求項２】
過去の２つ以上のタイミングで取得され同じ目標識別データと関連付けられた目標位置に基づき該目標識別データに対応する目標速度を求める速度情報算出手段をさらに備え、
前記マッチング処理手段は、前記目標速度に基づきあるタイミングにおける目標予測位置を求めて、該目標予測位置とあるタイミングで取得された前記目標位置との対応付けを行うことを特徴とする請求項１記載の目標追尾装置。
【請求項３】
前記マッチング処理手段は、複数の目標位置と複数の被マッチング目標位置に対してその距離の総和が最小となるような対応付けを行うことを特徴とする請求項１または２記載の目標追尾装置。
【請求項４】
前記マッチング処理手段は、複数の目標位置と複数の被マッチング目標位置に対してその距離及び目標の大きさの差異の乗算の総和が最小となるような対応付けを行うことを特徴とする請求項１または２記載の目標追尾装置。
【請求項５】
前記マッチング処理手段は、目標位置と被マッチング目標位置との距離が所定距離を超えている場合には、該目標位置と該被マッチング目標位置との対応付けは行わないことを特徴とする請求項３または４記載の目標追尾装置。
【請求項６】
前記所定距離は、前記目標信号より得られる目標の大きさに応じて変化させることを特徴とする請求項５記載の目標追尾装置。
【請求項７】
前記所定距離は、目標速度の大きさに応じて変化させることを特徴とする請求項５記載の目標追尾装置。
【請求項８】
前記マッチング処理手段は、
複数の目標位置と複数の被マッチング目標位置とのいずれか一方の位置（以下、第１目
標位置）に対して、その一方の位置から最も近い他方の位置（以下、第２目標位置）を仮対応付ける仮対応付け手段と、
複数の第１目標位置が仮対応付けられた第２目標位置に対して、一つの第１目標位置を仮対応付けて、別の第１目標位置を別の第２目標位置に仮対応付ける仮対応付け選択手段と、
全ての第２目標位置に１以下の第１目標位置が仮対応付けられたときに、仮対応付けられた第１目標位置と第２目標位置との対応付けを確定する確定手段と、
を備えることを特徴とする請求項１ないし７のいずれか１項に記載の目標追尾装置。
【請求項９】
前記仮対応付け選択手段は、
同じ第２目標位置に対して仮対応付けられた２つの第１目標位置Ｘ、Ｙに対して、それぞれ仮対応付けられた第２目標位置をＸ１、Ｙ１とし、Ｘ１、Ｙ１の次に近い第２目標位置をＸ２、Ｙ２としたときに、Ｘｄ＝｜Ｘ−Ｘ１｜＋｜Ｙ−Ｙ２｜とＹｄ＝｜Ｘ−Ｘ２｜＋｜Ｙ−Ｙ１｜（ここで、｜Ａ−Ｂ｜はＡ、Ｂ間の距離を表す）の大小を比較する手段と、
Ｘｄ≦Ｙｄである場合には、Ｘの仮対応点としてＸ１を選択し、Ｙの仮対応点としてＹ２を選択し、Ｘｄ＞Ｙｄである場合には、Ｙの仮対応点としてＹ１を選択し、Ｘの仮対応点としてＸ２を選択する手段と、
を備えることを特徴とする請求項８記載の目標追尾装置。
【請求項１０】
時系列的に取得された目標信号から得られる複数の目標位置に関するデータに基づいて目標の追尾を行う目標追尾方法において、
あるタイミングで取得された複数の目標位置と、過去のタイミングで取得された複数の目標位置または過去のタイミングで取得された目標位置に基づき求められたあるタイミングにおける複数の目標予測位置である複数の被マッチング目標位置との対応付けを行うマッチング処理工程と、
マッチング処理工程で被マッチング目標位置と対応付けられた目標位置を、被マッチング目標位置に対応する目標識別データと関連付けて格納する工程と、
過去のタイミングで取得された目標位置または過去のタイミングで取得された目標位置及び過去の２つ以上のタイミングで取得された目標位置に基づき求めた目標速度から目標の予測位置を求めて、該予測位置を中心としたウィンドウを追尾対象の目標毎に計算するウィンドウ計算工程と、
ウィンドウ計算工程で計算されたウィンドウ内にあるデータから目標を抽出する目標抽出工程と、を備え、
前記目標抽出工程は、ウィンドウ内に目標とするべきデータが２つ以上存在するとき、またはウィンドウ内に目標とするべきデータが存在しないときに、追尾対象の目標の目標識別データに関連付けて格納された目標位置を抽出することを特徴とする目標追尾方法。
【請求項１１】
前記マッチング処理工程は、
（ａ）  対象となる複数の目標位置と対象となる複数の被マッチング目標位置とのいずれか一方の位置（以下、第１目標位置）に対して、その一方の位置から最も近い他方の位置（以下、第２目標位置）を仮対応付ける工程と、
（ｂ）  第２目標位置に仮対応付けられた第１目標位置の数を求める工程と、
（ｃ）  同じ第２目標位置に仮対応付けられた第１目標位置の数が２以上である場合には、その数が２つである場合に、その２つの第１目標位置Ｘ、Ｙを選択し、その数が３以上である場合に、任意の２つの第１目標位置Ｘ、Ｙを選択する工程と、
（ｄ）  ２つの第１目標位置Ｘ、Ｙに対して、それぞれ仮対応付けられた第２目標位置をＸ１、Ｙ１とし、Ｘ１、Ｙ１の次に近い第２目標位置をＸ２、Ｙ２としたときに、Ｘｄ＝｜Ｘ−Ｘ１｜＋｜Ｙ−Ｙ２｜とＹｄ＝｜Ｘ−Ｘ２｜＋｜Ｙ−Ｙ１｜（ここで、｜Ａ−Ｂ｜はＡ、Ｂ間の距離を表す）の大小を比較する工程と、
（ｅ）  Ｘｄ≦Ｙｄである場合には、Ｘの仮対応点としてＸ１を選択し、Ｙの仮対応点としてＹ２を選択し、Ｘｄ＞Ｙｄである場合には、Ｙの仮対応点としてＹ１を選択し、Ｘの仮対応点としてＸ２を選択する工程と、
（ｆ）  （ｂ）、（ｄ）、（ｅ）の工程を、全ての第２目標位置に１以下の第１目標位置が仮対応付けられるまで繰り返す工程と、
（ｇ）  全ての第２目標位置に仮対応付けられた第１目標位置の数が１以下となったときに、仮対応付けられた第１目標位置と第２目標位置との対応付けを確定する工程と、
を備えることを特徴とする請求項１０記載の目標追尾方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、船舶等の目標の追尾を行う目標追尾装置及び目標追尾方法に関する。
【背景技術】
【０００２】
例えば、船舶等の衝突予防を行うＡＲＰＡ装置（Automatic Radar Plotting Aids）では、目標の追尾を行う目標追尾装置が設けられる。目標追尾装置では、レーダエコーの或るスキャンにおいて追尾している目標に対して、次のスキャンにおいて、その目標の周辺または過去の位置から算出したその目標の予測位置の周辺にウィンドウ（またはゲートともいう）を作成し、ウィンドウの中の目標の重心位置を求めて、それを目標の新しい位置として、順次、これを繰り返すことで、目標の追尾を行っている（例えば、特許文献１）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第３１１０５４９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来のウィンドウを設定する追尾においては次のような問題がある。
１）１つのウィンドウ内で１つの目標の追尾を行うことが原則となっているために、１つのウィンドウ内に複数のエコーが発生した場合には、そのウィンドウの中から目標のエコーを選択する必要がある。そのため、間違ったエコーを選択すると、乗り移りが発生する。特に、目標が輻輳している地帯では、１つのウィンドウ内に複数のエコーが発生することが頻繁に起こり、誤った選択がなされる可能性が高くなる。
２）追尾中の目標が急速な加減速や旋回を行った場合、ウィンドウからエコーが外れてしまう場合がある。一度ウィンドウからエコーが外れてしまうと、追尾を継続するのは非常に困難であり、目標を見失う（ロストする）結果となる。
【０００５】
ウィンドウの大きさが小さいとロストの可能性があり、逆に大きければ乗り移りの可能性があり、上記１）と２）の問題はトレードオフの関係にある。特に、近年船舶等の目標の高速化により、ウィンドウを大きくとる必要があり、そのため、乗り移りの可能性が高まっている。
【０００６】
本発明はかかる課題に鑑みなされたもので、その目的は、従来のウィンドウを用いた追尾の問題を解消することができる目標追尾装置及び目標追尾方法を提供することである。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明の請求項１記載の発明は、時系列的に取得された目標信号から得られる複数の目標位置に関するデータに基づいて目標の追尾を行う目標追尾装置において、
あるタイミングで取得された複数の目標位置と、過去のタイミングで取得された複数の目標位置または過去のタイミングで取得された目標位置に基づき求められたあるタイミングにおける複数の目標予測位置である複数の被マッチング目標位置との対応付けを行うマッチング処理手段と、
マッチング処理手段で被マッチング目標位置と対応付けられた目標位置を、被マッチング目標位置に対応する目標識別データと関連付けて格納するマッチングデータテーブルと、
過去のタイミングで取得された目標位置または過去のタイミングで取得された目標位置及び過去の２つ以上のタイミングで取得された目標位置に基づき求めた目標速度から目標の予測位置を求めて、該予測位置を中心としたウィンドウを追尾対象の目標毎に計算するウィンドウ計算工程と、
ウィンドウ計算工程で計算されたウィンドウ内にあるデータから目標を抽出する目標抽出工程と、を備え、
前記目標抽出工程は、ウィンドウ内に目標とするべきデータが２つ以上存在するとき、またはウィンドウ内に目標とするべきデータが存在しないときに、前記マッチングデータテーブルを参照して、追尾対象の目標の目標識別データに関連付けられた目標位置を抽出することを特徴とする。
【０００８】
請求項２記載の発明は、請求項１記載のものにおいて、過去の２つ以上のタイミングで取得され同じ目標識別データと関連付けられた目標位置に基づき該目標識別データに対応する目標速度を求める速度情報算出手段をさらに備え、
前記マッチング処理手段は、前記目標速度に基づきあるタイミングにおける目標予測位置を求めて、該目標予測位置とあるタイミングで取得された前記目標位置との対応付けを行うことを特徴とする。
【０００９】
請求項３記載の発明は、請求項１または２記載の前記マッチング処理手段が、複数の目標位置と複数の被マッチング目標位置に対してその距離の総和が最小となるような対応付けを行うことを特徴とする。
【００１０】
請求項４記載の発明は、請求項１または２記載の前記マッチング処理手段が、複数の目標位置と複数の被マッチング目標位置に対してその距離及び目標の大きさの差異の乗算の総和が最小となるような対応付けを行うことを特徴とする。
【００１１】
請求項５記載の発明は、請求項３または４記載の前記マッチング処理手段が、目標位置と被マッチング目標位置との距離が所定距離を超えている場合には、該目標位置と該被マッチング目標位置との対応付けは行わないことを特徴とする請求項３または４記載の目標追尾装置。
【００１２】
請求項６記載の発明は、請求項５記載の前記所定距離は、前記目標信号より得られる目標の大きさに応じて変化させることを特徴とする。
【００１３】
請求項７記載の発明は、請求項５記載の前記所定距離は、目標速度の大きさに応じて変化させることを特徴とする。
【００１５】
請求項８記載の発明は、請求項１ないし７のいずれか１項に記載の前記マッチング処理手段が、
複数の目標位置と複数の被マッチング目標位置とのいずれか一方の位置（以下、第１目標位置）に対して、その一方の位置から最も近い他方の位置（以下、第２目標位置）を仮対応付ける仮対応付け手段と、
複数の第１目標位置が仮対応付けられた第２目標位置に対して、一つの第１目標位置を仮対応付けて、別の第１目標位置を別の第２目標位置に仮対応付ける仮対応付け選択手段と、
全ての第２目標位置に１以下の第１目標位置が仮対応付けられたときに、仮対応付けられた第１目標位置と第２目標位置との対応付けを確定する確定手段と、
を備えることを特徴とする。
【００１６】
請求項９記載の発明は、請求項８記載の前記仮対応付け選択手段が、
同じ第２目標位置に対して仮対応付けられた２つの第１目標位置Ｘ、Ｙに対して、それぞれ仮対応付けられた第２目標位置をＸ１、Ｙ１とし、Ｘ１、Ｙ１の次に近い第２目標位置をＸ２、Ｙ２としたときに、Ｘｄ＝｜Ｘ−Ｘ１｜＋｜Ｙ−Ｙ２｜とＹｄ＝｜Ｘ−Ｘ２｜＋｜Ｙ−Ｙ１｜（ここで、｜Ａ−Ｂ｜はＡ、Ｂ間の距離を表す）の大小を比較する手段と、
Ｘｄ≦Ｙｄである場合には、Ｘの仮対応点としてＸ１を選択し、Ｙの仮対応点としてＹ２を選択し、Ｘｄ＞Ｙｄである場合には、Ｙの仮対応点としてＹ１を選択し、Ｘの仮対応点としてＸ２を選択する手段と、
を備えることを特徴とする。
【００１７】
請求項１０記載の発明は、時系列的に取得された目標信号から得られる複数の目標位置に関するデータに基づいて目標の追尾を行う目標追尾方法において、
あるタイミングで取得された複数の目標位置と、過去のタイミングで取得された複数の目標位置または過去のタイミングで取得された目標位置に基づき求められたあるタイミングにおける複数の目標予測位置である複数の被マッチング目標位置との対応付けを行うマッチング処理工程と、
マッチング処理工程で被マッチング目標位置と対応付けられた目標位置を、被マッチング目標位置に対応する目標識別データと関連付けて格納する工程と、
過去のタイミングで取得された目標位置または過去のタイミングで取得された目標位置及び過去の２つ以上のタイミングで取得された目標位置に基づき求めた目標速度から目標の予測位置を求めて、該予測位置を中心としたウィンドウを追尾対象の目標毎に計算するウィンドウ計算工程と、
ウィンドウ計算工程で計算されたウィンドウ内にあるデータから目標を抽出する目標抽出工程と、を備え、
前記目標抽出工程は、ウィンドウ内に目標とするべきデータが２つ以上存在するとき、またはウィンドウ内に目標とするべきデータが存在しないときに、追尾対象の目標の目標識別データに関連付けて格納された目標位置を抽出することを特徴とする。
【００１８】
請求項１１記載の発明は、請求項１０記載の方法において、
前記マッチング処理工程が、
（ａ）  対象となる複数の目標位置と対象となる複数の被マッチング目標位置とのいずれか一方の位置（以下、第１目標位置）に対して、その一方の位置から最も近い他方の位置（以下、第２目標位置）を仮対応付ける工程と、
（ｂ）  第２目標位置に仮対応付けられた第１目標位置の数を求める工程と、
（ｃ）  同じ第２目標位置に仮対応付けられた第１目標位置の数が２以上である場合には、その数が２つである場合に、その２つの第１目標位置Ｘ、Ｙを選択し、その数が３以上である場合に、任意の２つの第１目標位置Ｘ、Ｙを選択する工程と、
（ｄ）  ２つの第１目標位置Ｘ、Ｙに対して、それぞれ仮対応付けられた第２目標位置をＸ１、Ｙ１とし、Ｘ１、Ｙ１の次に近い第２目標位置をＸ２、Ｙ２としたときに、Ｘｄ＝｜Ｘ−Ｘ１｜＋｜Ｙ−Ｙ２｜とＹｄ＝｜Ｘ−Ｘ２｜＋｜Ｙ−Ｙ１｜（ここで、｜Ａ−Ｂ｜はＡ、Ｂ間の距離を表す）の大小を比較する工程と、
（ｅ）  Ｘｄ≦Ｙｄである場合には、Ｘの仮対応点としてＸ１を選択し、Ｙの仮対応点としてＹ２を選択し、Ｘｄ＞Ｙｄである場合には、Ｙの仮対応点としてＹ１を選択し、Ｘの仮対応点としてＸ２を選択する工程と、
（ｆ）  （ｂ）、（ｄ）、（ｅ）の工程を、全ての第２目標位置に１以下の第１目標位置が仮対応付けられるまで繰り返す工程と、
（ｇ）  全ての第２目標位置に仮対応付けられた第１目標位置の数が１以下となったときに、仮対応付けられた第１目標位置と第２目標位置との対応付けを確定する工程と、
を備えることを特徴とする。
【発明の効果】
【００１９】
本発明によれば、１つの目標にのみ着目してウィンドウを設定する従来と異なり、複数の目標位置と複数の被マッチング目標位置との対応付けを行うので、ウィンドウ内といった局所的な情報だけでは得られなかった周辺の目標の運動を考慮した目標追尾が可能となり、正確な追尾が可能となる。また、ウィンドウを設定する手段と組み合わせて精度を上げることも可能である。
【００２０】
これによって輻輳地域においても、乗り移りを抑制した正確な追尾が可能となる。また、追尾中の目標が急速な加減速や旋回を行った場合でも、全体として最適な組み合わせを求めるので、ロストを抑制した正確な追尾が可能となる。つまり、ある注目する目標が急速な加減速や旋回を行っても、周りの目標の追尾が継続できていれば、その注目した目標を容易に見つけることができるようになる。
【図面の簡単な説明】
【００２１】
【図１】本発明による目標追尾方法に従った目標追尾装置の第１実施形態のブロック図である。
【図２】本発明による目標追尾方法に従った目標追尾装置の第２実施形態のブロック図である。
【図３】本発明による目標追尾方法に従った目標追尾装置の第３実施形態のブロック図である。
【図４】マッチング処理の具体的なアルゴリズムを表すフローチャートである。
【図５】図４によるマッチング処理の具体例である。
【図６】図５の例において、従来のウィンドウ方式との違いを説明する説明図である。
【図７】目標の予測位置を求めるための一例を表す説明図である。
【図８】目標の位置を求めるための様々な例を表す説明図である。
【発明を実施するための形態】
【００２２】
以下、図面を用いて本発明の実施の形態を説明する。
【００２３】
(第１実施形態)
図１は、本発明による目標追尾方法に従った目標追尾装置１０の第１実施形態のブロック図を表している。この例の目標追尾装置１０は、ＡＲＰＡ装置に搭載されたものとしている。しかしながら、目標追尾装置１０は、周囲の海上または空中の船舶または飛行体といった目標を監視するために移動体自体に搭載する他に、海上または空中の船舶または飛行体といった目標を監視するために地上等に定置することができる。
【００２４】
図示したように、目標追尾装置１０は、大別して、信号受信ブロック１２、目標追尾ブロック１４、目標情報処理ブロック１６、画面表示部１８及びユーザインターフェース部２０と、を有している。
【００２５】
信号受信ブロック１２には、外部から目標エコー信号または目標信号が入力される。信号受信ブロック１２に入力される信号としては、目標位置検出手段であるレーダ装置から出力されたビデオ信号とすることができる。この場合、レーダアンテナが１つのパルスを送信し、それに対する反射波を受信するときの、その期間を１スイープと呼び、レーダアンテナが１回転する期間を１スキャンと呼ぶ。各パルスの真方位は、レーダアンテナから得られる方位（相対方位）に、ジャイロ等で得られるレーダアンテナのヘディング方位を加えることで求められる。
【００２６】
そして、信号受信ブロック１２は、入力されたビデオ信号を検波する信号受信部２２と、ビデオ信号に対してノイズ除去などの各種信号処理を行う信号処理部２４と、１スキャン毎に更新して１または複数スキャン分のエコーデータを格納するエコーデータメモリ２６を有する。
【００２７】
ここで、信号処理部２４は、エコーデータメモリ２６に格納する前処理として、複数スキャン分のエコーデータを格納することで、スキャン相関処理を行ってノイズ等を除去することも可能である。
【００２８】
尚、信号受信ブロック１２は、目標追尾装置１０が直接接続された目標位置検出手段であるレーダ装置からの信号を受信する他、目標追尾装置１０とは異なる移動体または目標追尾装置１０とは離れた地上に設置された目標位置検出手段であるレーダ装置からの信号を何らかの手段で受信してもよい。スレーブレーダ装置からの信号を直接に受信してもよく、また、レーザ装置またはその他の任意の目標位置検出手段から得られた信号をデジタルデータで受信してもよく、または、デジタルデータを任意の通信ネットワークを介して受信してもよい。いずれにしても、目標位置検出手段に対してある領域に着目したときに、その領域内に存在する目標を表す目標信号が時系列的に取得されて、その目標信号が信号受信ブロック１２に入力されればよい。また、信号受信ブロック１２に入力される前処理として、目標位置検出手段からの目標信号に対して、地図データを用いて地図データに基づき陸等の固定物を除去する処理を行った後に、信号受信ブロック１２に目標信号を入力するようにしてもよい。また、目標を検出する前に、ＡＩＳデータを用いて、追尾する必要のない既知の目標を除去する処理を行うようにしてもよい。
【００２９】
目標追尾ブロック１４は、目標検出部３２と、目標位置テーブル３４と、目標位置マッチング処理部３６と、マッチングデータテーブル４０と、を有する。
【００３０】
目標検出部３２は、エコーデータメモリ２６に格納された１スキャン分のデータから目標の位置を検出する。そのためにまず、目標エコーを抽出する。その抽出する方法としては、次の方法が考えられる。
【００３１】
（１）デジタルデータをある閾値を用いて二値化した後、距離方向、及び方位方向にある一定以上の連続性のある部分を検出して目標エコーとするか、または、ラベリングなどの画像処理手法によってエコーの塊の分離を行い、ある一定サイズ以上の塊を目標エコーとする。
（２）デジタルデータの値の２次元的な傾きを求め、極大となっている部分を目標位置とする。または、デジタルデータの値の変化量（傾き）が一定以上の部分を、目標エコーの境界として検出する。
（３）アナログ信号からそのまま目標を検出する方法として、スプリットゲート方式（吉田孝「改訂レーダ技術」、p.264、電子情報通信学会(1996)）を用い、アナログ信号を進みゲートと遅れゲートを通過させ、信号の中央部で両ゲートの差分が０となることを利用して目標エコーを検出する。
【００３２】
また、検出した目標エコーから目標位置を決定するには、次の方法が考えられる（図８参照）。
（１）目標エコーの塊の重心位置を目標位置とする（図８（ａ））。
（２）目標エコーの距離方向及び方位方向の中心位置を目標位置とする（図８（ｂ））。
（３）目標エコーの方位方向の中心位置、かつ距離方向の最もアンテナに近い位置を目標位置とする（図８（ｃ））。
【００３３】
また、信号受信ブロック１２に入力される前処理を行う代わりとして、目標検出部３２において目標を検出する際に、地図データを用いて地図データから陸等の固定物を除去する処理を行うようにしてもよい。尚、このように目標検出部３２によって１スキャン分のデータから全ての目標を自動的に抽出することができるので、ユーザが追尾目標を特定する必要はない。また、目標を検出する際に、ＡＩＳデータを用いて、追尾する必要のない既知の目標を除去する処理を行うようにしてもよい。
【００３４】
目標検出部３２で検出された１つ以上の目標位置は、目標位置テーブル３４に格納される。目標位置テーブル３４には、少なくとも最新の２スキャン分の目標位置が格納される。
【００３５】
目標位置マッチング処理部３６は、目標位置テーブル３４で格納された２スキャン分の目標位置を用いてマッチングを行う。具体的には、２スキャン間の目標位置に対して全ての目標の移動距離の和が最小になるようにマッチングを行う。または、和が最小となるような近似マッチング解を求めて、そのマッチングを行う。即ち、現スキャン（ｔ）の目標ｉの位置を（ｘi,t，ｙi,t）とし、その前のスキャンの目標位置を（ｘi,t-1，ｙi,t-1）としたときに、その位置の差ｄiの全目標の総和
【００３６】
【数１】
0005614969.tif 000002
が最小となるような組み合わせ（ｘit，ｙit）と（ｘi,t-1，ｙi,t-1）を見つけて、同じ目標にその位置情報を割り当てる。尚、３次元の場合でも同様に３次元に拡張すればよい。具体的な最小となるような組み合わせを近似的に見つけるマッチングアルゴリズム例については、後述する。
【００３７】
こうして、組み合わせを検出して対応付けると、マッチングデータテーブル４０において、目標を識別する目標識別データ及び時間に相当するそのスキャンタイミングと関連付けてその位置が格納される。ここで、目標識別データとしては、ある目標の履歴を他の目標の履歴と識別することができるものであれば、何でもよい。
【００３８】
目標情報処理ブロック１６は、マッチングデータテーブル４０で格納された目標の位置データから対象目標の針路・速力を計算する目標情報計算部（平滑部）４２と、その計算結果を格納する目標追尾情報テーブル４４と、目標データ検索部４６と、を有する。
【００３９】
目標情報計算部４２で行う針路・速力計算は、α−βフィルタ等を用いて計算することができる。
【００４０】
画面表示部１８には、前記目標追尾情報テーブル４４に格納された情報から目標が表示される。ユーザインターフェース部２０から、画面表示部１８に対して所望の目標が指定される、あるいは、予め指定された領域内に目標が侵入する、あるいは、目標情報処理ブロック１６により自動判断されることによって、目標データ検索部４６は、その目標について目標追尾情報テーブル４４で格納された情報、つまり針路・速力の情報を検索し、その情報を画面表示部１８に送ることで、画面表示部１８にその情報が表示される。
【００４１】
（第２実施形態）
図２は、本発明による目標追尾方法に従った目標追尾装置１０の第２実施形態のブロック図を表している。この第２実施形態は、第１実施形態と大部分が同じであるため、異なる部分のみ説明する。
【００４２】
即ち、この実施形態では、目標追尾ブロック１４が、速度情報算出部３８を更に備えている。
【００４３】
速度情報算出部３８は、マッチングデータテーブル４０の情報を用いて、各目標の速度（速力と針路）を算出する。即ち、速度は、過去の２つ以上のスキャンに亘る目標の位置を用いて求めることができる。必要に応じてα−βフィルタ等を用いてもよく、または、目標が一定の速度を維持するものとして、過去の２つのスキャンにおける目標の位置の差から単純に速度を求めてもよい。算出された速度は、目標識別データと関連付けてマッチングデータテーブル４０に格納される。
【００４４】
このように速度情報算出部３８を設けた目標位置マッチング処理部３６は、２スキャン間の目標位置に対して全ての目標の移動距離の和が最小になるようにマッチングを行う代わりに、現スキャン（ｔ）の各目標ｉの位置（ｘit，ｙit）と、速度情報算出部３８で算出した速度を用いた現スキャン（ｔ）における各目標の予測位置（ｘit＾，ｙ＾it）との距離の総和が最小になるようにマッチングを行う。この場合、各目標ｉの位置（ｘit，ｙit）と各目標の予測位置（ｘｓit＾，ｙ＾it）との位置の差ｄi＾は小さいと予測できるので、マッチングの精度を上げることができる。
【００４５】
以上のように第１及び第２実施形態による目標追尾装置においては、ウィンドウを設定せずに、その代わりに、全目標を対象として、過去の目標位置または予測目標位置と、現目標位置との間でのマッチングを行うようにしている。これに対して、従来は、個々の目標に対してそれぞれウィンドウを設定することで、個々の目標に対してそれぞれ独立的に追尾を行っている。そのため、１つのウィンドウ内に複数の目標が存在する場合もあり、またはウィンドウから目標エコーが外れてしまう場合もある。
【００４６】
しかしながら、本発明では、全ての目標を対象として追尾を行い、その移動量の和が最小となるように対応付けを行うので、ウィンドウ内といった局所的な情報だけでは得られなかった周辺の目標の運動を考慮した目標追尾が可能となり、正確な追尾が可能となる。これによって輻輳地域においても、乗り移りを抑制した正確な追尾が可能となる。また、追尾中の目標が急速な加減速や旋回を行った場合でも、全体として最適な組み合わせを求めるので、ロストを抑制した正確な追尾が可能となる。つまり、ある注目する目標が急速な加減速や旋回を行っても、周りの目標の追尾が継続できていれば、その注目した目標を容易に見つけることができる。
【００４７】
また、複数のスキャン間にまたがる信号を使用するために、信号処理部２４でスキャン相関処理を行うことで、ノイズ除去を行うことも可能になる。
【００４８】
尚、以上の第１及び第２実施形態では、１つのスキャン全体をマッチング対象としたが、これに限るものではなく、１つのスキャンを複数の領域に分けて、それぞれの領域においてマッチングを行うようにしてもよい。
【００４９】
(第３実施形態)
図３は、本発明による目標追尾方法に従った目標追尾装置１０の第３実施形態のブロック図を表している。第１または２実施形態の同一・同様の部材は同一の符号を付して、その詳細説明を省略する。
【００５０】
この実施形態では、ウィンドウを設定する方法とマッチングとを組み合わせて、追尾の精度を向上させるものである。
【００５１】
即ち、目標追尾装置１０は、大別して、信号受信ブロック５２、目標追尾ブロック５４、目標情報処理ブロック５６、画面表示部１８、ユーザインターフェース部２０、自動目標捕捉部６２及び補助追尾ブロック６４を有している。
【００５２】
信号受信ブロック５２は、信号受信ブロック１２と同様であり、外部からの目標エコー信号または目標信号が入力され、信号受信部２２と、信号処理部２４と、１スイープ分のエコーデータを格納するスイープデータメモリ６６とを有する。
【００５３】
目標追尾ブロック５４は、スイープデータ抽出部７０と、追尾ウィンドウメモリ７２と、目標抽出部７４と、目標位置算出部７６と、目標位置テーブル７８と、速度情報算出部８０と、ウィンドウ予測位置計算部８２と、ウィンドウ位置テーブル８４と、ウィンドウ初期位置計算部８６と、を有する。
【００５４】
スイープデータ抽出部７０は、後述のウィンドウ位置テーブル８４に従い、スイープデータメモリ６６からウィンドウ位置にあるエコーデータを抜き出し、追尾ウィンドウメモリ７２に格納する。
【００５５】
追尾ウィンドウメモリ７２に１つのウィンドウ分のデータの保存が完了すると、目標抽出部７４は、追尾ウィンドウメモリ７２に保存されたウィンドウ内の目標の抽出を行い、目標位置算出部７６で抽出した目標の位置を求める。目標抽出部７４は、必要に応じて後述の補助追尾ブロック６４で求めた目標を使用する。
【００５６】
求めた目標位置は、目標識別データ及び時間に相当するスキャンタイミングと関連付けて目標位置テーブル７８に格納される。速度情報算出部８０は、目標位置テーブル７８で格納された目標位置及び前回の目標位置から目標速度を求める。ウィンドウ予測位置計算部８２は、求めた目標位置及び速度情報算出部８０で求めた目標速度を用いて、次回の目標の予測位置を求め、さらにその予測位置を中心とするウィンドウを求める。また、ウィンドウ初期位置計算部８６は、初期捕捉時に、今回の目標位置を中心として、大きさの大きいウィンドウを求める。求めたウィンドウは、ウィンドウ位置テーブル８４に目標識別データと関連付けて格納される。
【００５７】
目標情報処理ブロック５６は、目標位置テーブル７８で格納された目標の位置データから対象目標の針路・速力を計算する目標情報計算部４２と、その計算結果を格納する目標追尾情報テーブル４４と、を有する。
【００５８】
追尾するべき目標の捕捉は、ユーザインターフェース部２０からユーザによって指定されるか、または自動目標捕捉部６２によってスイープデータメモリ６６からのエコーデータから、所定の条件を満足する目標が指定されることによって、行われる。
【００５９】
次に、補助追尾ブロック６４は、信号処理部２４からのエコーデータを、１スキャン毎に更新して１または複数スキャン分を格納するエコーデータメモリ９０と、目標検出部３２と、目標位置テーブル３４と、目標位置マッチング処理部３６と、必要に応じて速度情報算出部３８と、マッチングデータテーブル４０と、を有する。
【００６０】
補助追尾ブロック６４は、原則的に、第１実施形態または第２実施形態の目標追尾ブロック１４と同じであり、２スキャン分の目標位置のマッチングを行うか、または、現スキャンの各目標の位置と各目標の予測位置とのマッチングを行い、マッチングデータテーブル４０には、その目標識別データと目標位置とが格納される。
【００６１】
前記目標抽出部７４は、各目標識別データに対応して設定されたウィンドウ内で抽出された目標が１つであれば、そのまま抽出した目標の位置を目標識別データに対応するものとして求めればよい。しかしながら、所定の条件を満足する場合に、補助追尾ブロック６４のマッチングデータテーブル４０を参照して、目標の抽出を行うことができる。
【００６２】
所定の条件を満足する場合とは、例えば、ウィンドウ内に目標が２つ以上抽出された場合、またはウィンドウ内に目標が１つも抽出されない場合とすることができる。この場合に、補助追尾ブロック６４のマッチングデータテーブル４０において、マッチングされたデータを用いて問題とする目標識別データに対応する目標の位置を検出することができ、これによって、目標の乗り移りやロストを抑制することができる。
【００６３】
補助追尾ブロック６４は、目標追尾ブロック５４と独立的に演算処理を行うことができる。ここで、各ブロックにおいて、目標を識別するものとして使用される目標識別データは、共通としてもよく、または別個のものとしてもよい。別個のものとする場合には、互いの目標識別データの対応付けを、その他の情報、例えば、過去のスキャン時の目標位置の一致などから求めることができる。
【００６４】
そして、補助追尾ブロック６４は、常に、全ての目標に対してマッチングを行ってもよい。または、補助追尾ブロック６４の目標位置マッチング処理部３６によるマッチング処理は、限定された範囲でのマッチングを行うようにしてもよい。限定された範囲とは、前記目標抽出部７４において、前記所定の条件となる任意の目標のウィンドウ内に２つ以上の目標が抽出された場合、または該ウィンドウ内に全く目標が抽出されない場合に、そのウィンドウを含み該ウィンドウよりやや広めの領域に対してマッチングを行うようにしてもよい。また、マッチングを行う際に、前のスキャンの目標位置または、目標の速度を、目標位置テーブル７８から取得してもよい。
【００６５】
この実施形態によれば、ウィンドウを用いた追尾を用いることで、処理を簡素化すると共に、ウィンドウを用いた場合の欠点である乗り移り、ロストという問題が発生したときに、補助追尾ブロックによるマッチング結果を用いて目標を特定することで、ウィンドウを用いた場合の欠点を補うことができる。補助追尾ブロックにおいて、乗り移りやロストが発生している領域においてマッチング処理を行うようにすることで、処理負荷を低減させることもできる。
【００６６】
この実施形態では、目標抽出部７４でウィンドウ内の目標抽出を行う際にウィンドウ毎に抽出条件を変更することができ、目標の抽出をより確実にすることができる。例えば、大きな目標の影となってしまう領域においては、その信号は低くなるので、その場合には、該当するウィンドウにおける目標抽出の閾値を低く設定するなどのフレキシビリティを持たせることができる。
【００６７】
（マッチング処理の具体例）
次に、マッチング処理の具体例について以下に説明する。
ここで、マッチングは、ある時刻（タイミング）ｔ１にｎ個の点群が存在し、時刻（タイミング）ｔ２（ｔ２≠ｔ１）にｎ±α個の点群が存在するときに、時刻ｔ１の点とｔ２の点で、最も確からしい点を対応付けることである。
【００６８】
各点が情報として位置情報しか持たない場合には、マッチングの要素として持つのは、各点間の距離情報だけとなる。このとき、最も確からしいと思われるマッチングの解の一つは、対応付けられた点同士の距離ｄの総和（式（１）参照）が最小となる場合であると考えられる。
【００６９】
全ての組み合わせを計算して、総和が最小となる組み合わせを求めることも考えられるが、計算量が膨大となる欠点がある。そこで、この計算量を低減するアルゴリズムとして図４に示す次の例を提案する。
【００７０】
（１）  時刻ｔ＝ｔ１での点の集合Ｕ１＝｛Ａ１，Ａ２，・・・，Ａｎ｝と、時刻ｔ＝ｔ２での点の集合Ｕ２＝｛Ｂ１，Ｂ２，・・・，Ｂｎ｝を得る（ステップＳ１１）。
（２）  Ｕ２の各点それぞれについて、最も近いＵ１の点を一つ選びそれを仮対応付けとする（ステップＳ１２）。
（３）  Ｕ１の各点それぞれについて仮対応付けられたＵ２の点の数を調べる（ステップＳ１３）。
（４）  あるＵ１の点に対応付けられたＵ２の点が２つのときには、仮対応付けられたＵ２の２点を便宜上Ｘ、Ｙと称する（ステップＳ１４）。また、あるＵ１の点に対応付けられたＵ２の点が３つ以上のときには、３点以上の点から任意の２点を選び便宜上Ｘ、Ｙと称する（ステップＳ１５）。一方、Ｕ１の点に対応付けられたＵ２の点が１つである場合には、全ての仮対応付けを確定して、マッチングを終了する（ステップＳ２０）。
（５）  仮対応付けられたＵ２の２点Ｘ、Ｙに対して、点Ｘから最も近いＵ１の点をＸ１、その次に近いＵ１の点をＸ２とする。同様に、点Ｙから最も近いＵ１の点をＹ１、その次に近いＵ１の点をＹ２とする。ここで、Ｘ１とＹ１は同じ点である。尚、Ｕ１の点を選ぶときに、後述のステップＳ１８、１９で除外した点は除くものとする（ステップＳ１６）。
（６）  Ｘｄ＝｜Ｘ−Ｘ１｜＋｜Ｙ−Ｙ２｜とＹｄ＝｜Ｘ−Ｘ２｜＋｜Ｙ−Ｙ１｜の大小を比較する（ステップＳ１７）。ここで｜Ａ−Ｂ｜は、点Ａ、Ｂ２点間の距離を表す。
（７）  Ｘｄ≦Ｙｄである場合には、Ｘの仮対応点としてＸ１を選択し、Ｙの仮対応点としてＹ２を選択する。以降は、Ｙの対応点としてＹ１の点は除外する（ステップＳ１８）。
（８）  Ｘｄ＞Ｙｄである場合には、Ｙの仮対応点としてＹ１を選択し、Ｘの仮対応点としてＸ２を選択する。以降は、Ｘの対応点としてＸ１の点は除外する(ステップＳ１９)。
（９）  ステップＳ１３〜ステップＳ１９を、全てのＵ１の点に対応付けられたＵ２の点の数が１つとなるまで、繰り返す。
【００７１】
（１）〜（９）のアルゴリズムに従った具体的処理例を図５に基づき説明する。
図５（ａ）において、時刻ｔ２の白抜き点○Ａ〜Ｄを、時刻ｔ１の黒丸点●ア〜エにマッチングすることを考えると、まず、白抜き点○Ａ〜Ｄを、最も距離の近い時刻ｔ１の黒丸点●に仮対応付ける。●エに対して、○Ｂ、○Ｃ、○Ｄの３点が仮対応付けられるので、その中から任意の点として、○Ｂと○Ｃとを選択する。
【００７２】
そして、｜Ｂ−Ｂ１｜＋｜Ｃ−Ｃ２｜と｜Ｂ−Ｂ２｜＋｜Ｃ−Ｃ１｜との距離を比較する。後者の方が短いので、Ｂ２（●ウ）とＣ１（●エ）とを○Ｂ、○Ｃに仮対応付ける(図５（ｂ）)。ここで、○Ｂは●エの対応付けの候補としては落選したことになり、以降、●エは○Ｂの対応点としては除外される。
【００７３】
依然として、●エに対して○Ｃと○Ｄとが仮対応付けされているので、｜Ｃ−Ｃ１｜＋｜Ｄ−Ｄ２｜と｜Ｃ−Ｃ２｜＋｜Ｄ−Ｄ１｜との距離を比較する。後者の方が短いので、Ｃ２（●ウ）とＤ１（●エ）とを○Ｃと○Ｄに仮対応付ける(図５（ｃ）)。ここで、○Ｃは●エの対応付けの候補としては落選したことになり、以降、●エは○Ｃの対応点としては除外される。
【００７４】
すると、図５（ｄ）に示すように、●ウに対して、○Ｂ、○Ｃの２点が仮対応付けられるので、｜Ｂ−Ｂ１｜＋｜Ｃ−Ｃ２｜と｜Ｂ−Ｂ２｜＋｜Ｃ−Ｃ１｜との距離を比較する。ここで、最も近い点である●エは、前述のように対応点としては除外されている。後者の方が短いので、Ｂ２（●ア）とＣ１（●ウ）とを○Ｂ、○Ｃに仮対応付ける。以降、●ウは○Ｂの対応点としては除外される。
【００７５】
すると、図５（ｅ）に示すように、●アに対して、○Ａ、○Ｂの２点が仮対応付けられるので、｜Ａ−Ａ１｜＋｜Ｂ−Ｂ２｜と｜Ａ−Ａ２｜＋｜Ｂ−Ｂ１｜との距離を比較する。ここで、○Ｂについては、●エ、●ウは、前述のように対応点としては除外されている。前者の方が近いので、Ａ１（●ア）とＢ２（●イ）とを○Ａ、○Ｂに仮対応付ける。
【００７６】
こうして、複数対応付けられた黒丸点●がなくなったので、仮対応付けを確定し、処理を終了する。
【００７７】
以上のアルゴリズムで少ない計算量でマッチングを行うことができる。
【００７８】
図５で示した例の場合に、仮にウィンドウを設定して追尾を行ったとすると、初期捕捉時に、黒丸点●エを中心にウィンドウを設定した場合、図６に示すウィンドウＷ１ではウィンドウサイズが小さ過ぎるために、マッチングされた○Ｄが入らず、間違った追尾がなされてしまう。反対に、図６に示すウィンドウＷ２ではウィンドウサイズが大き過ぎるために、○Ａ、○Ｂ、○Ｃ、○Ｄが全て入ってしまい、通常は中心部に近い○Ｃが選択されてしまい、同様に間違った追尾がなされてしまう。マッチングによって、このような不具合を一掃することができ、正しい追尾を行うことができる。
【００７９】
以上の図４〜図５は、前回のスキャン時と今回のスキャン時における目標のマッチングを行う例であったが、今回のスキャン時の目標位置と目標の予測位置とのマッチングを行う場合は、例えば、図７に示すように、前々回のスキャン時の目標位置と前回のスキャン時の目標位置と、その速度が一定であることを前提として、予測位置ｔ１’を求めることができる。そして、予測位置と今回のスキャン時の目標位置との間でのマッチングについては、図４〜図５に示したアルゴリズムと同じ手法により（図５において●または○のいずれか一方を目標の予測位置、いずれか他方を今回のスキャン時の目標の位置とする）、計算量が少なく、且つより精度の高いマッチングを行うことができる。
【００８０】
以上の説明においては、前回のスキャン時と今回のスキャン時の目標数が一定であることを前提として行ったが、実際には、目標数が増減する。よって、そのような場合に備えて、マッチングする範囲を予め決めておき、その範囲内でなければマッチングを行わないように制約条件を付加することで、対応付けができない点を許すようにしてもよい。
【００８１】
具体的には、図４のアルゴリズムのステップＳ１２の仮対応付けにおいて、「最も近いＵ１の点が所定距離ｄ以上である場合には、仮対応付けしない」という条件を加え、ステップＳ１８及びＳ１９において、ＹとＹ２、ＸとＸ２との距離が所定距離ｄ以上である場合には、仮対応付けしない」という条件を加え、ステップＳ２０の仮対応付けを確定する条件としては、全てのＵ１の点に対応付けられたＵ２の点の数が１つ以下とすることができる。
【００８２】
これによって、明らかに離れた点同士が対応付けされてしまうことを防ぐことができ、時刻ｔ１とｔ２との数が同数でない場合でも確からしい対応付けを行うことができる。
【００８３】
目標数の増減の原因としては、次のような場合が考えられる。
・目標数が増加する場合
（１）  目標が処理範囲の外から進入してきた場合
（２）  ノイズやクラッタを目標と誤検出した場合
（３）  目標エコー信号が不安定となり、１つの目標が２つに見えてしまう場合
このような場合、マッチングされなかった新しい位置を、新規目標の位置として、追尾を開始するようにすればよい。ノイズの位置はランダムに近いので、ノイズが新規に捕捉されてしまっても、通常は数スキャンの間に追尾不可能になるはずである。同様に、１つの目標が一時的に２つに見えてしまった場合でも、片方は追尾が継続され、もう片方はロストとなる。しかしながら、追尾中の目標がノイズとマッチングしてしまうと、次回のスキャン時にロストになる可能性がある。
【００８４】
このようにノイズとマッチングすることを防ぐために、マッチングの条件として、移動距離を条件とするだけでなく、エコーサイズの大きさを条件とすることができる。即ち、現スキャン（ｔ）の目標ｉの位置を（ｘi,t，ｙi,t）、エコーの大きさをｓi,tとし、その前のスキャンの目標位置を（ｘi,t-1，ｙi,t-1）、エコーの大きさｓi,t-1としたときに、その位置の差ｄi及びサイズの差Δｓiとを乗算し、その全目標の総和
【００８５】
【数２】
0005614969.tif 000003
が最小となるようなマッチングを行うようにしてもよい。目標位置が３次元である場合も同様に３次元に拡張すればよい。
【００８６】
・目標数が減少する場合
（１）  大きな目標の影に隠れ、エコーが弱くなり、検出ができなくなった場合
（２）  他のエコーと一体化してしまった場合
（３）  目標が処理範囲の外へ出た場合
このような場合、マッチングされなかった前の目標を追尾不可能(ロスト)として扱えばよい。
【００８７】
但し、そのように扱うと、（１）、（２）の原因でも目標が検出できなかった場合にすぐにロストになってしまうので、一時的に検出できなかった場合、即ちマッチングできなかった場合であっても、数スキャンの間、仮追尾を行うことが考えられる。仮追尾は、図７に示すような予測位置を求めることで行うことができる。そうすることで、目標エコーが復帰したときに、継続して追尾を行うことができる。
【００８８】
また、以上の制約条件を一定に定める他に、マッチングを行う範囲を、追尾状態や船舶情報を用いて動的に変化させることも考えられる。例えば、初期捕捉時に速度情報が得られないときには、マッチング範囲を広くし（所定距離ｄを大きくする）、速度情報が得られた後は、マッチング範囲を狭く（所定距離ｄを小さくする）することが考えられる。または、速力が大きな目標は高速船であり高速旋回を行う可能性があるとして、マッチング範囲を広くしてもよい。または、エコーサイズが大きな目標は大型船であり高速旋回を行う可能性は低いとしてマッチング範囲を狭くしてもよい。
【符号の説明】
【００８９】
１０  目標追尾装置
２６  エコーデータメモリ（データメモリ）
３６  目標位置マッチング処理部（マッチング処理手段）
３８  速度情報算出部（速度情報算出手段）
４０  マッチングデータテーブル
７４  目標抽出部（目標抽出手段）
８２  ウィンドウ予測位置計算部（ウィンドウ計算手段）
８６  ウィンドウ初期位置計算部（ウィンドウ計算手段）
【図４】
0005614969.tif 000004
【図５】
0005614969.tif 000005
【図６】
0005614969.tif 000006
【図７】
0005614969.tif 000007
【図８】
0005614969.tif 000008
【図１】
0005614969.tif 000009
【図２】
0005614969.tif 000010
【図３】
0005614969.tif 000011