JP6363524B2 - Radar apparatus and radar signal processing method - Google Patents

Radar apparatus and radar signal processing method Download PDF

Info

Publication number
JP6363524B2
JP6363524B2 JP2015017275A JP2015017275A JP6363524B2 JP 6363524 B2 JP6363524 B2 JP 6363524B2 JP 2015017275 A JP2015017275 A JP 2015017275A JP 2015017275 A JP2015017275 A JP 2015017275A JP 6363524 B2 JP6363524 B2 JP 6363524B2
Authority
JP
Japan
Prior art keywords
range
doppler
target
δaz
δel
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.)
Active
Application number
JP2015017275A
Other languages
Japanese (ja)
Other versions
JP2016142578A (en
Inventor
晋一 竹谷
晋一 竹谷
正一郎 安達
正一郎 安達
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.)
Toshiba Corp
Toshiba Infrastructure Systems and Solutions Corp
Original Assignee
Toshiba Corp
Toshiba Infrastructure Systems and Solutions 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 Toshiba Corp, Toshiba Infrastructure Systems and Solutions Corp filed Critical Toshiba Corp
Priority to JP2015017275A priority Critical patent/JP6363524B2/en
Publication of JP2016142578A publication Critical patent/JP2016142578A/en
Application granted granted Critical
Publication of JP6363524B2 publication Critical patent/JP6363524B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Description

本実施形態は、レーダ装置及びレーダ信号処理方法に関する。   The present embodiment relates to a radar apparatus and a radar signal processing method.

レーダ装置にあっては、画像により目標を識別する方法として、目標の重心等をレンジ及びドップラの両軸で追跡して画像中心を取得し、レンジ圧縮及びAZ圧縮することにより、レンジ−ドップラ軸で目標を画像化するISAR(Inverse Synthetic Aperture Radar:逆合成開口レーダ)処理が知られている(非特許文献1参照)。従来のISAR処理では、剛体目標が移動している場合、その目標の回転運動により、目標上の各反射点におけるドップラ速度が異なる。このため、ISAR処理によって生成した画像が真の目標形状と異なり、形状寸法の絶対値とも異なるのが通常である。   In a radar apparatus, as a method for identifying a target by an image, the center of gravity of the target is tracked by both the range and Doppler axes, the center of the image is acquired, and the range compression and the AZ compression are performed. An ISAR (Inverse Synthetic Aperture Radar) process for imaging a target is known (see Non-Patent Document 1). In the conventional ISAR process, when a rigid target is moving, the Doppler speed at each reflection point on the target varies depending on the rotational movement of the target. For this reason, the image generated by the ISAR process is usually different from the true target shape and is also different from the absolute value of the shape dimension.

したがって、ISAR画像により目標を識別する方法では、目標の回転や振動によるドップラ速度の変化により、生成画像の形状や寸法が真の目標とは大きく異なってしまい、識別が困難になる場合があった。   Therefore, in the method for identifying a target using an ISAR image, the shape and size of the generated image may be significantly different from the true target due to a change in Doppler speed due to the rotation or vibration of the target, and it may be difficult to identify the target. .

また、ISAR画像を基に、大きさを持つ目標内の特定の点を検出して追跡する場合には、クロスレンジが絶対値と異なるため、誤差が大きく追跡できなくなるという問題があった。   In addition, when a specific point within a target having a size is detected and tracked based on an ISAR image, there is a problem that the error is large and tracking cannot be performed because the cross range is different from the absolute value.

また、ISAR画像を生成しないドップラレーダ等において、大きさを持つ目標を測角する場合には、測角点が広がった複数点のベクトル合成(複素合成)になるため、グリントノイズが大きくなって測角の精度も低下してしまうという問題があった。   In addition, when measuring a target having a size in a Doppler radar or the like that does not generate an ISAR image, vector synthesis (complex synthesis) of a plurality of points in which the angle measurement points are widened, resulting in increased glint noise. There was a problem that the accuracy of angle measurement was also lowered.

また、目標が被搭載機のフライト軸方向にある場合には、剛体目標の回転運動が小さいと、クロスレンジ方向の左右上下等の目標の反射点のドップラが同じになるため、真の目標形状とSAR画像の形状が大きく異なってしまい、識別できないという問題があった。   In addition, when the target is in the flight axis direction of the mounted machine, if the rotational movement of the rigid target is small, the Doppler of the target reflection point in the cross range direction, such as left and right and up and down, will be the same, so the true target shape There is a problem that the shapes of the SAR image and the SAR image are greatly different and cannot be identified.

また、合成開口長が短い場合には、剛体目標の回転運動が小さくても、画像セルの絶対位置精度が低いという問題があった。   Further, when the synthetic aperture length is short, there is a problem that the absolute position accuracy of the image cell is low even if the rotational motion of the rigid target is small.

SAR方式(ISAR)、吉田、‘改定レーダ技術’、電子情報通信学会、pp.280-283(1996)SAR method (ISAR), Yoshida, 'Revised radar technology', IEICE, pp.280-283 (1996) SAR方式(レンジ圧縮)、大内、“リモートセンシングのための合成開口レーダの基礎”、東京電機大学出版局、pp.131-149(2003)SAR method (range compression), Ouchi, “Basics of Synthetic Aperture Radar for Remote Sensing”, Tokyo Denki University Press, pp.131-149 (2003) SAR方式(AZ圧縮)、大内、“リモートセンシングのための合成開口レーダの基礎”、東京電機大学出版局、pp.171-178(2003)SAR method (AZ compression), Ouchi, “Basics of Synthetic Aperture Radar for Remote Sensing”, Tokyo Denki University Press, pp.171-178 (2003) SAR処理方式(ポーラフォーマット変換再構成処理)、MEHRDAD SOUMEKH,“Synthetic Aperture Radar Signal Processing”, JOHN WILEY & SONS,INC.,pp.319-325(1999)SAR processing method (Polar format conversion reconstruction process), MEHRDAD SOUMEKH, “Synthetic Aperture Radar Signal Processing”, JOHN WILEY & SONS, INC., Pp.319-325 (1999) 位相モノパルス(位相比較モノパルス)方式、電子情報通信学会、改訂レーダ技術、pp.262-264(1996)Phase monopulse (phase comparison monopulse) method, IEICE, revised radar technology, pp.262-264 (1996) 振幅モノパルス(振幅比較モノパルス)方式、電子情報通信学会、改訂レーダ技術、pp.260-262(1996)Amplitude monopulse (amplitude comparison monopulse) system, IEICE, revised radar technology, pp.260-262 (1996) ファーフィ−ルドレンジ、電子情報通信学会、改訂レーダ技術、pp.144-145(1996)Far field range, IEICE, Revised radar technology, pp.144-145 (1996)

特開2008−304321号公報JP 2008-304321 A

以上述べたように、従来のレーダ装置に適用されるISAR処理では、レーダ観測位置と目標の複数反射点の各点との相対速度が異なる場合に、生成画像の形状や寸法が真の目標形状や寸法と大きく異なってしまうという課題があった。この課題に伴い、目標の識別能力や目標内の特定点の位置観測能力が低下する、ドップラレーダ等にあってはグリントノイズが大きくなって測角の精度も低下してしまう、目標の位置を観測可能な被搭載機のフライト軸方向に目標がある場合には、真の目標形状とSAR画像の形状が大きく異なって識別できない、合成開口長が短い場合には、剛体目標の回転運動が小さくても、画像セルの絶対位置精度が低いという課題があった。   As described above, in the ISAR processing applied to the conventional radar apparatus, the shape and size of the generated image are true target shapes when the relative velocity between the radar observation position and each of the target multiple reflection points is different. There was a problem that it was greatly different from the dimensions. Along with this problem, target identification ability and position observation ability of a specific point within the target are reduced, and in Doppler radar etc., the glint noise becomes large and the angle measurement accuracy is also reduced. When there is a target in the flight axis direction of the observable aircraft, the true target shape and the shape of the SAR image are very different and cannot be distinguished. When the synthetic aperture length is short, the rotational motion of the rigid target is small. However, there is a problem that the absolute position accuracy of the image cell is low.

本実施形態は上記課題に鑑みなされたもので、レーダと目標の複数反射点の各点の相対速度が異なる場合でも、真の目標形状や寸法に近い画像を生成することができるレーダ装置とそのレーダ信号処理方法を提供することを目的とする。さらには、識別能力向上、目標内の特定点の位置観測精度の向上、ドップラレーダ等におけるグリントノイズを低減して高精度に測角し目標の位置を観測することができる自機のフライト軸方向に目標がある場合でも目標形状に近いSAR画像を生成し、また、剛体目標の回転運動が小さく、合成開口長が小さい場合でも絶対位置精度の高いSAR画像を得ることのできるレーダ装置とレーダ信号処理方法を提供することを目的とする。   The present embodiment has been made in view of the above problems, and a radar apparatus capable of generating an image close to the true target shape and dimensions even when the relative speeds of the radar and the target multiple reflection points are different. An object of the present invention is to provide a radar signal processing method. Furthermore, the flight axis direction of the aircraft that can detect the target position by measuring the angle with high accuracy by reducing the glint noise in the Doppler radar etc. A radar apparatus and a radar signal capable of generating a SAR image close to the target shape even when the target is present and obtaining a SAR image with high absolute position accuracy even when the rotational motion of the rigid target is small and the synthetic aperture length is small An object is to provide a processing method.

