JP2018136337A - Tracking processing device and tracking processing method - Google Patents
Tracking processing device and tracking processing method Download PDFInfo
- Publication number
- JP2018136337A JP2018136337A JP2018084387A JP2018084387A JP2018136337A JP 2018136337 A JP2018136337 A JP 2018136337A JP 2018084387 A JP2018084387 A JP 2018084387A JP 2018084387 A JP2018084387 A JP 2018084387A JP 2018136337 A JP2018136337 A JP 2018136337A
- Authority
- JP
- Japan
- Prior art keywords
- tracking
- target
- congestion
- congestion degree
- echo
- 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
Links
- 238000012545 processing Methods 0.000 title claims abstract description 110
- 238000003672 processing method Methods 0.000 title claims description 5
- 238000004364 calculation method Methods 0.000 claims description 45
- 238000000034 method Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 20
- 230000033001 locomotion Effects 0.000 description 35
- 238000001514 detection method Methods 0.000 description 21
- 238000012546 transfer Methods 0.000 description 19
- 238000009499 grossing Methods 0.000 description 18
- 239000013256 coordination polymer Substances 0.000 description 10
- 230000004048 modification Effects 0.000 description 8
- 238000012986 modification Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 239000000284 extract Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
- G01S13/726—Multiple target tracking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
- G01S13/44—Monopulse radar, i.e. simultaneous lobing
- G01S13/449—Combined with MTI or Doppler processing circuits
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/04—Display arrangements
- G01S7/06—Cathode-ray tube displays or other two dimensional or three-dimensional displays
- G01S7/22—Producing cursor lines and indicia by electronic means
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G3/00—Traffic control systems for marine craft
- G08G3/02—Anti-collision systems
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Ocean & Marine Engineering (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
本発明は、追尾対象を追尾する追尾処理装置及び追尾処理方法に関する。 The present invention relates to a tracking processing device and a tracking processing method for tracking a tracking target.
従来から知られている追尾処理装置として、例えば特許文献1に開示される追尾処理装置では、追尾対象としてのターゲットの運動状態が推定され、当該運動状態に基づいてターゲットが追尾される。これにより、複数の物標の中から、追尾したい物標を精度よく検出することができる。
As a conventionally known tracking processing device, for example, in the tracking processing device disclosed in
ところで、上述のように追尾対象となるターゲットを追尾する場合、ターゲットの周囲の環境によっては、ターゲットを正確に追尾できない場合がある。 By the way, when the target to be tracked is tracked as described above, the target may not be accurately tracked depending on the environment around the target.
本発明は、上記課題を解決するためのものであり、その目的は、周囲の環境によらず、追尾対象を正確に追尾することである。 The present invention is to solve the above-described problems, and an object thereof is to accurately track a tracking target regardless of the surrounding environment.
(1)上記課題を解決するため、本発明のある局面に係る追尾処理装置は、追尾対象を追尾する処理を行う追尾処理部と、前記追尾対象の予測位置が含まれるエリア内に存在する物標の混雑の度合いである混雑度を算出する混雑度算出部と、を備え、前記追尾処理部は、前記混雑度算出部によって算出された前記混雑度の値が低いほど前記追尾対象への追従性が高くなるようにゲインを設定することを特徴とする。 (1) In order to solve the above problem, a tracking processing device according to an aspect of the present invention includes a tracking processing unit that performs processing for tracking a tracking target, and an object that is present in an area including the predicted position of the tracking target. A congestion degree calculation unit that calculates a degree of congestion that is a degree of congestion of the target, and the tracking processing unit follows the tracking target as the value of the congestion degree calculated by the congestion degree calculation unit decreases. It is characterized in that the gain is set so as to increase the performance.
(2)好ましくは、前記追尾処理部は、前記混雑度算出部によって算出された前記混雑度の値が低いほど、段階的に前記追尾対象への追従性が高くなるようにゲインを設定することを特徴とする。 (2) Preferably, the tracking processing unit sets the gain so that the followability to the tracking target increases stepwise as the value of the congestion level calculated by the congestion level calculation unit decreases. It is characterized by.
(3)好ましくは、前記混雑度算出部は、前記エリア内に存在する前記物標の数としての物標数に基づいて前記混雑度を算出する。 (3) Preferably, the said congestion degree calculation part calculates the said congestion degree based on the target number as the number of the said target which exists in the said area.
(4)更に好ましくは、前記追尾処理装置は、前記エリアが分割された複数のセルのそれぞれに存在する前記物標数をカウントし、前記複数のセル毎にカウントされた前記物標数を、各前記セルに対応させて記憶するエコー分布生成部を更に備え、前記混雑度算出部は、前記追尾対象の前記予測位置が含まれる前記セル内に存在する前記物標の前記混雑度を算出する。 (4) More preferably, the tracking processing device counts the target number that exists in each of a plurality of cells into which the area is divided, and the target number counted for each of the plurality of cells, Further comprising an echo distribution generation unit that stores data corresponding to each cell, and the congestion level calculation unit calculates the congestion level of the target existing in the cell including the predicted position of the tracking target. .
(5)好ましくは、前記混雑度算出部は、複数のタイミングでの前記物標数に基づいて前記混雑度を算出する。 (5) Preferably, the said congestion degree calculation part calculates the said congestion degree based on the said target number in several timing.
(6)好ましくは、前記混雑度算出部は、複数のタイミングでの前記混雑度に基づいて平滑混雑度を算出し、前記追尾処理部は、前記混雑度算出部によって算出された前記平滑混雑度の値に応じて、前記追尾対象を追尾する処理を行う。 (6) Preferably, the congestion level calculation unit calculates a smooth congestion level based on the congestion levels at a plurality of timings, and the tracking processing unit calculates the smooth congestion level calculated by the congestion level calculation unit. In accordance with the value of, the tracking target is tracked.
(7)上記課題を解決するため、本発明のある局面に係る追尾処理方法は、追尾対象を追尾する処理を行うステップと、前記追尾対象の予測位置が含まれるエリア内に存在する物標の混雑の度合いである混雑度を算出するステップと、を含み、前記追尾する処理を行うステップは、前記混雑度を算出するステップによって算出された前記混雑度の値が低いほど前記追尾対象への追従性が高くなるようにゲインを設定することを特徴とする。 (7) In order to solve the above problem, a tracking processing method according to an aspect of the present invention includes a step of performing processing of tracking a tracking target, and a target existing in an area including the predicted position of the tracking target. Calculating the degree of congestion that is a degree of congestion, and performing the tracking process includes following the tracking target as the value of the degree of congestion calculated by the step of calculating the degree of congestion is lower. It is characterized in that the gain is set so as to increase the performance.
本発明によれば、周囲の環境によらず、追尾対象を正確に追尾することができる。 According to the present invention, it is possible to accurately track a tracking target regardless of the surrounding environment.
以下、本発明に係る追尾処理装置3の実施形態について図面を参照しつつ説明する。本発明は、追尾対象として選択された物標を追尾する、追尾処理装置として広く適用することができる。以下では、追尾対象として設定されたターゲットを「追尾物標」という。また、以下では、図中同一または相当部分には、同一符号を付してその説明は繰り返さない。
Hereinafter, an embodiment of a
図1は、本発明の実施形態に係る追尾処理装置3を含む、レーダ装置1のブロック図である。本実施形態のレーダ装置1は、例えば、漁船等の船舶に備えられる舶用レーダである。レーダ装置1は、主に他船等の物標の探知に用いられる。また、レーダ装置1は、追尾物標として選択された物標を追尾することが可能に構成されている。レーダ装置1は、複数の追尾物標を、同時に追尾可能に構成されている。レーダ装置1は、追尾物標の運動状態を推定するように構成されている。本実施形態では、レーダ装置1は、上記運動状態として、追尾物標の平滑速度を算出する。平滑速度とは、追尾物標について推定される、進行方向及び進行速度を示すベクトルである。レーダ装置1は、追尾物標の平滑速度を、画面に表示する。尚、以下では、レーダ装置1が備えられている船舶を「自船」という。
FIG. 1 is a block diagram of a
図1に示すように、レーダ装置1は、アンテナユニット2と、追尾処理装置3と、表示器4と、を備えている。
As shown in FIG. 1, the
アンテナユニット2は、アンテナ5と、受信部6と、A/D変換部7と、を含んでいる。
The
アンテナ5は、指向性の強いパルス状電波を送信可能なレーダアンテナである。また、アンテナ5は、物標からの反射波であるエコー信号を受信するように構成されている。即
ち、物標のエコー信号は、アンテナ5からの送信信号に対する、物標での反射波である。レーダ装置1は、パルス状電波を送信してからエコー信号を受信するまでの時間を測定する。これにより、レーダ装置1は、物標までの距離rを検出することができる。アンテナ5は、水平面上で360°回転可能に構成されている。アンテナ5は、パルス状電波の送信方向を変えながら(アンテナ角度を変えながら)、電波の送受信を繰り返し行うように構成されている。以上の構成で、レーダ装置1は、自船周囲の平面上の物標を、360°にわたり探知することができる。
The antenna 5 is a radar antenna capable of transmitting a pulsed radio wave having strong directivity. The antenna 5 is configured to receive an echo signal that is a reflected wave from a target. That is, the echo signal of the target is a reflected wave at the target with respect to the transmission signal from the antenna 5. The
なお、以下の説明では、パルス状電波を送信してから次のパルス状電波を送信するまでの動作を「スイープ」という。また、電波の送受信を行いながらアンテナを360°回転させる動作を「スキャン」と呼ぶ。以下では、あるスキャンのことを、「nスキャン」といい、nスキャンから1個前のスキャンのことを、「n−1スキャン」という。尚、nは自然数である。また、以下では、動作について特に説明なき場合、nスキャン時点におけるレーダ装置1の動作を説明する。
In the following description, an operation from transmission of a pulsed radio wave to transmission of the next pulsed radio wave is referred to as “sweep”. The operation of rotating the antenna 360 ° while transmitting / receiving radio waves is called “scan”. In the following, a certain scan is referred to as an “n scan”, and a scan immediately before the n scan is referred to as an “n−1 scan”. Note that n is a natural number. In the following, the operation of the
受信部6は、アンテナ5で受信したエコー信号を検波して増幅する。受信部6は、増幅したエコー信号を、A/D変換部7へ出力する。A/D変換部7は、アナログ形式のエコー信号をサンプリングし、複数ビットからなるデジタルデータ(エコーデータ)に変換する。ここで、上記エコーデータは、アンテナ5が受信したエコー信号の強度(信号レベル)を特定するデータを含んでいる。A/D変換部7は、エコーデータを、追尾処理装置3へ出力する。
The
追尾処理装置3は、複数の物標の中から選択された物標を追尾物標として特定し、当該追尾物標を観測して得られた観測位置に追尾フィルタ処理を施すことで、追尾物標を追尾する追尾処理を行うように構成されている。より具体的には、追尾処理装置3は、追尾物標の平滑速度、及び追尾物標の推定位置(平滑位置)等を算出するように構成されている。
The
追尾処理装置3は、CPU、RAM及びROM(図示せず)等を含むハードウェアを用いて構成されている。また、追尾処理装置3は、ROMに記憶された追尾処理プログラムを含む、ソフトウェアを用いて構成されている。
The
上記追尾処理プログラムは、本発明に係る追尾処理方法を、追尾処理装置3に実行させるためのプログラムである。上記ハードウェアとソフトウェアとは、協働して動作するように構成されている。これにより、追尾処理装置3を、信号処理部9、エコー検出部10、及び追尾処理部11等として機能させることができる。追尾処理装置3は、1スキャン毎に、以下に説明する処理を行うように構成されている。
The tracking processing program is a program for causing the
追尾処理装置3は、信号処理部9と、エコー検出部(検出部)10と、追尾処理部11と、混雑度算出処理部20と、を有している。
The
信号処理部9は、フィルタ処理等を施すことにより、エコーデータに含まれる干渉成分と、不要な波形データとを除去する。また、信号処理部9は、物標エコー像に関するエコーデータの特徴情報の検出を行うように構成されている。信号処理部9は、処理したエコーデータをエコー検出部10へ出力する。
The
エコー検出部10は、物標エコー像の検出と、物標エコー像に関するエコーデータの特徴情報の検出とを行うように構成されている。即ち、エコー検出部10は、物標エコー像検出部と、特徴情報抽出部とを含んでいる。信号処理部9及びエコー検出部10によって、物標の特徴情報を検出するための検出部が構成されている。
The
エコー検出部10は、信号処理部9からエコーデータを読み出すときの読出しアドレスに基づいて、当該エコーデータに対応する位置までの距離rを求める。また、アンテナ5からエコー検出部10へは、当該アンテナ5が現在どの方向を向いているか(アンテナ角度θ)を示すデータが出力されている。以上の構成で、エコー検出部10は、エコーデータを読み出す際には、当該エコーデータに対応する位置を、距離rとアンテナ角度θとの極座標で取得することができる。
The
エコー検出部10は、エコーデータに対応する位置に物標が存在するか否かを検出するように構成されている。エコー検出部10は、例えば、エコーデータに対応する位置の信号レベル、即ち、信号強度を判別する。エコー検出部10は、信号レベルが所定のしきいレベル値以上である位置には、物標が存在していると判別する。
The
次いで、エコー検出部10は、物標が存在している範囲を検出する。エコー検出部10は、例えば、物標が存在している一まとまりの領域を、物標エコー像が存在している領域として検出する。このようにして、エコー検出部10は、エコーデータを基に、物標エコー像を検出する。この物標エコー像の外郭形状は、物標の外郭形状と略合致する。但し、エコーデータに含まれるノイズ等に起因して、この物標エコー像の外郭形状と、物標の外郭形状とは、わずかに異なる。次に、エコー検出部10は、エコーデータを用いて、物標エコー像に関連する特徴情報を抽出する。
Next, the
図2は、自船100と、物標エコー像120との関係を説明するための模式的な平面図である。図2では、物標エコー像120は、矩形の像として例示されている。また、図2では、物標エコー像120によって特定される物標130が示されている。図2では、物標130の外郭形状は、物標エコー像120と合致した状態で表示されている。
FIG. 2 is a schematic plan view for explaining the relationship between the
図1及び図2に示すように、極座標系では、自船100の位置としての自船位置M1を基準として、自船位置M1からの直線距離が、距離rとして示され、自船位置M1周りの角度が、角度θとして示される。本実施形態では、自船位置M1は、アンテナ5の位置に相当する。エコー検出部10は、物標エコー像120の代表点Pの抽出に際しては、自船位置M1を中心とする、リング状部分の一部形状の像110を用いる。この像110は、第1直線111、第2直線112、第1円弧113、及び第2円弧114によって囲まれた領域の像である。
As shown in FIGS. 1 and 2, in the polar coordinate system, a straight line distance from the ship position M1 is indicated as a distance r with respect to the ship position M1 as the position of the
第1直線111は、物標エコー像120の後縁120aのうち自船位置M1に最も近い点と、自船位置M1と、を通る直線である。第2直線112は、物標エコー像120の前縁120bのうち自船位置M1に最も近い点と、自船位置M1と、を通る直線である。第1円弧113は、物標エコー像120のうち自船位置M1から最も近い部分120cを通る円弧である。第1円弧113の曲率中心点は、自船位置M1である。第2円弧114は、物標エコー像120のうち自船位置M1から最も遠い部分120dを通る円弧である。第2円弧114は、第1円弧113と同心である。
The first
図3は、物標エコー像120(物標130)に関して抽出されるデータを説明するためのデータ一覧表である。図2及び図3に示すように、本実施形態では、信号処理部9及びエコー検出部10は、協働して、物標エコー像120(物標130)について、下記12個のデータを、特徴情報データとして抽出する。即ち、エコー検出部10は、エコーデータ、及び物標エコー像120の画像データを基に、下記12個のテキストデータを抽出する。本実施形態では、12個のテキストデータとは、フラグのデータ201と、距離rpのデータ202と、終了角度θeのデータ203と、角度幅θwのデータ204と、最前縁距離rnのデータ205と、最後縁距離rfのデータ206と、面積ar(形状情報)
のデータ207と、代表点Pの座標データ208と、エコーレベルecのデータ209と、隣接距離(物標の周囲の状態を特定する情報)adのデータ210と、ドップラーシフト量dsのデータ211と、時刻tmのデータ212と、である。
FIG. 3 is a data list for explaining data extracted with respect to the target echo image 120 (target 130). As shown in FIGS. 2 and 3, in the present embodiment, the
上記のフラグは、例えば、物標エコー像120について、所定の状態であるか否か、等を特定するフラグである。このフラグは、「1」又は「0」等に設定されるように構成されている。
The above flag is a flag for specifying whether or not the
距離rpは、自船位置M1から物標エコー像120の代表点Pまでの直線距離である。本実施形態では、代表点Pは、像110の中心点である。終了角度θeは、物標エコー像120の検出が終わった時点における、前述のアンテナ角度θである。角度幅θwは、物標エコー像120について、自船位置M1回りの角度方向の幅である。角度幅θwは、第1直線111と第2直線112とがなす角度でもある。最前縁距離rnは、物標エコー像120の部分120cと、自船位置M1との距離である。最後縁距離rfは、物標エコー像120の部分120dと、自船位置M1との距離である。面積arは、リング状部分の一部形状の像110における面積であり、本実施形態では、物標エコー像120の面積として扱われる。
The distance rp is a linear distance from the ship position M1 to the representative point P of the
エコーレベルecは、物標エコー像120を特定するエコー信号の、強度を示している。この強度は、物標エコー像120を特定するエコー信号のピーク強度であってもよいし、当該エコー信号の強度の平均値であってもよい。本実施形態では、隣接距離adは、例えば、隣接する2つの物標エコー像120間の距離である。ドップラーシフト量dsは、例えば、アンテナ5から放射されたパルス信号の周波数と、物標エコー像120によって特定される物標130で反射したエコー信号の周波数と、の差である。ドップラーシフト量dsを基に、物標エコー像120によって特定される物標130と、自船100との相対速度を求めることが可能である。時刻tmは、物標エコー像120を検出した時点の時刻である。尚、エコー像120のデータは、予備のデータ領域を含んでいてもよい。本実施形態では、信号処理部9及びエコー検出部10は、各物標エコー像120について、上記12個の特徴情報を抽出する。これら12個の特徴情報は、何れも、数値によって表される情報である。
The echo level ec indicates the intensity of an echo signal that identifies the
エコー検出部10で検出された複数の物標エコー像120の一例を、図4に示している。図4は、エコー検出部10で検出された、物標エコー像120を示す模式的な平面図である。図4では、一例として、nスキャン時点における4つの物標エコー像120(121,122,123,124)を示している。図4では、物標エコー像120(121,122,123,124)の形状は、それぞれ、物標130(131,132,133,134)の形状と合致している。
An example of a plurality of target echo
物標エコー像121によって特定される物標131,物標エコー像122によって特定される物標132,及び物標エコー像123によって特定される物標133は、例えば、小型船舶である。物標エコー像124によって特定される物標134は、例えば、大型船舶である。エコー検出部10は、nスキャン時点において、代表点Pとして、物標エコー像121の代表点P1(n)と、物標エコー像122の代表点P2(n)と、物標エコー像123の代表点P3(n)と、物標エコー像124の代表点P4(n)と、を検出している。以下では、物標131が、追尾物標140である場合を例に説明する。
The target 131 specified by the target echo image 121, the target 132 specified by the target echo image 122, and the target 133 specified by the target echo image 123 are, for example, small ships. The target 134 specified by the target echo image 124 is, for example, a large ship. The
図1〜図4に示すように、エコー検出部10は、各物標エコー像120についての特徴情報データ201〜212を、追尾処理部11及び混雑度算出処理部20へ出力する。
As shown in FIGS. 1 to 4, the
[混雑度算出処理部の構成]
図5は、混雑度算出処理部20の動作を説明するためのフローチャートである。以下では、図1及び図5等を参照しながら、混雑度算出処理部20の構成及び動作を説明する。
[Configuration of Congestion Calculation Processing Unit]
FIG. 5 is a flowchart for explaining the operation of the congestion degree
混雑度算出処理部20は、追尾物標140の予測位置Xp(n)周辺に存在する物標の混雑の度合い(混雑度)を算出するように構成されている。混雑度算出処理部20は、エコー分布生成部21と、混雑度算出部22と、を有している。混雑度算出処理部20には、詳しくは後述するように追尾処理部11によって算出された追尾物標のnスキャン時における予測位置Xp(n)が入力される。
The congestion degree
図6は、エコー分布生成部21によって生成されたエコー分布を模式的に示す図である。以下では、図6に示すxy座標上における位置が(x,y)であるセルBにnスキャン時に含まれる物標の数(物標数)を、N(x,y,n)として説明する。例えば、図6を参照して、N(1,5,n)=3である。なお、図6では、カウントされた物標数を、黒のドットで図示している。また、以下では、座標系がxy座標である例を挙げて説明するが、これに限らず、rθ座標系を用いることもできる。
FIG. 6 is a diagram schematically illustrating the echo distribution generated by the echo
エコー分布生成部21は、エコー検出部10によって検出された各物標130の分布を生成する。具体的には、エコー分布生成部21では、図6を参照して、海上における所定のエリアAが格子状に分割された複数のセルBのそれぞれに存在する物標130の数が、セルB毎にカウントされる(図5のステップS1)。このとき、各セルBに対応して平滑物標数(詳しくは後述する)が記憶されていない場合(ステップS2のNo)、エコー分布生成部21は、セルB毎にカウントされた物標130の物標数を、平滑物標数として各セルBに対応させて記憶する(ステップS3)。なお、エリアAとしては、例えば一例として、一辺が64マイルの正方形状の領域を挙げることができる。また、セルBとしては、例えば一例として、一辺が0.5マイルの正方形状の領域を挙げることができる。なお、本実施形態では、セルB毎に物標数をカウントする例を挙げて説明しているが、これに限らず、例えば一例として、各セルBに含まれる物標の面積を合計した値を、上述した物標数の代わりに用いてもよい。
The echo
一方、各セルBに対応して平滑物標数が記憶されている場合(ステップS2のYes)、エコー分布生成部21は、各セルBに対応して記憶されている平滑物標数Ns(x,y,n−1)と、直近の物標数N(x,y,n)とに基づいて、平滑物標数Ns(x,y,n)を算出する(ステップS4)。具体的には、エコー分布生成部21は、以下で示す式(1)に基づいて、平滑物標数Ns(x,y,n)を算出する。
On the other hand, when the smooth target characteristic is stored corresponding to each cell B (Yes in step S2), the echo
Ns(x,y,n)=a1・N(x,y,n)+(1−a1)・Ns(x,y,n−1)…(1) Ns (x, y, n) = a 1 · N (x, y, n) + (1−a 1 ) · Ns (x, y, n−1) (1)
なお、a1は、0<a1≦1で表される係数であり、例えば一例として0.9である。 Here, a 1 is a coefficient represented by 0 <a 1 ≦ 1, for example 0.9 as an example.
そして、エコー分布生成部21は、それまでに各セルBに対応して記憶されていた平滑物標数Ns(x,y,n−1)を、直近で算出した平滑物標数Ns(x,y,t)に置き換えて記憶する(ステップS5)。
The echo
混雑度算出部22は、追尾物標140の予測位置Xp(n)が含まれるセルBの平滑物標数Ns(x,y,n)を用い、混雑度C(n)を算出する(ステップS6)。混雑度C(n)は、平滑物標数Ns(x,y,n)を変数とする任意の関数、すなわち、C(n)=f(Ns(x,y,n))と表すことができる。この式の一例として、例えば、C(n)=Ns(x,y,n)、を挙げることができる。そして、混雑度算出部22は、過去に算出された上記セルBの混雑度であって該セルBに対応して記憶されている平滑混雑度Cs(n−1)と、直近で算出した混雑度C(n)とに基づいて、平滑混雑度Cs(n)を算出する(ステップS7)。具体的には、混雑度算出部22は、以下の式(2)に基づいて、平滑混雑度Cs(n)を算出する。そして、混雑度算出部22は、それまで記憶されていた平滑混雑度Cs(n−1)を、直近で算出した平滑混雑度Cs(n)に置き換えて記憶する(ステップS8)。
The congestion
Cs(n)=a2・C(n)+(1−a2)・Cs(n−1)−1…(2) Cs (n) = a 2 · C (n) + (1−a 2 ) · Cs (n−1) −1 (2)
但し、この式の算出結果が0未満となる場合には、Cs(n)の値は0として設定される。また、a2は、0<a2≦1で表される係数であり、例えば一例として0.9である。 However, when the calculation result of this expression is less than 0, the value of Cs (n) is set to 0. Further, a 2 is a coefficient represented by 0 <a 2 ≦ 1, and is, for example, 0.9.
混雑度算出部22で算出された平滑混雑度Cs(n)は、追尾処理部11に通知される(ステップS9)。追尾処理部11は、詳しくは後述するが、混雑度算出部22から通知された平滑混雑度Cs(n)の値に応じた追尾処理を行う。
The smoothness degree Cs (n) calculated by the congestion
追尾処理部11は、複数の物標130のなかから、追尾物標140を特定し、且つ、当該追尾物標140の追尾処理を行うように構成されている。追尾物標140は、例えば、表示器4に表示された、複数の物標130を示すシンボル等を基に、オペレータによって選択される。オペレータによる、追尾物標140の選択指令は、例えば、オペレータが操作装置(図示せず)を操作することにより、発せられる。本実施形態では、追尾処理部11は、特に説明無き限り、X−Y座標系を基準に、処理を行うように構成されている。
The
追尾処理部11は、nスキャン時点(最新のスキャン時点)における、追尾物標140の平滑速度Vs(n)を算出するように構成されている。また、追尾処理部11は、平滑速度Vs(n)を表示器4に表示させるように、構成されている。
The
追尾処理部11は、特徴情報メモリ12と、選別部13と、関連付け部14と、運動推定部15と、を有している。
The
特徴情報メモリ12は、信号処理部9及びエコー検出部10から出力されたデータを蓄積するように構成されている。特徴情報メモリ12は、(n−T)スキャン時点から、nスキャン時点までの各時点における、全ての物標エコー像120について、特徴情報データ201〜212を蓄積している。尚、定数Tは、予め設定されている値であり、例えば、数十程度である。
The
選別部13は、エコー選別処理を行うように構成されている。具体的には、選別部13は、nスキャン時点(最新のスキャン時点)において、観測位置Xo(n)(追尾物標140の追尾代表点)が存在すると推定される領域を特定する。
The sorting
具体的には、選別部13は、追尾物標140について、運動推定部15で算出された予測位置Xp(n)の座標データを参照する。即ち、選別部13は、nスキャン時点において、追尾代表点Pが存在していると推定される予測位置Xp(n)の座標データを取得する。選別部13は、当該予測位置Xp(n)を中心とする選別領域S(n)を設定する。この選別領域S(n)は、例えば、上記予測位置Xp(n)を中心とする円状の領域である。関連付け部14は、この選別領域S(n)内を検索する。
Specifically, the sorting
ここで、選別部13は、混雑度算出部22によって算出された平滑混雑度Cs(n)の値に応じて、選別領域S(n)の範囲を決定する。具体的には、選別部13は、例えば一例として、平滑混雑度Cs(n)が所定の閾値N4以下の場合、選別領域S(n)の半径
をr1として設定する。一方、選別部13は、平滑混雑度Cs(n)が所定の閾値N4を超える場合、選別領域S(n)の半径を、r1よりも小さいr2として設定する。
Here, the
関連付け部14は、上記選別領域S(n)内に存在する1つ又は複数の物標エコー像の追尾代表点Pのなかから、追尾物標140の観測位置Xo(n)を、過去の追尾ターゲット情報(エコー面積、位置誤差等)から算出される各追尾代表点Pの尤度に基づき特定する。
The associating
具体的には、例えば、関連付け部14は、各追尾代表点Pについて複数の尤度Lhn(n=1,2,…)を算出し、これらの尤度Lhnを合成した結合尤度が最も高い追尾代表点Pを、追尾物標140の観測位置Xo(n)として運動推定部15へ出力している。なお、上述した尤度Lhnとは、各追尾代表点Pが、追尾物標の追尾代表点Pである尤もらしさを示す度合いであって、追尾物標の特徴量と、選別領域S(n)内に存在する各物標の特徴量とに基づいて算出される。ここでの特徴量としては、例えば一例として、物標のエコー面積が挙げられる。
Specifically, for example, the associating
関連付け部14では、各尤度Lhnに対応して閾値Nan(n=1,2,…)が設定されている。そして、関連付け部14は、追尾代表点Pの各尤度Lhnを対応する閾値Nanと比較し、その比較結果に応じて、その追尾代表点Pが、追尾物標の追尾代表点Pの候補となり得るか否かを判定する。具体的には、関連付け部14は、選別領域S(n)内に存在するある追尾代表点Pにおける尤度Lhnが閾値Nan以上の場合、当該追尾代表点Pを、追尾物標の追尾代表点Pの候補になり得る追尾代表点であると判定する。一方、前記のある追尾代表点Pにおける尤度Lhnが閾値Nan未満の場合、当該追尾代表点Pを、追尾物標の追尾代表点Pの候補から除外する。
In the associating
ここで、関連付け部14は、平滑混雑度Cs(n)の値に応じて上記閾値Nanの値を決定する。具体的には、関連付け部14は、例えば、平滑混雑度Cs(n)が所定の閾値N5以上の場合、閾値Nanの値を高くする。一方、平滑混雑度Cs(n)が所定の閾値N5未満の場合、閾値Nanの値を低くする。
Here, the
運動推定部15は、X−Y座標系において、追尾物標140の追尾処理を行うように構成されている。本実施形態では、運動推定部15は、地球の表面に対する向きが一定である座標系を用いて、追尾物標140の追尾処理を行う。
The
運動推定部15は、観測誤差(アンテナ5が回転する際に発生するぶれ等に起因する誤差)の影響を平滑化するために設けられている。本実施形態では、運動推定部15は、追尾フィルタ処理として、α−βフィルタ処理を行う。運動推定部15は、追尾物標140の予測位置Xp(n)と、平滑位置Xs(n)と、平滑速度Vs(n)と、を算出するように構成されている。
The
具体的には、運動推定部15は、下記式(3)、(4)、(5)を演算する。
予測位置Xp(n)=Xs(n−1)+T×Vs(n−1)…(3)
平滑位置Xs(n)=Xp(n)+α{Xo(n)−Xp(n)}…(4)
平滑速度Vs(n)=Vs(n−1)+(β/T){Xo(n)−Xp(n)}…(5)
Specifically, the
Predicted position Xp (n) = Xs (n−1) + T × Vs (n−1) (3)
Smooth position Xs (n) = Xp (n) + α {Xo (n) −Xp (n)} (4)
Smoothing speed Vs (n) = Vs (n−1) + (β / T) {Xo (n) −Xp (n)} (5)
尚、平滑位置Xs(n)は、nスキャン時点における、予測位置Xp(n)と、観測位置Xo(n)とに関連付けられた位置を示す。平滑位置Xs(n)は、nスキャン時点において、追尾物標140の追尾代表点P1が到達していると推定される位置を示す。また、平滑速度Vs(n)は、nスキャン時点における、追尾代表点P1の推定速度を示して
いる。
The smooth position Xs (n) indicates a position associated with the predicted position Xp (n) and the observation position Xo (n) at the n-scan time point. The smooth position Xs (n) indicates a position where the tracking representative point P1 of the tracking target 140 is estimated to have reached at the n scan time point. The smoothing speed Vs (n) indicates the estimated speed of the tracking representative point P1 at the n scan time point.
また、Tは、運動推定部15が前回の平滑処理を行ってから上記の平滑処理を行うまでに経過した時間であり、1スキャンに要する時間に相当する。また、αは、平滑位置Xs(n)を算出するために用いられるゲインである。βは、平滑速度Vs(n)を算出するために用いられるゲインである。
T is the time that elapses from when the
ゲインαは、X軸方向の成分αx、及びY軸方向の成分αyを有している。ゲインαは、ゲインα(αx,αy)としても表現できる。成分αxは、上記式(4)において、X軸方向の成分を算出する際に用いられる。成分αyは、上記式(4)において、Y軸方向の成分を算出する際に用いられる。 The gain α has a component αx in the X-axis direction and a component αy in the Y-axis direction. The gain α can also be expressed as a gain α (αx, αy). The component αx is used when calculating the component in the X-axis direction in the above equation (4). The component αy is used when calculating the component in the Y-axis direction in the above equation (4).
また、ゲインβは、X軸方向の成分βx、及びY軸方向の成分βyを有している。ゲインβは、ゲインβ(βx,βy)としても表現できる。成分βxは、上記式(5)において、X軸方向の成分を算出する際に用いられる。また、成分βyは、上記式(5)において、Y軸方向の成分を算出する際に用いられる。 The gain β has a component βx in the X-axis direction and a component βy in the Y-axis direction. The gain β can also be expressed as a gain β (βx, βy). The component βx is used when calculating the component in the X-axis direction in the above equation (5). The component βy is used when calculating the component in the Y-axis direction in the above equation (5).
平滑位置Xs(n)は、予測位置Xp(n)と、観測位置Xo(n)と、を結ぶ線分LS1上に位置する。 The smooth position Xs (n) is located on a line segment LS1 connecting the predicted position Xp (n) and the observation position Xo (n).
上記の構成によると、ゲインαが小さいほど、平滑位置Xs(n)は、予測位置Xp(n)に近くなる。また、ゲインβが小さいほど、平滑速度Vs(n)の変化量は、小さくなる。このため、運動推定部15での算出結果は、平滑化された度合いが大きくなり、追尾物標140の観測誤差に起因するばらつき量が小さくなる。但し、追尾物標140の追尾代表点P1の変針運動を、運動推定部15の演算結果に反映させる際の応答は、ゲインα,βが小さいほど、遅くなる。
According to the above configuration, the smoothed position Xs (n) is closer to the predicted position Xp (n) as the gain α is smaller. Further, the smaller the gain β, the smaller the amount of change in the smoothing speed Vs (n). For this reason, the degree of smoothing of the calculation result in the
一方、ゲインαが大きいほど、平滑位置Xs(n)は、観測位置Xo(n)に近くなる。また、ゲインβが大きいほど、平滑速度Vs(n)を平滑化する度合いは、小さくなる。このため、運動推定部15は、追尾物標140の追尾代表点P1の変針運動を、追尾フィルタ処理において、応答性よく反映できる。このため、運動推定部15は、追尾物標140の変針運動に対する追従性を高めることができる。尚、変針運動とは、追尾物標140の向きを変える運動をいう。但し、ゲインα,βが大きいほど、平滑位置Xs(n)及び平滑速度Vs(n)について、スキャン時点毎の変動量が大きくなってしまう。
On the other hand, as the gain α is larger, the smooth position Xs (n) is closer to the observation position Xo (n). Further, the degree of smoothing the smoothing speed Vs (n) decreases as the gain β increases. For this reason, the
そして、運動推定部15は、混雑度算出部22から通知された平滑混雑度Cs(n)の値に応じて、ゲインα,βの値を設定する。ゲインα,βは、閾値N1〜N3(0<N3<N2<N1)と平滑混雑度Cs(n)との関係性に応じて、例えば一例として、以下のように表すことができる。
Then, the
α=β=0(N1<Cs(n)の場合)
α=0.5α1、β=0.5β1(N2<Cs(n)≦N1の場合)
α=α1、β=β1(N3<Cs(n)≦N2の場合)
α=1.5α1、β=1.5β1(Cs(n)≦N3の場合)
但し、α1及びβ1は、所定の定数である。
α = β = 0 (when N 1 <Cs (n))
α = 0.5α 1 , β = 0.5β 1 (when N 2 <Cs (n) ≦ N 1 )
α = α 1 , β = β 1 (when N 3 <Cs (n) ≦ N 2 )
α = 1.5α 1 , β = 1.5β 1 (when Cs (n) ≦ N 3 )
However, α 1 and β 1 are predetermined constants.
上述のように、平滑混雑度Cs(n)に応じてゲインα,βの値を設定することにより、以下のような効果を得ることができる。 As described above, the following effects can be obtained by setting the values of the gains α and β in accordance with the smoothness congestion degree Cs (n).
平滑混雑度が非常に高い場合(本実施形態の場合、N1<Cs(n))、追尾物標の予
測位置付近に多数の物標が存在すると推測される。この場合、関連付け部14による関連付け処理を行う際に、追尾対象が、本来追尾したい物標から他の物標に乗り移ってしまう可能性(いわゆる乗り移りが生じる可能性)が高くなる。よって、このように乗り移りの可能性が高い状況下においては、ゲインα及びβを0に設定していわゆる予測追尾を行うことにより、追尾対象の乗り移りを防止できる。なお、予測追尾とは、平滑位置Xs(n)を算出する際に観測位置Xo(n)を考慮せず、予測位置Xp(n)を平滑位置として算出する追尾のことである。
When the smoothness congestion level is very high (in this embodiment, N 1 <Cs (n)), it is estimated that a large number of targets exist near the predicted position of the tracking target. In this case, when the association process by the
平滑混雑度がやや高い場合(本実施形態の場合、N2<Cs(n)≦N1)、追尾物標の予測位置付近には比較的多数の物標が存在すると推測される。この場合、関連付け部14による関連付け処理を行う際に、乗り移りが生じる可能性がやや高くなる。よって、このように乗り移りの可能性がやや高い状況下においては、ゲインα及びβをやや低い値(本実施形態の場合、α=0.5α1、β=0.5β1)に設定して、いわゆる予測追尾に近い追尾処理を行うことにより、追尾対象の乗り移りの可能性を抑えつつ、追尾物標の変針運動にもある程度対応可能な追尾処理を行うことができる。
When the smoothness congestion level is slightly high (in this embodiment, N 2 <Cs (n) ≦ N 1 ), it is estimated that a relatively large number of targets exist near the predicted position of the tracking target. In this case, when the association process by the
一方、平滑混雑度が0に近い値の場合(本実施形態の場合、Cs(n)≦N3)、追尾物標の予測位置付近にはほとんど物標が存在していない、又は追尾物標の予測位置付近には物標が存在していないと推測される。この場合には、上述のような乗り移りの危険性が非常に低いため、ゲインα及びβの値を通常時よりも大きな値(α=1.5α1、β=1.5β1)に設定することにより、追尾物標の変針運動に対する追従性を高める追尾処理を行うことができる。 On the other hand, when the smoothness congestion is a value close to 0 (in this embodiment, Cs (n) ≦ N 3 ), there is almost no target near the predicted position of the tracking target, or the tracking target It is estimated that there is no target near the predicted position. In this case, since the risk of transfer as described above is very low, the values of the gains α and β are set to values larger than normal (α = 1.5α 1 , β = 1.5β 1 ). Thus, it is possible to perform a tracking process that improves the tracking performance of the tracking target with respect to the changing needle movement.
なお、平滑混雑度がN3<Cs(n)≦N2の場合、ゲインα及びβの値は、追尾対象の乗り移りのリスク低減、及び追尾物標に対する追従性の確保、の双方のバランスのとれた値(α=α1、β=β1)に設定される。 Note that when the smoothness congestion level is N 3 <Cs (n) ≦ N 2 , the values of the gains α and β are a balance between the reduction of the risk of transfer of the tracking target and the ensuring of the tracking performance of the tracking target. The obtained value (α = α 1 , β = β 1 ) is set.
運動推定部15は、上述のように平滑混雑度に応じて適切に設定されたゲインα,βを用いて算出した各データを、選別部13、混雑度算出処理部20、及び表示器4へ出力する。具体的には、運動推定部15は、予測位置Xp(n)の座標データを、選別部13及び混雑度算出処理部20へ出力する。選別部13では、当該データは、(n+1)スキャン時点における選別処理に用いられる。一方、混雑度算出処理部20では、当該データは、混雑度の算出対象となるセルBを特定するために用いられる。また、運動推定部15は、平滑位置Xs(n)、及び平滑速度Vs(n)を特定するデータを、表示器4へ出力する。
The
表示器4は、例えばカラー表示可能な液晶ディスプレイである。表示器4は、各物標エコー像120の画像データを用いて、表示画面に、各物標エコー像120を表示する。また、表示器4は、平滑速度Vs(n)を、画像として表示する。これにより、表示器4の表示画面には、追尾物標140(物標エコー像121)について、平滑速度Vs(n)を示す画像が表示される。レーダ装置1のオペレータは、表示器4に表示されたレーダ映像を確認することにより、追尾物標140の運動状態を確認することができる。
The display 4 is a liquid crystal display capable of color display, for example. The display 4 displays each
[効果]
以上のように、本実施形態に係る追尾処理装置3では、追尾処理部が、追尾物標140の予測位置Xp(n)が含まれるエリア(本実施形態の場合、セルB)内に存在する物標の混雑度(本実施形態の場合、平滑混雑度Cs(n))に応じて、追尾処理を行っている。こうすると、追尾対象が本来追尾したい物標から他の物標に乗り移ってしまう可能性(乗り移りの可能性)に応じて、追尾処理を適切に制御することができる。具体的には、本実施形態の場合、乗り移りの可能性が高い状況下、すなわち平滑混雑度Cs(n)が高い
状況下においては予測追尾を行っているため、乗り移りを防止することができる。一方、乗り移りの可能性が低い状況下、すなわち平滑混雑度Cs(n)が低い状況下においては、追尾物標の変針運動に対する追従性が高い追尾処理を行っている。
[effect]
As described above, in the
従って、追尾処理装置3によれば、周囲の環境によらず、追尾対象を正確に追尾できる。
Therefore, according to the
また、追尾処理装置3では、追尾物標140の予測位置Xp(n)が含まれるエリア内に存在する物標の物標数に基づいて混雑度(平滑混雑度Cs(n))を算出しているため、混雑度を適切に算出することができる。
Further, the
また、追尾処理装置3では、所定のエリアAが分割された複数のセルBのうち、追尾対象の予測位置Xp(n)が含まれるセルBの混雑度(平滑混雑度Cs(n))が算出されている。これにより、追尾対象が含まれる比較的限られた領域を対象として算出された混雑度に基づいて追尾対象の追尾が行われるため、より正確に追尾対象を追尾できる。
Further, in the
また、追尾処理装置3によれば、複数のタイミングでの物標数に基づいて混雑度を算出することができる。こうすると、例えば一例として、そのエコー強度が物標として算出される閾値前後となるような物標、の物標数が時間的に平均化されるため、混雑度がタイミング毎に大きく変動してしまうことを抑制できる。
Further, according to the
また、追尾処理装置3によれば、複数のタイミングでの混雑度に基づいて平滑混雑度を算出することができる。こうすると、上述の場合と同様、そのエコー強度が物標として算出される閾値前後となるような物標、の物標数に基づいて算出される混雑度が時間的に平均化されるため、混雑度がタイミング毎に大きく変動してしまうことを抑制できる。
Further, according to the
また、追尾処理装置3では、追尾フィルタ処理を施す際に用いられるゲインα及びβが、平滑混雑度Cs(n)の値に応じて設定されている。具体的には、乗り移りが発生する可能性が高い場合(平滑混雑度Cs(n)が高い場合)には、ゲインα及びβが低い値に設定される。これにより、いわゆる予測追尾が行われるため、追尾対象が本来追尾したい物標から他の物標に乗り移ってしまう可能性を低減できる。一方、乗り移りが発生する可能性が低い場合(平滑混雑度Cs(n)が低い場合)には、ゲインα及びβが高い値に設定される。これにより、追尾物標の変針運動に対する追従性が高い追尾処理を行うことができる。
In the
また、追尾処理装置3では、選別領域S(n)の面積が、平滑混雑度Cs(n)の値に応じて設定される。具体的には、乗り移りが発生する可能性が高い場合(平滑混雑度Cs(n)が高い場合)には、選別領域S(n)の面積が小さく設定される。これにより、追尾対象となり得る物標の候補を絞り込むことができるため、乗り移りの危険性を低減できる。一方、乗り移りが発生する可能性が低い場合(平滑混雑度Cs(n)が低い場合)には、選別領域S(n)の面積が大きく設定される。これにより、大きく変針した追尾物標についても追尾対象となり得る物標の候補とすることができるため、追尾物標の変針運動に対する追従性が高い追尾処理を行うことができる。
Further, in the
また、追尾処理装置3では、関連付け部14で用いられる閾値Nanの値が、平滑混雑度Cs(n)の値に応じて設定される。乗り移りが発生する可能性が高い場合(平滑混雑度Cs(n)が高い場合)には、閾値Nanの値が大きく設定される。これにより、追尾対象となり得る物標の候補を絞り込むことができるため、乗り移りの危険性を低減できる。一方、乗り移りが発生する可能性が低い場合(平滑混雑度Cs(n)が低い場合)には、閾値Nanの値が小さく設定される。これにより、大きく変針した追尾物標についても
追尾対象となり得る物標の候補とすることができるため、追尾物標の変針運動に対する追従性が高い追尾処理を行うことができる。
Further, the
[変形例]
以上、本発明の実施形態について説明したが、本発明はこれらに限定されるものではなく、本発明の趣旨を逸脱しない限りにおいて種々の変更が可能である。
[Modification]
As mentioned above, although embodiment of this invention was described, this invention is not limited to these, A various change is possible unless it deviates from the meaning of this invention.
(1)上記実施形態では、あるスキャン時にカウントされた物標数を、各スキャンに対応する物標数としているが、これに限らない。具体的には、各スキャン時にカウントされた物標数と、当該スキャン前後の他のスキャン時にカウントされた物標数とを平均化した値を、前記各スキャンに対応する物標数としてもよい。 (1) In the above embodiment, the target number counted during a certain scan is the target number corresponding to each scan. However, the present invention is not limited to this. Specifically, a target value corresponding to each scan may be a value obtained by averaging the target count counted at each scan and the target count counted at other scans before and after the scan. .
(2)上記実施形態では、ゲインα,βの値を、平滑混雑度Cs(n)の値と閾値N1〜N3の値との関係性に応じて段階的に設定したが、この限りでない。具体的には、ゲインα,βの値が平滑混雑度Cs(n)の値に応じて線形的に変化するように、ゲインα,βの値を設定してもよい。ゲインα,βの値は、例えば一例として、以下のように表すことができる。 (2) In the above embodiment, the values of the gains α and β are set in stages according to the relationship between the value of the smoothness congestion level Cs (n) and the values of the thresholds N 1 to N 3. Not. Specifically, the values of the gains α and β may be set so that the values of the gains α and β change linearly according to the value of the smoothness congestion degree Cs (n). The values of the gains α and β can be expressed as follows, for example.
α=ν・α2
β=ν・β2
但し、N1<Cs(n)の場合、ν=0であり、N1≧Cs(n)の場合ν=1−(Cs(n)−N2)/(N1−N2)である。また、α2及びβ2は、所定の定数である。
α = ν · α 2
β = ν · β 2
However, when N 1 <Cs (n), ν = 0, and when N 1 ≧ Cs (n), ν = 1− (Cs (n) −N 2 ) / (N 1 −N 2 ). . Α 2 and β 2 are predetermined constants.
ゲインα,βを上述した式を用いて設定することにより、平滑混雑度Cs(n)の値がN1以上となる場合におけるゲインα,β値を、平滑混雑度Cs(n)の値に応じて線形的に変化させることができる。これにより、ゲインα,β値を、平滑混雑度Cs(n)の値に応じてより適切に設定することができる。 Gain alpha, by setting using the formula described above the beta, the gain alpha, beta value when the value of the smoothness congestion Cs (n) is the N 1 or more, the value of the smoothness congestion Cs (n) It can be changed linearly accordingly. Thereby, the gains α and β values can be set more appropriately according to the value of the smoothness congestion level Cs (n).
また、上記実施形態では、選別領域S(n)の範囲を、平滑混雑度Cs(n)の値と閾値N4との関係性に応じて段階的に設定したが、この限りでない。具体的には、上述した場合と同じように、選別領域S(n)の範囲が平滑混雑度Cs(n)の値に応じて線形的に変化するように、選別領域S(n)の範囲を設定してもよい。 In the above embodiment, the range of the selection area S (n), was set stepwise in accordance with the relationship between the value and the threshold value N 4 of smoothness congestion Cs (n), this shall not apply. Specifically, as in the case described above, the range of the selection region S (n) so that the range of the selection region S (n) changes linearly according to the value of the smoothness congestion degree Cs (n). May be set.
(3)図7は、変形例に係る追尾処理装置の混雑度算出部22によって算出される混雑度の算出方法について説明するための模式図である。図7では、複数のセルB(図6参照)のうち、追尾物標の予測位置Xp(n)付近のセルBを示している。上記実施形態では、追尾物標の予測位置Xp(n)が含まれるセルBにおける平滑物標数Ns(n)を、混雑度C(n)として算出したが、これに限らない。本変形例では、混雑度C(n)を、線形補完を用いて算出している。
(3) FIG. 7 is a schematic diagram for explaining a calculation method of the congestion degree calculated by the congestion
本変形例では、混雑度算出部22が混雑度C(n)を算出する際、複数のセルBに対応して記憶されている平滑物標数Ns(n)のうち、以下に示す条件を満たすセルBの平滑物標数Ns(n)が用いられる。具体的には、混雑度C(n)は、図7を参照して、各中心点CP1〜CP4で囲まれる正方形状の領域内に予測位置Xp(n)が含まれる4つのセルB(図7では、B1〜B4)のそれぞれに対応して記憶される平滑物標数Ns(n)(図7では、Ns1(n)〜Ns4(n))を用いて、以下の式(6)によって算出される。そして、本変形例の混雑度算出部22は、上記実施形態の場合と同様、それまで記憶されていた平滑混雑度C(n−1)を、直近で算出した平滑混雑度C(n)に置き換えて記憶する。
In the present modification, when the congestion
C(n)=Ns’(n)・ωy+Ns”(n)・(1−ωy)…(6) C (n) = Ns ′ (n) · ωy + Ns ″ (n) · (1−ωy) (6)
但し、Ns’(n)=Ns1(n)・(1−ωx)+Ns2(n)・ωx、Ns”(n)=Ns3(n)・(1−ωx)+Ns4(n)・ωx、である。なお、ωxは、CP1からCP2までの距離(すなわち、CP3からCP4までの距離)に対する、CP1からXp(n)までの距離のx軸方向成分(すなわち、CP3からXp(n)までの距離のx軸方向成分)の比、である。また、ωyは、CP1からCP3までの距離(すなわち、CP2からCP4までの距離)に対する、CP3からXp(n)までの距離のy軸方向成分(すなわち、CP4からXp(n)までの距離のy軸方向成分)の比、である。 However, Ns ′ (n) = Ns 1 (n) · (1−ωx) + Ns 2 (n) · ωx, Ns ″ (n) = Ns 3 (n) · (1−ωx) + Ns 4 (n) · Note that ωx is the x-axis direction component of the distance from CP 1 to Xp (n) with respect to the distance from CP 1 to CP 2 (ie, the distance from CP 3 to CP 4 ) (ie, (X-axis direction component of the distance from CP 3 to Xp (n)), and ωy is the CP to the distance from CP 1 to CP 3 (ie, the distance from CP 2 to CP 4 ). 3 is the ratio of the y-axis direction component of the distance from 3 to Xp (n) (that is, the y-axis direction component of the distance from CP 4 to Xp (n)).
上述のように算出された平滑混雑度C(n)は、上記実施形態の場合と同様、追尾処理部11に通知される。そして、追尾処理部11は、上記実施形態の場合と同様、平滑混雑度C(n)の値に応じた追尾処理を行う。
The smoothness degree C (n) calculated as described above is notified to the
以上のように、本変形例に係る追尾処理装置によれば、上述した実施形態に係る追尾処理装置の場合と同様、周囲の環境によらず、追尾対象を正確に追尾できる。 As described above, according to the tracking processing device according to the present modification, the tracking target can be accurately tracked regardless of the surrounding environment, as in the case of the tracking processing device according to the above-described embodiment.
そして、本変形例に係る追尾処理装置によれば、線形補完を用いて混雑度を算出しているため、物標の混雑度をより正確に算出することができる。 And according to the tracking processing apparatus concerning this modification, since the congestion degree is calculated using linear interpolation, the congestion degree of the target can be calculated more accurately.
(4)上記実施形態では、追尾物標の予測位置Xp(n)が含まれるセルBにおける平滑物標数Ns(x,y,n)を、混雑度C(n)として算出したが、これに限らない。具体的には、平滑物標数を算出せずに、nスキャン時における物標数N(x,y,n)を、混雑度C(n)として算出してもよい。 (4) In the above embodiment, the smooth target characteristic Ns (x, y, n) in the cell B including the predicted position Xp (n) of the tracking target is calculated as the congestion degree C (n). Not limited to. Specifically, the target number N (x, y, n) at the time of n scans may be calculated as the congestion degree C (n) without calculating the smooth target number.
(5)上記実施形態及び変形例では、平滑混雑度(又は混雑度)に応じて、選別領域S(n)の範囲、ゲインα,βの値、及び閾値Nanの値を制御していたが、これに限らず、これらのうちの少なくとも1つが、混雑度に応じた変数として制御されればよい。 (5) In the embodiment and the variants described above, in accordance with the smoothness congestion degree (or congestion), the range of selection area S (n), the gain alpha, has been controlled value of beta, and the threshold value Na n However, the present invention is not limited to this, and at least one of them may be controlled as a variable according to the degree of congestion.
1 レーダ装置
3 追尾処理装置
11 追尾処理部
22 混雑度算出部
DESCRIPTION OF
Claims (7)
前記追尾対象の予測位置が含まれるエリア内に存在する物標の混雑の度合いである混雑度を算出する混雑度算出部と、
を備え、
前記追尾処理部は、前記混雑度算出部によって算出された前記混雑度の値が低いほど前記追尾対象への追従性が高くなるようにゲインを設定することを特徴とする、追尾処理装置。 A tracking processing unit that performs processing for tracking a tracking target;
A congestion degree calculating unit for calculating a degree of congestion that is a degree of congestion of a target existing in an area including the predicted position of the tracking target;
With
The tracking processing device is characterized in that the tracking processing unit sets a gain such that the lower the value of the congestion degree calculated by the congestion degree calculating unit, the higher the followability to the tracking target.
前記追尾処理部は、前記混雑度算出部によって算出された前記混雑度の値が低いほど、段階的に前記追尾対象への追従性が高くなるようにゲインを設定することを特徴とする、追尾処理装置。 The tracking processing device according to claim 1,
The tracking processing unit sets the gain so that the tracking performance to the tracking target increases step by step as the congestion level calculated by the congestion level calculation unit is lower. Processing equipment.
前記混雑度算出部は、前記エリア内に存在する前記物標の数としての物標数に基づいて前記混雑度を算出することを特徴とする、追尾処理装置。 In the tracking processing device according to claim 1 or 2,
The tracking processing device, wherein the congestion degree calculating unit calculates the congestion degree based on a target number as the number of the target existing in the area.
前記エリアが分割された複数のセルのそれぞれに存在する前記物標数をカウントし、前記複数のセル毎にカウントされた前記物標数を、各前記セルに対応させて記憶するエコー分布生成部を更に備え、
前記混雑度算出部は、前記追尾対象の前記予測位置が含まれる前記セル内に存在する前記物標の前記混雑度を算出することを特徴とする、追尾処理装置。 In the tracking processing device according to claim 3,
An echo distribution generation unit that counts the target number present in each of a plurality of cells into which the area is divided and stores the target number counted for each of the plurality of cells in association with each cell. Further comprising
The tracking processing apparatus, wherein the congestion degree calculation unit calculates the congestion degree of the target existing in the cell including the predicted position of the tracking target.
前記混雑度算出部は、複数のタイミングでの前記物標数に基づいて前記混雑度を算出することを特徴とする、追尾処理装置。 In the tracking processing device according to claim 3 or 4,
The tracking processing device, wherein the congestion degree calculation unit calculates the congestion degree based on the target numbers at a plurality of timings.
前記混雑度算出部は、複数のタイミングでの前記混雑度に基づいて平滑混雑度を算出し、
前記追尾処理部は、前記混雑度算出部によって算出された前記平滑混雑度の値に応じて、前記追尾対象を追尾する処理を行うことを特徴とする、追尾処理装置。 In the tracking processing apparatus of any one of Claims 1-5,
The congestion degree calculation unit calculates a smooth congestion degree based on the congestion degree at a plurality of timings,
The tracking processing device, wherein the tracking processing unit performs processing for tracking the tracking target according to the value of the smoothness congestion level calculated by the congestion level calculation unit.
前記追尾対象の予測位置が含まれるエリア内に存在する物標の混雑の度合いである混雑度を算出するステップと、を含み、
前記追尾する処理を行うステップは、前記混雑度を算出するステップによって算出された前記混雑度の値が低いほど前記追尾対象への追従性が高くなるようにゲインを設定することを特徴とする、追尾処理方法。 Performing a process of tracking the tracking target;
Calculating a degree of congestion which is a degree of congestion of a target existing in an area including the predicted position of the tracking target,
The step of performing the tracking process is characterized in that a gain is set so that the followability to the tracking target becomes higher as the value of the congestion degree calculated by the step of calculating the congestion degree is lower. Tracking processing method.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014241908 | 2014-11-28 | ||
JP2014241908 | 2014-11-28 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016561449A Division JP6334730B2 (en) | 2014-11-28 | 2015-10-08 | Tracking processing apparatus and tracking processing method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018136337A true JP2018136337A (en) | 2018-08-30 |
JP6554205B2 JP6554205B2 (en) | 2019-07-31 |
Family
ID=56074076
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016561449A Active JP6334730B2 (en) | 2014-11-28 | 2015-10-08 | Tracking processing apparatus and tracking processing method |
JP2018084387A Active JP6554205B2 (en) | 2014-11-28 | 2018-04-25 | Tracking processing apparatus and tracking processing method |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016561449A Active JP6334730B2 (en) | 2014-11-28 | 2015-10-08 | Tracking processing apparatus and tracking processing method |
Country Status (5)
Country | Link |
---|---|
US (1) | US10585184B2 (en) |
EP (1) | EP3226025A4 (en) |
JP (2) | JP6334730B2 (en) |
CN (1) | CN107110967B (en) |
WO (1) | WO2016084498A1 (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10399650B2 (en) * | 2017-01-17 | 2019-09-03 | Harris Corporation | System for monitoring marine vessels and determining rendezvouses therebetween and related methods |
US10302769B2 (en) | 2017-01-17 | 2019-05-28 | Harris Corporation | System for monitoring marine vessels using fractal processing of aerial imagery and related methods |
US10200113B2 (en) | 2017-01-17 | 2019-02-05 | Harris Corporation | System for monitoring marine vessels providing expected passenger determination features and related methods |
US11500085B2 (en) * | 2017-09-22 | 2022-11-15 | Furuno Electric Co., Ltd. | Radar device and target tracking method |
AU2020205850A1 (en) * | 2019-01-11 | 2021-08-26 | Maerospace Corporation | System and method for tracking vessels |
KR102205478B1 (en) * | 2019-09-30 | 2021-01-20 | 엘아이지넥스원 주식회사 | Apparatus for Scanning Detection Region |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5400264A (en) * | 1993-04-29 | 1995-03-21 | International Business Machines Corporation | Suboptimal joint probabilistic data association |
JP2001242242A (en) * | 2000-02-29 | 2001-09-07 | Hitachi Ltd | Millimeter-wave radar device with function for improving detecting performance |
JP2001296354A (en) * | 2000-04-12 | 2001-10-26 | Mitsubishi Electric Corp | Apparatus for managing resource |
JP2007068008A (en) * | 2005-09-01 | 2007-03-15 | Matsushita Electric Ind Co Ltd | Automatic tracking apparatus |
JP2009052940A (en) * | 2007-08-24 | 2009-03-12 | Toshiba Corp | Moving object detection device and autonomous moving body |
JP2009271086A (en) * | 2009-08-17 | 2009-11-19 | Mitsubishi Electric Corp | Radar device |
JP2013152157A (en) * | 2012-01-25 | 2013-08-08 | Furuno Electric Co Ltd | Signal selection device, signal selection method, and rader system |
JP2014085274A (en) * | 2012-10-25 | 2014-05-12 | Furuno Electric Co Ltd | Gain setting method, gain setting program, gain setting apparatus, tracking processor, tracking processing method, tracking processing program and rader system |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7315274B2 (en) * | 2004-03-15 | 2008-01-01 | Kongsberg Seatex As | Method and system for determining the position of marine vessels and similar objects |
US8665113B2 (en) * | 2005-10-31 | 2014-03-04 | Wavetronix Llc | Detecting roadway targets across beams including filtering computed positions |
JP4802112B2 (en) * | 2007-02-08 | 2011-10-26 | 株式会社東芝 | Tracking method and tracking device |
US8836570B2 (en) * | 2011-10-26 | 2014-09-16 | Raytheon Canada Limited | Systems and methods for extending maritime domain awareness by sharing radar tracks between vessels |
CN102393901B (en) | 2011-11-18 | 2013-10-30 | 东南大学 | Traffic flow information perception method based on hybrid characteristic and system thereof |
JP6070689B2 (en) * | 2012-02-29 | 2017-02-01 | 日本電気株式会社 | Flow line information generation system, flow line information generation method, and flow line information generation program |
JP6084812B2 (en) * | 2012-10-29 | 2017-02-22 | 古野電気株式会社 | Tracking processing apparatus and tracking processing method |
JP6084810B2 (en) | 2012-10-29 | 2017-02-22 | 古野電気株式会社 | Tracking processing apparatus and tracking processing method |
JP5735673B2 (en) * | 2014-02-03 | 2015-06-17 | 東京計器株式会社 | Target tracking device and target tracking method |
-
2015
- 2015-10-08 EP EP15862592.1A patent/EP3226025A4/en active Pending
- 2015-10-08 JP JP2016561449A patent/JP6334730B2/en active Active
- 2015-10-08 WO PCT/JP2015/078625 patent/WO2016084498A1/en active Application Filing
- 2015-10-08 CN CN201580059386.6A patent/CN107110967B/en active Active
- 2015-10-08 US US15/528,700 patent/US10585184B2/en active Active
-
2018
- 2018-04-25 JP JP2018084387A patent/JP6554205B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5400264A (en) * | 1993-04-29 | 1995-03-21 | International Business Machines Corporation | Suboptimal joint probabilistic data association |
JP2001242242A (en) * | 2000-02-29 | 2001-09-07 | Hitachi Ltd | Millimeter-wave radar device with function for improving detecting performance |
JP2001296354A (en) * | 2000-04-12 | 2001-10-26 | Mitsubishi Electric Corp | Apparatus for managing resource |
JP2007068008A (en) * | 2005-09-01 | 2007-03-15 | Matsushita Electric Ind Co Ltd | Automatic tracking apparatus |
JP2009052940A (en) * | 2007-08-24 | 2009-03-12 | Toshiba Corp | Moving object detection device and autonomous moving body |
JP2009271086A (en) * | 2009-08-17 | 2009-11-19 | Mitsubishi Electric Corp | Radar device |
JP2013152157A (en) * | 2012-01-25 | 2013-08-08 | Furuno Electric Co Ltd | Signal selection device, signal selection method, and rader system |
JP2014085274A (en) * | 2012-10-25 | 2014-05-12 | Furuno Electric Co Ltd | Gain setting method, gain setting program, gain setting apparatus, tracking processor, tracking processing method, tracking processing program and rader system |
Also Published As
Publication number | Publication date |
---|---|
US20170261608A1 (en) | 2017-09-14 |
JPWO2016084498A1 (en) | 2017-09-07 |
US10585184B2 (en) | 2020-03-10 |
WO2016084498A1 (en) | 2016-06-02 |
CN107110967A (en) | 2017-08-29 |
JP6554205B2 (en) | 2019-07-31 |
JP6334730B2 (en) | 2018-05-30 |
EP3226025A4 (en) | 2018-07-18 |
EP3226025A1 (en) | 2017-10-04 |
CN107110967B (en) | 2020-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6554205B2 (en) | Tracking processing apparatus and tracking processing method | |
JP5658871B2 (en) | Signal processing apparatus, radar apparatus, signal processing program, and signal processing method | |
JP5697910B2 (en) | Threshold setting method, target detection method, threshold setting device, target detection device, threshold setting program, and target detection program | |
JP6084812B2 (en) | Tracking processing apparatus and tracking processing method | |
JP6084810B2 (en) | Tracking processing apparatus and tracking processing method | |
JP5398195B2 (en) | Radar equipment | |
JP6762726B2 (en) | Radar device and target tracking method | |
JP2008164545A (en) | Moving target detecting device, moving target detection method, and moving target detection program | |
EP3290947B1 (en) | Signal processing device and radar device | |
JP2012103197A (en) | Moving target adaptive type scan correlation method | |
JP5074718B2 (en) | Marine radar | |
JP6466702B2 (en) | Tracking processing device, radar device, and tracking processing method | |
JP5731859B2 (en) | Interference cancellation apparatus, radar apparatus, and interference cancellation method | |
JP2011242253A (en) | Echo signal processing device, radar device, echo signal processing program and echo signal processing method | |
JP6084811B2 (en) | Tracking processing apparatus and tracking processing method | |
JP2011191141A (en) | False image reducing device, radar device, false image reducing method, and false image reducing program | |
JP6180004B2 (en) | Tracking processing apparatus and tracking processing method | |
JP6192915B2 (en) | Gain setting method, gain setting program, and gain setting device | |
JP6235557B2 (en) | Radar device | |
JP2016121970A (en) | Echo signal processing device, radar device, echo signal processing method, and program | |
JP6339074B2 (en) | Sea state detection device, radar device, sea state detection method, and program | |
JP6466768B2 (en) | Radar equipment | |
JP2014092458A (en) | Radar receiver | |
JP6084813B2 (en) | Echo signal processing device, radar device, echo signal processing method, and echo signal processing program | |
JP7023633B2 (en) | Radar received signal processing equipment, programs and methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180517 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20181108 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20181120 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20181227 |
|
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: 20190702 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190705 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6554205 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 |