JPS6368980A - Image reproduction processing system for synthetic aperture radar - Google Patents

Image reproduction processing system for synthetic aperture radar

Info

Publication number
JPS6368980A
JPS6368980A JP61211589A JP21158986A JPS6368980A JP S6368980 A JPS6368980 A JP S6368980A JP 61211589 A JP61211589 A JP 61211589A JP 21158986 A JP21158986 A JP 21158986A JP S6368980 A JPS6368980 A JP S6368980A
Authority
JP
Japan
Prior art keywords
image
resampling
range
dimensional
rotation factor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP61211589A
Other languages
Japanese (ja)
Other versions
JPH07120405B2 (en
Inventor
Akira Maeda
章 前田
Fuminobu Furumura
文伸 古村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP61211589A priority Critical patent/JPH07120405B2/en
Publication of JPS6368980A publication Critical patent/JPS6368980A/en
Publication of JPH07120405B2 publication Critical patent/JPH07120405B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Image Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Image Analysis (AREA)

Abstract

PURPOSE:To attain reproduction high in picture quality by making one- dimensional sampling of reproduced complex image in azimuthal direction, correcting phase deviation in range direction multiplying a phase rotation factor and making one-dimensional resampling in range direction. CONSTITUTION:A synthetic aperture radar SAR provided with a radar sensor 1 and an antenna 2 is mounted on an artificial satellite etc., and photographs ground surface while moving on the flying path 3 in the direction of arrow mark 4. Photographing data from the radar SAR are received by a ground station 6. In such a case, phase deviation in azimuthal direction is corrected by multiplying complex image after azimuth compression by phase rotation factor exp(-2pijBt). Then, one-dimensional resampling is made in azimuthal direction. Image data resampled in azimuthal direction is multiplied by phase rotation factor exp(2pijBt) again, and phase deviation in range direction is corrected. Finally, one-dimensional resampling is made in range direction to complete two-dimensional resampling process, and a reproduced image is obtained.

Description

【発明の詳細な説明】 本発明は人工衛星あるいは航空機等に搭載される合成開
口レーダ(Synthetic Aperture R
adar。
DETAILED DESCRIPTION OF THE INVENTION The present invention is directed to a synthetic aperture radar (Synthetic Aperture R) mounted on an artificial satellite or an aircraft.
adar.

以下rSARJという)による撮像データから、人間が
理解できる画像を再生するためのディジタル処理システ
ムに係り、特に高画質を画像を再生するのに好適なSA
R画像再生処理方式に関する。
This relates to a digital processing system for reproducing human-understandable images from imaging data obtained by rSARJ (hereinafter referred to as rSARJ), and is particularly suitable for reproducing images with high image quality.
This invention relates to an R image reproduction processing method.

〔従来の技術〕[Conventional technology]

計算機によるディジタル処理でSAR画像を再生する方
式については、例えば「環境のリモートセンシングに関
する第13回国際シンポジウム予稿集337〜360頁
(1977年)におけるベネット及びカミングによる6
′ディジタルSAR画像生成:エアボーン及びサテライ
ト結果′″ (Bennette 。
Regarding the method of reproducing SAR images through digital processing by a computer, see, for example, ``6 by Bennett and Cumming in Proceedings of the 13th International Symposium on Environmental Remote Sensing, pp. 337-360 (1977).
``Digital SAR Image Generation: Airborne and Satellite Results'' (Bennette).

CCumm1n : ”Digital SARIma
ge Formation :Ajrbnne and
 5atellite Re5ults” 、 13t
hProceedings of Internati
onal Symposium onRemote S
ensing of Environment、 pp
337−360(1977) ) Jと題する文献にお
いて詳しく述べられている。また上記文献中には、SA
R再生画像の支間分解能を向上し、高画質な画像を再生
する方式としてアジマス圧縮後の複素画像に対してリサ
ンプリングを行ない、画像を拡大しておく方式%式% 〔発明が解決しようとする問題点〕 上記従来技術はアジマス圧縮後の複素画像の位相ずれを
考慮していないため上記サンプリングにより高周波雑音
成分の増加、縞状アーチファクトの発生などの端劣化を
招く。以下、これを説明する。
CCumm1n: ”Digital SARIma
ge Formation: Ajrbnne and
5atellite Re5ults”, 13t
hProceedings of International
onal Symposium onRemote S
Ensuring of Environment, pp.
337-360 (1977)) J. Also, in the above literature, SA
A method of enlarging the image by performing resampling on the complex image after azimuth compression as a method of improving the interval resolution of the R reproduced image and reproducing a high-quality image. Problem] Since the above-mentioned conventional technology does not take into account the phase shift of the complex image after azimuth compression, the above-mentioned sampling causes edge deterioration such as an increase in high-frequency noise components and the generation of striped artifacts. This will be explained below.

レンジ圧縮後のSAR画像においては、地表の各点の情
報はアジマス方向に1次元的に拡がって分布する。この
拡がりを表わす関数をアジマス方向点像パターンと呼び h(t) =exp[2πj (αt”+βt)]  
 −(1)と書ける。ここで、tはSARセンサのアジ
マス方向座標を表わし、jは虚数単位、α、βは再生パ
ラメータと呼ばれる量で、αをドツプラ変化率βをドツ
プラ中心周波数と呼ぶ。α、βは衛星と冷上目標点との
距離や相対速度に依存するため、レンジ方向にゆるやか
に変化する。したがってアジマス方向点像パターンは1
シーン中で一定ではない。
In the SAR image after range compression, information on each point on the earth's surface is distributed one-dimensionally in the azimuth direction. The function representing this spread is called the azimuth direction point spread pattern h(t) = exp[2πj (αt''+βt)]
- It can be written as (1). Here, t represents the coordinate in the azimuth direction of the SAR sensor, j is an imaginary unit, α and β are quantities called reproduction parameters, and α is the Doppler change rate β is called the Doppler center frequency. Since α and β depend on the distance and relative speed between the satellite and the cold target point, they change slowly in the range direction. Therefore, the point spread pattern in the azimuth direction is 1
Not constant throughout the scene.

実際のSAR画像においては、レンジ圧縮後の画像中の
アジマス方向点像パターンは、SARセンサと地表の各
点との相対距離変化に起因するレンジカーバチャ歪と呼
ばれる幾何学的歪を受けており、アジマス方向の直線上
ではなく、曲線上に分布する。以下ではこの歪はレンジ
カバーチャ補正により既に補正されていると考えるもの
とする。
In an actual SAR image, the azimuthal point image pattern in the image after range compression is subject to geometric distortion called range curvature distortion, which is caused by changes in the relative distance between the SAR sensor and each point on the ground surface. , is distributed not on a straight line in the azimuth direction but on a curve. In the following, it is assumed that this distortion has already been corrected by range coverage correction.

アジマス方向点像パターン(1)は、周波数領域でドツ
プラ中心周波数βを中心に、ある幅をもって分布する。
The azimuth direction point spread pattern (1) is distributed with a certain width around the Doppler center frequency β in the frequency domain.

したがって、レンジ圧縮後のSAR画像とアジマス方向
点像パターンの相関をとってアジマス圧縮した結果のデ
ータも、周波数領域上でドツプラ中心周波数βを中心に
分布する。すなわちアジマス圧縮後の複素画像において
は、アジマス方向に隣接する画素の間に位相のずれが存
在し、この複素画像に対してキュービックコンボリュー
ションなどの補間を行う従来方式では、空間分解能の劣
化を防ぐ目的に反して逆に画質の劣化を引き起す事にな
る。
Therefore, the data resulting from azimuth compression by correlating the SAR image after range compression with the azimuth direction point image pattern is also distributed around the Doppler center frequency β in the frequency domain. In other words, in a complex image after azimuth compression, there is a phase shift between adjacent pixels in the azimuth direction, and conventional methods that perform interpolation such as cubic convolution on this complex image prevent deterioration of spatial resolution. Contrary to the purpose, this will cause deterioration of image quality.