上記の課題を解決するために、本実施形態は、実開口アンテナの受信信号から取得される目標をレンジ及びドップラ周波数の両軸で追跡し、前記目標の画像中心を得てレンジ圧縮及びAZ圧縮することにより、レンジ−ドップラ周波数軸で目標を画像化するISAR(Inverse Synthetic Aperture Radar:逆合成開口レーダ)方式のレーダ装置において、ビーム生成手段により、前記実開口アンテナの受信信号からモノパルス出力であるΣ、ΔAZ及びΔELそれぞれのビームを生成してサンプリング出力し、画像生成手段により、前記Σ、ΔAZ及びΔELそれぞれのビームのサンプリング出力からレンジ−ドップラ周波数軸におけるレンジ−ドップラ画像を生成し、目標画像抽出手段により前記Σビームのレンジ−ドップラ画像の中で所定の目標画像を抽出し、前記目標画像取得手段により、前記抽出された目標画像について、前記Σビームによるレンジと前記Σ、ΔAZ及びΔELそれぞれのビームによるモノパルス測角値とに基づいて前記目標画像のドップラ−レンジセルの位置を算出して目標画像を得る。   In order to solve the above-described problem, the present embodiment tracks a target acquired from a received signal of a real aperture antenna on both the range and Doppler frequency axes, obtains the image center of the target, and performs range compression and AZ compression. In the ISAR (Inverse Synthetic Aperture Radar) type radar apparatus that images the target on the range-Doppler frequency axis, the beam generation means outputs a monopulse from the received signal of the real aperture antenna. A beam of Σ, ΔAZ, and ΔEL is generated and sampled and output, and a range-Doppler image in the range-Doppler frequency axis is generated from the sampling output of the beams of Σ, ΔAZ, and ΔEL by the image generation means, and a target image is generated. The extraction means extracts a predetermined target image from the range-Doppler image of the Σ beam. The position of the Doppler-range cell of the target image is extracted based on the range of the Σ beam and the monopulse angle measurement value of the Σ, ΔAZ, and ΔEL beams for the extracted target image. To obtain a target image.

第1の実施形態に係るレーダ装置の構成を示すブロック図。1 is a block diagram showing a configuration of a radar apparatus according to a first embodiment. 図1に示すレーダ装置において、狭帯域処理器の構成を示すブロック図。The block diagram which shows the structure of a narrow-band processor in the radar apparatus shown in FIG. 図1に示すレーダ装置において、合成開口におけるISAR処理の概要を説明するための概念図。The conceptual diagram for demonstrating the outline | summary of the ISAR process in a synthetic aperture in the radar apparatus shown in FIG. 図1に示すレーダ装置において、合成開口による取得データからPRI毎のデータ列を抽出する様子を示す図。The figure which shows a mode that the data sequence for every PRI is extracted from the acquisition data by a synthetic aperture in the radar apparatus shown in FIG. 図1に示すレーダ装置において、参照信号として生成される近似曲線を示す図。The figure which shows the approximated curve produced | generated as a reference signal in the radar apparatus shown in FIG. 図1に示すレーダ装置において、SAR画像を取得する手順を示す図。The figure which shows the procedure which acquires the SAR image in the radar apparatus shown in FIG. 図1に示すレーダ装置において、ポーラフォーマット変換により格子点のデータを生成する様子を示す図。The figure which shows a mode that the data of a lattice point are produced | generated by polar format conversion in the radar apparatus shown in FIG. 図1に示すレーダ装置において、レンジ・測角値に基づく座標変換による画像変換の様子を示す図。The figure which shows the mode of the image conversion by the coordinate conversion based on a range and a measured angle value in the radar apparatus shown in FIG. 第2の実施形態に係るレーダ装置の構成を示すブロック図。The block diagram which shows the structure of the radar apparatus which concerns on 2nd Embodiment. 図9に示すレーダ装置において、合成開口におけるISAR処理の概要を説明するための概念図。The conceptual diagram for demonstrating the outline | summary of the ISAR process in a synthetic aperture in the radar apparatus shown in FIG. 図9に示すレーダ装置において、ビームの円錐の曲線特性を示す概念図。The conceptual diagram which shows the curve characteristic of the cone of a beam in the radar apparatus shown in FIG. 第3の実施形態に係るレーダ装置の構成を示すブロック図。The block diagram which shows the structure of the radar apparatus which concerns on 3rd Embodiment. 図12に示すレーダ装置におけるサンプル手法を説明するための図。The figure for demonstrating the sample method in the radar apparatus shown in FIG. 第4の実施形態に係るレーダ装置の構成を示すブロック図。The block diagram which shows the structure of the radar apparatus which concerns on 4th Embodiment. (a)はレンジ1〜Nのビーム幅の広がり具合を示し、(b)はレンジ1〜Nの誤差電圧と目標レンジとの関係を示す図。(A) shows the extent of the beam width of ranges 1 to N, and (b) is a diagram showing the relationship between the error voltage of ranges 1 to N and the target range.

以下、実施形態について、図面を参照して説明する。尚、各実施形態の説明において、同一部分には同一符号を付して示し、重複する説明を省略する。   Hereinafter, embodiments will be described with reference to the drawings. In the description of each embodiment, the same portions are denoted by the same reference numerals, and redundant description is omitted.

(第1の実施形態)
以下、図1乃至図8を参照して、第1の実施形態に係るレーダ装置について説明する。
(First embodiment)
The radar apparatus according to the first embodiment will be described below with reference to FIGS.

図1は第1の実施形態に係るレーダ装置の構成を示すブロック図である。図1に示すレーダ装置において、アンテナ1は複数のアンテナ素子を配列して大開口アレイを形成してなるフェーズドアレイアンテナであり、送受信器2の送受信部21から特定の周期で繰り返し供給される特定周波数の送信パルス信号(以下、PRF(Pulse Repetition Frequency)信号)を指定方向に送出してその反射波を受信する。送受信器2は、送受信部21において、アンテナ1の複数のアンテナ素子でそれぞれ受信された信号をビーム制御部22からの指示に従って位相制御を施し合成することで、任意の方向に受信ビームを形成してPRF受信信号を取得する。ここで、ビーム制御部22は指定された目標方向の測角値に基づいてΣビーム、ΔAZビーム、ΔELビームを形成するように、送受信部21に対して各ビームに対応する位相制御を施す。これにより、送受信部21はΣ信号、ΔAZ信号、ΔEL信号を生成して信号処理器3へ出力する。   FIG. 1 is a block diagram showing the configuration of the radar apparatus according to the first embodiment. In the radar apparatus shown in FIG. 1, the antenna 1 is a phased array antenna formed by arranging a plurality of antenna elements to form a large aperture array. The antenna 1 is repeatedly supplied from the transmitting / receiving unit 21 of the transmitter / receiver 2 at a specific cycle. A transmission pulse signal having a frequency (hereinafter referred to as a PRF (Pulse Repetition Frequency) signal) is transmitted in a specified direction and the reflected wave is received. The transceiver 2 forms a received beam in an arbitrary direction by performing phase control on the signals received by the plurality of antenna elements of the antenna 1 in the transceiver unit 21 according to instructions from the beam controller 22 and combining them. To obtain a PRF received signal. Here, the beam control unit 22 performs phase control corresponding to each beam on the transmission / reception unit 21 so as to form a Σ beam, a ΔAZ beam, and a ΔEL beam based on the angle measurement value in the designated target direction. As a result, the transmission / reception unit 21 generates a Σ signal, a ΔAZ signal, and a ΔEL signal and outputs them to the signal processor 3.

上記信号処理器3に入力されたΣ信号、ΔAZ信号、ΔEL信号は、それぞれの系統において、AD(Analog-Digital)変換部311,312,313で系統別にデジタル信号に変換されて狭帯域処理部32に送られる。狭帯域処理部32は、入力したΣ信号、ΔAZ信号、ΔEL信号から目標方向を測角演算する。ここで得られた測角値(Σ、ΔAZ、ΔEL)はビーム制御部22に送られる。   The Σ signal, ΔAZ signal, and ΔEL signal input to the signal processor 3 are converted into digital signals for each system by AD (Analog-Digital) conversion units 311, 312, and 313 in each system, and narrowband processing unit 32. The narrowband processing unit 32 calculates the angle of the target direction from the input Σ signal, ΔAZ signal, and ΔEL signal. The measured angle values (Σ, ΔAZ, ΔEL) obtained here are sent to the beam controller 22.

