JP2009281865A - Multi-target tracking device - Google Patents

Multi-target tracking device Download PDF

Info

Publication number
JP2009281865A
JP2009281865A JP2008134327A JP2008134327A JP2009281865A JP 2009281865 A JP2009281865 A JP 2009281865A JP 2008134327 A JP2008134327 A JP 2008134327A JP 2008134327 A JP2008134327 A JP 2008134327A JP 2009281865 A JP2009281865 A JP 2009281865A
Authority
JP
Japan
Prior art keywords
tracking
hypothesis
mht
processing unit
correlation
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
JP2008134327A
Other languages
Japanese (ja)
Other versions
JP5197154B2 (en
Inventor
Yasushi Obata
康 小幡
Hiroshi Kameda
洋志 亀田
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 JP2008134327A priority Critical patent/JP5197154B2/en
Publication of JP2009281865A publication Critical patent/JP2009281865A/en
Application granted granted Critical
Publication of JP5197154B2 publication Critical patent/JP5197154B2/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 provide a multi-target tracking device can switch algorithm according to the degree of difficulty of tracking. <P>SOLUTION: This device includes: a MFA tracking processing part 4 generating the correlation solution with MFA processing while taking multi-dimensional correlation using the observation values over a plurality of flames from a sensor 1 when satisfying the first condition for carrying out the MHT tracking processing at the previous observation time and the second condition for determining that the tracking at the previous observation time is difficult; a MHT tracking processing part 2 generating new track while taking two-dimensional correlation between the latest observation value from the sensor 1 and the existing tracking track and generating a plurality of assumptions when either of the first and the second conditions is not satisfied; a MHT assumption monitoring part 3 determining whether the tracking is difficult or not from the situation of the generated assumption during carrying out the MHT tracking processing; and a MHT assumption setting part 5 replacing the correlation solution generated by the MFA tracking processing part 4 as an assumption of the MHT tracking processing part 2. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

この発明は、レーダ等のセンサシステムにおいて、目標位置を観測するセンサから得られる観測値を使って航跡を生成する追尾技術に関し、特に、この追尾に関して相関の問題を多次元的に解く技術において、演算時間を削減し、かつ良好の追尾性能を実現するための相関解を導出する多目標追尾装置に関するものである。   The present invention relates to a tracking technique for generating a wake using an observation value obtained from a sensor for observing a target position in a sensor system such as a radar, and more particularly, in a technique for solving a correlation problem in a multidimensional manner with respect to this tracking. The present invention relates to a multi-target tracking apparatus for deriving a correlation solution for reducing calculation time and realizing good tracking performance.

現在、センサにより目標を観測して得られた観測値を使って目標を追尾する技術については既に多くの論文、特許等の文献で取り挙げられており、その装置および方法については様々な提案がなされている。   Currently, many papers, patents, etc. have already mentioned the technology for tracking a target using observation values obtained by observing the target with a sensor, and various proposals have been made for the device and method. Has been made.

近接多目標を追尾する場合や、クラッタ等の不要信号環境下で目標を追尾する場合、相関の決定方法は追尾性能を左右する重要な部分である。相関とは、図12に示す様な、各観測時刻で複数の観測値が得られた場合、「観測時刻t1から観測時刻t3に至るどの観測値を同じ目標と見なすか」を示すものである。   When tracking a multi-proximity target or tracking a target in an unnecessary signal environment such as clutter, the correlation determination method is an important part that affects the tracking performance. The correlation indicates “what observation value from observation time t1 to observation time t3 is regarded as the same target” when a plurality of observation values are obtained at each observation time as shown in FIG. .

この相関決定の方式としては、図13に示す様に、各観測時刻毎に航跡を生成し、1観測時刻前の航跡と最新観測時刻の観測値の対応付けを行う方式が用いられてきた。この航跡と観測値の2次元の相関を決定する方法については様々なものがあるが、一つの例としては、各々の既存航跡に対して可能な相関全てを考慮して仮説を構成し、各々の相関の尤もらしさを基に仮説の信頼度を計算する多重仮説追尾法(MHT:Multiple Hypothesis Tracking)がある(例えば、非特許文献1参照)。   As a method of determining this correlation, as shown in FIG. 13, a method of generating a wake at each observation time and associating the wake one observation time before the observation value at the latest observation time has been used. There are various methods for determining the two-dimensional correlation between the track and the observed value, but one example is to construct a hypothesis considering all possible correlations for each existing track, There is a multiple hypothesis tracking (MHT) method that calculates the reliability of a hypothesis based on the likelihood of correlation (for example, see Non-Patent Document 1).

上記の技術は既存の航跡と最新の観測値の割り当てを決定する方式であるが、さらなる相関性能の向上を目指し、図14の様に、複数観測時刻に跨る観測値の相関をまとめて計算して航跡を生成するアルゴリズムとして複数フレーム割当法(MFA:Multiple Frame Assignment)がある(例えば、非特許文献2参照)。なお、ここでフレームとは、1観測時刻において得られた観測値の集合である。   The above-mentioned technique is a method for determining the allocation of the existing wake and the latest observation value. To further improve the correlation performance, as shown in Fig. 14, the correlation between observation values over multiple observation times is calculated and calculated together. There is a multiple frame assignment (MFA) method as an algorithm for generating a wake (see, for example, Non-Patent Document 2). Here, the frame is a set of observation values obtained at one observation time.

この複数フレーム割当法では、N(自然数)観測時刻に跨る観測値の相関決定をN次元の割り当て問題に変換して解く。定式化された割り当て問題を以下に示す。   In this multiple frame allocation method, the correlation determination of observation values over N (natural number) observation times is converted into an N-dimensional allocation problem and solved. The formulated assignment problem is shown below.

Figure 2009281865
Figure 2009281865

ここで、コスト行列ci1,…,iNは、観測値i,…,iから航跡を生成するのにかかるコストであり、通常、航跡尤度比の逆数を使用する。 Here, cost matrix c i1, ..., iN is observations i 1, ..., a i N a cost to generate the wake typically use the inverse of track likelihood ratio.

Figure 2009281865
Figure 2009281865

上記の割り当て問題は、コスト行列の各次元から要素を一つずつ選択して、選択した要素を最小にするという問題である。簡単な例として、N=2の場合の観測値と対応するコスト行列の例を図15に示す。このコスト行列から得られる解の例と対応する航跡の構成を図16に示す。この解によるコストは、「1」が設定された部分に対応するコスト行列の要素の和をとって、c11+c24+c32+c43となる。また、制約条件は、「1」と設定される要素の列が重複してはならないことを示す。例えば、図17の様な、相関決定は、第2観測時刻における第2要素に重複して「1」が設定されており、制約条件を満たさない。これは追尾において、1つの観測値を重複して使用していることに相当する。 The above assignment problem is a problem of selecting one element from each dimension of the cost matrix and minimizing the selected element. As a simple example, FIG. 15 shows an example of the cost matrix corresponding to the observed value when N = 2. An example of the solution obtained from this cost matrix and the configuration of the corresponding wake are shown in FIG. The cost of this solution is c 11 + c 24 + c 32 + c 43 by taking the sum of the elements of the cost matrix corresponding to the part where “1” is set. The constraint condition indicates that the column of elements set to “1” must not overlap. For example, as shown in FIG. 17, in the correlation determination, “1” is set to overlap the second element at the second observation time, and the constraint condition is not satisfied. This corresponds to the fact that one observation value is used redundantly in tracking.

以上に説明したN次元の割り当てで実現可能な(制約条件を満たす)解の数は、観測値の数に応じて指数的に増加する。そのため、コストを最小とする最適解を得るためには、処理時間が観測値の数や処理する次元数Nに応じて指数的に増加するという問題がある。これを解決するために、上記の非特許文献2で示されている方式ではラグランジェ(Lagrange)緩和法を使って、N次元の割り当て問題を2次元に緩和しながら準最適な解を得る。   The number of solutions that can be realized by the N-dimensional assignment described above (satisfying the constraint condition) increases exponentially according to the number of observation values. Therefore, in order to obtain an optimal solution that minimizes the cost, there is a problem that the processing time increases exponentially according to the number of observed values and the number N of dimensions to be processed. In order to solve this problem, the method shown in Non-Patent Document 2 described above uses a Lagrange relaxation method to obtain a sub-optimal solution while relaxing the N-dimensional assignment problem in two dimensions.

非特許文献2で示された従来の複数フレーム割当法による多目標追尾装置について図18及び図19を参照しながら説明する。図18は、複数フレーム割当法による従来の多目標追尾装置の構成を示す図である。また、図19は、複数フレーム割当法による従来の多目標追尾装置の動作を示すフローチャートである。   A conventional multi-target tracking device based on the multi-frame allocation method shown in Non-Patent Document 2 will be described with reference to FIGS. FIG. 18 is a diagram illustrating a configuration of a conventional multi-target tracking apparatus using a multiple frame allocation method. FIG. 19 is a flowchart showing the operation of the conventional multi-target tracking apparatus based on the multiple frame allocation method.

図17において、従来の複数フレーム割当法による多目標追尾装置は、目標位置を観測するセンサ1と、MFA追尾処理部4と、航跡表示部6とが設けられている。また、MFA追尾処理部4は、コスト行列計算部41と、Lagrange乗数設定部42と、Lagrange緩和解算出部43と、実現可能解算出部44と、緩和解とコスト記憶部45と、実現可能解とコスト記憶部46とが設けられている。   In FIG. 17, the conventional multi-target tracking apparatus based on the multiple frame allocation method is provided with a sensor 1 for observing a target position, an MFA tracking processing unit 4, and a track display unit 6. In addition, the MFA tracking processing unit 4 can be implemented as a cost matrix calculation unit 41, a Lagrange multiplier setting unit 42, a Lagrange relaxation solution calculation unit 43, a feasible solution calculation unit 44, a relaxation solution and cost storage unit 45. A solution and a cost storage unit 46 are provided.

まず、ステップ901において、コスト行列計算部41は、センサ1から得られた過去から最新までのN観測時刻分の観測値を入力する。   First, in step 901, the cost matrix calculation unit 41 inputs observation values for N observation times from the past to the latest obtained from the sensor 1.

次に、ステップ902において、コスト行列計算部41は、各観測時刻から一つずつ観測値を選択し、その観測値の組み合わせによって形成される航跡の尤度比を計算し、その逆数を使ってコスト行列を計算する。   Next, in step 902, the cost matrix calculation unit 41 selects observation values one by one from each observation time, calculates the likelihood ratio of the wake formed by the combination of the observation values, and uses the reciprocal thereof. Calculate the cost matrix.

次に、ステップ903において、Lagrange乗数設定部42は、制約条件を緩和するためのラグランジェ乗数の設定を行う。このラグランジェ乗数は、第3観測時刻から第N観測時刻に至る観測値全てについて設定される。   Next, in step 903, the Lagrange multiplier setting unit 42 sets a Lagrange multiplier for relaxing the constraint condition. This Lagrange multiplier is set for all observed values from the third observation time to the Nth observation time.

観測値jに関するラグランジェ乗数をujrとする。このラグランジェ乗数を設定することにより、N次元のコスト行列の第1次元〜第N次元の観測値に関する制約を2次元に緩和することができる。この緩和によりコスト関数は2次元のコスト行列に変形され、その要素は以下の式で表せる。 Let u jr be the Lagrange multiplier for the observed value j r . By setting this Lagrangian multiplier, it is possible to relax the restrictions on the first to N-th observed values of the N-dimensional cost matrix in two dimensions. By this relaxation, the cost function is transformed into a two-dimensional cost matrix, and its elements can be expressed by the following equations.

Figure 2009281865
Figure 2009281865

このステップ903における処理の一例として、初期のラグランジェ乗数を全て0に設定する。   As an example of the processing in step 903, all initial Lagrangian multipliers are set to zero.

次に、ステップ904おいて、Lagrange緩和解算出部43は、乗数によって緩和された2次元の割り当て問題の解を算出する。この2次元解の算出方法は、ハンガリー法等の高速な解導出方式があり、そのうちの何れかを使用する。導出された緩和解のコストをq(u)とし、緩和解とコスト記憶部45に記憶する。   Next, in Step 904, the Lagrange relaxation solution calculation unit 43 calculates a solution of the two-dimensional assignment problem that has been relaxed by the multiplier. This two-dimensional solution calculation method includes a high-speed solution derivation method such as the Hungarian method, and one of them is used. The derived relaxation solution cost is set to q (u) and stored in the relaxation solution and cost storage unit 45.

次に、ステップ905において、実現可能解算出部44は、上記の緩和解を修正して、全ての制約条件を満たす実現可能解を導出する。導出された実現可能解のコストをν(オーバーラインz)とし、実現可能解とコスト記憶部46に記憶する。なお、オーバーラインzは、zの上に−があることを表す。   Next, in step 905, the feasible solution calculation unit 44 corrects the above relaxation solution and derives a feasible solution that satisfies all the constraint conditions. The derived feasible solution cost is set as ν (overline z) and stored in the feasible solution and cost storage unit 46. Note that the overline z indicates that-is present on z.

次に、ステップ906において、実現可能解算出部44は、ラグランジェ緩和解のコストq(u)と実現可能解のコストν(オーバーラインz)の差をq(u)で割った値を算出し、それがあらかじめ設定した閾値パラメータ以下であったら、ステップ908へ進む。また、閾値条件を満たさない場合は、ステップ907へ進む。   Next, in Step 906, the feasible solution calculation unit 44 calculates a value obtained by dividing the difference between the cost q (u) of the Lagrangian relaxation solution and the cost ν (overline z) of the feasible solution by q (u). If it is below the preset threshold parameter, the process proceeds to step 908. If the threshold condition is not satisfied, the process proceeds to step 907.

ステップ907において、閾値条件を満たさない場合には、Lagrange乗数設定部42は、ラグランジェ乗数の再設定を行う。この再設定では、最新の緩和解における制約条件の逸脱具合に応じた乗数設定を行う。   If the threshold condition is not satisfied in step 907, the Lagrange multiplier setting unit 42 resets the Lagrange multiplier. In this resetting, a multiplier is set according to the deviation of the constraint condition in the latest relaxation solution.

このラグランジェ乗数再設定ステップ907を実行した場合は、以降、ラグランジェ緩和解算出ステップ904と、実現可能解算出ステップ905が、閾値条件を満たすまで繰り返される。   When this Lagrange multiplier resetting step 907 is executed, the Lagrangee relaxation solution calculation step 904 and the feasible solution calculation step 905 are repeated until the threshold condition is satisfied.

ステップ908において、航跡表示部6は、実現可能解をこのN次元割り当て問題の準最適解として出力する。   In step 908, the wake display unit 6 outputs the feasible solution as a suboptimal solution of this N-dimensional assignment problem.

上記の複数フレーム割当法MFAは、相関性能が高いものの、演算負荷が重いという問題があった。この問題を解決するために、複数フレーム割当法による追尾方式と2次元相関による追尾方式を併用する技術が提案されている(例えば、特許文献1参照)。   Although the above-described multiple frame allocation method MFA has high correlation performance, it has a problem of heavy computation load. In order to solve this problem, a technique has been proposed in which a tracking method based on the multiple frame allocation method and a tracking method based on two-dimensional correlation are used in combination (for example, see Patent Document 1).

図20は、この併用技術の構成を示す図である。この併用技術では、観測時刻を観測間隔で割った値(処理回数)があらかじめパラメータとして設定されていた値nの倍数である場合には複数フレーム割り当て処理を行い、そうでない場合には2次元相関算出処理を行う。   FIG. 20 is a diagram showing the configuration of this combined technique. In this combined technique, when the value obtained by dividing the observation time by the observation interval (number of processing times) is a multiple of the value n set in advance as a parameter, multiple frame allocation processing is performed, and otherwise, two-dimensional correlation is performed. Perform the calculation process.

特開2007−212244号公報JP 2007-212244 A Donald B. Reid “An Algorithm for Tracking Multiple Targets”(IEEE Transactions on Automatic Control, Vol.AC-24, No6, December,1979)Donald B. Reid “An Algorithm for Tracking Multiple Targets” (IEEE Transactions on Automatic Control, Vol.AC-24, No6, December, 1979) A. B. Poore and A. J. Robertson III “A New Multi-dimensional Data Association Algorithm for Multisensor-Multitarget Tracking”(Proc. Of SPIE vol.2561,1995)A. B. Poore and A. J. Robertson III “A New Multi-dimensional Data Association Algorithm for Multisensor-Multitarget Tracking” (Proc. Of SPIE vol.2561,1995)

追尾における相関処理は、センサの探知状況が良好で周囲に他の目標や不要信号が存在しない場合は正解を導くのが容易であるが、目標の探知抜けが起きたり、周囲に不要信号が発生したりする場合には困難となる。容易な場合には従来の2次元相関による追尾方式で十分であり、困難な場合にはMFA追尾方式を適用した方が、一般に正解が得られる確率が高い。よって、追尾が容易な場合に2次元相関を適用し、困難な場合にMFA追尾処理を適用する切り替え方式が理想的である。しかし、上記の特許文献1の従来技術は、二つの追尾アルゴリズムの切り替え時期が固定的であり、追尾の難易度によってアルゴリズムの切り替えを行うことができないという問題点があった。   Correlation processing in tracking is easy to find the correct answer when the detection status of the sensor is good and there are no other targets or unnecessary signals in the surroundings, but target detection is lost or unnecessary signals are generated in the surroundings It becomes difficult to do. If easy, the conventional tracking method using two-dimensional correlation is sufficient, and if difficult, the probability of obtaining a correct answer is generally higher when the MFA tracking method is applied. Therefore, a switching method in which two-dimensional correlation is applied when tracking is easy and MFA tracking processing is applied when difficult is ideal. However, the prior art disclosed in Patent Document 1 has a problem that the switching time of the two tracking algorithms is fixed, and the algorithm cannot be switched depending on the difficulty of tracking.

この発明は、上述のような課題を解決するためになされたもので、その目的は、追尾の難易度に応じてアルゴリズムの切り換えを行うことができ、演算負荷の削減と正確な相関解の導出を両立させることができる多目標追尾装置を得るものである。   The present invention has been made to solve the above-described problems, and its purpose is to switch algorithms according to the difficulty of tracking, to reduce the calculation load and to derive an accurate correlation solution. A multi-target tracking device capable of satisfying both of the above is obtained.

この発明に係る多目標追尾装置は、前観測時刻で多重仮説追尾法による追尾処理を実行する第1の条件、及び前観測時刻で追尾が困難であると判定する第2の条件を満たす場合には、センサからの複数のフレームに跨る観測値を用いて多次元の相関を取りながら、複数フレーム割当法による処理で相関解を生成するMFA追尾処理部と、前記第1及び第2の条件の何れかの条件を満たさない場合には、センサからの最新観測値と既存追尾航跡の2次元の相関を取り、複数の仮説を生成しながら新たな航跡を生成するMHT追尾処理部と、前記MHT追尾処理部による追尾処理実行中に、生成された仮説の状況から追尾が困難であるかどうかを判定するMHT仮説監視部と、前記MFA追尾処理部で生成された相関解を前記MHT追尾処理部の仮説として置き換えるMHT仮説設定部とを設けたものである。   The multi-target tracking device according to the present invention satisfies the first condition for executing the tracking process by the multiple hypothesis tracking method at the previous observation time and the second condition for determining that tracking is difficult at the previous observation time. Includes an MFA tracking processing unit that generates a correlation solution by a process based on a multi-frame allocation method while obtaining a multidimensional correlation using observation values across a plurality of frames from the sensor, and the first and second conditions. If any of the conditions is not satisfied, the MHT tracking processing unit that generates a new track while generating a plurality of hypotheses by taking a two-dimensional correlation between the latest observed value from the sensor and the existing tracking track, and the MHT An MHT hypothesis monitoring unit that determines whether tracking is difficult from the status of the generated hypothesis during execution of the tracking processing by the tracking processing unit, and the correlation solution generated by the MFA tracking processing unit is the MHT tracking processing unit. of It is provided with a and MHT hypothesis setting unit to replace a theory.

この発明に係る多目標追尾装置は、追尾の難易度に応じてアルゴリズムの切り換えを行うことができ、演算負荷の削減と正確な相関解の導出を両立させることができるという効果を奏する。   The multi-target tracking device according to the present invention can perform algorithm switching according to the tracking difficulty level, and has the effect of being able to achieve both reduction in calculation load and accurate derivation of a correlation solution.

実施の形態1.
この発明の実施の形態1に係る多目標追尾装置について図1から図7までを参照しながら説明する。図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 diagram showing the configuration of the multi-target tracking device according to Embodiment 1 of the present invention. In addition, in each figure, the same code | symbol shows the same or equivalent part.

図1において、この発明の実施の形態1に係る多目標追尾装置は、目標位置を観測するセンサ1と、センサ1からの最新観測値と既存追尾航跡の2次元の相関を取り、複数の仮説を生成しながら新たな航跡を生成するMHT追尾処理部2と、MHT追尾処理部2による追尾処理実行中に、生成された仮説の状況から追尾が困難であるかどうかを判定するMHT仮説監視部3と、センサ1からの複数のフレームに跨る観測値を用いて多次元の相関を取りながら、複数フレーム割当法による処理で相関解を生成するMFA追尾処理部4と、MFA追尾処理部4で生成された相関解をMHT追尾処理部2の仮説として置き換えるMHT仮説設定部5と、航跡をユーザに表示する航跡表示部6と、航跡を記憶する航跡記憶部21と、仮説を記憶する仮説記憶部22と、Nフレーム分の観測値、航跡を記憶するNフレームデータ記憶部23とが設けられている。   In FIG. 1, a multi-target tracking device according to Embodiment 1 of the present invention takes a two-dimensional correlation between a sensor 1 for observing a target position, the latest observed value from the sensor 1 and an existing tracking track, and a plurality of hypotheses. MHT tracking processing unit 2 that generates a new track while generating MHT, and an MHT hypothesis monitoring unit that determines whether tracking is difficult from the status of the generated hypothesis during execution of tracking processing by the MHT tracking processing unit 2 3 and an MFA tracking processing unit 4 that generates a correlation solution by a process using a multiple frame allocation method while obtaining a multidimensional correlation using observation values over a plurality of frames from the sensor 1. An MHT hypothesis setting unit 5 that replaces the generated correlation solution as a hypothesis of the MHT tracking processing unit 2, a wake display unit 6 that displays the wake to the user, a wake storage unit 21 that stores the wake, and a temporary storage that stores the hypothesis. A storage unit 22, the observed value of the N frames, and the N-frame data storing unit 23 for storing the track is provided.

つぎに、この実施の形態1に係る多目標追尾装置の動作について図面を参照しながら説明する。図2は、この発明の実施の形態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.

最初に、ステップ101において、センサ1から観測された観測値が装置に入力され、それを観測値記憶部(図示せず)に蓄積しておく。   First, in step 101, the observation value observed from the sensor 1 is input to the apparatus and stored in an observation value storage unit (not shown).

次に、ステップ102において、以下の条件(a)及び(b)を共に満たす場合には、複数フレーム割当法(MFA)による追尾処理を選択し、ステップ106に移行する。
(a)前観測時刻で多重仮説追尾法(MHT)による追尾処理が実行された場合。
(b)前観測時刻で追尾が困難であると判定された場合。
Next, in step 102, when both of the following conditions (a) and (b) are satisfied, the tracking process by the multiple frame allocation method (MFA) is selected, and the process proceeds to step 106.
(A) When tracking processing by multiple hypothesis tracking method (MHT) is executed at the previous observation time.
(B) When it is determined that tracking is difficult at the previous observation time.

上記の何れかの条件を満たさない場合には、多重仮説追尾法による追尾処理を選択し、ステップ103の処理を行う。なお、最初の観測時刻では無条件に多重仮説追尾法による追尾処理を選択する。   If any of the above conditions is not satisfied, the tracking process by the multiple hypothesis tracking method is selected, and the process of step 103 is performed. At the first observation time, tracking processing by the multiple hypothesis tracking method is selected unconditionally.

多重仮説追尾法による処理を選択した場合、まず、ステップ103において、MHT追尾処理部2は、従来の多重仮説追尾法による処理と同じ処理を実行する。その処理内容の詳細は、上述したように、非特許文献1で示されている。   When processing by the multiple hypothesis tracking method is selected, first, in step 103, the MHT tracking processing unit 2 executes the same processing as the processing by the conventional multiple hypothesis tracking method. The details of the processing contents are shown in Non-Patent Document 1, as described above.

その後、ステップ104において、MHT追尾処理部2は、ある観測時刻における最良仮説に含まれる航跡について、その子孫をその観測時刻からN回の観測時刻が経過する間、全て生成し、航跡記憶部21に記憶する。また、仮説についても、仮説記憶部22に記憶する。   Thereafter, in step 104, the MHT tracking processing unit 2 generates all descendants of the wake included in the best hypothesis at a certain observation time during the N observation times from the observation time, and the wake storage unit 21. To remember. The hypothesis is also stored in the hypothesis storage unit 22.

ここで、航跡、仮説の親子関係について説明する。ある観測時刻において生成された航跡Tを基にして、その後得られた何れかの観測値を用いて追尾フィルタ計算によって運動諸元推定値を更新することにより生成された航跡は全てTの子孫である。ある観測時刻の航跡Tの例を図3(a)に示す。その後、2観測時刻経過した後で得られる航跡例を図3(b)に示すが、4つの航跡のうち、T11とT12はTの子孫であり、T21とT22はTの子孫ではない。すなわち、航跡を組合せて生成する仮説については信頼度に応じて取捨選択するが、航跡については全て維持する点で、従来の多重仮説追尾法による追尾処理と異なる。   Here, the wake and hypothesis parent-child relationship will be described. Based on the wake T generated at a certain observation time, all the wakes generated by updating the motion specification estimates by tracking filter calculation using any of the observation values obtained thereafter are descendants of T. is there. An example of the wake T at a certain observation time is shown in FIG. Thereafter, an example of a wake obtained after two observation times have elapsed is shown in FIG. 3B. Of the four wakes, T11 and T12 are T descendants, and T21 and T22 are not T descendants. That is, the hypothesis generated by combining the wakes is selected according to the reliability, but all the wakes are maintained, which is different from the tracking process by the conventional multiple hypothesis tracking method.

多重仮説追尾法による追尾処理を行った場合、ステップ105において、MHT仮説監視部3は、次の観測時刻で多重仮説追尾法の追尾が困難であるか否かの判定を行う。この判定は以下の(1)〜(3)の何れかによって行う。   When tracking processing by the multiple hypothesis tracking method is performed, in step 105, the MHT hypothesis monitoring unit 3 determines whether or not tracking of the multiple hypothesis tracking method is difficult at the next observation time. This determination is performed by any one of the following (1) to (3).

(1)MHT追尾処理部2が生成した仮説のうち、第1位の仮説H1と第2位の仮説H2の信頼度の差が閾値th以上の場合、すなわち、   (1) Among the hypotheses generated by the MHT tracking processing unit 2, when the difference in reliability between the first hypothesis H1 and the second hypothesis H2 is greater than or equal to the threshold th,

Figure 2009281865
Figure 2009281865

である場合は、多重仮説追尾法による追尾が困難でないとする。反対に、   In this case, it is assumed that tracking by the multiple hypothesis tracking method is not difficult. Conversely,

Figure 2009281865
Figure 2009281865

である場合には、多重仮説追尾法による追尾が困難であるとする。図4に、第1位の仮説H1と第2位の仮説H2の信頼度の推移の例を示す。観測時刻(a)、(b)では、二つの仮説の信頼度の差は閾値thを上回るため、多重仮説追尾法の追尾が困難ではないと判定される。観測時刻(c)では、二つの仮説の信頼度の差が閾値thを下回るため、多重仮説追尾法の追尾が困難であると判定される。   In this case, it is assumed that tracking by the multiple hypothesis tracking method is difficult. FIG. 4 shows an example of the transition of the reliability of the first hypothesis H1 and the second hypothesis H2. At the observation times (a) and (b), since the difference in reliability between the two hypotheses exceeds the threshold th, it is determined that tracking of the multiple hypothesis tracking method is not difficult. At the observation time (c), since the difference in reliability between the two hypotheses is below the threshold th, it is determined that tracking by the multiple hypothesis tracking method is difficult.

(2)MHT追尾処理部2では、最終的に維持する仮説数m(自然数)を事前設定のパラメータとする場合がある。このパラメータが指定されている場合に、従来の多重仮説追尾法では、第m位の仮説Hmと第m+1位の仮説Hm+1の信頼度の差が閾値th以上の場合、すなわち、   (2) The MHT tracking processing unit 2 may use the hypothesis number m (natural number) to be finally maintained as a preset parameter. When this parameter is specified, in the conventional multiple hypothesis tracking method, when the reliability difference between the m-th hypothesis Hm and the m + 1-th hypothesis Hm + 1 is greater than or equal to the threshold th, that is,

Figure 2009281865
Figure 2009281865

である場合には、多重仮説追尾法による追尾が困難でないとする。反対に、   In this case, it is assumed that tracking by the multiple hypothesis tracking method is not difficult. Conversely,

Figure 2009281865
Figure 2009281865

である場合には、多重仮説追尾法による追尾が困難であるとする。図5に、第m位の仮説Hmと第m+1位の仮説Hm+1の信頼度の推移の例を示す。観測時刻(a)、(b)では、二つの仮説の信頼度の差は閾値thを上回るため、多重仮説追尾法の追尾が困難ではないと判定される。観測時刻(c)では、二つの仮説の信頼度の差が閾値thを下回るため、多重仮説追尾法の追尾が困難であると判定される。   In this case, it is assumed that tracking by the multiple hypothesis tracking method is difficult. FIG. 5 shows an example of transition of reliability of the m-th hypothesis Hm and the m + 1-th hypothesis Hm + 1. At the observation times (a) and (b), since the difference in reliability between the two hypotheses exceeds the threshold th, it is determined that tracking of the multiple hypothesis tracking method is not difficult. At the observation time (c), since the difference in reliability between the two hypotheses is below the threshold th, it is determined that tracking by the multiple hypothesis tracking method is difficult.

(3)MHT追尾処理部2が生成した仮説のうち、第1位と第2位の仮説の子孫の順位が前観測時刻から入れ替わったら多重仮説追尾法による追尾が困難であるとする。なお、ステップ104で航跡の親子関係について説明したが、仮説についてもその参照する航跡の親子関係に従って、親子関係を定義することができる。例えば、ある観測時刻t1における航跡T1を採択する仮説H1があり、その後の観測時刻t2において航跡T1の子孫であるT11を採択する仮説H11がある場合は、仮説H11は仮説H1の子孫である。ある観測時刻において第1位の仮説H1と第2位の仮説H2があったとし、各々の子孫の信頼度の推移の例を図6に示す。観測時刻(a)、(b)では、二つの子孫の間で信頼度の順位は入れ替わっていないため、多重仮説追尾法の追尾が困難ではないと判定される。一方、観測時刻(c)では、仮説H2の子孫の信頼度が仮説H1の子孫の信頼度を上回るため、多重仮説追尾法の追尾が困難であると判定される。   (3) Of the hypotheses generated by the MHT tracking processing unit 2, it is assumed that tracking by the multiple hypothesis tracking method is difficult if the ranks of descendants of the first and second hypotheses are switched from the previous observation time. In addition, although the parent-child relationship of a wake was demonstrated in step 104, a parent-child relationship can be defined also about a hypothesis according to the parent-child relationship of the wake referred to. For example, if there is a hypothesis H1 that adopts a wake T1 at a certain observation time t1, and there is a hypothesis H11 that adopts T11 that is a descendant of the wake T1 at a subsequent observation time t2, the hypothesis H11 is a descendant of the hypothesis H1. FIG. 6 shows an example of transition of reliability of each descendant assuming that there is a first hypothesis H1 and a second hypothesis H2 at a certain observation time. At the observation times (a) and (b), the reliability rankings are not interchanged between the two descendants, so that it is determined that tracking of the multiple hypothesis tracking method is not difficult. On the other hand, at the observation time (c), since the reliability of the descendants of the hypothesis H2 exceeds the reliability of the descendants of the hypothesis H1, it is determined that tracking of the multiple hypothesis tracking method is difficult.

続いて、上記のステップ102において説明した2条件(a)及び(b)を共に満たし、複数フレーム割当法による追尾処理を選択した場合について説明する。   Next, a case where both of the two conditions (a) and (b) described in step 102 above are satisfied and the tracking process by the multiple frame allocation method is selected will be described.

まず、ステップ106において、MFA追尾処理部4は、従来の複数フレーム割当法による追尾処理を実行する。この処理内容については様々な技術が適用可能であるが、上述した非特許文献2でその一例が説明されている。なお、複数フレーム割当法による追尾処理を行う場合、複数フレーム割当法で相関処理を行う最新Nフレーム分の観測値であるスライディングウィンドウの第1フレームは、観測値でなく航跡の平滑値とする必要がある。この平滑値として、前観測時刻において多重仮説追尾法による処理によって生成された最も信頼度が高い仮説に含まれる航跡のスライディングウィンドウの最初の時刻の平滑値を適用する。この第1フレームの平滑値を実現する航跡の子孫は、MHT追尾処理部2側によって全て生成されているため、MFA追尾処理部4側ではそれらの航跡を用いて、Nフレーム分の相関のコスト行列を生成する。   First, in step 106, the MFA tracking processing unit 4 executes tracking processing by a conventional multiple frame allocation method. Various techniques can be applied to the processing content, but an example is described in Non-Patent Document 2 described above. In addition, when tracking processing by the multiple frame allocation method is performed, the first frame of the sliding window, which is an observation value for the latest N frames subjected to correlation processing by the multiple frame allocation method, needs to be a smooth value of a wake instead of an observation value. There is. As the smooth value, the smooth value of the first time of the sliding window of the wake included in the hypothesis having the highest reliability generated by the processing by the multiple hypothesis tracking method at the previous observation time is applied. Since all descendants of the wake that realize the smooth value of the first frame are generated by the MHT tracking processing unit 2 side, the MFA tracking processing unit 4 side uses these wakes to calculate the correlation cost for N frames. Generate a matrix.

次に、ステップ107において、MHT仮説設定部5は、MFA追尾処理部4が生成した相関解をMHT追尾処理部2が管理する第1位の仮説とする。また、MHT追尾処理部2が前観測時刻に生成した第1位〜第m−1位の仮説を最も単純な相関処理で延長する。ただし、MFA追尾処理部4が生成した相関解が第1位となるため、延長仮説は前観測時刻から順位を一つずつ下げ、第2位〜第m位の仮説とする。この例を図7に示す。   Next, in Step 107, the MHT hypothesis setting unit 5 sets the correlation solution generated by the MFA tracking processing unit 4 as the first hypothesis managed by the MHT tracking processing unit 2. In addition, the hypothesis of the 1st place to the (m-1) th place generated by the MHT tracking processing unit 2 at the previous observation time is extended by the simplest correlation process. However, since the correlation solution generated by the MFA tracking processing unit 4 is first, the extended hypothesis is lowered from the previous observation time one by one to be the second to mth hypotheses. An example of this is shown in FIG.

この単純な相関処理による仮説の延長方法を説明する。各仮説に含まれる航跡の子孫で、かつMFA追尾処理部4が生成した航跡のうち最も尤度が高い航跡を選び、延長仮説に含める。この結果、一つの延長仮説に含まれる複数の航跡間で、使用する観測値に重複が生じた場合、尤度の高い航跡を残す。尤度の低い航跡については別の子孫の航跡を選択する。また、MFA追尾処理部4が生成した全航跡の中に、MHT追尾処理部2の仮説に含まれる航跡の子航跡が存在しない場合は、延長仮説中でメモリトラック航跡(現観測時刻において探知抜けを起こした航跡)とする。   A hypothesis extension method based on this simple correlation process will be described. A track that is descendant of a track included in each hypothesis and that has the highest likelihood is selected from the tracks generated by the MFA tracking processing unit 4 and included in the extended hypothesis. As a result, when there is an overlap in observation values to be used among a plurality of wakes included in one extension hypothesis, a wake having a high likelihood is left. For a track with a low likelihood, a track of another descendant is selected. Further, if there is no child track of the track included in the hypothesis of the MHT tracking processing unit 2 among all the tracks generated by the MFA tracking processing unit 4, the memory track track (not detected at the current observation time) is included in the extended hypothesis. Wake).

なお、第2位以下の仮説は全てMFA追尾処理部4で計算済みの航跡を使って構成されているため、MFA追尾処理部4が最も尤度の高い相関解を導出する限り信頼度の順位が入れ替わることはないが、第2位の仮説と同一内容となる可能性がある。また、演算時間が十分でない場合、さらに下位の仮説と同一内容となる可能性がある。そのため、後処理として、まずMFA追尾処理部4の下位とMHT追尾処理部2の仮説を比較し、一致した場合はMFA追尾処理部4の解の方を削除する。また、信頼度の上下関係が仮説の順位と一致しない場合には、仮説の順位を入れ換える。なお、複数フレーム割当法による追尾処理を選択した観測時刻の次の観測時刻では必ず多重仮説追尾法による処理を選択する。   Since all hypotheses of the second and lower ranks are configured using wakes already calculated by the MFA tracking processing unit 4, as long as the MFA tracking processing unit 4 derives the most likely correlation solution, the reliability ranking May not be replaced, but may have the same content as the second hypothesis. In addition, if the calculation time is not sufficient, the content may be the same as that of a lower hypothesis. Therefore, as post-processing, first, hypotheses of the MFA tracking processing unit 4 and the hypothesis of the MHT tracking processing unit 2 are compared. If they match, the solution of the MFA tracking processing unit 4 is deleted. In addition, when the hierarchical relationship of reliability does not match the hypothesis rank, the hypothesis rank is switched. Note that the processing based on the multiple hypothesis tracking method is always selected at the observation time next to the observation time at which the tracking processing by the multiple frame allocation method is selected.