一方、周波数領域における平行移動または実時間領域に
おける位相回転因子の乗算により、上記位相ずれを補正
する方式も考えられる。この方式ではアジマス圧縮にお
いて、各ライン毎にexp(−2πjβt〕なる位相回
転因子を乗算することと等価であるが、βがレンジ方向
に変化する量であるため、ライン毎に異なる位相回転因
子が乗算される結果、レンジ方向に隣接する画素間に位
相ずれが生じる。
On the other hand, a method of correcting the above-mentioned phase shift by parallel movement in the frequency domain or multiplication by a phase rotation factor in the real time domain is also conceivable. This method is equivalent to multiplying each line by a phase rotation factor of exp (-2πjβt) in azimuth compression, but since β is the amount that changes in the range direction, a different phase rotation factor is applied to each line. As a result of the multiplication, a phase shift occurs between pixels adjacent in the range direction.

以上述べた様に、アジマス圧縮後の複素画像にはアジマ
スまたはレンジ方向のいずれかの方向に位相ずれの発生
が避けられず、この位相ずれを考慮せずにサンプリング
を行なうと、前述した如く画質劣化が生ずるという問題
点があった。
As mentioned above, it is inevitable that a complex image after azimuth compression will have a phase shift in either the azimuth or range direction, and if sampling is performed without taking this phase shift into consideration, the image quality will deteriorate as described above. There was a problem that deterioration occurred.

本発明では以上の如き問題点を解決し、両画質のSAR
画像を再生する処理方式を提供する事にある。
The present invention solves the above problems and provides SAR of both image quality.
The objective is to provide a processing method for reproducing images.

〔問題点を解決するための手段〕[Means for solving problems]

上記目的は、アジマス圧縮後のSAR複素画像に対する
2次元リサンプリングを、レンジ方向とアジマス方向の
各々1次元リサンプリングに分割して行ない、レンジ方
向またはアジマス方向いずれかの方向に1次元リサンプ
リングを行なった後の複素画像に対して位相回転因子を
乗算して位相ずれを補正する事により達成される。
The above purpose is to perform two-dimensional resampling on the SAR complex image after azimuth compression by dividing it into one-dimensional resampling in the range direction and one-dimensional resampling in the azimuth direction, and perform one-dimensional resampling in either the range direction or the azimuth direction. This is achieved by multiplying the complex image by a phase rotation factor to correct the phase shift.

〔作用〕[Effect]

一般に2次元リサンプリングはたて・よこ各々1次元リ
サンプリングに分割できる。したがってSARの場合、
まずアジマス方向のみの1次元リサンプリング、次いで
レンジ方向のみの1次元リサンプリングに分割する事が
できる。したがってアジマス方向1次元リサンプリング
の際には再生複素画像のアジマス方向の位相ずれを補正
しておき、次いでアジマス方向リサンプリング後の画像
に位相回転因子を乗算してレンジ方向の位相ずれを補正
しておき、しかる後にレンジ方向1次元リサンプリング
を行なう。これによってレンジ・アジマス各方向の1次
元リサンプリングは、つねに位相ずれの補正されたデー
タに対してなされるため、リサンプリングによる画質の
劣化を防ぐ事ができる。また、レンジ方向1次元リサン
プリングから先に行なう場合も同様である。
Generally, two-dimensional resampling can be divided into vertical and horizontal one-dimensional resampling. Therefore, in the case of SAR,
First, it can be divided into one-dimensional resampling only in the azimuth direction, and then one-dimensional resampling only in the range direction. Therefore, when performing one-dimensional resampling in the azimuth direction, the phase shift in the azimuth direction of the reproduced complex image is corrected, and then the phase shift in the range direction is corrected by multiplying the image after the azimuth direction resampling by a phase rotation factor. After that, one-dimensional resampling in the range direction is performed. As a result, one-dimensional resampling in each direction of range and azimuth is always performed on data whose phase shift has been corrected, so that deterioration in image quality due to resampling can be prevented. The same applies to the case where one-dimensional resampling in the range direction is performed first.

〔実施例〕〔Example〕

以下、本発明の一実施例を第1図および第2図により説
明する。
An embodiment of the present invention will be described below with reference to FIGS. 1 and 2.

第2図にSARの全体システムを示した。レーダ・セン
サ1.アンテナ2を有するSARは1人工衛星等に搭載
されて飛行経路3上を矢印4方向に移動しつつ表地の撮
像を行う。SARからの撮像データは地上局6で受信さ
れ、データ・プロセッサ6により処理されて映像フィル
ム7の作成、データ記憶用磁気テープ8の作成等が行わ
れる。
Figure 2 shows the entire SAR system. Radar sensor 1. A SAR having an antenna 2 is mounted on an artificial satellite or the like and images the surface while moving on a flight path 3 in the direction of an arrow 4. Imaging data from the SAR is received by a ground station 6 and processed by a data processor 6 to create a video film 7, a magnetic tape 8 for data storage, and the like.

なお、9は分解能セルを、10はSARで採取されるデ
ータの地表上のレンジ方向を、11は同アジマス方向を
、12はアンテナ・ビームをそして13は刈り幅をそれ
ぞれ示している。
Note that 9 indicates a resolution cell, 10 indicates a range direction on the ground surface of data collected by SAR, 11 indicates an azimuth direction, 12 indicates an antenna beam, and 13 indicates a cutting width.

第1−図は発明方式によるSAR画像再生処理中の複素
画像リサンプリング処理のフローチャートである。
FIG. 1 is a flowchart of complex image resampling processing during SAR image reproduction processing according to the invention method.

アジマス圧縮処理直後のSAR複素画像は、レンジ方向
に位相が揃っているものとする。このとき前述した様に
、アジマス方向には exp[2πjβt〕だけの位相ずれが存在する。
It is assumed that the SAR complex image immediately after the azimuth compression process has the same phase in the range direction. At this time, as described above, there is a phase shift of exp[2πjβt] in the azimuth direction.

まずステップ21において、アジマス圧縮後の複素画像
25に対してexp(−2πjβt〕なる位相回転因子
を乗算し、アジマス方向の位相ずれを補正する。次いで
ステップ22においてアジマス方向1次元リサンプリン
グを行なう。次にステンプ23においてアジマス方向に
リサンプルした画像データに再びexp[2πjβt〕
なる位相回転因子を乗算し、レンジ方向の位相ずれを補
正する。
First, in step 21, the complex image 25 after azimuth compression is multiplied by a phase rotation factor exp (-2πjβt) to correct the phase shift in the azimuth direction.Next, in step 22, one-dimensional resampling in the azimuth direction is performed. Next, in the step 23, exp[2πjβt] is again applied to the image data resampled in the azimuth direction.
The phase shift in the range direction is corrected by multiplying by a phase rotation factor.

最後にステップ24においてレンジ方向1次元リサンプ
リング行なう事により2次元リサンプリング処理を終了
し再生画像26を得る。
Finally, in step 24, one-dimensional resampling in the range direction is performed to complete the two-dimensional resampling process and obtain a reproduced image 26.

以上、本実施例ではアジマス方向リサンプリングを先に
実行するものとしたが、先にレンジ方向リサンプリング
を行なう場合も全く同様に、レンジ方向リサンプリング
後の画像に位相回転因子を −乗算し、アジマス方向の
位置ずれを補正する事により同様の効果が得られる事は
言うまでもない。
As described above, in this embodiment, the azimuth direction resampling is performed first, but when performing the range direction resampling first, the image after the range direction resampling is -multiplied by the phase rotation factor, It goes without saying that the same effect can be obtained by correcting the positional deviation in the azimuth direction.

〔発明の効果〕〔Effect of the invention〕

以上、本発明によれば、常に位相ずれの補正が行なわれ
たデータに対してリサンプリング処理がなされるため、
位相ずれによる画質の劣化を防ぐ事ができ、SAR画像
の高画質再生が可能になるという効果がある。
As described above, according to the present invention, resampling processing is always performed on data whose phase shift has been corrected.
This has the effect of preventing deterioration of image quality due to phase shift and enabling high-quality reproduction of SAR images.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は本発明の一実施例によるSAR複素画像のりサ
ンプリング処理のフロー図、第2図は第 / 囮
FIG. 1 is a flow diagram of SAR complex image sampling processing according to an embodiment of the present invention, and FIG.

Claims (1)

【特許請求の範囲】[Claims] 1、合成開口レーダによる撮像データから画像を再生す
る画像処理システムにおいて、再生された複素画像に対
してレンジまたはアジマス方向に1次元リサンプリング
を行なつた画像に位相回転因子を乗算する事を特徴とす
る合成開口レーダの画像再生処理方式。
1. In an image processing system that reproduces an image from imaging data obtained by synthetic aperture radar, a feature is that the reproduced complex image is subjected to one-dimensional resampling in the range or azimuth direction and the image is multiplied by a phase rotation factor. Image reproduction processing method for synthetic aperture radar.
JP61211589A 1986-09-10 1986-09-10 Image reproduction processing method for synthetic aperture radar Expired - Lifetime JPH07120405B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP61211589A JPH07120405B2 (en) 1986-09-10 1986-09-10 Image reproduction processing method for synthetic aperture radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP61211589A JPH07120405B2 (en) 1986-09-10 1986-09-10 Image reproduction processing method for synthetic aperture radar

Publications (2)

Publication Number Publication Date
JPS6368980A true JPS6368980A (en) 1988-03-28
JPH07120405B2 JPH07120405B2 (en) 1995-12-20

Family

ID=16608267

Family Applications (1)

Application Number Title Priority Date Filing Date
JP61211589A Expired - Lifetime JPH07120405B2 (en) 1986-09-10 1986-09-10 Image reproduction processing method for synthetic aperture radar

Country Status (1)

Country Link
JP (1) JPH07120405B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08152469A (en) * 1994-11-29 1996-06-11 Toshiba Corp Radar apparatus
WO2015008554A1 (en) * 2013-07-19 2015-01-22 国立大学法人東北大学 Synthetic aperture processing sensor, processing method for sensor, and program
CN105629232A (en) * 2015-12-24 2016-06-01 中国科学院电子学研究所 Arrangement method of multi-ground calibration station used for GEOSAR phase calibration

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08152469A (en) * 1994-11-29 1996-06-11 Toshiba Corp Radar apparatus
WO2015008554A1 (en) * 2013-07-19 2015-01-22 国立大学法人東北大学 Synthetic aperture processing sensor, processing method for sensor, and program
JPWO2015008554A1 (en) * 2013-07-19 2017-03-02 国立大学法人東北大学 Sensor with synthetic aperture processing, method for processing the sensor, and program
CN105629232A (en) * 2015-12-24 2016-06-01 中国科学院电子学研究所 Arrangement method of multi-ground calibration station used for GEOSAR phase calibration

Also Published As

Publication number Publication date
JPH07120405B2 (en) 1995-12-20

Similar Documents

Publication Publication Date Title
US4616227A (en) Method of reconstructing synthetic aperture radar image
Rifman Digital rectification of ERTS multispectral imagery
Ran et al. An autofocus algorithm for estimating residual trajectory deviations in synthetic aperture radar
CN108318880B (en) Polar coordinate SAR imaging method with parametric motion error
CN108205135B (en) Radar video imaging method based on non-interpolation fusion fast backward projection
JPH0646405B2 (en) Synthetic aperture radar image reproduction processing method
Xiao et al. An N 2 logN back-projection algorithm for SAR image formation
Jia et al. Simplified real‐time imaging flow for high‐resolution FMCW SAR
CN109188432B (en) Rapid BP imaging method of parallel bistatic bunching SAR
Qiu et al. An Omega-K algorithm with phase error compensation for bistatic SAR of a translational invariant case
CN114910911A (en) Satellite-borne multi-base SAR imaging method based on multi-phase center reconstruction
US9146312B1 (en) Pre-processing SAR image stream to facilitate compression for transport on bandwidth-limited-link
CN111220979B (en) Imaging method for curved synthetic aperture radar
CN114325704A (en) Rapid time domain imaging method of synthetic aperture radar based on wave number spectrum splicing
JP2000276045A (en) Method and device for making map using photographed picture, and method of correcting distortion of photographed picture
JPS6368980A (en) Image reproduction processing system for synthetic aperture radar
CN112558070A (en) Frequency domain imaging method and device of circular scanning foundation SAR
JPS6280767A (en) Reproducing processing system for synthetic aperture radar image
CN114384519B (en) Ultrahigh-resolution satellite-borne synthetic aperture radar imaging method and device
CN115685202A (en) Baseline estimation method and device of space-borne distributed interferometric synthetic aperture radar
JPS61187079A (en) Reproducing processing system of synthetic aperture radar picture image
CN112305539B (en) ArcSAR polar coordinate format imaging method based on spherical wave decomposition
JPH0799392B2 (en) Synthetic aperture radar image reproduction processing method
KR102070974B1 (en) Method and apparatus for correcting geometric distortion
CN113552564B (en) SAR time domain rapid imaging method, system, terminal and application for complex terrain scene

Legal Events

Date Code Title Description
EXPY Cancellation because of completion of term