図2に上記狭帯域処理部32の具体的な構成を示す。図2において、デジタル化されたΣ信号、ΔAZ信号、ΔEL信号は、それぞれの系統において、データ抽出部3211〜3213に送られる。データ抽出部3211〜3213は、サンプリングレートを下げて処理規模を削減するために、入力されたΣ信号、ΔAZ信号、ΔEL信号を所定の周波数フィルタに通した後、PRI(Pulse repetition Interval)内のデータをレンジセル単位で取得する。各データ抽出部3211〜3213で得られたΣ、ΔAZ、ΔELそれぞれのレンジセル信号は、PRI軸FFT(Fast Fourier Transformation)処理部3221〜3223によってPRI軸の周波数領域信号に変換された後、レンジ(時間)圧縮部3231〜3233によってレンジ軸上でパルス圧縮されてレンジ圧縮信号となる。   FIG. 2 shows a specific configuration of the narrowband processing unit 32. In FIG. 2, the digitized Σ signal, ΔAZ signal, and ΔEL signal are sent to the data extraction units 3211 to 213 in each system. In order to reduce the sampling rate and reduce the processing scale, the data extraction units 3211 to 3213 pass the input Σ signal, ΔAZ signal, and ΔEL signal through a predetermined frequency filter, and then within a PRI (Pulse repetition Interval) Acquire data by range cell. The range cell signals of Σ, ΔAZ, and ΔEL obtained by the data extraction units 3211 to 213 are converted into PRI axis frequency domain signals by the PRI axis FFT (Fast Fourier Transformation) processing units 3221 to 223, and then the range ( Time) Pulse compression is performed on the range axis by the compression units 3231 to 233 to form a range compressed signal.

このうち、Σ系統のレンジ圧縮信号については、CFAR(Constant False Alarm Rate)処理部324に送られ、そのCFAR処理によって複数の極値を持つレンジセルが検出される。ここで検出されたレンジセルに対し、スレショルド検出部325は、閾値以上の極値を持つレンジセルを検出してそのレンジ圧縮信号を出力する。また、ΔAZ系統、ΔEL系統のレンジ圧縮信号については、それぞれセル検出部3261,3262に送られ、それぞれΣ系統のスレショルド検出部325で検出されたレンジセルに対応するレンジセルの信号が抽出されて、Σ系統の検出セルの信号と共に測角部327に送られ、ここで目標方向の測角演算が行われる。   Among them, the Σ-system range compressed signal is sent to a CFAR (Constant False Alarm Rate) processing unit 324, and a range cell having a plurality of extreme values is detected by the CFAR processing. For the range cell detected here, the threshold detection unit 325 detects a range cell having an extreme value equal to or greater than the threshold and outputs the range compression signal. Further, the range compression signals of the ΔAZ system and the ΔEL system are sent to the cell detection units 3261 and 3262, respectively, and a range cell signal corresponding to the range cell detected by the threshold detection unit 325 of the Σ system is extracted, and the Σ Together with the signal of the detection cell of the system, it is sent to the angle measuring unit 327, where angle measurement in the target direction is performed.

一方、上記AD変換部311でデジタル化されたΣ系統の受信信号は、データ抽出部33に送られる。このデータ抽出部33では、目標重心等の代表点のドップラ−レンジを観測してビームを目標方向に指向させ、ドップラ−レンジの情報を基にして、AZ圧縮用の際の参照信号を生成するためのデータを抽出する。ここで抽出されたデータは、PRI軸FFT処理部34にてCPI(Coherent Processing Interval)毎にPRI軸の周波数領域信号に変換され、レンジ圧縮部35にてレンジ軸上でパルス圧縮される。このパルス圧縮信号は、CFAR処理部36に送られ、そのCFAR処理によって複数の極値を持つレンジセルが検出される。ここで検出されたレンジセルに対し、スレショルド検出部37は、閾値以上の極値を持つレンジセルが検出する。レンジ・ドップラ算出部38は、スレショルド検出部37で検出されたレンジセルの信号からCPI毎に目標のレンジとドップラ周波数を算出する。参照信号生成部39は、CPI毎に検出したレンジとドップラ周波数からパルス圧縮のための参照信号を生成する。   On the other hand, the received signal of the Σ system digitized by the AD conversion unit 311 is sent to the data extraction unit 33. The data extraction unit 33 observes the Doppler range of a representative point such as the target center of gravity, directs the beam in the target direction, and generates a reference signal for AZ compression based on the Doppler range information. Extract data for The data extracted here is converted into a frequency domain signal of the PRI axis for each CPI (Coherent Processing Interval) by the PRI axis FFT processing unit 34, and is pulse-compressed on the range axis by the range compression unit 35. This pulse compression signal is sent to the CFAR processing unit 36, and a range cell having a plurality of extreme values is detected by the CFAR processing. For the range cell detected here, the threshold detection unit 37 detects a range cell having an extreme value equal to or greater than a threshold value. The range / Doppler calculation unit 38 calculates a target range and Doppler frequency for each CPI from the signal of the range cell detected by the threshold detection unit 37. The reference signal generation unit 39 generates a reference signal for pulse compression from the range and Doppler frequency detected for each CPI.

また、上記AD変換部311〜313でデジタル化されたΣ、ΔAZ、ΔELそれぞれの受信信号は、データ抽出部3A1,3A2,3A3にてサンプリング周波数を満足するための所定の周波数フィルタに通される。そして、レンジ圧縮部3B1,3B2,3B3によってレンジ圧縮され、さらにAZ圧縮部3C1,3C2,3C3によって参照信号に基づいてクロスレンジ圧縮(AZ圧縮)された後、ポーラフォーマット変換による画像化部3D1,3D2,3D3に送られる。画像化部3D1,3D2,3D3は、レンジ・クロスレンジ圧縮された受信信号に対して2次元FFT処理により周波数領域の信号に変換した後、ポーラフォーマット変換による画像処理を行い、2次元逆FFT処理により時間領域に信号に戻して画像化する。尚、画像化部3D1,3D2,3D3は、画像品質を向上する場合に用いればよく、省略可能である。   The received signals of Σ, ΔAZ, and ΔEL digitized by the AD conversion units 311 to 313 are passed through predetermined frequency filters for satisfying the sampling frequency in the data extraction units 3A1, 3A2, and 3A3. . Then, range compression is performed by the range compression units 3B1, 3B2, and 3B3, and further, cross-range compression (AZ compression) is performed based on the reference signal by the AZ compression units 3C1, 3C2, and 3C3, and then the imaging unit 3D1 by polar format conversion is performed. Sent to 3D2 and 3D3. The imaging units 3D1, 3D2, and 3D3 convert the received signals subjected to range / cross-range compression into frequency domain signals by two-dimensional FFT processing, perform image processing by polar format conversion, and perform two-dimensional inverse FFT processing. Thus, the signal is converted back to the time domain and imaged. Note that the imaging units 3D1, 3D2, and 3D3 may be used to improve image quality and can be omitted.

その後、Σ系統では、セル選定部3Eにおいて、画像化部3D1,3D2,3D3で生成された画像の中で振幅が所定のスレショルドを超えるセルを選定する。また、ΔAZ,ΔEL系統では、セル選定部3Eで選定されたセルをセル抽出部3F1,3F2に通知し、選定セルに対応するΔAZ,ΔELのセルを抽出する。測角部3Gは、セル選定部3E及びセル抽出部3F1,3F2で得られたΣ、ΔAZ、ΔELのセルの信号から目標方向を測角し、画像化部3Hは測角部3Gで測角された目標を画像化し出力する。   Thereafter, in the Σ system, the cell selection unit 3E selects a cell whose amplitude exceeds a predetermined threshold among the images generated by the imaging units 3D1, 3D2, and 3D3. In the ΔAZ, ΔEL system, the cells selected by the cell selection unit 3E are notified to the cell extraction units 3F1, 3F2, and cells of ΔAZ, ΔEL corresponding to the selected cells are extracted. The angle measuring unit 3G measures the target direction from the Σ, ΔAZ, and ΔEL cell signals obtained by the cell selecting unit 3E and the cell extracting units 3F1 and 3F2, and the imaging unit 3H measures the angle by the angle measuring unit 3G. The target is imaged and output.

上記構成によるレーダ装置において、以下、図3乃至図8を参照してそのISAR処理を説明する。   In the radar apparatus having the above configuration, the ISAR process will be described below with reference to FIGS.

まず、図3を参照してレーダ装置の合成開口の概要を説明する。飛翔体搭載のレーダ装置では、実開口ビームを目標に常に照射するように向けて、合成開口時間(1サイクル)内にPRI(Pulse Repetition Interval)間隔でパルスを送信し、その反射波を受信してパルス毎にPRI内のデータをレンジセル単位で取得する。この取得データを用いてISAR処理を実施して目標のISAR画像を得る。尚、図3は飛翔体搭載のレーダ装置の場合を示しているが、ISAR画像を得られれば、レーダ装置は固定の場合でもよい。   First, the outline of the synthetic aperture of the radar apparatus will be described with reference to FIG. In a radar device mounted on a flying object, a pulse is transmitted at a PRI (Pulse Repetition Interval) interval within a synthetic aperture time (1 cycle) and a reflected wave is received so that the target is always irradiated with a real aperture beam. For each pulse, data in the PRI is acquired in units of range cells. Using this acquired data, an ISAR process is performed to obtain a target ISAR image. Although FIG. 3 shows the case of a radar apparatus mounted on a flying object, the radar apparatus may be fixed if an ISAR image can be obtained.

図1において、ビーム制御部22により画像化したい範囲にビームを指向させて送受信し、その受信信号をΣ、ΔAZ、ΔELそれぞれの系統振り分けてデジタル信号に変換する。ビーム指向方向は、狭帯域処理部32により、目標方向を測角した方向とする。この測角については、狭帯域処理を用いなくても、広帯域のモノパルス測角の結果を用いてもよい。   In FIG. 1, the beam control unit 22 directs a beam to a range to be imaged and transmits / receives the received signal, and the received signal is divided into systems of Σ, ΔAZ, and ΔEL and converted into a digital signal. The beam directing direction is the direction in which the narrow band processing unit 32 measures the target direction. For this angle measurement, the result of broadband monopulse angle measurement may be used without using narrowband processing.