最終のステップ108では、航跡表示部6は、MHT追尾処理部2が管理する仮説群の中の第1位の仮説に含まれる航跡をユーザに表示する。   In the final step 108, the wake display unit 6 displays the wake included in the first hypothesis in the hypothesis group managed by the MHT tracking processing unit 2 to the user.

以上のように、この実施の形態1によれば、多重仮説追尾法(MHT)による処理と複数フレーム割当法(MFA)による処理を使い分け、さらに多重仮説追尾法による処理における仮説の状況によって追尾の難易を判定し、追尾が困難である場合には複数フレーム割当法による処理を行うため、処理負荷の軽減と正確な相関解の導出が可能となる。   As described above, according to the first embodiment, the processing based on the multiple hypothesis tracking method (MHT) and the processing based on the multiple frame allocation method (MFA) are selectively used, and the tracking is performed according to the hypothesis status in the processing based on the multiple hypothesis tracking method. When difficulty is determined and tracking is difficult, processing by the multiple frame allocation method is performed, so that the processing load can be reduced and an accurate correlation solution can be derived.

実施の形態2.
この発明の実施の形態2に係る多目標追尾装置について図8及び図9を参照しながら説明する。図8は、この発明の実施の形態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. 8 is a diagram showing the configuration of the multi-target tracking device according to Embodiment 2 of the present invention.

