JP2005156276A - Radar system for target identification - Google Patents

Radar system for target identification Download PDF

Info

Publication number
JP2005156276A
JP2005156276A JP2003393328A JP2003393328A JP2005156276A JP 2005156276 A JP2005156276 A JP 2005156276A JP 2003393328 A JP2003393328 A JP 2003393328A JP 2003393328 A JP2003393328 A JP 2003393328A JP 2005156276 A JP2005156276 A JP 2005156276A
Authority
JP
Japan
Prior art keywords
target
polarization
feature
unit
received
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
JP2003393328A
Other languages
Japanese (ja)
Other versions
JP4096867B2 (en
Inventor
Masato Imamura
誠人 今村
Hiroshi Suwa
啓 諏訪
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 JP2003393328A priority Critical patent/JP4096867B2/en
Publication of JP2005156276A publication Critical patent/JP2005156276A/en
Application granted granted Critical
Publication of JP4096867B2 publication Critical patent/JP4096867B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To improve target identification performance by narrowing down reference data, in a radar system for target identification for identifying observed targets. <P>SOLUTION: The radar system for target identification is provided with receiver 106 for receiving a plurality of polarization signals, through the use of a plurality of aerials having polarization characteristics difference from one another; a transmitter 101 for generating and transmitting wide-band pulses; a target detecting part 108 for detecting the targets, on the basis of a plurality of received polarization signals received by the receiver 106; and a target identifying part 109 for identifying the targets by combining the feature quantities. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

この発明は、レーダの分野における目標のレンジプロフィールと偏波の情報を用いて目標を識別するレーダ装置に関するものである。   The present invention relates to a radar apparatus for identifying a target using information on a target range profile and polarization in the field of radar.

目標識別を行うレーダ装置に関する従来の技術として、高い距離分解能をもつレーダ装置において観測される目標のレンジプロフィールを利用して目標を識別するものがあり、識別の対象とする候補目標のレンジプロフィールを参照用レンジプロフィールとして予め蓄積しておくライブラリと、観測された目標レンジプロフィールとライブラリに蓄積された参照用レンジプロフィールの相関を計算する手段とを備え、観測された目標レンジプロフィールと相関の高い参照用レンジプロフィールがある場合、観測された目標はその参照用レンジプロフィールに対応する目標であると識別することが開示されている(例えば、非特許文献1参照。)。   As a conventional technique related to a radar apparatus for performing target identification, there is a technique for identifying a target using a target range profile observed in a radar apparatus having a high range resolution. A library pre-stored as a reference range profile and means for calculating the correlation between the observed target range profile and the reference range profile stored in the library, and a reference highly correlated with the observed target range profile In the case where there is a range profile, it is disclosed that an observed target is identified as a target corresponding to the reference range profile (see Non-Patent Document 1, for example).

“Correlation Filters for Aircraft Identification From Radar Range Profiles”,IEEE Transactions on Aerospace and Electronic Systems,vol.29,no.3,pp.741-748,July 1993“Correlation Filters for Aircraft Identification From Radar Range Profiles”, IEEE Transactions on Aerospace and Electronic Systems, vol.29, no.3, pp.741-748, July 1993

従来の目標識別用レーダ装置は以上のように構成されているが、レンジプロフィールの類似した目標が複数ある場合や、ミリ波帯などの目標の大きさに対して波長の短い周波数帯を用いた場合、これらの目標を識別できないという問題が発生する。   The conventional target identification radar device is configured as described above. However, when there are multiple targets with similar range profiles, a frequency band with a shorter wavelength than the target size such as a millimeter wave band is used. If so, the problem arises that these goals cannot be identified.

この発明は、上記のような課題を解決するためになされたもので、レンジプロフィールの類似した目標が複数ある場合や、ミリ波帯などの目標の大きさに対して波長の短い周波数帯を用いた場合にも目標識別性能を向上することができる目標識別用レーダ装置を得ることを目的とする。   The present invention has been made to solve the above-described problems. When there are a plurality of targets having similar range profiles, a frequency band having a short wavelength with respect to a target size such as a millimeter wave band is used. An object of the present invention is to obtain a target identification radar apparatus that can improve the target identification performance even in the case of such a situation.

この発明に係る目標識別用レーダ装置は、広帯域パルスを発生して送信する送信機と、相互に偏波特性の異なる複数のアンテナを用いて複数の偏波信号を受信する受信機と、上記受信機によって受信した複数の受信偏波信号から目標を検出する目標検出部と、上記目標検出部からの信号に基づき特徴量を算出する特徴量抽出部と、上記特徴量抽出部によって抽出された特徴量を用いて目標の識別を行う目標識別部とを備えたものである。   A radar device for target identification according to the present invention includes a transmitter that generates and transmits a broadband pulse, a receiver that receives a plurality of polarization signals using a plurality of antennas having different polarization characteristics, and A target detection unit that detects a target from a plurality of received polarization signals received by the receiver, a feature amount extraction unit that calculates a feature amount based on a signal from the target detection unit, and a feature amount extraction unit And a target identifying unit that identifies a target using a feature amount.

この発明によれば、各偏波チャネルのレンジプロフィールの大まかな形状を示す複数種類の指標と、複数偏波チャネルのレンジプロフィール間の関係を示す複数種類の指標を算出し、特徴量を組み合わせて目標識別を行うことにより、目標を正しく識別する確率を向上できるという効果がある。   According to the present invention, a plurality of types of indicators indicating the rough shape of the range profile of each polarization channel and a plurality of types of indicators indicating the relationship between the range profiles of the plurality of polarization channels are calculated, and the feature amounts are combined. By performing target identification, there is an effect that the probability of correctly identifying the target can be improved.

実施の形態1.
図1はこの発明の実施の形態1による目標識別用レーダ装置の構成を示すブロック図である。目標のレンジプロフィールは目標のアスペクト角(目標の中心線とレーダの視線方向のなす角)に大きく依存するため、各候補目標に関して、あらゆるアスペクト角で観測したレンジプロフィールを参照データとして用意し、これらを全て候補目標レンジプロフィールのライブラリに蓄積しておく必要がある。
Embodiment 1 FIG.
FIG. 1 is a block diagram showing a configuration of a target identifying radar apparatus according to Embodiment 1 of the present invention. Since the target range profile greatly depends on the target aspect angle (the angle between the target center line and the radar line-of-sight direction), for each candidate target, the range profile observed at any aspect angle is prepared as reference data. Must be stored in a library of candidate target range profiles.

しかしながらこの方法では、各候補目標に関してアスペクト角を非常に細かい刻みで変化させて参照データを作成して蓄積することが要求されるため、候補目標レンジプロフィールのライブラリが膨大なものとなってしまう問題が発生する。候補目標レンジプロフィールのライブラリを可能な限り小さくするためには、参照データとして蓄積しておくアスペクト角の刻みを大きくすることが望ましいが、今度は、候補目標レンジプロフィールのライブラリが、観測時の目標のアスペクト角と完全に一致したデータを含まない場合に、目標の識別性能が劣化する問題が発生する。そこでこの実施の形態では、特徴量を組み合わせて目標識別を行うことによって問題を解決する。   However, with this method, it is required to create and store reference data by changing the aspect angle in very fine steps for each candidate target, so that the library of candidate target range profiles becomes enormous. Occurs. In order to make the candidate target range profile library as small as possible, it is desirable to increase the increment of the aspect angle that is stored as reference data, but this time the candidate target range profile library is If the data does not include data that completely matches the aspect angle of the target, there arises a problem that the target identification performance deteriorates. Therefore, in this embodiment, the problem is solved by performing target identification by combining feature amounts.

以下図1について説明する。
101はパルス信号を送信する送信機、102は送受信を切り換える送受切換器、103は偏波切換器、104は第2偏波送受信アンテナ105の偏波特性と直交する偏波特性を有する第1偏波送受信アンテナ、105は第1偏波送受信アンテナ104の偏波特性と直交する偏波特性を有する第2偏波送受信アンテナである。
Hereinafter, FIG. 1 will be described.
101 is a transmitter for transmitting a pulse signal, 102 is a transmission / reception switch for switching transmission / reception, 103 is a polarization switch, and 104 is a first polarization characteristic orthogonal to the polarization characteristic of the second polarization transmission / reception antenna 105. A single polarization transmission / reception antenna 105 is a second polarization transmission / reception antenna having a polarization characteristic orthogonal to the polarization characteristic of the first polarization transmission / reception antenna 104.

なお、第1偏波送受信アンテナ104及び第2偏波送受信アンテナ105における偏波特性が直交する組み合わせとして、例えば、垂直偏波と水平偏波の組み合わせや、右旋円偏波と左旋円偏波の組み合わせなどが考えられる。   Note that as the combinations in which the polarization characteristics of the first polarization transmitting / receiving antenna 104 and the second polarization transmitting / receiving antenna 105 are orthogonal, for example, a combination of vertical polarization and horizontal polarization, right-handed circular polarization and left-handed circular polarization, and the like. A combination of waves can be considered.

106は第1偏波送受信アンテナ104及び第2偏波送受信アンテナ105により受信された受信信号に対する位相検波とA/D変換処理を実施して、その偏波信号の振幅と位相を示すディジタルの受信偏波信号(後述する散乱行列レンジプロフィール)を出力する受信機である。   A digital reception unit 106 performs phase detection and A / D conversion processing on the reception signals received by the first polarization transmission / reception antenna 104 and the second polarization transmission / reception antenna 105, and indicates the amplitude and phase of the polarization signal. It is a receiver that outputs a polarization signal (scattering matrix range profile described later).

107は受信機106から出力された散乱行列レンジプロフィールを蓄積する受信偏波信号記憶部、108は受信偏波信号記憶部より送られた散乱行列レンジプロフィールより目標を検出し、検出された目標位置近傍の散乱行列レンジプロフィールの切り出しを行う目標検出部、109は目標検出部108によって切り出された散乱行列レンジプロフィールから複数の特徴量を算出し、ここで得られた特徴量を、事前に蓄積された既知の目標に関する特徴量のライブラリと比較することで目標を識別する目標識別部、110は目標の検出結果と識別結果を表示する表示部である。   Reference numeral 107 denotes a received polarization signal storage unit that accumulates the scattering matrix range profile output from the receiver 106, and 108 denotes a target position detected from the scattering matrix range profile sent from the received polarization signal storage unit. A target detection unit 109 that cuts out a nearby scatter matrix range profile, 109 calculates a plurality of feature amounts from the scatter matrix range profile cut out by the target detection unit 108, and the feature amounts obtained here are accumulated in advance. A target identification unit 110 for identifying a target by comparing it with a library of feature quantities relating to known targets, and a display unit 110 for displaying the detection result and identification result of the target.

なお、109aは目標検出部108より送られた散乱行列レンジプロフィールに対して偏波基底変換処理を行い、複数偏波チャネルのレンジプロフィールを生成する偏波基底変換処理部、109bは109aによって得られた複数偏波チャネルレンジプロフィールより、識別に利用する複数の特徴量を算出する特徴量算出部、109cは偏波基底変換処理部109aと特徴量算出部109bを備える特徴量抽出部、109dは既知の目標に関する特徴量を蓄積している参照用目標特徴量ライブラリ、109eは特徴量抽出部109cで抽出された特徴量と参照用目標特徴量ライブラリ109dに蓄積された特徴量とから、事前に定められた規則に従って目標識別をおこなう判定部である。   109a is a polarization basis conversion processing unit that performs polarization basis conversion processing on the scattering matrix range profile sent from the target detection unit 108 to generate a range profile of a plurality of polarization channels, and 109b is obtained by 109a. 109c is a feature quantity extraction unit that calculates a plurality of feature quantities used for identification from a plurality of polarization channel range profiles, 109c is a feature quantity extraction unit that includes a polarization basis conversion processing unit 109a and a feature quantity calculation unit 109b, and 109d is known 109e is a reference target feature quantity library 109e that stores feature quantities related to the target, and is determined in advance from the feature quantities extracted by the feature quantity extraction unit 109c and the feature quantities accumulated in the reference target feature quantity library 109d. It is a determination unit that performs target identification in accordance with established rules.

図3はレンジプロフィールの概念を示す説明図である。
図3において、波長の短い電磁波のパルスを送信し、目標によって反射された反射波を受信して時間の関数として記録したものがレンジプロフィールである。レンジプロフィールは目標の反射係数の空間分布をレーダの視線方向に投影したものであると見なすことができる。
FIG. 3 is an explanatory diagram showing the concept of a range profile.
In FIG. 3, a range profile is obtained by transmitting a pulse of an electromagnetic wave having a short wavelength, receiving a reflected wave reflected by a target, and recording it as a function of time. The range profile can be regarded as a projection of the spatial distribution of the target reflection coefficient in the direction of the line of sight of the radar.

時刻 t1に送信されたパルスは、時刻 t2 に図3に示すようなパルス幅 w の短冊状の領域(距離分解能セル)に入射し、そこで反射される。また異なる時刻 t3 には、送信パルスは図に示すごとく別のレンジビンに入射し、そこで反射される。したがって、目標からの反射波を時間の関数として受信し、記録すると、各距離分解能セルからの反射波をアンテナに近い方から順に記録することができる。このようにして目標のレンジプロフィールが観測される。 The pulse transmitted at time t 1 is incident on a strip-shaped region (distance resolution cell) having a pulse width w as shown in FIG. 3 at time t 2 and reflected there. At a different time t 3 , the transmission pulse enters another range bin as shown in the figure and is reflected there. Therefore, when the reflected wave from the target is received and recorded as a function of time, the reflected wave from each distance resolution cell can be recorded in order from the side closer to the antenna. In this way, the target range profile is observed.

次に、散乱行列レンジプロフィールについて説明する。
目標への入射波や、目標からの散乱波は、Maxwellの方程式から明らかなように、本来空間内のベクトルとして表現される。特に、目標が自由空間中で十分遠方に存在する場合には、これらの電磁波は平面波とみなすことができることから、電界(磁界)を進行方向に直交した平面上の二次元ベクトルとして取り扱うことができる。平面波におけるこれら電界(磁界)ベクトルの時間変化の様子は、波の偏り、つまり、いわゆる偏波の概念として理解、分類されている。電磁波をベクトルで表現した場合には、観測対象の散乱特性も、レーダ断面積(Radar Cross Section : RCS)のようなスカラー値ではなく、式(1)、式(2)の散乱行列(Scattering Matrix)として表現される。
Next, the scattering matrix range profile will be described.
The incident wave to the target and the scattered wave from the target are originally expressed as a vector in space, as is clear from Maxwell's equation. In particular, when the target is sufficiently far away in free space, these electromagnetic waves can be regarded as plane waves, so that the electric field (magnetic field) can be treated as a two-dimensional vector on a plane orthogonal to the traveling direction. . The state of time variation of these electric field (magnetic field) vectors in a plane wave is understood and classified as a concept of wave bias, that is, so-called polarization. When electromagnetic waves are expressed as vectors, the scattering characteristics of the observation target are not scalar values such as the radar cross section (Radar Cross Section: RCS), but the scattering matrix (Scattering Matrix) of Equations (1) and (2). ).

Figure 2005156276
Figure 2005156276

Figure 2005156276
Figure 2005156276

時刻t、位置ベクトルrにおける周波数f、波数ベクトルkの平面波E (r, t)は次式のように表される。   The plane wave E (r, t) of the time t, the frequency f at the position vector r, and the wave vector k is expressed by the following equation.

Figure 2005156276
Figure 2005156276

ここで、E0は複素の電界ベクトルであり、水平方向(horizontal:h) の電界成分Ehと垂直方向(vertical:v)の電界成分Evを用いて次式のように表すことができる。 Here, E 0 is a complex electric field vector, which can be expressed by the following equation using a horizontal (horizontal: h) electric field component Eh and a vertical (vertical: v) electric field component Ev.

Figure 2005156276
Figure 2005156276

式(4)において、平面波の偏波状態(polarization state)を特徴づけるのは、ベクトル[ 1,ρ]T(上付きのTは転置)である。そこで、[ 1,ρ]Tのユークリッドノルムを1とするベクトルEJを次式で定義する。 In equation (4), the vector [1, ρ] T (superscript T is transposed) characterizes the polarization state of the plane wave. Therefore, a vector E J in which the Euclidean norm of [1, ρ] T is 1 is defined by the following equation.

Figure 2005156276
Figure 2005156276

ここで上付きの*は複素共役を表す。
以下では、式(4)の電界ベクトルの表現形式をJones Vector形式と呼ぶことにする。
目標への入射波の偏波状態、すなわち送信アンテナの偏波状態をJones Vector形式の複素電界ベクトルEtで表現する。この場合の散乱波の偏波状態を表すベクトルEsは次式で与えられる。
Here, the superscript * represents a complex conjugate.
Hereinafter, the expression format of the electric field vector of Equation (4) is referred to as Jones Vector format.
Polarization state of the incident wave to the target, i.e. to represent the polarization state of transmit antennas at the complex electric field vector E t of Jones Vector format. The vector E s representing the polarization state of the scattered wave in this case is given by the following equation.

Figure 2005156276
Figure 2005156276

上式において、Sは観測対象の散乱特性を表す散乱行列(scattering matrix)であり次式で表される。   In the above equation, S is a scattering matrix representing the scattering characteristics of the observation target, and is represented by the following equation.

Figure 2005156276
Figure 2005156276

ここで、
Shhは入射波の偏波がHの時の散乱波のH成分、
Shvは入射波の偏波がHの時の散乱波のV成分、
Svvは入射波の偏波がVの時の散乱波のV成分、
Svhは入射波の偏波がVの時の散乱波のH成分を表す。
here,
Shh is the H component of the scattered wave when the polarization of the incident wave is H,
Shv is the V component of the scattered wave when the polarization of the incident wave is H,
Svv is the V component of the scattered wave when the polarization of the incident wave is V,
Svh represents the H component of the scattered wave when the polarization of the incident wave is V.

この散乱波を、その偏波状態がJones Vector形式の複素電界ベクトルErで与えられる受信アンテナで受信した場合の受信電圧Vsは次式で与えられる。 The received voltage Vs when this scattered wave is received by a receiving antenna whose polarization state is given by a complex electric field vector Er in the Jones Vector format is given by the following equation.

Figure 2005156276
Figure 2005156276

よって、この場合の受信電力Psは次のように表される。   Therefore, the received power Ps in this case is expressed as follows.

Figure 2005156276
Figure 2005156276

ここで、式(8)と式(9)においては、偏波状態に影響されない定数項は省いている。
なお、送信アンテナと受信アンテナがモノスタティックで構成されている場合は次式が成り立つことが知られている。
Here, in Equations (8) and (9), constant terms that are not affected by the polarization state are omitted.
It is known that the following equation holds when the transmission antenna and the reception antenna are monostatic.

Figure 2005156276
Figure 2005156276

さて、観測対象の散乱行列を観測することが出来るならば、観測対象のレンジプロフィールも、散乱行列のレンジプロフィールとして観測できる。   If the scattering matrix of the observation target can be observed, the range profile of the observation target can also be observed as the range profile of the scattering matrix.

図4は散乱行列レンジプロフィールの概念を示す説明図であり、図2はこの実施の形態1による目標識別方法を示すフローチャートである。
すなわち図4に示すように、散乱行列の4つの要素Shh,Svh,Shv,Svvのそれぞれについて、図3で示したものと同様なレンジプロフィールを得ることができる。これら4つのレンジプロフィールをまとめてここでは散乱行列レンジプロフィールと呼ぶこととする。あるいは、距離分解能セル毎の散乱行列が距離方向に並んだものを散乱行列レンジプロフィールと考えても同義である。なお、図4においては図示の都合上、レンジプロフィールは電界強度の形で示されているが、散乱行列レンジプロフィールの各要素は本来複素数である。
FIG. 4 is an explanatory diagram showing the concept of the scattering matrix range profile, and FIG. 2 is a flowchart showing the target identification method according to the first embodiment.
That is, as shown in FIG. 4, a range profile similar to that shown in FIG. 3 can be obtained for each of the four elements Shh, Svh, Shv, Svv of the scattering matrix. These four range profiles are collectively referred to herein as a scattering matrix range profile. Or it is synonymous even if it considers that the scattering matrix for every distance resolution cell arranged in the distance direction as a scattering matrix range profile. In FIG. 4, for convenience of illustration, the range profile is shown in the form of electric field strength, but each element of the scattering matrix range profile is inherently a complex number.

次に、動作について説明する。
送信機101が広帯域パルスとしてパルス信号を生成すると、送受切換器102が当該パルス信号を偏波切換器103に送り、偏波切換器103が、第1偏波送受信アンテナ104を駆動することにより、そのパルス信号を第1偏波送受信アンテナ104から空間に放射させる。なお、以下では第1偏波送受信アンテナ104と第2偏波送受信アンテナ105の偏波特性が、それぞれ水平偏波と垂直偏波である場合について述べる。
Next, the operation will be described.
When the transmitter 101 generates a pulse signal as a broadband pulse, the transmission / reception switch 102 sends the pulse signal to the polarization switch 103, and the polarization switch 103 drives the first polarization transmitting / receiving antenna 104, The pulse signal is radiated from the first polarization transmitting / receiving antenna 104 to the space. Hereinafter, a case will be described in which the polarization characteristics of the first polarization transmitting / receiving antenna 104 and the second polarization transmitting / receiving antenna 105 are horizontal polarization and vertical polarization, respectively.

第1偏波送受信アンテナ104から空間に放射されたパルス信号は観測対象によって散乱される。
偏波切換器103は、第1偏波送受信アンテナ104と第2偏波送受信アンテナ105の双方を駆動することにより、第1偏波送受信アンテナ104と第2偏波送受信アンテナ105が観測対象によって散乱された散乱波をそれぞれ受信すると、第1偏波送受信アンテナ104と第2偏波送受信アンテナ105が受信信号を受信機106に送り、受信機106が第1偏波送受信アンテナ104と第2偏波送受信アンテナ105において受信された受信信号のそれぞれに対する位相検波処理とA/D変換処理を実施して、それぞれの受信信号の振幅と位相を示すディジタル受信信号S11(m),S21(m)を出力する(ステップST1)。
The pulse signal radiated into the space from the first polarization transmitting / receiving antenna 104 is scattered by the observation target.
The polarization switching unit 103 drives both the first polarization transmission / reception antenna 104 and the second polarization transmission / reception antenna 105, so that the first polarization transmission / reception antenna 104 and the second polarization transmission / reception antenna 105 are scattered by the observation target. When the scattered waves are received, the first polarization transmitting / receiving antenna 104 and the second polarization transmitting / receiving antenna 105 send the received signal to the receiver 106, and the receiver 106 receives the first polarization transmitting / receiving antenna 104 and the second polarization. Digital detection signals S 11 (m) and S 21 (m) indicating the amplitude and phase of each received signal by performing phase detection processing and A / D conversion processing on each received signal received by the transmitting / receiving antenna 105. Is output (step ST1).

なお、Sij(m)は、第j偏波送受信アンテナで送信して第i偏波送受信アンテナで受信された受信信号のm番目(m=1,2,・・・,M)のサンプル値である。ここで、Mはサンプル数である。同様に、送信機101で生成した広帯域パルスを、送受切換器102を介して偏波切換器103に送り、これを第2偏波送受信アンテナ105から目標に照射して第1偏波送受信アンテナ105と第2偏波送受信アンテナ106において受信された受信信号同様の処理を繰り返すことにより、受信信号S12(m),S22(m)を得る(ステップST1)。 S ij (m) is the m-th (m = 1, 2,..., M) sample value of the received signal transmitted by the j-th polarization transmission / reception antenna and received by the i-th polarization transmission / reception antenna. It is. Here, M is the number of samples. Similarly, the broadband pulse generated by the transmitter 101 is sent to the polarization switching unit 103 via the transmission / reception switch 102, and the target is irradiated from the second polarization transmitting / receiving antenna 105 to the first polarization transmitting / receiving antenna 105. By repeating the same processing as the received signal received by the second polarization transmitting / receiving antenna 106, the received signals S 12 (m) and S 22 (m) are obtained (step ST1).

図5は第1偏波送受信アンテナと第2偏波送受信アンテナの各時刻における動作モードを示す説明図である。
第1偏波送受信アンテナ104と第2偏波送受信アンテナ105の動作モードを示し、インターバルが、受信信号S11(m),S21(m),S12(m),S21(m)の一組を得るのに要する処理のひとまとめである。
FIG. 5 is an explanatory diagram showing operation modes at each time of the first polarization transmitting / receiving antenna and the second polarization transmitting / receiving antenna.
The operation modes of the first polarization transmission / reception antenna 104 and the second polarization transmission / reception antenna 105 are shown, and the intervals of the received signals S 11 (m), S 21 (m), S 12 (m), S 21 (m) are shown. It is a batch of processing required to obtain a set.

受信偏波信号記憶部107は、得られた受信信号S11(m),S21(m),S12(m),S22(m)を一時蓄積し、蓄積された信号は、観測対象に関する各分解能セルに散乱行列Sm(m=1,2,・・・,M)の値をもつ受信散乱行列レンジプロフィールであり、式(11)のように表される。
なお、以下では、送信アンテナと受信アンテナの偏波状態の組み合わせのことを偏波チャネルと呼ぶことがある。
The received polarization signal storage unit 107 temporarily stores the obtained received signals S 11 (m), S 21 (m), S 12 (m), and S 22 (m). Is a reception scattering matrix range profile having a value of the scattering matrix Sm (m = 1, 2,..., M) in each resolution cell, and is expressed as in Expression (11).
Hereinafter, the combination of the polarization states of the transmission antenna and the reception antenna may be referred to as a polarization channel.

Figure 2005156276
Figure 2005156276

次に、目標検出部108は受信散乱行列レンジプロフィールから、目標を検出する。
ここでは、例えばPWF(Polarimetric Whitening Filter)や、PNF(Polarimetric Notch Filter)等の方式によってクラッタを抑圧することで、目標検出性能の向上を図っても良い。さらに、目標検出部108は検出された目標位置近傍のサンプルSm(m=M0,M0+1,・・・,M0+MT-1)を切り出す(ステップST2)。
Next, the target detection unit 108 detects a target from the received scattering matrix range profile.
Here, for example, the target detection performance may be improved by suppressing clutter by a method such as PWF (Polarimetric Whitening Filter) or PNF (Polarimetric Notch Filter). Further, the target detection unit 108 cuts out the sample Sm (m = M 0 , M 0 +1,..., M 0 + M T −1) near the detected target position (step ST2).

ここで切り出すサンプルの数MTは、例えば観測対象として想定している目標物の大きさなどの情報から事前に決定することができる。
基底変換処理部109aは、目標検出部108において切り出された目標の受信散乱行列レンジプロフィールSm(m=M0,M0+1,・・・,M0+MT-1)に対し、まず式(12)で示す偏波基底変換あるいは式(13)に例示する偏波成分分解を行う(ステップST3)。
The number M T of samples to be cut out here can be determined in advance from information such as the size of a target assumed as an observation target.
First, the basis conversion processing unit 109a applies to the target reception scattering matrix range profile Sm (m = M 0 , M 0 +1,..., M 0 + M T −1) extracted by the target detection unit 108. Polarization basis conversion shown in equation (12) or polarization component decomposition exemplified in equation (13) is performed (step ST3).

Figure 2005156276
Figure 2005156276

ここで、Vはユニタリ行列であり、上付きの*Tは行列の共役転置を表す。
目的に応じて任意のユニタリ行列を設定できるものとする。
式(12)で表される偏波基底変換によって、第1偏波送受信アンテナ104と第2偏波送受信アンテナ105の直交する偏波の組み合わせとは異なる二つの直交偏波の組み合わせで同様な観測を行った場合の散乱行列レンジプロフィールUm(m=M0,M0+1,・・・,M0+MT-1)を求める。
Here, V is a unitary matrix, and the superscript * T represents the conjugate transpose of the matrix.
An arbitrary unitary matrix can be set according to the purpose.
The same observation is performed with a combination of two orthogonal polarizations different from the orthogonal polarization combination of the first polarization transmission / reception antenna 104 and the second polarization transmission / reception antenna 105 by the polarization basis conversion represented by the equation (12). Scattering matrix range profile Um (m = M 0 , M 0 +1,..., M 0 + M T −1) is obtained.

あるいは、例えば式(13)で示す偏波成分分解を用いて、受信波のうち、目標物において奇数回反射された後に送受信アンテナに到達した奇数回反射成分Ko(m)、目標物において偶数回反射された後に送受信アンテナに到達した偶数回反射成分Ke(m)、目標で反射される際に生じた、入射波と直交する偏波成分(交差成分)Kc(m)を求める。 Alternatively, for example, by using the polarization component decomposition expressed by Equation (13), the odd-numbered reflection component K o (m) that has reached the transmitting / receiving antenna after being reflected an odd number of times in the target of the received wave, and the even number in the target The even-numbered reflection component K e (m) that reaches the transmission / reception antenna after being reflected once, and the polarization component (cross component) K c (m) orthogonal to the incident wave that is generated when reflected by the target are obtained.

Figure 2005156276
Figure 2005156276

以上のように、種々の偏波基底変換や偏波成分分解を行うことにより、任意の種類の偏波チャネルについてレンジプロフィールを得ることが可能である。
ここで得られたレンジプロフィールの組を、以下では目標の複数偏波チャネルレンジプロフィールと呼ぶ。なお、ここでは、「偏波チャネル」は送信偏波と受信偏波の組み合わせだけでなく、Ko,Ke,Kcなどの偏波成分のことも含むこととする(以下同様)。
次に、特徴量算出部109bは、上記の基底変換処理によって得られた複数偏波チャネルレンジプロフィールから、以下に定義する特徴量を算出する(ステップST4)。
As described above, it is possible to obtain a range profile for any type of polarization channel by performing various polarization basis conversions and polarization component decomposition.
The set of range profiles obtained here is hereinafter referred to as a target multi-polarization channel range profile. Here, the “polarization channel” includes not only a combination of transmission polarization and reception polarization but also polarization components such as K o , K e , and K c (the same applies hereinafter).
Next, the feature quantity calculation unit 109b calculates a feature quantity defined below from the multiple polarization channel range profile obtained by the above-described base conversion process (step ST4).

Figure 2005156276
Figure 2005156276

式(14)において、P(m)はある偏波チャネルの信号の電力値を表しており、例えば、水平偏波送信水平偏波受信の偏波チャネルについては、式(15)のように定義される。rgは重心であり、式(16)のように定義される。Xpは偏波チャネルpのレンジプロフィールのフーリエスペクトルを表す。 In Expression (14), P (m) represents the power value of a signal of a certain polarization channel. For example, the polarization channel of horizontal polarization transmission horizontal polarization reception is defined as Expression (15). Is done. r g is the gravity center, is defined by the equation (16). Xp represents the Fourier spectrum of the range profile of the polarization channel p.

式(14)において、2次と3次のクロスモーメントスペクトルは、それぞれ偏波チャネル1と2と偏波チャネル1と2と3とで例示してあるが、基底変換処理部109aで得られた全ての偏波チャネルの組み合わせについて計算することが可能である。   In the equation (14), the second-order and third-order cross moment spectra are exemplified by the polarization channels 1 and 2, and the polarization channels 1, 2, and 3, respectively, but obtained by the basis conversion processing unit 109a. It is possible to calculate for all combinations of polarization channels.

Figure 2005156276
Figure 2005156276

図6は特徴量算出部で算出される特徴量のうち、平均、標準偏差、最大値および重心周りのモーメントの定義を示す説明図である。
特徴量は特徴量算出部109bで算出され、平均値は各偏波チャネルの平均的な反射強度、最大値は最大の反射強度、標準偏差は各偏波チャネルのレンジプロフィール形状の複雑さ、重心周りのモーメントは主にレンジプロフィールの長さを表す指標となっていることがわかる。
FIG. 6 is an explanatory diagram showing definitions of the average, standard deviation, maximum value, and moment around the center of gravity among the feature amounts calculated by the feature amount calculation unit.
The feature amount is calculated by the feature amount calculation unit 109b, the average value is the average reflection intensity of each polarization channel, the maximum value is the maximum reflection intensity, the standard deviation is the complexity of the range profile shape of each polarization channel, the center of gravity It can be seen that the surrounding moment is an index mainly representing the length of the range profile.

さらに、式(14)より、正規化された重心周りのモーメントとバイスペクトルはレンジプロフィール上での散乱強度の強い点の分布の様子、クロスモーメントスペクトルは複数偏波チャネルのレンジプロフィール間の関係をそれぞれ示す指標となっていると言える。以上をまとめると、式(14)に示した特徴量は、各偏波チャネルのレンジプロフィールの大まかな形状を示す指標と、複数偏波チャネルのレンジプロフィール間の関係を示す指標に大別できる。   Furthermore, from equation (14), the normalized moment around the center of gravity and bispectrum show the distribution of points with strong scattering intensity on the range profile, and the cross moment spectrum shows the relationship between the range profiles of multiple polarization channels. It can be said that each is an indicator. Summarizing the above, the feature quantity shown in Equation (14) can be broadly classified into an index indicating the rough shape of the range profile of each polarization channel and an index indicating the relationship between the range profiles of a plurality of polarization channels.

参照用目標特徴量ライブラリ109dには、同様にして事前に観測したN種類の目標の散乱行列レンジプロフィールより抽出された特徴量が格納されている。Smnは事前の観測の他に、理論計算(例えば、GTD:Geometrical Theory of Diffraction)などから得ることも可能である。   In the reference target feature amount library 109d, feature amounts extracted from N types of target scattering matrix range profiles observed in advance in the same manner are stored. Smn can be obtained from theoretical calculations (for example, GTD: Geometrical Theory of Diffraction) in addition to prior observations.

判定部109eは、参照用目標特徴量ライブラリ109dに蓄積されたN種類の目標に関する特徴量と比較を行い、あらかじめ定められた規則に従って目標の識別を行う(ステップST5)。   The determination unit 109e compares the N-type target feature quantities stored in the reference target feature quantity library 109d, and identifies the target according to a predetermined rule (step ST5).

例えば、識別に利用する特徴量の数がLの場合、各目標に関するL個の特徴量をL次元の特徴量ベクトルの形で参照用目標特徴量ライブラリ109dに蓄積しておき、受信信号から抽出したL個の特徴量からL次元の特徴量ベクトルを構成した後に、受信信号のL次元特徴量ベクトルと各参照用目標の特徴量ベクトルとのユークリッド距離をそれぞれ計算し、最も距離の短い参照用目標と同じ種別であると判定する方法が考えられる。
最後に、表示部110は、目標検出部108、判定部109eの出力を表示する。
For example, when the number of feature quantities used for identification is L, L feature quantities relating to each target are accumulated in the reference target feature quantity library 109d in the form of L-dimensional feature quantity vectors and extracted from the received signal. After constructing an L-dimensional feature quantity vector from the L feature quantities, the Euclidean distance between the L-dimensional feature quantity vector of the received signal and the feature quantity vector of each reference target is calculated, and the reference distance is the shortest. A method of determining that the type is the same as the target is conceivable.
Finally, the display unit 110 displays the outputs of the target detection unit 108 and the determination unit 109e.

この実施の形態1によれば、複数偏波チャネルのレンジプロフィール間の関係を示す複数種類の指標を算出し、該特徴量を組み合わせて目標識別を行う構成としたので、目標を正しく識別する確率を向上することができることと、目標の識別に散乱行列レンジプロフィールから、目標の距離方向の移動とアスペクト角の変化による変動の小さい特徴量として、各偏波チャネルのレンジプロフィールの大まかな形状を示す複数種類の指標を使用することで、ライブラリメモリ量、計算処理量を減らすことができるという効果がある。   According to the first embodiment, since a plurality of types of indexes indicating relationships between range profiles of a plurality of polarization channels are calculated, and the target identification is performed by combining the feature amounts, the probability of correctly identifying the target From the scattering matrix range profile for target identification, the rough shape of the range profile of each polarization channel is shown as a feature quantity with small variation due to movement of the target in the distance direction and change of the aspect angle By using a plurality of types of indexes, there is an effect that the amount of library memory and the amount of calculation processing can be reduced.

実施の形態2.
図7はこの発明の実施の形態2による目標識別用レーダ装置の構成を示すブロック図であり、109fは特徴量抽出部109cにおいて抽出された特徴量に、事前の観測結果に基づいて定めた重み付け係数をかける重み付け部である。
Embodiment 2. FIG.
FIG. 7 is a block diagram showing the configuration of the target identifying radar apparatus according to Embodiment 2 of the present invention, and 109f is a weight determined on the basis of a previous observation result for the feature quantity extracted by the feature quantity extraction unit 109c. It is a weighting unit for applying a coefficient.

一般に識別の対象とする目標物によって、角度依存性の低下に有効な特徴量の組み合わせは異なることが予想されるため、式(14)で示した特徴量Ci(i=1,2,・・・,L:Lは特徴量の総数)の中から適当なものを選択して組み合わせる、あるいは適当な重み付けを行って判別することが望ましい。 In general, it is expected that the combination of feature quantities effective in reducing the angle dependency varies depending on the target to be identified. Therefore, the feature quantities C i (i = 1, 2,. .., L: It is desirable to select and combine appropriate ones from among L (the total number of feature values), or perform appropriate weighting.

識別の対象とする目標が既知である場合は、事前に観測を行って特徴量を抽出する際に、最も正しく識別する確率が高くなる重み付けの係数wi(i=1,2,・・・,L)を決定する。重み付け部109fは特徴量抽出部109cにおいて抽出された特徴量に対し、次式に従って重み付けを行う。なお、重み付け係数の値が0であれば、対応する特徴量を選択しないことを意味する。 If the target to be identified is known, the weighting coefficient w i (i = 1, 2,...) That increases the probability of the most correct identification when performing feature extraction by observation in advance. , L). The weighting unit 109f weights the feature amount extracted by the feature amount extraction unit 109c according to the following equation. If the value of the weighting coefficient is 0, it means that the corresponding feature amount is not selected.

なおここで、特徴量とは、各偏波チャネル信号の平均電力、各偏波チャネル信号の標準偏差、各偏波チャネル信号の最大値、各偏波チャネル信号の重心周りのモーメント、各偏波チャネル信号の重心周りの正規化されたモーメント、各偏波チャネル信号の高次スペクトル、上記複数の偏波チャネル信号のクロスモーメントスペクトル等を表す。   Here, the feature amount means the average power of each polarization channel signal, the standard deviation of each polarization channel signal, the maximum value of each polarization channel signal, the moment around the center of gravity of each polarization channel signal, and each polarization. It represents a normalized moment around the center of gravity of the channel signal, a higher-order spectrum of each polarization channel signal, a cross moment spectrum of the plurality of polarization channel signals, and the like.

Figure 2005156276
Figure 2005156276

なお、重み付け係数wi(i=1,2,・・・,L)の組み合わせを、対象とする目標の種別に応じて複数用意しておき、動作時には、オペレータが適した重み付け係数の組み合わせを選択することによって、汎用性を確保することができる。 It should be noted that a plurality of combinations of weighting coefficients w i (i = 1, 2,..., L) are prepared in accordance with the target type to be processed, and an operator selects a combination of weighting coefficients suitable for the operation. By selecting, versatility can be ensured.

この実施の形態2によれば、識別の対象とする目標物によって、適当な重み付けを特徴量に加えることができるため、目標を正しく識別する確率を向上することができる効果を奏する。   According to the second embodiment, since an appropriate weight can be added to the feature amount depending on the target to be identified, there is an effect that the probability of correctly identifying the target can be improved.

実施の形態3.
図8はこの発明の実施の形態3による目標識別用レーダ装置の構成を示すブロック図であり、109gは各偏波チャネルの信号対雑音電力比に応じて識別に利用する偏波チャネルを選択する偏波チャネル選択部である。
Embodiment 3 FIG.
FIG. 8 is a block diagram showing a configuration of a target identifying radar apparatus according to Embodiment 3 of the present invention. 109g selects a polarization channel to be used for identification according to the signal-to-noise power ratio of each polarization channel. It is a polarization channel selector.

偏波チャネル選択部109gでは、各偏波チャネルの平均電力を求めた後、これを既知の受信機雑音レベルと比較して、平均電力が受信機雑音レベルのα倍よりも低い偏波チャネルのデータを、特徴量算出部109bへ送らないようにする。   The polarization channel selection unit 109g obtains the average power of each polarization channel, and then compares this with the known receiver noise level, and the average power of the polarization channel whose average power is lower than α times the receiver noise level. The data is not sent to the feature amount calculation unit 109b.

この実施の形態3によれば、信号対雑音比の低い偏波チャネルのデータを排除できるため、雑音の影響による目標の誤識別確率を低減する効果を奏する。   According to the third embodiment, since data of a polarization channel having a low signal-to-noise ratio can be excluded, there is an effect of reducing the target misidentification probability due to the influence of noise.

実施の形態4.
図9、図10、図11はこの発明の実施の形態4による目標識別用レーダ装置の構成を示すブロック図である。
Embodiment 4 FIG.
9, FIG. 10 and FIG. 11 are block diagrams showing the configuration of a target identifying radar apparatus according to Embodiment 4 of the present invention.

図9、図10、図11において、109hは複数命中の観測値が得られた場合に、各命中信号における識別結果を比較して、最も多い結果を最終的な識別結果として出力する多数決部、109iは複数命中の観測値が得られた場合に、各命中信号で得られた散乱行列レンジプロフィールの平均処理を行うレンジプロフィール平均処理部、109jは複数命中の観測値が得られた場合に、特徴量抽出部109cにおいて命中毎に抽出された特徴量から特徴量毎に平均値を算出する特徴量平均処理部であり、111は複数ヒット観測値が得られた場合に、各命中信号における散乱行列レンジプロフィールを蓄積しておく複数命中の受信偏波信号記憶部である。   9, 10, and 11, when a plurality of hit observation values are obtained, 109h compares the discrimination results in the respective hit signals, and outputs the most result as the final discrimination result, 109i is a range profile average processing unit that averages the scattering matrix range profile obtained with each hit signal when multiple hit observations are obtained, and 109j is a multiple hit observation value obtained. A feature amount average processing unit that calculates an average value for each feature amount from the feature amounts extracted for each hit in the feature amount extraction unit 109c. 111 is a scatter in each hit signal when multiple hit observation values are obtained. It is a received polarization signal storage unit for multiple hits that stores matrix range profiles.

複数命中の観測値が得られた場合、大別して二通りの方法が可能である。判定部109eにおいて命中毎に判定を行った後、もっとも妥当と思われる判定結果を選択する方法と、判定部109eへの入力値を予め平均化する方法である。
ここで、図9は前者、図10と図11は後者の方法をそれぞれ具体化した構成を示している。
When multiple hits are obtained, there are two main methods. This is a method of selecting a determination result that seems to be most appropriate after determining each hit in the determination unit 109e and a method of averaging the input values to the determination unit 109e in advance.
Here, FIG. 9 shows a configuration embodying the former, and FIGS. 10 and 11 show configurations embodying the latter method.

まず、図9の構成において多数決部109hの動作を説明する。観測に用いた命中数がHの場合、判定部109eはH個の目標識別結果を出力する。
これらの識別結果が全て等しくない場合に、多数決部109hは最も数の多かったものを最終的な目標識別結果として出力する。あるいは、事前に定めた規則に従って、識別結果と、その結果に対する確度を同時に出力するようにしてもよい。
First, the operation of the majority decision unit 109h in the configuration of FIG. 9 will be described. When the number of hits used for observation is H, the determination unit 109e outputs H target identification results.
When these identification results are not all equal, the majority decision unit 109h outputs the largest number as the final target identification result. Alternatively, according to a predetermined rule, the identification result and the accuracy for the result may be output at the same time.

例えば、10命中の観測を行った場合に、判定部109eの判定結果が、目標A5回、目標B3回、目標C2回の場合に、多数決部109hは『目標はAである』と出力してもよいし、『目標Aである確率50%、目標Bである確率30%、目標Cである確率20%』と出力してもよい。   For example, when 10 hits are observed, if the determination result of the determination unit 109e is the target A5 times, the target B3 times, and the target C2 times, the majority decision unit 109h outputs “the target is A” Alternatively, “50% probability of target A, 30% probability of target B, 20% probability of target C” may be output.

次に図10の構成において、レンジプロフィール平均処理部109iの動作を説明する。
観測に用いた命中数がHの場合、複数ヒット受信偏波信号記憶部111はH個の散乱行列レンジプロフィールSm,h(m=1,2,・・・,M;h=1,2,・・・,H)を蓄積する。
Next, the operation of the range profile average processing unit 109i in the configuration of FIG. 10 will be described.
When the hit number used for the observation is H, the multiple hit received polarization signal storage unit 111 stores H scattering matrix range profiles Sm, h (m = 1, 2,..., M; h = 1, 2, ..., H) are accumulated.

レンジプロフィール平均処理部109iは、複数ヒット受信偏波信号記憶部111に蓄積されたこれらの散乱行列レンジプロフィールを次式に従って平均化し、得られた平均化散乱行列レンジプロフィールを特徴量抽出部109cに送る。   The range profile average processing unit 109i averages these scattering matrix range profiles accumulated in the multiple hit received polarization signal storage unit 111 according to the following equation, and the obtained averaged scattering matrix range profile is sent to the feature amount extraction unit 109c. send.

Figure 2005156276
Figure 2005156276

ここで、θh(h=1,2,・・・,H)は、プラットフォームの移動による位相の変化量を補正するための補正項である。
次に図11の構成において、特徴量平均処理部109jの動作を説明する。特徴量平均処理部109jは、特徴量抽出部109cにおいて得られた、命中毎のL個の特徴量Ci,h(i=1,2,・・・,L:Lは特徴量の総数; h=1,2,・・・,H)を次式に従って平均化し、得られた平均化特徴量を判定部109eに送る。
Here, θ h (h = 1, 2,..., H) is a correction term for correcting the amount of phase change due to the movement of the platform.
Next, the operation of the feature amount average processing unit 109j in the configuration of FIG. 11 will be described. The feature quantity average processing unit 109j obtains L feature quantities C i, h (i = 1, 2,..., L: L is the total number of feature quantities) for each hit obtained by the feature quantity extraction unit 109c; h = 1, 2,..., H) are averaged according to the following equation, and the obtained average feature amount is sent to the determination unit 109e.