狭帯域処理部32では、図2に示すように、図1の系統のΣ、ΔAZ及びΔELのAD変換出力を入力とし、データ抽出部3211〜3213により、サンプリングレートを下げて処理規模を削減するために、周波数フィルタで所定の帯域を抽出した後、サンプリングする。そのうち、Σ系統の信号をPRI軸FFT処理部34でPRI軸の周波数領域信号に変換して、レンジ圧縮部35でレンジ方向にパルス圧縮し、CFAR処理部36で複数の極値を持つレンジセルを検出し、スレショルド検出部37で閾値以上の極値を持つレンジセルを検出する。この検出したセルに対応するΔAZとΔELのセルをセル抽出部3261,3262で抽出し、測角部327でΣ、ΔAZ及びΔEL信号を用いて測角する。この測角値を用いて、図1のビーム制御部22によりビームを目標方向に指向させる。   As shown in FIG. 2, the narrowband processing unit 32 receives AD conversion outputs of Σ, ΔAZ, and ΔEL of the system shown in FIG. 1, and the data extraction units 3211 to 213 reduce the sampling rate to reduce the processing scale. Therefore, sampling is performed after a predetermined band is extracted by the frequency filter. Among them, the Σ system signal is converted into a PRI-axis frequency domain signal by the PRI-axis FFT processing unit 34, pulse-compressed in the range direction by the range compression unit 35, and a range cell having a plurality of extreme values is converted by the CFAR processing unit 36. Then, the range detection unit 37 detects a range cell having an extreme value greater than or equal to the threshold value. The cells of ΔAZ and ΔEL corresponding to the detected cells are extracted by the cell extraction units 3261 and 3262, and the angle measuring unit 327 performs angle measurement using the Σ, ΔAZ, and ΔEL signals. Using this angle measurement value, the beam controller 22 in FIG. 1 directs the beam in the target direction.

次に、図1のΣ系のデータ抽出部33により、必要なサンプルデータを抽出する。データ抽出33では、目標重心等の代表点のドップラ−レンジを観測し、ビームを目標方向に指向させたときのドップラ−レンジ情報をもとにして、AZ圧縮用の参照信号を生成するためのデータを抽出する。すなわち、図4に示すように、合成開口のために取得したデータNallの中から、重複を含めてNヒットのPRIデータをPRI1〜PRINまで抽出する。ここで、CPI(m)(m=1〜M)とすると、CPI毎にPRI軸FFT処理部34でPRI軸の周波数領域信号に変換し、レンジ圧縮部35でレンジ方向にパルス圧縮し、CFAR処理部36で複数の極値を持つセルを検出し、スレショルド検出部37でスレショルド以上の極値を持つセルを検出し、レンジ・ドップラ算出部38により、目標のレンジ及びドップラ周波数を算出する。次に、参照信号生成部39により、CPI毎に算出したレンジとドップラ周波数により、図5に示す参照信号Refの近似曲線を算出する。これを定式化すると次の通りである。

Figure 0006363524
Next, necessary sample data is extracted by the Σ-system data extraction unit 33 of FIG. In the data extraction 33, a Doppler range of a representative point such as a target center of gravity is observed, and a reference signal for AZ compression is generated based on Doppler range information when the beam is directed in the target direction. Extract data. That is, as shown in FIG. 4, N hit PRI data including duplicates is extracted from PRI1 to PRIN from the data Nall acquired for the synthetic aperture. Here, assuming CPI (m) (m = 1 to M), the PRI axis FFT processing unit 34 converts each CPI into a PRI axis frequency domain signal, the range compression unit 35 performs pulse compression in the range direction, and CFAR. A cell having a plurality of extreme values is detected by the processing unit 36, a cell having an extreme value equal to or higher than the threshold is detected by the threshold detection unit 37, and a target range and Doppler frequency are calculated by the range / Doppler calculation unit 38. Next, the reference signal generation unit 39 calculates an approximate curve of the reference signal Ref shown in FIG. 5 based on the range and Doppler frequency calculated for each CPI. This is formulated as follows.
Figure 0006363524

一方、アンテナ1のモノパルス出力ΣとΔAZとΔELの出力は、それぞれAD変換部311,312,313によりデジタル信号に変換され、データ抽出部3A1,3A2,3A3によりサンプリング出力され、レンジ圧縮部3B1,3B2,3B3により変調信号の波形に対応するレンジ圧縮用の参照信号と相関演算された後、AZ圧縮部3C1,3C2,3C3において(1)式で生成された参照信号との相関演算によりAZ圧縮される。これを定式化すると次の通りである。   On the other hand, the monopulse outputs Σ, ΔAZ, and ΔEL of the antenna 1 are converted into digital signals by the AD conversion units 311, 312, and 313, sampled and output by the data extraction units 3 A 1, 3 A 2, and 3 A 3, respectively, and the range compression units 3 B 1, 3 B 1 After the correlation calculation with the reference signal for range compression corresponding to the waveform of the modulation signal is performed by 3B2 and 3B3, the AZ compression is performed by the correlation calculation with the reference signal generated by the expression (1) in the AZ compression units 3C1, 3C2, and 3C3 Is done. This is formulated as follows.

まずレンジ圧縮について述べる(非特許文献2参照)。レンジ圧縮は、入力信号とレンジ圧縮用信号の相関処理であり、これを周波数領域で行う場合について定式化すると次の通りである。

Figure 0006363524
First, range compression will be described (see Non-Patent Document 2). Range compression is a correlation process between an input signal and a range compression signal. The case of performing this in the frequency domain is formulated as follows.
Figure 0006363524

また、参照信号sref(線形チャープ信号)は次式で表現できる。

Figure 0006363524
The reference signal sref (linear chirp signal) can be expressed by the following equation.
Figure 0006363524

この参照信号sref(t)のサンプル長を入力信号に合わせて0埋めした信号に置き換える。

Figure 0006363524
The sample length of this reference signal sref (t) is replaced with a zero-padded signal in accordance with the input signal.
Figure 0006363524

時間軸上の信号にするには、このsを逆フーリエ変換すればよいが、この後クロスレンジ圧縮(AZ圧縮、非特許文献3参照)を行うために、信号sの(ω,u)軸のままとする。次にクロスレンジ圧縮を行うが、その参照信号は(1)式を用い、ω軸については同じ値を入れたものとする。

Figure 0006363524
In order to obtain a signal on the time axis, the s may be subjected to inverse Fourier transform. However, in order to perform cross range compression (AZ compression, see Non-Patent Document 3), the (ω, u) axis of the signal s Leave as it is. Next, cross range compression is performed, and the reference signal uses the equation (1), and the same value is input for the ω axis.
Figure 0006363524

(5)式と(6−3)式を乗算して信号csを得る。

Figure 0006363524
The signal cs is obtained by multiplying the equation (5) and the equation (6-3).
Figure 0006363524

これを用いて、u軸でFFT処理して周波数領域の信号fcs(ω,ku)を得る。

Figure 0006363524
Using this, FFT processing is performed on the u axis to obtain a frequency domain signal fcs (ω, ku).
Figure 0006363524

FFT画像出力は、信号fcsのω軸に関する逆FFTにより算出することができる。

Figure 0006363524
The FFT image output can be calculated by inverse FFT with respect to the ω axis of the signal fcs.
Figure 0006363524

次に、より精度の高い画像生成手法として、図6及び図7を参照してポーラフォーマット変換による画像化部3D1,3D2,3D3のSAR画像処理について述べる(非特許文献4参照)。   Next, as a more accurate image generation method, SAR image processing of the imaging units 3D1, 3D2, and 3D3 by polar format conversion will be described with reference to FIGS. 6 and 7 (see Non-Patent Document 4).

図6は第1の実施形態に係るレーダ装置でSAR画像を取得する手順を示している。まず、サンプリング時間tに対する合成開口長のサンプル点のデータuを入力信号sig(t,u)として取り込んで(S11)、t軸に対するフーリエ変換(FFTx(ω,u))を行い(S12)、その結果Sin(ω,u)に参照信号のt軸に対するフーリエ変換(FFTx(ω,u))の結果を乗算してレンジ圧縮を行い(S13)、レンジ圧縮の出力sを得る。このレンジ圧縮信号sを用いて、u軸に関してフーリエ変換する(S14)。

Figure 0006363524
FIG. 6 shows a procedure for acquiring a SAR image by the radar apparatus according to the first embodiment. First, sample point data u of the synthetic aperture length with respect to the sampling time t is taken in as an input signal sig (t, u) (S11), and Fourier transform (FFTx (ω, u)) with respect to the t-axis is performed (S12). As a result, the range compression is performed by multiplying Sin (ω, u) by the result of Fourier transform (FFTx (ω, u)) with respect to the t-axis of the reference signal (S13), and an output s of the range compression is obtained. Using this range compression signal s, Fourier transform is performed with respect to the u-axis (S14).
Figure 0006363524

次にAZ圧縮用の参照信号を生成する。