図8において、この発明の実施の形態2に係る多目標追尾装置は、センサ1と、MHT追尾処理部2と、MFA追尾処理部4と、航跡表示部6と、MHT追尾処理部2によって生成された第1位の仮説とMFA追尾処理部4によって生成された相関解の尤度を比較する解比較部7と、解比較部7による比較結果として、MFA追尾処理部4の相関解の尤度が大きい場合には、MHT追尾処理部2の第1位の仮説をMFA追尾処理部4の相関解で置き換えるMHT仮説修正部8と、解比較部7による比較結果として、MHT追尾処理部2の第1位の仮説の尤度が大きい場合には、MFA追尾処理部4の相関解をMHT追尾処理部2の第1位の仮説で置き換えるMFA相関解修正部9と、航跡記憶部21と、仮説記憶部22と、Nフレームデータ記憶部23とが設けられている。   In FIG. 8, the multi-target tracking device according to Embodiment 2 of the present invention is generated by the sensor 1, the MHT tracking processing unit 2, the MFA tracking processing unit 4, the wake display unit 6, and the MHT tracking processing unit 2. As a comparison result by the solution comparison unit 7 that compares the first hypothesis and the likelihood of the correlation solution generated by the MFA tracking processing unit 4 and the comparison result by the solution comparison unit 7, the likelihood of the correlation solution of the MFA tracking processing unit 4 When the degree is high, the MHT hypothesis correcting unit 8 that replaces the first hypothesis of the MHT tracking processing unit 2 with the correlation solution of the MFA tracking processing unit 4 and the comparison result by the solution comparing unit 7 are used as the comparison result of the MHT tracking processing unit 2. When the likelihood of the first hypothesis is large, the MFA correlation solution correcting unit 9 that replaces the correlation solution of the MFA tracking processing unit 4 with the first hypothesis of the MHT tracking processing unit 2, the track storage unit 21, , Hypothesis storage unit 22 and N frame data And 憶部 23 is provided.

つぎに、この実施の形態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. 9 is a flowchart showing the operation of the multitarget tracking apparatus according to Embodiment 2 of the present invention.

最初に、ステップ201において、センサ1から観測された観測値が装置に入力され、それを観測値記憶部(図示せず)に蓄積しておく。   First, in step 201, an observation value observed from the sensor 1 is input to the apparatus and is stored in an observation value storage unit (not shown).

次に、ステップ202及び203において、多重仮説追尾法による追尾処理と、複数フレーム割当法による追尾処理を並行して実行する。   Next, in steps 202 and 203, the tracking process by the multiple hypothesis tracking method and the tracking process by the multiple frame allocation method are executed in parallel.

ステップ202では、MHT追尾処理部2は、従来の多重仮説追尾法による処理と同じ処理を実行する。その処理内容の詳細は、上述したように、非特許文献1で示されている。   In step 202, the MHT tracking processing unit 2 executes the same processing as the processing by the conventional multiple hypothesis tracking method. The details of the processing contents are shown in Non-Patent Document 1, as described above.

また、ステップ203では、MFA追尾処理部4は、Nフレーム分の多次元相関による追尾処理を行う。この処理内容については様々な技術が適用可能であるが、上述したように、非特許文献2でその一例が示されている。   In step 203, the MFA tracking processing unit 4 performs tracking processing based on multi-dimensional correlation for N frames. Various techniques can be applied to the processing content, but as described above, Non-Patent Document 2 shows an example.

次に、ステップ204において、解比較部7は、多重仮説追尾法による処理で生成された第1位の仮説の尤度と、複数フレーム割当法による処理で生成された相関解の尤度を比較する。   Next, in step 204, the solution comparison unit 7 compares the likelihood of the first hypothesis generated by the process based on the multiple hypothesis tracking method and the likelihood of the correlation solution generated by the process based on the multiple frame allocation method. To do.

比較した結果として、複数フレーム割当法の相関解の尤度が多重仮説追尾法の第1位の仮説の尤度を上回り、かつその差が閾値thを上回る場合には、ステップ205へ移行する。 As a result of the comparison, when the likelihood of the correlation solution of the multiple frame allocation method exceeds the likelihood of the first hypothesis of the multiple hypothesis tracking method and the difference exceeds the threshold th 1 , the process proceeds to step 205. .

ステップ205において、MHT仮説修正部8は、多重仮説追尾法の第1位の仮説を複数フレーム割当法の相関解と同一となる様に置き換えを行う。まず、MFA追尾処理部4が生成した相関解をMHT追尾処理部2が管理する第1位の仮説とする。また、MHT追尾処理部2が前観測時刻に生成した第1位〜第m−1位の仮説については順位を一つずつ下げて第2位〜第m位の仮説とするか、全てを破棄する。   In step 205, the MHT hypothesis correcting unit 8 replaces the first hypothesis of the multiple hypothesis tracking method so as to be the same as the correlation solution of the multiple frame allocation method. First, the correlation solution generated by the MFA tracking processing unit 4 is set as the first hypothesis managed by the MHT tracking processing unit 2. In addition, the first to m-1th hypotheses generated at the previous observation time by the MHT tracking processing unit 2 are lowered one by one to be the second to mth hypotheses or all are discarded. To do.

解比較部7の比較結果として、多重仮説追尾法の第1位の仮説の尤度が複数フレーム割当法の相関解の尤度を上回り、かつその差が閾値−thを上回る場合には、ステップ206へ移行する。 As a comparison result of the solution comparison unit 7, when the likelihood of the first hypothesis of the multiple hypothesis tracking method exceeds the likelihood of the correlation solution of the multiple frame allocation method and the difference exceeds the threshold −th 2 , The process proceeds to step 206.

ステップ206において、MFA相関解修正部9は、MFA追尾処理部4の相関解を多重仮説追尾法の第1位の仮説と同一となる様に置き換えを行う。   In step 206, the MFA correlation solution correction unit 9 replaces the correlation solution of the MFA tracking processing unit 4 so as to be the same as the first hypothesis of the multiple hypothesis tracking method.

最終のステップ207において、航跡表示部6は、MHT追尾処理部2の相関解に含まれる航跡をユーザに表示する。   In the final step 207, the wake display unit 6 displays the wake included in the correlation solution of the MHT tracking processing unit 2 to the user.

以上のように、この実施の形態2によれば、多重仮説追尾法(MHT)による処理と複数フレーム割当法(MFA)による処理を並行して実行し、さらに両者が生成した仮説および相関解をその尤度の大小に応じて互いに置き換えるため、何れかの処理で追尾性能が劣化する場合においても常に正確な相関解の導出が可能となる。   As described above, according to the second embodiment, the processing based on the multiple hypothesis tracking method (MHT) and the processing based on the multiple frame allocation method (MFA) are executed in parallel, and the hypothesis and the correlation solution generated by both are executed. Since they are replaced with each other in accordance with the magnitude of the likelihood, an accurate correlation solution can always be derived even when the tracking performance deteriorates in any process.

実施の形態3.
この発明の実施の形態3に係る多目標追尾装置について図10及び図11を参照しながら説明する。図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 diagram showing a configuration of a multi-target tracking device according to Embodiment 3 of the present invention.

図10において、この発明の実施の形態3に係る多目標追尾装置は、センサ1と、MHT追尾処理部2と、複数のMFA追尾処理部4と、航跡表示部6と、MHT追尾処理部2によって生成された上位の複数の仮説について各仮説に含まれる航跡の平滑値を複数のMFA追尾処理部4のそれぞれの入力データの第1フレームとするMFA初期解生成部10と、複数のMFA追尾処理部4によって生成された複数の相関解の中から複数の相関解の尤度に応じて1つを選択する解調整部11と、Nフレームデータ記憶部23とが設けられている。また、MFA追尾処理部4は、m個のMFA追尾処理4,4,…,4が設けられている。 10, the multi-target tracking device according to Embodiment 3 of the present invention includes a sensor 1, an MHT tracking processing unit 2, a plurality of MFA tracking processing units 4, a wake display unit 6, and an MHT tracking processing unit 2. MFA initial solution generation unit 10 that uses the smooth value of the wake included in each hypothesis as the first frame of the respective input data of the plurality of MFA tracking processing units 4 and the plurality of MFA tracking A solution adjustment unit 11 that selects one of the plurality of correlation solutions generated by the processing unit 4 according to the likelihood of the plurality of correlation solutions, and an N frame data storage unit 23 are provided. The MFA tracking processing unit 4 is provided with m MFA tracking processes 4 1 , 4 2 ,..., 4 m .

つぎに、この実施の形態3に係る多目標追尾装置の動作について図面を参照しながら説明する。図11は、この発明の実施の形態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.

最初に、ステップ301において、センサ1から観測された観測値が装置に入力され、それを観測値記憶部(図示せず)に蓄積しておく。   First, in step 301, the observation value observed from the sensor 1 is input to the apparatus, and is stored in an observation value storage unit (not shown).

次に、ステップ302において、センサ1の初観測値時刻から最新観測時刻に到る全観測値がN+1フレーム以上か判断し、観測開始時点からNフレーム分が経過するまでは、常にステップ303へ移行し、経過している場合には、ステップ305へ移行する。   Next, in step 302, it is determined whether all observation values from the first observation value time of the sensor 1 to the latest observation time are N + 1 frames or more, and the process always proceeds to step 303 until N frames have elapsed from the observation start time. If it has passed, the process proceeds to step 305.

ステップ303において、MHT追尾処理部2は、従来の多重仮説追尾法による処理と同じ処理を実行する。その処理内容の詳細は、非特許文献1で示されている。   In step 303, the MHT tracking processing unit 2 executes the same processing as the processing based on the conventional multiple hypothesis tracking method. The details of the processing contents are shown in Non-Patent Document 1.

次に、多重仮説追尾法による追尾処理を実行した場合は、ステップ304において、MFA初期解生成部10は、MHT追尾処理部2が生成した仮説群の上位m個について、各仮説に含まれる航跡の平滑値を複数フレーム割当法による追尾処理の入力データの第1フレームとする。   Next, when the tracking process by the multiple hypothesis tracking method is executed, in step 304, the MFA initial solution generation unit 10 includes the tracks included in each hypothesis for the top m hypotheses of the hypothesis group generated by the MHT tracking processing unit 2. Is the first frame of the input data of the tracking process by the multiple frame allocation method.

ステップ305において、MFA追尾処理部4は、Nフレーム分の多次元相関による追尾処理を行う。この処理内容については様々な技術が適用可能であるが、上述した非特許文献2でその一例が示されている。   In step 305, the MFA tracking processing unit 4 performs tracking processing by multidimensional correlation for N frames. Various techniques can be applied to this processing content, but an example is shown in Non-Patent Document 2 described above.

なお、ここでは、m個のMFA追尾処理4,4,…,4が並行して動作する。各MFA追尾処理部に入力され相関処理の対象となるNフレーム分のデータについて、その第1フレームの内容をMHT追尾処理部2が生成したm個の仮説に含まれる航跡の平滑値とする。すなわち、各MFA追尾処理部で第1フレームの平滑値が異なり残りのフレームの観測値が共通であるNフレーム分のデータを処理することになる。よって、相関解も合計m個生成される。 Here, m MFA tracking processes 4 1 , 4 2 ,..., 4 m operate in parallel. For the data of N frames that are input to each MFA tracking processing unit and are subject to correlation processing, the contents of the first frame are set as smooth values of wakes included in the m hypotheses generated by the MHT tracking processing unit 2. That is, each MFA tracking processing unit processes data for N frames in which the smooth value of the first frame is different and the observation values of the remaining frames are common. Therefore, a total of m correlation solutions are generated.

次に、ステップ306において、解調整部11は、MFA追尾処理部4によって生成されたm個の相関解のうち最も尤度の高い相関解を選択する。   Next, in step 306, the solution adjustment unit 11 selects the most likely correlation solution among the m correlation solutions generated by the MFA tracking processing unit 4.

最終のステップ307において、航跡表示部6は、解調整部11によって選択された相関解に含まれる航跡をユーザに表示する。   In the final step 307, the wake display unit 6 displays the wake included in the correlation solution selected by the solution adjustment unit 11 to the user.

以上のように、この実施の形態3によれば、多重仮説追尾法(MHT)による処理によって生成された初期のm個の仮説を基に、複数の複数フレーム割当法(MFA)による処理を並行して実行し、さらに複数フレーム割当法(MFA)による処理によって生成された複数の相関解をその尤度の大小に応じて選択してユーザに表示するため、単一の複数フレーム割当法(MFA)による処理では追尾性能が劣化する場合においても正確な相関解の導出が可能となる。   As described above, according to the third embodiment, processing by a plurality of multiple frame allocation methods (MFA) is performed in parallel based on the initial m hypotheses generated by processing by the multiple hypothesis tracking method (MHT). In addition, a plurality of correlation solutions generated by processing by the multiple frame allocation method (MFA) are selected according to the likelihood and displayed to the user. ) Makes it possible to derive an accurate correlation solution even when the tracking performance deteriorates.

この発明の実施の形態1に係る多目標追尾装置の構成を示す図である。It is a figure 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. ある観測時刻の航跡Tの例と、その後、2観測時刻経過した後で得られる航跡例を示す図である。It is a figure which shows the example of the wake T of a certain observation time, and the wake example obtained after two observation time passes after that. 第1位の仮説H1と第2位の仮説H2の信頼度の推移の例を示す図である。It is a figure which shows the example of transition of the reliability of 1st hypothesis H1 and 2nd hypothesis H2. 第m位の仮説Hmと第m+1位の仮説Hm+1の信頼度の推移の例を示す図である。It is a figure which shows the example of the transition of the reliability of the hypothesis Hm of the m-th rank and the hypothesis Hm + 1 of the m + 1-th rank. 第1位の仮説H1の子孫と第2位の仮説H2の子孫の信頼度の推移の例を示す図である。It is a figure which shows the example of transition of the reliability of the descendant of the 1st hypothesis H1, and the descendant of the 2nd hypothesis H2. この発明の実施の形態1に係る多目標追尾装置の仮説置き換えの例を示す図である。It is a figure which shows the example of hypothesis replacement | exchange of the multi-target tracking apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態2に係る多目標追尾装置の構成を示す図である。It is a figure 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. この発明の実施の形態3に係る多目標追尾装置の構成を示す図である。It is a figure 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. 多目標追尾における相関を説明するための図である。It is a figure for demonstrating the correlation in multitarget tracking. 2次元相関方式における相関決定の例を示す図である。It is a figure which shows the example of the correlation determination in a two-dimensional correlation system. MFA相関方式における相関決定の例を示す図である。It is a figure which shows the example of the correlation determination in a MFA correlation system. 2次元(2観測時刻)の観測値と対応するコスト行列の例を示す図である。It is a figure which shows the example of the cost matrix corresponding to the observation value of two dimensions (2 observation time). コスト行列から得られる解の例と対応する航跡の構成を示す図である。It is a figure which shows the structure of the wake corresponding to the example of the solution obtained from a cost matrix. 制約条件を満たさない相関決定の例を示す図である。It is a figure which shows the example of the correlation determination which does not satisfy | fill a constraint condition. 複数フレーム割当法による従来の多目標追尾装置の構成を示す図である。It is a figure which shows the structure of the conventional multi target tracking apparatus by the multiple frame allocation method. 複数フレーム割当法による従来の多目標追尾装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the conventional multi-target tracking apparatus by the multiple frame allocation method. 複数フレーム割当法による追尾方式と2次元相関による追尾方式を併用する従来の多目標追尾装置の構成を示す図である。It is a figure which shows the structure of the conventional multi-target tracking apparatus which uses together the tracking system by the multiple frame allocation method and the tracking system by two-dimensional correlation.

符号の説明Explanation of symbols

1 センサ、2 MHT追尾処理部、3 MHT仮説監視部、4 MFA追尾処理部、5 MHT仮説設定部、6 航跡表示部、7 解比較部、8 MHT仮説修正部、9 MFA相関解修正部、10 MFA初期解生成部、11 解調整部、21 航跡記憶部、22 仮説記憶部、23 Nフレームデータ記憶部。   1 sensor, 2 MHT tracking processing unit, 3 MHT hypothesis monitoring unit, 4 MFA tracking processing unit, 5 MHT hypothesis setting unit, 6 track display unit, 7 solution comparison unit, 8 MHT hypothesis correction unit, 9 MFA correlation solution correction unit, 10 MFA initial solution generation unit, 11 solution adjustment unit, 21 track storage unit, 22 hypothesis storage unit, 23 N frame data storage unit.

Claims (8)

前観測時刻で多重仮説追尾法による追尾処理を実行する第1の条件、及び前観測時刻で追尾が困難であると判定する第2の条件を満たす場合には、センサからの複数のフレームに跨る観測値を用いて多次元の相関を取りながら、複数フレーム割当法による処理で相関解を生成するMFA追尾処理部と、
前記第1及び第2の条件の何れかの条件を満たさない場合には、センサからの最新観測値と既存追尾航跡の2次元の相関を取り、複数の仮説を生成しながら新たな航跡を生成するMHT追尾処理部と、
前記MHT追尾処理部による追尾処理実行中に、生成された仮説の状況から追尾が困難であるかどうかを判定するMHT仮説監視部と、
前記MFA追尾処理部で生成された相関解を前記MHT追尾処理部の仮説として置き換えるMHT仮説設定部と
を備えたことを特徴とする多目標追尾装置。
When the first condition for executing the tracking process by the multiple hypothesis tracking method at the previous observation time and the second condition for determining that the tracking is difficult at the previous observation time are satisfied, a plurality of frames from the sensor are straddled. An MFA tracking processing unit that generates a correlation solution by processing using a multi-frame allocation method while obtaining multi-dimensional correlation using observation values;
If either of the first and second conditions is not met, a two-dimensional correlation between the latest observation value from the sensor and the existing tracking track is taken, and a new track is generated while generating a plurality of hypotheses. An MHT tracking processing unit,
An MHT hypothesis monitoring unit that determines whether tracking is difficult from the status of the generated hypothesis during execution of the tracking process by the MHT tracking processing unit;
A multi-target tracking device comprising: an MHT hypothesis setting unit that replaces the correlation solution generated by the MFA tracking processing unit as a hypothesis of the MHT tracking processing unit.
前記MHT仮説監視部は、前記MHT追尾処理部が生成した第1位の仮説と第2位の仮説の信頼度の差によって追尾が困難であるかどうかを判定する
ことを特徴とする請求項1記載の多目標追尾装置。
The MHT hypothesis monitoring unit determines whether tracking is difficult based on a difference in reliability between the first hypothesis and the second hypothesis generated by the MHT tracking processing unit. The multi-target tracking device described.
前記MHT仮説監視部は、前記MHT追尾処理部の仮説選択方法が信頼度に応じた上位m(自然数)個の仮説を残すという方法を採用していることを前提として、前記MHT追尾処理部が生成した第m位の仮説と第m+1位の仮説の信頼度の差によって追尾が困難であるかどうかを判定する
ことを特徴とする請求項1記載の多目標追尾装置。
The MHT hypothesis monitoring unit assumes that the MHT tracking processing unit adopts a method in which the hypothesis selection method of the MHT tracking processing unit retains the top m (natural number) hypotheses according to reliability. The multi-target tracking device according to claim 1, wherein it is determined whether or not tracking is difficult based on a difference in reliability between the generated m-th hypothesis and m + 1-th hypothesis.
前記MHT仮説監視部は、前記MHT追尾処理部が生成したある観測時刻における第1位の仮説と第2位の仮説の子孫の仮説の順位の変動の有無によって追尾が困難であるかどうかを判定する
ことを特徴とする請求項1記載の多目標追尾装置。
The MHT hypothesis monitoring unit determines whether or not tracking is difficult depending on whether or not the rank of the hypothesis of the first hypothesis and the hypothesis of the second hypothesis at a certain observation time generated by the MHT tracking processing unit is changed. The multi-target tracking device according to claim 1, wherein:
センサからの最新観測値と既存追尾航跡の2次元の相関を取り、複数の仮説を生成しながら新たな航跡を生成するMHT追尾処理部と、
前記MHT追尾処理部と並行して、センサからの複数のフレームに跨る観測値を用いて多次元の相関を取りながら、複数フレーム割当法による処理で相関解を生成するMFA追尾処理部と、
前記MHT追尾処理部によって生成された第1位の仮説と前記MFA追尾処理部によって生成された相関解の尤度を比較する解比較部と、
前記解比較部による比較結果として、前記MFA追尾処理部の相関解の尤度が大きい場合には、前記MHT追尾処理部の第1位の仮説を前記MFA追尾処理部の相関解で置き換えるMHT仮説修正部と、
前記解比較部による比較結果として、前記MHT追尾処理部の第1位の仮説の尤度が大きい場合には、前記MFA追尾処理部の相関解を前記MHT追尾処理部の第1位の仮説で置き換えるMFA相関解修正部と
を備えたことを特徴とする多目標追尾装置。
A MHT tracking processing unit that takes a two-dimensional correlation between the latest observation value from the sensor and the existing tracking track, and generates a new track while generating a plurality of hypotheses;
In parallel with the MHT tracking processing unit, an MFA tracking processing unit that generates a correlation solution by processing using a multiple frame allocation method while obtaining a multidimensional correlation using observation values across a plurality of frames from a sensor;
A solution comparison unit that compares the first hypothesis generated by the MHT tracking processing unit with the likelihood of the correlation solution generated by the MFA tracking processing unit;
When the likelihood of the correlation solution of the MFA tracking processing unit is large as a comparison result by the solution comparison unit, the MHT hypothesis that replaces the first hypothesis of the MHT tracking processing unit with the correlation solution of the MFA tracking processing unit Correction part,
If the likelihood of the first hypothesis of the MHT tracking processing unit is large as a comparison result by the solution comparison unit, the correlation solution of the MFA tracking processing unit is expressed as the first hypothesis of the MHT tracking processing unit. A multi-target tracking device comprising: a replacement MFA correlation solution correction unit.
センサからの最新観測値と既存追尾航跡の2次元の相関を取り、複数の仮説を生成しながら新たな航跡を生成するMHT追尾処理部と、
センサからの複数のフレームに跨る観測値を用いて多次元の相関を取りながら、複数フレーム割当法による処理で相関解をそれぞれ生成する複数のMFA追尾処理部と、
前記MHT追尾処理部によって生成された上位の複数の仮説について各仮説に含まれる航跡の平滑値を前記複数のMFA追尾処理部のそれぞれの入力データの第1フレームとするMFA初期解生成部と、
前記複数のMFA追尾処理部によって生成された複数の相関解の中から前記複数の相関解の尤度に応じて1つを選択する解調整部と
を備えたことを特徴とする多目標追尾装置。
A MHT tracking processing unit that takes a two-dimensional correlation between the latest observation value from the sensor and the existing tracking track, and generates a new track while generating a plurality of hypotheses;
A plurality of MFA tracking processing units each generating a correlation solution by processing by a multi-frame allocation method while obtaining multi-dimensional correlation using observation values across a plurality of frames from a sensor;
An MFA initial solution generator that uses a smooth value of a wake included in each hypothesis for a plurality of upper hypotheses generated by the MHT tracking processor as a first frame of input data of each of the plurality of MFA tracking processors;
A multi-target tracking apparatus comprising: a solution adjusting unit that selects one of the plurality of correlation solutions generated by the plurality of MFA tracking processing units according to the likelihood of the plurality of correlation solutions. .
センサからの最新観測値と既存追尾航跡の2次元の相関を取り、複数の仮説を生成しながら新たな航跡を生成するMHT追尾処理部と、
前記MHT追尾処理部による追尾処理実行中に、生成された仮説の状況から追尾が困難であるかどうかを判定するMHT仮説監視部とを備え、
前記MHT仮説監視部により追尾が困難でないと判断された場合には、前記MHT追尾処理部による処理を継続する
ことを特徴とする多目標追尾装置。
A MHT tracking processing unit that takes a two-dimensional correlation between the latest observation value from the sensor and the existing tracking track, and generates a new track while generating a plurality of hypotheses;
An MHT hypothesis monitoring unit that determines whether tracking is difficult from the generated hypothesis status during the execution of the tracking process by the MHT tracking processing unit;
When the MHT hypothesis monitoring unit determines that tracking is not difficult, the multi-target tracking device continues processing by the MHT tracking processing unit.
前観測時刻で多重仮説追尾法による追尾処理を実行する第1の条件、及び前観測時刻で前記MHT仮説監視部により追尾が困難であると判定される第2の条件を満たす場合には、センサからの複数のフレームに跨る観測値を用いて多次元の相関を取りながら、複数フレーム割当法による処理で相関解を生成するMFA追尾処理部と、
前記MFA追尾処理部で生成された相関解を前記MHT追尾処理部の仮説として置き換えるMHT仮説設定部とをさらに備えた
ことを特徴とする請求項7記載の多目標追尾装置。
When the first condition for executing the tracking process by the multiple hypothesis tracking method at the previous observation time and the second condition for determining that tracking is difficult by the MHT hypothesis monitoring unit at the previous observation time are satisfied, An MFA tracking processing unit that generates a correlation solution by a process using a multi-frame allocation method while obtaining a multi-dimensional correlation using observation values over a plurality of frames from
The multi-target tracking device according to claim 7, further comprising an MHT hypothesis setting unit that replaces the correlation solution generated by the MFA tracking processing unit as a hypothesis of the MHT tracking processing unit.
JP2008134327A 2008-05-22 2008-05-22 Multi-target tracking device Active JP5197154B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008134327A JP5197154B2 (en) 2008-05-22 2008-05-22 Multi-target tracking device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008134327A JP5197154B2 (en) 2008-05-22 2008-05-22 Multi-target tracking device