Figure 2005156276
Figure 2005156276

この実施の形態4によれば、複数回観測を行うことでデータ量を増加することができるため、目標を正しく識別する確率を向上できる効果を奏する。   According to the fourth embodiment, since the amount of data can be increased by performing observation a plurality of times, there is an effect that the probability of correctly identifying the target can be improved.

実施の形態5.
図12はこの発明の実施の形態5による目標識別用レーダ装置の構成を示すブロック図であり、115は目標移動軌跡算出部であり、1〜110は実施の形態1の説明と同じものであり、実施の形態1と同じ部分についての説明は実施の形態1の説明によるものとする。
Embodiment 5 FIG.
FIG. 12 is a block diagram showing the configuration of the target identifying radar apparatus according to the fifth embodiment of the present invention. Reference numeral 115 denotes a target movement trajectory calculation unit. Reference numerals 1-110 are the same as those described in the first embodiment. The description of the same parts as those of the first embodiment is based on the description of the first embodiment.

目標移動軌跡算出部115は目標検出部108によって受信偏波信号から目標を検出した信号を受け、時々刻々と変化する目標の移動ベクトルを航跡として捕らえ目標の姿勢を算出し参照すべき目標の特徴量を絞り込み、一方109は目標検出部108によって切り出された散乱行列レンジプロフィールから複数の特徴量を算出し、ここで得られた特徴量を、事前に蓄積された既知の目標に関する特徴量のライブラリと比較することで目標を識別する目標識別部を持っている。   The target movement trajectory calculation unit 115 receives the signal obtained by detecting the target from the received polarization signal by the target detection unit 108, captures the target movement vector that changes every moment as a wake, calculates the attitude of the target, and refers to the target characteristics On the other hand, 109 calculates a plurality of feature amounts from the scattering matrix range profile extracted by the target detection unit 108, and the feature amounts obtained here are a library of feature amounts relating to known targets accumulated in advance. It has a target identification unit that identifies the target by comparing with.

例えば、識別に利用する特徴量の数がLの場合、各目標に関するL個の特徴量をL次元の特徴量ベクトルの形で参照用目標特徴量ライブラリ109dに蓄積しておき、受信信号から抽出したL個の特徴量からL次元の特徴量ベクトルを構成した後に、受信信号のL次元特徴量ベクトルと各参照用目標の特徴量ベクトルとのユークリッド距離をそれぞれ計算し、最も距離の短い参照用目標と同じ種別であると判定する方法が考えられる。
最後に、表示部110は、目標検出部108、判定部109eの出力を表示する。
For example, when the number of feature quantities used for identification is L, L feature quantities relating to each target are accumulated in the reference target feature quantity library 109d in the form of L-dimensional feature quantity vectors and extracted from the received signal. After constructing an L-dimensional feature quantity vector from the L feature quantities, the Euclidean distance between the L-dimensional feature quantity vector of the received signal and the feature quantity vector of each reference target is calculated, and the reference distance is the shortest. A method of determining that the type is the same as the target is conceivable.
Finally, the display unit 110 displays the outputs of the target detection unit 108 and the determination unit 109e.

なおここで、特徴量とは各偏波チャネル信号の平均電力、各偏波チャネル信号の標準偏差、各偏波チャネル信号の最大値、各偏波チャネル信号の重心周りのモーメント等である。   Here, the feature amount includes an average power of each polarization channel signal, a standard deviation of each polarization channel signal, a maximum value of each polarization channel signal, a moment around the center of gravity of each polarization channel signal, and the like.

この実施の形態5によれば、目標の識別に散乱行列レンジプロフィールから、特徴量として、各偏波チャネルのレンジプロフィールの大まかな形状を示す複数種類の指標と、複数の偏波チャネルのレンジプロフィール間の関係を示す複数種類の指標を算出し、該特徴量を組み合わせて目標識別を行う構成としたので、目標を正しく識別する確率を向上することができる効果を奏する。   According to the fifth embodiment, from the scattering matrix range profile for target identification, a plurality of types of indicators indicating the rough shape of the range profile of each polarization channel, and the range profiles of the plurality of polarization channels are used as features. Since a plurality of types of indices indicating the relationship between them are calculated and the target identification is performed by combining the feature amounts, the probability of correctly identifying the target can be improved.

実施の形態6.
図13はこの発明の実施の形態6による目標識別用レーダ装置の構成を示すブロック図であり、112は目標検出部108によって目標の受信偏波信号を偏波特性として画像化した後、目標の姿勢を算出し参照すべき目標の特徴量を絞り込み、一方109は目標検出部108によって切り出された散乱行列レンジプロフィールから複数の特徴量を算出し、ここで得られた特徴量を、事前に蓄積された既知の目標に関する特徴量のライブラリと比較することで目標を識別するものである。
Embodiment 6 FIG.
FIG. 13 is a block diagram showing the configuration of a target identifying radar apparatus according to Embodiment 6 of the present invention. Reference numeral 112 shows a target received polarization signal imaged as a polarization characteristic by the target detector 108, and then the target is detected. The target feature quantity to be referred to is narrowed down, while 109 calculates a plurality of feature quantities from the scattering matrix range profile extracted by the target detection unit 108, and the obtained feature quantities are calculated in advance. The target is identified by comparing it with a library of feature quantities related to the accumulated known target.