Figure 0006363524
Next, a reference signal for AZ compression is generated.
Figure 0006363524

また、kx,kyは次式により求めることができる。

Figure 0006363524
Further, kx and ky can be obtained by the following equations.
Figure 0006363524

画像中心については、レンジ、ドップラ、角度を観測して、次式により算出する位置とする。

Figure 0006363524
For the image center, the range, Doppler, and angle are observed and set as the position calculated by the following equation.
Figure 0006363524

ここで、観測距離、観測AZ角度、観測EL角度については、例えば、図4の各CPIにおける観測値の平均値とすればよい。 Here, the observation distance, the observation AZ angle, and the observation EL angle may be, for example, the average value of the observation values at each CPI in FIG.

このAZ圧縮参照信号を用いて、次式によりfsmを算出する(S15)。

Figure 0006363524
Using this AZ compression reference signal, fsm is calculated by the following equation (S15).
Figure 0006363524

fsmを用いて、図7に示すポーラフォーマット変換を行い、kx軸、ky軸で格子点のデータF(kx,ky)を生成する(S16)。ポーラフォーマット変換は、取得データを用いて、(kx,ky)軸の格子点の補間後データを内挿手法等を用いて算出する手法であり、細部については非特許文献4の通りであり割愛する。このfsmを用いて、2次元逆FFTによりポーラフォーマット変換を用いた画像fsarを出力する(S17)。

Figure 0006363524
Polar format conversion shown in FIG. 7 is performed using fsm, and lattice point data F (kx, ky) is generated on the kx axis and ky axis (S16). Polar format conversion is a method of calculating the interpolated data of the (kx, ky) axis grid points using the acquired data using an interpolation method, etc., and details are omitted as in Non-Patent Document 4. To do. Using this fsm, an image fsar using polar format conversion is output by two-dimensional inverse FFT (S17).
Figure 0006363524

このポーラフォーマット変換による画像処理は、常に必要ではなく、画像品質を向上する場合に用いればよい。 Image processing by this polar format conversion is not always necessary, and may be used to improve image quality.

次に、Σ系の画像出力を用いて、X−Y軸(ポーラフォーマット変換を用いない場合はドップラ−レンジ軸に対応)の画像の中で、振幅が所定のスレショルドを超えるセル(P個)を選定する。このP個のセルの各々に対応するΔAZ及びΔEL系の画像出力のセルをセル抽出部3F1,3F2で抽出し、測角部3Gにおいて、Σ、ΔAZ及びΔELのセル信号(複素信号)を用いてモノパルス測角する(非特許文献5参照)。この測角値と距離を用いて、画像化部3Hにおいて、次式により画像化する。

Figure 0006363524
Next, using the Σ-system image output, cells whose amplitude exceeds a predetermined threshold in an XY axis image (corresponding to a Doppler range axis when polar format conversion is not used) (P cells) Is selected. The ΔAZ and ΔEL image output cells corresponding to each of the P cells are extracted by the cell extraction units 3F1 and 3F2, and the angle measurement unit 3G uses Σ, ΔAZ, and ΔEL cell signals (complex signals). Then, monopulse angle measurement is performed (see Non-Patent Document 5). Using this angle measurement value and distance, the imaging unit 3H forms an image according to the following equation.
Figure 0006363524

この画像変換の様子を図8に示す。図8において、(a)はレンジ−ドップラ画像であり、(b)は(a)の画像のレンジ・測角値による座標変換を行った変換画像である。図8から明らかなように、レンジ・測角値による座標変換により、真値の位置に近い画像を生成することができる。 The state of this image conversion is shown in FIG. In FIG. 8, (a) is a range-Doppler image, and (b) is a converted image obtained by performing coordinate conversion using the range / angle-measured value of the image of (a). As is apparent from FIG. 8, an image close to the true value position can be generated by coordinate conversion based on the range / angle measurement value.

尚、第1の実施形態の手法によれば、被搭載機のフライト方向に目標があり、目標が回転運動を持たず、クロスレンジ方向にアンビギュイティを持ち、SAR画像が実際の目標形状と大きく異なる場合でも、実開口の測角はクロスレンジ方向の区別ができるため、実際の形状に近い画像を得ることができる。   According to the method of the first embodiment, there is a target in the flight direction of the mounted machine, the target has no rotational motion, has ambiguity in the cross range direction, and the SAR image has the actual target shape. Even if they differ greatly, the angle of measurement of the actual aperture can be distinguished in the cross-range direction, so that an image close to the actual shape can be obtained.

(第2の実施形態)
第1の実施形態では、実開口の測角精度が低い場合がある。第2の実施形態は実開口の測角精度を改善する手法について述べる。
(Second Embodiment)
In the first embodiment, the angle measurement accuracy of the actual aperture may be low. The second embodiment describes a method for improving the angle measurement accuracy of the actual aperture.

図9は第2の実施形態に係るレーダ装置の構成を示すブロック図である。図9において、図1に示した構成と異なる点は、ΔSARウェイト乗算部3I、AZ圧縮部3C4、ポーラフォーマット変換による画像化部3D4、セル抽出部3F3を備え、図1に示した測角部3Gに代わって、SAR測角部3J、実開口測角部3K、測角相関部3Lを備えるようにしたことにある。   FIG. 9 is a block diagram showing a configuration of a radar apparatus according to the second embodiment. 9 is different from the configuration shown in FIG. 1 in that a ΔSAR weight multiplication unit 3I, an AZ compression unit 3C4, an imaging unit 3D4 by polar format conversion, and a cell extraction unit 3F3 are provided, and the angle measuring unit shown in FIG. Instead of 3G, a SAR angle measurement unit 3J, an actual aperture angle measurement unit 3K, and an angle measurement correlation unit 3L are provided.

第1の実施形態では、実開口アンテナのモノパルス測角値を用いたが、測角精度を向上するために、図10に示すような合成開口を用いる。すなわち、目標が飛翔経路軸付近にある場合は、合成開口によるリニアアレイのエンドファイア方向(リニアアレイの軸方向)にビームが指向される。この場合、ブロ−ドサイド方向(リニアアレイの軸に垂直方向)にビームを形成する場合に比べて合成開口ビーム幅は広がるが、合成開口長が大きい場合には、実開口アンテナによるビーム幅よりも十分狭くなる。このエンドファイアアレイのPRI毎のNサンプルのデータを使った和信号であるΣ信号と、1〜N/2とN/2+1〜Nの開口2分割した差信号であるΔ信号を用いてモノパルス測角(非特許文献5参照)を行う。これによってSAR測角値が得られる(特許文献1参照)。   In the first embodiment, the monopulse angle measurement value of the real aperture antenna is used. However, in order to improve the angle measurement accuracy, a synthetic aperture as shown in FIG. 10 is used. That is, when the target is in the vicinity of the flight path axis, the beam is directed in the endfire direction of the linear array (the axial direction of the linear array) by the synthetic aperture. In this case, the synthetic aperture beam width is wider than when the beam is formed in the broad side direction (perpendicular to the axis of the linear array), but when the synthetic aperture length is large, it is larger than the beam width of the real aperture antenna. It becomes narrow enough. Monopulse measurement using the Σ signal, which is a sum signal using data of N samples for each PRI of this endfire array, and the Δ signal, which is a difference signal obtained by dividing the apertures 1 to N / 2 and N / 2 + 1 to N by two. A corner (see Non-Patent Document 5) is performed. Thereby, the SAR angle measurement value is obtained (see Patent Document 1).

このビ−ムは、エンドファイアアレイの軸を中心に図10及び図11に示すような円錐状のビ−ムであり、この円錐状に目標がある場合には、一意に位置を特定できない。この対策のために、第1の実施形態と同じように実開口アンテナのモノパルス測角値(AZ角とEL角)を用いる。具体的には、目標距離RとエンドファイアによるSAR測角値θSARとにより、図10に示すビームの円錐の曲線を描き、(16)式で決まる実開口アンテナによる目標位置から最も近い円錐状の点を抽出し、第1の実施形態の変換位置と置き換えればよい。   This beam is a conical beam as shown in FIGS. 10 and 11 centering on the axis of the endfire array. If the target is in this conical shape, the position cannot be uniquely identified. For this countermeasure, the monopulse angle measurement value (AZ angle and EL angle) of the real aperture antenna is used as in the first embodiment. Specifically, the conical curve of the beam shown in FIG. 10 is drawn by the target distance R and the SAR angle measurement value θSAR by the endfire, and the conical shape closest to the target position by the real aperture antenna determined by the equation (16) is drawn. A point may be extracted and replaced with the conversion position of the first embodiment.

以上の処理を図9に示す系統図で説明する。第1の実施形態と同様の部分の説明は省略する。レンジ圧縮部3B1でΣ信号のレンジ圧縮した出力のうち、上記のように開口2分割した信号に対してΔ信号用のサイドローブを低減するためのウェイトをΔSARウェイト乗算部3Iにより乗算し、AZ圧縮部3C4、ポーラフォーマット変換による画像化部3D1でそれぞれAZ圧縮、ポーラフォーマット変換を行って画像化する。さらに、セル選定部3Eにて選定されたセルをセル抽出部3F3に通知し、Σの選定セルに対応するΔSARの信号を抽出する。これにより得られたΣ及びΔSARの信号を用いて、次式の誤差電圧εを演算する。

