JP2009109457A - Radar device - Google Patents

Radar device Download PDF

Info

Publication number
JP2009109457A
JP2009109457A JP2007284860A JP2007284860A JP2009109457A JP 2009109457 A JP2009109457 A JP 2009109457A JP 2007284860 A JP2007284860 A JP 2007284860A JP 2007284860 A JP2007284860 A JP 2007284860A JP 2009109457 A JP2009109457 A JP 2009109457A
Authority
JP
Japan
Prior art keywords
frequency
doppler frequency
unit
pulse signal
weather clutter
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
JP2007284860A
Other languages
Japanese (ja)
Other versions
JP5196959B2 (en
Inventor
Yuichiro Fukuma
雄一郎 福間
Toshio Wakayama
俊夫 若山
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 JP2007284860A priority Critical patent/JP5196959B2/en
Publication of JP2009109457A publication Critical patent/JP2009109457A/en
Application granted granted Critical
Publication of JP5196959B2 publication Critical patent/JP5196959B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide a radar device that accurately measures Doppler frequency of a weather clutter and appropriately suppresses the weather clutter. <P>SOLUTION: The device for detecting a target by transmission and reception of a pulse signal includes a weather clutter detecting section 6 for detecting the weather clutter from a low frequency reception pulse signal to output the Doppler frequency and strength of the weather clutter, a first return corrective section 8 for executing return correction processing to the Doppler frequency for each beam scan to output the Doppler frequency after the first correction, a second return corrective section 10 for executing return correction processing to the Doppler frequency after the first correction to output the Doppler frequency after the second correction, and an averaging processing section 12 for executing averaging processing to the Doppler frequency after the second correction to output the Doppler frequency average value. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

この発明は、例えば艦船等に搭載され、電磁波を送受信して目標物までの距離等を測定するレーダ装置に関する。   The present invention relates to a radar device that is mounted on, for example, a ship and that measures the distance to a target by transmitting and receiving electromagnetic waves.

従来から、電磁波の送受信によって目標物までの距離を測定する装置として、レーダ装置が知られている。この種のレーダ装置では、パルス信号を空間に送信し、目標物で反射したパルス信号を受信して、受信したパルス信号の空間伝搬時間を測定することにより、目標物までの距離を測定している。
受信したパルス信号には、目標物で反射したパルス信号の他に、例えば、雨粒や雲等で反射した信号(気象クラッタ)が含まれている。目標物で反射したパルス信号がこの気象クラッタに埋もれると、目標物までの距離を測定することが困難になる。
Conventionally, a radar device is known as a device that measures the distance to a target by transmitting and receiving electromagnetic waves. This type of radar device measures the distance to the target by transmitting a pulse signal to space, receiving the pulse signal reflected by the target, and measuring the spatial propagation time of the received pulse signal. Yes.
The received pulse signal includes, for example, a signal (weather clutter) reflected by raindrops, clouds, etc., in addition to the pulse signal reflected by the target. If the pulse signal reflected by the target is buried in the weather clutter, it becomes difficult to measure the distance to the target.

そこで、気象クラッタを抑圧するために、フィルタが利用されている。このフィルタの中心周波数は、気象クラッタのドップラー周波数によって決まる値であり、フィルタの特性を決める重要なパラメータの1つである。すなわち、フィルタの中心周波数を正確に決定するためには、気象クラッタのドップラー周波数を高精度に測定する必要がある。   Therefore, filters are used to suppress weather clutter. The center frequency of this filter is a value determined by the Doppler frequency of the weather clutter, and is one of important parameters that determine the characteristics of the filter. That is, in order to accurately determine the center frequency of the filter, it is necessary to measure the Doppler frequency of the weather clutter with high accuracy.

気象クラッタのドップラー周波数を測定する方法としては、気象レーダで降雨量などの測定に広く利用されているパルスペア方式を利用する方法がある。パルスペア方式では、同一距離から反射した2つの連続するパルス信号を受信し、パルス信号の位相差を検出してドップラー周波数を測定している(例えば、非特許文献1参照)。   As a method for measuring the Doppler frequency of the weather clutter, there is a method using a pulse pair method widely used for measuring rainfall and the like by a weather radar. In the pulse pair method, two continuous pulse signals reflected from the same distance are received, the phase difference between the pulse signals is detected, and the Doppler frequency is measured (for example, see Non-Patent Document 1).

また、ドップラー周波数を高精度に測定するために、レーダ覆域中の距離方向および方位方向について、気象クラッタのドップラー周波数の測定値を平均化する処理が実行されている。
ドップラー周波数の測定値を平均化するものとして、従来のレーダ信号処理装置は、クラッタのドップラー周波数を距離方向に測定した後、ドップラー周波数の測定値に対してレンジブロック単位で平均化処理を実行する平均処理部を備えている(例えば、特許文献1参照)。
Further, in order to measure the Doppler frequency with high accuracy, a process of averaging measured values of the Doppler frequency of the weather clutter is performed for the distance direction and the azimuth direction in the radar coverage.
In order to average the measured values of the Doppler frequency, the conventional radar signal processing apparatus measures the Doppler frequency of the clutter in the distance direction, and then performs an averaging process on the measured value of the Doppler frequency in units of range blocks. An average processing unit is provided (see, for example, Patent Document 1).

しかしながら、パルスペア方式を用いて気象クラッタのドップラー周波数を測定する場合には、パルス繰り返し周波数に対応してドップラー周波数の測定値の符号が反転し、距離方向または方位方向で値が不連続となる折り返し現象が発生する恐れがある。折り返し現象が発生すると、ドップラー周波数の測定値を平均化した場合に、符号が反転した箇所では、正確なドップラー周波数の平均値を得ることができず、ドップラー周波数を高精度に測定することができないという問題点があった。   However, when measuring the Doppler frequency of meteorological clutter using the pulse pair method, the sign of the Doppler frequency measurement value is inverted corresponding to the pulse repetition frequency, and the value becomes discontinuous in the distance direction or azimuth direction. The phenomenon may occur. When the aliasing phenomenon occurs, when the measured value of Doppler frequency is averaged, the accurate average value of Doppler frequency cannot be obtained at the place where the sign is reversed, and the Doppler frequency cannot be measured with high accuracy. There was a problem.

そこで、ドップラー周波数を高精度に測定するために、ドップラー周波数の測定値に対する平均化処理に先立って、ドップラー周波数の測定値に対する折り返し補正処理が実行されている。折り返し補正処理では、距離方向または方位方向で値が連続となるように、符号が反転したドップラー周波数の測定値に対して補正値を加算している(例えば、非特許文献1参照)。   Therefore, in order to measure the Doppler frequency with high accuracy, the aliasing correction processing for the measured value of Doppler frequency is performed prior to the averaging processing for the measured value of Doppler frequency. In the aliasing correction process, the correction value is added to the measured value of the Doppler frequency whose sign is inverted so that the value is continuous in the distance direction or the azimuth direction (see, for example, Non-Patent Document 1).

なお、上記従来のレーダ装置では、気象クラッタのドップラー周波数の測定値に対する平均化処理および折り返し補正処理は、レーダ装置のビーム走査毎(レーダ装置のビームがレーダ覆域を1周する毎)に実行されている。   In the above-described conventional radar apparatus, the averaging process and the aliasing correction process for the measured value of the Doppler frequency of the weather clutter are executed for each beam scanning of the radar apparatus (every time the radar apparatus beam makes one round of the radar coverage). Has been.

気象研究所編集、「気象研究所技術報告(第19号)ドップラーレーダによる気象・海象の研究」、気象研究所、1986年3月31日、p.10‐13、p.63‐64Edited by Meteorological Research Institute, "Meteorological Research Institute Technical Report (No. 19) Meteorological / Sea Research by Doppler Radar", Meteorological Research Institute, March 31, 1986, p. 10-13, p. 63-64 特開2001−201568号公報JP 2001-151568 A

例えば、気象レーダ等のレーダ装置では、低速でアンテナを回転させながら各方位についてパルス信号を送受信しているので、1回のビーム走査で大量のパルス信号を送受信することができる。そのため、気象クラッタのドップラー周波数の測定値に対する平均化処理および折り返し補正処理を実行して、ドップラー周波数を高精度に測定することができる。   For example, in a radar device such as a weather radar, a pulse signal is transmitted / received in each direction while rotating an antenna at a low speed, so that a large amount of pulse signals can be transmitted / received by one beam scanning. Therefore, the Doppler frequency can be measured with high accuracy by executing the averaging process and the aliasing correction process on the measured value of the Doppler frequency of the weather clutter.

これに対して、例えば艦船等に搭載される捜索用のレーダ装置では、高速でアンテナを回転させながら各方位についてパルス信号を送受信している。そのため、ビーム走査1回あたりの時間が短くなり、送受信されるパルス信号の数(パルスヒット数)が気象レーダの場合よりも減少する。したがって、気象クラッタのドップラー周波数の測定値に対する平均化処理および折り返し補正処理を実行した場合であっても、ドップラー周波数を高精度に測定することができないという問題点があった。   On the other hand, for example, in a search radar device mounted on a ship or the like, pulse signals are transmitted and received in each direction while rotating an antenna at high speed. For this reason, the time per one beam scan is shortened, and the number of pulse signals transmitted and received (number of pulse hits) is smaller than in the case of weather radar. Therefore, there is a problem that the Doppler frequency cannot be measured with high accuracy even when the averaging process and the aliasing correction process are performed on the measured value of the Doppler frequency of the weather clutter.

そこで、艦船等に搭載されるレーダ装置において気象クラッタのドップラー周波数を高精度に測定するための方法として、複数回のビーム走査によって測定されたドップラー周波数どうしを平均化すること、すなわちビーム走査間(複数回のビーム走査の間)でドップラー周波数の測定値を平均化することが考えられる。しかしながら、ビーム走査間でドップラー周波数の測定値を平均化する場合には、ビーム走査毎に実行される折り返し補正処理による折り返し方向が、ビーム走査間で必ずしも同じ方向にならない。そのため、ビーム走査間で正確なドップラー周波数の平均値を得ることができず、ドップラー周波数を高精度に測定することができないという問題点があった。   Therefore, as a method for measuring the Doppler frequency of meteorological clutter with high accuracy in a radar device mounted on a ship or the like, averaging Doppler frequencies measured by a plurality of beam scans, that is, between beam scans ( It is conceivable to average the Doppler frequency measurements between multiple beam scans). However, when the measured Doppler frequency values are averaged between the beam scans, the return direction by the return correction process executed for each beam scan is not necessarily the same direction between the beam scans. Therefore, there is a problem that an accurate average value of Doppler frequency cannot be obtained between beam scans, and the Doppler frequency cannot be measured with high accuracy.

この発明は、上記のような課題を解決するためになされたものであって、その目的は、気象クラッタのドップラー周波数を高精度に測定し、気象クラッタを適切に抑圧することができるレーダ装置を提供することにある。   The present invention has been made to solve the above-described problems, and an object of the present invention is to provide a radar device that can measure the Doppler frequency of the weather clutter with high accuracy and appropriately suppress the weather clutter. It is to provide.

この発明に係るレーダ装置は、周波数変換を行うための局部発振信号を発生して出力する局部発振部と、低周波送信パルス信号を発生するとともに、局部発振信号を用いて低周波送信パルス信号を高周波送信パルス信号に周波数変換して出力する送信部と、高周波送信パルス信号を空間に放射するとともに、目標物で反射した高周波送信パルス信号を高周波受信パルス信号として受信する空中線部と、局部発振信号を用いて高周波受信パルス信号を低周波受信パルス信号に周波数変換するとともに、低周波受信パルス信号をAD変換して出力する受信部と、デジタル化された低周波受信パルス信号から不要な信号である気象クラッタを検出し、気象クラッタのドップラー周波数および強度を出力する気象クラッタ検出部と、ビームがレーダ覆域を走査するビーム走査毎に、ドップラー周波数に対する折り返し補正処理を実行し、第1補正後ドップラー周波数を出力する第1折り返し補正部と、複数回のビーム走査の間であるビーム走査間で、第1補正後ドップラー周波数に対する折り返し補正処理を実行し、第2補正後ドップラー周波数を出力する第2折り返し補正部と、第2補正後ドップラー周波数に対する平均化処理を実行して、ドップラー周波数平均値を出力する平均処理部と、第1補正後ドップラー周波数およびドップラー周波数平均値のいずれか一方に基づいて、デジタル化された低周波受信パルス信号に含まれる気象クラッタを抑圧するためのフィルタ係数を算出するフィルタ設定部と、フィルタ係数に基づいて、デジタル化された低周波受信パルス信号に対するフィルタ処理を実行し、気象クラッタを抑圧した抑圧後受信パルス信号を出力するクラッタ抑圧処理部と、抑圧後受信パルス信号に基づいて目標物を検出する目標検出部とを備えたものである。   A radar apparatus according to the present invention generates a local oscillation signal for performing frequency conversion and outputs a local oscillation unit, a low-frequency transmission pulse signal, and a low-frequency transmission pulse signal using the local oscillation signal. A transmitter that converts the frequency into a high-frequency transmission pulse signal and outputs it, an antenna that radiates the high-frequency transmission pulse signal to space and receives the high-frequency transmission pulse signal reflected by the target as a high-frequency reception pulse signal, and a local oscillation signal This is an unnecessary signal from the receiving unit that converts the frequency of the high-frequency received pulse signal into the low-frequency received pulse signal using the A / D, and AD-converts and outputs the low-frequency received pulse signal, and the digitized low-frequency received pulse signal A weather clutter detector that detects weather clutter and outputs the Doppler frequency and intensity of the weather clutter, and a beam that covers the radar The first aliasing correction unit that executes the aliasing correction process for the Doppler frequency and outputs the first corrected Doppler frequency for each beam scanning, and the beam scanning that is between a plurality of beam scannings. Executes aliasing correction processing for the corrected Doppler frequency, outputs a second aliasing correction unit that outputs the second corrected Doppler frequency, executes averaging processing for the second corrected Doppler frequency, and outputs the Doppler frequency average value. Filter setting for calculating a filter coefficient for suppressing weather clutter included in the digitized low frequency received pulse signal based on either the average processing unit and the first corrected Doppler frequency or the Doppler frequency average value And a filter for the digitized low frequency received pulse signal based on the filter coefficients. Run the data processing, but with the clutter suppression processing section for outputting a received pulse signal after suppression that suppresses the weather clutter, and a target detector for detecting a target object on the basis of suppression after receiving the pulse signal.

この発明のレーダ装置によれば、第1折り返し補正部は、ビーム走査毎に気象クラッタのドップラー周波数に対する折り返し補正処理を実行し、第1補正後ドップラー周波数を出力する。第2折り返し補正部は、ビーム走査間で第1補正後ドップラー周波数に対する折り返し補正処理を実行し、第2補正後ドップラー周波数を出力する。また、平均処理部は、第2補正後ドップラー周波数に対する平均化処理を実行して、ドップラー周波数平均値を出力する。クラッタ抑圧処理部は、ドップラー周波数平均値に基づいて算出されたフィルタ係数に応じて、デジタル化された低周波受信パルス信号に対するフィルタ処理を実行し、気象クラッタを抑圧する。
そのため、気象クラッタのドップラー周波数を高精度に測定し、気象クラッタを適切に抑圧することができる。
According to the radar apparatus of the present invention, the first aliasing correction unit executes the aliasing correction process for the Doppler frequency of the weather clutter for each beam scan, and outputs the first corrected Doppler frequency. The second aliasing correction unit executes aliasing correction processing for the first corrected Doppler frequency between beam scans, and outputs the second corrected Doppler frequency. The average processing unit executes an averaging process on the second corrected Doppler frequency, and outputs an average Doppler frequency value. The clutter suppression processing unit executes filter processing on the digitized low-frequency received pulse signal according to the filter coefficient calculated based on the average Doppler frequency value, and suppresses weather clutter.
For this reason, the Doppler frequency of the weather clutter can be measured with high accuracy, and the weather clutter can be appropriately suppressed.

以下、この発明の各実施の形態について図に基づいて説明するが、各図において同一、または相当する部材、部位については、同一符号を付して説明する。   Hereinafter, embodiments of the present invention will be described with reference to the drawings. In the drawings, the same or corresponding members and parts will be described with the same reference numerals.

実施の形態1.
図1は、この発明の実施の形態1に係るレーダ装置を示すブロック構成図である。
図1において、このレーダ装置は、局部発振部1と、送信部2と、サーキュレータ3と、空中線部4と、受信部5と、気象クラッタ検出部6と、気象クラッタマップ作成部7と、第1折り返し補正部8と、切り換えスイッチ9と、第2折り返し補正部10と、データ記憶部11と、平均処理部12と、フィルタ設定部13と、クラッタ抑圧処理部14と、目標検出部15とを備えている。
Embodiment 1 FIG.
1 is a block diagram showing a radar apparatus according to Embodiment 1 of the present invention.
In FIG. 1, this radar apparatus includes a local oscillator 1, a transmitter 2, a circulator 3, an antenna unit 4, a receiver 5, a weather clutter detector 6, a weather clutter map generator 7, 1 folding correction unit 8, changeover switch 9, second folding correction unit 10, data storage unit 11, average processing unit 12, filter setting unit 13, clutter suppression processing unit 14, target detection unit 15, It has.

局部発振部1は、周波数変換を行うための局部発振信号を発生して、送信部2および受信部5に出力する。
送信部2は、低周波の正弦波信号をパルス変調して低周波送信パルス信号を発生する。また、送信部2は、局部発振部1からの局部発振信号を用いて、低周波送信パルス信号をギガヘルツ帯の高周波送信パルス信号に周波数変換して出力する。
The local oscillation unit 1 generates a local oscillation signal for performing frequency conversion and outputs it to the transmission unit 2 and the reception unit 5.
The transmitter 2 performs pulse modulation on the low frequency sine wave signal to generate a low frequency transmission pulse signal. Further, the transmission unit 2 uses the local oscillation signal from the local oscillation unit 1 to frequency-convert the low frequency transmission pulse signal into a high frequency transmission pulse signal in the gigahertz band and outputs the signal.

サーキュレータ3は、送信部2からの高周波送信パルス信号を空中線部4に出力するとともに、空中線部4からの高周波受信パルス信号を受信部5に出力する。
空中線部4は、ビームを電子的に走査しながら、サーキュレータ3からの高周波送信パルス信号を空間に放射する。また、空中線部4は、空間中に存在する目標物で反射した高周波送信パルス信号を高周波受信パルス信号として受信する。
The circulator 3 outputs the high frequency transmission pulse signal from the transmission unit 2 to the antenna unit 4 and also outputs the high frequency reception pulse signal from the antenna unit 4 to the reception unit 5.
The antenna unit 4 emits a high-frequency transmission pulse signal from the circulator 3 to the space while electronically scanning the beam. The antenna unit 4 receives a high-frequency transmission pulse signal reflected by a target existing in the space as a high-frequency reception pulse signal.

受信部5は、局部発振部1からの局部発振信号を用いて、空中線部4からの高周波受信パルス信号を低周波帯(数十〜数百メガヘルツ)の低周波受信パルス信号に周波数変換する。また、受信部5は、低周波受信パルス信号をAD変換(アナログ−デジタル変換)し、デジタル化された低周波受信パルス信号を出力する。   The receiving unit 5 uses the local oscillation signal from the local oscillation unit 1 to frequency-convert the high frequency reception pulse signal from the antenna unit 4 into a low frequency reception pulse signal in a low frequency band (several tens to several hundreds of megahertz). The receiving unit 5 performs AD conversion (analog-digital conversion) on the low-frequency reception pulse signal and outputs a digitized low-frequency reception pulse signal.

以下、図2を参照しながら、受信部5による低周波受信パルス信号のAD変換処理について説明する。
図2は、この発明の実施の形態1に係る低周波受信パルス信号のAD変換処理を示す説明図である。
図2において、低周波送信パルス信号は、パルス繰り返し時間間隔で発生され、高周波送信パルス信号に周波数変換されて送信される。また、目標物で反射して受信された高周波受信パルス信号は、低周波受信パルス信号に周波数変換される。
Hereinafter, the AD conversion processing of the low-frequency received pulse signal by the receiving unit 5 will be described with reference to FIG.
FIG. 2 is an explanatory diagram showing AD conversion processing of the low-frequency received pulse signal according to Embodiment 1 of the present invention.
In FIG. 2, a low frequency transmission pulse signal is generated at a pulse repetition time interval, and is converted into a high frequency transmission pulse signal and transmitted. Further, the high frequency received pulse signal reflected and received by the target is frequency-converted into a low frequency received pulse signal.

ここで、パルス繰り返し時間は、10の区間に分割されており、1区間がサンプリング間隔に相当する。また、これら10区間のうち、低周波送信パルス信号の発生タイミングを除く9区間は、レーダ装置からの距離(第1距離〜第9距離)に対応している。
すなわち、低周波受信パルス信号をAD変換することにより、レーダからの距離毎の情報が得られる。また、AD変換した低周波受信パルス信号をパルス繰り返し時間毎に並べ替えることにより、同一距離におけるパルス繰り返し時間毎の時系列データが、各距離について作成される。
図2には、例として、デジタル化された低周波受信パルス信号がパルス繰り返し時間毎に並べ替えられて作成された第5距離における時系列データを示す。この時系列データのサンプリング時間間隔は、パルス繰り返し時間と一致する。
Here, the pulse repetition time is divided into 10 sections, and one section corresponds to the sampling interval. Of these 10 sections, 9 sections excluding the generation timing of the low-frequency transmission pulse signal correspond to the distance from the radar device (the first distance to the ninth distance).
That is, information for each distance from the radar can be obtained by AD converting the low-frequency received pulse signal. Further, by rearranging the AD-converted low-frequency received pulse signals for each pulse repetition time, time series data for each pulse repetition time at the same distance is created for each distance.
FIG. 2 shows, as an example, time-series data at the fifth distance created by rearranging digitized low-frequency received pulse signals for each pulse repetition time. The sampling time interval of this time series data coincides with the pulse repetition time.

気象クラッタ検出部6は、受信部5でデジタル化されてパルス繰り返し時間毎に並べ替えられた低周波送信パルス信号から、雨粒や雲等で反射した不要な信号である気象クラッタを検出し、気象クラッタのドップラー周波数および強度を出力する。
気象クラッタのドップラー周波数は、前述したパルスペア方式を用いて、次式(1)のように算出される。また、気象クラッタのドップラー周波数は、レーダ装置からの各距離について算出される。
The weather clutter detection unit 6 detects weather clutter, which is an unnecessary signal reflected by raindrops, clouds, and the like, from the low-frequency transmission pulse signals digitized by the reception unit 5 and rearranged for each pulse repetition time. Outputs the Doppler frequency and intensity of the clutter.
The Doppler frequency of the weather clutter is calculated as the following equation (1) using the pulse pair method described above. Further, the Doppler frequency of the weather clutter is calculated for each distance from the radar device.

Figure 2009109457
Figure 2009109457

式(1)において、f(x)は第x距離における気象クラッタのドップラー周波数、Tはパルス繰り返し時間、argは複素数値の位相角、Mはパルスヒット数、v(x,k)は第x距離におけるkヒット目の低周波受信パルス信号、*は複素共役を示している。
また、気象クラッタの強度は、各距離について低周波受信パルス信号の電力または電圧を測定することによって求められる。
なお、気象クラッタのドップラー周波数および強度は、ビームが走査範囲(ビーム覆域)を走査するビーム走査毎に、各距離および各方位について測定される。
In equation (1), f d (x) is the Doppler frequency of the weather clutter at the x-th distance, T s is the pulse repetition time, arg is the complex-valued phase angle, M is the number of pulse hits, and v (x, k) is The low-frequency received pulse signal of the kth hit at the xth distance, * indicates the complex conjugate.
In addition, the intensity of the weather clutter is obtained by measuring the power or voltage of the low-frequency received pulse signal for each distance.
In addition, the Doppler frequency and intensity | strength of a weather clutter are measured about each distance and each azimuth | direction for every beam scanning in which a beam scans a scanning range (beam coverage).

気象クラッタマップ作成部7は、気象クラッタ検出部6から出力された気象クラッタのドップラー周波数および強度のデータに基づいて、ドップラー周波数および強度の距離特性を示す気象クラッタマップを、方位毎に作成する。
以下、図3を参照しながら、気象クラッタマップ作成部7で作成される気象クラッタマップについて説明する。
The weather clutter map creation unit 7 creates a weather clutter map indicating the distance characteristics of the Doppler frequency and intensity for each direction based on the Doppler frequency and intensity data of the weather clutter output from the weather clutter detection unit 6.
The weather clutter map created by the weather clutter map creation unit 7 will be described below with reference to FIG.

図3は、この発明の実施の形態1に係る気象クラッタマップを示す説明図である。図3では、例として、i番目の方位(第i方位)における気象クラッタマップを示す。
図3において、気象クラッタマップは、距離方向を横軸とし、ドップラー周波数および強度をそれぞれ縦軸として方位毎に作成されている。また、距離方向の各レンジセルは、図2で示した第1距離〜第9距離にそれぞれ対応している。
FIG. 3 is an explanatory diagram showing a weather clutter map according to Embodiment 1 of the present invention. FIG. 3 shows a weather clutter map in the i-th azimuth (i-th azimuth) as an example.
In FIG. 3, the weather clutter map is created for each azimuth with the distance direction as the horizontal axis and the Doppler frequency and intensity as the vertical axis. Each range cell in the distance direction corresponds to each of the first distance to the ninth distance shown in FIG.

第1折り返し補正部8は、気象クラッタマップ作成部7で作成された気象クラッタマップに基づいて、気象クラッタのドップラー周波数に対する折り返し補正処理を距離方向および方位方向について実行し、第1補正後ドップラー周波数および気象クラッタの強度を出力する。第1折り返し補正部8は、ビーム走査毎にこの折り返し補正処理を実行する。また、第1補正後ドップラー周波数および気象クラッタの強度は、ビーム走査毎にデータ記憶部11に記憶される。   The first aliasing correction unit 8 executes the aliasing correction process for the Doppler frequency of the weather clutter based on the weather clutter map created by the weather clutter map creation unit 7 for the distance direction and the azimuth direction, and the first corrected Doppler frequency And the intensity of weather clutter. The first aliasing correction unit 8 executes this aliasing correction process for each beam scan. The first corrected Doppler frequency and the intensity of the weather clutter are stored in the data storage unit 11 for each beam scan.

気象クラッタの折り返しには、距離の折り返しおよびドップラー周波数の折り返しがある。距離の折り返しは、前に送信された高周波送信パルス信号が非常に遠い物体で反射し、次に送信された高周波送信パルス信号に混入して受信される場合に発生するものである。
ここで、距離の折り返しを抑制するために、パルス繰り返し時間を大きくすることが考えられるが、この場合には、ドップラー周波数の折り返しが発生する。
Weather clutter folding includes distance folding and Doppler frequency folding. The return of the distance occurs when the previously transmitted high-frequency transmission pulse signal is reflected by a very far object and mixed with the next transmitted high-frequency transmission pulse signal.
Here, in order to suppress the return of the distance, it is conceivable to increase the pulse repetition time, but in this case, the return of the Doppler frequency occurs.

ドップラー周波数の折り返しは、次式(2)で表される最大折り返し周波数fmaxよりも高いドップラー周波数が正確に認識されないために発生するものである。
式(2)において、Tはパルス繰り返し時間を示している。
The Doppler frequency aliasing occurs because a Doppler frequency higher than the maximum aliasing frequency f max expressed by the following equation (2) is not accurately recognized.
In equation (2), T s indicates the pulse repetition time.

Figure 2009109457
Figure 2009109457

すなわち、図4に示すように、真のドップラー周波数に対して、ドップラー周波数の測定値は、−fmax〜fmaxの値をとる。そのため、最大折り返し周波数fmaxよりも大きなドップラー周波数を直接測定することができない。
ドップラー周波数の折り返しが発生すると、真のドップラー周波数に最大折り返し周波数fmaxを偶数倍した値が加算または減算されるので、折り返しの発生前後で距離方向に不連続な領域が生じる。
That is, as shown in FIG. 4, the measured value of the Doppler frequency takes a value of −f max to f max with respect to the true Doppler frequency. Therefore, it is impossible to directly measure a Doppler frequency higher than the maximum folding frequency f max .
When the Doppler frequency aliasing occurs, a value obtained by multiplying the true Doppler frequency by an even multiple of the maximum aliasing frequency f max is added or subtracted, so that a discontinuous region occurs in the distance direction before and after the aliasing occurs.

そこで、第1折り返し補正部8は、ドップラー周波数の折り返しを検出するとともに、ドップラー周波数に対する折り返し補正処理を距離方向について実行する。
以下、図5を参照しながら、第1折り返し補正部8によるドップラー周波数に対する距離方向の折り返し補正処理について説明する。
Therefore, the first aliasing correction unit 8 detects the aliasing of the Doppler frequency and executes the aliasing correction process for the Doppler frequency in the distance direction.
Hereinafter, the folding correction process in the distance direction with respect to the Doppler frequency by the first folding correction unit 8 will be described with reference to FIG.

図5は、この発明の実施の形態1に係るドップラー周波数に対する距離方向の折り返し補正処理を示す説明図である。
図5において、まず、第1折り返し補正部8は、距離毎に気象クラッタの強度とあらかじめ設定されたしきい値とを比較し、気象クラッタの強度がしきい値よりも小さい場合に、その気象クラッタのドップラー周波数および強度をノイズとして除去する。また、第1折り返し補正部8は、残った気象クラッタの強度のうち、強度が最大の距離(レンジセル)を検出し、この距離におけるドップラー周波数を基準(真)値とする。以下、このドップラー周波数の基準値を基準ドップラー周波数と表記する。
FIG. 5 is an explanatory diagram showing the aliasing correction processing in the distance direction with respect to the Doppler frequency according to Embodiment 1 of the present invention.
In FIG. 5, first, the first aliasing correction unit 8 compares the intensity of the weather clutter with a preset threshold value for each distance, and if the intensity of the weather clutter is smaller than the threshold value, the weather Remove the Doppler frequency and intensity of the clutter as noise. Further, the first aliasing correction unit 8 detects a distance (range cell) having the maximum intensity from the intensity of the remaining weather clutter, and sets the Doppler frequency at this distance as a reference (true) value. Hereinafter, the reference value of the Doppler frequency is referred to as a reference Doppler frequency.

続いて、第1折り返し補正部8は、着目する距離におけるドップラー周波数fと、ドップラー周波数fに隣接するドップラー周波数fとを比較する。第1折り返し補正部8は、両者の差の絶対値が最大折り返し周波数fmaxよりも大きくなる場合、すなわち次式(3)が成り立つ場合に、ドップラー周波数の折り返しが発生していると判断する。 Subsequently, the first aliasing correction unit 8 compares the Doppler frequency f n at the target distance and the Doppler frequency f m adjacent to the Doppler frequency f n . The first aliasing correction unit 8 determines that the aliasing of the Doppler frequency has occurred when the absolute value of the difference between the two is greater than the maximum aliasing frequency fmax , that is, when the following equation (3) is satisfied.

|f−f|>fmax ・・・(3) | F n −f m |> f max (3)

第1折り返し補正部8は、最初に上記基準ドップラー周波数と、基準ドップラー周波数に隣接するドップラー周波数とを比較する。このとき、着目する距離におけるドップラー周波数fが基準ドップラー周波数となり、ドップラー周波数fに隣接するドップラー周波数fが基準ドップラー周波数に隣接するドップラー周波数となる。 The first aliasing correction unit 8 first compares the reference Doppler frequency with a Doppler frequency adjacent to the reference Doppler frequency. In this case, the Doppler frequency f n at the distance of interest is a reference Doppler frequency, the Doppler frequency Doppler frequency f m which is adjacent to the Doppler frequency f n is adjacent to a reference Doppler frequency.

第1折り返し補正部8は、ドップラー周波数の折り返しが発生していないと判断した場合には、ドップラー周波数fを補正しない。
また、第1折り返し補正部8は、ドップラー周波数の折り返しが発生していると判断した場合には、ドップラー周波数fがドップラー周波数fに最も近くなるように、ドップラー周波数fに最大折り返し周波数fmaxを偶数倍した値を加算または減算して、次式(4)で表されるドップラー周波数fに補正する。
式(4)において、sは1または−1である。
The first reflection correcting unit 8, when wrapping Doppler frequency is determined not to occur does not correct the Doppler frequency f m.
The first reflection correcting unit 8, when wrapping Doppler frequency is determined to have occurred, as the Doppler frequency f m is closest to the Doppler frequency f n, the maximum folding frequency in the Doppler frequency f m the f max by adding or subtracting even multiple values, corrects the Doppler frequency f p represented by the following formula (4).
In the formula (4), s is 1 or -1.

=f+2sfmax ・・・(4) f p = f m + 2sf max (4)

次に、第1折り返し補正部8は、基準ドップラー周波数に隣接するドップラー周波数と、基準ドップラー周波数に隣接するドップラー周波数にさらに隣接する(すなわち、基準ドップラー周波数の2つ隣の)ドップラー周波数とを比較する。このとき、着目する距離におけるドップラー周波数fは、基準ドップラー周波数に隣接するドップラー周波数(fまたはf)となる。第1折り返し補正部8は、以下、同様にして次々と隣接する距離方向のドップラー周波数を補正する。 Next, the first aliasing correction unit 8 compares the Doppler frequency adjacent to the reference Doppler frequency with the Doppler frequency further adjacent to the Doppler frequency adjacent to the reference Doppler frequency (that is, adjacent to the reference Doppler frequency). To do. In this case, the Doppler frequency f n at the distance of interest is a Doppler frequency that is adjacent to a reference Doppler frequency (f m or f p). Hereinafter, the first aliasing correction unit 8 similarly corrects the Doppler frequencies in the distance direction adjacent to each other in the same manner.

なお、気象クラッタのドップラー周波数および強度の測定値には、ノイズ等による誤差が含まれているので、ドップラー周波数および強度は、不連続な値をとることがある。この不連続な値を、ドップラー周波数の折り返しが発生していると誤って判断すると、ドップラー周波数を適切に補正することができない。
そのため、ドップラー周波数および強度の測定値に対して、既存の処理であるメジアンフィルタ処理等のノイズ除去処理を事前に施す必要がある。
The measured values of the Doppler frequency and intensity of the weather clutter include errors due to noise and the like, so the Doppler frequency and intensity may take discontinuous values. If this discontinuous value is erroneously determined that the Doppler frequency aliasing has occurred, the Doppler frequency cannot be corrected appropriately.
For this reason, it is necessary to apply noise removal processing such as median filter processing, which is an existing processing, to the Doppler frequency and intensity measurement values in advance.

以上、図5を参照しながら、ドップラー周波数に対する距離方向の折り返し補正処理について説明した。しかしながら、距離方向の折り返し補正処理を実行した場合であっても、隣接する方位で折り返しの向き(プラス方向側およびマイナス方向側のいずれかの向き)が一致しないので、方位方向でドップラー周波数が不連続な値をとる可能性がある。   As described above, the aliasing correction processing in the distance direction with respect to the Doppler frequency has been described with reference to FIG. However, even when the aliasing correction processing in the distance direction is executed, the aliasing direction (either the plus direction side or the minus direction side) does not match in the adjacent orientation, so the Doppler frequency does not match in the orientation direction. Possible to take continuous values.

そこで、第1折り返し補正部8は、距離方向についてドップラー周波数に対する折り返し補正処理を実行した後に、方位方向についてドップラー周波数に対する折り返し補正処理を実行する。
以下、図6を参照しながら、第1折り返し補正部8によるドップラー周波数に対する方位方向の折り返し補正処理について説明する。
Therefore, the first aliasing correction unit 8 executes the aliasing correction process for the Doppler frequency in the azimuth direction after executing the aliasing correction process for the Doppler frequency in the distance direction.
Hereinafter, the aliasing correction processing in the azimuth direction with respect to the Doppler frequency by the first aliasing correction unit 8 will be described with reference to FIG.

図6は、この発明の実施の形態1に係るドップラー周波数に対する方位方向の折り返し補正処理を示す説明図である。
図6において、まず、第1折り返し補正部8は、着目する第i方位において気象クラッタの強度が最大となる距離におけるドップラー周波数fと、第i方位に隣接する第i+1方位において気象クラッタの強度が最大となる距離におけるドップラー周波数fi+1とを比較する。第1折り返し補正部8は、両者の差の絶対値が最大折り返し周波数fmaxよりも大きくなる場合、すなわち次式(5)が成り立つ場合に、ドップラー周波数の折り返しが発生していると判断する。
FIG. 6 is an explanatory diagram showing a fold-back correction process in the azimuth direction with respect to the Doppler frequency according to Embodiment 1 of the present invention.
6, first, the first reflection correcting unit 8, and the Doppler frequency f i at the distance which the intensity of the weather clutter in the i direction of interest is maximized, the intensity of the weather clutter in the (i + 1) orientation adjacent to the i-th orientation Is compared with the Doppler frequency f i + 1 at the distance where becomes the maximum. The first aliasing correction unit 8 determines that the aliasing of the Doppler frequency has occurred when the absolute value of the difference between the two is greater than the maximum aliasing frequency f max , that is, when the following equation (5) is satisfied.

|fi+1−f|>fmax ・・・(5) | F i + 1 −f i |> f max (5)

第1折り返し補正部8は、ドップラー周波数の折り返しが発生していないと判断した場合には、ドップラー周波数fi+1を補正しない。
また、第1折り返し補正部8は、ドップラー周波数の折り返しが発生していると判断した場合には、第i+1方位のドップラー周波数が第i方位のドップラー周波数に最も近くなるように、第i+1方位のドップラー周波数を補正する。具体的には、第i+1方位における各距離のドップラー周波数のそれぞれに、最大折り返し周波数fmaxを偶数倍した値を加算または減算する。第i+1方位における各距離の補正後のドップラー周波数fi+1(x)は、次式(6)で表される。
When determining that no Doppler frequency aliasing has occurred, the first aliasing correction unit 8 does not correct the Doppler frequency fi + 1 .
In addition, when the first aliasing correction unit 8 determines that the aliasing of the Doppler frequency has occurred, the first aliasing correction unit 8 has the i + 1 orientation so that the Doppler frequency in the i + 1 orientation is closest to the Doppler frequency in the i orientation. Correct the Doppler frequency. Specifically, a value obtained by multiplying the maximum folding frequency f max by an even number is added to or subtracted from each of the Doppler frequencies at each distance in the (i + 1) th azimuth. The Doppler frequency f i + 1 (x) after correction of each distance in the ( i + 1 ) th azimuth is expressed by the following equation (6).

i+1(x)=fi+1,x+2tfmax ・・・(6) f i + 1 (x) = f i + 1, x + 2tf max (6)

式(6)において、fi+1,xは第i+1方位における第x距離のドップラー周波数を示している。また、式(6)において、tは1または−1である。
次に、第1折り返し補正部8は、第i+1方位において気象クラッタの強度が最大となる距離におけるドップラー周波数fi+1と、第i+1方位に隣接する第i+2方位において気象クラッタの強度が最大となる距離におけるドップラー周波数fi+2とを比較する。第1折り返し補正部8は、以下、同様にして次々と隣接する方位方向のドップラー周波数を補正する。
In Equation (6), f i + 1, x represents the Doppler frequency of the xth distance in the (i + 1) th azimuth. Moreover, in Formula (6), t is 1 or -1.
Next, the first aliasing correction unit 8 performs the Doppler frequency f i + 1 at the distance where the intensity of the weather clutter is maximum in the i + 1 direction and the distance where the intensity of the weather clutter is maximum in the i + 2 direction adjacent to the i + 1 direction. And the Doppler frequency f i + 2 at. The first aliasing correction unit 8 corrects Doppler frequencies in the adjacent azimuth directions one after another in the same manner.

切り換えスイッチ9は、第1折り返し補正部8から出力された第1補正後ドップラー周波数の進路を切り換えて、第2折り返し補正部10およびフィルタ設定部13のいずれか一方に出力する。
切り換えスイッチ9をフィルタ設定部13側に切り換えた場合には、ビーム走査毎に折り返し補正処理が実行された気象クラッタのドップラー周波数、および気象クラッタの強度に基づいて、気象クラッタを抑圧するためのフィルタ係数が算出される。
一方、切り換えスイッチ9を第2折り返し補正部10側に切り換えた場合には、ビーム走査間で折り返し補正処理が実行され、平均化処理された気象クラッタのドップラー周波数(ドップラー周波数平均値)、および平均化処理された気象クラッタの強度に基づいて、フィルタ係数が算出される。
The changeover switch 9 switches the course of the first post-correction Doppler frequency output from the first aliasing correction unit 8 and outputs it to one of the second aliasing correction unit 10 and the filter setting unit 13.
When the changeover switch 9 is switched to the filter setting unit 13 side, a filter for suppressing weather clutter based on the Doppler frequency of the weather clutter for which the aliasing correction processing is executed for each beam scan and the intensity of the weather clutter. A coefficient is calculated.
On the other hand, when the changeover switch 9 is switched to the second aliasing correction unit 10 side, aliasing correction processing is executed between the beam scans, and the averaged Doppler frequency (Doppler frequency average value) of the weather clutter is averaged. The filter coefficient is calculated based on the intensity of the weather clutter that has been processed.

第2折り返し補正部10は、ビーム走査間で、第1折り返し補正部8から出力された第1補正後ドップラー周波数に対する折り返し補正処理を実行し、第2補正後ドップラー周波数および気象クラッタの強度を出力する。
前述の第1折り返し補正部8は、気象クラッタのドップラー周波数に対する距離方向および方位方向の折り返し補正処理をビーム走査毎に実行している。しかしながら、ビーム走査毎に折り返しの向きがプラス方向側またはマイナス方向側となり、ビーム走査間で見ると、第1補正後ドップラー周波数の折り返しが発生している可能性がある。
The second aliasing correction unit 10 performs aliasing correction processing on the first corrected Doppler frequency output from the first aliasing correction unit 8 between beam scans, and outputs the second corrected Doppler frequency and the intensity of weather clutter. To do.
The first aliasing correction unit 8 performs the aliasing correction processing in the distance direction and the azimuth direction with respect to the Doppler frequency of the weather clutter for each beam scan. However, the turn-back direction is positive or negative for each beam scan, and there is a possibility that the first corrected Doppler frequency is turned back when viewed between the beam scans.

そこで、第2折り返し補正部10は、ビーム走査間で第1補正後ドップラー周波数に対する折り返し補正処理を実行する。
以下、図7を参照しながら、第2折り返し補正部10による第1補正後ドップラー周波数に対する折り返し補正処理について説明する。
図7は、この発明の実施の形態1に係る第1補正後ドップラー周波数に対する折り返し補正処理を示す説明図である。
Therefore, the second aliasing correction unit 10 executes aliasing correction processing for the first post-correction Doppler frequency between beam scans.
Hereinafter, the aliasing correction process for the first corrected Doppler frequency by the second aliasing correction unit 10 will be described with reference to FIG.
FIG. 7 is an explanatory diagram showing the aliasing correction process for the first corrected Doppler frequency according to the first embodiment of the present invention.

図7において、まず、第2折り返し補正部10は、j回目のビーム走査における第1補正後ドップラー周波数と、j+1回目のビーム走査における第1補正後ドップラー周波数とを方位毎に比較する。具体的には、j回目のビーム走査において第i方位の気象クラッタの強度が最大となる距離における第1補正後ドップラー周波数fj,iと、j+1回目のビーム走査において第i方位の気象クラッタの強度が最大となる距離における第1補正後ドップラー周波数fj+1,iとを比較する。
第2折り返し補正部10は、両者の差の絶対値が最大折り返し周波数fmaxよりも大きくなる場合、すなわち次式(7)が成り立つ場合に、第1補正後ドップラー周波数の折り返しが発生していると判断する。
In FIG. 7, first, the second folding correction unit 10 compares the first corrected Doppler frequency in the j-th beam scanning and the first corrected Doppler frequency in the j + 1-th beam scanning for each direction. Specifically, the first corrected Doppler frequency f j, i at a distance at which the intensity of the weather clutter in the i-th direction is maximum in the j-th beam scan, and the weather clutter in the i-th direction in the j + 1-th beam scan. The first corrected Doppler frequency f j + 1, i at the distance where the intensity is maximum is compared.
The second aliasing correction unit 10 causes aliasing of the first post-correction Doppler frequency when the absolute value of the difference between the two becomes larger than the maximum aliasing frequency f max , that is, when the following equation (7) holds: Judge.

|fj+1,i−fj,i|>fmax ・・・(7) | F j + 1, i −f j, i |> f max (7)

第2折り返し補正部10は、第1補正後ドップラー周波数の折り返しが発生していないと判断した場合には、第1補正後ドップラー周波数fj+1,iを補正しない。
また、第2折り返し補正部10は、第1補正後ドップラー周波数の折り返しが発生していると判断した場合には、j+1回目のビーム走査における第i方位の第1補正後ドップラー周波数が、j回目のビーム走査における第i方位の第1補正後ドップラー周波数に最も近くなるように、j+1回目のビーム走査における第i方位の第1補正後ドップラー周波数を補正する。具体的には、j+1回目のビーム走査の第i方位における各距離の第1補正後ドップラー周波数のそれぞれに、最大折り返し周波数fmaxを偶数倍した値を加算または減算する。j+1回目のビーム走査の第i方位における各距離の補正後の第1補正後ドップラー周波数fj+1,i(x)は、次式(8)で表される。
The second aliasing correction unit 10 does not correct the first corrected Doppler frequency f j + 1, i when determining that the aliasing of the first corrected Doppler frequency has not occurred.
Further, when the second aliasing correction unit 10 determines that aliasing of the first corrected Doppler frequency has occurred, the first corrected Doppler frequency in the i-th azimuth in the j + 1-th beam scanning is the jth time. The first corrected Doppler frequency in the i-th azimuth in the j + 1-th beam scanning is corrected so as to be closest to the first corrected Doppler frequency in the i-th direction in the beam scanning. Specifically, a value obtained by multiplying the maximum folding frequency f max by an even number is added to or subtracted from each first corrected Doppler frequency at each distance in the i-th azimuth of the (j + 1) th beam scanning. The first corrected Doppler frequency f j + 1, i (x) after correction of each distance in the i-th azimuth of the j + 1-th beam scanning is expressed by the following equation (8).

i+1,i(x)=fj+1,i,x+2ufmax ・・・(8) f i + 1, i (x) = f j + 1, i, x + 2uf max (8)

式(8)において、fj+1,i,xは、j+1回目のビーム走査における第i方位、第x距離の第1補正後ドップラー周波数を示している。また、式(8)において、uは1または−1である。
なお、第2折り返し補正部10によるビーム走査間の折り返し補正処理は、ビーム走査回数が2回以上になると実行することができる。すなわち、第2折り返し補正部10は、ビーム走査回数が2回以上になると、データ記憶部11に記憶された第1補正後ドップラー周波数を読み出す。また、第2折り返し補正部10は、第1補正後ドップラー周波数に対する折り返し補正処理を実行した後、第2補正後ドップラー周波数としてデータ記憶部11に格納する。
In the equation (8), f j + 1, i, x indicate the first corrected Doppler frequency of the i-th azimuth and the x-th distance in the j + 1-th beam scanning. Moreover, in Formula (8), u is 1 or -1.
Note that the folding correction process between beam scans by the second folding correction unit 10 can be executed when the number of beam scans is two or more. That is, the second aliasing correction unit 10 reads the first post-correction Doppler frequency stored in the data storage unit 11 when the number of beam scans is two or more. The second aliasing correction unit 10 executes the aliasing correction process on the first corrected Doppler frequency, and then stores the second aliasing correction frequency in the data storage unit 11 as the second corrected Doppler frequency.

データ記憶部11は、第1折り返し補正部8から出力された第1補正後ドップラー周波数および気象クラッタの強度、および第2折り返し補正部10から出力された第2補正後ドップラー周波数を記憶する。   The data storage unit 11 stores the first corrected Doppler frequency and weather clutter intensity output from the first aliasing correction unit 8 and the second corrected Doppler frequency output from the second aliasing correction unit 10.

平均処理部12は、第2折り返し補正部10から出力された第2補正後ドップラー周波数および気象クラッタの強度に対する平均化処理を実行して、ドップラー周波数平均値および強度平均値を出力する。平均処理部12は、ビーム走査間での同一距離および同一方位における第2補正後ドップラー周波数どうし、および強度どうしを平均する。
なお、平均化処理に用いるビーム走査回数は、固定の値であるとする。平均化処理に用いるビーム走査回数が多いほど、信号処理に用いるパルスヒット数が増加するので、気象クラッタのドップラー周波数の測定精度を向上させることができる。
The average processing unit 12 performs an averaging process on the second corrected Doppler frequency output from the second aliasing correction unit 10 and the intensity of the weather clutter, and outputs the Doppler frequency average value and the intensity average value. The average processing unit 12 averages the second post-correction Doppler frequencies and the intensities at the same distance and the same direction between the beam scans.
It is assumed that the number of beam scans used for the averaging process is a fixed value. As the number of beam scans used for averaging processing increases, the number of pulse hits used for signal processing increases, so that the measurement accuracy of the Doppler frequency of the weather clutter can be improved.

フィルタ設定部13は、受信部5からのデジタル化された低周波受信パルス信号に含まれる気象クラッタを抑圧するためのフィルタ係数を算出する。
フィルタ設定部13は、第1折り返し補正部8からの第1補正後ドップラー周波数および気象クラッタの強度が入力された場合には、第1補正後ドップラー周波数に減衰特性を有し、気象クラッタの強度以上の減衰量特性を有する帯域制限フィルタのフィルタ係数を算出する。
The filter setting unit 13 calculates a filter coefficient for suppressing weather clutter included in the digitized low frequency reception pulse signal from the reception unit 5.
When the first corrected Doppler frequency and the intensity of the weather clutter are input from the first aliasing correction unit 8, the filter setting unit 13 has an attenuation characteristic in the first corrected Doppler frequency, and the intensity of the weather clutter The filter coefficient of the band limiting filter having the above attenuation characteristic is calculated.

また、フィルタ設定部13は、平均処理部12からのドップラー周波数平均値および強度平均値が入力された場合には、ドップラー周波数平均値に減衰特性を有し、強度平均値以上の減衰量特性を有する帯域制限フィルタのフィルタ係数を算出する。
ここで、例えば公知のFIR型またはIIR型のデジタルフィルタにおけるフィルタ係数を設定することにより、気象クラッタを抑圧するフィルタが設定される。
Further, when the Doppler frequency average value and the intensity average value are input from the average processing unit 12, the filter setting unit 13 has an attenuation characteristic in the Doppler frequency average value, and has an attenuation characteristic that is equal to or greater than the intensity average value. A filter coefficient of the band limiting filter that is included is calculated.
Here, for example, a filter for suppressing weather clutter is set by setting a filter coefficient in a known FIR type or IIR type digital filter.

クラッタ抑圧処理部14は、フィルタ設定部13で算出された帯域制限フィルタのフィルタ係数に基づいて、受信部5からのデジタル化された低周波受信パルス信号に対するフィルタ処理を実行し、気象クラッタを抑圧した抑圧後受信パルス信号を出力する。   The clutter suppression processing unit 14 performs filter processing on the digitized low-frequency received pulse signal from the reception unit 5 based on the filter coefficient of the band limiting filter calculated by the filter setting unit 13 to suppress weather clutter. The post-suppression received pulse signal is output.

すなわち、切り換えスイッチ9が第2折り返し補正部10側に切り換えられている場合、クラッタ抑圧処理部14には、n回目までの過去の複数回のビーム走査に対してビーム走査間で平均化された第2補正後ドップラー周波数および気象クラッタの強度に応じた帯域制限フィルタが設定されている。また、クラッタ抑圧処理部14は、この帯域制限フィルタを用いて、n+1回目のビーム走査で受信部5から出力されるデジタル化された低周波受信パルス信号に対するフィルタ処理を実行する。この結果、デジタル化された低周波受信パルス信号から気象クラッタ成分が抑圧される。   That is, when the changeover switch 9 is switched to the second aliasing correction unit 10 side, the clutter suppression processing unit 14 averages the beam scans over a plurality of past beam scans up to the nth time. A band limiting filter according to the second corrected Doppler frequency and the intensity of the weather clutter is set. Further, the clutter suppression processing unit 14 performs filter processing on the digitized low-frequency reception pulse signal output from the reception unit 5 in the (n + 1) th beam scanning using the band limiting filter. As a result, the weather clutter component is suppressed from the digitized low frequency received pulse signal.

目標検出部15は、クラッタ抑圧処理部14から出力された抑圧後受信パルス信号に基づいて、目標物を検出する。
目標検出部15は、例えば公知のCFAR(Constant False Alarm Rate)処理等の方法によって、目標物を検出する。
The target detection unit 15 detects a target based on the received pulse signal after suppression output from the clutter suppression processing unit 14.
The target detection unit 15 detects a target by a method such as a known CFAR (Constant False Alarm Rate) process.

この発明の実施の形態1に係るレーダ装置によれば、第1折り返し補正部8は、気象クラッタマップ作成部7からの気象クラッタマップに基づいて、ビーム走査毎に気象クラッタのドップラー周波数に対する折り返し補正処理を実行し、第1補正後ドップラー周波数を出力する。第2折り返し補正部10は、ビーム走査間で第1補正後ドップラー周波数に対する折り返し補正処理を実行し、第2補正後ドップラー周波数を出力する。また、平均処理部12は、第2補正後ドップラー周波数に対する平均化処理を実行して、ドップラー周波数平均値を出力する。クラッタ抑圧処理部14は、ドップラー周波数平均値に基づいてフィルタ設定部13で算出されたフィルタ係数に応じて、デジタル化された低周波受信パルス信号に対するフィルタ処理を実行し、気象クラッタを抑圧する。
そのため、ビーム走査間で正確なドップラー周波数平均値を得ることができるので、気象クラッタのドップラー周波数を高精度に測定することができる。
また、このドップラー周波数平均値に基づいて帯域制限フィルタのフィルタ係数を正確に算出することにより、気象クラッタを適切に抑圧することができる。
According to the radar apparatus according to Embodiment 1 of the present invention, the first aliasing correction unit 8 performs the aliasing correction for the Doppler frequency of the weather clutter for each beam scanning based on the weather clutter map from the weather clutter map creation unit 7. The process is executed and the first corrected Doppler frequency is output. The second aliasing correction unit 10 performs aliasing correction processing for the first corrected Doppler frequency between beam scans, and outputs the second corrected Doppler frequency. Moreover, the average process part 12 performs the averaging process with respect to the 2nd corrected Doppler frequency, and outputs a Doppler frequency average value. The clutter suppression processing unit 14 performs a filtering process on the digitized low-frequency received pulse signal according to the filter coefficient calculated by the filter setting unit 13 based on the Doppler frequency average value, and suppresses weather clutter.
Therefore, since an accurate average Doppler frequency value can be obtained between beam scans, the Doppler frequency of weather clutter can be measured with high accuracy.
Further, by accurately calculating the filter coefficient of the band limiting filter based on this Doppler frequency average value, weather clutter can be appropriately suppressed.

実施の形態2.
上記実施の形態1では、平均処理部12において、第2補正後ドップラー周波数および気象クラッタの強度に対する平均化処理に用いるビーム走査回数は、固定の値であるとした。しかしながら、これに限定されず、平均化処理に用いるビーム走査回数は、目標検出部15における目標物の検出結果に応じて任意に設定されてもよい。
以下、目標検出部15における目標物の検出結果に応じて平均化処理に用いるビーム走査回数を可変設定する処理について説明する。
Embodiment 2. FIG.
In Embodiment 1 described above, the number of beam scans used for the averaging process for the second corrected Doppler frequency and the intensity of the weather clutter in the averaging processing unit 12 is a fixed value. However, the present invention is not limited to this, and the number of beam scans used for the averaging process may be arbitrarily set according to the detection result of the target in the target detection unit 15.
Hereinafter, a process of variably setting the number of beam scans used for the averaging process according to the target detection result in the target detection unit 15 will be described.

図8は、この発明の実施の形態2に係るレーダ装置を示すブロック構成図である。
図8において、このレーダ装置は、図1に示した平均処理部12に代えて、平均処理部12Aを備えている。また、このレーダ装置は、目標検出部15における目標物の検出結果に応じて、第2補正後ドップラー周波数および気象クラッタの強度に対する平均化処理に用いるビーム走査回数を設定し、平均処理部12Aに出力する平均回数制御部16を備えている。
その他の構成については、前述の実施の形態1と同様であり、その説明を省略する。
FIG. 8 is a block diagram showing a radar apparatus according to Embodiment 2 of the present invention.
In FIG. 8, the radar apparatus includes an average processing unit 12A instead of the average processing unit 12 shown in FIG. In addition, the radar apparatus sets the number of beam scans used for the averaging process with respect to the second corrected Doppler frequency and the intensity of the weather clutter according to the detection result of the target in the target detection unit 15, and the average processing unit 12A An average number control unit 16 for output is provided.
Other configurations are the same as those of the first embodiment, and the description thereof is omitted.

一般的に、気象クラッタを抑制するためのフィルタ係数は、気象クラッタのドップラー周波数および強度の測定精度が高いほど正確な値となる。気象クラッタのドップラー周波数および強度の測定精度を向上させるための方法として、平均化処理に用いるパルスヒット数を増加させることが考えられる。パルスヒット数を増加させるためには、平均化処理に用いるビーム走査回数を増加させればよい。   In general, the filter coefficient for suppressing weather clutter becomes more accurate as the measurement accuracy of the Doppler frequency and intensity of the weather clutter is higher. As a method for improving the measurement accuracy of the Doppler frequency and intensity of the weather clutter, it is conceivable to increase the number of pulse hits used in the averaging process. In order to increase the number of pulse hits, the number of beam scans used in the averaging process may be increased.

そこで、平均回数制御部16は、目標検出部15における目標物の検出結果、具体的には、クラッタ抑圧処理部14からの抑圧後受信パルス信号における気象クラッタの消え残りの大きさに応じて、平均化処理に用いるビーム走査回数を設定する。平均処理部12Aは、平均回数制御部16で設定されたビーム走査回数に応じて、第2補正後ドップラー周波数および気象クラッタの強度に対する平均化処理を実行する。   Therefore, the average number control unit 16 determines the target detection result in the target detection unit 15, specifically, according to the magnitude of the remaining uncleared weather clutter in the received pulse signal after suppression from the clutter suppression processing unit 14. The number of beam scans used for the averaging process is set. The average processing unit 12A performs an averaging process on the second corrected Doppler frequency and the intensity of the weather clutter according to the number of beam scans set by the average number control unit 16.

このため、抑圧後受信パルス信号における所望信号対クラッタ信号の電力比で表される気象クラッタの消え残りの大きさを、平均化処理に用いるビーム走査回数と対応付けてあらかじめ設定しておく。平均回数制御部16は、所望信号対クラッタ信号の電力比を実際に測定した結果に基づいて、平均化処理に用いるビーム走査回数を設定する。なお、所望信号とは、目標物を示す信号であり、クラッタ信号とは、気象クラッタを示す信号である。   For this reason, the magnitude of the unerased weather clutter represented by the power ratio of the desired signal to the clutter signal in the received pulse signal after suppression is set in advance in association with the number of beam scans used in the averaging process. The average number control unit 16 sets the number of beam scans used for the averaging process based on the result of actually measuring the power ratio of the desired signal to the clutter signal. The desired signal is a signal indicating a target, and the clutter signal is a signal indicating weather clutter.

ここで、抑圧後受信パルス信号における所望信号対クラッタ信号の電力比と、平均化処理に用いるビーム走査回数との関係を図9に例示する。
図9において、所望信号対クラッタ信号の電力比が小さいほど、所望信号が気象クラッタに埋もれているので、平均化処理に用いるビーム走査回数が大きく設定されている。
なお、所望信号対クラッタ信号の電力比と、平均化処理に用いるビーム走査回数との関係は、図9に示されたものに限定されず、レーダ装置の仕様や電波環境等に応じて変更される。
Here, FIG. 9 illustrates the relationship between the power ratio of the desired signal to the clutter signal in the received pulse signal after suppression and the number of beam scans used for the averaging process.
In FIG. 9, as the power ratio of the desired signal to the clutter signal is smaller, the desired signal is buried in the weather clutter, so the number of beam scans used for the averaging process is set larger.
Note that the relationship between the power ratio of the desired signal to the clutter signal and the number of beam scans used for the averaging process is not limited to that shown in FIG. 9, but is changed according to the specifications of the radar device, the radio wave environment, and the like. The

この発明の実施の形態2に係るレーダ装置によれば、平均回数制御部16は、目標検出部15における目標物の検出結果、すなわち気象クラッタの消え残りの大きさを示す所望信号対クラッタ信号の電力比に応じて、平均化処理に用いるビーム走査回数を設定する。
そのため、気象クラッタのドップラー周波数および強度の測定精度を向上させることができるとともに、気象クラッタを抑圧するための帯域フィルタのフィルタ係数をより正確に算出することができる。
According to the radar apparatus according to the second embodiment of the present invention, the average number control unit 16 detects the target signal in the target detection unit 15, that is, the desired signal vs. the clutter signal indicating the magnitude of the remaining weather clutter. The number of beam scans used for the averaging process is set according to the power ratio.
Therefore, the measurement accuracy of the Doppler frequency and intensity of the weather clutter can be improved, and the filter coefficient of the band filter for suppressing the weather clutter can be calculated more accurately.

実施の形態3.
上記実施の形態1では、第2折り返し補正部10は、隣接する2回のビーム走査の間で第1補正後ドップラー周波数の折り返しの発生を判断し、第1補正後ドップラー周波数に対する折り返し補正処理を実行していた。しかしながら、これに限定されず、2回以上の複数回のビーム走査の間で第1補正後ドップラー周波数の折り返し方向を一度に判定し、この折り返し方向に基づいて、第1補正後ドップラー周波数に対する折り返し補正処理を実行してもよい。
以下、複数回のビーム走査の間で第1補正後ドップラー周波数の折り返し方向を判定し、第1補正後ドップラー周波数に対する折り返し補正処理を実行する処理について説明する。
Embodiment 3 FIG.
In the first embodiment, the second aliasing correction unit 10 determines the occurrence of aliasing of the first post-correction Doppler frequency between two adjacent beam scans, and performs the aliasing correction processing for the first post-correction Doppler frequency. Was running. However, the present invention is not limited to this, and the folding direction of the first corrected Doppler frequency is determined at a time between two or more beam scans, and the folding with respect to the first corrected Doppler frequency is performed based on the folding direction. Correction processing may be executed.
Hereinafter, a process of determining the folding direction of the first corrected Doppler frequency between a plurality of beam scans and executing the folding correction process on the first corrected Doppler frequency will be described.

図10は、この発明の実施の形態3に係るレーダ装置を示すブロック構成図である。
図10において、このレーダ装置は、図1に示した第2折り返し補正部10に代えて、第2折り返し補正部10Bを備えている。また、このレーダ装置は、折り返し方向判定部17を備えている。
その他の構成については、前述の実施の形態1と同様であり、その説明を省略する。
FIG. 10 is a block diagram showing a radar apparatus according to Embodiment 3 of the present invention.
In FIG. 10, the radar apparatus includes a second folding correction unit 10B instead of the second folding correction unit 10 shown in FIG. The radar apparatus also includes a turn-back direction determination unit 17.
Other configurations are the same as those of the first embodiment, and the description thereof is omitted.

折り返し方向判定部17は、複数回のビーム走査の間で第1補正後ドップラー周波数の折り返し方向を判定する。また、折り返し方向判定部17は、基準のビーム走査として決定したビーム走査の回数番号、および第1補正後ドップラー周波数に対する折り返し補正処理を実行する必要のあるビーム走査の回数番号を第2折り返し補正部10Bに出力する。第2折り返し補正部10Bは、基準として決定したビーム走査の値を基準にして、折り返し補正処理が必要と判断されたビーム走査についてのみ、第1補正後ドップラー周波数に対する折り返し補正処理を実行する。   The folding direction determination unit 17 determines the folding direction of the first post-correction Doppler frequency between a plurality of beam scans. The folding direction determination unit 17 also sets the number of times of beam scanning determined as the reference beam scanning and the number of times of beam scanning that needs to be subjected to folding correction processing for the first post-correction Doppler frequency to the second folding correction unit. Output to 10B. The second aliasing correction unit 10B executes the aliasing correction process for the first post-correction Doppler frequency only for the beam scanning determined to require the aliasing correction process with reference to the beam scanning value determined as the reference.

以下、図11を参照しながら、折り返し方向判定部17の動作について詳細に説明する。
図11は、この発明の実施の形態3に係る折り返し方向判定部17の動作を示す説明図である。ここでは、3回のビーム走査の間で第1補正後ドップラー周波数の折り返し方向を判定し、基準のビーム走査および折り返し補正処理を実行する必要のあるビーム走査を決定する処理について説明する。
Hereinafter, the operation of the folding direction determination unit 17 will be described in detail with reference to FIG.
FIG. 11 is an explanatory diagram showing the operation of the folding direction determination unit 17 according to the third embodiment of the present invention. Here, the process of determining the folding direction of the first post-correction Doppler frequency between the three beam scannings and determining the beam scanning that needs to perform the reference beam scanning and the folding correction process will be described.

図11において、まず、折り返し方向判定部17には、データ記憶部11から、複数回のビーム走査における第1補正後ドップラー周波数および気象クラッタの強度が入力される。続いて、折り返し方向判定部17は、各々のビーム走査における気象クラッタの強度を比較し、最大の強度が含まれているビーム走査を基準のビーム走査として決定する。
具体的には、j−1回目、j回目、j+1回目のビーム走査のうち、強度が最大の気象クラッタを測定したビーム走査を基準のビーム走査とする。ここでは、j−1回目の第p方位において気象クラッタの強度が最大になったとし、j−1回目のビーム走査を基準のビーム走査と決定する。
In FIG. 11, first, the first corrected Doppler frequency and the intensity of weather clutter in a plurality of beam scans are input from the data storage unit 11 to the folding direction determination unit 17. Subsequently, the return direction determination unit 17 compares the intensity of the weather clutter in each beam scan, and determines the beam scan including the maximum intensity as the reference beam scan.
Specifically, of the j−1th, jth, and j + 1th beam scans, the beam scan that measures the weather clutter having the maximum intensity is set as the reference beam scan. Here, it is assumed that the intensity of the weather clutter is maximized in the (j−1) th p-th direction, and the j−1th beam scan is determined as the reference beam scan.

続いて、折り返し方向判定部17は、基準のビーム走査において気象クラッタの強度が最大である方位と同じ方位に限定して、基準のビーム走査と他のビーム走査との間で、気象クラッタの強度が最大となる距離における第1補正後ドップラー周波数を比較する。折り返し方向判定部17は、互いの差の絶対値が前述した式(2)に示す最大折り返し周波数fmaxよりも大きくなる場合に、他のビーム走査を、折り返し補正処理を実行する必要のあるビーム走査と決定する。 Subsequently, the folding direction determination unit 17 limits the intensity of the weather clutter between the reference beam scan and another beam scan by limiting to the same direction as the direction in which the intensity of the weather clutter is maximum in the reference beam scan. Compare the first corrected Doppler frequency at a distance where becomes the maximum. When the absolute value of the difference between the folding direction determination unit 17 is larger than the maximum folding frequency f max shown in the above-described equation (2), the beam that needs to perform the folding correction process for another beam scan. Determine to scan.

具体的には、j−1回目のビーム走査において第p方位の気象クラッタの強度が最大となる距離における第1補正後ドップラー周波数fj−1,pと、j回目およびj+1回目のビーム走査において第p方位の気象クラッタの強度が最大となる距離における第1補正後ドップラー周波数fj,p、fj+1,pとをそれぞれ比較する。ここでは、j回目のビーム走査において第1補正後ドップラー周波数の折り返しが発生しているので、j回目のビーム走査を、折り返し補正処理を実行する必要のあるビーム走査と決定する。 Specifically, the first corrected Doppler frequency f j−1, p at the distance at which the intensity of the weather clutter in the p-th direction is the maximum in the j− 1th beam scan, and the jth and j + 1th beam scans. The first corrected Doppler frequencies f j, p , f j + 1, p at the distance at which the intensity of the weather clutter in the p-th direction is maximum are respectively compared. Here, since aliasing of the first post-correction Doppler frequency occurs in the j-th beam scanning, the j-th beam scanning is determined to be a beam scanning that needs to be subjected to aliasing correction processing.

この発明の実施の形態3に係るレーダ装置によれば、折り返し方向判定部17は、複数回のビーム走査の間で第1補正後ドップラー周波数の折り返し方向を判定する。また、折り返し方向判定部17は、基準のビーム走査として決定したビーム走査の回数番号、および第1補正後ドップラー周波数に対する折り返し補正処理を実行する必要のあるビーム走査の回数番号を第2折り返し補正部10Bに出力する。
そのため、2回以上の複数回のビーム走査について第1補正後ドップラー周波数に対する折り返し補正処理を実行する場合に、第1補正後ドップラー周波数の折り返し方向および折り返し補正処理を実行する必要のあるビーム走査が分かった上で、折り返し補正処理を実行することができる。
したがって、折り返し方向を逐一決定して折り返し補正処理を実行する場合と比較して、手順を簡素化することができる。
According to the radar apparatus according to Embodiment 3 of the present invention, the folding direction determination unit 17 determines the folding direction of the first post-correction Doppler frequency between a plurality of beam scans. The folding direction determination unit 17 also sets the number of times of beam scanning determined as the reference beam scanning and the number of times of beam scanning that needs to be subjected to folding correction processing for the first post-correction Doppler frequency to the second folding correction unit. Output to 10B.
Therefore, when the folding correction processing for the first post-correction Doppler frequency is executed for a plurality of beam scans of two or more times, the beam scanning that needs to perform the folding direction and the folding correction processing of the first post-correction Doppler frequency is performed. After knowing, the aliasing correction processing can be executed.
Therefore, the procedure can be simplified as compared with the case where the return direction is determined and the return correction process is executed.

実施の形態4.
上記実施の形態1では、フィルタ設定部13において、ドップラー周波数平均値に減衰特性を有し、強度平均値以上の減衰量特性を有する帯域制限フィルタのフィルタ係数が算出されるとした。しかしながら、これに限定されず、気象クラッタの周波数幅を測定して、気象クラッタを抑圧するフィルタの帯域幅を設定してもよい。
以下、気象クラッタの周波数幅を測定して、気象クラッタを抑圧するフィルタの帯域幅を設定する処理について説明する。
Embodiment 4 FIG.
In the first embodiment, the filter setting unit 13 calculates the filter coefficient of the band limiting filter having the attenuation characteristic in the average Doppler frequency and the attenuation characteristic higher than the average intensity value. However, the present invention is not limited to this, and the bandwidth of a filter that suppresses weather clutter may be set by measuring the frequency width of weather clutter.
Hereinafter, a process for measuring the frequency width of the weather clutter and setting the bandwidth of the filter for suppressing the weather clutter will be described.

図12は、この発明の実施の形態4に係るレーダ装置を示すブロック構成図である。
図12において、このレーダ装置は、図1に示した気象クラッタ検出部6、気象クラッタマップ作成部7および平均処理部12に代えて、気象クラッタ検出部6C、気象クラッタマップ作成部7Cおよび平均処理部12Cを備えている。また、このレーダ装置は、フィルタ形状算出部19を備えている。
その他の構成については、前述の実施の形態1と同様であり、その説明を省略する。
FIG. 12 is a block diagram showing a radar apparatus according to Embodiment 4 of the present invention.
In FIG. 12, this radar apparatus replaces the weather clutter detection unit 6, the weather clutter map creation unit 7 and the average processing unit 12 shown in FIG. 1 with a weather clutter detection unit 6C, a weather clutter map creation unit 7C and an average process. 12C is provided. The radar apparatus also includes a filter shape calculation unit 19.
Other configurations are the same as those of the first embodiment, and the description thereof is omitted.

気象クラッタ検出部6Cは、受信部5でデジタル化されてパルス繰り返し時間毎に並べ替えられた低周波送信パルス信号から、気象クラッタを検出し、気象クラッタのドップラー周波数、強度および周波数幅を出力する。
ここで、気象クラッタのドップラー周波数および強度は、前述した実施の形態1で示した方法によって算出される。一方、気象クラッタの周波数幅は、低周波送信パルス信号を使用して、前述したパルスペア方式により、次式(9)のように算出される。なお、気象クラッタの周波数幅は、ビーム走査毎に、各距離および各方位について測定される。
The weather clutter detection unit 6C detects the weather clutter from the low-frequency transmission pulse signals digitized by the reception unit 5 and rearranged for each pulse repetition time, and outputs the Doppler frequency, intensity, and frequency width of the weather clutter. .
Here, the Doppler frequency and intensity of the weather clutter are calculated by the method described in the first embodiment. On the other hand, the frequency width of the weather clutter is calculated as shown in the following equation (9) by the above-described pulse pair method using the low-frequency transmission pulse signal. Note that the frequency width of the weather clutter is measured for each distance and each direction for each beam scan.

Figure 2009109457
Figure 2009109457

式(9)において、σ(x)は第x距離における気象クラッタの周波数幅、Tはパルス繰り返し時間、R(x,T)は第x距離においてパルス繰り返し時間Tをタイムラグとした受信パルス列の自己相関、R(x,0)は第x距離における受信パルス列の自己相関を示している。また、記号| |は絶対値を示している。
なお、式(9)におけるR(x,T)およびR(x,0)は、次式(10)および次式(11)でそれぞれ表される。
In Equation (9), σ f (x) is the frequency width of the weather clutter at the x-th distance, T s is the pulse repetition time, and R (x, T s ) is the pulse lag time T s at the x-th distance. The autocorrelation of the received pulse train, R (x, 0), indicates the autocorrelation of the received pulse train at the xth distance. The symbol || represents an absolute value.
In addition, R (x, T s ) and R (x, 0) in the formula (9) are represented by the following formula (10) and the following formula (11), respectively.

Figure 2009109457
Figure 2009109457
Figure 2009109457
Figure 2009109457

式(10)および式(11)において、Tはパルス繰返し時間、argは複素数値の位相角、Mはパルスヒット数、v(x,k)は第x距離におけるkヒット目の低周波受信パルス信号、*は複素共役を示している。
気象クラッタマップ作成部7Cは、気象クラッタ検出部6Cから出力された気象クラッタのドップラー周波数、強度および周波数幅のデータに基づいて、ドップラー周波数、強度および周波数幅の距離特性を示す気象クラッタマップを、方位毎に作成する。
ここで、第i方位における周波数幅の距離特性を示す気象クラッタマップを図13に例示する。
In equations (10) and (11), T s is the pulse repetition time, arg is the complex-valued phase angle, M is the number of pulse hits, and v (x, k) is the low frequency reception of the k hit at the xth distance. A pulse signal, * indicates a complex conjugate.
The weather clutter map creation unit 7C generates a weather clutter map indicating the distance characteristics of the Doppler frequency, intensity, and frequency width based on the Doppler frequency, intensity, and frequency width data of the weather clutter output from the weather clutter detection unit 6C. Create for each direction.
Here, FIG. 13 illustrates a weather clutter map showing the distance characteristics of the frequency width in the i-th direction.

また、平均処理部12Cは、第2折り返し補正部10から出力された第2補正後ドップラー周波数および気象クラッタの強度に対する平均化処理を実行するとともに、複数のビーム走査間で気象クラッタの周波数幅に対する平均化処理を実行し、周波数幅平均値を出力する。   The average processing unit 12C performs an averaging process on the second corrected Doppler frequency and the intensity of the weather clutter output from the second aliasing correction unit 10, and performs the frequency processing on the frequency width of the weather clutter between a plurality of beam scans. Averaging processing is executed and the average frequency width is output.

フィルタ形状算出部19は、平均処理部12Cから出力されるドップラー周波数平均値、強度平均値および周波数幅平均値に基づいて、気象クラッタを抑圧するためのフィルタ形状を算出する。すなわち、ドップラー周波数平均値に一致するようにフィルタの中心周波数を設定し、強度平均値以上となるように減衰量を設定し、周波数幅平均値以上となるようにフィルタの帯域幅を設定する。なお、フィルタの設定は、各方位について距離毎に行われる。   The filter shape calculation unit 19 calculates a filter shape for suppressing weather clutter based on the Doppler frequency average value, the intensity average value, and the frequency width average value output from the average processing unit 12C. That is, the center frequency of the filter is set so as to match the Doppler frequency average value, the attenuation amount is set so as to be equal to or greater than the intensity average value, and the filter bandwidth is set so as to be equal to or greater than the frequency width average value. The filter is set for each distance in each direction.

また、フィルタ設定部13は、前述した実施の形態1で示した方法で、帯域制限フィルタのフィルタ係数を算出する。すなわち、例えば公知のFIR型またはIIR型のデジタルフィルタにおけるフィルタ係数を設定し、気象クラッタを抑圧するフィルタを設定する。   Further, the filter setting unit 13 calculates the filter coefficient of the band limiting filter by the method described in the first embodiment. That is, for example, a filter coefficient in a known FIR type or IIR type digital filter is set, and a filter that suppresses weather clutter is set.

この発明の実施の形態4に係るレーダ装置によれば、気象クラッタ検出部6Cは、気象クラッタの周波数幅を測定する。また、フィルタ形状算出部19は、平均処理部12Cで平均化された周波数幅平均値を用いて、気象クラッタを抑圧するフィルタの帯域幅を設定する。
そのため、気象クラッタを抑圧するフィルタの形状を正確に設定することができる。
According to the radar apparatus according to Embodiment 4 of the present invention, the weather clutter detection unit 6C measures the frequency width of the weather clutter. The filter shape calculation unit 19 sets the bandwidth of the filter that suppresses the weather clutter using the average frequency width averaged by the average processing unit 12C.
Therefore, the shape of the filter that suppresses the weather clutter can be set accurately.

実施の形態5.
上記実施の形態4では、フィルタ形状算出部19は、気象クラッタの強度および周波数幅について、平均処理部12Cで平均化された強度平均値および周波数幅平均値を用いてフィルタ形状を算出するとした。しかしながら、これに限定されず、気象クラッタの強度および周波数幅の最大値を検出し、これらの最大値を用いてフィルタ形状を算出してもよい。
以下、気象クラッタの強度および周波数幅の最大値を用いてフィルタ形状を算出する処理について説明する。
Embodiment 5 FIG.
In the fourth embodiment, the filter shape calculation unit 19 calculates the filter shape using the intensity average value and the frequency width average value averaged by the average processing unit 12C for the intensity and frequency width of the weather clutter. However, the present invention is not limited to this, and the maximum values of the intensity and frequency width of the weather clutter may be detected, and the filter shape may be calculated using these maximum values.
Hereinafter, the process of calculating the filter shape using the maximum value of the intensity and frequency width of the weather clutter will be described.

図14は、この発明の実施の形態5に係るレーダ装置を示すブロック構成図である。
図14において、このレーダ装置は、図12に示した平均処理部12Cおよびフィルタ形状算出部19に代えて、平均処理部12Dおよびフィルタ形状算出部19Dを備えている。また、このレーダ装置は、最大値検出部20を備えている。
その他の構成については、前述の実施の形態4と同様であり、その説明を省略する。
FIG. 14 is a block diagram showing a radar apparatus according to Embodiment 5 of the present invention.
In FIG. 14, the radar apparatus includes an average processing unit 12D and a filter shape calculation unit 19D instead of the average processing unit 12C and the filter shape calculation unit 19 shown in FIG. In addition, this radar apparatus includes a maximum value detection unit 20.
Other configurations are the same as those of the fourth embodiment described above, and a description thereof is omitted.

平均処理部12Dは、第2折り返し補正部10から出力された第2補正後ドップラー周波数に対する平均化処理を実行して、ドップラー周波数平均値を出力する。
最大値検出部20は、複数のビーム走査(ビーム走査間で平均化処理を実行しようとした場合のビーム走査回数)で、各方位について距離毎に気象クラッタの強度および周波数幅の最大値を検出し、強度最大値および周波数幅最大値を出力する。
The average processing unit 12D performs an averaging process on the second corrected Doppler frequency output from the second aliasing correction unit 10 and outputs an average Doppler frequency value.
The maximum value detection unit 20 detects the maximum values of the intensity and frequency width of the weather clutter for each distance in each direction in a plurality of beam scans (the number of beam scans when the averaging process is to be executed between beam scans). The maximum intensity value and the maximum frequency width value are output.

フィルタ形状算出部19は、平均処理部12Dから出力されるドップラー周波数平均値と、最大値検出部20から出力される強度最大値および周波数幅最大値とに基づいて、気象クラッタを抑圧するためのフィルタ形状を算出する。すなわち、ドップラー周波数平均値に一致するようにフィルタの中心周波数を設定し、強度最大値以上となるように減衰量を設定し、周波数幅最大値以上となるようにフィルタの帯域幅を設定する。なお、フィルタの設定は、各方位について距離毎に行われる。   The filter shape calculation unit 19 suppresses weather clutter based on the Doppler frequency average value output from the average processing unit 12D and the maximum intensity value and frequency width maximum value output from the maximum value detection unit 20. The filter shape is calculated. That is, the center frequency of the filter is set so as to coincide with the average Doppler frequency, the attenuation is set so as to be equal to or greater than the maximum intensity value, and the bandwidth of the filter is set so as to be equal to or greater than the maximum frequency width. The filter is set for each distance in each direction.

この発明の実施の形態5に係るレーダ装置によれば、フィルタ形状算出部19Dは、最大値検出部20で検出された強度最大値および周波数幅最大値を用いて、気象クラッタを抑圧するフィルタの減衰量および帯域幅を設定する。
そのため、フィルタの減衰特性を十分に確保することができ、気象クラッタの抑圧性能を向上させることができる。
According to the radar apparatus according to Embodiment 5 of the present invention, the filter shape calculation unit 19D uses a maximum intensity value and a maximum frequency width value detected by the maximum value detection unit 20 to use a filter that suppresses weather clutter. Set attenuation and bandwidth.
Therefore, it is possible to sufficiently secure the attenuation characteristic of the filter and improve the suppression performance of the weather clutter.

なお、フィルタ形状算出部19は、気象クラッタの時間的な変動に応じて、強度平均値および周波数幅平均値と、強度最大値および周波数幅最大値との何れかを選択してフィルタ形状を算出してもよい。
この場合も、上記実施の形態4または5と同様の効果を奏することができる。
The filter shape calculation unit 19 calculates the filter shape by selecting either the intensity average value and the frequency width average value, the intensity maximum value, and the frequency width maximum value according to the temporal variation of the weather clutter. May be.
Also in this case, the same effect as in the fourth or fifth embodiment can be obtained.

なお、上記実施の形態1〜5では、レンジセルの数を9としているが、これに限定されない。レンジセルの数は、任意の値に設定することができる。
この場合も、上記実施の形態1〜5と同様の効果を奏することができる。
In the first to fifth embodiments, the number of range cells is nine. However, the present invention is not limited to this. The number of range cells can be set to an arbitrary value.
Also in this case, the same effects as those of the first to fifth embodiments can be obtained.

この発明の実施の形態1に係るレーダ装置を示すブロック構成図である。It is a block block diagram which shows the radar apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る低周波受信パルス信号のAD変換処理を示す説明図である。It is explanatory drawing which shows the AD conversion process of the low frequency received pulse signal which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る気象クラッタマップを示す説明図である。It is explanatory drawing which shows the weather clutter map which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る真のドップラー周波数とドップラー周波数の測定値との関係を示す説明図である。It is explanatory drawing which shows the relationship between the true Doppler frequency which concerns on Embodiment 1 of this invention, and the measured value of a Doppler frequency. この発明の実施の形態1に係るドップラー周波数に対する距離方向の折り返し補正処理を示す説明図である。It is explanatory drawing which shows the folding correction process of the distance direction with respect to the Doppler frequency which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係るドップラー周波数に対する方位方向の折り返し補正処理を示す説明図である。It is explanatory drawing which shows the folding correction process of the azimuth | direction direction with respect to the Doppler frequency which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る第1補正後ドップラー周波数に対する折り返し補正処理を示す説明図である。It is explanatory drawing which shows the aliasing correction process with respect to the 1st post-correction Doppler frequency which concerns on Embodiment 1 of this invention. この発明の実施の形態2に係るレーダ装置を示すブロック構成図である。It is a block block diagram which shows the radar apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る所望信号対クラッタ信号の電力比と、平均化処理に用いるビーム走査回数との関係を示す説明図である。It is explanatory drawing which shows the relationship between the power ratio of the desired signal and the clutter signal which concerns on Embodiment 2 of this invention, and the frequency | count of beam scanning used for an averaging process. この発明の実施の形態3に係るレーダ装置を示すブロック構成図である。It is a block block diagram which shows the radar apparatus which concerns on Embodiment 3 of this invention. この発明の実施の形態3に係る折り返し方向判定部の動作を示す説明図である。It is explanatory drawing which shows operation | movement of the folding | turning direction determination part which concerns on Embodiment 3 of this invention. この発明の実施の形態4に係るレーダ装置を示すブロック構成図である。It is a block block diagram which shows the radar apparatus which concerns on Embodiment 4 of this invention. この発明の実施の形態4に係る周波数幅の距離特性を示す説明図である。It is explanatory drawing which shows the distance characteristic of the frequency width which concerns on Embodiment 4 of this invention. この発明の実施の形態5に係るレーダ装置を示すブロック構成図である。It is a block block diagram which shows the radar apparatus concerning Embodiment 5 of this invention.

符号の説明Explanation of symbols

1 局部発振部、2 送信部、4 空中線部、5 受信部、6、6C 気象クラッタ検出部、8 第1折り返し補正部、10、10B 第2折り返し補正部、12、12A、12C、12D 平均処理部、13 フィルタ設定部、14 クラッタ抑圧処理部、15 目標検出部、16 平均回数制御部、17 折り返し方向判定部、19、19D フィルタ形状算出部、20 最大値検出部。   DESCRIPTION OF SYMBOLS 1 Local oscillation part, 2 Transmission part, 4 Antenna part, 5 Reception part, 6, 6C Weather clutter detection part, 8 1st return correction part, 10, 10B 2nd return correction part, 12, 12A, 12C, 12D Average process Unit, 13 filter setting unit, 14 clutter suppression processing unit, 15 target detection unit, 16 average number control unit, 17 folding direction determination unit, 19, 19D filter shape calculation unit, 20 maximum value detection unit.

Claims (5)

周波数変換を行うための局部発振信号を発生して出力する局部発振部と、
低周波送信パルス信号を発生するとともに、前記局部発振信号を用いて前記低周波送信パルス信号を高周波送信パルス信号に周波数変換して出力する送信部と、
前記高周波送信パルス信号を空間に放射するとともに、目標物で反射した前記高周波送信パルス信号を高周波受信パルス信号として受信する空中線部と、
前記局部発振信号を用いて前記高周波受信パルス信号を低周波受信パルス信号に周波数変換するとともに、前記低周波受信パルス信号をAD変換して出力する受信部と、
デジタル化された前記低周波受信パルス信号から不要な信号である気象クラッタを検出し、前記気象クラッタのドップラー周波数および強度を出力する気象クラッタ検出部と、
ビームがレーダ覆域を走査するビーム走査毎に、前記ドップラー周波数に対する折り返し補正処理を実行し、第1補正後ドップラー周波数を出力する第1折り返し補正部と、
複数回の前記ビーム走査の間であるビーム走査間で、前記第1補正後ドップラー周波数に対する折り返し補正処理を実行し、第2補正後ドップラー周波数を出力する第2折り返し補正部と、
前記第2補正後ドップラー周波数に対する平均化処理を実行して、ドップラー周波数平均値を出力する平均処理部と、
前記第1補正後ドップラー周波数および前記ドップラー周波数平均値のいずれか一方に基づいて、前記デジタル化された低周波受信パルス信号に含まれる前記気象クラッタを抑圧するためのフィルタ係数を算出するフィルタ設定部と、
前記フィルタ係数に基づいて、前記デジタル化された低周波受信パルス信号に対するフィルタ処理を実行し、前記気象クラッタを抑圧した抑圧後受信パルス信号を出力するクラッタ抑圧処理部と、
前記抑圧後受信パルス信号に基づいて前記目標物を検出する目標検出部と、
を備えたことを特徴とするレーダ装置。
A local oscillator for generating and outputting a local oscillation signal for frequency conversion;
A transmitter that generates a low-frequency transmission pulse signal, converts the frequency of the low-frequency transmission pulse signal into a high-frequency transmission pulse signal using the local oscillation signal, and outputs the high-frequency transmission pulse signal;
An antenna unit that radiates the high-frequency transmission pulse signal into space and receives the high-frequency transmission pulse signal reflected by the target as a high-frequency reception pulse signal;
Using the local oscillation signal, the high-frequency reception pulse signal is frequency-converted into a low-frequency reception pulse signal, and the low-frequency reception pulse signal is AD-converted and output;
Detecting a weather clutter that is an unnecessary signal from the digitized low frequency received pulse signal, and outputting a Doppler frequency and intensity of the weather clutter;
A first aliasing correction unit that performs aliasing correction processing on the Doppler frequency and outputs the first corrected Doppler frequency each time the beam scans the radar coverage;
A second aliasing correction unit that performs aliasing correction processing on the first corrected Doppler frequency between beam scannings that are a plurality of times of the beam scanning, and outputs a second corrected Doppler frequency;
An averaging process unit that performs an averaging process on the second corrected Doppler frequency and outputs an averaged Doppler frequency value;
A filter setting unit that calculates a filter coefficient for suppressing the weather clutter included in the digitized low-frequency received pulse signal based on one of the first corrected Doppler frequency and the average Doppler frequency When,
A clutter suppression processing unit that performs filtering on the digitized low-frequency received pulse signal based on the filter coefficient, and outputs a post-suppression received pulse signal that suppresses the weather clutter;
A target detection unit for detecting the target based on the received pulse signal after suppression;
A radar apparatus comprising:
前記目標検出部における前記目標物の検出結果に応じて、前記平均化処理に用いる前記ビーム走査の回数を設定する平均回数制御部を備え、
前記平均処理部は、前記平均回数制御部で設定された前記ビーム走査の回数に応じて、前記平均化処理を実行することを特徴とする請求項1に記載のレーダ装置。
According to the detection result of the target in the target detection unit, comprising an average number control unit for setting the number of times of the beam scanning used for the averaging process,
The radar apparatus according to claim 1, wherein the average processing unit executes the averaging process according to the number of beam scans set by the average number control unit.
複数回の前記ビーム走査による前記第1補正後ドップラー周波数に基づいて、前記第1補正後ドップラー周波数の折り返し方向を判定する折り返し方向判定部を備え、
前記第2折り返し補正部は、前記折り返し方向判定部で判定された前記折り返し方向に応じて、前記第1補正後ドップラー周波数に対する折り返し補正処理を実行することを特徴とする請求項1に記載のレーダ装置。
A folding direction determination unit that determines a folding direction of the first corrected Doppler frequency based on the first corrected Doppler frequency by the plurality of beam scans;
2. The radar according to claim 1, wherein the second aliasing correction unit executes a aliasing correction process for the first corrected Doppler frequency according to the aliasing direction determined by the aliasing direction determination unit. apparatus.
フィルタ形状算出部を備え、
前記気象クラッタ検出部は、前記気象クラッタの周波数幅を出力し、
前記平均処理部は、前記ビーム走査間で前記周波数幅に対する平均化処理を実行して周波数幅平均値を出力し、
前記フィルタ形状算出部は、前記周波数幅平均値を用いて、前記気象クラッタを抑圧するためのフィルタの形状を算出することを特徴とする請求項1に記載のレーダ装置。
A filter shape calculator,
The weather clutter detection unit outputs a frequency width of the weather clutter,
The average processing unit performs an averaging process on the frequency width between the beam scans and outputs an average frequency width value,
The radar apparatus according to claim 1, wherein the filter shape calculation unit calculates a shape of a filter for suppressing the weather clutter using the average frequency width value.
最大値検出部を備え、
前記気象クラッタ検出部は、前記気象クラッタの強度を出力し、
前記最大値検出部は、前記気象クラッタの強度および前記気象クラッタの周波数幅の最大値を検出して強度最大値および周波数幅最大値を出力し、
前記フィルタ形状算出部は、前記強度最大値および前記周波数幅最大値を用いて、前記気象クラッタを抑圧するためのフィルタの形状を算出することを特徴とする請求項4に記載のレーダ装置。
Equipped with a maximum value detector,
The weather clutter detection unit outputs the intensity of the weather clutter,
The maximum value detecting unit detects the intensity of the weather clutter and the maximum value of the frequency width of the weather clutter and outputs the maximum intensity value and the maximum frequency width value,
The radar apparatus according to claim 4, wherein the filter shape calculation unit calculates a shape of a filter for suppressing the weather clutter using the maximum intensity value and the maximum frequency width value.
JP2007284860A 2007-11-01 2007-11-01 Radar equipment Expired - Fee Related JP5196959B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007284860A JP5196959B2 (en) 2007-11-01 2007-11-01 Radar equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007284860A JP5196959B2 (en) 2007-11-01 2007-11-01 Radar equipment

Publications (2)

Publication Number Publication Date
JP2009109457A true JP2009109457A (en) 2009-05-21
JP5196959B2 JP5196959B2 (en) 2013-05-15

Family

ID=40778076

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007284860A Expired - Fee Related JP5196959B2 (en) 2007-11-01 2007-11-01 Radar equipment

Country Status (1)

Country Link
JP (1) JP5196959B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013024775A (en) * 2011-07-22 2013-02-04 Mitsubishi Electric Corp Radar device
CN111458701A (en) * 2020-04-12 2020-07-28 西安电子工程研究所 Meteorological track restraining method based on track characteristic iterative updating
CN112105952A (en) * 2018-05-09 2020-12-18 古野电气株式会社 Weather radar device, weather observation method, and weather observation program

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62116278A (en) * 1985-11-15 1987-05-27 Nec Corp Weather clutter suppressing circuit
JP2000230976A (en) * 1999-02-12 2000-08-22 Mitsubishi Electric Corp Device and method for doppler velocity detection
JP2001091645A (en) * 1999-09-22 2001-04-06 Toshiba Corp Doppler radar device
JP2001201568A (en) * 2000-01-18 2001-07-27 Mitsubishi Electric Corp Radar signal processor
JP2007322331A (en) * 2006-06-02 2007-12-13 Mitsubishi Electric Corp Radar device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62116278A (en) * 1985-11-15 1987-05-27 Nec Corp Weather clutter suppressing circuit
JP2000230976A (en) * 1999-02-12 2000-08-22 Mitsubishi Electric Corp Device and method for doppler velocity detection
JP2001091645A (en) * 1999-09-22 2001-04-06 Toshiba Corp Doppler radar device
JP2001201568A (en) * 2000-01-18 2001-07-27 Mitsubishi Electric Corp Radar signal processor
JP2007322331A (en) * 2006-06-02 2007-12-13 Mitsubishi Electric Corp Radar device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013024775A (en) * 2011-07-22 2013-02-04 Mitsubishi Electric Corp Radar device
CN112105952A (en) * 2018-05-09 2020-12-18 古野电气株式会社 Weather radar device, weather observation method, and weather observation program
CN111458701A (en) * 2020-04-12 2020-07-28 西安电子工程研究所 Meteorological track restraining method based on track characteristic iterative updating
CN111458701B (en) * 2020-04-12 2022-12-27 西安电子工程研究所 Meteorological track inhibition method based on track characteristic iterative update

Also Published As

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

Similar Documents

Publication Publication Date Title
JP6150474B2 (en) Signal processing apparatus and signal processing method in wind profiler
KR101135982B1 (en) Synchronization method of radar systems for the rejection of interference in FMCW radars
US20140210661A1 (en) Systems and methods for calibrating dual polarization radar systems
JP2007322331A (en) Radar device
EP2290396A2 (en) Doppler radar apparatus and method of calculating Doppler velocity
JP6676151B2 (en) Wave height calculating device, radar device, and wave height calculating method
US9689979B2 (en) Pulse signal setting device, radar apparatus, method and program of setting pulse signal
JP4999592B2 (en) Radar equipment
EP3006956B1 (en) Surface tidal-current estimation device, radar device, surface tidal-current estimation method and surface tidal-current estimation program
JP5196959B2 (en) Radar equipment
JP4836496B2 (en) Radar apparatus and signal processing method thereof
JP5355321B2 (en) Interference wave detection apparatus and interference wave detection method
JP5574907B2 (en) Radar equipment
JP4836497B2 (en) Radar apparatus and signal processing method thereof
JP5925264B2 (en) Radar equipment
JP6147617B2 (en) Radar apparatus and signal processing method thereof
US11927666B2 (en) Signal processing apparatus and signal processing method
JP6132576B2 (en) Signal processing apparatus, radar apparatus, and signal processing method
JP2013113723A (en) Radar system
JP5491788B2 (en) Weather radar system and its precipitation intensity calculation method
RU102391U1 (en) RADAR SYSTEM
US11892535B2 (en) Signal processing apparatus and signal processing method
JP2007212332A (en) Radar system and its signal processing method
JP2013205268A (en) Radar signal processor
JP4999586B2 (en) Radar equipment

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100913

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120425

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120508

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120628

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130108

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130205

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

Free format text: PAYMENT UNTIL: 20160215

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees