JP2007292553A - Multiple target tracking system - Google Patents

Multiple target tracking system Download PDF

Info

Publication number
JP2007292553A
JP2007292553A JP2006119462A JP2006119462A JP2007292553A JP 2007292553 A JP2007292553 A JP 2007292553A JP 2006119462 A JP2006119462 A JP 2006119462A JP 2006119462 A JP2006119462 A JP 2006119462A JP 2007292553 A JP2007292553 A JP 2007292553A
Authority
JP
Japan
Prior art keywords
observation
track
value
temporary
unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2006119462A
Other languages
Japanese (ja)
Other versions
JP5047534B2 (en
Inventor
Yasushi Obata
康 小幡
Takashi Matsuzaki
貴史 松崎
Masayoshi Ito
正義 系
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2006119462A priority Critical patent/JP5047534B2/en
Publication of JP2007292553A publication Critical patent/JP2007292553A/en
Application granted granted Critical
Publication of JP5047534B2 publication Critical patent/JP5047534B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To realize a multiple target tracking system capable of obtaining an accurately tracking result in the allowable processing time. <P>SOLUTION: The multiple target tracking system comprises an observation value selecting section 1 for determining a combination of observed values and temporary tracks used for performing predicting computation; a predicting section 22 for computing prediction of a temporary track, based on the observed value acquired time about the combination of observed values and temporary tracks determined by the observation value selecting section 1; a gate inside/outside determining section 23 for determining as to whether the temporary tracks can be associated with the observed values, by comparing the observed values with the predicted values of the temporary tracks computed by the predicting section 22; a track likelihood evaluating section 24 for computing the likelihood of the temporary track, when it is associated with the observed value; a correlation determining section 25 for determining the association of the temporary tracks and the observed values, by using the likelihood as an index which is computed by the track likelihood evaluating section 24; a smoothing section 26 for computing smoothing parameters, by using a predicted parameter of the temporary tracks and the observed value associated by the correlation determining section 25; and a track determining section 27 which determines a temporary track, having a high likelihood among target tracks. <P>COPYRIGHT: (C)2008,JPO&INPIT

Description

捜索レーダ等のセンサシステムにおいて、一定領域内にある目標位置を観測するセンサから目標の観測値を使って仮航跡を生成する追尾開始技術が必要となる。この発明は、この追尾開始において、センサのサーチ内で観測時刻にずれが生じる場合に、現実的な処理時間で目標航跡を追尾開始する多目標追尾装置に関するものである。   In a sensor system such as a search radar, a tracking start technique for generating a temporary track using a target observation value from a sensor that observes a target position within a certain area is required. The present invention relates to a multi-target tracking device that starts tracking a target track in a realistic processing time when there is a difference in observation time within a sensor search at the start of tracking.

センサにより一定領域内をサーチして得られた観測値を使って目標の追尾開始を行う技術については、すでに多くの論文、特許等の文献で取り挙げられており、その装置および方法については様々な提案がなされている。   There are many papers, patents, and other literature already on the technology to start tracking a target using observation values obtained by searching a certain area with a sensor. Various devices and methods are available. Proposals have been made.

目標追尾では、仮航跡の予測位置と観測位置を比較することによって仮航跡と観測値の対応付けが可能か否か(相関の有無)を判定する。通常の仮航跡の予測位置の計算では、仮航跡の最新時刻における平滑諸元を観測値の時刻へ外挿する計算を行う。捜索レーダによる観測の場合、図23の例の様に1サーチ内でもビームが到達する時刻によって観測値が得られる時刻は異なる。この相違のため、仮航跡の予測処理における予測時刻は観測値の数だけ発生する。   In target tracking, it is determined whether or not it is possible to associate the temporary track with the observed value by comparing the predicted position of the temporary track and the observed position (whether there is a correlation). In the calculation of the predicted position of a normal tentative track, a calculation is performed by extrapolating the smoothing parameters at the latest time of the tentative track to the time of the observed value. In the case of observation by a search radar, the time at which an observation value is obtained differs depending on the time when the beam arrives even within one search as in the example of FIG. Because of this difference, the predicted time in the temporary track prediction process is generated by the number of observation values.

この様な観測値を使って追尾開始処理を行う場合の従来の多目標追尾装置について図21から図24までを参照しながら説明する。図21は、従来の多目標追尾装置の構成を示すブロック図である。図22は、従来の多目標追尾装置の動作を示すフローチャートである。図23は、従来の多目標追尾装置におけるサーチ領域を示す図である。図24は、従来の多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。   A conventional multi-target tracking device in the case where tracking start processing is performed using such observation values will be described with reference to FIGS. FIG. 21 is a block diagram showing a configuration of a conventional multi-target tracking device. FIG. 22 is a flowchart showing the operation of the conventional multi-target tracking device. FIG. 23 is a diagram showing a search area in a conventional multi-target tracking device. FIG. 24 is a timing chart showing the relationship between the observation time of the conventional multi-target tracking device and the predicted time of the temporary track.

図21において、従来の多目標追尾装置20は、時刻、予測諸元を含む仮航跡を複数個格納する仮航跡ファイル21と、予測処理部22と、ゲート内外判定部23と、航跡尤度評価部24と、相関決定部25と、平滑処理部26と、航跡決定部27と、目標航跡を格納する目標航跡ファイル28とが設けられている。   In FIG. 21, a conventional multi-target tracking device 20 includes a temporary track file 21 that stores a plurality of temporary tracks including time and prediction specifications, a prediction processing unit 22, a gate inside / outside determination unit 23, and a track likelihood evaluation. A unit 24, a correlation determining unit 25, a smoothing processing unit 26, a track determining unit 27, and a target track file 28 for storing the target track are provided.

「観測値入力」ステップ900において、捜索レーダ(センサ)10から観測値が入力されると、まず、「予測処理」ステップ901において、予測処理部22は、仮航跡について観測値の時刻における予測位置と予測誤差共分散の計算を行う。   In the “observation value input” step 900, when an observation value is input from the search radar (sensor) 10, first, in a “prediction process” step 901, the prediction processing unit 22 predicts the predicted position at the time of the observation value for the temporary track. And calculate the prediction error covariance.

次に、「ゲート内外判定」ステップ902において、ゲート内外判定部23は、仮航跡と観測値が対応付けられる可能性を、仮航跡の予測誤差共分散と観測値の観測誤差共分散を使ったカイ平方検定によって判定する。   Next, in the “gate inside / outside determination” step 902, the gate inside / outside determination unit 23 uses the prediction error covariance of the temporary track and the observation error covariance of the observation value to determine the possibility that the temporary track and the observed value are associated with each other. Determined by chi-square test.

次に、「航跡尤度評価」ステップ903において、航跡尤度評価部24は、仮航跡と前ステップでゲート内にあると判定された観測値の組み合わせに対して、その組み合わせにより観測値を仮航跡に割り当てた場合の仮航跡の尤度を計算する。   Next, in the “wake likelihood evaluation” step 903, the wake likelihood evaluation unit 24 tentatively calculates an observation value based on the combination of the tentative track and the observation value determined to be in the gate in the previous step. Calculate the likelihood of a temporary wake when assigned to a wake.

次に、「相関決定」ステップ904において、相関決定部25は、観測値と仮航跡の組み合わせのどれを最終的に選択するか、前ステップで計算された尤度を基に決定する。   Next, in the “correlation determination” step 904, the correlation determination unit 25 determines which of the combination of the observed value and the provisional track is finally selected based on the likelihood calculated in the previous step.

次に、「平滑処理」ステップ905において、平滑処理部26は、前ステップで決定された相関結果を基に、仮航跡の予測諸元と組み合わされた観測値を使って平滑諸元を計算する。すなわち、仮航跡ファイル21に格納されている仮航跡の時刻を最新時刻の値に更新する。   Next, in the “smoothing process” step 905, the smoothing processing unit 26 calculates the smooth specification using the observation value combined with the prediction specification of the temporary track based on the correlation result determined in the previous step. . That is, the temporary track time stored in the temporary track file 21 is updated to the latest time value.

そして、「航跡決定」ステップ906において、航跡決定部27は、仮航跡の尤度や観測値との相関回数を基に、仮航跡を目標航跡と見做すかどうか、決定する。決定した目標航跡を目標航跡ファイル28に格納する。   Then, in a “track determination” step 906, the track determination unit 27 determines whether or not to consider the temporary track as the target track based on the likelihood of the temporary track and the number of correlations with the observed value. The determined target track is stored in the target track file 28.

また、捜索レーダを想定した1サーチ内での観測値の時刻ずれに起因する演算負荷増大の問題を解消するための従来の他の目標追尾装置が提案されている(例えば、特許文献1参照)。この目標追尾装置では、観測領域を事前に分割し、分割領域毎に追尾処理を行うことによって演算負荷の削減を図っている。   In addition, another conventional target tracking device for solving the problem of an increase in calculation load caused by a time lag of observation values in one search assuming a search radar has been proposed (see, for example, Patent Document 1). . In this target tracking device, the observation region is divided in advance, and the tracking processing is performed for each divided region to reduce the calculation load.

特開2004−233136号公報JP 2004-233136 A

上述したような従来の多目標追尾装置では、時刻ずれがある1サーチ分の観測値を処理しようとすると、1サーチの追尾開始処理において必要な予測処理回数は「仮航跡数×観測値数」となる。図23に示す観測例では、図24に示すタイミングチャートの様に計25回の予測処理が必要となる。予測処理に起因する演算負荷は、目標数が多い程高くなり、現実的な処理時間内での追尾開始が困難となる。   In the conventional multi-target tracking device as described above, if an attempt is made to process observation values for one search with a time lag, the number of prediction processes required in the tracking start processing of one search is “the number of temporary tracks × the number of observation values”. It becomes. In the observation example shown in FIG. 23, a total of 25 prediction processes are required as in the timing chart shown in FIG. The calculation load resulting from the prediction process increases as the target number increases, making it difficult to start tracking within a realistic processing time.

また、時刻ずれ観測値の処理時間短縮を図った従来の他の目標追尾装置では、目標航跡が分割領域の境界を跨って運動する場合に正確な追尾を行うことができないという問題点があった。   In addition, other conventional target tracking devices designed to reduce the processing time of the time lag observation values have had the problem that accurate tracking cannot be performed when the target track moves across the boundary of the divided areas. .

この発明は、上述のような課題を解決するためになされたもので、その目的は、許容可能な処理時間で、正確な追尾結果を得ることができる多目標追尾装置を得るものである。   The present invention has been made to solve the above-described problems, and an object thereof is to obtain a multi-target tracking device capable of obtaining an accurate tracking result in an allowable processing time.

この発明に係る多目標追尾装置は、センサで観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、予測計算を行う仮航跡と観測値の組み合わせを決定する観測値選別部と、前記観測値選別部により決定された仮航跡と観測値の組み合わせについて、仮航跡を観測値の時刻によって予測計算する予測処理部と、仮航跡と観測値が対応付け可能かを、前記予測処理部が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部と、観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部と、前記航跡尤度評価部が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部と、仮航跡の予測諸元と前記相関決定部により対応付けられた観測値を使って平滑諸元を計算する平滑処理部と、仮航跡のうち尤度が高いものを目標航跡とする航跡決定部とを設けたものである。   The multi-target tracking device according to the present invention is a multi-target tracking device that creates a target track using observation values obtained by searching an observation area with a sensor, and is a combination of a temporary track and an observation value for performing prediction calculation An observation value selection unit that determines the tentative value, a prediction processing unit that predicts and calculates the temporary track according to the time of the observation value for the combination of the temporary track and the observation value determined by the observation value selection unit, and the provisional track and the observation value A gate internal / external determination unit that determines whether the predicted value calculated by the prediction processing unit is compared with an observed value, and a track that calculates the likelihood of the temporary track when the observed value is associated A likelihood evaluation unit, a correlation determination unit that determines which temporary track is associated with which observation value, using the likelihood calculated by the track likelihood evaluation unit as an index, a prediction parameter of the temporary track, and the correlation determination unit Associated The observed values with the smoothing processing unit for calculating a smoothed specifications was one in which those likelihood of tentative track is high is provided a track determiner that the target track.

この発明に係る多目標追尾装置は、許容可能な処理時間で、正確な追尾結果を得ることができるという効果を奏する。   The multi-target tracking device according to the present invention has an effect that an accurate tracking result can be obtained in an allowable processing time.

実施の形態1.
この発明の実施の形態1に係る多目標追尾装置について図1から図4までを参照しながら説明する。図1は、この発明の実施の形態1に係る多目標追尾装置の構成を示すブロック図である。なお、以降では、各図中、同一符号は同一又は相当部分を示す。
Embodiment 1 FIG.
A multi-target tracking device according to Embodiment 1 of the present invention will be described with reference to FIGS. FIG. 1 is a block diagram showing the configuration of the multitarget tracking apparatus according to Embodiment 1 of the present invention. In the following, in each figure, the same reference numerals indicate the same or corresponding parts.

図1において、この実施の形態1に係る多目標追尾装置20は、観測値選別部1と、時刻、予測諸元を含む仮航跡を複数個格納する仮航跡ファイル21と、予測処理部22と、ゲート内外判定部23と、航跡尤度評価部24と、相関決定部25と、平滑処理部26と、航跡決定部27と、目標航跡を格納する目標航跡ファイル28とが設けられている。   In FIG. 1, the multi-target tracking device 20 according to the first embodiment includes an observation value selection unit 1, a temporary track file 21 that stores a plurality of temporary tracks including time and prediction specifications, and a prediction processing unit 22. A gate inside / outside determination unit 23, a track likelihood evaluation unit 24, a correlation determination unit 25, a smoothing processing unit 26, a track determination unit 27, and a target track file 28 for storing the target track are provided.

つぎに、この実施の形態1に係る多目標追尾装置の動作について図面を参照しながら説明する。図2は、この発明の実施の形態1に係る多目標追尾装置の動作を示すフローチャートである。図3は、この発明の実施の形態1に係る多目標追尾装置におけるサーチ領域を示す図である。図4は、この発明の実施の形態1に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。   Next, the operation of the multi-target tracking device according to the first embodiment will be described with reference to the drawings. FIG. 2 is a flowchart showing the operation of the multitarget tracking apparatus according to Embodiment 1 of the present invention. FIG. 3 is a diagram showing a search area in the multi-target tracking device according to Embodiment 1 of the present invention. FIG. 4 is a timing chart showing the relationship between the observation time of the multi-target tracking device according to Embodiment 1 of the present invention and the predicted time of the temporary track.

「観測値入力」ステップ100において、捜索レーダ(センサ)10から観測値が入力されると、まず「観測値選別処理」ステップ101において、観測値選別部1は、仮航跡に対して、予測処理を行う観測値を選別する。この選別では、仮航跡と観測値の全ての組み合わせについて、その組み合わせを相関の候補とするか否かを判定する。ある仮航跡とある観測値の組み合わせについての判定は、以下の4方式の何れかとする。   In the “observation value input” step 100, when an observation value is input from the search radar (sensor) 10, first, in the “observation value selection process” step 101, the observation value selection unit 1 performs a prediction process on the temporary track. Select observations to be performed. In this selection, it is determined whether or not all combinations of the temporary track and the observed value are candidates for correlation. The determination of a combination of a certain temporary track and a certain observed value is any one of the following four methods.

(1−1)観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差が、閾値(第1の閾値)を超えない場合にその観測値を候補とする。すなわち、以下の不等式(1)が成立するときに候補とする。   (1-1) If the difference between the observation time of the observation value and the smoothing time in the previous search of the temporary track does not exceed the threshold value (first threshold value), the observation value is set as a candidate. That is, it is determined as a candidate when the following inequality (1) holds.

Figure 2007292553
Figure 2007292553

ここで、ΔTは観測時刻と仮航跡の平滑時刻の差であり、Tsearchはサーチにかかる時間、thは事前に設定する閾値(第1の閾値)である。すなわち、観測値選別部1は、入力される観測値に含まれる観測時刻と、仮航跡ファイル21から読み出した仮航跡に含まれる平滑時刻とからΔTを求める。また、求めたΔTと、捜索レーダ(センサ)10から入手したTsearchから、観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差を求める。次に、求めた観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差が、閾値th未満の場合には、その観測値を選別する。 Here, ΔT is the difference between the observation time and the smooth time of the temporary track, T search is the time required for the search, and th T is a threshold value (first threshold value) set in advance. That is, the observation value selection unit 1 obtains ΔT from the observation time included in the input observation value and the smoothing time included in the temporary track read from the temporary track file 21. Further, from the obtained ΔT and T search obtained from the search radar (sensor) 10, the difference between the observation time of the observed value and the smooth time in the previous search of the temporary track is obtained. Next, when the difference between the observation time of the obtained observation value and the smoothing time in the previous search of the temporary track is less than the threshold th T , the observation value is selected.

(1−2)図3に示す様に、観測値の時刻と仮航跡の前サーチにおける平滑時刻の差が、サーチの角速度とサーチにかかる時間と目標の最大角速度から求めた許容時間差を越えない場合に、その観測値を候補とする。この場合の判定式(2)は以下のようになる。   (1-2) As shown in FIG. 3, the difference between the time of the observed value and the smoothing time in the search before the temporary track does not exceed the allowable time difference obtained from the angular velocity of the search, the time required for the search, and the target maximum angular velocity. In this case, the observed value is a candidate. The judgment formula (2) in this case is as follows.

Figure 2007292553
Figure 2007292553

ここで、ΔTは観測時刻と仮航跡の平滑時刻の差であり、Tsearchはサーチにかかる時間、vAz_maxは事前に設定する目標の最大角速度、vsearchはサーチの角速度である。すなわち、観測値選別部1は、入力される観測値に含まれる観測時刻と、仮航跡ファイル21から読み出した仮航跡に含まれる平滑時刻とからΔTを求める。また、求めたΔTと、捜索レーダ(センサ)10から入手したTsearchから、観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差を求める。さらに、目標の最大角速度と、捜索レーダ(センサ)10から入手したTsearch、vsearchから、式(2)の右辺の値を求める。次に、求めた観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差が、式(2)の右辺の値未満の場合には、その観測値を選別する。 Here, ΔT is the difference between the observation time and the smooth time of the tentative track, T search is the time required for the search , v Az_max is the target maximum angular velocity set in advance, and v search is the angular velocity of the search . That is, the observation value selection unit 1 obtains ΔT from the observation time included in the input observation value and the smoothing time included in the temporary track read from the temporary track file 21. Further, from the obtained ΔT and T search obtained from the search radar (sensor) 10, the difference between the observation time of the observed value and the smooth time in the previous search of the temporary track is obtained. Further, the value on the right side of the equation (2) is obtained from the target maximum angular velocity and T search and v search obtained from the search radar (sensor) 10. Next, if the difference between the observed time of the observed value and the smoothing time in the previous search of the temporary track is less than the value on the right side of Equation (2), the observed value is selected.

(1−3)上記の(1−2)の判定に加え、さらに観測値と仮航跡の距離差も判定に加える。すなわち、式(2)の条件に加えて、以下の式(3)の条件も成立する場合に限って、その観測値を候補とする。   (1-3) In addition to the determination of (1-2) above, the distance difference between the observed value and the temporary track is also added to the determination. That is, the observation value is a candidate only when the following equation (3) is satisfied in addition to the equation (2).

Figure 2007292553
Figure 2007292553

ここで、Rは仮航跡の前サーチの平滑時刻からTsearch(1サーチにかかる時間)後の位置から計算したセンサからの距離、Rは観測値のセンサからの距離であり、thは事前に設定する閾値(第2の閾値)である。すなわち、観測値選別部1は、仮航跡ファイル21から読み出した仮航跡に含まれる前サーチの平滑時刻と、捜索レーダ(センサ)10から入手したTsearchとからRを求める。また、求めたRと、入力される観測値に含まれるRをから、観測値と仮航跡の距離差を求める。次に、求めた観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差が、式(2)の右辺の値未満の場合で、かつ、観測値と仮航跡の距離差が、閾値thR未満の場合には、その観測値を選別する。 Here, R T is the distance from the sensor calculated from the position after T search (time required for one search) from the smoothing time of the previous search of the tentative track, R O is the distance from the sensor of the observed value, and th R Is a threshold (second threshold) set in advance. That is, the observed value selecting unit 1 obtains a smoothing time of the search before included in the tentative track read from the tentative track file 21, the R T and a T search obtained from search radar (sensor) 10. Further, from the R T obtained, the R O contained in the observation value input, obtains the distance difference observations and tentative track. Next, when the difference between the observed time of the observed value and the smoothing time in the previous search for the temporary track is less than the value on the right side of Equation (2), and the distance difference between the observed value and the temporary track is the threshold th If it is less than R , the observed value is selected.

(1−4)上記の判定式(2)において、vAz_maxを事前設定のパラメータでなく、前サーチで計算された仮航跡の速度の中の最大の値とする。すなわち、観測値選別部1は、仮航跡ファイル21から読み出した仮航跡に含まれる速度の中の最大値を求める。また、入力される観測値に含まれる観測時刻と、仮航跡ファイル21から読み出した仮航跡に含まれる平滑時刻とからΔTを求める。また、求めたΔTと、捜索レーダ(センサ)10から入手したTsearchから、観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差を求める。さらに、求めた仮航跡の速度の中の最大の値と、捜索レーダ(センサ)10から入手したTsearch、vsearchから、式(2)の右辺の値を求める。次に、求めた観測値の観測時刻と仮航跡の前サーチにおける平滑時刻の差が、式(2)の右辺の値未満の場合には、その観測値を選別する。 (1-4) In the above-described determination formula (2), v Az_max is not the preset parameter, but the maximum value of the speed of the temporary track calculated in the previous search. That is, the observation value selection unit 1 obtains the maximum value among the speeds included in the temporary track read from the temporary track file 21. Further, ΔT is obtained from the observation time included in the input observation value and the smoothing time included in the temporary track read from the temporary track file 21. Further, from the obtained ΔT and T search obtained from the search radar (sensor) 10, the difference between the observation time of the observed value and the smooth time in the previous search of the temporary track is obtained. Further, the value on the right side of the equation (2) is obtained from the maximum value of the obtained tentative track speed and the T search and v search obtained from the search radar (sensor) 10. Next, if the difference between the observed time of the observed value and the smoothing time in the previous search of the temporary track is less than the value on the right side of Equation (2), the observed value is selected.

次に、「予測処理」ステップ102において、予測処理部22は、仮航跡について観測値の観測時刻における予測位置と予測誤差共分散の計算を行う。ただし、ある仮航跡と、その仮航跡について前ステップで候補と判定されなかった観測値に対する予測計算は行わない。   Next, in the “prediction process” step 102, the prediction processing unit 22 calculates the predicted position and the prediction error covariance at the observation time of the observed value for the temporary track. However, prediction calculation is not performed for a certain tentative track and an observed value that is not determined as a candidate in the previous step for the tentative track.

次に、「ゲート内外判定」ステップ103において、ゲート内外判定部23は、仮航跡と観測値が対応付けられる可能性を、仮航跡の予測誤差共分散と観測値の観測誤差共分散を使ったカイ平方検定によって判定する。   Next, in the “gate inside / outside determination” step 103, the gate inside / outside determination unit 23 uses the prediction error covariance of the temporary track and the observation error covariance of the observation value as a possibility of associating the temporary track with the observation value. Determined by chi-square test.

次に、「航跡尤度評価」ステップ104において、航跡尤度評価部24は、仮航跡と前ステップでゲート内にあると判定された観測値の組み合わせに対して、その組み合わせにより観測値を仮航跡に割り当てた場合の仮航跡の尤度を計算する。   Next, in the “wake likelihood evaluation” step 104, the wake likelihood evaluation unit 24 tentatively determines an observation value based on the combination of the temporary wake and the observation value determined to be in the gate in the previous step. Calculate the likelihood of a temporary wake when assigned to a wake.

次に、「相関決定」ステップ105において、相関決定部25は、観測値と仮航跡の組み合わせのどれを最終的に選択するか、前ステップで計算された尤度を基に決定する。   Next, in a “correlation determination” step 105, the correlation determination unit 25 determines which of the combination of the observed value and the provisional track is finally selected based on the likelihood calculated in the previous step.

次に、「平滑処理」ステップ106において、平滑処理部26は、前ステップで決定された相関結果を基に、仮航跡の予測諸元と組み合わされた観測値を使って平滑諸元を計算する。すなわち、仮航跡ファイル21に格納されている仮航跡の時刻を最新時刻の値に更新する。   Next, in the “smoothing process” step 106, the smoothing processing unit 26 calculates the smooth specification using the observation value combined with the prediction specification of the temporary track based on the correlation result determined in the previous step. . That is, the temporary track time stored in the temporary track file 21 is updated to the latest time value.

そして、「航跡決定」ステップ107において、航跡決定部27は、仮航跡の尤度や観測値との相関回数を基に、仮航跡を目標航跡と見做すかどうか、決定する。決定した目標航跡を目標航跡ファイル28に格納する。   In “track determination” step 107, the track determination unit 27 determines whether or not to consider the temporary track as the target track based on the likelihood of the temporary track and the number of correlations with the observed value. The determined target track is stored in the target track file 28.

本実施の形態1に係る多目標追尾装置を図23に示す観測例に適用した場合(図上、左から1番目〜5番目の仮航跡、1番目〜5番目の観測値とする。)、1番目の仮航跡について1番目〜3番目の観測値が上記の判定式(2)を満たし、2番目の仮航跡について2番目〜4番目の観測値が上記の判定式(2)を満たし、3番目の仮航跡について3番目〜4番目の観測値が上記の判定式(2)を満たし、4番目の仮航跡について4番目〜5番目の観測値が上記の判定式(2)を満たし、5番目の仮航跡について3番目〜4番目の観測値が上記の判定式(2)を満たしていると、図4のタイムチャートに示す様に計12回の予測計算が必要となる。これは従来方式で必要な予測計算回数25回よりも少ない。   When the multi-target tracking device according to the first embodiment is applied to the observation example shown in FIG. 23 (the first to fifth temporary tracks and the first to fifth observation values from the left in the figure). The first to third observation values for the first temporary track satisfy the above-described determination formula (2), and the second to fourth observation values for the second temporary track satisfy the above-described determination formula (2). The third to fourth observation values for the third temporary track satisfy the above-described determination formula (2), and the fourth to fifth observation values for the fourth temporary track satisfy the above-described determination formula (2). If the third to fourth observation values for the fifth temporary track satisfy the above-described determination formula (2), a total of twelve prediction calculations are required as shown in the time chart of FIG. This is less than the required number of prediction calculations of 25 in the conventional method.

すなわち、この実施の形態1に係る多目標追尾装置は、センサ10で観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、予測計算を行う仮航跡と観測値の組み合わせを決定する観測値選別部1と、前記観測値選別部1により決定された仮航跡と観測値の組み合わせについて、仮航跡を観測値の時刻によって予測計算する予測処理部22と、仮航跡と観測値が対応付け可能かを、前記予測処理部22が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部23と、観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部24と、前記航跡尤度評価部24が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部25と、仮航跡の予測諸元と前記相関決定部25により対応付けられた観測値を使って平滑諸元を計算する平滑処理部26と、仮航跡のうち尤度が高いものを目標航跡とする航跡決定部27とを設けたものである。   That is, the multi-target tracking device according to the first embodiment is a multi-target tracking device that creates a target wake using observation values obtained by searching an observation region with the sensor 10, and is a temporary target that performs prediction calculation. An observation value selection unit 1 that determines a combination of a wake and an observation value, and a prediction processing unit 22 that predicts and calculates the temporary track according to the time of the observation value for the combination of the tentative track and the observation value determined by the observation value selection unit 1. And the gate internal / external determination unit 23 that determines whether the temporary track and the observed value can be associated with each other by comparing the predicted value of the temporary track calculated by the prediction processing unit 22 with the observed value. The track likelihood evaluation unit 24 for calculating the likelihood of the temporary track in the case, and the correlation determination unit 25 for determining which temporary track and which observation value are associated with each other using the likelihood calculated by the track likelihood evaluation unit 24 as an index. When A smoothing processing unit 26 that calculates a smooth specification using the predicted value of the provisional wake and the observation value associated with the correlation determination unit 25, and a wake determination that uses a tentative wake with a high likelihood as a target wake. A portion 27 is provided.

以上の様に、この実施の形態1に係る多目標追尾装置によれば、各仮航跡に対して予測処理を行う観測値を選別する前処理を行うので、処理負荷が軽減される。よって、目標数が多い場合でも対処可能となる。   As described above, according to the multi-target tracking device according to the first embodiment, the preprocessing for selecting the observation value for which the prediction process is performed for each temporary track is performed, so that the processing load is reduced. Therefore, even when the target number is large, it can be dealt with.

実施の形態2.
この発明の実施の形態2に係る多目標追尾装置について図5から図9までを参照しながら説明する。図5は、この発明の実施の形態2に係る多目標追尾装置の構成を示すブロック図である。
Embodiment 2. FIG.
A multi-target tracking device according to Embodiment 2 of the present invention will be described with reference to FIGS. FIG. 5 is a block diagram showing the configuration of the multi-target tracking device according to Embodiment 2 of the present invention.

図5において、この実施の形態2に係る多目標追尾装置20は、観測値グループ分け部2と、グループ重心計算部3と、時刻、予測諸元を含む仮航跡を複数個格納する仮航跡ファイル21と、予測処理部22と、ゲート内外判定部23と、航跡尤度評価部24と、相関決定部25と、平滑処理部26と、航跡決定部27と、目標航跡を格納する目標航跡ファイル28とが設けられている。   In FIG. 5, the multi-target tracking device 20 according to the second embodiment includes an observation value grouping unit 2, a group center of gravity calculation unit 3, and a temporary track file that stores a plurality of temporary tracks including time and prediction specifications. 21, prediction processing unit 22, gate inside / outside determination unit 23, track likelihood evaluation unit 24, correlation determination unit 25, smoothing processing unit 26, track determination unit 27, and target track file that stores the target track 28 are provided.

つぎに、この実施の形態2に係る多目標追尾装置の動作について図面を参照しながら説明する。図6は、この発明の実施の形態2に係る多目標追尾装置の動作を示すフローチャートである。図7は、この発明の実施の形態2に係る多目標追尾装置におけるサーチ領域を示す図である。図8は、この発明の実施の形態2に係る多目標追尾装置の観測値グループ分け部の同一グループ判定を説明するための図である。図9は、この発明の実施の形態2に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。   Next, the operation of the multi-target tracking device according to the second embodiment will be described with reference to the drawings. FIG. 6 is a flowchart showing the operation of the multitarget tracking apparatus according to Embodiment 2 of the present invention. FIG. 7 is a diagram showing a search area in the multi-target tracking device according to Embodiment 2 of the present invention. FIG. 8 is a diagram for explaining the same group determination of the observation value grouping unit of the multitarget tracking apparatus according to Embodiment 2 of the present invention. FIG. 9 is a timing chart showing the relationship between the observation time of the multi-target tracking device according to Embodiment 2 of the present invention and the predicted time of the temporary track.

「観測値入力」ステップ200において、捜索レーダ(センサ)10から観測値が入力されると、まず、「観測値グループ分け」ステップ201において、観測値グループ分け部2は、距離の近い観測値同士でグループ化する。ある観測値と別の観測値の組み合わせについての同一グループ判定は以下の3方式の何れかとする。   In the “observation value input” step 200, when an observation value is input from the search radar (sensor) 10, first, in the “observation value grouping” step 201, the observation value grouping unit 2 determines that observation values having a short distance are Group by. The same group determination for a combination of a certain observation value and another observation value is any one of the following three methods.

(2−1)2つの観測値の観測時刻の差が、閾値(第3の閾値)を超えない場合にその観測値を候補とする。すなわち、以下の不等式(4)が成立するときに候補とする。   (2-1) If an observation time difference between two observation values does not exceed a threshold value (third threshold value), the observation value is set as a candidate. That is, it is determined as a candidate when the following inequality (4) holds.

Figure 2007292553
Figure 2007292553

ここで、ΔTは2つの観測値の観測時刻の差であり、thは事前に設定する閾値(第3の閾値)である。すなわち、観測値グループ分け部2は、入力される観測値に含まれる観測時刻から、2つの観測値の観測時刻の差を求める。次に、求めた2つの観測値の観測時刻の差が、閾値th未満の場合には、これらの2つの観測値を同一グループと判定する。 Here, ΔT is a difference between observation times of two observation values, and th T is a threshold value (third threshold value) set in advance. That is, the observation value grouping unit 2 obtains the difference between the observation times of the two observation values from the observation times included in the input observation values. Next, when the difference between the observed times of the two observed values is less than the threshold th T , these two observed values are determined as the same group.

(2−2)図8に示す様に、2つの観測値の観測時刻の差が、サーチの角速度とサーチにかかる時間と目標の最大角速度から求めた許容時間差を越えない場合に、その観測値を候補とする。この場合の判定式(5)は以下のようになる。   (2-2) As shown in FIG. 8, when the difference between the observation times of the two observation values does not exceed the allowable time difference obtained from the angular velocity of the search, the time required for the search, and the target maximum angular velocity, the observation value Is a candidate. The judgment formula (5) in this case is as follows.

Figure 2007292553
Figure 2007292553

ここで、ΔTは2つの観測値の時刻の差であり、Tsearchはサーチにかかる時間、vAz_maxは事前に設定する目標の最大角速度、vsearchはサーチの角速度である。すなわち、観測値グループ分け部2は、入力される観測値に含まれる観測時刻から、2つの観測値の観測時刻の差を求める。また、目標の最大角速度と、捜索レーダ(センサ)10から入手したTsearch、vsearchから、式(5)の右辺を求める。次に、求めた2つの観測値の観測時刻の差が、式(5)の右辺の値未満の場合には、これらの2つの観測値を同一グループと判定する。 Here, ΔT is the time difference between the two observation values, T search is the time required for the search , v Az_max is the target maximum angular velocity set in advance, and v search is the angular velocity of the search . That is, the observation value grouping unit 2 obtains the difference between the observation times of the two observation values from the observation times included in the input observation values. Further, the right side of Expression (5) is obtained from the maximum angular velocity of the target and T search and v search obtained from the search radar (sensor) 10. Next, when the difference between the observation times of the two obtained observation values is less than the value on the right side of Equation (5), these two observation values are determined to be the same group.

(2−3)上記の(2−2)の判定に加え、さらに、2つの観測値の距離差も判定に加える。すなわち、式(5)の条件に加えて、以下の式(6)の条件も成立する場合に限って、その観測値を候補とする。   (2-3) In addition to the determination of (2-2) above, the distance difference between the two observation values is also added to the determination. That is, the observation value is a candidate only when the following equation (6) is satisfied in addition to the equation (5).

Figure 2007292553
Figure 2007292553

ここで、ΔRは2つの観測値の距離差、thは事前に設定する閾値(第4の閾値)である。すなわち、観測値グループ分け部2は、入力される観測値に含まれる観測時刻から、2つの観測値の観測時刻の差を求める。また、目標の最大角速度と、捜索レーダ(センサ)10から入手したTsearch、vsearchから、式(5)の右辺を求める。さらに、入力される観測値に含まれる距離から、2つの観測値の距離差を求める。次に、求めた2つの観測値の観測時刻の差が、式(5)の右辺の値未満の場合で、かつ、求めた2つの観測値の距離差が、閾値th未満の場合には、これらの2つの観測値を同一グループと判定する。 Here, ΔR is a distance difference between two observation values, and th R is a threshold (fourth threshold) set in advance. That is, the observation value grouping unit 2 obtains the difference between the observation times of the two observation values from the observation times included in the input observation values. Further, the right side of Expression (5) is obtained from the maximum angular velocity of the target and T search and v search obtained from the search radar (sensor) 10. Further, the distance difference between the two observation values is obtained from the distance included in the input observation value. Next, when the difference in observation time between the two observed values is less than the value on the right side of Equation (5) and the distance difference between the two observed values is less than the threshold th R These two observation values are determined as the same group.

次に、「グループ重心計算」ステップ202において、グループ重心計算部3は、グループに所属する観測値群の重心の位置に相当する観測時刻を計算する(以下、重心時刻と呼ぶ)。グループ重心計算部3は、入力される観測値に含まれる位置、観測時刻から、平均などにより、同じグループに所属する観測値群の重心の位置を求め、重心の位置に対応する観測時刻、つまり重心時刻を求める。   Next, in the “group centroid calculation” step 202, the group centroid calculation unit 3 calculates the observation time corresponding to the position of the centroid of the observation value group belonging to the group (hereinafter referred to as the centroid time). The group center-of-gravity calculation unit 3 obtains the position of the center of gravity of the observation value group belonging to the same group from the position and the observation time included in the input observation value by averaging, and the observation time corresponding to the position of the center of gravity, that is, Find the centroid time.

次に、「予測処理」ステップ203において、予測処理部22は、仮航跡について観測値が所属するグループの重心時刻における予測位置と予測誤差共分散の計算を行う。   Next, in the “prediction process” step 203, the prediction processing unit 22 calculates the predicted position and the prediction error covariance at the barycentric time of the group to which the observation value belongs for the temporary track.

次に、「ゲート内外判定」ステップ204において、ゲート内外判定部23は、仮航跡と観測値が対応付けられる可能性を、仮航跡の予測誤差共分散と観測値の観測誤差共分散を使ったカイ平方検定によって判定する。   Next, in the “gate inside / outside determination” step 204, the gate inside / outside determination unit 23 uses the prediction error covariance of the temporary track and the observation error covariance of the observation value as a possibility of associating the temporary track with the observation value. Determined by chi-square test.

次に、「航跡尤度評価」ステップ205において、航跡尤度評価部24は、仮航跡と前ステップでゲート内にあると判定された観測値の組み合わせに対して、その組み合わせにより観測値を仮航跡に割り当てた場合の仮航跡の尤度を計算する。   Next, in the “track likelihood evaluation” step 205, the track likelihood evaluation unit 24, based on the combination of the temporary track and the observation value determined to be in the gate in the previous step, temporarily calculates the observation value. Calculate the likelihood of a temporary wake when assigned to a wake.

次に、「相関決定」ステップ206において、相関決定部25は、観測値と仮航跡の組み合わせのどれを最終的に選択するか、前ステップで計算された尤度を基に決定する。   Next, in a “correlation determination” step 206, the correlation determination unit 25 determines which of the combination of the observed value and the provisional track is finally selected based on the likelihood calculated in the previous step.

次に、「平滑処理」ステップ207において、平滑処理部26は、前ステップで決定された相関結果を基に、仮航跡の予測諸元と組み合わされた観測値を使って平滑諸元を計算する。すなわち、仮航跡ファイル21に格納されている仮航跡の時刻を最新時刻の値に更新する。   Next, in the “smoothing process” step 207, the smoothing processing unit 26 calculates the smooth specification using the observation value combined with the prediction specification of the temporary track based on the correlation result determined in the previous step. . That is, the temporary track time stored in the temporary track file 21 is updated to the latest time value.

そして、「航跡決定」ステップ208において、航跡決定部27は、仮航跡の尤度や観測値との相関回数を基に、仮航跡を目標航跡と見做すかどうか、決定する。決定した目標航跡を目標航跡ファイル28に格納する。   Then, in a “track determination” step 208, the track determination unit 27 determines whether to consider the temporary track as the target track based on the likelihood of the temporary track and the number of correlations with the observed value. The determined target track is stored in the target track file 28.

本実施の形態2に係る多目標追尾装置を図23に示す観測例に適用した場合、図7の様にグループが2つできるとすると、図9のタイムチャートに示す様に、仮航跡数×グループ数(5×2)の計10回の予測計算が必要となる。これは従来方式で必要な予測計算回数25回よりも少ない。   When the multi-target tracking device according to the second embodiment is applied to the observation example shown in FIG. 23, assuming that two groups can be formed as shown in FIG. 7, the number of provisional tracks ×× as shown in the time chart of FIG. A total of 10 prediction calculations for the number of groups (5 × 2) are required. This is less than the required number of prediction calculations of 25 in the conventional method.

すなわち、この実施の形態2に係る多目標追尾装置は、センサ10で観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、1サーチで得られた観測値群をグループ分けする観測値グループ分け部2と、グループ毎にグループの重心位置に相当する観測時刻である重心時刻を計算するグループ重心計算部3と、仮航跡と観測値の組み合わせについて、仮航跡を観測値が所属するグループの前記重心時刻によって予測計算する予測処理部22と、仮航跡と観測値が対応付け可能かを、前記予測処理部22が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部23と、観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部24と、前記航跡尤度評価部24が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部25と、仮航跡の予測諸元と前記相関決定部25により対応付けられた観測値を使って平滑諸元を計算する平滑処理部26と、仮航跡のうち尤度が高いものを目標航跡とする航跡決定部27とを設けたものである。   That is, the multi-target tracking device according to the second embodiment is a multi-target tracking device that creates a target wake using observation values obtained by searching an observation region with the sensor 10, and is obtained in one search. The observed value grouping unit 2 for grouping the observed value groups, the group centroid calculating unit 3 for calculating the centroid time which is the observation time corresponding to the centroid position of the group for each group, and the combination of the temporary track and the observed value The prediction processing unit 22 that predicts and calculates the temporary track based on the centroid time of the group to which the observation value belongs, and the prediction value of the temporary track that the prediction processing unit 22 calculates whether the temporary track and the observation value can be associated with each other. Gate inside / outside determination unit 23 that is determined by comparing observation values, track likelihood evaluation unit 24 that calculates the likelihood of a temporary track when the observation values are associated, and the track likelihood evaluation unit 2 The correlation determination unit 25 that determines which temporary track and which observation value are associated with each other using the likelihood calculated by the above as an index, and smoothing using the predicted values of the temporary track and the observation values associated with the correlation determination unit 25 A smoothing processing unit 26 that calculates specifications and a track determination unit 27 that uses a temporary track having a high likelihood as a target track are provided.

以上の様に、この実施の形態2に係る多目標追尾装置によれば、観測値をグループ化してその重心時刻で各仮航跡に対する予測処理を行うので、処理負荷が軽減される。よって、目標数が多い場合でも対処可能となる。   As described above, according to the multi-target tracking device according to the second embodiment, the observation values are grouped and the prediction process for each temporary track is performed at the center of gravity time, so the processing load is reduced. Therefore, even when the target number is large, it can be dealt with.

実施の形態3.
この発明の実施の形態3に係る多目標追尾装置について図10から図13までを参照しながら説明する。図10は、この発明の実施の形態3に係る多目標追尾装置の構成を示すブロック図である。
Embodiment 3 FIG.
A multi-target tracking device according to Embodiment 3 of the present invention will be described with reference to FIGS. FIG. 10 is a block diagram showing the configuration of the multitarget tracking apparatus according to Embodiment 3 of the present invention.

図10において、この実施の形態3に係る多目標追尾装置20は、重心計算部4と、誤差共分散拡大部5と、時刻、予測諸元を含む仮航跡を複数個格納する仮航跡ファイル21と、予測処理部22と、ゲート内外判定部23と、航跡尤度評価部24と、相関決定部25と、平滑処理部26と、航跡決定部27と、目標航跡を格納する目標航跡ファイル28とが設けられている。   In FIG. 10, the multi-target tracking device 20 according to the third embodiment includes a centroid calculating unit 4, an error covariance expanding unit 5, and a temporary track file 21 that stores a plurality of temporary tracks including time and prediction specifications. A prediction processing unit 22, a gate inside / outside determination unit 23, a track likelihood evaluation unit 24, a correlation determination unit 25, a smoothing processing unit 26, a track determination unit 27, and a target track file 28 that stores a target track. And are provided.

つぎに、この実施の形態3に係る多目標追尾装置の動作について図面を参照しながら説明する。図11は、この発明の実施の形態3に係る多目標追尾装置の動作を示すフローチャートである。図12は、この発明の実施の形態3に係る多目標追尾装置におけるサーチ領域を示す図である。図13は、この発明の実施の形態3に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。   Next, the operation of the multi-target tracking device according to the third embodiment will be described with reference to the drawings. FIG. 11 is a flowchart showing the operation of the multitarget tracking apparatus according to Embodiment 3 of the present invention. FIG. 12 is a diagram showing search areas in the multitarget tracking apparatus according to Embodiment 3 of the present invention. FIG. 13 is a timing chart showing the relationship between the observation time of the multi-target tracking device according to Embodiment 3 of the present invention and the predicted time of the temporary track.

「観測値入力」ステップ300において、捜索レーダ(センサ)10から観測値が入力されると、まず、「重心計算」ステップ301において、重心計算部4は、得られた全観測値の重心の位置に相当する観測時刻を計算する(以下、全体重心時刻と呼ぶ)。重心計算部4は、入力される観測値に含まれる位置、観測時刻から、平均などにより、全観測値の重心の位置を求め、重心の位置に対応する観測時刻、つまり全体重心時刻を求める。   When an observation value is input from the search radar (sensor) 10 in the “observation value input” step 300, first, in a “centroid calculation” step 301, the centroid calculation unit 4 determines the position of the centroid of all the obtained observation values. Is calculated (hereinafter referred to as the total centroid time). The center-of-gravity calculation unit 4 obtains the position of the center of gravity of all the observed values from the position and the observation time included in the input observation value by averaging, and obtains the observation time corresponding to the position of the center of gravity, that is, the total center-of-gravity time.

次に、「予測処理」ステップ302において、予測処理部22は、仮航跡について全体重心時刻における予測位置と予測誤差共分散の計算を行う。   Next, in the “prediction process” step 302, the prediction processing unit 22 calculates the predicted position and the prediction error covariance at the total center-of-gravity time for the temporary track.

次に、「誤差共分散拡大」ステップ303において、誤差共分散拡大部5は、全体重心時刻で予測処理したことによる誤差分を考慮して、予測誤差共分散行列を拡大する。この拡大分の計算は以下のいずれかとする。   Next, in an “error covariance expansion” step 303, the error covariance expansion unit 5 expands the prediction error covariance matrix in consideration of the error due to the prediction processing performed at the total centroid time. The calculation for this enlargement is either:

(3−1)事前に設定されたパラメータである目標の最大角速度をvmax、全体重心時刻と各観測値の観測時刻の差の最大値をΔtmaxとすると、予測誤差共分散の拡大分をvmaxΔtmaxとする。すなわち、誤差共分散拡大部5は、入力される観測値に含まれる観測時刻から、ステップ301で求めた全体重心時刻と各観測値の観測時刻の差の最大値を求める。次に、目標の最大角速度vmaxと最大値Δtmaxの積から予測誤差共分散の拡大分を求める。 (3-1) Assuming that the target maximum angular velocity, which is a preset parameter, is v max , and the maximum value of the difference between the overall centroid time and the observation time of each observation value is Δt max , the expansion of the prediction error covariance is Let v max Δt max . That is, the error covariance enlarging unit 5 obtains the maximum value of the difference between the total centroid time obtained in step 301 and the observation time of each observation value from the observation time included in the input observation value. Next, an enlargement of the prediction error covariance is obtained from the product of the target maximum angular velocity v max and the maximum value Δt max .

(3−2)仮航跡の速度の最大値をvTmax、全体重心時刻と各観測値の観測時刻の差の最大値をΔtmaxとすると、予測誤差共分散の拡大分をvTmaxΔtmaxとする。すなわち、誤差共分散拡大部5は、仮航跡ファイル21から読み出した仮航跡に含まれる速度から、仮航跡の速度の最大値を求める。次に、仮航跡の速度の最大値vTmaxと最大値Δtmaxの積から予測誤差共分散の拡大分を求める。 (3-2) Assuming that the maximum value of the speed of the temporary track is v Tmax and the maximum value of the difference between the overall center of gravity time and the observation time of each observation value is Δt max , the expanded portion of the prediction error covariance is v Tmax Δt max To do. That is, the error covariance enlarging unit 5 obtains the maximum value of the speed of the temporary track from the speed included in the temporary track read from the temporary track file 21. Next, an enlargement of the prediction error covariance is obtained from the product of the maximum value v Tmax and the maximum value Δt max of the speed of the temporary track.

次に、「ゲート内外判定」ステップ304において、ゲート内外判定部23は、仮航跡と観測値が対応付けられる可能性を、仮航跡の予測誤差共分散と観測値の観測誤差共分散を使ったカイ平方検定によって判定する。   Next, in the “gate inside / outside determination” step 304, the gate inside / outside determination unit 23 uses the prediction error covariance of the temporary track and the observation error covariance of the observation value to determine whether the temporary track and the observed value are associated with each other. Determined by chi-square test.

次に、「航跡尤度評価」ステップ305において、航跡尤度評価部24は、仮航跡と前ステップでゲート内にあると判定された観測値の組み合わせに対して、その組み合わせにより観測値を仮航跡に割り当てた場合の仮航跡の尤度を計算する。   Next, in “Wake Likelihood Evaluation” step 305, the wake likelihood evaluation unit 24 tentatively determines an observation value based on the combination of the tentative track and the observation value determined to be in the gate in the previous step. Calculate the likelihood of a temporary wake when assigned to a wake.

次に、「相関決定」ステップ306において、相関決定部25は、観測値と仮航跡の組み合わせのどれを最終的に選択するか、前ステップで計算された尤度を基に決定する。   Next, in the “correlation determination” step 306, the correlation determination unit 25 determines which of the combination of the observed value and the provisional track is to be finally selected based on the likelihood calculated in the previous step.

次に、「平滑処理」ステップ307において、平滑処理部26は、前ステップで決定された相関結果を基に、仮航跡の予測諸元と組み合わされた観測値を使って平滑諸元を計算する。すなわち、仮航跡ファイル21に格納されている仮航跡の時刻を最新時刻の値に更新する。   Next, in the “smoothing process” step 307, the smoothing processing unit 26 calculates the smooth specification using the observation value combined with the prediction specification of the temporary track based on the correlation result determined in the previous step. . That is, the temporary track time stored in the temporary track file 21 is updated to the latest time value.

そして、「航跡決定」ステップ308において、航跡決定部27は、仮航跡の尤度や観測値との相関回数を基に、仮航跡を目標航跡と見做すかどうか、決定する。決定した目標航跡を目標航跡ファイル28に格納する。   In a “wake determination” step 308, the wake determination unit 27 determines whether the temporary wake is regarded as the target wake based on the likelihood of the tentative wake and the number of correlations with the observed value. The determined target track is stored in the target track file 28.

本実施の形態3に係る多目標追尾装置を図23に示す観測例に適用した場合、図12の様に誤差楕円が拡大され、図13のタイムチャートに示す様に仮航跡数の計5回の予測計算が必要となる。これは従来方式で必要な予測計算回数25回よりも少ない。   When the multi-target tracking device according to the third embodiment is applied to the observation example shown in FIG. 23, the error ellipse is enlarged as shown in FIG. 12, and the total number of tentative tracks is five times as shown in the time chart of FIG. Prediction calculation is required. This is less than the required number of prediction calculations of 25 in the conventional method.

すなわち、この実施の形態3に係る多目標追尾装置は、センサ10で観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、1サーチで得られた全観測値の重心の位置に相当する観測時刻である全体重心時刻を計算する重心計算部4と、仮航跡と観測値の組み合わせについて、仮航跡を前記全体重心時刻によって予測計算する予測処理部22と、前記全体重心時刻で予測処理したことによる誤差分を考慮して、予測誤差共分散行列を拡大する誤差共分散拡大部5と、仮航跡と観測値が対応付け可能かを、前記予測処理部22が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部23と、観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部24と、前記航跡尤度評価部24が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部25と、仮航跡の予測諸元と前記相関決定部25により対応付けられた観測値を使って平滑諸元を計算する平滑処理部26と、仮航跡のうち尤度が高いものを目標航跡とする航跡決定部27とを設けたものである。   That is, the multi-target tracking device according to the third embodiment is a multi-target tracking device that creates a target wake using observation values obtained by searching an observation region with the sensor 10, and is obtained in one search. A center-of-gravity calculation unit 4 for calculating the total center-of-gravity time, which is an observation time corresponding to the position of the center of gravity of all observed values, and a prediction processing unit for predicting and calculating the temporary track based on the total center-of-gravity time for a combination of the temporary track and the observed value 22 and the error covariance expansion unit 5 that expands the prediction error covariance matrix in consideration of the error due to the prediction processing at the total center-of-gravity time, and whether the temporary track and the observation value can be associated with each other A gate inside / outside determination unit 23 that determines by comparing the predicted value of the temporary track calculated by the processing unit 22 and the observed value, and a track likelihood evaluation unit 24 that calculates the likelihood of the temporary track when the observed value is associated. When, Correlation determining unit 25 that determines which temporary track and which observation value are associated with each other using the likelihood calculated by recorded track likelihood evaluation unit 24 as an index, and the prediction parameters of temporary track and the correlation determining unit 25 are associated with each other. The smoothing processing unit 26 that calculates the smooth specification using the observed values, and the wake determination unit 27 that uses the tentative track having a high likelihood as the target track is provided.

以上の様に、この実施の形態3に係る多目標追尾装置によれば、観測値全体の全体重心時刻で各仮航跡に対する予測処理を行い、さらに時刻近似によって発生する共分散行列の誤差を補償するので、追尾性能を劣化させることなく処理負荷が軽減される。よって、目標数が多い場合でも対処可能となる。   As described above, according to the multi-target tracking device according to the third embodiment, prediction processing is performed for each temporary track at the entire center of gravity time of the entire observation value, and further, the error of the covariance matrix generated by time approximation is compensated. Therefore, the processing load is reduced without degrading the tracking performance. Therefore, even when the target number is large, it can be dealt with.

実施の形態4.
この発明の実施の形態4に係る多目標追尾装置について図14から図16までを参照しながら説明する。図14は、この発明の実施の形態4に係る多目標追尾装置の構成を示すブロック図である。
Embodiment 4 FIG.
A multi-target tracking device according to Embodiment 4 of the present invention will be described with reference to FIGS. FIG. 14 is a block diagram showing the configuration of the multitarget tracking apparatus according to Embodiment 4 of the present invention.

図14において、この実施の形態4に係る多目標追尾装置20は、観測値グループ分け部2と、グループ重心計算部3と、誤差共分散拡大部5と、時刻、予測諸元を含む仮航跡を複数個格納する仮航跡ファイル21と、予測処理部22と、ゲート内外判定部23と、航跡尤度評価部24と、相関決定部25と、平滑処理部26と、航跡決定部27と、目標航跡を格納する目標航跡ファイル28とが設けられている。   In FIG. 14, the multi-target tracking device 20 according to the fourth embodiment includes an observation value grouping unit 2, a group center-of-gravity calculation unit 3, an error covariance expansion unit 5, and a temporary track including time and prediction specifications. Tentative track file 21, a prediction processing unit 22, a gate inside / outside determination unit 23, a track likelihood evaluation unit 24, a correlation determination unit 25, a smoothing processing unit 26, a track determination unit 27, A target track file 28 for storing the target track is provided.

つぎに、この実施の形態4に係る多目標追尾装置の動作について図面を参照しながら説明する。図15は、この発明の実施の形態4に係る多目標追尾装置の動作を示すフローチャートである。図16は、この発明の実施の形態4に係る多目標追尾装置におけるサーチ領域を示す図である。   Next, the operation of the multi-target tracking device according to the fourth embodiment will be described with reference to the drawings. FIG. 15 is a flowchart showing the operation of the multitarget tracking apparatus according to Embodiment 4 of the present invention. FIG. 16 is a diagram showing search areas in the multi-target tracking device according to Embodiment 4 of the present invention.

「観測値入力」ステップ400において、捜索レーダ(センサ)10から観測値が入力されると、まず、「観測値グループ分け」ステップ401において、観測値グループ分け部2は、距離の近い観測値同士でグループ化する。ある観測値と別の観測値の組み合わせについての同一グループ判定については、上記実施の形態2の同名のステップと同じ方式とする。   When an observation value is input from the search radar (sensor) 10 in the “observation value input” step 400, first, in an “observation value grouping” step 401, the observation value grouping unit 2 determines that observation values having a short distance are Group by. About the same group determination about the combination of a certain observation value and another observation value, it is set as the same method as the step of the same name of the said Embodiment 2. FIG.

次に、「グループ重心計算」ステップ402において、グループ重心計算部3は、グループに所属する観測値群の重心の位置に相当する観測時刻を計算する(以下、重心時刻と呼ぶ)。重心時刻の計算方法については、上記実施の形態2の同名のステップと同じ方式とする。   Next, in a “group centroid calculation” step 402, the group centroid calculation unit 3 calculates an observation time corresponding to the position of the centroid of the observation value group belonging to the group (hereinafter referred to as the centroid time). The calculation method of the centroid time is the same as the step of the same name in the second embodiment.

次に、「予測処理」ステップ403において、予測処理部22は、仮航跡について観測値が所属するグループの重心時刻における予測位置と予測誤差共分散の計算を行う。   Next, in a “prediction process” step 403, the prediction processing unit 22 calculates the predicted position and the prediction error covariance at the centroid time of the group to which the observation value belongs for the temporary track.

次に、「誤差共分散拡大」ステップ404において、誤差共分散拡大部5は、重心時刻で予測処理したことによる誤差分を考慮して、予測誤差共分散行列を拡大する。この拡大分の計算については、上記実施の形態3の同名のステップと同じ方式とする。但し、全体重心時刻を重心時刻とする。   Next, in the “error covariance expansion” step 404, the error covariance expansion unit 5 expands the prediction error covariance matrix in consideration of the error due to the prediction processing at the centroid time. The enlargement is calculated in the same manner as the step of the same name in the third embodiment. However, the center of gravity time is defined as the center of gravity time.

次に、「ゲート内外判定」ステップ405において、ゲート内外判定部23は、仮航跡と観測値が対応付けられる可能性を、仮航跡の予測誤差共分散と観測値の観測誤差共分散を使ったカイ平方検定によって判定する。   Next, in the “gate inside / outside determination” step 405, the gate inside / outside determination unit 23 uses the prediction error covariance of the temporary track and the observation error covariance of the observation value to determine whether the temporary track and the observed value are associated with each other. Determined by chi-square test.

次に、「航跡尤度評価」ステップ406において、航跡尤度評価部24は、仮航跡と前ステップでゲート内にあると判定された観測値の組み合わせに対して、その組み合わせにより観測値を仮航跡に割り当てた場合の仮航跡の尤度を計算する。   Next, in the “wake likelihood evaluation” step 406, the wake likelihood evaluation unit 24 tentatively determines an observation value based on the combination of the temporary track and the observation value determined to be in the gate in the previous step. Calculate the likelihood of a temporary wake when assigned to a wake.

次に、「相関決定」ステップ407において、相関決定部25は、観測値と仮航跡の組み合わせのどれを最終的に選択するか、前ステップで計算された尤度を基に決定する。   Next, in the “correlation determination” step 407, the correlation determination unit 25 determines which of the combination of the observed value and the provisional track is finally selected based on the likelihood calculated in the previous step.

次に、「平滑処理」ステップ408において、平滑処理部26は、前ステップで決定された相関結果を基に、仮航跡の予測諸元と組み合わされた観測値を使って平滑諸元を計算する。すなわち、仮航跡ファイル21に格納されている仮航跡の時刻を最新時刻の値に更新する。   Next, in the “smoothing process” step 408, the smoothing processing unit 26 calculates the smooth specification using the observation value combined with the prediction specification of the temporary track based on the correlation result determined in the previous step. . That is, the temporary track time stored in the temporary track file 21 is updated to the latest time value.

そして、「航跡決定」ステップ409において、航跡決定部27は、仮航跡の尤度や観測値との相関回数を基に、仮航跡を目標航跡と見做すかどうか、決定する。決定した目標航跡を目標航跡ファイル28に格納する。   Then, in a “track determination” step 409, the track determination unit 27 determines whether or not to consider the temporary track as the target track based on the likelihood of the temporary track and the number of correlations with the observed value. The determined target track is stored in the target track file 28.

本実施の形態4に係る多目標追尾装置を図23に示す観測例に適用した場合、図16の様に観測値のグループが2つできるとすると、上記の実施の形態2の図9と同様のタイムチャートとなり、仮航跡数×グループ数(5×2)の計10回の予測計算が必要となる。これは従来方式で必要な予測計算回数25回よりも少ない。   When the multi-target tracking device according to the fourth embodiment is applied to the observation example shown in FIG. 23, if two observation value groups can be formed as shown in FIG. 16, the same as in FIG. 9 of the second embodiment described above. Thus, a total of 10 prediction calculations of the number of temporary tracks × the number of groups (5 × 2) are required. This is less than the required number of prediction calculations of 25 in the conventional method.

すなわち、この実施の形態4に係る多目標追尾装置は、センサ10で観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、1サーチで得られた観測値群をグループ分けする観測値グループ分け部2と、グループ毎にグループの重心位置に相当する観測時刻である重心時刻を計算するグループ重心計算部3と、仮航跡と観測値の組み合わせについて、仮航跡を観測値が所属するグループの前記重心時刻によって予測計算する予測処理部22と、前記重心時刻で予測処理したことによる誤差分を考慮して、予測誤差共分散行列を拡大する誤差共分散拡大部5と、仮航跡と観測値が対応付け可能かを、前記予測処理部22が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部23と、観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部24と、前記航跡尤度評価部24が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部25と、仮航跡の予測諸元と前記相関決定部25により対応付けられた観測値を使って平滑諸元を計算する平滑処理部26と、仮航跡のうち尤度が高いものを目標航跡とする航跡決定部27とを設けたものである。   That is, the multi-target tracking device according to the fourth embodiment is a multi-target tracking device that creates a target wake using observation values obtained by searching an observation region with the sensor 10, and is obtained in one search. The observed value grouping unit 2 for grouping the observed value groups, the group centroid calculating unit 3 for calculating the centroid time which is the observation time corresponding to the centroid position of the group for each group, and the combination of the temporary track and the observed value The prediction processing unit 22 that predicts and calculates the tentative track based on the centroid time of the group to which the observation value belongs, and the error co-variance matrix that expands the prediction error covariance matrix in consideration of the error caused by the prediction processing at the centroid time. In the gate, the dispersion enlarging unit 5 determines whether the temporary track and the observed value can be associated with each other by comparing the predicted value and the observed value of the temporary track calculated by the prediction processing unit 22 The determination unit 23, the track likelihood evaluation unit 24 for calculating the likelihood of the temporary track when the observation value is associated, and the temporary track and which of the temporary track, using the likelihood calculated by the track likelihood evaluation unit 24 as an index Correlation determining unit 25 for determining whether to associate observation values, smoothing processing unit 26 for calculating smooth specifications using predicted values of temporary track and observation values associated with correlation determining unit 25, Of these, a track determination unit 27 is provided which sets a target track having a high likelihood.

以上の様に、この実施の形態4に係る多目標追尾装置によれば、観測値をグループ化し、グループの重心時刻で各仮航跡に対する予測処理を行い、さらに時刻近似によって発生する共分散行列の誤差を補償するので、追尾性能を劣化させることなく処理負荷が軽減される。よって、目標数が多い場合でも対処可能となる。   As described above, according to the multi-target tracking device according to the fourth embodiment, the observation values are grouped, the prediction process is performed for each temporary track at the center of gravity time of the group, and the covariance matrix generated by time approximation is further calculated. Since the error is compensated, the processing load is reduced without degrading the tracking performance. Therefore, even when the target number is large, it can be dealt with.

実施の形態5.
この発明の実施の形態5に係る多目標追尾装置について図17から図20までを参照しながら説明する。図17は、この発明の実施の形態5に係る多目標追尾装置の構成を示すブロック図である。
Embodiment 5 FIG.
A multi-target tracking device according to Embodiment 5 of the present invention will be described with reference to FIGS. FIG. 17 is a block diagram showing a configuration of a multitarget tracking apparatus according to Embodiment 5 of the present invention.

図17において、この実施の形態5に係る多目標追尾装置20は、観測値グループ分け部2と、グループ相関計算部6と、同一グループの仮航跡を記述する仮航跡グループ情報を格納する航跡グループ情報ファイル7と、航跡グループ分け部8と、時刻、予測諸元を含む仮航跡を複数個格納する仮航跡ファイル21と、予測処理部22と、ゲート内外判定部23と、航跡尤度評価部24と、相関決定部25と、平滑処理部26と、航跡決定部27と、目標航跡を格納する目標航跡ファイル28とが設けられている。   In FIG. 17, the multi-target tracking device 20 according to the fifth embodiment includes an observation value grouping unit 2, a group correlation calculation unit 6, and a wake group that stores tentative track group information describing a tentative track of the same group. Information file 7, track grouping unit 8, temporary track file 21 for storing a plurality of temporary tracks including time and prediction specifications, prediction processing unit 22, gate inside / outside determination unit 23, track likelihood evaluation unit 24, a correlation determining unit 25, a smoothing processing unit 26, a track determining unit 27, and a target track file 28 for storing the target track.

つぎに、この実施の形態5に係る多目標追尾装置の動作について図面を参照しながら説明する。図18は、この発明の実施の形態5に係る多目標追尾装置の動作を示すフローチャートである。図19は、この発明の実施の形態5に係る多目標追尾装置におけるサーチ領域を示す図である。図20は、この発明の実施の形態5に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。   Next, the operation of the multi-target tracking device according to the fifth embodiment will be described with reference to the drawings. FIG. 18 is a flowchart showing the operation of the multitarget tracking apparatus according to Embodiment 5 of the present invention. FIG. 19 is a diagram showing search areas in the multi-target tracking device according to Embodiment 5 of the present invention. FIG. 20 is a timing chart showing the relationship between the observation time of the multi-target tracking device according to Embodiment 5 of the present invention and the predicted time of the temporary track.

「観測値入力」ステップ500において、捜索レーダ(センサ)10から観測値が入力されると、まず、「観測値グループ分け」ステップ501において、観測値グループ分け部2は、距離の近い観測値同士でグループ化する。ある観測値と別の観測値の組み合わせについての同一グループ判定については、上記実施の形態2の同名のステップと同じ方式とする。   When an observation value is input from the search radar (sensor) 10 in the “observation value input” step 500, first, in the “observation value grouping” step 501, the observation value grouping unit 2 determines that observation values having a short distance are Group by. About the same group determination about the combination of a certain observation value and another observation value, it is set as the same method as the step of the same name of the said Embodiment 2. FIG.

次に、「グループ相関処理」ステップ502において、グループ相関計算部6は、前ステップで求めた観測値のグループと、航跡グループ情報ファイル7から読み出した仮航跡のグループの相関処理を行う。これは観測値の位置と観測誤差共分散による分布を混合することによって計算された分布と、仮航跡の位置と誤差共分散行列による分布を混合することによって計算された分布を元に、「グループが同一である」という仮説を検定することによって相関の有無を決定する。   Next, in a “group correlation process” step 502, the group correlation calculation unit 6 performs a correlation process between the observed value group obtained in the previous step and the temporary track group read from the track group information file 7. This is based on the distribution calculated by mixing the distribution of the observed value and the distribution due to the observation error covariance and the distribution calculated by mixing the distribution based on the temporary track position and the error covariance matrix. The presence or absence of correlation is determined by testing the hypothesis that “is the same”.

次に、「予測処理」ステップ503において、予測処理部22は、前ステップで相関があると判定されたグループの仮航跡と観測値の組み合わせについて、観測時刻における仮航跡の予測位置と予測誤差共分散の計算を行う。   Next, in the “prediction process” step 503, the prediction processing unit 22 determines the predicted position of the temporary track and the prediction error at the observation time for the combination of the temporary track and the observed value of the group determined to be correlated in the previous step. Calculate the variance.

次に、「ゲート内外判定」ステップ504において、ゲート内外判定部23は、仮航跡と観測値が対応付けられる可能性を、仮航跡の予測誤差共分散と観測値の観測誤差共分散を使ったカイ平方検定によって判定する。   Next, in the “gate inside / outside determination” step 504, the gate inside / outside determination unit 23 uses the prediction error covariance of the temporary track and the observation error covariance of the observation value as a possibility of associating the temporary track with the observation value. Determined by chi-square test.

次に、「航跡尤度評価」ステップ505において、航跡尤度評価部24は、仮航跡と前ステップでゲート内にあると判定された観測値の組み合わせに対して、その組み合わせにより観測値を仮航跡に割り当てた場合の仮航跡の尤度を計算する。   Next, in the “wake likelihood evaluation” step 505, the wake likelihood evaluation unit 24 tentatively calculates an observation value for the combination of the temporary track and the observation value determined to be in the gate in the previous step. Calculate the likelihood of a temporary wake when assigned to a wake.

次に、「相関決定」ステップ506において、相関決定部25は、観測値と仮航跡の組み合わせのどれを最終的に選択するか、前ステップで計算された尤度を基に決定する。   Next, in the “correlation determination” step 506, the correlation determination unit 25 determines which of the combination of the observed value and the provisional track is finally selected based on the likelihood calculated in the previous step.

次に、「平滑処理」ステップ507において、平滑処理部26は、前ステップで決定された相関結果を基に、仮航跡の予測諸元と組み合わされた観測値を使って平滑諸元を計算する。すなわち、仮航跡ファイル21に格納されている仮航跡の時刻を最新時刻の値に更新する。   Next, in a “smoothing process” step 507, the smoothing processing unit 26 calculates a smooth specification using the observation value combined with the predicted specification of the temporary track based on the correlation result determined in the previous step. . That is, the temporary track time stored in the temporary track file 21 is updated to the latest time value.

次に、「航跡グループ分け」ステップ508において、航跡グループ分け部8は、平滑処理に利用した観測値が所属するグループによって、仮航跡をグループ分けする。その後、仮航跡グループ情報として、航跡グループ情報ファイル7に格納する。   Next, in a “track grouping” step 508, the track grouping unit 8 groups the temporary tracks according to the group to which the observation value used for the smoothing process belongs. Then, it is stored in the wake group information file 7 as temporary wake group information.

そして、「航跡決定」ステップ509において、仮航跡の尤度や観測値との相関回数を基に、仮航跡を目標航跡と見做すかどうか、決定する。決定した目標航跡を目標航跡ファイル28に格納する。   In a “wake determination” step 509, it is determined whether to consider the temporary track as the target track based on the likelihood of the temporary track and the number of correlations with the observed value. The determined target track is stored in the target track file 28.

本実施の形態5に係る多目標追尾装置を図23に示す観測例に適用した場合、図19の様に観測値のグループが2つ、仮航跡のグループが2つあり、観測値グループ1と航跡グループ1に相関があり、かつ観測値グループ2と航跡グループ2に相関があるとする。このとき図20のタイムチャートに示す様に、相関のあるグループ内で予測処理を行うためグループ1については仮航跡数×観測値数(3×3)の9回、グループ2については仮航跡数×観測値数(2×2)の4回の予測計算が必要となり、予測計算の回数は計13(=9+4)回となる。これは従来方式で必要な予測計算回数25回よりも少ない。   When the multi-target tracking device according to the fifth embodiment is applied to the observation example shown in FIG. 23, there are two observation value groups and two temporary track groups as shown in FIG. Assume that wake group 1 has a correlation, and observation value group 2 and wake group 2 have a correlation. At this time, as shown in the time chart of FIG. 20, in order to perform the prediction process within the correlated group, the number of temporary tracks for group 1 × the number of observed values (3 × 3) is 9 times, and the number of temporary tracks for group 2 X Prediction calculation of 4 times of the number of observed values (2 x 2) is required, and the total number of prediction calculations is 13 (= 9 + 4). This is less than the required number of prediction calculations of 25 in the conventional method.

すなわち、この実施の形態5に多目標追尾装置は、センサ10で観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、1サーチで得られた観測値群をグループ分けする観測値グループ分け部2と、観測値のグループと仮航跡のグループ間で相関処理を行うグループ相関計算部6と、仮航跡と観測値の組み合わせについて、仮航跡を観測値の観測時刻によって予測計算する予測処理部22と、仮航跡と観測値が対応付け可能かを、前記予測処理部22が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部23と、観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部24と、前記航跡尤度評価部24が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部25と、仮航跡の予測諸元と前記相関決定部25により対応付けられた観測値を使って平滑諸元を計算する平滑処理部26と、平滑処理に利用した観測値が所属するグループによって、仮航跡のグループ分けをする航跡グループ分け部8と、仮航跡のうち尤度が高いものを目標航跡とする航跡決定部27とを設けたものである。   That is, the multi-target tracking device according to the fifth embodiment is a multi-target tracking device that creates a target wake using observation values obtained by searching an observation region with the sensor 10, and is obtained in one search. Observation group grouping unit 2 for grouping observation groups, group correlation calculation unit 6 for performing correlation processing between the observation group and the temporary track group, and observing the temporary track for the combination of the temporary track and the observation value A prediction processing unit 22 that performs prediction calculation according to the observation time of the value, and a gate that determines whether the temporary track and the observation value can be associated by comparing the predicted value of the temporary track calculated by the prediction processing unit 22 with the observation value The inside / outside determination unit 23, the track likelihood evaluation unit 24 for calculating the likelihood of the temporary track when the observation values are associated, and the temporary track with which the likelihood calculated by the track likelihood evaluation unit 24 is used as an index Which A correlation determining unit 25 for determining whether to associate the measurement values, a smoothing processing unit 26 for calculating a smooth specification using the predicted values of the temporary track and the observation values associated with the correlation determining unit 25, and smoothing processing According to the group to which the observation value used belongs, there is provided a wake grouping unit 8 for grouping the tentative track, and a track determination unit 27 for setting a target track having a high likelihood among the tentative tracks.

以上の様に、この実施の形態5に係る多目標追尾装置によれば、観測値と仮航跡をグループ化し、グループ同士の相関を行ってから、相関のあるグループ内の仮航跡と観測値の組に対して予測処理を行うので、処理負荷が軽減される。よって、目標数が多い場合でも対処可能となる。   As described above, according to the multi-target tracking device according to the fifth embodiment, the observation value and the temporary track are grouped and the correlation between the groups is performed, and then the temporary track and the observation value in the correlated group are correlated. Since the prediction process is performed on the set, the processing load is reduced. Therefore, even when the target number is large, it can be dealt with.

この発明の実施の形態1に係る多目標追尾装置の構成を示すブロック図である。It is a block diagram which shows the structure of the multi-target tracking apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る多目標追尾装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the multi target tracking apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る多目標追尾装置におけるサーチ領域を示す図である。It is a figure which shows the search area | region in the multi-target tracking apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。It is a timing chart which shows the relationship between the observation time of the multitarget tracking device concerning Embodiment 1 of this invention, and the prediction time of a temporary track. この発明の実施の形態2に係る多目標追尾装置の構成を示すブロック図である。It is a block diagram which shows the structure of the multi-target tracking apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る多目標追尾装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the multi-target tracking apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る多目標追尾装置におけるサーチ領域を示す図である。It is a figure which shows the search area | region in the multi-target tracking apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る多目標追尾装置の観測値グループ分け部の同一グループ判定を説明するための図である。It is a figure for demonstrating the same group determination of the observation value grouping part of the multitarget tracking device which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。It is a timing chart which shows the relationship between the observation time of the multitarget tracking device concerning Embodiment 2 of this invention, and the prediction time of a temporary track. この発明の実施の形態3に係る多目標追尾装置の構成を示すブロック図である。It is a block diagram which shows the structure of the multi-target tracking apparatus which concerns on Embodiment 3 of this invention. この発明の実施の形態3に係る多目標追尾装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the multi-target tracking apparatus which concerns on Embodiment 3 of this invention. この発明の実施の形態3に係る多目標追尾装置におけるサーチ領域を示す図である。It is a figure which shows the search area | region in the multi-target tracking apparatus which concerns on Embodiment 3 of this invention. この発明の実施の形態3に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。It is a timing chart which shows the relationship between the observation time of the multitarget tracking device concerning Embodiment 3 of this invention, and the prediction time of a temporary track. この発明の実施の形態4に係る多目標追尾装置の構成を示すブロック図である。It is a block diagram which shows the structure of the multi-target tracking apparatus which concerns on Embodiment 4 of this invention. この発明の実施の形態4に係る多目標追尾装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the multi-target tracking apparatus which concerns on Embodiment 4 of this invention. この発明の実施の形態4に係る多目標追尾装置におけるサーチ領域を示す図である。It is a figure which shows the search area | region in the multi-target tracking apparatus which concerns on Embodiment 4 of this invention. この発明の実施の形態5に係る多目標追尾装置の構成を示すブロック図である。It is a block diagram which shows the structure of the multi-target tracking apparatus which concerns on Embodiment 5 of this invention. この発明の実施の形態5に係る多目標追尾装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the multitarget tracking device which concerns on Embodiment 5 of this invention. この発明の実施の形態5に係る多目標追尾装置におけるサーチ領域を示す図である。It is a figure which shows the search area | region in the multi-target tracking apparatus which concerns on Embodiment 5 of this invention. この発明の実施の形態5に係る多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。It is a timing chart which shows the relationship between the observation time of the multitarget tracking device concerning Embodiment 5 of this invention, and the prediction time of a temporary track. 従来の多目標追尾装置の構成を示すブロック図である。It is a block diagram which shows the structure of the conventional multi-target tracking apparatus. 従来の多目標追尾装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the conventional multi-target tracking apparatus. 従来の多目標追尾装置におけるサーチ領域を示す図である。It is a figure which shows the search area | region in the conventional multi-target tracking apparatus. 従来の多目標追尾装置の観測時刻と仮航跡の予測時刻の関係を示すタイミングチャートである。It is a timing chart which shows the relationship between the observation time of the conventional multi-target tracking apparatus and the prediction time of a temporary track.

符号の説明Explanation of symbols

1 観測値選別部、2 観測値グループ分け部、3 グループ重心計算部、4 重心計算部、5 誤差共分散拡大部、6 グループ相関計算部、7 航跡グループ情報ファイル、8 航跡グループ分け部、10 捜索レーダ(センサ)、20 多目標追尾装置、21 仮航跡ファイル、22 予測処理部、23 ゲート内外判定部、24 航跡尤度評価部、25 相関決定部、26 平滑処理部、27 航跡決定部、28 目標航跡ファイル。   1 observation value selection unit, 2 observation value grouping unit, 3 group centroid calculation unit, 4 centroid calculation unit, 5 error covariance expansion unit, 6 group correlation calculation unit, 7 track group information file, 8 track grouping unit, 10 Search radar (sensor), 20 multi-target tracking device, 21 temporary track file, 22 prediction processing unit, 23 gate inside / outside determination unit, 24 track likelihood evaluation unit, 25 correlation determination unit, 26 smoothing processing unit, 27 track determination unit, 28 Target track file.

Claims (16)

センサで観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、
予測計算を行う仮航跡と観測値の組み合わせを決定する観測値選別部と、
前記観測値選別部により決定された仮航跡と観測値の組み合わせについて、仮航跡を観測値の時刻によって予測計算する予測処理部と、
仮航跡と観測値が対応付け可能かを、前記予測処理部が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部と、
観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部と、
前記航跡尤度評価部が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部と、
仮航跡の予測諸元と前記相関決定部により対応付けられた観測値を使って平滑諸元を計算する平滑処理部と、
仮航跡のうち尤度が高いものを目標航跡とする航跡決定部と
を備えたことを特徴とする多目標追尾装置。
A multi-target tracking device that creates a target wake using observations obtained by searching an observation area with a sensor,
An observation value selector for determining a combination of the temporary track and the observation value for which the prediction calculation is performed;
For the combination of the temporary track and the observation value determined by the observation value selection unit, a prediction processing unit that predicts and calculates the temporary track according to the time of the observation value;
A gate inside / outside determination unit that determines whether the temporary track and the observed value can be associated with each other by comparing the predicted value and the observed value of the temporary track calculated by the prediction processing unit,
A wake likelihood evaluation unit that calculates the likelihood of a tentative wake when an observation value is associated,
Using the likelihood calculated by the track likelihood evaluation unit as an index, a correlation determination unit that determines which temporary track and which observation value are associated with each other,
A smoothing processing unit for calculating a smoothing specification using the predicted values of the provisional wake and the observation values associated with the correlation determination unit;
A multi-target tracking device, comprising: a track determination unit that sets a target track with a high likelihood among temporary tracks.
前記観測値選別部は、観測値の時刻と仮航跡の前サーチにおける平滑時刻の差が、事前に設定された第1の閾値を超えない場合に前記観測値を候補とする
ことを特徴とする請求項1記載の多目標追尾装置。
The observation value selection unit is characterized in that the observation value is a candidate when the difference between the time of the observation value and the smoothing time in the previous search of the temporary track does not exceed the first threshold value set in advance. The multi-target tracking device according to claim 1.
前記観測値選別部は、観測値の時刻と仮航跡の前サーチにおける平滑時刻の差が、サーチの角速度とサーチにかかる時間と目標の最大角速度から求めた許容時間差を越えない場合に前記観測値を候補とする
ことを特徴とする請求項1記載の多目標追尾装置。
The observed value selection unit determines the observed value when the difference between the time of the observed value and the smoothing time in the preliminary search before the temporary track does not exceed the allowable time difference obtained from the angular velocity of the search, the time required for the search, and the target maximum angular velocity. The multi-target tracking device according to claim 1, wherein
前記観測値選別部は、観測値の時刻と仮航跡の前サーチにおける平滑時刻の差が、サーチの角速度とサーチにかかる時間と目標の最大角速度から求めた許容時間差を越えない場合で、かつ、観測値と仮航跡の距離差が、事前に設定された第2の閾値を超えない場合に前記観測値を候補とする
ことを特徴とする請求項1記載の多目標追尾装置。
The observation value selection unit is the case where the difference between the time of the observation value and the smoothing time in the previous search of the temporary track does not exceed the allowable time difference obtained from the angular velocity of the search and the time required for the search and the target maximum angular velocity, and The multi-target tracking device according to claim 1, wherein the observation value is set as a candidate when a difference in distance between the observation value and the provisional track does not exceed a preset second threshold value.
前記目標の最大角速度は、前サーチで計算された仮航跡の角速度の中の最大値である
ことを特徴とする請求項3又は4記載の多目標追尾装置。
5. The multi-target tracking device according to claim 3, wherein the maximum angular velocity of the target is a maximum value among the angular velocities of the tentative track calculated in the previous search.
センサで観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、
1サーチで得られた観測値群をグループ分けする観測値グループ分け部と、
グループ毎にグループの重心位置に相当する観測時刻である重心時刻を計算するグループ重心計算部と、
仮航跡と観測値の組み合わせについて、仮航跡を観測値が所属するグループの前記重心時刻によって予測計算する予測処理部と、
仮航跡と観測値が対応付け可能かを、前記予測処理部が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部と、
観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部と、
前記航跡尤度評価部が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部と、
仮航跡の予測諸元と前記相関決定部により対応付けられた観測値を使って平滑諸元を計算する平滑処理部と、
仮航跡のうち尤度が高いものを目標航跡とする航跡決定部と
を備えたことを特徴とする多目標追尾装置。
A multi-target tracking device that creates a target wake using observations obtained by searching an observation area with a sensor,
An observation value grouping unit for grouping observation value groups obtained in one search;
A group centroid calculation unit for calculating a centroid time which is an observation time corresponding to the centroid position of the group for each group;
About the combination of the temporary track and the observed value, a prediction processing unit that predicts and calculates the temporary track based on the centroid time of the group to which the observed value belongs,
A gate inside / outside determination unit that determines whether the temporary track and the observed value can be associated with each other by comparing the predicted value and the observed value of the temporary track calculated by the prediction processing unit,
A wake likelihood evaluation unit that calculates the likelihood of a tentative wake when an observation value is associated,
Using the likelihood calculated by the track likelihood evaluation unit as an index, a correlation determination unit that determines which temporary track and which observation value are associated with each other,
A smoothing processing unit for calculating a smoothing specification using the predicted values of the provisional wake and the observation values associated with the correlation determination unit;
A multi-target tracking device, comprising: a track determination unit that sets a target track with a high likelihood among temporary tracks.
センサで観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、
1サーチで得られた全観測値の重心の位置に相当する観測時刻である全体重心時刻を計算する重心計算部と、
仮航跡と観測値の組み合わせについて、仮航跡を前記全体重心時刻によって予測計算する予測処理部と、
前記全体重心時刻で予測処理したことによる誤差分を考慮して、予測誤差共分散行列を拡大する誤差共分散拡大部と、
仮航跡と観測値が対応付け可能かを、前記予測処理部が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部と、
観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部と、
前記航跡尤度評価部が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部と、
仮航跡の予測諸元と前記相関決定部により対応付けられた観測値を使って平滑諸元を計算する平滑処理部と、
仮航跡のうち尤度が高いものを目標航跡とする航跡決定部と
を備えたことを特徴とする多目標追尾装置。
A multi-target tracking device that creates a target wake using observations obtained by searching an observation area with a sensor,
A center-of-gravity calculation unit that calculates an overall center-of-gravity time that is an observation time corresponding to the position of the center of gravity of all observation values obtained in one search;
For a combination of the temporary track and the observed value, a prediction processing unit that predicts and calculates the temporary track based on the total centroid time,
An error covariance expansion unit that expands the prediction error covariance matrix in consideration of the error due to the prediction processing at the overall centroid time,
A gate inside / outside determination unit that determines whether the temporary track and the observed value can be associated with each other by comparing the predicted value and the observed value of the temporary track calculated by the prediction processing unit,
A wake likelihood evaluation unit that calculates the likelihood of a tentative wake when an observation value is associated,
Using the likelihood calculated by the track likelihood evaluation unit as an index, a correlation determination unit that determines which temporary track and which observation value are associated with each other,
A smoothing processing unit for calculating a smoothing specification using the predicted values of the provisional wake and the observation values associated with the correlation determination unit;
A multi-target tracking device, comprising: a track determination unit that sets a target track with a high likelihood among temporary tracks.
センサで観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、
1サーチで得られた観測値群をグループ分けする観測値グループ分け部と、
グループ毎にグループの重心位置に相当する観測時刻である重心時刻を計算するグループ重心計算部と、
仮航跡と観測値の組み合わせについて、仮航跡を観測値が所属するグループの前記重心時刻によって予測計算する予測処理部と、
前記重心時刻で予測処理したことによる誤差分を考慮して、予測誤差共分散行列を拡大する誤差共分散拡大部と、
仮航跡と観測値が対応付け可能かを、前記予測処理部が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部と、
観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部と、
前記航跡尤度評価部が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部と、
仮航跡の予測諸元と前記相関決定部により対応付けられた観測値を使って平滑諸元を計算する平滑処理部と、
仮航跡のうち尤度が高いものを目標航跡とする航跡決定部と
を備えたことを特徴とする多目標追尾装置。
A multi-target tracking device that creates a target wake using observations obtained by searching an observation area with a sensor,
An observation value grouping unit for grouping observation value groups obtained in one search;
A group centroid calculation unit for calculating a centroid time which is an observation time corresponding to the centroid position of the group for each group;
About the combination of the temporary track and the observed value, a prediction processing unit that predicts and calculates the temporary track based on the centroid time of the group to which the observed value belongs,
In consideration of the error due to the prediction processing at the centroid time, an error covariance expansion unit that expands the prediction error covariance matrix;
A gate inside / outside determination unit that determines whether the temporary track and the observed value can be associated with each other by comparing the predicted value and the observed value of the temporary track calculated by the prediction processing unit,
A wake likelihood evaluation unit that calculates the likelihood of a tentative wake when an observation value is associated,
Using the likelihood calculated by the track likelihood evaluation unit as an index, a correlation determination unit that determines which temporary track and which observation value are associated with each other,
A smoothing processing unit for calculating a smoothing specification using the predicted values of the provisional wake and the observation values associated with the correlation determination unit;
A multi-target tracking device, comprising: a track determination unit that sets a target track with a high likelihood among temporary tracks.
センサで観測領域をサーチすることによって得られた観測値を使って目標航跡を作る多目標追尾装置であって、
1サーチで得られた観測値群をグループ分けする観測値グループ分け部と、
観測値のグループと仮航跡のグループ間で相関処理を行うグループ相関計算部と、
仮航跡と観測値の組み合わせについて、仮航跡を観測値の観測時刻によって予測計算する予測処理部と、
仮航跡と観測値が対応付け可能かを、前記予測処理部が計算した仮航跡の予測値と観測値を比較することによって判定するゲート内外判定部と、
観測値を対応付けた場合の仮航跡の尤度を計算する航跡尤度評価部と、
前記航跡尤度評価部が計算した尤度を指標として、どの仮航跡とどの観測値を対応付けるか決定する相関決定部と、
仮航跡の予測諸元と前記相関決定部により対応付けられた観測値を使って平滑諸元を計算する平滑処理部と、
平滑処理に利用した観測値が所属するグループによって、仮航跡のグループ分けをする航跡グループ分け部と、
仮航跡のうち尤度が高いものを目標航跡とする航跡決定部と
を備えたことを特徴とする多目標追尾装置。
A multi-target tracking device that creates a target wake using observations obtained by searching an observation area with a sensor,
An observation value grouping unit for grouping observation value groups obtained in one search;
A group correlation calculation unit that performs correlation processing between the observation value group and the temporary track group;
A prediction processing unit that predicts and calculates the temporary track based on the observation time of the observation value for the combination of the temporary track and the observation value;
A gate inside / outside determination unit that determines whether the temporary track and the observed value can be associated with each other by comparing the predicted value and the observed value of the temporary track calculated by the prediction processing unit,
A wake likelihood evaluation unit that calculates the likelihood of a tentative wake when an observation value is associated,
Using the likelihood calculated by the track likelihood evaluation unit as an index, a correlation determination unit that determines which temporary track and which observation value are associated with each other,
A smoothing processing unit for calculating a smoothing specification using the predicted values of the provisional wake and the observation values associated with the correlation determination unit;
A track grouping unit for grouping temporary tracks according to the group to which the observation values used for smoothing belong,
A multi-target tracking device, comprising: a track determination unit that sets a target track with a high likelihood among temporary tracks.
前記観測値グループ分け部は、2つの観測値の観測時刻の差が、事前に設定された第3の閾値を超えない場合に前記2つの観測値を同一グループとする
ことを特徴とする請求項6、8又は9記載の多目標追尾装置。
The observation value grouping unit makes the two observation values the same group when a difference between observation times of the two observation values does not exceed a preset third threshold value. The multi-target tracking device according to 6, 8, or 9.
前記観測値グループ分け部は、2つの観測値の観測時刻の差が、サーチの角速度とサーチにかかる時間と目標の最大角速度から求めた許容時間差を越えない場合に前記2つの観測値を同一グループとする
ことを特徴とする請求項6、8又は9記載の多目標追尾装置。
The observation value grouping unit assigns the two observation values to the same group when the difference between the observation times of the two observation values does not exceed the allowable time difference obtained from the angular velocity of the search, the time required for the search, and the target maximum angular velocity. The multi-target tracking device according to claim 6, 8 or 9.
前記観測値グループ分け部は、2つの観測値の観測時刻の差が、サーチの角速度とサーチにかかる時間と目標の最大角速度から求めた許容時間差を越えない場合で、かつ、前記2つの観測値の距離差が、事前に設定された第4の閾値を超えない場合に前記2つの観測値を同一グループとする
ことを特徴とする請求項6、8又は9記載の多目標追尾装置。
The observation value grouping unit is configured so that the difference between the observation times of the two observation values does not exceed the allowable time difference obtained from the angular velocity of the search, the time required for the search, and the target maximum angular velocity, and the two observation values The multi-target tracking device according to claim 6, wherein the two observation values are set to be in the same group when the difference in distance does not exceed a preset fourth threshold value.
前記誤差共分散拡大部は、事前に設定された目標の最大角速度と、前記全体重心時刻及び各観測値の時刻の差の最大値との積を予測誤差共分散行列の拡大分とする
ことを特徴とする請求項7記載の多目標追尾装置。
The error covariance enlarging unit sets the product of the target maximum angular velocity set in advance and the maximum value of the difference between the overall centroid time and each observed value as an enlargement of the prediction error covariance matrix. The multi-target tracking device according to claim 7, wherein:
前記誤差共分散拡大部は、仮航跡の速度の最大値と、前記全体重心時刻及び各観測値の時刻の差の最大値との積を予測誤差共分散行列の拡大分とする
ことを特徴とする請求項7記載の多目標追尾装置。
The error covariance enlarging unit is characterized in that a product of the maximum value of the speed of the tentative track and the maximum value of the difference between the total centroid time and the time of each observation value is an expansion of the prediction error covariance matrix. The multi-target tracking device according to claim 7.
前記誤差共分散拡大部は、事前に設定された目標の最大角速度と、前記重心時刻及び各観測値の時刻の差の最大値との積を予測誤差共分散行列の拡大分とする
ことを特徴とする請求項8記載の多目標追尾装置。
The error covariance expansion unit uses a product of a target maximum angular velocity set in advance and the maximum value of the difference between the center of gravity time and each observation value as an expansion of the prediction error covariance matrix. The multi-target tracking device according to claim 8.
前記誤差共分散拡大部は、仮航跡の速度の最大値と、前記重心時刻及び各観測値の時刻の差の最大値との積を予測誤差共分散行列の拡大分とする
ことを特徴とする請求項8記載の多目標追尾装置。
The error covariance expansion unit is characterized in that a product of the maximum value of the speed of the tentative track and the maximum value of the difference between the center of gravity time and the time of each observation value is an expansion of the prediction error covariance matrix. The multi-target tracking device according to claim 8.
JP2006119462A 2006-04-24 2006-04-24 Multi-target tracking device Active JP5047534B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2006119462A JP5047534B2 (en) 2006-04-24 2006-04-24 Multi-target tracking device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006119462A JP5047534B2 (en) 2006-04-24 2006-04-24 Multi-target tracking device

Publications (2)

Publication Number Publication Date
JP2007292553A true JP2007292553A (en) 2007-11-08
JP5047534B2 JP5047534B2 (en) 2012-10-10

Family

ID=38763305

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006119462A Active JP5047534B2 (en) 2006-04-24 2006-04-24 Multi-target tracking device

Country Status (1)

Country Link
JP (1) JP5047534B2 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010085164A (en) * 2008-09-30 2010-04-15 Mitsubishi Electric Corp Image radar system
JP2010230404A (en) * 2009-03-26 2010-10-14 Mitsubishi Electric Corp Target tracking apparatus
JP2011169751A (en) * 2010-02-18 2011-09-01 Mitsubishi Electric Corp Tracking device
JP2012247272A (en) * 2011-05-26 2012-12-13 Mitsubishi Electric Corp Tracker, computer program and tracking method
JP2015064229A (en) * 2013-09-24 2015-04-09 富士通株式会社 Tracking device, tracking method, and tracking program
JP2018054520A (en) * 2016-09-30 2018-04-05 株式会社Subaru Tracking target setting device, tracking target setting method and tracking target setting program
CN111458702A (en) * 2020-03-13 2020-07-28 四川九洲防控科技有限责任公司 Maneuvering target tracking method, system, electronic equipment and storage medium
JP2021143856A (en) * 2020-03-10 2021-09-24 三菱電機株式会社 Radar device
CN116758110A (en) * 2023-08-15 2023-09-15 中国科学技术大学 Robust multi-target tracking method under complex motion scene

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09264955A (en) * 1996-03-29 1997-10-07 Mazda Motor Corp Obstacle sensing apparatus
JP2004037262A (en) * 2002-07-03 2004-02-05 Mitsubishi Electric Corp Target tracking system
JP2004233136A (en) * 2003-01-29 2004-08-19 Mitsubishi Electric Corp Target tracking device
JP2005227184A (en) * 2004-02-13 2005-08-25 Mitsubishi Electric Corp Multi-target tracking apparatus
JP2005257386A (en) * 2004-03-10 2005-09-22 Mitsubishi Electric Corp Parallel multiple target tracking system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09264955A (en) * 1996-03-29 1997-10-07 Mazda Motor Corp Obstacle sensing apparatus
JP2004037262A (en) * 2002-07-03 2004-02-05 Mitsubishi Electric Corp Target tracking system
JP2004233136A (en) * 2003-01-29 2004-08-19 Mitsubishi Electric Corp Target tracking device
JP2005227184A (en) * 2004-02-13 2005-08-25 Mitsubishi Electric Corp Multi-target tracking apparatus
JP2005257386A (en) * 2004-03-10 2005-09-22 Mitsubishi Electric Corp Parallel multiple target tracking system

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010085164A (en) * 2008-09-30 2010-04-15 Mitsubishi Electric Corp Image radar system
JP2010230404A (en) * 2009-03-26 2010-10-14 Mitsubishi Electric Corp Target tracking apparatus
JP2011169751A (en) * 2010-02-18 2011-09-01 Mitsubishi Electric Corp Tracking device
JP2012247272A (en) * 2011-05-26 2012-12-13 Mitsubishi Electric Corp Tracker, computer program and tracking method
JP2015064229A (en) * 2013-09-24 2015-04-09 富士通株式会社 Tracking device, tracking method, and tracking program
JP2018054520A (en) * 2016-09-30 2018-04-05 株式会社Subaru Tracking target setting device, tracking target setting method and tracking target setting program
JP2021143856A (en) * 2020-03-10 2021-09-24 三菱電機株式会社 Radar device
JP7337008B2 (en) 2020-03-10 2023-09-01 三菱電機株式会社 radar equipment
CN111458702A (en) * 2020-03-13 2020-07-28 四川九洲防控科技有限责任公司 Maneuvering target tracking method, system, electronic equipment and storage medium
CN111458702B (en) * 2020-03-13 2023-06-20 四川九洲防控科技有限责任公司 Maneuvering target tracking method, maneuvering target tracking system, electronic equipment and storage medium
CN116758110A (en) * 2023-08-15 2023-09-15 中国科学技术大学 Robust multi-target tracking method under complex motion scene
CN116758110B (en) * 2023-08-15 2023-11-17 中国科学技术大学 Robust multi-target tracking method under complex motion scene

Also Published As

Publication number Publication date
JP5047534B2 (en) 2012-10-10

Similar Documents

Publication Publication Date Title
JP5047534B2 (en) Multi-target tracking device
US9453911B2 (en) Target tracking system and target tracking method
CN104424648B (en) Method for tracing object and equipment
US20190113603A1 (en) Method for predicting a motion of an object
JP2007024773A (en) Track correlation integrating device
JP4777152B2 (en) Tracking allocation device
JP5495581B2 (en) Sensor control device
CN109087331A (en) A kind of motion forecast method based on KCF algorithm
JP2010223813A (en) Radar system
JP2014044131A5 (en)
JP4801566B2 (en) Data stream monitoring device, data stream monitoring method, program thereof, and recording medium
JP4166651B2 (en) Target tracking device
JP5764867B2 (en) Target tracking device, target tracking method, and program
JP5516429B2 (en) Target tracking processing apparatus and target tracking processing method
JP5623253B2 (en) Multi-target tracking device
JP2014153162A (en) Track-to-track association device
JP2003149328A (en) Target correlation device of radar
JP3616029B2 (en) Tracking device
JP4786738B2 (en) Optimal position calculation device, optimal position calculation method, and program
JP6351497B2 (en) Target classifier
JP6303254B2 (en) Tracking device
JP4519617B2 (en) Tracking device
JP2004233136A (en) Target tracking device
JP3859908B2 (en) Target tracking device
JP2012093286A (en) Target detection device and target detection method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090224

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110719

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110720

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110914

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120313

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120405

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120717

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120718

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150727

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 5047534

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250