Figure 0006363524
The above processing will be described with reference to the system diagram shown in FIG. A description of the same parts as those in the first embodiment will be omitted. Of the outputs of the Σ signal that have been subjected to range compression by the range compression unit 3B1, a signal for reducing the side lobe for the Δ signal is multiplied by the ΔSAR weight multiplication unit 3I with respect to the signal divided into two apertures as described above, and AZ The compression unit 3C4 and the polar format conversion imaging unit 3D1 perform AZ compression and polar format conversion, respectively, to form an image. Further, the cell selection unit 3E notifies the cell extraction unit 3F3 of the cell selected, and extracts the ΔSAR signal corresponding to the selected cell of Σ. Using the Σ and ΔSAR signals thus obtained, an error voltage ε of the following equation is calculated.
Figure 0006363524

測角部3Gにおいて、上記誤差電圧εと予め作成した誤差電圧テーブルとの比較により測角すれば、円錐状の角度を算出することができる。このSAR測角値を用いて、前述の処理により、高精度な目標位置を算出すればよい。   If the angle measuring unit 3G measures the angle by comparing the error voltage ε with an error voltage table created in advance, a conical angle can be calculated. Using this SAR angle measurement value, a highly accurate target position may be calculated by the above-described processing.

なお、SAR測角について、ΣとΔsar信号による位相モノパルス測角の場合について述べたが、Δsar信号の代わりに、Σ信号ビームからスクイント角分ずらせたΣ2sarビームを用いて、スクイント測角してもよいのは言うまでもない。この場合は、誤差電圧εは、次式となる。

Figure 0006363524
As for SAR angle measurement, the case of phase monopulse angle measurement using Σ and Δsar signals has been described. However, instead of the Δsar signal, Σ2sar beam obtained by dividing the Σ signal beam by the squint angle may be used for squint angle measurement. Needless to say, it is good. In this case, the error voltage ε is as follows.
Figure 0006363524

この誤差電圧εと予め作成した誤差電圧テーブルを用いて測角することができる。また、振幅のみを用いた振幅モノパルスの場合でも測角できる(非特許文献6参照)。

Figure 0006363524
The angle can be measured using this error voltage ε and a previously created error voltage table. Further, the angle can be measured even in the case of an amplitude monopulse using only the amplitude (see Non-Patent Document 6).
Figure 0006363524

この測角値と距離を用いて、画像化部3Hにおいて、(16)式により3次元位置を画像化する。これにより、第2の実施形態によれば、目標が飛翔経路軸付近にある場合でも、実開口の測角精度を改善して真値に近い画像を取得することができる。   Using the angle measurement value and the distance, the imaging unit 3H images the three-dimensional position according to equation (16). Thereby, according to 2nd Embodiment, even when a target exists in the vicinity of a flight path axis | shaft, the angle measurement precision of a real opening can be improved and the image close | similar to a true value can be acquired.

なお、第2の実施形態においては、合成開口長を開口2分割したモノパルス測角手法の主旨で説明しているが、合成開口長内のサンプル列のFFTを用いる手法であり、ドップラ周波数軸におけるモノパルス手法とも言える。これより、剛体目標の回転運動が大きい場合には、そのドップラ周波数成分による位相誤差が生じるため、クロスレンジ方向に画像が歪むのと同様に、測角値も誤差を持つことになる。したがって、本手法により絶対位置精度の向上を図る条件としては剛体目標の回転運動が低い場合か、剛体目標の回転運動を観測し、そのドップラによる位相成分を補正できる場合等となり、モノパルス測角手法を用いる分だけ、クロスレンジの精度向上となる。   In the second embodiment, the gist of the monopulse angle measurement method in which the synthetic aperture length is divided into two apertures has been described. However, this is a method that uses the FFT of the sample sequence within the synthetic aperture length, and is based on the Doppler frequency axis. It can also be said to be a monopulse technique. As a result, when the rotational motion of the rigid target is large, a phase error due to the Doppler frequency component occurs, so that the angle measurement value also has an error in the same manner as the image is distorted in the cross range direction. Therefore, the conditions for improving the absolute position accuracy with this method are when the rotational motion of the rigid target is low, or when the rotational motion of the rigid target is observed and the phase component due to the Doppler can be corrected. As much as is used, the accuracy of the cross range is improved.

(第3の実施形態)
第1の実施形態及び第2の実施形態では、Σの他にΔAZとΔEL系統が必要であり、処理規模が増える。第3の実施形態は、Σ、ΔAZ、ΔELの各系統の共用を図り、処理規模の増加を抑えるものである。
(Third embodiment)
In the first embodiment and the second embodiment, ΔAZ and ΔEL systems are required in addition to Σ, and the processing scale increases. In the third embodiment, each of the Σ, ΔAZ, and ΔEL systems is shared, and an increase in processing scale is suppressed.

図12は第3の実施形態に係るレーダ装置の構成を示すブロック図、図13は第3の実施形態におけるサンプル手法を説明するための図である。第1及び第2の実施形態と同様の部分の説明は省略する。   FIG. 12 is a block diagram showing a configuration of a radar apparatus according to the third embodiment, and FIG. 13 is a diagram for explaining a sampling method in the third embodiment. Description of the same parts as those in the first and second embodiments is omitted.

第1及び第2の実施形態では、Σ、ΔAZ及びΔELの3系統分を同時にサンプルしていたが、第3の実施形態では、送受信器2において、送受信部21で得られるΣ、ΔAZ、ΔEL信号をサンプル切り替え部23で各信号をサンプリングし、AD変換を1系統分としてAD変換部31のみとしている。   In the first and second embodiments, three systems of Σ, ΔAZ, and ΔEL are simultaneously sampled. However, in the third embodiment, Σ, ΔAZ, and ΔEL obtained by the transmitter / receiver 21 in the transceiver 2. Each signal is sampled by the sample switching unit 23, and only one AD conversion unit 31 is used for AD conversion.

すなわち、送受信器2のサンプル切り替え部23は、図13(a),(b)に示すように、各系統を時分割に切り替える。この場合、3系統分を順にまとめて同時にサンプルする方式とすると、時間差が生じるため、角度誤差が大きくなる。この対策のために、3系統分のサンプル順を擬似ランダム化して、間引きサンプルを実施する。AD変換部を1chとするためには、送受信器2において、サンプル切り替え部23により、擬似ランダム順に送受信部21からの3系統の順番を切り替える。以降の処理は第1及び第2の実施形態と同様である。   That is, the sample switching unit 23 of the transceiver 2 switches each system to time division as shown in FIGS. 13 (a) and 13 (b). In this case, if the system for sampling three systems in order is sampled at the same time, a time difference is generated, and the angle error becomes large. For this countermeasure, the sampling order for the three systems is pseudo-randomized and thinned samples are implemented. In order to set the AD conversion unit to 1ch, in the transceiver 2, the sample switching unit 23 switches the order of the three systems from the transmission / reception unit 21 in a pseudo-random order. Subsequent processing is the same as in the first and second embodiments.

(第4の実施形態)
レーダと目標の距離は遠距離の場合と近距離の場合があるが、第1及び第2の実施形態では、特に近距離の場合に、ファーフィールドレンジ(非特許文献7参照)の条件を満足せず、ビーム幅が広がる場合が想定される。この場合に、遠距離の場合と同じ測角テーブルを用いると、誤差電圧との対応がとれずに、測角値に誤差が生じる。第4の実施形態ではこの対策を施している。図14は、第4の実施形態に係るレーダ装置の構成を示すブロック図である。
(Fourth embodiment)
The distance between the radar and the target may be a long distance or a short distance. In the first and second embodiments, the far field range (see Non-Patent Document 7) is satisfied, particularly in the short distance. Without, it is assumed that the beam width is widened. In this case, if the same angle measurement table as in the case of a long distance is used, the error voltage is not matched and an error occurs in the angle measurement value. In the fourth embodiment, this countermeasure is taken. FIG. 14 is a block diagram illustrating a configuration of a radar apparatus according to the fourth embodiment.

測角曲線としては、一般にファーフィールド条件を満足する遠距離(レンジ1)はビーム幅が狭く、ファーフィールド条件を満足しない近距離(レンジNに近い方)はビーム幅が広くなる。この様子を図15に示す。図15(a)はレンジ1〜Nのビーム幅の広がり具合を示し、図15(b)はレンジ1〜Nの誤差電圧と目標レンジとの関係を示す図である。図15から明らかなように、観測レンジに対応したΣとΔビームより、測角用の誤差電圧を予めテーブル化しておき、観測した目標レンジに応じて測角テーブルを切り替えればよい。   As the angle measurement curve, generally, a far distance (range 1) that satisfies the far field condition has a narrow beam width, and a short distance that does not satisfy the far field condition (closer to the range N) has a wide beam width. This is shown in FIG. FIG. 15A shows the extent of the beam width in the ranges 1 to N, and FIG. 15B shows the relationship between the error voltage in the ranges 1 to N and the target range. As is apparent from FIG. 15, the angle measurement error voltage may be tabulated in advance from the Σ and Δ beams corresponding to the observation range, and the angle measurement table may be switched according to the observed target range.

