JP2010060448A - Radar image reproduction apparatus - Google Patents
Radar image reproduction apparatus Download PDFInfo
- Publication number
- JP2010060448A JP2010060448A JP2008226965A JP2008226965A JP2010060448A JP 2010060448 A JP2010060448 A JP 2010060448A JP 2008226965 A JP2008226965 A JP 2008226965A JP 2008226965 A JP2008226965 A JP 2008226965A JP 2010060448 A JP2010060448 A JP 2010060448A
- Authority
- JP
- Japan
- Prior art keywords
- range
- azimuth
- compression
- frequency
- cell migration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
この発明は、人工衛星や航空機等の移動プッラットフォームに搭載される合成開口レーダ等のレーダ画像再生装置に関する。 The present invention relates to a radar image reproducing apparatus such as a synthetic aperture radar mounted on a mobile platform such as an artificial satellite or an aircraft.
合成開口レーダ観測では、アジマス方向の高分解能化とレンジ方向の観測領域の拡大を同時に行うことは難しい。これはそれぞれを実現するために必要なパルス繰り返し周波数(PRF:Pulse Repetition Frequency)に対する要求が相反するためである。アジマス方向の高分解能化には、ドップラー帯域以上の高いPRFが要求される。一方、レンジ方向の観測領域の拡大には、観測領域の拡大に伴うエコーの受信時間増加に応じて、低いPRFが要求される。そこで、レンジ方向の観測領域拡大とアジマス方向の高分解能化を同時に実現すべく、複数のレーダビームを用いる合成開口レーダの観測方式がある(以下、複数ビーム観測方式と呼ぶ)。例えば、図1に示すように、1つの送信アンテナと2つの受信アンテナを用いる複数ビーム観測方式では、物理的なPRFは送信アンテナでのパルス照射間隔(PRI:Pulse Repetition Interval)の逆数となるが、実効的なPRFをその2倍(つまり、PRIは半分)とすることができる。物理的なPRFを低く設定できるためレンジ領域の観測領域を狭める必要がない。そして、等価的に幅の広い1つのレーダビームを形成することでアジマス分解能を改善している。 In synthetic aperture radar observation, it is difficult to simultaneously increase the resolution in the azimuth direction and expand the observation area in the range direction. This is because the requirements for the pulse repetition frequency (PRF) necessary to realize each of these are contradictory. Higher resolution in the azimuth direction requires a higher PRF than the Doppler band. On the other hand, in order to expand the observation area in the range direction, a low PRF is required in accordance with an increase in echo reception time accompanying the expansion of the observation area. Therefore, there is a synthetic aperture radar observation method using a plurality of radar beams (hereinafter referred to as a “multiple beam observation method”) in order to simultaneously realize an observation region expansion in the range direction and a high resolution in the azimuth direction. For example, as shown in FIG. 1, in the multiple beam observation method using one transmission antenna and two reception antennas, the physical PRF is the reciprocal of the pulse irradiation interval (PRI: Pulse Repetition Interval) at the transmission antenna. Effective PRF can be doubled (ie, PRI is half). Since the physical PRF can be set low, there is no need to narrow the observation area in the range area. The azimuth resolution is improved by forming one radar beam that is equivalently wide.
この複数レーダビーム観測方式のための従来のレーダ画像再生装置は、複数の受信アンテナで得た観測信号をスペクトル上で合成し、等価的な1つの受信アンテナで得た観測信号を生成した後で、画像再生処理を行っていた(例えば、非特許文献1)。各受信アンテナで得られた観測信号のアジマス方向のスペクトルは、図2に示すように、物理的なPRFで決まるサンプリング帯域幅が受信アンテナのビームの幅で決まるアジマス周波数(ドップラー周波数)帯域よりも狭いため、エイリアシング誤差を持つ。そこで、スペクトル上での合成において、このエイリアシング誤差を相殺して解消する復元フィルタ(Reconstruction Filter)を適用し、エイリアシング誤差を含まないスペクトルを得るようにしている。 The conventional radar image reproducing apparatus for the multiple radar beam observation method synthesizes observation signals obtained by a plurality of reception antennas on a spectrum and generates an observation signal obtained by an equivalent single reception antenna. The image reproduction processing was performed (for example, Non-Patent Document 1). As shown in FIG. 2, the spectrum in the azimuth direction of the observation signal obtained by each receiving antenna is more than the azimuth frequency (Doppler frequency) band in which the sampling bandwidth determined by the physical PRF is determined by the beam width of the receiving antenna. Because it is narrow, it has aliasing errors. Therefore, in the synthesis on the spectrum, a reconstruction filter that cancels and eliminates the aliasing error is applied to obtain a spectrum that does not include the aliasing error.
従来のレーダ画像再生装置は、レンジ圧縮後の2次元信号上で復元フィルタを適用しエイリアシング誤差を解消していた。レンジドップラー領域において、レンジ圧縮後の1散乱点の信号は、レンジセルマイグレーション(RCM:Range Cell Migration)を持ちレンジ方向に広がっている。このレンジ方向の範囲がM点、アジマス周波数方向の範囲がN点であるとすると、L個の受信アンテナがある場合、2次元信号上で復元フィルタを適用するためには、1レンジ当りM×N×Lの演算量が必要であり、計算負荷が高いという問題点があった。さらに、復元フィルタを適用して合成した信号を格納するためにM×N×Lの配列を主記憶容量中に確保する必要があり、計算機要求が高いという問題があった。
この発明は上記のような問題点を解決するためになされたもので、演算量が少なく、かつ、大きな配列を必要としない、計算負荷および計算要求の小さなレーダ画像再生装置を得ることを目的とする。
A conventional radar image reproducing apparatus applies an restoration filter on a two-dimensional signal after range compression to eliminate an aliasing error. In the range Doppler region, the signal at one scattering point after range compression has range cell migration (RCM) and spreads in the range direction. Assuming that the range in the range direction is M points and the range in the azimuth frequency direction is N points, when there are L receiving antennas, M × There is a problem that the calculation amount of N × L is required and the calculation load is high. Furthermore, in order to store the signal synthesized by applying the restoration filter, it is necessary to secure an M × N × L array in the main storage capacity, and there is a problem that the computer requirement is high.
The present invention has been made to solve the above-described problems, and an object of the present invention is to obtain a radar image reproducing device with a small calculation load and a small calculation load and a low calculation requirement. To do.
この発明に係るレーダ画像再生装置は、複数の受信アンテナで同時受信された観測信号をそれぞれ格納する複数の格納部と、各受信アンテナで得られ、格納部に格納された信号をそれぞれレンジ圧縮する複数のレンジ圧縮部と、各レンジ圧縮部でレンジ圧縮された信号をそれぞれレンジセルマイグレーション補正する複数のレンジセルマイグレーション補正部と、各レンジセルマイグレーション補正結果をそれぞれコヒーレントに加算する複数のコヒーレント加算部と、各コヒーレント加算部の結果にそれぞれ復元フィルタを適用する複数の復元フィルタ部と、各受信アンテナで得られた信号に対するそれぞれの復元フィルタの複数の適用結果を合成しアジマス方向のエイリアシングエラーを相殺させ解消する合成部と、前記合成部の出力に対しアジマス圧縮を行うアジマス圧縮部を備える。 The radar image reproducing apparatus according to the present invention includes a plurality of storage units each storing observation signals simultaneously received by a plurality of reception antennas, and a range compression of the signals obtained by each reception antenna and stored in the storage unit. A plurality of range compression units, a plurality of range cell migration correction units that perform range cell migration correction on signals that have been range compressed by each range compression unit, and a plurality of coherent addition units that add each range cell migration correction result coherently, respectively And multiple restoration filter units that apply the restoration filter to the result of each coherent adder unit, and multiple application results of each restoration filter for the signal obtained by each receiving antenna to cancel aliasing errors in the azimuth direction To the output of the synthesis unit Comprising were azimuth compression unit to perform azimuth compression.
この発明によるレーダ画像再生装置によれば、復元フィルタを適用する前に、エイリアシングを考慮したレンジセルマイグレーション補償を行い、信号のレンジ方向の広がりを補償しているので、復元フィルタを適用するデータ量を(1/Mに)削減することができ、このデータ量の削減により、復元フィルタ処理のための演算量を(1/Mに)減らし、処理を高速化することができる。
また、復元フィルタを適用するデータ量の削減により、復元フィルタ適用後の合成処理に必要な配列サイズを小さくすることができ、これに伴い配列の確保に必要な主記憶容量を抑圧でき、計算機要求を緩和できる。
According to the radar image reproducing apparatus of the present invention, before applying the restoration filter, the range cell migration compensation is performed in consideration of aliasing, and the spread in the range direction of the signal is compensated. Therefore, the amount of data to which the restoration filter is applied (1 / M) can be reduced, and by reducing the amount of data, the amount of computation for restoration filter processing can be reduced (to 1 / M), and the processing speed can be increased.
In addition, by reducing the amount of data to which the restoration filter is applied, the array size required for the composition processing after the restoration filter is applied can be reduced, and as a result, the main storage capacity necessary to secure the array can be suppressed, and the computer request Can be relaxed.
実施の形態1.
図3はこの発明の実施の形態1によるレーダ画像再生装置を示すブロック構成図である。なお、ここでは便宜上、受信アンテナが2つある場合を想定して説明する。複数受信アンテナがある場合は、各受信アンテナに応じた処理ブロックを受信アンテナ個数分用意すればよい。
3 is a block diagram showing a radar image reproducing apparatus according to
図3において、Rx1格納部1は受信アンテナ1で得られた観測信号を格納し、Rx2格納部2は受信アンテナ2で得られた観測信号を格納する。出力格納部3は、本装置の出力を格納する。レンジ圧縮部11とレンジ圧縮部12は、それぞれ受信アンテナ1と2で得られた観測信号をレンジ圧縮および2次レンジ圧縮する。折り返しを考慮したRCM補正部13と折り返しを考慮したRCM補正部14は、それぞれ対応するアンテナのレンジ圧縮された信号に対しレンジセルマイグレーション補正をする。コヒーレント加算部15とコヒーレント加算部16は、それぞれ対応するアンテナのレンジセルマイグレーション補正後の信号における、エイリアシングした部分とエイリアシングしていない部分をコヒーレントに加算する。復元フィルタ部P1(fa;R0)17と復元フィルタ部P2(fa;R0)18は、それぞれ対応するアンテナの信号のアジマススペクトルに対し復元フィルタを適用する。合成部19は、復元フィルタ適用後の各受信アンテナの信号をアジマススペクトル上で合成する。アジマス圧縮部20は、合成後の信号に対しアジマス圧縮処理を施す。
In FIG. 3, the
なお、この明細書の説明において、部位という語は、専用の電子回路または素子を意味しているが、汎用的な中央演算装置(CPU:Central Processing Unit)を搭載したコンピュータに所定の処理を実行させるコンピュータプログラムモジュールの形で構成するようにしても良い。 In the description of this specification, the term “part” means a dedicated electronic circuit or element, but performs predetermined processing on a computer equipped with a general-purpose central processing unit (CPU). You may make it comprise in the form of the computer program module to be made.
次に動作について説明する。
図4は、実施の形態1によるレーダ画像再生装置の処理の流れを示すフローチャートである。
まずステップST1において、レンジ圧縮部11とレンジ圧縮部12が、それぞれ受信アンテナ1と2で得られた観測信号に対しレンジ圧縮を行う。このレンジ圧縮は、レンジ方向の周波数領域で、式1で与えられる参照関数を乗ずることで実現する。
Next, the operation will be described.
FIG. 4 is a flowchart showing a processing flow of the radar image reproducing apparatus according to the first embodiment.
First, in step ST1, the
ここで、frはレンジ方向周波数を表し、Krは送信パルスのFM(Frequency Modulation)レートを表す。
なお、式1は周波数領域におけるレンジ方向のマッチドフィルタを示すが、レンジ周波数領域に限るものではなく、レンジ時間領域でマッチドフィルタを生成し、レンジ圧縮しても良い。
Here, f r represents the range direction frequency, K r denotes the FM (Frequency Modulation) rate of the transmitted pulse.
In addition, although Formula 1 shows the matched filter of the range direction in a frequency domain, it is not restricted to a range frequency domain, A matched filter may be produced | generated in a range time domain, and range compression may be carried out.
このレンジ圧縮の結果、レンジドップラー領域におけるレンジR0に位置する散乱点の信号の軌跡は、図4中100に示すようになる。信号の軌跡の内、アジマス周波数 faの範囲が(−PRF/2,PRF/2) である成分は双曲線状の弧を描くが、fa<−PRF/2 の成分とfa>PRF/2 の成分はエイリアシングして現れる。この内fa<−PRF/2 の成分は PRF/2 側に現れ、fa>PRF/2 の成分は−PRF/2 側に現れる。そして、これらエイリアシングした部分には、エイリアシング誤差が発生する。ここで、図4中100に示す、MはレンジRCMによるデータのレンジ方向の広がった範囲(データ点数)を示し、Nはアジマス方向の範囲(データ点数)に相当する。 As a result of the range compression, the signal locus of the scattering point located in the range R0 in the range Doppler region is as indicated by 100 in FIG. Among trajectory signals, the range of azimuth frequency f a is the component which is (-PRF / 2, PRF / 2 ) draws a hyperbolic arc but, f a <-PRF / 2 component and f a> PRF / The second component appears with aliasing. Among these, the component of f a <−PRF / 2 appears on the PRF / 2 side, and the component of f a > PRF / 2 appears on the −PRF / 2 side. An aliasing error occurs in these aliased portions. Here, M shown in 100 in FIG. 4 indicates a range (data points) in the range direction of the data by the range RCM, and N corresponds to a range (data points) in the azimuth direction.
次に、ST2において、再生処理対象となるデータのレンジに対するループ処理に入る。このループ処理では、後述する2次レンジ圧縮処理、復元フィルタ処理およびアジマス圧縮処理におけるレンジR0をループ毎に更新しながら設定する。このレンジの更新間隔はレンジビン間隔にするのが適当であるが、これに限られるものではなく、ユーザが任意に設定してよい。 Next, in ST2, loop processing for the range of data to be reproduced is entered. In this loop processing, the range R0 in the secondary range compression processing, restoration filter processing, and azimuth compression processing described later is set while being updated for each loop. The range update interval is suitably the range bin interval, but is not limited to this, and the user may arbitrarily set it.
そして、ステップST3において、データのアジマス周波数の区画に対するループ処理に入る。このループ処理は、データのアジマス周波数の範囲をループ毎に更新して設定する。本実施の形態1における説明のように、受信アンテナが2つの場合は、設定されるアジマス周波数faの範囲は(−3PRF/2,−PRF/2),(−PRF/2, PRF/2)および(PRF/2,3PRF/2)の3区画である。なお、受信アンテナがL個の場合は、アジマス周波数の範囲は式2のように(2L-1)区画となる。
Then, in step ST3, a loop process for the azimuth frequency section of the data is entered. In this loop processing, the azimuth frequency range of data is updated and set for each loop. As described in the first embodiment, when the receiving antenna is two, the range of azimuth frequency f a is set (-3PRF / 2, -PRF / 2 ), (- PRF / 2, PRF / 2 ) And (PRF / 2, 3PRF / 2). When the number of receiving antennas is L, the azimuth frequency range is (2L-1) sections as shown in
ここでは、便宜上アジマス中心周波数(ドップラー中心周波数)fdcを0として記述したが、これに限るものではない。ビームスクイントでアジマス中心周波数fdcが非ゼロとなる場合は、上記アジマス周波数の範囲にfdcをオフセットとして加えれば良い。 Here, the azimuth center frequency (Doppler center frequency) f dc is described as 0 for convenience, but the present invention is not limited to this. If azimuth center frequency f dc becomes nonzero Beams Quint, may be added to f dc as an offset to the range of the azimuth frequency.
ステップST4において、レンジ圧縮部11とレンジ圧縮部12がそれぞれ受信アンテナ1と2の信号を2次元周波数空間で2次レンジ圧縮処理する。この2次レンジ圧縮処理は式3で与えられる関数を乗じることで実現する。
In step ST4, the
ここで、faはアジマス方向周波数を示す。また、R0は処理対象レンジのレンジを示す。そして、KSTC(R0,fa)は式4で与えられる変数である。 Here, f a indicates the azimuth direction frequency. Also, R 0 represents a range of the process target range. K STC (R 0 , f a ) is a variable given by Equation 4.
ここで、V はレーダの速度、fcはレーダ送信周波数、cは光の速度、λはレーダ波長である。
なお、この2次レンジ圧縮処理は必須ではない。ビームスクイント角が小さい、あるいは、アジマス方向分解能が小さく、レンジとアジマス方向のカップリングが問題とならない場合は、この2次レンジ圧縮処理を省略しても良い。
Here, V is the speed of the radar, is f c radar transmission frequency, c is the speed of light, lambda is the radar wavelength.
This secondary range compression process is not essential. When the beam squint angle is small, or the resolution in the azimuth direction is small and the coupling between the range and the azimuth direction does not cause a problem, the secondary range compression process may be omitted.
ステップST5において、折り返しを考慮したRCM補正部13と折り返しを考慮したRCM補正部14が、それぞれ2次レンジ圧縮後の受信アンテナ1と2のデータに対しRCM補正を行う。ここでは、先の2次レンジ圧縮で設定した3つのアジマス周波数範囲において、式5を用いてアジマス周波数faに対するレンジを算出する。
In step ST5, the
この操作を、図5を用いて説明する。同図は、式5で算出されるアジマス周波数faに応じたレンジR(fa,R0)を示す模式図である。同図に示すように、式5を用いて算出したアジマス周波数faに対するレンジは、アジマス周波数の範囲に応じて異なる。そして、このレンジの軌跡は、レンジドップラー領域におけるレンジ圧縮および2次レンジ圧縮後のレンジR0に位置する散乱点の信号の軌跡に一致する。そこで、式5に基づき、2次レンジ圧縮時と同じアジマス周波数の範囲の設定でレンジR(fa,R0)に位置するデータを集めることで、折り返しを考慮したRCM補正が可能となる。
この結果、図4中の200に示すように、アジマス周波数faの範囲 (−3PRF/2,−PRF/2),(−PRF/2, PRF/2) および (PRF/2,3PRF/2)のそれぞれについてRCM補正結果を得る。
This operation will be described with reference to FIG. This figure is a schematic diagram showing a range R (f a , R 0 ) corresponding to the azimuth frequency f a calculated by equation (5). As shown in the figure, the range with respect to the azimuth frequency f a calculated using Equation 5 differs depending on the range of the azimuth frequency. The locus of this range corresponds to the locus of the signal of scattered points located range R 0 after range compression and secondary range compression in the range Doppler region. Therefore, RCM correction considering aliasing can be performed by collecting data located in the range R (f a , R 0 ) with the same azimuth frequency range setting as in the secondary range compression based on Equation 5.
As a result, as shown in 200 in FIG. 4, the range of azimuth frequency f a (-3PRF / 2, -PRF / 2), (- PRF / 2, PRF / 2) and (PRF / 2,3PRF / 2 ) To obtain the RCM correction result.
ステップST6において、全てのアジマス周波数の区画に対し、ステップST4からST5までの処理が適用された場合は、アジマス周波数の区画に対するループを終了する。そうでない場合は、ループを継続する。
ステップST7において、コヒーレント加算部15とコヒーレンス加算部16が、それぞれ受信アンテナ1と2の信号に対し、前述したアジマス周波数faの3つの区画におけるRMC補正結果をコヒーレントに加算する。この結果、図4中300に示すように、RCM補正前にM×N点に広がっていたデータの範囲はRCM補正後に1×N点となる。
ステップST8において、復元フィルタ部P1(fa;R0)17と復元フィルタ部P2(fa;R0)18が、RCM補正されたそれぞれの受信アンテナの信号に対し復元フィルタを適用する。この受信アンテナ1と2に対する、復元フィルタは例えば、式6で与えられる。
In step ST6, when the processing from step ST4 to ST5 is applied to all the azimuth frequency sections, the loop for the azimuth frequency section is terminated. If not, continue the loop.
In step ST7, the
In step ST8, the restoration filter unit P1 (fa; R0) 17 and the restoration filter unit P2 (fa; R0) 18 apply the restoration filter to the signals of the respective reception antennas subjected to the RCM correction. The restoration filter for the receiving
なお、ここでは2つの受信アンテナを用いる場合の復元フィルタを式6で与えたが、これに限られるものではなく、同等の効果を持つ式で与えられる復元フィルタを適用しても良い。 Here, the restoration filter in the case of using two receiving antennas is given by Expression 6, but the present invention is not limited to this, and a restoration filter given by an expression having an equivalent effect may be applied.
ステップST9において、合成部19が、アジマススペクトル上で復元フィルタ適用後の各受信アンテナの信号を合成し、エイリアシング誤差を相殺することで解消する。この結果、図4中400に示すように、合成後の信号の範囲は1×2N点(受信アンテナがL個の場合は1×N×L点)となる。
In step ST9, the combining
ステップST10において、全てのレンジビンに対し、ステップST3からステップST9までの処理が適用された場合は、データのレンジに対するループを終了する。そうでない場合は、全レンジビンに対して一連の処理が終わるまでループを続ける。
ステップST11において、アジマス圧縮部20が、合成後の信号に対しアジマス圧縮処理を施す。アジマス圧縮処理は、レンジドップラー領域において式7で与えられる関数を乗ずることで実現する。
If the processing from step ST3 to step ST9 is applied to all the range bins in step ST10, the loop for the data range is terminated. Otherwise, the loop is continued until a series of processing is completed for all range bins.
In step ST11, the
なお、式7は周波数領域におけるアジマス方向のマッチドフィルタを示すが、アジマス圧縮処理はアジマス周波数領域のマッチドフィルタ処理に限るものではなく、アジマス時間領域でマッチドフィルタを生成し、アジマス時間領域でアジマス圧縮しても良い。
また、サブアパーチャアルゴリズムによるアジマス圧縮処理を導入し、データをアジマス方向で合成開口長以下に分割、アジマス圧縮処理した上で結合する手順をとっても良い。
そして、このアジマス圧縮処理によりレーダ画像再生が完了し、レーダ画像を得る。
Note that Equation 7 shows a matched filter in the azimuth direction in the frequency domain, but the azimuth compression processing is not limited to the matched filter processing in the azimuth frequency domain. You may do it.
Further, it is possible to adopt a procedure in which azimuth compression processing using a sub-aperture algorithm is introduced, data is divided in the azimuth direction to be equal to or shorter than the synthetic aperture length, and after azimuth compression processing is combined.
The radar image reproduction is completed by this azimuth compression processing, and a radar image is obtained.
以上のように、復元フィルタを適用する前に、エイリアシングを考慮したレンジセルマイグレーション補償を行い、信号のレンジ方向の広がりを補償しているので、復元フィルタを適用するデータ量を(1/Mに)削減することができる。
このデータ量の削減により、復元フィルタ処理のための演算量を(1/Mに)減らし、処理を高速化することができる。
また、復元フィルタを適用するデータ量を削減しているので、復元フィルタ適用後の合成処理に必要な配列サイズを小さくすることができ、これに伴い配列の確保に必要な主記憶容量を抑圧でき、計算機要求を緩和できる。
As described above, before applying the restoration filter, the range cell migration compensation considering aliasing is performed to compensate for the spread in the range direction of the signal, so the data amount to which the restoration filter is applied is reduced to 1 / M. ) Can be reduced.
By reducing the amount of data, the amount of computation for the restoration filter processing can be reduced (to 1 / M), and the processing speed can be increased.
In addition, since the amount of data to which the restoration filter is applied is reduced, the array size required for the composition processing after the restoration filter is applied can be reduced, and the main storage capacity necessary for securing the array can be suppressed accordingly. Can ease the computer requirements.
実施の形態2.
以上の実施の形態1は、2次レンジ圧縮処理をレンジR0を更新しながら適用して高精度な2次レンジ圧縮を行う形態であるが、次に2次レンジ圧縮処理のレンジR0を観測対象領域中央のレンジに固定することで演算量をさらに削減した実施の形態2を示す。
以下、この発明の実施の形態2を図に基づいて説明するが、図において実施の形態1と同一、または、相当する部分については同一記号を付して説明を省略する。
図6は、実施の形態2によるレーダ画像再生装置を示すブロック構成図である。実施の形態1とは、レンジに対するループ処理を表すフィードバックループの戻り位置が変更されている。
The first embodiment described above is a mode in which the secondary range compression processing is applied while updating the range R 0 to perform highly accurate secondary range compression. Next, the range R 0 of the secondary range compression processing is set.
The second embodiment of the present invention will be described below with reference to the drawings. In the figure, the same or corresponding parts as those of the first embodiment are denoted by the same reference numerals and the description thereof is omitted.
FIG. 6 is a block diagram showing a radar image reproducing apparatus according to the second embodiment. In the first embodiment, the return position of the feedback loop representing the loop processing for the range is changed.
次に動作について説明する。
図7は、実施の形態2によるレーダ画像再生装置の処理の流れを示すフローチャートである。
ステップST1_4において、レンジ圧縮部11とレンジ圧縮部12がそれぞれ受信アンテナ1と2のデータを2次元周波数空間で2次レンジ圧縮処理する。この2次レンジ圧縮処理は式8で与えられる関数を乗じることで実現する。
Next, the operation will be described.
FIG. 7 is a flowchart showing a processing flow of the radar image reproducing apparatus according to the second embodiment.
In step ST1_4, the
以上のように、2次レンジ圧縮処理を、後段のレンジに対するループ処理の前にレンジ圧縮処理とともに1回の演算で行っているので、レンジに対するループ処理内での2次レンジ圧縮処理の演算がなくなり、処理を高速化することができる。 As described above, since the secondary range compression processing is performed by one operation together with the range compression processing before the loop processing for the subsequent range, the calculation of the secondary range compression processing within the loop processing for the range is performed. The processing can be speeded up.
この発明は、人工衛星や航空機などの移動プッラットフォームに搭載され地表や海面等の高分解能画像を得る合成開口レーダに適用される。 The present invention is applied to a synthetic aperture radar that is mounted on a moving platform such as an artificial satellite or an aircraft and obtains a high-resolution image of the ground surface or the sea surface.
1;Rx1格納部、2;Rx2格納部、3;出力格納部、11,12;レンジ圧縮部、13,14;RCM補正部、15,16;コヒーレント加算部、17;復元フィルタ部P1(fa;R0)、18;復元フィルタ部P2(fa;R0)、19;合成部、20;アジマス圧縮部。 1; Rx1 storage unit, 2; Rx2 storage unit, 3; Output storage unit, 11, 12; Range compression unit, 13, 14; RCM correction unit, 15, 16; Coherent addition unit, 17; Restoration filter unit P1 (fa R0), 18; restoration filter part P2 (fa; R0), 19; synthesis part, 20; azimuth compression part.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008226965A JP5106323B2 (en) | 2008-09-04 | 2008-09-04 | Radar image playback device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008226965A JP5106323B2 (en) | 2008-09-04 | 2008-09-04 | Radar image playback device |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010060448A true JP2010060448A (en) | 2010-03-18 |
JP5106323B2 JP5106323B2 (en) | 2012-12-26 |
Family
ID=42187391
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008226965A Expired - Fee Related JP5106323B2 (en) | 2008-09-04 | 2008-09-04 | Radar image playback device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5106323B2 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015001428A (en) * | 2013-06-14 | 2015-01-05 | 三菱電機株式会社 | Radar signal processor and radar signal processing method |
WO2016027326A1 (en) * | 2014-08-20 | 2016-02-25 | 三菱電機株式会社 | Signal processing device |
CN105572648A (en) * | 2016-02-01 | 2016-05-11 | 中国科学院电子学研究所 | Synthetic aperture radar echo data range cell migration correction method and device |
KR20190018346A (en) * | 2017-08-14 | 2019-02-22 | 국방과학연구소 | An information processing apparatus and method for processing a reflection signal for a target object |
KR20190018347A (en) * | 2017-08-14 | 2019-02-22 | 국방과학연구소 | An information processing apparatus and method for processing a reflection signal for a target object |
KR20190143307A (en) * | 2018-06-20 | 2019-12-30 | 고려대학교 산학협력단 | Improved azimuth signal reconstruction method in time domain for multichannel sar systems at highy nonuniform sampling |
JP2020016585A (en) * | 2018-07-26 | 2020-01-30 | 日本無線株式会社 | Synthetic aperture radar device, synthetic aperture radar signal processor, and synthetic aperture radar signal processing program |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5840914A (en) * | 1981-09-02 | 1983-03-10 | Nec Corp | Band dividing and synthesizing filter |
JPH06237183A (en) * | 1992-11-30 | 1994-08-23 | American Teleph & Telegr Co <Att> | Decoding method of encoded signal |
JPH0737075A (en) * | 1993-04-05 | 1995-02-07 | General Electric Co <Ge> | Method for generating images having high resolution for radiation imaging system, x-ray imaging system and incident radiation and method for enhancing decomposable space frequency of digital radiation imaging system |
JPH11125673A (en) * | 1997-10-21 | 1999-05-11 | Mitsubishi Electric Corp | Radar device |
JPH11248833A (en) * | 1998-03-04 | 1999-09-17 | Nec Corp | Sar device and its signal-processing method |
JPH11258342A (en) * | 1998-03-10 | 1999-09-24 | Mitsubishi Electric Corp | Radar signal processor |
JP2000207538A (en) * | 1998-11-06 | 2000-07-28 | General Electric Co <Ge> | Method for generating high resolution image of incident radiation and radiation imaging system |
JP2006162518A (en) * | 2004-12-09 | 2006-06-22 | Mitsubishi Electric Corp | Radar system |
JP2007256058A (en) * | 2006-03-23 | 2007-10-04 | Mitsubishi Electric Corp | Radar image processing apparatus |
JP2008203228A (en) * | 2007-02-23 | 2008-09-04 | Nec Corp | Synthetic aperture radar and the synthetic aperture radar image reproduction processing method |
-
2008
- 2008-09-04 JP JP2008226965A patent/JP5106323B2/en not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5840914A (en) * | 1981-09-02 | 1983-03-10 | Nec Corp | Band dividing and synthesizing filter |
JPH06237183A (en) * | 1992-11-30 | 1994-08-23 | American Teleph & Telegr Co <Att> | Decoding method of encoded signal |
JPH0737075A (en) * | 1993-04-05 | 1995-02-07 | General Electric Co <Ge> | Method for generating images having high resolution for radiation imaging system, x-ray imaging system and incident radiation and method for enhancing decomposable space frequency of digital radiation imaging system |
JPH11125673A (en) * | 1997-10-21 | 1999-05-11 | Mitsubishi Electric Corp | Radar device |
JPH11248833A (en) * | 1998-03-04 | 1999-09-17 | Nec Corp | Sar device and its signal-processing method |
JPH11258342A (en) * | 1998-03-10 | 1999-09-24 | Mitsubishi Electric Corp | Radar signal processor |
JP2000207538A (en) * | 1998-11-06 | 2000-07-28 | General Electric Co <Ge> | Method for generating high resolution image of incident radiation and radiation imaging system |
JP2006162518A (en) * | 2004-12-09 | 2006-06-22 | Mitsubishi Electric Corp | Radar system |
JP2007256058A (en) * | 2006-03-23 | 2007-10-04 | Mitsubishi Electric Corp | Radar image processing apparatus |
JP2008203228A (en) * | 2007-02-23 | 2008-09-04 | Nec Corp | Synthetic aperture radar and the synthetic aperture radar image reproduction processing method |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015001428A (en) * | 2013-06-14 | 2015-01-05 | 三菱電機株式会社 | Radar signal processor and radar signal processing method |
WO2016027326A1 (en) * | 2014-08-20 | 2016-02-25 | 三菱電機株式会社 | Signal processing device |
JPWO2016027326A1 (en) * | 2014-08-20 | 2017-04-27 | 三菱電機株式会社 | Signal processing device |
EP3185034A4 (en) * | 2014-08-20 | 2018-03-14 | Mitsubishi Electric Corporation | Signal processing device |
US10754006B2 (en) | 2014-08-20 | 2020-08-25 | Mitsubishi Electric Corporation | Signal processing device |
CN105572648A (en) * | 2016-02-01 | 2016-05-11 | 中国科学院电子学研究所 | Synthetic aperture radar echo data range cell migration correction method and device |
KR20190018347A (en) * | 2017-08-14 | 2019-02-22 | 국방과학연구소 | An information processing apparatus and method for processing a reflection signal for a target object |
KR101996025B1 (en) | 2017-08-14 | 2019-07-03 | 국방과학연구소 | An information processing apparatus and method for processing a reflection signal for a target object |
KR102009402B1 (en) * | 2017-08-14 | 2019-08-09 | 국방과학연구소 | An information processing apparatus and method for processing a reflection signal for a target object |
KR20190018346A (en) * | 2017-08-14 | 2019-02-22 | 국방과학연구소 | An information processing apparatus and method for processing a reflection signal for a target object |
KR20190143307A (en) * | 2018-06-20 | 2019-12-30 | 고려대학교 산학협력단 | Improved azimuth signal reconstruction method in time domain for multichannel sar systems at highy nonuniform sampling |
KR102100905B1 (en) * | 2018-06-20 | 2020-04-14 | 고려대학교 산학협력단 | Improved azimuth signal reconstruction method in time domain for multichannel sar systems at highy nonuniform sampling |
JP2020016585A (en) * | 2018-07-26 | 2020-01-30 | 日本無線株式会社 | Synthetic aperture radar device, synthetic aperture radar signal processor, and synthetic aperture radar signal processing program |
WO2020022468A1 (en) * | 2018-07-26 | 2020-01-30 | 国立大学法人千葉大学 | Synthetic aperture radar device, synthetic aperture radar signal processing device, synthetic aperture radar signal processing program, and synthetic aperture radar observation method |
JP7160268B2 (en) | 2018-07-26 | 2022-10-25 | 日本無線株式会社 | Synthetic aperture radar signal processing device and synthetic aperture radar signal processing program |
Also Published As
Publication number | Publication date |
---|---|
JP5106323B2 (en) | 2012-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5106323B2 (en) | Radar image playback device | |
CN111443339B (en) | Bistatic SAR space-variant correction imaging method, device, equipment and storage medium | |
KR102451696B1 (en) | Phase calibration of a stepped―chirp signal for a synthetic aperature radar | |
JP6066934B2 (en) | System for generating a plurality of SAR images on a satellite or an aerial platform, a satellite comprising the system, and a method for generating a synthetic aperture radar (SAR) image on a satellite or an aerial platform | |
US7397418B1 (en) | SAR image formation with azimuth interpolation after azimuth transform | |
CN104007440B (en) | One accelerated decomposition rear orientation projection spot beam SAR formation method | |
CN110515050B (en) | Satellite-borne SAR real-time echo simulator based on GPU | |
KR101839041B1 (en) | FMCW-SAR system using correction continuous motion effect and method for reconstructing SAR image using FMCW-SAR system | |
JP2011169869A (en) | Apparatus for processing radar signal | |
US7551119B1 (en) | Flight path-driven mitigation of wavefront curvature effects in SAR images | |
JP5219897B2 (en) | Image radar device | |
WO2009102926A1 (en) | Radar imaging of buildings using model-based focusing and data-based focusing | |
JP2012093257A (en) | Radar image processing device | |
JP2006266775A (en) | Radar device | |
JP2015129694A (en) | Synthetic aperture radar device and image processing method thereof | |
JP5901767B2 (en) | Radar system, observation device and data processing device | |
CN111273292A (en) | Synthetic aperture radar high-frequency vibration compensation method and device, electronic equipment and medium | |
JP2007225442A (en) | Radar device | |
JP2010066131A (en) | Image radar device and image radar system | |
CN109343056A (en) | A kind of RD imaging method and device for nonlinear frequency modulation SAR | |
CN117129994A (en) | Improved backward projection imaging method based on phase compensation nuclear GNSS-SAR | |
Renga et al. | Multi-Platform Image Synthesis for Distributed Synthetic Aperture Radar in Long Baseline Bistatic Configurations | |
JPWO2020136871A1 (en) | Radar signal processing device and radar signal processing method | |
CN104749575B (en) | A kind of improved geostationary orbit SAR frequency domain imaging methods | |
US7876256B2 (en) | Antenna back-lobe rejection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110609 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120817 |
|
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: 20120904 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20121002 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20151012 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |