JP2011247593A - Image radar device - Google Patents
Image radar device Download PDFInfo
- Publication number
- JP2011247593A JP2011247593A JP2010117656A JP2010117656A JP2011247593A JP 2011247593 A JP2011247593 A JP 2011247593A JP 2010117656 A JP2010117656 A JP 2010117656A JP 2010117656 A JP2010117656 A JP 2010117656A JP 2011247593 A JP2011247593 A JP 2011247593A
- Authority
- JP
- Japan
- Prior art keywords
- range
- phase change
- hit
- unnecessary
- unnecessary phase
- 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.)
- Pending
Links
Images
Abstract
Description
この発明は、観測対象となる目標の電波画像を取得する画像レーダ装置に関し、特に目標と画像レーダ装置との間の距離の時間変化の影響で発生する画像のぼけを、受信信号そのものに基づいて補償する技術に関するものである。 The present invention relates to an image radar apparatus that acquires a radio wave image of a target to be observed, and in particular, blurring of an image that occurs due to the temporal change in the distance between the target and the image radar apparatus is based on the received signal itself. It relates to the technology to compensate.
従来の画像レーダ装置においては、送信機から発生した高周波信号を、送受切換器および送受信アンテナを介して目標に照射し、目標で散乱された高周波信号の一部を、送受信アンテナおよび送受切換器を介して受信機で受信し、受信信号のレンジ方向の分解能をレンジ圧縮により向上させてレンジプロフィールを得る、という一連の処理を、目標と画像レーダ装置との間の相対位置関係を変えながら繰り返し実行し、得られたレンジプロフィールの時間履歴（以下、「レンジヒストリ」という）を画像再生処理することにより目標のレーダ画像を得ている。 In a conventional image radar apparatus, a target is irradiated with a high-frequency signal generated from a transmitter via a transmission / reception switch and a transmission / reception antenna, and a part of the high-frequency signal scattered by the target is transmitted to the transmission / reception antenna and the transmission / reception switch. A series of processing to obtain the range profile by improving the resolution in the range direction of the received signal by range compression and repeatedly changing the relative positional relationship between the target and the image radar device The target radar image is obtained by performing image reproduction processing on the time history of the obtained range profile (hereinafter referred to as “range history”).
レーダ画像を得る代表的なレーダ装置としては、合成開口レーダ（Ｓｙｎｔｈｅｔｉｃ Ａｐｅｒｔｕｒｅ Ｒａｄａｒ：ＳＡＲ）と、逆合成開口レーダ（Ｉｎｖｅｒｓｅ ＳＡＲ：ＩＳＡＲ）とが挙げられる。 Typical radar apparatuses that obtain radar images include a synthetic aperture radar (SAR) and an inverse synthetic aperture radar (Inverse SAR: ISAR).
ＳＡＲは、たとえば地表面などの観測対象を、レーダ装置（レーダプラットフォーム）が空間内を移動しながら観測し、得られたレンジプロフィールを適切に合成することにより、あたかも大口径のアンテナを空間に配置したような効果を得て、高分解能なレーダ画像を得るシステムである。 For example, SAR observes an observation target such as the ground surface while a radar device (radar platform) moves in space, and synthesizes the obtained range profile appropriately, as if placing a large-diameter antenna in space. This is a system that obtains a high-resolution radar image with the above effects.
逆に、ＩＳＡＲは、運動する目標をレーダ装置で観測することにより、ＳＡＲと同様の効果を得て、高分解能なレーダ画像を得るシステムである。この場合、たとえば、目標に固定された座標系内でのレーダ位置の変化を考えれば、ＳＡＲと同じ効果が得られることが分かる。 On the other hand, ISAR is a system that obtains a high-resolution radar image by observing a moving target with a radar device and obtaining the same effect as SAR. In this case, for example, if the change of the radar position in the coordinate system fixed to the target is considered, it can be seen that the same effect as the SAR can be obtained.
ＳＡＲ／ＩＳＡＲの画像再生処理としては、多数の種類が提案されているが、以下では、一例としてレンジドップラー画像再生の概要について説明する。
レンジドップラー画像再生においては、目標上の反射点を、レーダ装置からの距離を表すレンジと、目標とレーダ装置との間の相対運動で発生するドップラー周波数差とで、分離してレーダ画像を得ている。
Many types of SAR / ISAR image reproduction processing have been proposed, but an outline of range Doppler image reproduction will be described below as an example.
In range Doppler image reproduction, the reflection point on the target is separated into a range representing the distance from the radar device and a Doppler frequency difference generated by the relative motion between the target and the radar device to obtain a radar image. ing.
上記レーダ画像は、すでに得られたレンジヒストリを、観測の繰り返し方向にフーリエ変換することにより得られる。
以下、観測繰り返しの順番を「ヒット」といい、第ｈ本目のレンジプロフィールを得たときの観測順番を、「第ｈヒット」という。また、繰り返し方向を、単に「ヒット方向」という。
The radar image is obtained by Fourier transforming the already obtained range history in the direction of repeated observation.
Hereinafter, the observation repetition order is referred to as “hit”, and the observation order when the h-th range profile is obtained is referred to as “h-th hit”. The repeat direction is simply referred to as “hit direction”.
なお、等時間間隔でデータを観測する場合においては、ヒットは、それぞれの観測開始時刻に比例する。
以下の説明においては、等時間間隔の観測を想定している。
When observing data at equal time intervals, hits are proportional to the respective observation start times.
In the following description, observation at equal time intervals is assumed.
この処理は、レンジに直交するクロスレンジ方向の分解能を向上させることから、「クロスレンジ圧縮」と呼ばれる。
ただし、このクロスレンジ圧縮処理においては、目標上の反射点が、目標とレーダ装置との間の相対運動の影響によって、観測中に以下の条件（ａ）、（ｂ）となった場合に、レーダ画像にぼけが生じることが知られている。
This process is called “cross range compression” because it improves the resolution in the cross range direction orthogonal to the range.
However, in this cross-range compression processing, when the reflection point on the target becomes the following conditions (a) and (b) during observation due to the influence of the relative motion between the target and the radar device, It is known that the radar image is blurred.
（ａ）目標上の反射点が、レンジ分解能セルを越えて移動する。
（ｂ）目標上の反射点の位相が、（ドップラー周波数の時間変化を引き起こす）時間に対する２次以上の変化成分を含む。
(A) The reflection point on the target moves beyond the range resolution cell.
(B) The phase of the reflection point on the target includes a second-order or higher-order change component with respect to time (which causes a time change of the Doppler frequency).
そこで、レーダ画像再生処理においては、上記条件（ａ）、（ｂ）を補償するための運動補償処理が必要になる。
特に、上記条件（ａ）、（ｂ）がセンサで計測不可能な場合（たとえば、ＳＡＲにおいて、レーダプラットフォームに搭載する運動センサの精度が低い場合、または、ＩＳＡＲにおいて、運動が未知の目標を観測する場合）には、受信信号そのものから上記条件（ａ）、（ｂ）の補償量を推定して補償するオートフォーカスが必要になる。
Therefore, in the radar image reproduction process, a motion compensation process for compensating the conditions (a) and (b) is required.
In particular, when the above conditions (a) and (b) cannot be measured by a sensor (for example, in SAR, when a motion sensor mounted on a radar platform is low in accuracy, or in ISAR, a target whose motion is unknown is observed. In this case, autofocus is required to estimate and compensate the compensation amounts of the above conditions (a) and (b) from the received signal itself.
画像レーダにおけるオートフォーカスは、以下の補償（Ａ）、（Ｂ）に大別され、各補償（Ａ）、（Ｂ）ごとに、それぞれ様々な方法が提案されている。
（Ａ）条件（ａ）を推定して補償するレンジ補償。
（Ｂ）条件（ｂ）を推定して補償する位相補償。
Autofocus in image radar is roughly classified into the following compensations (A) and (B), and various methods have been proposed for each compensation (A) and (B).
(A) Range compensation for estimating and compensating for condition (a).
(B) Phase compensation for estimating and compensating for the condition (b).
基本的には、まずレンジ補償（Ａ）を行うことにより、各反射点のレンジ分解能セルを越えた移動を除去した後に、１つ以上のレンジ分解能セル上のヒット方向に並ぶ受信信号に基づいて位相補償（Ｂ）を行う。
位相補償（Ｂ）の有力な方法の１つとして、ＰＧＡ（Ｐｈａｓｅ Ｇｒａｄｉｅｎｔ Ａｕｔｏｆｏｃｕｓ）法が挙げられる（たとえば、特許文献１、特許文献２参照）
Basically, the range compensation (A) is first performed to remove the movement of each reflection point beyond the range resolution cell, and then based on the received signals arranged in the hit direction on one or more range resolution cells. Phase compensation (B) is performed.
As one of the effective methods of phase compensation (B), there is a PGA (Phase Gradient Autofocus) method (for example, refer to
ＰＧＡ法においては、注目するレンジ分解能セルに存在する孤立反射点に生じた「位相誤差の一次微分」を計算し、位相誤差の一次微分値を積分することにより、２次以上の位相変化（以下、「不要位相変化」という）を推定する。
特に、複数のレンジ分解能セル孤立反射点の位相誤差の一次微分値を加算することにより、雑音の影響は低減される。
In the PGA method, a “first-order differential of phase error” generated at an isolated reflection point existing in the target range resolution cell is calculated, and a first-order differential value of the phase error is integrated to calculate a second-order or higher-order phase change (hereinafter referred to as “the first-order differential”). , “Unnecessary phase change”).
In particular, the influence of noise is reduced by adding the first derivative of the phase error of a plurality of range resolution cell isolated reflection points.
従来の画像レーダ装置では、特許文献１、特許文献２に記載された位相補償法（ＰＧＡ法）によれば、雑音の影響が低減されるものの、受信信号の性質（たとえば、注目するレンジセルの反射点の数や位置関係、不要位相変化の次数など）によっては、不要位相変化の推定誤差（たとえば、全ヒットにわたる推定誤差や、或るヒット付近のみで増大する局所的な推定誤差）が増大する可能性があるという課題があった。
In the conventional image radar apparatus, according to the phase compensation method (PGA method) described in
この発明は、上記のような課題を解決するためになされたものであり、ＰＧＡ法を用いて不要位相変化を推定する際に、受信信号の条件によっては発生する可能性がある「全ヒットにわたる推定誤差」、または「或るヒット付近のみで増大する局所的な推定誤差」の影響を低減することにより、位相補償精度の向上、ひいては、再生されたレーダ画像の分解能向上を実現した画像レーダ装置を得ることを目的とする。 The present invention has been made in order to solve the above-described problems. When estimating an unnecessary phase change using the PGA method, it may occur depending on the condition of the received signal. An image radar apparatus that realizes an improvement in phase compensation accuracy and, in turn, an improvement in resolution of a reproduced radar image by reducing the influence of “estimation error” or “local estimation error that increases only near a certain hit” The purpose is to obtain.
この発明に係る画像レーダ装置は、目標のレンジプロフィールを取得するレーダ観測器と、レンジプロフィールのレンジヒストリに基づいてレーダ画像を生成するレーダ画像化器とからなる画像レーダ装置であって、レーダ観測器は、高周波信号を発生する送信機と、高周波信号を目標に照射するとともに目標で散乱された高周波信号の一部を受信する送受信アンテナと、送受信アンテナを介して目標からの受信信号を取り込む受信機と、送信機からの高周波信号を送受信アンテナに送出するとともに送受信アンテナからの受信信号を受信機に入力する送受切換器と、受信信号のレンジ方向の分解能をレンジ圧縮により向上させてレンジプロフィールを取得するレンジ圧縮器と、を備え、レンジプロフィールを取得する一連の処理を、目標とレーダ観測器との間の相対位置関係を変えながら繰り返し実行してレンジヒストリを取得し、レーダ画像化器は、レンジヒストリをクロスレンジ圧縮することにより、目標上の反射強度分布をレンジドップラー分布として表現したレーダ画像を得るために、最終的にクロスレンジ圧縮に用いるヒット幅以上のヒット幅のレンジヒストリに対して、目標とレーダ観測器との間の相対運動で発生する観測中の各反射点のレンジ分解能セルの移動を推定して補償するレンジ補償器と、レンジ補償器によるレンジ補償後のレンジヒストリ上の代表的な複数の反射点の位相変化を調べ、レーダ画像のドップラー周波数方向のぼけの原因となるヒットに対する２次以上の不要位相変化の、最終的にクロスレンジ圧縮に用いるヒット幅の範囲での値を推定する不要位相変化推定器と、不要位相変化推定器で得られた不要位相変化に基づき、レンジ補償後のレンジヒストリから、最終的にクロスレンジ圧縮に用いるヒット幅分だけ切出したレンジヒストリに含まれる不要位相変化成分を補償する位相補償回路と、を備えたものである。 An image radar apparatus according to the present invention is an image radar apparatus that includes a radar observation device that acquires a target range profile and a radar imager that generates a radar image based on a range history of the range profile. A transmitter that generates a high-frequency signal; a transmission / reception antenna that irradiates the target with the high-frequency signal and receives a portion of the high-frequency signal scattered by the target; and a receiver that captures a reception signal from the target via the transmission / reception antenna A transmission / reception switch that sends a high-frequency signal from the transmitter to the transmission / reception antenna and inputs a reception signal from the transmission / reception antenna to the receiver, and improves the range direction resolution of the reception signal by range compression. A range compressor for acquiring, and a series of processes for acquiring a range profile as a target The radar imager obtains the range history by repeatedly executing the relative position relationship with the radar observer, and the radar imager converts the range history to the Doppler distribution by cross-compressing the range history. In order to obtain a radar image expressed as, each reflection during observation that occurs due to relative movement between the target and the radar observer for a range history with a hit width that is greater than or equal to the hit width used for cross-range compression A range compensator that estimates and compensates for the movement of a point's range resolution cell, and examines the phase change of multiple representative reflection points on the range history after range compensation by the range compensator, in the Doppler frequency direction of the radar image. Estimate the value of the second or higher-order unnecessary phase change for hits that cause blurs in the range of the hit width used for cross-range compression. Based on the unnecessary phase change estimator and the unnecessary phase change obtained by the unnecessary phase change estimator, it is not included in the range history that is finally extracted from the range history after range compensation by the hit width used for cross range compression. And a phase compensation circuit for compensating for the phase change component.
この発明によれば、位相補償精度を向上させて、再生されたレーダ画像の分解能を向上させることができる。 According to the present invention, it is possible to improve the phase compensation accuracy and improve the resolution of the reproduced radar image.
実施の形態１．
図１はこの発明の実施の形態１に係る画像レーダ装置の全体構成を示すブロック図である。
図１において、画像レーダ装置は、レーダ観測器１と、レーダ画像化器２と、により構成されている。
1 is a block diagram showing the overall configuration of an image radar apparatus according to
In FIG. 1, the image radar apparatus includes a
レーダ観測器１は、目標Ｔのレーダ画像Ｕ（ｒ，ｆ）を生成するために必要なレンジヒストリＳ（レンジプロフィールの時間履歴）を得るための観測を行う。
レーダ画像化器２は、レーダ観測器１で得られたレンジヒストリＳに基づいて、レーダ画像Ｕ（ｒ，ｆ）を生成する。
The
The
図２は図１内のレーダ観測器１の構成を示すブロック図である。
図２において、レーダ観測器１は、高周波信号を発生する送信機３と、送受切換器４と、送受信アンテナ５と、受信信号を増幅および検波する受信機６と、レンジ圧縮器７と、を備えている。
FIG. 2 is a block diagram showing the configuration of the
In FIG. 2, the
送受切換器４は、送信時と受信時とで、信号の流れる方向を切り換える。
送受信アンテナ５は、目標Ｔに向けて電波を照射するとともに、目標Ｔからの散乱電波を受信する。
レンジ圧縮器７は、受信機６で得られた受信信号のレンジ分解能を、送信波形の情報に基づき高分解能化し、高分解能化したレンジプロフィールを取得してレーダ画像化器２に入力する。
The transmission /
The transmission / reception antenna 5 emits radio waves toward the target T and receives scattered radio waves from the target T.
The
図３は図１内のレーダ画像化器２の構成を示すブロック図である。
図３において、レーダ画像化器２は、レンジヒストリ切出器１１と、レンジ補償器１２と、位相補償器１３と、クロスレンジ圧縮器１４と、を備えている。
FIG. 3 is a block diagram showing the configuration of the
In FIG. 3, the
レンジヒストリ切出器１１は、レーダ観測器１から入力されたレンジヒストリＳに基づき、指定された切出し開始ヒットおよび指定されたヒット幅Ｈ０のレンジヒストリＳ０（ｒ，ｈ）を切出す。
レンジ補償器１２は、レンジヒストリ切出器１１で切出されたレンジヒストリＳ０（ｒ，ｈ）のレンジ補償を行う。具体的には、最終的にクロスレンジ圧縮に用いるヒット幅Ｈ以上のヒット幅Ｈ０のレンジヒストリＳ０（ｒ，ｈ）に対して、目標Ｔとレーダ観測器１との間の相対運動で発生する、観測中の各反射点のレンジ分解能セルの移動を推定して補償する。
The
The
位相補償器１３は、レンジヒストリ切出器１１およびレンジ補償器１２を介して入力されたレンジヒストリＳ１（ｒ，ｈ）に基づき、最終的に画像化に用いるヒット範囲のレンジヒストリＳ１（ｒ，ｈ）の位相補償を行う。
クロスレンジ圧縮器１４は、位相補償器１３で得られた位相補償後のレンジヒストリＳ３（ｒ，ｈ）のクロスレンジ圧縮を行い、レーダ画像Ｕ（ｒ，ｆ）を生成する。
The
The
図４は図３内の位相補償器１３の構成を示すブロック図である。
図４において、位相補償器１３は、不要位相変化推定器２１および位相補償回路２２を備えている。
FIG. 4 is a block diagram showing the configuration of the
In FIG. 4, the
不要位相変化推定器２１は、レンジ補償器１２で得られたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化を推定する。
具体的には、不要位相変化推定器２１は、レンジ補償器１２から得られるレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）上の代表的な複数の反射点の位相変化を調べ、レーダ画像Ｕ（ｒ，ｆ）のドップラー周波数方向のぼけの原因となる不要位相変化（ヒットに対する２次以上の変化成分）の、最終的にクロスレンジ圧縮に用いるヒット幅Ｈの範囲での不要位相変化φ２（ｈ）の値を推定して、位相補償回路２２に入力する。
The unnecessary
Specifically, the unnecessary
位相補償回路２２は、不要位相変化推定器２１から得られた不要位相変化φ２（ｈ）に基づき、レンジ補償器１２から得られたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）のうちの、最終的に画像化に用いるヒット範囲のレンジヒストリＳ３（ｒ，ｈ）の不要位相変化を補償する。
The
具体的には、位相補償回路２２は、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）から、最終的にクロスレンジ圧縮に用いるヒット幅Ｈ分だけ切出したレンジヒストリＳ３（ｒ，ｈ）に含まれる不要位相変化成分を補償して、クロスレンジ圧縮器１４に入力する。
Specifically, the
図５は図４内の不要位相変化推定器２１の構成を示すブロック図である。
図５において、不要位相変化推定器２１は、ＰＧＡ推定器３１およびヒット範囲切出器１１１を備えている。
FIG. 5 is a block diagram showing the configuration of the unnecessary
In FIG. 5, the unnecessary
ＰＧＡ推定器３１は、レンジ補償器１２から入力されたレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）を、位相補償法の一種であるＰＧＡ法の原理にしたがって推定し、ヒット範囲切出器１１１に入力する。
ヒット範囲切出器１１１は、ＰＧＡ推定器３１から得られた不要位相変化φ１（ｈ）に基づき、最終的に画像化に用いるヒット範囲のレンジヒストリに対応するヒット範囲のデータを切出して、不要位相変化φ２（ｈ）として位相補償回路２２に入力する。
The
Based on the unnecessary phase change φ1 (h) obtained from the
次に、図６を参照しながら、図１〜図５に示したこの発明の実施の形態１による動作について説明する。
図６は位相補償器１３の処理内容を示す説明図であり、図６（ｂ）〜図６（ｅ）において、共通の横軸はヒットｈを表し、縦軸はレンジｒ、不要位相変化レベルを表している。
Next, the operation according to the first embodiment of the present invention shown in FIGS. 1 to 5 will be described with reference to FIG.
FIG. 6 is an explanatory diagram showing the processing contents of the
図６（ａ）はレンジヒストリＳ１（ｒ，ｈ）を格納する配列イメージを図式的に表している。
図６（ｂ）はＰＧＡ推定器３１で得られた不要位相変化φ１（ｈ）を格納する配列イメージを図式的に表しており、ヒットｈごとに１つの不要位相変化φ１（ｈ）が得られる状態を示している。
FIG. 6A schematically shows an array image storing the range history S1 (r, h).
FIG. 6B schematically shows an arrangement image storing the unnecessary phase change φ1 (h) obtained by the
図６（ｃ）は不要位相変化φ１（ｈ）、φ２（ｈ）の値の具体的な変化イメージを表しており、細実線は不要位相変化φ１（ｈ）、太実線は不要位相変化φ２（ｈ）を示している。
図６（ｄ）は不要位相変化φ２（ｈ）を格納する配列イメージを図式的に表しており、 図６（ｅ）はレンジヒストリＳ０（ｒ，ｈ）を格納する配列イメージを図式的に表している。
FIG. 6C shows a specific change image of the values of the unnecessary phase changes φ1 (h) and φ2 (h), the thin solid line indicates the unnecessary phase change φ1 (h), and the thick solid line indicates the unnecessary phase change φ2 ( h).
FIG. 6D schematically shows an array image for storing the unnecessary phase change φ2 (h), and FIG. 6E schematically shows an array image for storing the range history S0 (r, h). ing.
まず、レーダ観測器１（図２）においては、従来装置と同様に、送信機３から高周波信号を発生し、送受切換器４および送受信アンテナ５を介して、高周波信号を目標Ｔに照射する。
また、レーダ観測器１は、目標Ｔで散乱された高周波信号の一部を、送受信アンテナ５および送受切換器４を介して、受信機６で受信する。
First, in the radar observation device 1 (FIG. 2), a high-frequency signal is generated from the
Further, the
レンジ圧縮器７は、受信信号のレンジ方向の分解能を、一般的なパルス圧縮などで向上させて、レンジプロフィールを得る。
以上のレーダ観測器１による一連処理は、目標Ｔとレーダ観測器１との間の相対位置関係を変えながら繰り返し実行され、これにより、レンジヒストリＳ（レンジプロフィールの時間履歴）が得られる。
The
The series of processes by the
レンジヒストリＳは、ヒットｈと、レンジセル番号（以下、単に「レンジ」、「レンジセル」ともいう）ｒと、を軸とする２次元配列で与えられる。各配列要素には、各ヒットｈにおける、各レンジセルｒの複数の反射点の振幅および位相情報を含む複素数値が格納される。 The range history S is given as a two-dimensional array with a hit h and a range cell number (hereinafter simply referred to as “range” or “range cell”) r as axes. Each array element stores a complex value including amplitude and phase information of a plurality of reflection points of each range cell r in each hit h.
続いて、レーダ画像化器２（図３）においては、レーダ観測器１で得られたレンジヒストリＳを処理して、レーダ画像Ｕ（ｒ，ｆ）を生成する。
まず、レンジヒストリ切出器１１は、レーダ観測器１から入力されたレンジヒストリＳに基づき、指定された切出し開始位置、および、入力されたヒット幅以下の指定されたヒット幅Ｈ０のレンジヒストリＳ０（ｒ，ｈ）を切出して、レンジ補償器１２に入力する。
Subsequently, the radar imager 2 (FIG. 3) processes the range history S obtained by the
First, the
すなわち、レンジヒストリ切出器１１は、レーダ画像再生および運動補償（レーダ画像再生の前処理）で用いるヒット位置およびヒット幅Ｈ０のレンジヒストリＳ０（ｒ，ｈ）を、レーダ観測器１から切出す。
That is, the
ここでは、レンジセル数をＲで表し、レンジｒ＝０、１、・・・、Ｒ−１、とする。
また、図６のように、切出したヒット幅Ｈ０を偶数とし、ヒットｈ＝−Ｈ０／２、−Ｈ０／２＋１、・・・、０、・・・、Ｈ０／２−２、Ｈ０／２−１、と番号付ける。
Here, the number of range cells is represented by R, and the ranges r = 0, 1,..., R−1.
Further, as shown in FIG. 6, the extracted hit width H0 is an even number, and hits h = −H0 / 2, −H0 / 2 + 1,..., 0,..., H0 / 2-2, H0 / 2−. Number one.
なお、ヒット幅Ｈ０が奇数の場合は、ヒットｈ＝−（Ｈ０−１）／２、・・・、０、・・・、（Ｈ０−１）／２とすればよいが、以下では、偶数／奇数の場合分けが煩雑になるので、ヒット幅Ｈ０が偶数の場合のみについて説明する。ヒット幅Ｈ０が奇数の場合も、以下の処理はそのまま適用できるのは言うまでもない。 When the hit width H0 is an odd number, the hit h = − (H0-1) / 2,..., 0,..., (H0-1) / 2 may be used. / Since the case of odd numbers becomes complicated, only the case where the hit width H0 is an even number will be described. It goes without saying that the following processing can be applied as it is even when the hit width H0 is an odd number.
また、レーダ画像化器２（図３）においては、画像再生や運動補償に用いるヒット幅Ｈ以上の幅でレンジヒストリＳを収集し、収集されたレンジヒストリＳから、必要なヒット位置およびヒット幅のデータを切出して画像化するような運用（たとえば、切出し位置を連続的に変えて動画像を得るような場合）を想定して、レーダ画像化器２の初段にレンジヒストリ切出器１１が配置されている。
In the radar imager 2 (FIG. 3), the range history S is collected with a width equal to or larger than the hit width H used for image reproduction and motion compensation, and the required hit position and hit width are collected from the collected range history S. Assuming an operation in which the data is cut out and imaged (for example, when moving images are obtained by continuously changing the cut-out position), the range history cut-out
なお、レーダ観測器１において、画像再生や運動補償に用いるのと同じヒット幅Ｈ０のレンジヒストリのみを収集するような運用も当然考えられ、この場合には、切出処理そのものは不要となる。
ただし、この場合も、レンジヒストリ切出器１１による切出処理は、入力されたレンジヒストリＳと同じヒット幅Ｈ０のレンジヒストリＳ０（ｒ，ｈ）を切出す処理と等価なので、汎用性をもたせる意味で、レーダ画像化器２はレンジヒストリ切出器１１を含むものとする。
It should be noted that the
However, in this case as well, the extraction process by the
続いて、レンジ補償器１２は、レンジヒストリ切出器１１で得られた切出し後のレンジヒストリＳ０（ｒ，ｈ）に対して、既存のレンジ補償処理を適用して、目標Ｔとレーダ観測器１との間の相対運動の影響で発生した、各反射点のレンジセル移動成分を補償する。
なお、具体的なレンジ補償処理は、種々の公知文献（特許第３３６０５６２号公報、特開２００６−３４３２９０号公報など）において参照することができる。
Subsequently, the
Specific range compensation processing can be referred to in various known documents (Japanese Patent No. 3360562, Japanese Patent Laid-Open No. 2006-343290, etc.).
レンジ補償器１２から得られるレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）は、ヒットｈ＝−Ｈ０／２、−Ｈ０／２＋１、・・・、０、・・・、Ｈ０／２−２、Ｈ０／２−１、であり、レンジｒ＝０、１、・・・、Ｒ−１、である。
レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）においては、各反射点のレンジ移動成分が補償されているので、図６（ａ）のように、各反射点は、全ヒットｈにわたって、それぞれ同じレンジセルｒ内に留まる。
The range history S1 (r, h) after the range compensation obtained from the
In the range history S1 (r, h) after the range compensation, the range movement component of each reflection point is compensated, so that each reflection point is the same over the entire hit h as shown in FIG. Stay in range cell r.
次に、位相補償器１３（図４）は、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）上に残存する「各反射点の位相の時間ｔに対する（または、ヒットｈに対する）２次以上の変化成分（不要位相変化φ２（ｈ））」を推定して補償する。 Next, the phase compensator 13 (FIG. 4) remains on the range history S1 (r, h) after the range compensation, “second order or higher with respect to time t of the phase of each reflection point (or hit h)” The change component (unnecessary phase change φ2 (h)) ”is estimated and compensated.
一般に、或る時刻ｔにおける受信信号の位相φ（ｔ）［ｒａｄ］と、瞬時のドップラー周波数ｆｄ（ｔ）［Ｈｚ］との間には、以下の式（１）で表される周知の関係がある。 In general, a known relationship represented by the following expression (1) between the phase φ (t) [rad] of a received signal at a certain time t and the instantaneous Doppler frequency fd (t) [Hz]. There is.
式（１）から明らかなように、ドップラー周波数ｆｄ（ｔ）は、位相φ（ｔ）に２次以上の時間変化成分が含まれれば、１次以上の時間変化が生じるので、各反射点を、レンジとドップラー差とで分離して目標Ｔを画像化する場合、各反射点の像は、ドップラー周波数ｆｄ（ｔ）方向に画像ぼけが発生することになる。
したがって、上記画像ぼけを解消するためには、図３、図４のように、位相補償器１３が必要となる。
As is clear from the equation (1), the Doppler frequency fd (t) has a first-order or higher-order time change if the phase φ (t) includes a second-order or higher-order time change component. When the target T is imaged separately by the range and the Doppler difference, the image of each reflection point is blurred in the direction of the Doppler frequency fd (t).
Therefore, in order to eliminate the image blur, the
なお、不要位相変化φ１（ｈ）は、前述のように、目標Ｔとレーダ観測器１との間の相対位置の変化の影響によって発生するが、相対位置の変化によって発生する見込角の変化が小さい場合には、レーダ観測器１と目標Ｔの重心との間の距離変化（時間に対する２次以上の成分）によって発生するものと見なすことができる。また、レーダ観測器１から目標Ｔまでの距離に比べて、目標Ｔの大きさが十分に小さい場合には、距離変化は全反射点に共通なものと見なすことができる。
Note that the unnecessary phase change φ1 (h) occurs due to the influence of the change in the relative position between the target T and the
したがって、近似的に、不要位相変化φ１（ｈ）も、全反射点で一致するものと見なすことができる。
ここでは、上記近似条件が成立する範囲での適用を想定しており、前述のＰＧＡ法（位相補償法）も、上記近似条件が成立する範囲での適用が想定されている。
Therefore, approximately, the unnecessary phase change φ1 (h) can be regarded as matching at the total reflection point.
Here, it is assumed that the approximation condition is satisfied, and the above-described PGA method (phase compensation method) is also assumed to be applied in the range where the approximation condition is satisfied.
位相補償器１３（図４）に入力されたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）は、不要位相変化推定器２１および位相補償回路２２の両方に入力される。
不要位相変化推定器２１（図５）内のＰＧＡ推定器３１は、公知のＰＧＡ法の原理に基づいて、レンジヒストリＳ１（ｒ，ｈ）の不要位相変化φ１（ｈ）を推定する。
ＰＧＡ法に基づく不要位相変化φ１（ｈ）の推定方法については、たとえば前述の特許文献１に詳述されているので、以下に要点のみを述べる。
The range history S1 (r, h) after the range compensation input to the phase compensator 13 (FIG. 4) is input to both the unnecessary
The
Since the estimation method of the unnecessary phase change φ1 (h) based on the PGA method is described in detail in, for example, the above-mentioned
一般に、目標Ｔ上の或る反射点に着目した場合、その反射点の受信信号ｇ０（ｈ）の位相には、反射点の位相やドップラー位置に関する１次以下の位相変化に画像ぼけの原因となる不要位相変化φｅ（ｈ）が含まれる。 In general, when attention is paid to a certain reflection point on the target T, the phase of the reception signal g0 (h) at the reflection point may cause image blurring due to the phase change of the reflection point or the first-order phase related to the Doppler position. The unnecessary phase change φe (h) is included.
この場合、反射点がぼけているので、反射点のドップラー位置を正確に推定することは困難であるが、反射点のドップラー位置の情報に基づき、受信信号ｇ０（ｈ）から１次の位相変化をできるだけ補償した結果をｇ（ｈ）とする。
ここで、理想の（一次の位相変化が完全に除去された）状態においては、不要位相変化φｅ（ｈ）と補償後の受信信号ｇ（ｈ）との間には、以下の式（２）で示す関係が成立する。
In this case, since the reflection point is blurred, it is difficult to accurately estimate the Doppler position of the reflection point. However, based on the information of the Doppler position of the reflection point, the first-order phase change from the received signal g0 (h). Let g (h) be the result of compensating as much as possible.
Here, in an ideal state (primary phase change is completely removed), the following equation (2) is present between the unnecessary phase change φe (h) and the received signal g (h) after compensation. The relationship indicated by is established.
ＰＧＡ法においては、式（２）の関係を利用して、不要位相変化の一次微分値を取得して、これを積分することにより、最終的な不要位相変化φ１（ｈ）を得る。
特に、受信機６の雑音の影響を低減するために、式（２）の関係を、複数のレンジセルｒについて取得して、これを平均処理することができる。
In the PGA method, a first-order differential value of an unnecessary phase change is obtained using the relationship of the expression (2) and integrated to obtain a final unnecessary phase change φ1 (h).
In particular, in order to reduce the influence of the noise of the
また、同じレンジセルｒに複数の反射点が存在する場合を想定して、周波数軸上でウィンドウをかける処理や、推定された不要位相変化φ１（ｈ）に基づいてレンジヒストリを位相補償したうえで、再度の推定処理を行うシーケンス処理を繰り返すことにより、段階的に不要位相変化φ１（ｈ）の推定精度を向上させる方法なども含まれる。 In addition, assuming that there are a plurality of reflection points in the same range cell r, the phase history is compensated for the phase based on the process of applying a window on the frequency axis and the estimated unnecessary phase change φ1 (h). Also included is a method of improving the estimation accuracy of the unnecessary phase change φ1 (h) step by step by repeating the sequence processing for performing the estimation processing again.
ＰＧＡ推定器３１は、上記公知のＰＧＡ法を用いて、入力されたレンジヒストリＳ１（ｒ，ｈ）に対する不要位相変化φ１（ｈ）を得る。
図６（ｂ）のように、レンジヒストリＳ１（ｒ，ｈ）の不要位相変化φ１（ｈ）については、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）（図６（ａ））と同じように、ヒットｈ＝−Ｈ０／２、−Ｈ０／２＋１、・・・、０、・・・、Ｈ０／２−２、Ｈ０／２−１、のヒット範囲で値が得られる。
The
As shown in FIG. 6B, the unnecessary phase change φ1 (h) of the range history S1 (r, h) is the same as the range history S1 (r, h) after range compensation (FIG. 6A). Further, values are obtained in hit ranges of hits h = −H0 / 2, −H0 / 2 + 1,..., 0,..., H0 / 2-2, H0 / 2-1.
従来装置のＰＧＡ処理では、以上のように得られた不要位相変化φ１（ｈ）を用いて、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）の位相補償が行われている。
公知のＰＧＡ処理は、有力な方法ではあるが、受信信号の性質（たとえば、注目するレンジセルｒの反射点の数や位置関係、不要位相変化φ１（ｈ）の次数など）によっては、不要位相変化φ１（ｈ）の推定誤差（たとえば全ヒットにわたる推定誤差や、或るヒット付近のみで増大する局所的な推定誤差）が増大する可能性がある。
In the PGA processing of the conventional apparatus, the phase compensation of the range history S1 (r, h) after the range compensation is performed using the unnecessary phase change φ1 (h) obtained as described above.
The known PGA process is a powerful method, but depending on the nature of the received signal (for example, the number and positional relationship of the reflection points of the target range cell r, the order of the unnecessary phase change φ1 (h), etc.) There is a possibility that an estimation error of φ1 (h) (for example, an estimation error over all hits or a local estimation error that increases only in the vicinity of a certain hit) increases.
そこで、この発明の実施の形態１（図１〜図６）においては、特に、ヒットｈの中央付近（−Ｈ／２〜Ｈ／２−１）では、不要位相変化φ１（ｈ）の推定精度が比較的高いものの、ヒットｈの端部付近（−Ｈ０／２〜−Ｈ／２−１、Ｈ／２〜Ｈ０／２−１）では、不要位相変化φ１（ｈ）の推定精度が比較的低くなる、という状況を想定する。 Therefore, in the first embodiment (FIGS. 1 to 6) of the present invention, particularly in the vicinity of the center of the hit h (−H / 2 to H / 2-1), the estimation accuracy of the unnecessary phase change φ1 (h) is estimated. However, in the vicinity of the end of the hit h (-H0 / 2 to -H / 2-1, H / 2 to H0 / 2-1), the estimation accuracy of the unnecessary phase change φ1 (h) is relatively high. Assume a situation of lowering.
たとえば、図６（ｂ）において、ヒットの端部付近では、相対運動の影響とは異なる位相変化が発生しているような状況をイメージしている。
このような場合に、全ヒットｈの不要位相変化φ１（ｈ）（推定値）を用いて位相補償を行うと、不要位相変化の推定誤差の影響で画像ぼけが残存し、場合によっては画像ぼけが増大する可能性もある。
For example, in FIG. 6B, a situation is imagined in which a phase change different from the effect of relative motion occurs near the end of the hit.
In such a case, if phase compensation is performed using the unnecessary phase change φ1 (h) (estimated value) of all hits h, an image blur remains due to the effect of an estimation error of the unnecessary phase change, and in some cases, the image blur is caused. May also increase.
したがって、不要位相変化推定器２１は、上記のような推定誤差の特性を仮定して、最終的に画像再生に用いるヒット幅Ｈ以上の広いヒット幅Ｈ０（≧Ｈ）でレンジ補償を行い、ＰＧＡ法による不要位相変化φ１（ｈ）を推定したうえで、不要位相変化φ１（ｈ）のうち、誤差が発生している可能性がある不要な端部の推定結果を捨てて、位相補償に本来必要な不要位相変化φ２（ｈ）を最終的に取得する。
Therefore, the unnecessary
なお、処理順序の都合上、ヒット幅Ｈ０でレンジ補償を行い、不要位相変化φ１（ｈ）を推定するという説明を先にしたが、本来的には、まず、画像化に必要なヒット幅Ｈが定まり、その後に、ヒット幅Ｈ以上の適当なヒット幅Ｈ０が定まることになる。 For convenience of processing order, the explanation has been given that range compensation is performed with the hit width H0 and the unnecessary phase change φ1 (h) is estimated. However, first, the hit width H necessary for imaging is essentially the first. After that, an appropriate hit width H0 greater than or equal to the hit width H is determined.
具体的には、ヒット範囲切出器１１１は、ヒット幅Ｈの不要位相変化φ２（ｈ）として、不要位相変化φ１（ｈ）のうちの中央付近（ｈ＝−Ｈ／２〜Ｈ／２−１）を、以下の式（３）のように抽出して、位相補償回路２２に入力する。
Specifically, the
ヒット幅Ｈの不要位相変化φ２（ｈ）を格納する配列のイメージは、図６（ｄ）に示した通りである。
また、不要位相変化φ２（ｈ）の値の変化のイメージは、図６（ｃ）内の太実線で示した通りである。
An image of an array for storing the unnecessary phase change φ2 (h) of the hit width H is as shown in FIG.
The image of the change in the value of the unnecessary phase change φ2 (h) is as shown by the thick solid line in FIG.
続いて、位相補償回路２２は、不要位相変化推定器２１から最終的に得られた不要位相変化φ２（ｈ）を用いて、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）の位相補償を行う。
Subsequently, the
具体的には、まず、レンジ補償後のヒット幅Ｈ０のレンジヒストリＳ１（ｒ，ｈ）から、最終的に画像再生に用いるヒット幅ＨのレンジヒストリＳ２（ｒ，ｈ）を、以下の式（４）のように得る。 Specifically, first, from the range history S1 (r, h) of the hit width H0 after the range compensation, the range history S2 (r, h) of the hit width H finally used for image reproduction is expressed by the following formula ( It is obtained as in 4).
ヒット幅ＨのレンジヒストリＳ２（ｒ，ｈ）を格納する配列のイメージは、図６（ｅ）に示した通りである。
また、位相補償回路２２は、以下の式（５）のように、レンジヒストリＳ２（ｒ，ｈ）の位相補償を行い、位相補償後のレンジヒストリＳ３（ｒ，ｈ）（ｈ＝−Ｈ／２、−Ｈ／２＋１、・・・、０、・・・、Ｈ／２−２、Ｈ／２−１；ｒ＝０、１、・・・、Ｒ−１）を得る。
An image of an array for storing the range history S2 (r, h) of the hit width H is as shown in FIG.
Further, the
これにより、クロスレンジ軸方向の画像ぼけの要因である不要位相変化φ２（ｈ）を補償することができる。
位相補償後のレンジヒストリＳ３（ｒ，ｈ）は、クロスレンジ圧縮器１４に入力される。
Thereby, it is possible to compensate for the unnecessary phase change φ2 (h) that is a cause of image blur in the cross range axis direction.
The range history S3 (r, h) after phase compensation is input to the
最後に、レーダ画像化器２（図３）内のクロスレンジ圧縮器１４は、公知のクロスレンジ圧縮処理、すなわちヒットｈ方向の離散フーリエ変換を行うことにより、レンジヒストリＳ３（ｒ，ｈ）のクロスレンジ方向の分解能を向上させたレーダ画像Ｕ（ｒ，ｆ）（ｆ＝０、１、・・・、Ｈ／２−１；ｒ＝０、１、・・・、Ｒ−１）を、以下の式（６）のように取得する。
Finally, the
なお、切出処理のヒット幅Ｈ０が最終的な画像再生用のヒット幅Ｈと一致する場合は、通常のＰＧＡ法を用いた場合の位相補償と一致するので、この発明の実施の形態１による上記方法は、ＰＧＡ法の自然な拡張方法の１つと位置付けられる。 Note that, when the hit width H0 of the clipping process matches the final image reproduction hit width H, it matches the phase compensation in the case of using the normal PGA method, and therefore, according to the first embodiment of the present invention. The above method is positioned as one of natural extension methods of the PGA method.
以上のように、この発明の実施の形態１（図１〜図６）に係る画像レーダ装置は、目標Ｔのレンジプロフィールを取得するレーダ観測器１と、レンジプロフィールのレンジヒストリＳに基づいてレーダ画像Ｕ（ｒ，ｆ）を生成するレーダ画像化器２と、により構成されている。
As described above, the image radar apparatus according to Embodiment 1 (FIGS. 1 to 6) of the present invention is based on the
レーダ観測器１（図２）は、高周波信号を発生する送信機３と、高周波信号を目標Ｔに照射するとともに目標Ｔで散乱された高周波信号の一部を受信する送受信アンテナ５と、送受信アンテナ５を介して目標Ｔからの受信信号を取り込む受信機６と、送信機３からの高周波信号を送受信アンテナ５に送出するとともに送受信アンテナ５からの受信信号を受信機６に入力する送受切換器４と、受信信号のレンジ方向の分解能をレンジ圧縮により向上させてレンジプロフィールを取得するレンジ圧縮器７と、を備えている。
The radar observation device 1 (FIG. 2) includes a
レーダ観測器１は、レンジプロフィールを取得する一連の処理を、目標Ｔとレーダ観測器１との間の相対位置関係を変えながら繰り返し実行してレンジヒストリＳを取得する。
レーダ画像化器２（図３）は、レンジヒストリＳをクロスレンジ圧縮することにより、目標Ｔ上の反射強度分布をレンジドップラー分布として表現したレーダ画像Ｕ（ｒ，ｆ）を得る。
The
The radar imager 2 (FIG. 3) obtains a radar image U (r, f) expressing the reflection intensity distribution on the target T as a range Doppler distribution by cross-range compressing the range history S.
レーダ画像化器２は、最終的にクロスレンジ圧縮に用いるヒット幅Ｈ以上のヒット幅Ｈ０のレンジヒストリＳ０（ｒ，ｈ）に対して、目標Ｔとレーダ観測器１との間の相対運動で発生する観測中の各反射点のレンジ分解能セルの移動を推定して補償するレンジ補償器１２を備えている。
The
また、レーダ画像化器２は、レンジ補償器１２によるレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）上の代表的な複数の反射点の位相変化を調べ、レーダ画像のドップラー周波数方向のぼけの原因となるヒットｈに対する２次以上の不要位相変化φ２（ｈ）の、最終的にクロスレンジ圧縮に用いるヒット幅Ｈの範囲での値を推定する不要位相変化推定器２１（図４、図５）を備えている。
Further, the
さらに、レーダ画像化器２は、不要位相変化推定器２１で得られた不要位相変化φ２（ｈ）に基づき、レンジ補償後のレンジヒストリから、最終的にクロスレンジ圧縮に用いるヒット幅分だけ切出したレンジヒストリに含まれる不要位相変化成分を補償する位相補償回路２２（図４）を備えている。
Further, the
不要位相変化推定器２１（図５）は、何らかの方法で一度推定されたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）について、不要位相変化φ１（ｈ）のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅Ｈの不要位相変化φ２（ｈ）のみを抽出するヒット範囲切出器１１１を含む。
The unnecessary phase change estimator 21 (FIG. 5) uses the unnecessary phase change φ1 (h) for the unnecessary phase change φ1 (h) included in the range history S1 (r, h) after the range compensation once estimated by some method. This includes a
これにより、一般的なＰＧＡ法で位相補償する場合に発生する補償誤差、すなわち、不要位相変化φ１（ｈ）の推定誤差を低減することができる。
すなわち、ＰＧＡ法を用いて不要位相変化φ１（ｈ）を推定する際に、受信信号の条件によっては発生する可能性がある「全ヒットにわたる推定誤差」、または「或るヒット付近のみで増大する局所的な推定誤差」の影響を低減して、位相補償精度を向上させることができる。
As a result, it is possible to reduce a compensation error that occurs when phase compensation is performed by a general PGA method, that is, an estimation error of the unnecessary phase change φ1 (h).
That is, when the unnecessary phase change φ1 (h) is estimated using the PGA method, “estimation error over all hits” that may occur depending on the condition of the received signal, or “increases only in the vicinity of a certain hit”. It is possible to improve the phase compensation accuracy by reducing the influence of “local estimation error”.
特に、ＰＧＡ法において、受信信号の性質（たとえば、注目するレンジセルｒの反射点の数や位置関係、不要位相変化φ１（ｈ）の次数など）によって、ヒットｈの中央付近では不要位相変化の推定精度が比較的高いが、端部付近では不要位相変化φ１（ｈ）の推定精度が悪くなるような状況を想定して、あらかじめ画像再生に用いるヒット幅Ｈよりも広いヒット幅Ｈ０のレンジヒストリＳ０（ｒ，ｈ）のレンジ補償および不要位相変化推定を施したうえで、ここから誤差が発生している可能性がある不要な端部の推定結果を捨てて、位相補償に本来必要な不要位相変化φ２（ｈ）を最終的に得ることにより、仮定した不要位相変化φ１（ｈ）の推定誤差の影響を低減することができる。
この結果、再生されたレーダ画像Ｕ（ｒ，ｆ）の分解能を向上させることができる。
In particular, in the PGA method, an unnecessary phase change is estimated in the vicinity of the center of the hit h depending on the nature of the received signal (for example, the number and positional relationship of the reflection points of the target range cell r, the order of the unnecessary phase change φ1 (h), etc.). A range history S0 having a hit width H0 wider than the hit width H used in advance for image reproduction is assumed assuming that the estimation accuracy of the unnecessary phase change φ1 (h) is poor near the end, although the accuracy is relatively high. After performing (r, h) range compensation and unnecessary phase change estimation, discard the unnecessary end estimation result that may have caused an error from here, and unnecessary phase originally required for phase compensation By finally obtaining the change φ2 (h), it is possible to reduce the influence of the estimated error of the assumed unnecessary phase change φ1 (h).
As a result, the resolution of the reproduced radar image U (r, f) can be improved.
実施の形態２．
なお、上記実施の形態１（図５）では、不要位相変化推定器２１内にヒット範囲切出器１１１を設けたが、図７のように、ヒット範囲切出器１１１に代えて、位相フィッティング器１２１を設けてもよい。
図７はこの発明の実施の形態２による不要位相変化推定器２１Ａの構成を示すブロック図である。また、この発明の実施の形態２に係る画像レーダ装置の他の構成は、図１〜図４に示した通りである。
In the first embodiment (FIG. 5), the
FIG. 7 is a block diagram showing a configuration of an unnecessary
図７において、不要位相変化推定器２１Ａは、ＰＧＡ推定器３１の後段に、前述（図５）のヒット範囲切出器１１１に代えて、位相フィッティング器１２１を備えている。
位相フィッティング器１２１は、入力した不要位相変化φ１（ｈ）に対し、あらかじめ設定した次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化φ３（ｈ）として抽出し、位相補償回路２２に入力する。
In FIG. 7, the unnecessary
The phase
次に、図１〜図４、図６および図７を参照しながら、この発明の実施の形態１による動作について説明する。
まず、レーダ観測器１（図１、図２）は、前述と同様に、レンジヒストリＳの取得処理を行い、レンジヒストリ切出器１１（図３）は、運動補償および画像再生に用いるレンジヒストリＳ０（ｒ，ｈ）の切出処理を行う。
Next, the operation according to the first embodiment of the present invention will be described with reference to FIGS. 1 to 4, 6 and 7.
First, the radar observation device 1 (FIGS. 1 and 2) performs the process of acquiring the range history S as described above, and the range history extractor 11 (FIG. 3) uses the range history used for motion compensation and image reproduction. The extraction process of S0 (r, h) is performed.
また、前述と同様に、レンジ補償器１２（図３、図６）は、切出したレンジヒストリＳ０（ｒ，ｈ）のレンジ補償を行い、レンジヒストリＳ１（ｒ，ｈ）を取得し、不要位相変化推定器２１Ａ（図７）内のＰＧＡ推定器３１は、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に基づき、不要位相変化φ１（ｈ）を取得する処理を行う。
In the same manner as described above, the range compensator 12 (FIGS. 3 and 6) performs range compensation on the extracted range history S0 (r, h), acquires the range history S1 (r, h), and obtains an unnecessary phase. The
以下、この発明の実施の形態２による不要位相変化推定器２１Ａ（図７）内の位相フィッティング器１２１は、最小２乗法を適用して、不要位相変化φ３（ｈ）を抽出し、位相補償器１３（図４）内の位相補償回路２２に入力する。
Hereinafter, the phase
最後に、前述と同様に、位相補償回路２２は、不要位相変化推定器２１Ａからの不要位相変化φ３（ｈ）に基づき位相補償後のレンジヒストリＳ３（ｒ，ｈ）を取得し、クロスレンジ圧縮器１４は、レンジヒストリＳ３（ｒ，ｈ）をクロスレンジ圧縮して、レーダ画像Ｕ（ｒ，ｆ）の取得処理を行う。
Finally, as described above, the
ここでは、前述と同様に、ＰＧＡ法を用いて不要位相変化φ１（ｈ）を推定した際に、受信信号の性質によっては、不要位相変化φ１（ｈ）の推定誤差が増大する可能性があることを想定し、この影響を低減することを目的としている。
また、特に、同じレンジｒに存在する複数の反射点間の干渉の影響や、雑音の影響で、ヒットｈごとに細かく変動する不要位相変化φ１（ｈ）の推定誤差が存在するような状況に対処することを目的としている。
Here, as described above, when the unnecessary phase change φ1 (h) is estimated using the PGA method, the estimation error of the unnecessary phase change φ1 (h) may increase depending on the nature of the received signal. This is intended to reduce this effect.
In particular, there is an estimation error of the unnecessary phase change φ1 (h) that varies finely for each hit h due to the influence of interference between a plurality of reflection points existing in the same range r or the influence of noise. The purpose is to deal with it.
上記推定誤差を低減するために、ＰＧＡ推定器３１の後段に設けられた位相フィッティング器１２１は、ＰＧＡ推定器３１で得られた不要位相変化φ１（ｈ）に対し、あらかじめ設定した次数の一般的な最小２乗法を適用し、２乗誤差を最小とする不要位相変化φ３（ｈ）を得る。
In order to reduce the estimation error, the phase
このように、最小２乗法を適用することにより、最小２乗法の低域通過フィルタとしての性質を利用して、ヒットｈごとに細かく変動する不要位相変化φ１（ｈ）の推定誤差を低減することができる。 In this way, by applying the least square method, the estimation error of the unnecessary phase change φ1 (h) that varies finely for each hit h can be reduced by utilizing the property of the least square method as a low-pass filter. Can do.
なお、前述と同様に、レンジヒストリ切出器１１（図３）で切出すヒット幅Ｈ０を、最終的に画像再生に用いるヒット幅Ｈと一致させることにより、不要位相変化φ１（ｈ）、φ３（ｈ）のヒット幅がいずれもＨとなる。
また、切出しヒット幅Ｈ＝Ｈ０の条件により、Ｓ１（ｒ，ｈ）と完全に一致する場合、位相補償器１３（図４）内の位相補償回路２２において、最終的に画像再生に用いるヒット幅Ｈのレンジ補償後のレンジヒストリＳ２（ｒ，ｈ）を、不要位相変化φ３（ｈ）を用いて位相補償する処理についても、前述と同様である。
Similarly to the above, unnecessary phase changes φ1 (h), φ3 are obtained by making the hit width H0 cut out by the range history extractor 11 (FIG. 3) finally coincide with the hit width H used for image reproduction. Both hit widths in (h) are H.
When the cut hit width H = H0 completely matches S1 (r, h), the
以上のように、この発明の実施の形態２（図７）による不要位相変化推定器２１Ａは、一度推定されたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化φ３（ｈ）の推定値として出力する位相フィッティング器１２１を備えている。
As described above, the unnecessary
これにより、ＰＧＡ法を用いて不要位相変化を推定した際に、受信信号の性質によっては増大する可能性がある不要位相変化φ１（ｈ）の推定誤差のうちの、特に同じレンジｒに存在する複数の反射点間の干渉の影響や雑音の影響で、ヒットｈごとに細かく変動する不要位相変化φ１（ｈ）の推定誤差を、最小２乗法の低域通過フィルタとしての性質を利用して低減することができるので、不要位相変化の推定精度を向上させることができる。 As a result, when the unnecessary phase change is estimated using the PGA method, the estimation error of the unnecessary phase change φ1 (h), which may increase depending on the properties of the received signal, exists particularly in the same range r. Reduces the estimation error of the unnecessary phase change φ1 (h) that fluctuates finely for each hit h due to the influence of interference between multiple reflection points and the noise, utilizing the characteristics of a low-pass filter of the least square method Therefore, the estimation accuracy of the unnecessary phase change can be improved.
実施の形態３．
なお、上記実施の形態１、２（図５、図７）では、不要位相変化推定器２１、２１Ａ内にヒット範囲切出器１１１、位相フィッティング器１２１の一方を設けたが、図８のように、不要位相変化推定器２１Ｂ内にヒット範囲切出器１１１および位相フィッティング器１２１の両方を設けてもよい。
In the first and second embodiments (FIGS. 5 and 7), one of the
図８はこの発明の実施の形態３による不要位相変化推定器２１Ｂの構成を示すブロック図である。また、この発明の実施の形態３に係る画像レーダ装置の他の構成は、図１〜図４に示した通りである。
図８において、不要位相変化推定器２１Ｂは、ＰＧＡ推定器３１の後段に、位相フィッティング器１２１を備え、さらに、位相フィッティング器１２１の後段に、ヒット範囲切出器１１１を備えている。
FIG. 8 is a block diagram showing a configuration of an unnecessary phase change estimator 21B according to
In FIG. 8, the unnecessary
次に、図１〜図４、図６および図８を参照しながら、この発明の実施の形態３による動作について説明する。
まず、レーダ観測器１（図１、図２）は、前述と同様に、レンジヒストリＳの取得処理を行い、レンジヒストリ切出器１１（図３）は、運動補償および画像再生に用いるレンジヒストリＳ０（ｒ，ｈ）の切出処理を行う。
Next, the operation according to the third embodiment of the present invention will be described with reference to FIGS. 1 to 4, 6 and 8.
First, the radar observation device 1 (FIGS. 1 and 2) performs the process of acquiring the range history S as described above, and the range history extractor 11 (FIG. 3) uses the range history used for motion compensation and image reproduction. The extraction process of S0 (r, h) is performed.
また、前述と同様に、レンジ補償器１２（図３、図６）は、切出したレンジヒストリＳ０（ｒ，ｈ）のレンジ補償を行い、レンジヒストリＳ１（ｒ，ｈ）を取得し、不要位相変化推定器２１Ｂ（図８）内のＰＧＡ推定器３１は、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に基づき、不要位相変化φ１（ｈ）を取得する処理を行う。
In the same manner as described above, the range compensator 12 (FIGS. 3 and 6) performs range compensation on the extracted range history S0 (r, h), acquires the range history S1 (r, h), and obtains an unnecessary phase. The
以下、この発明の実施の形態３による不要位相変化推定器２１Ｂ（図８）内の位相フィッティング器１２１は、最小２乗法を適用して、不要位相変化φ３（ｈ）を抽出し、ヒット範囲切出器１１１に入力する。
ヒット範囲切出器１１１は、不要位相変化φ３（ｈ）から不要位相変化φ２（ｈ）を切出して位相補償器１３（図４）内の位相補償回路２２に入力する。
Hereinafter, the phase
The
位相補償回路２２は、前述と同様に、不要位相変化推定器２１Ｂからの不要位相変化φ２（ｈ）に基づき、最終的に画像再生に用いるヒット幅Ｈのレンジ補償後のレンジヒストリＳ２（ｒ，ｈ）を補償して、位相補償後のレンジヒストリＳ３（ｒ，ｈ）を取得する。
最後に、クロスレンジ圧縮器１４は、レンジヒストリＳ３（ｒ，ｈ）をクロスレンジ圧縮して、レーダ画像Ｕ（ｒ，ｆ）の取得処理を行う。
As described above, the
Finally, the
この発明の実施の形態３においても、レーダ画像化器２（図３）内のレンジヒストリ切出器１１により切出すヒット幅Ｈ０は、最終的に画像再生に用いるヒット幅Ｈ以上に設定されており、レンジヒストリＳ１（ｒ，ｈ）および不要位相変化φ１（ｈ）のヒット幅はＨ０であり、不要位相変化φ２（ｈ）のヒット幅はＨである。
Also in the third embodiment of the present invention, the hit width H0 extracted by the
また、前述と同様に、ＰＧＡ法を用いて不要位相変化を推定した際に、受信信号の性質によっては、不要位相変化の推定誤差が増大する可能性があることを想定し、この影響を低減することを目的としている。 Also, as described above, when estimating the unnecessary phase change using the PGA method, it is assumed that the estimation error of the unnecessary phase change may increase depending on the nature of the received signal, and this effect is reduced. The purpose is to do.
さらに、特に、同じレンジｒに存在する複数の反射点間の干渉の影響や、雑音の影響で、ヒットｈごとに細かく変動する不要位相変化の推定誤差が存在するような場合や、ヒットｈの中央付近では不要位相変化の推定精度が比較的高いが、端部付近では不要位相変化の推定精度が低いような場合が混在して発生するような状況に対処することを目的としている。 Furthermore, in particular, when there is an estimation error of an unnecessary phase change that varies finely for each hit h due to the influence of interference between a plurality of reflection points existing in the same range r or the influence of noise, The object of the present invention is to deal with a situation in which the estimation accuracy of the unnecessary phase change is relatively high near the center but the estimation accuracy of the unnecessary phase change is low near the end.
ＰＧＡ推定器３１の後段の位相フィッティング器１２１は、ＰＧＡ推定器３１で得られた不要位相変化φ１（ｈ）に対して、あらかじめ設定した次数の一般的な最小２乗法を適用して、２乗誤差を最小とする不要位相変化φ３（ｈ）（ｈ＝−Ｈ０／２、−Ｈ０／２＋１、・・・、０、・・・、Ｈ０／２−２、Ｈ０／２−１）を得る。
The phase
不要位相変化φ３（ｈ）においては、最小２乗法の低域通過フィルタとしての性質により、ヒットｈごとに細かく変動する不要位相変化の推定誤差が低減されているものと期待されるが、端部付近で増大する不要位相変化φ１（ｈ）の推定誤差の影響が含まれている可能性がある。 In the unnecessary phase change φ3 (h), it is expected that the estimation error of the unnecessary phase change that varies finely for each hit h is reduced due to the property as a low-pass filter of the least square method. There is a possibility that the influence of the estimation error of the unnecessary phase change φ1 (h) that increases in the vicinity is included.
そこで、ヒット範囲切出器１１１は、不要位相変化φ３（ｈ）のうちの、推定誤差が増大している可能性が高く、かつ不要な、端部付近のデータを捨てて、最終的に画像再生に用いるヒット幅Ｈの不要位相変化φ２（ｈ）を取得する。
この処理は、前述の式（３）のように、不要位相変化φ１（ｈ）に基づき不要位相変化φ２（ｈ）を取得する処理と同様である。
これにより、端部付近で増大する不要位相変化の推定誤差を低減することができる。
Therefore, the
This process is the same as the process of obtaining the unnecessary phase change φ2 (h) based on the unnecessary phase change φ1 (h) as in the above-described equation (3).
Thereby, the estimation error of the unnecessary phase change that increases near the end can be reduced.
以上のように、この発明の実施の形態３（図８）による不要位相変化推定器２１Ｂは、一度推定されたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化φ３（ｈ）の推定値として出力する位相フィッティング器１２１を備えている。
As described above, the unnecessary phase change estimator 21B according to the third embodiment (FIG. 8) of the present invention uses the unnecessary phase change φ1 (h) included in the range history S1 (r, h) after the range compensation once estimated. ), A phase
また、不要位相変化推定器２１Ｂは、位相フィッティング器１２１からの不要位相変化φ３（ｈ）のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅Ｈの不要位相変化φ２（ｈ）を抽出するヒット範囲切出器１１１を備えている。
Further, the unnecessary phase change estimator 21B discards the estimated values near the hit ends of the unnecessary phase change φ3 (h) from the phase
これにより、前述と同様の効果を奏するとともに、同じレンジｒに存在する複数の反射点間の干渉の影響や雑音の影響で、ヒットｈごとに細かく変動する不要位相変化の推定誤差が存在するような場合や、ヒットｈの中央付近では不要位相変化の推定精度が比較的高いが、端部付近では不要位相変化の推定精度が低いような場合が混在して発生するような状況でも、不要位相変化φ２（ｈ）の推定誤差を低減して、高精度に位相補償を行うことができ、再生されたレーダ画像Ｕ（ｒ，ｆ）のぼけの発生を低減することができる。 As a result, the same effect as described above can be obtained, and there is an estimation error of unnecessary phase change that varies finely for each hit h due to the influence of interference between a plurality of reflection points existing in the same range r or the influence of noise. Even in situations where the estimation accuracy of the unnecessary phase change is relatively high near the center of the hit h, and the estimation accuracy of the unnecessary phase change is low near the edge, The estimation error of the change φ2 (h) can be reduced, phase compensation can be performed with high accuracy, and the occurrence of blurring of the reproduced radar image U (r, f) can be reduced.
実施の形態４．
なお、上記実施の形態３（図８）では、ヒット範囲切出器１１１の前段に位相フィッティング器１２１を設けたが、図９のように、ヒット範囲切出器１１１の後段に位相フィッティング器１２１を設けてもよい。
In the third embodiment (FIG. 8), the phase
図９はこの発明の実施の形態４による不要位相変化推定器２１Ｃの構成を示すブロック図である。また、この発明の実施の形態４に係る画像レーダ装置の他の構成は、図１〜図４に示した通りである。
図９において、不要位相変化推定器２１Ｃは、ＰＧＡ推定器３１の後段に、ヒット範囲切出器１１１を備え、さらに、ヒット範囲切出器１１１の後段に、位相フィッティング器１２１を備えている。
FIG. 9 is a block diagram showing a configuration of an unnecessary
In FIG. 9, the unnecessary
次に、図１〜図４、図６および図９を参照しながら、この発明の実施の形態４による動作について説明する。
この場合、不要位相変化推定器２１Ｃ内のヒット範囲切出器１１１および位相フィッティング器１２１の処理順序が入れ替わったのみであり、他の処理は前述と同様である。
Next, the operation according to the fourth embodiment of the present invention will be described with reference to FIGS. 1 to 4, 6 and 9.
In this case, only the processing order of the
すなわち、レーダ観測器１（図２）によるレンジヒストリＳの取得処理、レーダ画像化器２内のレンジヒストリ切出器１１（図３）によるレンジヒストリＳ０（ｒ，ｈ）の切出処理、レンジ補償器１２（図３、図６）によるレンジヒストリＳ１（ｒ，ｈ）の取得処理、不要位相変化推定器２１Ｃ内のＰＧＡ推定器３１による不要位相変化φ１（ｈ）の取得処理、位相補償回路２２（図４）によるレンジヒストリＳ３（ｒ，ｈ）の取得処理、および、クロスレンジ圧縮器１４によるレーダ画像Ｕ（ｒ，ｆ）の取得処理は、前述と同様である。
That is, the range history S acquisition process by the radar observer 1 (FIG. 2), the range history S0 (r, h) extraction process by the range history extractor 11 (FIG. 3) in the
また、この発明の実施の形態４においても、レンジヒストリ切出器１１で切出すヒット幅Ｈ０は、最終的に画像再生に用いるヒット幅Ｈ以上に設定されており、レンジヒストリＳ１（ｒ，ｈ）および不要位相変化φ１（ｈ）のヒット幅はＨ０であり、不要位相変化φ２（ｈ）のヒット幅はＨである。
Also in
また、前述と同様に、ＰＧＡ法を用いて不要位相変化を推定した際に、受信信号の性質によっては不要位相変化の推定誤差が増大する可能性があることを想定し、この影響を低減することを目的としている。 Further, as described above, when an unnecessary phase change is estimated using the PGA method, it is assumed that an estimation error of the unnecessary phase change may increase depending on the characteristics of the received signal, and this influence is reduced. The purpose is that.
さらに、同じレンジに存在する複数の反射点間の干渉の影響や、雑音の影響で、ヒットごとに細かく変動する不要位相変化の推定誤差が存在するような場合や、ヒットの中央付近では不要位相変化の推定精度が比較的高いが、端部付近では不要位相変化の推定精度が低いような場合が混在して発生するような状況に対処することを目的としている。 In addition, there is an unnecessary phase change estimation error that varies finely from hit to hit due to the influence of interference between multiple reflection points in the same range, or the influence of noise, or an unnecessary phase near the center of the hit. The object is to cope with a situation in which the estimation accuracy of the change is relatively high but the estimation accuracy of the unnecessary phase change is low near the end portion.
不要位相変化推定器２１Ｃにおいて、ＰＧＡ推定器３１の後段のヒット範囲切出器１１１は、前述と同様に、不要位相変化φ１（ｈ）から、最終的に画像再生に用いるヒット幅Ｈの不要位相変化φ２（ｈ）（ｈ＝−Ｈ／２、−Ｈ／２＋１、・・・、０、・・・、Ｈ／２−２、Ｈ／２−１）を抽出する。
なお、ヒット幅Ｈは、前述の実施の形態３（図８）における不要位相変化φ３（ｈ）のヒット幅とは異なる。
In the unnecessary
The hit width H is different from the hit width of the unnecessary phase change φ3 (h) in the above-described third embodiment (FIG. 8).
不要位相変化φ２（ｈ）においては、端部付近で増大する不要位相変化の推定誤差が低減されているものと期待されるが、ヒットごとに細かく変動する不要位相変化の推定誤差は残存している可能性が高い。 In the unnecessary phase change φ2 (h), it is expected that the estimation error of the unnecessary phase change that increases in the vicinity of the end is reduced, but the estimation error of the unnecessary phase change that varies finely for each hit remains. There is a high possibility.
そこで、ヒット範囲切出器１１１の後段の位相フィッティング器１２１は、不要位相変化φ２（ｈ）に対して、あらかじめ設定した次数の最小２乗法を適用して、２乗誤差を最小とする不要位相変化φ３（ｈ）（ｈ＝−Ｈ／２、−Ｈ／２＋１、・・・、０、・・・、Ｈ／２−２、Ｈ／２−１）を取得する。
Therefore, the phase
不要位相変化φ３（ｈ）においては、不要位相変化φ２（ｈ）に対して、最小２乗法の低域通過フィルタとしての性質を適用することにより、ヒットｈごとに細かく変動する不要位相変化の推定誤差を低減できているものと期待される。 In the unnecessary phase change φ3 (h), by applying the property as a low-pass filter of the least square method to the unnecessary phase change φ2 (h), estimation of the unnecessary phase change that varies finely for each hit h. It is expected that the error can be reduced.
以上のように、この発明の実施の形態４（図９）による不要位相変化推定器２１Ｃは、入力された不要位相変化φ１（ｈ）に対して、不要位相変化φ１（ｈ）のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅Ｈの不要位相変化φ２（ｈ）を抽出するヒット範囲切出器１１１と、入力された不要位相変化φ２（ｈ）に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化φ３（ｈ）の推定値として出力する位相フィッティング器１２１と、を備えている。
As described above, the unnecessary
すなわち、ヒット範囲切出器１１１は、一度推定されたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）について、不要位相変化φ１（ｈ）のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅Ｈの不要位相変化φ２（ｈ）を抽出する。
位相フィッティング器１２１は、入力されたヒット範囲切出器１１１からの不要位相変化φ２（ｈ）に対して、設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化φ３（ｈ）の推定値として出力する。
That is, the
The phase
これにより、前述の実施の形態１、２と同様の効果を奏するとともに、実施の形態３と同様に、同じレンジに存在する複数の反射点間の干渉の影響や、雑音の影響で、ヒットごとに細かく変動する不要位相変化の推定誤差が存在するような場合や、ヒットの中央付近では不要位相変化の推定精度が比較的高いが、端部付近では不要位相変化の推定精度が低いような場合が混在して発生するような状況でも、不要位相変化の推定誤差を低減して、高精度に位相補償を行うことができ、換言すれば、再生画像のぼけの発生を低減することができる。 As a result, the same effects as those of the first and second embodiments described above can be obtained, and, similarly to the third embodiment, each hit may be affected by the influence of interference between a plurality of reflection points existing in the same range or the influence of noise. If there is an estimation error of the unnecessary phase change that fluctuates slightly, or if the estimation accuracy of the unnecessary phase change is relatively high near the center of the hit, but the estimation accuracy of the unnecessary phase change is low near the edge Even in a situation where both occur, the estimation error of the unnecessary phase change can be reduced and the phase compensation can be performed with high accuracy. In other words, the occurrence of blurring of the reproduced image can be reduced.
ところで、前述の実施の形態３（図８）とこの発明の実施の形態４（図９）とを比較すると、構成的には、位相フィッティング器１２１とヒット範囲切出器１１１の順序が入れ替わっている点のみが相違しているが、この相違点による効果の違いについて述べる。
By the way, when the above-described third embodiment (FIG. 8) and the fourth embodiment (FIG. 9) of the present invention are compared, the order of the phase
前述（図８）の構成では、ＰＧＡ推定器３１の直後に位相フィッティング器１２１が配置されているので、位相フィッティング器１２１が出力する不要位相変化φ３（ｈ）においては、端部付近の大きな推定誤差の影響が広いヒット範囲に伝播する可能性がある。これは特に、最小２乗法の次数を低くした場合に顕著になる。
つまり、図８の場合、位相フィッティング器１２１の後段のヒット範囲切出器１１１で抽出するヒット幅にも上記誤差の影響が漏れ込んでくる可能性がある。
In the configuration of FIG. 8 (FIG. 8), since the phase
That is, in the case of FIG. 8, there is a possibility that the influence of the error leaks into the hit width extracted by the
これに対し、図９の構成によれば、ヒット範囲切出器１１１により先に端部のデータを捨てた後に、位相フィッティング器１２１により最小２乗法を適用するので、図８の構成による上記誤差伝播の発生を低減することができる。
On the other hand, according to the configuration of FIG. 9, the least square method is applied by the phase
実施の形態５．
なお、上記実施の形態３、４（図８、図９）では、不要位相変化推定器２１Ｂ、２１Ｃ内に単一のヒット範囲切出器１１１を設けたが、図１０のように、不要位相変化推定器２１Ｄ内において、位相フィッティング器１２１の前後に、第１ヒット範囲切出器１５１および第２ヒット範囲切出器１５２を設けてもよい。
Embodiment 5 FIG.
In the third and fourth embodiments (FIGS. 8 and 9), the single
図１０はこの発明の実施の形態５による不要位相変化推定器２１Ｄの構成を示すブロック図であり、図１１は不要位相変化推定器２１Ｄの処理内容を示す説明図である。また、この発明の実施の形態５に係る画像レーダ装置の他の構成は、図１〜図４に示した通りである。
FIG. 10 is a block diagram showing a configuration of an unnecessary
図１０において、不要位相変化推定器２１Ｄは、ＰＧＡ推定器３１の後段に、第１ヒット範囲切出器１５１を備え、さらに、ヒット範囲切出器１１１の後段に、位相フィッティング器１２１を介して、第２ヒット範囲切出器１５２を備えている。
In FIG. 10, the unnecessary
第１ヒット範囲切出器１５１は、運動補償に用いるレンジヒストリと同じヒット幅Ｈ０のデータ列（不要位相変化φ１（ｈ））から、ヒット幅Ｈ０以下で、かつ、最終的に画像再生に用いるヒット幅Ｈ以上の中間的なヒット幅Ｈ１のデータ列（不要位相変化φ２１（ｈ））を抽出する。
The first
第２ヒット範囲切出器１５２は、位相フィッティング器１２１を介した中間的なヒット幅Ｈ１のデータ列（不要位相変化φ３（ｈ））から、最終的に画像再生に用いるヒット幅Ｈのデータ列（不要位相変化φ２２（ｈ））を抽出する。
The second
次に、図１〜図４、図１０および図１１を参照しながら、この発明の実施の形態５による動作について説明する。
この場合、不要位相変化推定器２１Ｄ内において、位相フィッティング器１２１の前後に第１ヒット範囲切出器１５１および第１ヒット範囲切出器１５２が設けられた点を除けば、他の処理は前述と同様である。
Next, the operation according to the fifth embodiment of the present invention will be described with reference to FIGS. 1 to 4, 10 and 11. FIG.
In this case, in the unnecessary
すなわち、レーダ観測器１（図２）によるレンジヒストリＳの取得処理、レーダ画像化器２内のレンジヒストリ切出器１１（図３）によるレンジヒストリＳ０（ｒ，ｈ）の切出処理、レンジ補償器１２（図３、図１１）によるレンジヒストリＳ１（ｒ，ｈ）の取得処理、不要位相変化推定器２１Ｄ内のＰＧＡ推定器３１による不要位相変化φ１（ｈ）の取得処理、位相補償回路２２（図４）によるレンジヒストリＳ３（ｒ，ｈ）の取得処理、および、クロスレンジ圧縮器１４によるレーダ画像Ｕ（ｒ，ｆ）の取得処理は、前述と同様である。
That is, the range history S acquisition process by the radar observer 1 (FIG. 2), the range history S0 (r, h) extraction process by the range history extractor 11 (FIG. 3) in the
また、この場合も、レンジヒストリ切出器１１で切出すヒット幅Ｈ０は、最終的に画像再生に用いるヒット幅Ｈ以上に設定されており、レンジヒストリＳ１（ｒ，ｈ）および不要位相変化φ１（ｈ）のヒット幅はＨ０であり、不要位相変化φ２（ｈ）のヒット幅はＨである。
Also in this case, the hit width H0 cut out by the
また、前述と同様に、ＰＧＡ法を用いて不要位相変化を推定した際に、受信信号の性質によっては不要位相変化の推定誤差が増大する可能性があることを想定し、この影響を低減することを目的としている。 Further, as described above, when an unnecessary phase change is estimated using the PGA method, it is assumed that an estimation error of the unnecessary phase change may increase depending on the characteristics of the received signal, and this influence is reduced. The purpose is that.
さらに、同じレンジに存在する複数の反射点間の干渉の影響や、雑音の影響で、ヒットごとに細かく変動する不要位相変化の推定誤差が存在するような場合や、ヒットの中央付近では不要位相変化の推定精度が比較的高いが、端部付近では不要位相変化の推定精度が低いような場合が混在して発生するような状況に対処することを目的としている。 In addition, there is an unnecessary phase change estimation error that varies finely from hit to hit due to the influence of interference between multiple reflection points in the same range, or the influence of noise, or an unnecessary phase near the center of the hit. The object is to cope with a situation in which the estimation accuracy of the change is relatively high but the estimation accuracy of the unnecessary phase change is low near the end portion.
このため、まず、ＰＧＡ推定器３１の後段の第１ヒット範囲切出器１５１は、不要位相変化φ１（ｈ）から、以下の式（７）のように、中間的なヒット幅Ｈ１（Ｈ≦Ｈ１≦Ｈ０）の不要位相変化φ２１（ｈ）を抽出する。
Therefore, first, the first
なお、レンジヒストリＳ１（ｒ，ｈ）および不要位相変化φ１（ｈ）、φ２１（ｈ）のイメージは、図１１（ａ）〜図１１（ｃ）に示した通りである。
これにより、不要位相変化φ２１（ｈ）においては、端部付近で増大する不要位相変化の推定誤差が低減されているものと期待される。ただし、ヒットごとに細かく変動する不要位相変化の推定誤差は残存している可能性が高い。
The images of the range history S1 (r, h) and unnecessary phase changes φ1 (h), φ21 (h) are as shown in FIGS. 11 (a) to 11 (c).
Thereby, in the unnecessary phase change φ21 (h), it is expected that the estimation error of the unnecessary phase change that increases near the end portion is reduced. However, there is a high possibility that an estimation error of an unnecessary phase change that varies finely for each hit remains.
そこで、第１ヒット範囲切出器１５１の後段の位相フィッティング器１２１は、不要位相変化φ２１（ｈ）に対して、あらかじめ設定した次数の最小２乗法を適用し、２乗誤差を最小とする不要位相変化φ３（ｈ）（ｈ＝−Ｈ１／２、−Ｈ１／２＋１、・・・、０、・・・、Ｈ１／２−２、Ｈ１／２−１）を取得する。
Therefore, the phase
不要位相変化φ３（ｈ）のイメージは、図１１（ｄ）に示した通りである。
不要位相変化φ３（ｈ）においては、不要位相変化φ２１（ｈ）に対して、さらに、最小２乗法の低域通過フィルタとしての性質を適用することにより、ヒットｈごとに細かく変動する不要位相変化の推定誤差を低減できているものと期待される。
An image of the unnecessary phase change φ3 (h) is as shown in FIG.
In the unnecessary phase change φ3 (h), the property as a low-pass filter of the least square method is further applied to the unnecessary phase change φ21 (h), so that the unnecessary phase change finely varies for each hit h. It is expected that the estimation error can be reduced.
ただし、最小２乗法を適用する際に、たまたま不要位相変化φ２１（ｈ）のヒット端部付近で、受信信号の性質の影響で誤差が増大していた場合には、端部付近で最小２乗法のフィッティングに失敗して推定誤差も増大する可能性がある。 However, when the least square method is applied, if the error increases due to the influence of the nature of the received signal near the hit end of the unnecessary phase change φ21 (h), the least square method is used near the end. There is a possibility that the estimation error also increases due to failure of fitting.
そこで、端部付近での最小２乗法のフィッティング失敗の可能性を考慮して、位相フィッティング器１２１の後段に、さらに第２ヒット範囲切出器１５２が設けられている。
第２ヒット範囲切出器１５２は、不要位相変化φ３（ｈ）のヒット端部の付近のデータを捨てて、以下の式（８）のように、最終的に画像再生に用いるヒット幅Ｈの不要位相変化φ２２（ｈ）を取得する。
Therefore, in consideration of the possibility of failure of the least-squares fitting near the end, a second
The second
実際には、まず、最終的なヒット幅Ｈに対して少し広げたヒット幅Ｈ１が設定され、その後、さらに広げたヒット幅Ｈ０が設定される。
ヒット幅Ｈの不要位相変化φ２２（ｈ）のイメージは、図１１（ｅ）に示した通りである。
Actually, a hit width H1 that is slightly expanded with respect to the final hit width H is set first, and then a further expanded hit width H0 is set.
An image of the unnecessary phase change φ22 (h) of the hit width H is as shown in FIG.
以上のように、この発明の実施の形態５（図１０、図１１）による不要位相変化推定器２１Ｄは、第１ヒット範囲切出器１５１と、位相フィッティング器１２１と、第２ヒット範囲切出器１５２と、を備えている。
As described above, the unnecessary
第１ヒット範囲切出器１５１は、一度推定されたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）について、不要位相変化φ１（ｈ）のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅Ｈ以上でレンジ補償後のレンジヒストリのヒット幅Ｈ０以下の中間的なヒット幅Ｈ１（Ｈ≦Ｈ１≦Ｈ０）の不要位相変化φ２１（ｈ）を出力する。
For the unnecessary phase change φ1 (h) included in the range history S1 (r, h) after the range compensation once estimated, the first
位相フィッティング器１２１は、第１ヒット範囲切出器１５１からの不要位相変化φ２１（ｈ）に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を、中間的なヒット幅Ｈ１の不要位相変化φ３（ｈ）として出力する。
第２ヒット範囲切出器１５２は、位相フィッティング器１２１からの不要位相変化φ３（ｈ）のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅Ｈの不要位相変化φ２２（ｈ）を出力する。
The phase
The second
これにより、前述の実施の形態１〜４と同等の効果を奏するとともに、最小２乗法を適用した際にヒット端部で大きくなる可能性がある推定誤差を低減することが可能となる。 As a result, the same effects as those of the first to fourth embodiments described above can be obtained, and estimation errors that can increase at the hit end when the least square method is applied can be reduced.
実施の形態６．
なお、上記実施の形態５（図１０）では、位相フィッティング器１２１の前後に、第１ヒット範囲切出器１５１および第２ヒット範囲切出器１５２を設けたが、図１２のように、不要位相変化推定器２１Ｅ内において、位相フィッティング器１２１の前段に分割ヒット範囲切出器１６１を設け、位相フィッティング器１２１の後段にヒット範囲切出器１１１を設けてもよい。
In the fifth embodiment (FIG. 10), the first
図１２はこの発明の実施の形態６による不要位相変化推定器２１Ｅの構成を示すブロック図であり、図１３は不要位相変化推定器２１Ｅの処理内容を示す説明図である。また、この発明の実施の形態６に係る画像レーダ装置の他の構成は、図１〜図４に示した通りである。
FIG. 12 is a block diagram showing a configuration of an unnecessary
図１２において、不要位相変化推定器２１Ｅは、ＰＧＡ推定器３１の後段に、分割ヒット範囲切出器１６１を備え、さらに、分割ヒット範囲切出器１６１の後段に、位相フィッティング器１２１を介して、ヒット範囲切出器１１１を備えている。
分割ヒット範囲切出器１６１は、ＰＧＡ推定器３１で得られた不要位相変化φ１（ｈ）から、範囲の分割も許容した切出範囲から、不要位相変化φｃ（ｈ）を抽出する。
In FIG. 12, the unnecessary
The division hit
次に、図１〜図４、図１２および図１３を参照しながらこの発明の実施の形態６による動作について説明する。
この場合、不要位相変化推定器２１Ｅ内において、位相フィッティング器１２１の前段に分割ヒット範囲切出器１６１が設けられた点を除けば、他の構成は前述の実施の形態３（図８）と同様である。
Next, the operation according to the sixth embodiment of the present invention will be described with reference to FIGS. 1 to 4, FIG. 12 and FIG.
In this case, except for the point that the split hit
すなわち、レーダ観測器１（図２）によるレンジヒストリＳの取得処理、レーダ画像化器２内のレンジヒストリ切出器１１（図３）によるレンジヒストリＳ０（ｒ，ｈ）の切出処理、レンジ補償器１２（図３、図１１）によるレンジヒストリＳ１（ｒ，ｈ）の取得処理、不要位相変化推定器２１Ｄ内のＰＧＡ推定器３１による不要位相変化φ１（ｈ）の取得処理、位相補償回路２２（図４）によるレンジヒストリＳ３（ｒ，ｈ）の取得処理、および、クロスレンジ圧縮器１４によるレーダ画像Ｕ（ｒ，ｆ）の取得処理は、前述と同様である。
That is, the range history S acquisition process by the radar observer 1 (FIG. 2), the range history S0 (r, h) extraction process by the range history extractor 11 (FIG. 3) in the
また、この場合も、レンジヒストリ切出器１１で切出すヒット幅Ｈ０は、最終的に画像再生に用いるヒット幅Ｈ以上に設定されており、レンジヒストリＳ１（ｒ，ｈ）および不要位相変化φ１（ｈ）のヒット幅はＨ０であり、不要位相変化φ２（ｈ）のヒット幅はＨである。
Also in this case, the hit width H0 cut out by the
また、前述と同様に、ＰＧＡ法を用いて不要位相変化を推定した際に、受信信号の性質によっては不要位相変化の推定誤差が増大する可能性があることを想定し、この影響を低減することを目的としている。
ここでは、特に、或る複数の区分的なヒット範囲において不要位相変化φ１（ｈ）の推定誤差が増大するような場合を想定し、この場合に不要位相変化の推定精度を向上させることを目的としている。
Further, as described above, when an unnecessary phase change is estimated using the PGA method, it is assumed that an estimation error of the unnecessary phase change may increase depending on the characteristics of the received signal, and this influence is reduced. The purpose is that.
Here, in particular, a case is assumed in which the estimation error of the unnecessary phase change φ1 (h) increases in a plurality of piecewise hit ranges. In this case, the purpose is to improve the estimation accuracy of the unnecessary phase change. It is said.
図１３において、図１３（ａ）は不要位相変化φ１（ｈ）のイメージであり、複数の区分的なヒット範囲において不要位相変化φ１（ｈ）の推定誤差が増大するような状況を示している。 In FIG. 13, FIG. 13 (a) is an image of the unnecessary phase change φ1 (h), and shows a situation where the estimation error of the unnecessary phase change φ1 (h) increases in a plurality of piecewise hit ranges. .
図１３（ａ）において、特定ヒット範囲Ａ、Ｂでは、推定誤差が増大していると考えられ、このような不要位相変化の推定誤差が大きい特定ヒット範囲Ａ、Ｂのデータで位相補償を行うと、画像のぼけが残存することになるうえ、場合によっては画像のぼけ量が位相補償前の状態よりも大きくなる可能性がある。 In FIG. 13A, it is considered that the estimation error is increased in the specific hit ranges A and B, and phase compensation is performed with the data of the specific hit ranges A and B having a large estimation error of such unnecessary phase change. As a result, blurring of the image remains, and in some cases, the blurring amount of the image may be larger than that before the phase compensation.
そこで、分割ヒット範囲切出器１６１は、このような推定誤差が大きくなっている可能性がある特定ヒット範囲Ａ、Ｂのデータを捨てて、推定誤差が小さいと期待されるヒット範囲のデータのみを切出して、不要位相変化φｃ（ｈ）として位相フィッティング器１２１に入力する。
Therefore, the divided hit
ただし、不要位相変化φｃ（ｈ）におけるデータ存在ヒット範囲は、これまでのように連続的な範囲ではなく、図１３（ｂ）に示すように、飛び飛びの範囲になることも有り得る。また、図１３（ｂ）内の特定ヒット範囲Ａのように、データ欠落範囲が、最終的な画像再生に用いるヒット幅Ｈ内に含まれることも許容する。 However, the data presence hit range in the unnecessary phase change φc (h) is not a continuous range as described above, but may be a jump range as shown in FIG. Further, as in the specific hit range A in FIG. 13B, the data missing range is allowed to be included in the hit width H used for final image reproduction.
続いて、位相フィッティング器１２１は、図１３（ｃ）に示すように、不要位相変化φｃ（ｈ）に対して、あらかじめ設定した次数の最小２乗法を適用して、２乗誤差を最小とした不要位相変化φ３（ｈ）を取得する。
Subsequently, as shown in FIG. 13C, the phase
最小２乗法においては、推定された各次元の係数に基づいて、任意のヒットｈの値を計算することができるので、不要位相変化φ３（ｈ）においては、不要位相変化φｃ（ｈ）でデータが欠落した特定ヒット範囲Ａ、Ｂの領域のデータ（実線参照）も、当然ながら取得することが可能となる。 In the least square method, the value of an arbitrary hit h can be calculated based on the estimated coefficient of each dimension. Therefore, in the unnecessary phase change φ3 (h), data is represented by the unnecessary phase change φc (h). Naturally, it is also possible to acquire the data (refer to the solid line) of the specific hit ranges A and B in which is missing.
すなわち、不要位相変化φ３（ｈ）は、ｈ＝−Ｈ０／２、・・・、０、・・・、Ｈ０／２−１のヒット範囲で得られる。
以下、前述と同様に、ヒット範囲切出器１１１は、不要位相変化φ３（ｈ）から、最終的に画像再生に使うヒット幅Ｈの不要位相変化φ２（ｈ）を抽出する。
That is, the unnecessary phase change φ3 (h) is obtained in the hit range of h = −H0 / 2,..., 0,.
Thereafter, as described above, the
以上のように、この発明の実施の形態６（図１２、図１３）による不要位相変化推定器２１Ｅは、分割ヒット範囲切出器１６１と、位相フィッティング器１２１と、ヒット範囲切出器１１１と、を備えている。
As described above, the unnecessary
分割ヒット範囲切出器１６１は、一度推定されたレンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）について、不要位相変化φ１（ｈ）から、不要位相変化φ１（ｈ）の推定精度の予測結果を踏まえて設定された特定ヒット範囲Ａ、Ｂのデータを捨てて、特定ヒット範囲Ａ、Ｂ以外のヒット範囲の不要位相変化φｃ（ｈ）を抽出する。
The division hit
位相フィッティング器１２１は、特定ヒット範囲Ａ、Ｂでデータが欠落する分割ヒット範囲切出器１６１からの不要位相変化φｃ（ｈ）に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化φ３（ｈ）の推定値として出力する
ヒット範囲切出器１１１は、位相フィッティング器１２１からの不要位相変化φ３（ｈ）から、最終的に必要なヒット幅Ｈの不要位相変化φ２（ｈ）を抽出して、位相補償回路２２に入力する。
The phase
これにより、前述の実施の形態１〜４と同様の効果を奏するとともに、最小２乗法を適用した際にヒット端部で大きくなる可能性がある推定誤差を低減することができる。
特に、複数の区分的な特定ヒット範囲Ａ、Ｂで不要位相変化φ１（ｈ）の推定誤差が増大するような場合を想定して、特定ヒット範囲Ａ、Ｂのデータを捨てることにより、不要位相変化の推定精度を向上させることができる。
As a result, the same effects as those of the first to fourth embodiments described above can be achieved, and estimation errors that may increase at the hit end when the least square method is applied can be reduced.
In particular, assuming that the estimation error of the unnecessary phase change φ1 (h) increases in a plurality of piecewise specific hit ranges A and B, the unnecessary phase is discarded by discarding the data of the specific hit ranges A and B. The estimation accuracy of the change can be improved.
なお、上記実施の形態１〜６による不要位相変化推定器２１、２１Ａ〜２１Ｅは、レンジ補償後のレンジヒストリＳ１（ｒ，ｈ）に含まれる不要位相変化φ１（ｈ）を、公知のＰＧＡ（Ｐｈａｓｅ Ｇｒａｄｉｅｎｔ Ａｕｔｏｆｏｃｕｓ）法の原理を利用して推定するＰＧＡ推定器３１を備えている。
これにより、一般的なＰＧＡ法で位相補償する場合に発生する補償誤差（換言すれば、不要位相変化の推定誤差）を低減することができる。
Note that the unnecessary
As a result, it is possible to reduce a compensation error (in other words, an estimation error of an unnecessary phase change) that occurs when phase compensation is performed by a general PGA method.
１ レーダ観測器、２ レーダ画像化器、３ 送信機、４ 送受切換器、５ 送受信アンテナ、６ 受信機、７ レンジ圧縮器、１１ レンジヒストリ切出器、１２ レンジ補償器、１３ 位相補償器、１４ クロスレンジ圧縮器、２１、２１Ａ〜２１Ｅ 不要位相変化推定器、２２ 位相補償回路、３１ ＰＧＡ推定器、１１１ ヒット範囲切出器、１２１ 位相フィッティング器、１５１ 第１ヒット範囲切出器、１５２ 第２ヒット範囲切出器、１６１ 分割ヒット範囲切出器、Ａ、Ｂ 特定ヒット範囲、Ｈ 最適的なヒット幅、Ｈ０ 切出後のヒット幅、Ｈ１ 中間的なヒット幅、ｒ レンジ（レンジセル）、Ｓ、Ｓ０（ｒ，ｈ）、Ｓ１（ｒ，ｈ）、Ｓ２（ｒ，ｈ）、Ｓ３（ｒ，ｈ） レンジヒストリ、Ｔ 目標、Ｕ（ｒ，ｆ） レーダ画像、φ１（ｈ）、φ２（ｈ）、φ２１（ｈ）、φ２２（ｈ）、φ３（ｈ）、φｃ 不要位相変化。 1 radar observation device, 2 radar imager, 3 transmitter, 4 transmission / reception switching device, 5 transmission / reception antenna, 6 receiver, 7 range compressor, 11 range history extractor, 12 range compensator, 13 phase compensator, 14 cross range compressor, 21, 21A-21E unnecessary phase change estimator, 22 phase compensation circuit, 31 PGA estimator, 111 hit range extractor, 121 phase fitting unit, 151 first hit range extractor, 152 2 hit range extractor, 161 split hit range extractor, A, B specific hit range, H optimal hit width, H0 hit width after extraction, H1 intermediate hit width, r range (range cell), S, S0 (r, h), S1 (r, h), S2 (r, h), S3 (r, h) Range history, T target, U (r, f) Radar image, φ1 (h , Φ2 (h), φ21 (h), φ22 (h), φ3 (h), φc unnecessary phase change.
Claims (9)
前記レーダ観測器は、
高周波信号を発生する送信機と、
前記高周波信号を前記目標に照射するとともに前記目標で散乱された前記高周波信号の一部を受信する送受信アンテナと、
前記送受信アンテナを介して前記目標からの受信信号を取り込む受信機と、
前記送信機からの高周波信号を前記送受信アンテナに送出するとともに前記送受信アンテナからの受信信号を前記受信機に入力する送受切換器と、
前記受信信号のレンジ方向の分解能をレンジ圧縮により向上させて前記レンジプロフィールを取得するレンジ圧縮器と、を備え、
前記レンジプロフィールを取得する一連の処理を、前記目標と前記レーダ観測器との間の相対位置関係を変えながら繰り返し実行して前記レンジヒストリを取得し、
前記レーダ画像化器は、
前記レンジヒストリをクロスレンジ圧縮することにより、前記目標上の反射強度分布をレンジドップラー分布として表現した前記レーダ画像を得るために、
最終的にクロスレンジ圧縮に用いるヒット幅以上のヒット幅のレンジヒストリに対して、前記目標と前記レーダ観測器との間の相対運動で発生する観測中の各反射点のレンジ分解能セルの移動を推定して補償するレンジ補償器と、
前記レンジ補償器によるレンジ補償後のレンジヒストリ上の代表的な複数の反射点の位相変化を調べ、前記レーダ画像のドップラー周波数方向のぼけの原因となるヒットに対する２次以上の不要位相変化の、前記最終的にクロスレンジ圧縮に用いるヒット幅の範囲での値を推定する不要位相変化推定器と、
前記不要位相変化推定器で得られた不要位相変化に基づき、前記レンジ補償後のレンジヒストリから前記最終的にクロスレンジ圧縮に用いるヒット幅分だけ切出したレンジヒストリに含まれる不要位相変化成分を補償する位相補償回路と、
を備えたことを特徴とする画像レーダ装置。 An image radar apparatus comprising a radar observer that obtains a target range profile and a radar imager that generates a radar image based on a range history of the range profile,
The radar observer is
A transmitter that generates a high-frequency signal;
A transmitting and receiving antenna that irradiates the target with the high-frequency signal and receives a portion of the high-frequency signal scattered by the target;
A receiver that captures a received signal from the target via the transmit / receive antenna;
A transmission / reception switch for sending a high-frequency signal from the transmitter to the transmission / reception antenna and receiving a reception signal from the transmission / reception antenna;
A range compressor that obtains the range profile by improving the resolution in the range direction of the received signal by range compression; and
A series of processing for acquiring the range profile is repeatedly executed while changing the relative positional relationship between the target and the radar observation device, and the range history is acquired.
The radar imager is
In order to obtain the radar image expressing the reflection intensity distribution on the target as a range Doppler distribution by cross-range compressing the range history,
Finally, for the range history with a hit width greater than the hit width used for cross-range compression, move the range resolution cell of each reflection point under observation generated by relative movement between the target and the radar observer. A range compensator to estimate and compensate;
The phase change of a plurality of representative reflection points on the range history after the range compensation by the range compensator is examined, and the second or higher-order unnecessary phase change with respect to a hit that causes blur in the Doppler frequency direction of the radar image. An unnecessary phase change estimator for estimating a value in a range of hit widths used for the final cross range compression;
Based on the unnecessary phase change obtained by the unnecessary phase change estimator, the unnecessary phase change component included in the range history extracted from the range history after the range compensation by the hit width used for the final cross range compression is compensated. A phase compensation circuit to
An image radar apparatus comprising:
一度推定された前記レンジ補償後のレンジヒストリに含まれる不要位相変化について、前記不要位相変化のヒット両端付近の推定値を捨てて、前記最終的にクロスレンジ圧縮に必要なヒット幅の不要位相変化のみを抽出するヒット範囲切出器を含むことを特徴とする請求項１に記載の画像レーダ装置。 The unnecessary phase change estimator is
For the unnecessary phase change included in the range history after the range compensation once estimated, the estimated value near the both ends of the hit of the unnecessary phase change is discarded, and the unnecessary phase change of the hit width finally required for the cross range compression is performed. The image radar apparatus according to claim 1, further comprising a hit range extractor that extracts only the signal.
一度推定された前記レンジ補償後のレンジヒストリに含まれる不要位相変化に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化の推定値として出力する位相フィッティング器を含むことを特徴とする請求項１に記載の画像レーダ装置。 The unnecessary phase change estimator is
By applying a predetermined order least square method to the unnecessary phase change included in the range history after the range compensation once estimated, the phase change that minimizes the square error is estimated as the unnecessary phase change. The image radar apparatus according to claim 1, further comprising a phase fitting unit that outputs a value.
入力された前記不要位相変化に対して、前記不要位相変化のヒット両端付近の推定値を捨てて、前記最終的にクロスレンジ圧縮に必要なヒット幅の不要位相変化を抽出するヒット範囲切出器と、
入力された前記不要位相変化に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化の推定値として出力する位相フィッティング器と、
を含むことを特徴とする請求項１に記載の画像レーダ装置。 The unnecessary phase change estimator is
A hit range extractor that discards the estimated values near both ends of the hit of the unnecessary phase change and extracts the unnecessary phase change of the hit width necessary for the cross range compression with respect to the input unnecessary phase change. When,
A phase fitting device that applies a least-square method of a preset order to the input unnecessary phase change and outputs a phase change that minimizes a square error as an estimated value of the unnecessary phase change;
The image radar apparatus according to claim 1, comprising:
一度推定された前記レンジ補償後のレンジヒストリに含まれる不要位相変化に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化の推定値として出力する位相フィッティング器と、
位相フィッティング器からの不要位相変化のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅の不要位相変化を抽出するヒット範囲切出器と、
を含むことを特徴とする請求項４に記載の画像レーダ装置。 The unnecessary phase change estimator is
By applying a predetermined order least square method to the unnecessary phase change included in the range history after the range compensation once estimated, the phase change that minimizes the square error is estimated as the unnecessary phase change. A phase fitting device to output as a value;
A hit range extractor that discards the estimated values near both ends of the hit of the unnecessary phase change from the phase fitting device, and finally extracts the unnecessary phase change of the hit width required for cross range compression,
The image radar apparatus according to claim 4, comprising:
一度推定された前記レンジ補償後のレンジヒストリに含まれる不要位相変化について、前記不要位相変化のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅の不要位相変化を抽出するヒット範囲切出器と、
前記ヒット範囲切出器からの不要位相変化に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化の推定値として出力する位相フィッティング器と、
を含むことを特徴とする請求項４に記載の画像レーダ装置。 The unnecessary phase change estimator is
For the unnecessary phase change included in the range history after the range compensation once estimated, the estimated values near the both ends of the hit of the unnecessary phase change are discarded, and finally the unnecessary phase change of the hit width necessary for the cross range compression is performed. A hit range extractor to extract,
Phase fitting that outputs a phase change that minimizes a square error as an estimated value of an unnecessary phase change by applying a least-square method of a preset order to the unnecessary phase change from the hit range extractor And
The image radar apparatus according to claim 4, comprising:
一度推定された前記レンジ補償後のレンジヒストリに含まれる不要位相変化について、前記不要位相変化のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅以上でレンジ補償後のレンジヒストリのヒット幅以下の中間的なヒット幅の不要位相変化を出力する第１ヒット範囲切出器と、
前記第１ヒット範囲切出器からの不要位相変化に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を、中間的なヒット幅の不要位相変化として出力する位相フィッティング器と、
前記位相フィッティング器からの不要位相変化のヒット両端付近の推定値を捨てて、最終的にクロスレンジ圧縮に必要なヒット幅の不要位相変化を出力する第２ヒット範囲切出器と、
を含むことを特徴とする請求項１に記載の画像レーダ装置。 The unnecessary phase change estimator is
For the unnecessary phase change included in the range history after the range compensation once estimated, the estimated values near both ends of the hit of the unnecessary phase change are discarded, and after the range compensation at the final hit width required for cross range compression A first hit range extractor that outputs an unnecessary phase change with an intermediate hit width equal to or smaller than the hit width of the range history of
A phase change that minimizes the square error is applied to the unnecessary phase change from the first hit range extractor by applying a least-square method of a preset order, and an unnecessary phase having an intermediate hit width is obtained. A phase fitting device that outputs as a change;
A second hit range extractor that discards the estimated value near both ends of the hit of the unnecessary phase change from the phase fitting device and finally outputs the unnecessary phase change of the hit width required for cross range compression;
The image radar apparatus according to claim 1, comprising:
一度推定された前記レンジ補償後のレンジヒストリに含まれる不要位相変化について、前記不要位相変化から、前記不要位相変化の推定精度の予測結果を踏まえて設定された特定ヒット範囲のデータを捨てて、前記特定ヒット範囲以外のヒット範囲の不要位相変化を抽出する分割ヒット範囲切出器と、
前記特定ヒット範囲でデータが欠落する前記分割ヒット範囲切出器からの不要位相変化に対して、あらかじめ設定された次数の最小２乗法を適用して、２乗誤差を最小とする位相変化を不要位相変化の推定値として出力する位相フィッティング器と、
前記位相フィッティング器からの不要位相変化から、最終的に必要なヒット幅の不要位相変化を抽出するヒット範囲切出器と、
を含むことを特徴とする請求項１に記載の画像レーダ装置。 The unnecessary phase change estimator is
For the unnecessary phase change included in the range history after the range compensation once estimated, discard the data of the specific hit range set based on the prediction result of the estimation accuracy of the unnecessary phase change from the unnecessary phase change, A split hit range extractor for extracting unnecessary phase changes in hit ranges other than the specific hit range;
Applying the least-squares method of a preset order to the unnecessary phase change from the divided hit range extractor in which data is missing in the specific hit range, eliminating the phase change that minimizes the square error A phase fitting device that outputs as an estimate of the phase change;
A hit range extractor for extracting an unnecessary phase change of a finally required hit width from an unnecessary phase change from the phase fitting device;
The image radar apparatus according to claim 1, comprising:
前記レンジ補償後のレンジヒストリに含まれる不要位相変化を、特許文献１や特許文献２に記載のＰＧＡ（Ｐｈａｓｅ Ｇｒａｄｉｅｎｔ Ａｕｔｏｆｏｃｕｓ）法の原理を利用して推定するＰＧＡ推定器を含むことを特徴とする請求項１から請求項８までのいずれか１項に記載の画像レーダ装置。 The unnecessary phase change estimator is
It includes a PGA estimator that estimates an unnecessary phase change included in the range history after the range compensation using a principle of a PGA (Phase Gradient Autofocus) method described in Patent Document 1 or Patent Document 2. The image radar device according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010117656A JP2011247593A (en) | 2010-05-21 | 2010-05-21 | Image radar device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010117656A JP2011247593A (en) | 2010-05-21 | 2010-05-21 | Image radar device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2011247593A true JP2011247593A (en) | 2011-12-08 |
Family
ID=45413053
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010117656A Pending JP2011247593A (en) | 2010-05-21 | 2010-05-21 | Image radar device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2011247593A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014020820A (en) * | 2012-07-13 | 2014-02-03 | Mitsubishi Electric Corp | Radar device |
| WO2014092052A1 (en) * | 2012-12-14 | 2014-06-19 | 三菱電機株式会社 | Radar apparatus |
| WO2016148104A1 (en) * | 2015-03-16 | 2016-09-22 | 三菱電機株式会社 | Radar moving image creation device and method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11118913A (en) * | 1997-10-21 | 1999-04-30 | Mitsubishi Electric Corp | Phase compensation circuit, phase compensation method, and radar equipment |
| JP2001194456A (en) * | 2000-01-14 | 2001-07-19 | Mitsubishi Electric Corp | Composite aperture radar |
| JP2004157048A (en) * | 2002-11-07 | 2004-06-03 | Mitsubishi Electric Corp | Radar installation with synthetic aperture |
| JP2007333583A (en) * | 2006-06-15 | 2007-12-27 | Mitsubishi Electric Corp | Image radar device |
| JP2009505068A (en) * | 2005-08-10 | 2009-02-05 | レイセオン カンパニー | Efficient autofocus method for swath SAR |
-
2010
- 2010-05-21 JP JP2010117656A patent/JP2011247593A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11118913A (en) * | 1997-10-21 | 1999-04-30 | Mitsubishi Electric Corp | Phase compensation circuit, phase compensation method, and radar equipment |
| JP2001194456A (en) * | 2000-01-14 | 2001-07-19 | Mitsubishi Electric Corp | Composite aperture radar |
| JP2004157048A (en) * | 2002-11-07 | 2004-06-03 | Mitsubishi Electric Corp | Radar installation with synthetic aperture |
| JP2009505068A (en) * | 2005-08-10 | 2009-02-05 | レイセオン カンパニー | Efficient autofocus method for swath SAR |
| JP2007333583A (en) * | 2006-06-15 | 2007-12-27 | Mitsubishi Electric Corp | Image radar device |
Non-Patent Citations (3)
| Title |
|---|
| JPN6013056829; 山本 和彦 Ｋａｚｕｈｉｋｏ ＹＡＭＡＭＯＴＯ，岩本 雅史 Ｍａｓａｆｕｍｉ ＩＷＡＭＯＴＯ，小菅 義夫 Ｙｏｓｈｉｏ ＫＯＳＵＧＥ: '"ドップラー微係数の時間変化を利用した画像レーダのオートフォーカス法 Ａｎ Ａｕｔｏｆｏｃｕｓ Ａｌｇｏｒｉｔｈｍ ｆｏｒ ｔｈｅ' 電子情報通信学会論文誌Ｂ通信 Ｖｏｌ．Ｊ９０-Ｂ，Ｎｏ．４, 20070401, 社団法人電子情報通信学会 * |
| JPN6013056831; 山本 和彦 Ｋａｚｕｈｉｋｏ ＹＡＭＡＭＯＴＯ，岩本 雅史 Ｍａｓａｆｕｍｉ ＩＷＡＭＯＴＯ，坂爪 俊之 Ｔｏｓｈｉｙｕｋｉ ＳＡＫＡＺＵＭＥ，志村 誠一: '"ドップラーを利用したレーダ画像再生のための位相補償アルゴリズム Ａ Ｎｅｗ Ｐｈａｓｅ Ｃｏｍｐｅｎｓａｔｉｏｎ Ａｌｇｏｒｉｔｈｍ' 電子情報通信学会論文誌Ｂ通信 Ｖｏｌ．Ｊ８３-Ｂ，Ｎｏ．１０, 20001025, 社団法人電子情報通信学会 * |
| JPN6013056833; 山本 和彦 Ｋａｚｕｈｉｋｏ ＹＡＭＡＭＯＴＯ，岩本 雅史 Ｍａｓａｆｕｍｉ ＩＷＡＭＯＴＯ，小菅 義夫 Ｙｏｓｈｉｏ ＫＯＳＵＧＥ: '"二次以上のレンジ変化を考慮したＳＡＲオートフォーカス法 Ａ ＳＡＲ Ａｕｔｏｆｏｃｕｓ Ａｌｇｏｒｉｔｈｍ ｆｏｒ Ｈｉｇｈ Ｏｒｄｅｒ Ｒ' 電子情報通信学会論文誌Ｂ通信 Ｖｏｌ．Ｊ８９-Ｂ，Ｎｏ．７, 20060701, ｐ．１０５４-１０６２, 社団法人電子情報通信学会 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014020820A (en) * | 2012-07-13 | 2014-02-03 | Mitsubishi Electric Corp | Radar device |
| WO2014092052A1 (en) * | 2012-12-14 | 2014-06-19 | 三菱電機株式会社 | Radar apparatus |
| JP5823062B2 (en) * | 2012-12-14 | 2015-11-25 | 三菱電機株式会社 | Radar equipment |
| WO2016148104A1 (en) * | 2015-03-16 | 2016-09-22 | 三菱電機株式会社 | Radar moving image creation device and method |
| JPWO2016148104A1 (en) * | 2015-03-16 | 2017-04-27 | 三菱電機株式会社 | Radar moving image creating apparatus and method |
| US10495749B2 (en) | 2015-03-16 | 2019-12-03 | Mitsubishi Electric Corporation | Radar video creation apparatus and method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8619908B2 (en) | Wireless ranging system and related methods | |
| CN110221268B (en) | Phase-amplitude calibration filter and method for improving short pulse non-coherent radar signal processing performance | |
| JP6080582B2 (en) | Image radar device | |
| CN102520396A (en) | Time-varying narrow-band interference suppression method based on complex empirical mode decomposition | |
| US20140232590A1 (en) | Sar image formation | |
| CN113812979A (en) | System and method for ultrasound elastography with continuous transducer vibration | |
| JP5606097B2 (en) | Passive radar device | |
| WO2009102926A1 (en) | Radar imaging of buildings using model-based focusing and data-based focusing | |
| EP3771921A1 (en) | Method for processing non-interrupted phase synchronization signal of bistatic sar based on coded signal | |
| JP2007040953A (en) | Correlation reception processor | |
| JP4046422B2 (en) | Motion compensation circuit and radar device | |
| JP3645133B2 (en) | Radar equipment | |
| US20150145718A1 (en) | Radar weather data signal processing method and signal processing module | |
| JP4722574B2 (en) | Image radar device | |
| KR20110036430A (en) | Method of generating inverse synthetic aperture radar image compensating rotational motion of target and apparatus thereof | |
| JP5979868B2 (en) | Image radar device | |
| WO2012006226A1 (en) | Processing radar return signals to detect targets | |
| JP2014020820A (en) | Radar device | |
| JP2011247593A (en) | Image radar device | |
| JP4943065B2 (en) | Image radar device | |
| JP2010127771A (en) | Synthetic aperture sonar, and method and program for correcting phase error of synthetic aperture sonar | |
| KR101569673B1 (en) | Method for Reducing Sidelobe In Ultrasound Imaging | |
| JP2015532848A (en) | Image tissue motion estimation | |
| JP4982229B2 (en) | Moving target detection device | |
| KR20180012590A (en) | A Method of ISAR rotational motion compensation for highly focused ISAR images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20121119 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131023 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131119 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140107 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20140805 |