図14に示すブロック図を参照して第4の実施形態を説明する。測角テーブルとしては、Σ系統のセル選定部3Eとセル抽出部3F3でそれぞれ選定・抽出されたセルに基づきSAR測角値を観測するSAR可変測角部3Mと、ΔAZ、ΔEL系統のセル抽出結果に基づいて実開口の可変測角を行う実開口可変測角部3Nがあり、各測角部3M,3NはΣ系統のレンジドップラ算出部38で算出したレンジによりテーブルを切り替える。   The fourth embodiment will be described with reference to the block diagram shown in FIG. As the angle measurement table, the SAR variable angle measurement unit 3M that observes the SAR angle measurement value based on the cells selected and extracted by the cell selection unit 3E and the cell extraction unit 3F3 of the Σ system, and the cell extraction of the ΔAZ and ΔEL systems, respectively. There is an actual aperture variable angle measuring unit 3N that performs variable angle measurement of the actual aperture based on the result, and each angle measuring unit 3M, 3N switches the table according to the range calculated by the range Doppler calculating unit 38 of the Σ system.

尚、第1乃至第4の実施形態において、処理規模を低減するために狭帯域処理部32により、ビ−ム方向を制御する手法について述べたが、広帯域処理(SAR処理)や、他のセンサからの情報等によりビーム方向がわかれば、狭帯域処理を行わなくてもよい。   In the first to fourth embodiments, the method of controlling the beam direction by the narrowband processing unit 32 in order to reduce the processing scale has been described. However, wideband processing (SAR processing) and other sensors are described. If the beam direction is known from the information from, etc., narrowband processing may not be performed.

その他、本実施形態は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。更に、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。   In addition, the present embodiment is not limited to the above-described embodiment as it is, and can be embodied by modifying the components without departing from the scope of the invention in the implementation stage. In addition, various inventions can be formed by appropriately combining a plurality of components disclosed in the embodiment. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, constituent elements over different embodiments may be appropriately combined.

1…アンテナ、11…送信アンテナ素子、2…送受信器、21…送受信器、22…ビーム制御器、23…サンプル切り替え部、3…信号処理器、31,311,312,313…AD変換部、32…狭帯域処理部、3211〜3213…データ抽出部、3221〜3223…PRI軸FFT処理部、3231〜3233…レンジ(時間)圧縮部、324…CFAR処理部、325…スレショルド検出部、3261,3262…セル検出部、327…測角部、33,3A1,3A2,3A3…データ抽出部、34…PRI軸FFT処理部、35,3B1,3B2,3B3…レンジ圧縮部、36…CFAR処理部、37…スレショルド検出部、38…レンジ・ドップラ算出部、39…参照信号生成部、3C1,3C2,3C3,3C4…AZ圧縮部、3D1,3D2,3D3,3D4…ポーラフォーマット変換による画像化部、3E…セル選定部、3F1,3F2,3F3…セル抽出部、3G…測角部、3H…画像化部、3I…ΔSARウェイト算出部、3J…SAR測角部、3K…実開口測角部、3L…測角相関部、3M…SAR可変測角部、3N…実開口可変測角部。 DESCRIPTION OF SYMBOLS 1 ... Antenna, 11 ... Transmitting antenna element, 2 ... Transmitter / receiver, 21 ... Transmitter / receiver, 22 ... Beam controller, 23 ... Sample switching part, 3 ... Signal processor, 31, 311, 312, 313 ... AD converter, 32 ... Narrow band processing unit, 3211 to 2213 ... Data extraction unit, 3221 to 2223 ... PRI axis FFT processing unit, 3231 to 2233 ... Range (time) compression unit, 324 ... CFAR processing unit, 325 ... Threshold detection unit, 3261 3262 ... Cell detection unit, 327 ... Angle measurement unit, 33, 3A1, 3A2, 3A3 ... Data extraction unit, 34 ... PRI axis FFT processing unit, 35, 3B1, 3B2, 3B3 ... Range compression unit, 36 ... CFAR processing unit, 37: Threshold detection unit, 38: Range / Doppler calculation unit, 39: Reference signal generation unit, 3C1, 3C2, 3C3, 3C4 ... AZ compression 3D1, 3D2, 3D3, 3D4 ... Polar format conversion imaging unit, 3E ... cell selection unit, 3F1,3F2,3F3 ... cell extraction unit, 3G ... angle measuring unit, 3H ... imaging unit, 3I ... ΔSAR weight calculation 3J: SAR angle measuring unit, 3K: actual aperture angle measuring unit, 3L: angle measurement correlating unit, 3M: SAR variable angle measuring unit, 3N: actual aperture variable angle measuring unit.

Claims (8)