この実施の形態6によれば、識別の対象とする目標の偏波特性を画像化するため、対象物そのものの形状から姿勢角を把握できるため、より適切な重み付けを特徴量に加えること、もしくは選定できるため、目標を正しく識別する確率を向上することができる効果を奏する。   According to the sixth embodiment, in order to image the target polarization characteristics to be identified, the posture angle can be grasped from the shape of the object itself, so that more appropriate weighting is added to the feature amount. Or since it can select, the effect which can improve the probability which identifies a target correctly is produced.

実施の形態7.
図14はこの発明の実施の形態7による目標識別用レーダ装置の構成を示すブロック図であり、113は目標検出部108によって受信偏波信号を用いて目標を検出した後、検出された近傍のレンジプロファイルを切り出した結果から、目標のRCS(レーダ反射面積)データを算出し、そのRCSデータを元に、事前の観測結果に基づいて定めたRCSデータとパターンマッティング方式により目標の姿勢を求め参照すべき目標の特徴量を絞り込み、一方109は目標検出部108によって切り出された散乱行列レンジプロフィールから複数の特徴量を算出し、ここで得られた特徴量を、事前に蓄積された既知の目標に関する特徴量のライブラリと比較することで目標を識別する。
Embodiment 7 FIG.
FIG. 14 is a block diagram showing the configuration of a target identifying radar apparatus according to Embodiment 7 of the present invention. Reference numeral 113 denotes a target detected by the target detecting unit 108 using the received polarization signal, and the detected neighborhood. The target RCS (radar reflection area) data is calculated from the result of cutting out the range profile, and based on the RCS data, the target posture is obtained by the RCS data determined based on the previous observation result and the pattern matting method. The target feature amount to be referred to is narrowed down, while 109 calculates a plurality of feature amounts from the scattering matrix range profile cut out by the target detection unit 108, and the obtained feature amounts are stored in advance as known features. The target is identified by comparing it with a library of feature values for the target.

この実施の形態7によれば、目標の姿勢を算出するための参照データは容易に確保できるため、計算負荷は低く、あらゆる目標に対して、特徴量を用いた目標識別の確率を向上することができる効果を奏する。   According to the seventh embodiment, since reference data for calculating a target posture can be easily secured, the calculation load is low, and the probability of target identification using feature quantities is improved for every target. There is an effect that can.

実施の形態8.
図15はこの発明の実施の形態8による目標識別用レーダ装置の構成を示すブロック図であり、109は目標検出部108によって切り出された散乱行列レンジプロフィールから複数の特徴量を算出する目標識別部であり、一方、114は基底変換処理部109aによって切り出された受信偏波信号に偏波基底変換処理を施し、複数の偏波チャンネル信号を得て、1回反射成分、2回反射成分及び多重反射成分の3種類の基本構成に分解し、それぞれの寄与率から目標の姿勢を算出し、参照すべき目標の特徴量を絞り込み、目標識別部109で得られた特徴量を、事前に蓄積された既知の目標に関する特徴量のライブラリと比較することで目標を識別する基本構成3成分寄与率算出部である。
Embodiment 8 FIG.
FIG. 15 is a block diagram showing the configuration of a target identification radar apparatus according to Embodiment 8 of the present invention. Reference numeral 109 denotes a target identification unit for calculating a plurality of feature amounts from the scattering matrix range profile extracted by the target detection unit 108. On the other hand, 114 performs polarization basis conversion processing on the received polarization signal cut out by the basis conversion processing unit 109a to obtain a plurality of polarization channel signals to obtain a one-time reflection component, two-time reflection component, and multiplexing. It is decomposed into three basic components of the reflection component, the target posture is calculated from each contribution rate, the target feature amount to be referred to is narrowed down, and the feature amount obtained by the target identification unit 109 is accumulated in advance. The basic component three-component contribution rate calculation unit for identifying a target by comparing it with a library of feature amounts related to a known target.