Publications (2)

Publication Number Publication Date
JP2009281865A true JP2009281865A (en) 2009-12-03
JP5197154B2 JP5197154B2 (en) 2013-05-15

Family

ID=41452469

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008134327A Active JP5197154B2 (en) 2008-05-22 2008-05-22 Multi-target tracking device

Country Status (1)

Country Link
JP (1) JP5197154B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011145096A (en) * 2010-01-12 2011-07-28 Mitsubishi Electric Corp Multi-target tracking apparatus
JP2012112807A (en) * 2010-11-25 2012-06-14 Mitsubishi Electric Corp Multi-target tracking device
JP2012127890A (en) * 2010-12-17 2012-07-05 Mitsubishi Electric Corp Estimation device, computer program, and estimation method
JP2015031656A (en) * 2013-08-06 2015-02-16 株式会社Ihi Method and device for estimating object position
JP2017067624A (en) * 2015-09-30 2017-04-06 三菱電機株式会社 Target trailing device and radar device

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103759732B (en) * 2014-01-14 2016-06-22 北京航空航天大学 A kind of centralized multisensor multiple hypotheis tracking method of angle information auxiliary
RU2657005C1 (en) * 2017-05-05 2018-06-08 Акционерное общество "НИИ измерительных приборов - Новосибирский завод имени Коминтерна" (АО "НПО НИИИП-НЗиК") Method of target tracking by surveillance radar station (options)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11118922A (en) * 1997-10-14 1999-04-30 Mitsubishi Electric Corp Target tracking apparatus
JP2003057337A (en) * 2001-08-20 2003-02-26 Nec Corp Tracking processing device, tracking processing method, and tracking processing program
JP2007212244A (en) * 2006-02-08 2007-08-23 Mitsubishi Electric Corp Multiple target tracking device
JP2009162615A (en) * 2008-01-07 2009-07-23 Mitsubishi Electric Corp Track allocation device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11118922A (en) * 1997-10-14 1999-04-30 Mitsubishi Electric Corp Target tracking apparatus
JP2003057337A (en) * 2001-08-20 2003-02-26 Nec Corp Tracking processing device, tracking processing method, and tracking processing program
JP2007212244A (en) * 2006-02-08 2007-08-23 Mitsubishi Electric Corp Multiple target tracking device
JP2009162615A (en) * 2008-01-07 2009-07-23 Mitsubishi Electric Corp Track allocation device
JP5046962B2 (en) * 2008-01-07 2012-10-10 三菱電機株式会社 Wake assignment device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN6012042583; Samuel Blackman, et al.: "Design and Analysis of Modern Tracking Systems" , 1999, p.403-423, Attech House *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011145096A (en) * 2010-01-12 2011-07-28 Mitsubishi Electric Corp Multi-target tracking apparatus
JP2012112807A (en) * 2010-11-25 2012-06-14 Mitsubishi Electric Corp Multi-target tracking device
JP2012127890A (en) * 2010-12-17 2012-07-05 Mitsubishi Electric Corp Estimation device, computer program, and estimation method
JP2015031656A (en) * 2013-08-06 2015-02-16 株式会社Ihi Method and device for estimating object position
JP2017067624A (en) * 2015-09-30 2017-04-06 三菱電機株式会社 Target trailing device and radar device

Also Published As

Publication number Publication date
JP5197154B2 (en) 2013-05-15

Similar Documents

Publication Publication Date Title
JP5197154B2 (en) Multi-target tracking device
Bouwmans et al. Robust PCA via principal component pursuit: A review for a comparative evaluation in video surveillance
JP6610278B2 (en) Machine learning apparatus, machine learning method, and machine learning program
Norena et al. Bayesian analysis of inflation. III. Slow roll reconstruction using model selection
US10108884B2 (en) Learning user preferences for photo adjustments
WO2020119318A1 (en) Self-adaptive selection and design method for convolutional-layer hardware accelerator
US11961298B2 (en) Memory-guided video object detection
KR101852116B1 (en) Denoiser, and control method thereof
US11829869B2 (en) Multiple task transfer learning
WO2017197330A1 (en) Two-stage training of a spoken dialogue system
JP2018134199A5 (en)
US8290884B2 (en) Method for approximating user task representations by document-usage clustering
JP5014637B2 (en) Multi-target tracking device
JP5034469B2 (en) Information processing apparatus, information processing method, and program
Srivastava et al. An asynchronous distributed expectation maximization algorithm for massive data: the DEM algorithm
Mukherjee et al. Distribution-dependent and time-uniform bounds for piecewise iid bandits
WO2013126336A1 (en) Active acquisition of priviledged information
JP2020038485A (en) Learning apparatus, determining apparatus, learning method, determining method, learning program and determining program
JP2014110020A (en) Image processor, image processing method and image processing program
JP5161172B2 (en) Frame rate conversion method, frame rate conversion device, and frame rate conversion program
Castro et al. Compressed sensing vs. active learning
Darkhovsky et al. A method of generating random vectors with a given probability density function
CN106815858A (en) A kind of Extracting of Moving Object and device
Hwang et al. Overcoming Overconfidence for Active Learning
CN105651284B (en) The method and device of raising experience navigation interior joint efficiency of selection

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110401

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120725

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120814

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20121015

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: 20130108

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130205

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

Free format text: PAYMENT UNTIL: 20160215

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5197154

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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