実開口アンテナの受信信号から取得される目標をレンジ及びドップラ周波数の両軸で追跡し、前記目標の画像中心を得てレンジ圧縮及びAZ圧縮することにより、レンジ−ドップラ周波数軸で目標を画像化するISAR(Inverse Synthetic Aperture Radar:逆合成開口レーダ)方式のレーダ装置において、
前記実開口アンテナの受信信号からモノパルス出力であるΣ、ΔAZ及びΔELそれぞれのビームを生成してサンプリング出力するビーム生成手段と、
前記Σ、ΔAZ及びΔELそれぞれのビームのサンプリング出力からレンジ−ドップラ周波数軸におけるレンジ−ドップラ画像を生成する画像生成手段と、
前記Σビームのレンジ−ドップラ画像の中で所定の目標画像を抽出する目標画像抽出手段と、
前記抽出された目標画像について、前記Σビームによるレンジと前記Σ、ΔAZ及びΔELそれぞれのビームによるモノパルス測角値とに基づいて前記目標画像のドップラ−レンジセルの位置を算出して目標画像を得る目標画像取得手段と
を具備するレーダ装置。
The target acquired from the received signal of the real aperture antenna is tracked on both the range and Doppler frequency axes, and the target is imaged on the range-Doppler frequency axis by obtaining the target image center and performing range compression and AZ compression. In the ISAR (Inverse Synthetic Aperture Radar) type radar device,
Beam generating means for generating and sampling each beam of Σ, ΔAZ and ΔEL, which are monopulse outputs, from the received signal of the real aperture antenna;
Image generation means for generating a range-Doppler image on the range-Doppler frequency axis from the sampling outputs of the beams of Σ, ΔAZ, and ΔEL;
A target image extracting means for extracting a predetermined target image from the range-Doppler image of the Σ beam;
A target for obtaining a target image by calculating a position of a Doppler-range cell of the target image based on a range by the Σ beam and a monopulse angle value by the beams of Σ, ΔAZ, and ΔEL for the extracted target image. A radar apparatus comprising image acquisition means.
前記目標画像取得手段は、前記実開口アンテナの合成開口長におけるCPI(Coherent Processing Interval:コヒーレント処理間隔)の任意の複数ヒットを2分割し、前記Σビームによる和信号ΣSARによりドップラ−レンジセルを選定し、選定したセルに対応するΣSARと前記Σビームによる差信号ΔSARを用いたモノパルス処理により円錐角θSARを測角し、測角した円錐角θSARの中で最も前記ΔAZ及びΔELのビームそれぞれの測角値θAZ、θELに近い角度θAZcalとθELcalにより、前記ドップラ−レンジセルの位置を算出する請求項1記載のレーダ装置。   The target image acquisition means divides an arbitrary plurality of hits of CPI (Coherent Processing Interval) in the synthetic aperture length of the real aperture antenna into two, and selects a Doppler-range cell by the sum signal ΣSAR by the Σ beam The cone angle θSAR is measured by monopulse processing using the ΣSAR corresponding to the selected cell and the difference signal ΔSAR by the Σ beam, and the angle of each of the above-mentioned ΔAZ and ΔEL beams among the measured cone angles θSAR is measured. The radar apparatus according to claim 1, wherein the position of the Doppler-range cell is calculated from angles θAZcal and θELcal close to values θAZ and θEL. 前記ビーム生成手段は、前記実開口アンテナの合成開口長におけるCPI(Coherent Processing Interval:コヒーレント処理間隔)の任意の複数ヒット内で、疑似ランダム順に前記Σビーム、ΔAZビーム及びΔELビームを切り替えてサンプリング出力する請求項1または2に記載のレーダ装置。   The beam generating means performs sampling output by switching the Σ beam, ΔAZ beam, and ΔEL beam in pseudo-random order within an arbitrary plurality of hits of CPI (Coherent Processing Interval) in the synthetic aperture length of the real aperture antenna. The radar device according to claim 1 or 2. 前記目標画像取得手段は、前記実開口アンテナの合成開口長におけるCPI(Coherent Processing Interval:コヒーレント処理間隔)の任意の複数ヒットを2分割し、前記Σ、ΔAZ、ΔELの各ビームの出力から第1のモノパルス測角曲線を算出し、前記Σビームによる和信号ΣSARと前記Σビームによる差信号ΔSARから第2のモノパルス測角曲線を算出し、所定のレンジ範囲毎にテーブル化して、観測した目標レンジにより前記第1、第2のモノパルス測角曲線を切り替えて、切り替えられたモノパルス測角曲線に基づいて前記目標画像のドップラ−レンジセルの位置を算出する請求項1記載のレーダ装置。   The target image acquisition means divides an arbitrary plurality of hits of CPI (Coherent Processing Interval) in the synthetic aperture length of the real aperture antenna into two, and outputs the first from the outputs of the Σ, ΔAZ, and ΔEL beams. The monopulse angle curve is calculated, the second monopulse angle curve is calculated from the sum signal ΣSAR by the Σ beam and the difference signal ΔSAR by the Σ beam, and is tabulated for each predetermined range range, and the observed target range The radar apparatus according to claim 1, wherein the first and second monopulse angle measurement curves are switched by calculating a Doppler-range cell position of the target image based on the switched monopulse angle measurement curves. 実開口アンテナの受信信号から取得される目標をレンジ及びドップラ周波数の両軸で追跡し、前記目標の画像中心を得てレンジ圧縮及びAZ圧縮することにより、レンジ−ドップラ周波数軸で目標を画像化するISAR(Inverse Synthetic Aperture Radar:逆合成開口レーダ)方式のレーダ信号処理方法において、
前記実開口アンテナの受信信号からモノパルス出力であるΣ、ΔAZ及びΔELそれぞれのビームを生成してサンプリング出力し、
前記Σ、ΔAZ及びΔELそれぞれのビームのサンプリング出力からレンジ−ドップラ周波数軸におけるレンジ−ドップラ画像を生成し、
前記Σビームのレンジ−ドップラ画像の中で所定の目標画像を抽出し、
前記抽出された目標画像について、前記Σビームによるレンジと前記Σ、ΔAZ及びΔELそれぞれのビームによるモノパルス測角値とに基づいて前記目標画像のドップラ−レンジセルの位置を算出して目標画像を得るレーダ信号処理方法。
The target acquired from the received signal of the real aperture antenna is tracked on both the range and Doppler frequency axes, and the target is imaged on the range-Doppler frequency axis by obtaining the target image center and performing range compression and AZ compression. In the ISAR (Inverse Synthetic Aperture Radar) type radar signal processing method,
Generate and output Σ, ΔAZ and ΔEL beams, which are monopulse outputs, from the received signal of the real aperture antenna,
A range-Doppler image in the range-Doppler frequency axis is generated from the sampling outputs of the beams of Σ, ΔAZ, and ΔEL,
A predetermined target image is extracted from the range-Doppler image of the Σ beam,
A radar that obtains a target image by calculating the position of a Doppler-range cell of the target image based on the extracted target image based on the range of the Σ beam and the monopulse angle value of the Σ, ΔAZ, and ΔEL beams. Signal processing method.
前記目標画像の取得は、前記実開口アンテナの合成開口長におけるCPI(Coherent Processing Interval:コヒーレント処理間隔)の任意の複数ヒットを2分割し、前記Σビームによる和信号ΣSARによりドップラ−レンジセルを選定し、選定したセルに対応するΣSARと前記Σビームによる差信号ΔSARを用いたモノパルス処理により円錐角θSARを測角し、測角した円錐角θSARの中で最も前記ΔAZ及びΔELのビームそれぞれの測角値θAZ、θELに近い角度θAZcalとθELcalにより、前記ドップラ−レンジセルの位置を算出する請求項5記載のレーダ信号処理方法。   The target image is acquired by dividing a plurality of hits of CPI (Coherent Processing Interval) in the synthetic aperture length of the real aperture antenna into two, and selecting a Doppler-range cell by the sum signal ΣSAR using the Σ beam. The cone angle θSAR is measured by monopulse processing using the ΣSAR corresponding to the selected cell and the difference signal ΔSAR by the Σ beam, and the angle of each of the above-mentioned ΔAZ and ΔEL beams among the measured cone angles θSAR is measured. The radar signal processing method according to claim 5, wherein the position of the Doppler-range cell is calculated based on angles θAZcal and θELcal close to the values θAZ and θEL. 前記ビーム生成は、前記実開口アンテナの合成開口長におけるCPI(Coherent Processing Interval:コヒーレント処理間隔)の任意の複数ヒット内で、疑似ランダム順に前記Σビーム、ΔAZビーム及びΔELビームを切り替えてサンプリング出力する請求項5に記載のレーダ信号処理方法。   The beam generation is performed by sampling and outputting the Σ beam, ΔAZ beam, and ΔEL beam in pseudo-random order within an arbitrary plurality of hits of the CPI (Coherent Processing Interval) in the synthetic aperture length of the real aperture antenna. The radar signal processing method according to claim 5. 前記目標画像の取得は、前記実開口アンテナの合成開口長におけるCPI(Coherent Processing Interval:コヒーレント処理間隔)の任意の複数ヒットを2分割し、前記Σ、ΔAZ、ΔELの各ビームの出力から第1のモノパルス測角曲線を算出し、前記Σビームによる和信号ΣSARと前記Σビームによる差信号ΔSARから第2のモノパルス測角曲線を算出し、所定のレンジ範囲毎にテーブル化して、観測した目標レンジにより前記第1、第2のモノパルス測角曲線を切り替えて、切り替えられたモノパルス測角曲線に基づいて前記目標画像のドップラ−レンジセルの位置を算出する請求項5記載のレーダ信号処理方法。   The acquisition of the target image is performed by dividing an arbitrary plurality of hits of CPI (Coherent Processing Interval) in the synthetic aperture length of the real aperture antenna into two and firstly outputting from the outputs of the beams of Σ, ΔAZ, and ΔEL. The monopulse angle curve is calculated, the second monopulse angle curve is calculated from the sum signal ΣSAR by the Σ beam and the difference signal ΔSAR by the Σ beam, and is tabulated for each predetermined range range, and the observed target range The radar signal processing method according to claim 5, wherein the first and second monopulse angle measurement curves are switched by calculating the position of the Doppler-range cell of the target image based on the switched monopulse angle measurement curves.
JP2015017275A 2015-01-30 2015-01-30 Radar apparatus and radar signal processing method Active JP6363524B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015017275A JP6363524B2 (en) 2015-01-30 2015-01-30 Radar apparatus and radar signal processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015017275A JP6363524B2 (en) 2015-01-30 2015-01-30 Radar apparatus and radar signal processing method

Publications (2)

Publication Number Publication Date
JP2016142578A JP2016142578A (en) 2016-08-08
JP6363524B2 true JP6363524B2 (en) 2018-07-25

Family

ID=56570211

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015017275A Active JP6363524B2 (en) 2015-01-30 2015-01-30 Radar apparatus and radar signal processing method

Country Status (1)

Country Link
JP (1) JP6363524B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6989663B1 (en) * 2020-07-16 2022-01-05 株式会社東芝 Radar device
CN115079151A (en) * 2022-05-25 2022-09-20 北京集创北方系统技术有限公司 Detection system and detection method based on Doppler radar

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4855747A (en) * 1987-08-17 1989-08-08 Trustees Of The University Of Pennsylvania Method of target imaging and identification
JP2938728B2 (en) * 1993-09-27 1999-08-25 防衛庁技術研究本部長 Radar equipment
JPH1078481A (en) * 1996-09-02 1998-03-24 Mitsubishi Electric Corp Aircraft loaded radar
JP2008304321A (en) * 2007-06-07 2008-12-18 Toshiba Corp Tracking device
JP4976968B2 (en) * 2007-09-20 2012-07-18 株式会社東芝 Radar equipment
JP5728418B2 (en) * 2012-03-19 2015-06-03 株式会社東芝 Target angle detection device, target angle detection method, and guidance device

Also Published As

Publication number Publication date
JP2016142578A (en) 2016-08-08

Similar Documents

Publication Publication Date Title
JP6271032B2 (en) Antenna specification estimating device and radar device
CN108885254B (en) Object detection device
JP6352688B2 (en) Radar apparatus and radar signal processing method thereof
JP6301749B2 (en) Doppler radar apparatus and radar signal processing method thereof
JP6470152B2 (en) Radar apparatus and radar signal processing method
JP2014182010A (en) Radar apparatus
JP2018205174A (en) Radar device and radar signal processing method thereof
JP6346082B2 (en) Pulse compression radar apparatus and radar signal processing method thereof
JP6466263B2 (en) Radar apparatus and radar signal processing method
CN108776342A (en) A kind of high speed platform SAR moving-target detection and speed estimation method at a slow speed
JP6367143B2 (en) Synthetic aperture radar apparatus and radar signal processing method thereof
JP6324327B2 (en) Passive radar equipment
JP5044358B2 (en) Radar equipment
JP2019168255A (en) Pulse compression radar device and radar signal processing method therefor
JP6546109B2 (en) Radar equipment
Li et al. Radar and inverse scattering
JP2017106799A (en) Synthetic-aperture radar device and radar signal processing method thereof
JP5247077B2 (en) Moving target detection device
JP6251087B2 (en) Target detection apparatus and target detection method
JP6178244B2 (en) Synthetic aperture radar apparatus and image processing method thereof
JP6178222B2 (en) Synthetic aperture radar apparatus and image processing method thereof
JP6363524B2 (en) Radar apparatus and radar signal processing method
JP6367134B2 (en) Radar apparatus and radar signal processing method
JP6400340B2 (en) Pulse compression radar apparatus and radar signal processing method thereof
JP2015052549A (en) Synthetic aperture radar system and method of image processing thereof

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170815

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20170911

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20170911

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180523

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20180529

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180628

R150 Certificate of patent or registration of utility model

Ref document number: 6363524

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150