この実施の形態8によれば、識別の対象とする目標の偏波特性の精度向上がはかれ、その結果が姿勢角算出精度に反映されるため、より適切な重み付けを特徴量に加えること、もしくは選定できるため、目標を正しく識別する確率を向上することができる効果を奏する。   According to the eighth embodiment, the accuracy of the target polarization characteristic to be identified is improved, and the result is reflected in the attitude angle calculation accuracy. Therefore, more appropriate weighting is added to the feature amount. Or, since it can be selected, it is possible to improve the probability of correctly identifying the target.

この発明の実施の形態1による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 1 of this invention. この発明の実施の形態1による目標識別方法を示すフローチャートである。It is a flowchart which shows the target identification method by Embodiment 1 of this invention. レンジプロフィールの概念を示す説明図である。It is explanatory drawing which shows the concept of a range profile. 散乱行列レンジプロフィールの概念を示す説明図である。It is explanatory drawing which shows the concept of a scattering matrix range profile. 第1偏波送受信アンテナと第2偏波送受信アンテナの各時刻における動作モードを示す説明図である。It is explanatory drawing which shows the operation mode in each time of a 1st polarized wave transmission / reception antenna and a 2nd polarization transmission / reception antenna. 特徴量算出部で算出される特徴量のうち、平均、標準偏差、最大値および重心周りのモーメントの定義を示す説明図である。It is explanatory drawing which shows the definition of an average, a standard deviation, a maximum value, and the moment around a gravity center among the feature-values calculated in a feature-value calculation part. この発明の実施の形態2による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 2 of this invention. この発明の実施の形態3による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 3 of this invention. この発明の実施の形態4による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 4 of this invention. この発明の実施の形態4による他の目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for another target identification by Embodiment 4 of this invention. この発明の実施の形態4による図9または図10以外の目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification other than FIG. 9 or FIG. 10 by Embodiment 4 of this invention. この発明の実施の形態5による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 5 of this invention. この発明の実施の形態6による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 6 of this invention. この発明の実施の形態7による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 7 of this invention. この発明の実施の形態8による目標識別用レーダ装置の構成を示すブロック図である。It is a block diagram which shows the structure of the radar apparatus for target identification by Embodiment 8 of this invention.

符号の説明Explanation of symbols

101 送信機、 102 送受切換器、 103 偏波切換器、 104 第1偏波送受信アンテナ、 105 第2偏波送受信アンテナ、 106 受信機、 107 受信偏波信号記憶部、 108 目標検出部、 109 目標識別部、 109a 基底変換処理部、 109b 特徴量算出部、 109c 特徴量抽出部、 109d 参照用目標特徴量ライブラリ、 109e 判定部、 109f 重み付け部、 109g 偏波チャネル選択部、 109h 多数決部、 109i レンジプロフィール平均処理部、 109j 特徴量平均処理部、 110 表示部、 111 複数ヒット受信偏波信号記憶部、 112 目標画像生成部、 113・RCS算出部、 114・基本構成3成分寄与率算出部、 115 目標移動航跡算出部。   DESCRIPTION OF SYMBOLS 101 Transmitter, 102 Transmission / reception switch, 103 Polarization switch, 104 1st polarization transmission / reception antenna, 105 2nd polarization transmission / reception antenna, 106 receiver, 107 Reception polarization signal memory | storage part, 108 Target detection part, 109 Target Identification unit, 109a basis conversion processing unit, 109b feature quantity calculation unit, 109c feature quantity extraction unit, 109d reference target feature quantity library, 109e determination unit, 109f weighting unit, 109g polarization channel selection unit, 109h majority decision unit, 109i range Profile average processing unit, 109j feature quantity average processing unit, 110 display unit, 111 multiple hit received polarization signal storage unit, 112 target image generation unit, 113 · RCS calculation unit, 114 · basic configuration three component contribution rate calculation unit, 115 Target movement track calculation unit.

Claims (13)

広帯域パルスを発生して送信する送信機と、
相互に偏波特性の異なる複数のアンテナを用いて複数の偏波信号を受信する受信機と、上記受信機によって受信した複数の受信偏波信号から目標を検出する目標検出部と、
上記目標検出部からの信号に基づき特徴量を算出する特徴量抽出部と、
上記特徴量抽出部によって抽出された特徴量を用いて目標の識別を行う目標識別部と
を備えた目標識別用レーダ装置。
A transmitter that generates and transmits wideband pulses; and
A receiver that receives a plurality of polarization signals using a plurality of antennas having different polarization characteristics from each other; a target detection unit that detects a target from a plurality of received polarization signals received by the receiver;
A feature amount extraction unit that calculates a feature amount based on a signal from the target detection unit;
A target identifying radar apparatus comprising: a target identifying unit that identifies a target using the feature amount extracted by the feature amount extracting unit.
上記特徴量抽出部として、受信した複数の受信偏波信号から複数の偏波チャネル信号を算出し、特徴量を算出することを特徴とする請求項1記載の目標識別用レーダ装置。 2. The target identifying radar apparatus according to claim 1, wherein the feature quantity extraction unit calculates a feature quantity by calculating a plurality of polarization channel signals from a plurality of received polarization signals received. 上記特徴量抽出部として、目標の距離方向の移動とアスペクト角の変化による変動の小さい特徴量を算出することを特徴とする請求項1または請求項2記載の目標識別用レーダ装置。 3. The target identifying radar apparatus according to claim 1, wherein the feature amount extracting unit calculates a feature amount that is small in variation due to movement of the target in the distance direction and change in the aspect angle. 上記特徴量抽出部によって抽出された特徴量を、事前に蓄積された既知の目標に関する特徴量のライブラリと比較することで目標を識別することを特徴とする請求項1〜請求項3の何れか1項記載の目標識別用レーダ装置。 4. The target is identified by comparing the feature quantity extracted by the feature quantity extraction unit with a library of feature quantities relating to known targets accumulated in advance. The radar device for target identification according to 1. 上記目標識別部として、事前の観測結果に基づいて決定された上記特徴量の信頼度に基づいて、各特徴量に適当な重みづけを加えた後、目標の識別を行うことを特徴とする請求項1〜請求項4の何れか1項記載の目標識別用レーダ装置。 The target identification unit is configured to identify a target after adding an appropriate weight to each feature amount based on the reliability of the feature amount determined based on a previous observation result as the target identification unit. The radar device for target identification according to any one of claims 1 to 4. 上記目標識別部として、事前に既知の雑音電力に対して平均電力が一定レベル以下の信号を除去した偏波チャネルの信号を用いて目標の識別を行うことを特徴とする請求項1〜請求項5の何れか1項記載の目標識別用レーダ装置。 The target identification is performed by using a signal of a polarization channel obtained by removing a signal whose average power is equal to or lower than a certain level with respect to known noise power in advance as the target identification unit. The radar device for target identification according to any one of 5. 上記特徴量抽出部として、受信した複数命中信号から上記特徴量を算出すること、および、上記目標識別部は、命中信号毎に上記特徴量抽出部によって算出された特徴量の全てあるいは一部を用いて、目標の識別を行った後、命中信号毎の識別結果を集計して最も多かった識別結果を最終的な識別結果として出力する多数決部を有することを特徴とする請求項1〜請求項6の何れか1項記載の目標識別用レーダ装置。 The feature amount extraction unit calculates the feature amount from the received multiple hit signals, and the target identification unit obtains all or part of the feature amounts calculated by the feature amount extraction unit for each hit signal. Claims 1 to claim, further comprising: a majority voting unit that aggregates the identification results for each hit signal and outputs the most frequent identification result as a final identification result after identifying the target. The radar device for target identification according to any one of 6. 上記特徴量抽出部として、受信した複数命中信号を命中方向に平均して平均受信偏波信号を算出した後、該平均受信偏波信号から上記特徴量を算出することを特徴とする請求項1〜請求項7の何れか1項記載の目標識別用レーダ装置。 2. The feature quantity extracting unit, after averaging a plurality of received hit signals in the hit direction to calculate an average received polarization signal, calculating the feature quantity from the average received polarization signal. The target identifying radar apparatus according to claim 7. 上記特徴量抽出部として、受信した複数命中信号から上記特徴量を算出した後、各特徴量を命中信号方向に平均化することを特徴とする請求項1〜請求項8の何れか1項記載の目標識別用レーダ装置。 The said feature-value extraction part calculates the said feature-value from the received multiple hit signal, Then, averages each feature-value in the hit signal direction, The any one of Claims 1-8 characterized by the above-mentioned. Target identification radar device. 検出した目標の移動航跡を算出し目標の姿勢を決定する目標移動航跡算出部を具備し、
上記特徴量抽出部として、受信した複数の受信偏波信号から複数の偏波チャネル信号を算出し、特徴量を算出することを特徴とする請求項1記載の目標識別用レーダ装置。
A target movement track calculation unit for calculating a movement track of the detected target and determining a posture of the target;
2. The target identifying radar apparatus according to claim 1, wherein the feature quantity extraction unit calculates a feature quantity by calculating a plurality of polarization channel signals from a plurality of received polarization signals received.
上記目標識別部として、受信した複数の受信偏波信号を画像化し、目標の姿勢を算出する目標画像生成部を備え、事前の観測結果に基づいて決定された特徴量に対して適切な比較特徴量を取捨選択し、目標の識別能力を向上することを特徴とする請求項1および請求項3〜請求項10の何れか1項記載の目標識別用レーダ装置。 The target identification unit includes a target image generation unit that images a plurality of received polarization signals received and calculates a target attitude, and is a comparative feature suitable for a feature amount determined based on a previous observation result 11. The target identifying radar apparatus according to claim 1, wherein the target identifying ability is improved by selecting a quantity. 上記目標識別部として、受信した複数の受信偏波信号からRCS(レーダ反射面積)を算出し、目標の姿勢を決定するRCS算出部を備え、事前の観測結果に基づいて決定された特徴量に対して適切な比較特徴量を取捨選択し、目標の識別能力を向上することを特徴とする請求項1および請求項3〜請求項11の何れか1項記載の目標識別用レーダ装置。 As the target identification unit, an RCS (radar reflection area) is calculated from a plurality of received polarization signals received, and an RCS calculation unit for determining a target attitude is provided. The target identification radar apparatus according to any one of claims 1 and 3 to 11, wherein an appropriate comparison feature amount is selected to improve target identification capability. 上記目標識別部として、受信した複数の受信偏波信号から基底変換処理を行い、3種類の基本構成の成分に分解し、且つそれぞれの寄与率を算出し目標の姿勢を算出する3種類の基本構成の成分の寄与率を算出する基本構成3成分寄与率算出部を備え、事前の観測結果に基づいて決定された特徴量に対して適切な比較特徴量を取捨選択し、目標の識別能力を向上することを特徴とする請求項1〜請求項12の何れか1項記載の目標識別用レーダ装置。 As the target identification unit, three types of basics for performing base conversion processing from a plurality of received polarization signals received, decomposing the components into three types of basic components, calculating respective contribution rates, and calculating a target attitude A basic configuration 3 component contribution rate calculation unit that calculates the contribution rate of the component of the configuration is selected, and appropriate comparison feature amounts are selected with respect to the feature amounts determined based on the previous observation results, and the target discrimination ability is increased. 13. The target identifying radar apparatus according to claim 1, wherein the target identifying radar apparatus is improved.
JP2003393328A 2003-11-25 2003-11-25 Radar device for target identification Expired - Fee Related JP4096867B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003393328A JP4096867B2 (en) 2003-11-25 2003-11-25 Radar device for target identification

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003393328A JP4096867B2 (en) 2003-11-25 2003-11-25 Radar device for target identification

Publications (2)

Publication Number Publication Date
JP2005156276A true JP2005156276A (en) 2005-06-16
JP4096867B2 JP4096867B2 (en) 2008-06-04

Family

ID=34719721

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003393328A Expired - Fee Related JP4096867B2 (en) 2003-11-25 2003-11-25 Radar device for target identification

Country Status (1)

Country Link
JP (1) JP4096867B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007218689A (en) * 2006-02-15 2007-08-30 Mitsubishi Electric Corp Target classifier
JP2008002818A (en) * 2006-06-20 2008-01-10 Alpine Electronics Inc Vehicle identification system
JP2010230462A (en) * 2009-03-26 2010-10-14 Nec Corp Polarimetric sar system and target identification method using polarimetric sar data
JP2011520125A (en) * 2008-05-07 2011-07-14 コロラド ステート ユニバーシティー リサーチ ファウンデーション Dual Polarization Radar Processing System Using Time Domain Method
JP2019519768A (en) * 2016-05-27 2019-07-11 ロンバス システムズ グループ, インコーポレイテッドRhombus Systems Group, Inc. Radar system for tracking low-flying unmanned aerial vehicles and objects

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007218689A (en) * 2006-02-15 2007-08-30 Mitsubishi Electric Corp Target classifier
JP2008002818A (en) * 2006-06-20 2008-01-10 Alpine Electronics Inc Vehicle identification system
JP2011520125A (en) * 2008-05-07 2011-07-14 コロラド ステート ユニバーシティー リサーチ ファウンデーション Dual Polarization Radar Processing System Using Time Domain Method
US8665144B2 (en) 2008-05-07 2014-03-04 Colorado State University Research Foundation Dual-polarization radar processing system using time domain method
JP2010230462A (en) * 2009-03-26 2010-10-14 Nec Corp Polarimetric sar system and target identification method using polarimetric sar data
JP2019519768A (en) * 2016-05-27 2019-07-11 ロンバス システムズ グループ, インコーポレイテッドRhombus Systems Group, Inc. Radar system for tracking low-flying unmanned aerial vehicles and objects
JP7036744B2 (en) 2016-05-27 2022-03-15 ロンバス システムズ グループ, インコーポレイテッド Radar system for tracking low-flying unmanned aerial vehicles and objects
US11294048B2 (en) 2016-05-27 2022-04-05 Rhombus Systems Group, Inc. Radar system to track low flying unmanned aerial vehicles and objects
US11656354B2 (en) 2016-05-27 2023-05-23 Rhombus Systems Group, Inc. Radar system to track low flying unmanned aerial vehicles and objects

Also Published As

Publication number Publication date
JP4096867B2 (en) 2008-06-04

Similar Documents

Publication Publication Date Title
Conte et al. Statistical analysis of real clutter at different range resolutions
US20090021422A1 (en) Device and method for estimating the number of arrival signals
CN104155648B (en) The high-frequency ground wave radar single snap MUSIC direction-finding method reset based on array data
US20120313809A1 (en) Signal processing unit and method
CA2411615A1 (en) Surface wave radar
JP4972852B2 (en) Radar equipment
US20050242995A1 (en) Method for the blind wideband localization of one or more transmitters from a carrier that is passing by
JP2010038744A (en) Radar system for target identification
JP4096867B2 (en) Radar device for target identification
JP2010286404A (en) Moving target detector
JP5093781B2 (en) Radar equipment
JP2010175383A (en) Target detection apparatus and target detection method
JP2005291816A (en) Radar system
KR102099388B1 (en) Method of estimating direction of arrival of radar signal based on antenna array extrapolation and apparatus for the same
JP2008164494A (en) Guidance device
JP2014153088A (en) Rader system and tracking processor
JP4698941B2 (en) Target classification device
JP2005189171A (en) Radar apparatus
JP2012237705A (en) Azimuth and elevation angle measuring system
JP3659577B2 (en) Clutter suppression device and clutter suppression method
JP3905369B2 (en) Target detection apparatus and method
CN110940977B (en) Constant false alarm detector adaptive to sea condition change and constant false alarm detection method
JP4275425B2 (en) Radar equipment
JP5724548B2 (en) Direction measuring apparatus, direction measuring method and direction measuring program
JP3783842B2 (en) Target detection apparatus and target detection method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070821

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20071022

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20071120

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080121

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080303

R151 Written notification of patent or utility model registration

Ref document number: 4096867

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

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

Free format text: PAYMENT UNTIL: 20110321

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20110321

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20120321

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20130321

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20130321

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20140321

Year of fee payment: 6

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

LAPS Cancellation because of no payment of annual fees