JPS6165366A - 合成開口レ−ダ画像再生処理方法 - Google Patents

合成開口レ−ダ画像再生処理方法

Info

Publication number
JPS6165366A
JPS6165366A JP59186317A JP18631784A JPS6165366A JP S6165366 A JPS6165366 A JP S6165366A JP 59186317 A JP59186317 A JP 59186317A JP 18631784 A JP18631784 A JP 18631784A JP S6165366 A JPS6165366 A JP S6165366A
Authority
JP
Japan
Prior art keywords
image
processing
executed
synthetic
apperture
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
JP59186317A
Other languages
English (en)
Other versions
JPH0646405B2 (ja
Inventor
Akira Maeda
章 前田
Akira Tsuboi
坪井 晃
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 JP59186317A priority Critical patent/JPH0646405B2/ja
Priority to CA000489007A priority patent/CA1257370A/en
Priority to US06/767,046 priority patent/US4758838A/en
Publication of JPS6165366A publication Critical patent/JPS6165366A/ja
Publication of JPH0646405B2 publication Critical patent/JPH0646405B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004SAR image acquisition techniques
    • G01S13/9011SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Processing Or Creating Images (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔発明の利用分野〕 本発明は人工衛星あるいは航空機等に搭載される合成開
口レーダ(Synthetic Aperture R
adar。
以下rsARJという)による撮像データから、人間が
理解できる画像を再生するためのディジタル処理システ
ムに係り、特に高画質画像を高速に再生するのに好適な
画像再生処理方式に関する。
〔発明の背景〕
人工衛星あるいは航空機等を用いたりモートセンシング
の分野では、地表を撮像するためのセンサとして、雲を
透過するマイクロ波帯で高解像の画像が得られるSAR
が注目されている。
第1図にSARの全体システムを示した。レーダ・セン
サ1.アンテナ2を有するSARは、人工衛星等に搭載
されて飛行経路3−ヒを矢印4方向に移動しつつ地表の
撮像を行う。SARからの撮像データは地上局5で受信
され、データ・プロセッサ6により処理されて映像フィ
ルム7の作成、データ記憶用磁気テーブル8の作成等が
行われる。
なお、9は分解セルを、10はSARで採取されるデー
タの地表上のレンジ方向を、11は同アジマス方向を、
12はアンテナ・ビームをそして13は刈り幅をそれぞ
れ示している。
以下、SARで採取されたデータの処理の概要を述べる
。なお、詳補については、例えば〔安藤″合成開ロレー
ダと間接計測技術″計測と制御。
Vo12.22.翫2〕 を参照されたい。
SARの受信画像中においては、原画像上の1点が点像
パターンh(x+y)の広がりをもって分布しており、
このままでは利用できない。ここでXは前記レンジ方向
を、yは前記アジマス方向を示している。前記受信画像
中で広がっている情報は、まずレンジ方向に圧縮され、
次にアジマス方向に圧縮される。この様子を第2図に示
す。第2図のAは地表にマイクロ波反射点が2点だけ存
在したときの受信画像を模式的に示したものであるが、
2方向に圧縮処理を行えば、Bのようにもとの地表パタ
ーンを得ることができる。前記レンジ圧縮処理は画像デ
ータ1ラインごとの点像パターンデータとのコリレーシ
ョン処理によって行う。
但し、コリレーション処理をそのまま実行すると、莫大
な処理時間がかかるため、高速フーリエ変換(以下、r
FFTJという)、複素乗算、高速逆フーリエ変換(以
下、rIFFTJという。)を用いて高速化が図られる
。FFTを用いてコリレーション処理を行なうには、ま
ず点像パターンを計算機によるディジタル処理で生成し
、点像パターンと画像データ1ラインの両方のFFTを
計算する。2個のデータのコリレーションはFFTを計
算した後の周波数領域では単なる乗算になることから、
上記2個のデータのFFT計算結果の積をとり、それを
IFFTすることにより、1ライン分のコリレーション
結果が得られる。
レンジ圧縮処理における点像パターンは、SARセンサ
システムに固有なものであり、少なくとも1シーンの中
では全く変化しないとみなせるので、上記処理における
点像パターンの生成とそのFETは1シーンの処理で高
々1回行なえばよく、その処理時間のSAR再生処理時
間全体に占める割合は無視できる程度に小さい。
次いで行なわれるアジマス圧縮処理も同様にFFT使用
した高速化がなされるが、アジマス圧線処理における点
像パターンは、SARセンサと目標との距離や相対速度
に依存するため、1ラインごとに異なったものとなる。
したがって従来の方法ではアジマス圧縮処理において、
アジマス方向点像パターンの生成、FFTを毎ライン繰
り返さねばならずその演算量はSAR再生処理全体の演
算量の約20%にも達し、SAR再生処理全体の演算量
が増大する、といった欠点があった。
〔発明の目的〕
本発明の目的は、上記の如き欠点を解消し、SAR再生
画像の画質を劣化させることなく、SAR再生処理の演
算量を減少させ、高画像のSAR画像を高速に再生する
処理方式を提供することにある。
〔発明の概要〕
オ 以下では、アジマス吟向圧縮処理に限定して説明する。
4を アジマス方向点線パターンは次の様に表わされる。
h、=exp (2tc j Ca n”+βn)) 
  −(1)ユニでnは整数で、サンプル番号を表わす
。α。
βは再生パラメータと呼ばれる量で、αをドツプラ変化
率、βをドツプラ中心周波数と呼ぶ。α。
βは衛星と目標との距離や相対速度に依存するため、ラ
イン方向(レンジ方向)にゆるやかに変化する。
したがって、本発明では、α、βというパラメータがゆ
るやかに変化することから各ライン毎に式(1)に従っ
て点像パターンを生成、FFTするのではなく、数ライ
ンに1回に省略することにより、演算量を減少させる点
に特徴がある。
すなわち、本発明方式では、ゆるやかに変化するα、β
パラメータを階段状に近似することに相当する。第3図
は、この近似を説明するもので曲線31が正確なαまた
はβの値、階段状の折れ線32が近似した結果である。
ところが、このままだと、βの値が不連続的に変化する
場所で再生画像に不連続な位置ずれ、位相ずれが生じ、
再生画像の画質を劣化させる原因となる。以下、この現
象を説明する。
地表に電波を反射する点が1点だけ存在するとき、ある
ラインにおけるSARの観測データをg、とすると、 g、=exp (2πj ((E n2+βn))  
−(2)となる。一般の場合は同様のパターンの線型重
ね合わせになるが、ここでは簡単のため地表には1夜 点しか存柑しないものとして説明する。式(2)と式(
1)の点線パターンとのコリレーションをとって圧縮し
た結果をf、とすると !、=ΣFg、、hhh” =Σexp (2πj (a (n+k)”+β(n+
k)−rt k2−βk))=A(n)exp (2π
j(rzn”+βn))−−(3)ただし A(n)=Σexp (2π、i(2αnk))   
”・(4)である。
次に再生パラメータβに誤差δβを含む場合、点像パタ
ーンは、 h、=exp (2tc j (a n”+(β+δβ
)n)) −(5)となる。式(5)の点像パターンで
式(2)の観測データg (n)を圧縮した結果をり、
とするとfk=Σg。ek)lk” =Σexp (2tt j ((2an−δβ)k、+
αn”十βn)]2α 式(6)が示すように、圧縮された点像は、再生パラメ
ータβに含まれる誤差δβに比例した位置ずれを生じる
この位置ずれを補正するため、本発明方式では周波数領
域で位相回転因子を乗算する。アジマス圧縮におけるI
FFTは と書ける。ここでFkは点像パターンと観測データをF
FTにした結果を乗じたもの、Nはデータ点数である。
乗じIFFTL/た結果をJ、′とすれば=f1.◆ 
         ・・・・・・ (8)となる。これ
はf、に位置ずれφを加えたものである。
したがって、周波数領域で適当な位相回転因子を乗する
ことにより、(6)式の位置ずれを補正できる。その結
果をf、とすると 2α 生じていることが分かる。位置ずれ同様、この位相ずれ
も再生画像画質劣化の要因となるので補正する必要があ
る。ただしこの位相ずれはサンプル番号nに依存しない
ライン毎の定数なので、(6)式の位置ずれ補正の際に
、同時に補正することが可能である。
したがって周波数領域で ターン生成を省略したことによる位置ずれ、位相ずれ共
に補正され、再生画像の画質劣化を防ぐ事イjζ二 ができる。また、位相ずれの補正は装置ずれ補正とは別
に、実時間領域における位相回転因子の乗算により補正
することも可能である。
〔発明の実施例〕
以下、本発明の一実施例を第4図、第5図により説明す
る。
第4図はSAR画像再生処理の全体フローである。
観測データはまずレンジ圧縮処理4]によりレンジ方向
に圧縮される。次にアジマス圧縮に先立ち衛星の軌道・
姿勢パラメータから再生パラメータを求める再生パラメ
ータ計算処理42を行ない、そこで求めた再生パラメー
タを用いてアジマス圧縮処理43を行ない、アジマス方
向に圧縮する。
次に圧縮された結果の歪を補正する歪補正リサンプリン
グ44を行なう。このリサンプリングは、再生画像の周
波数帯域を考慮して複素数のデータに対して行なう。そ
の後、ノルム化処理45を施すことによりデータを実数
化し、雑音低減処理であるマルチルツクメ加算処理46
を行なった後、量子化処理47を行ない最終画像を得る
第5図は、第4図のフローで本発明に関わるアジマス圧
縮処理43の1ラインに対する処理を詳細に示したもの
である。
まずレンジ圧縮後の観測データにFFT51をジ剰 施し、周祈数領域に変換した後、衛星と目標の間の距離
変化を補正するレンジカーパチャ補正52を行な嵐。処
理53では、点像パターンの生成を必要とするか否かの
判定を行う。例えば、16ラインに1回だけ点像パター
ンを生成するものとすれば、処理53では現在処理中の
ライン番号を16で割った剰余を計算し、その値で判定
を行なう。ここで点像パターン生成が必要と判定された
場合は点像パターン生成処理54、そのFF7本 55が実行され、そ労以外の場合は以前に生成されたデ
ータを使用して処理56の実行に移る。処理56ではレ
ンジカーバチャ補正52の結果と、FFTされた点像パ
ターンの複素乗算が行なわれる。
次いで、処理57において、現在処理中のラインに対応
する再生パラメータα、βと、点像パターン生成に使用
した再生パラメータβ′から、位値ずれ、位相ずれを補
正する位相回転因子−β、Nはデータ点数)を式(10
)に従って生成し、処理56の結果に乗する。最後にI
FFT58を施して、アジマス圧縮結果を得る・以上述
べた様に、本実施例によれば、点像パタできる、という
効果がある。
〔発明の効果〕
本発明によれば、SAR画像再生処理においてアジマス
方向点像パターン生成とそのFFTを再生画像の画質を
劣化させることなく大部分省略することができるので、
SAR画像再生処理の演算量を減らし、高速化できると
いう効果がある。
【図面の簡単な説明】
第1図はSARシステムの全体構成図、第2図はSAR
画像再生の原理を示す図、第3図は再生パラメータの近
似法を説明する図、第4図は5ARPii像再生処理像
再生処理−チャート、第5図は、第4図の処理43の詳
細なフローチャートである。 1・・・レーダ・センサ、2・・・アンテナ、3・・・
飛行経路、5・・・地上局、6・・・データ・プロセッ
サ、10V2図 竿牛目 多5 目

Claims (1)

  1. 【特許請求の範囲】 1、合成開口レーダによる撮像データから原画像を再生
    する画像再生処理システムにおいて、原画像上の各点に
    対応して点像パータンを生成し、生成された点像パター
    ンと1ライン分の撮像データとを高速フーリエ変換する
    処理を複数ラインおきにおこなうことを特徴とする合成
    開口レーダ画像再生処理方式。 2、上記点像関数はドップラ中心周波数を再生パラメー
    タとして含み、該ドップラ中心周波数に含まれる誤差に
    もとづく再生画像上の位置ずれと位相ずれを、高速フー
    リエ変換された点像パターンと1ライン分の画像データ
    との積に位相回転因子を乗算して補正することを特徴と
    する第1項の合成開口レーダ画像再生処理方式。
JP59186317A 1984-09-07 1984-09-07 合成開口レ−ダ画像再生処理方法 Expired - Lifetime JPH0646405B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP59186317A JPH0646405B2 (ja) 1984-09-07 1984-09-07 合成開口レ−ダ画像再生処理方法
CA000489007A CA1257370A (en) 1984-09-07 1985-08-19 Method of reconstructing images from synthetic aperture radar's data
US06/767,046 US4758838A (en) 1984-09-07 1985-08-19 Method of reconstructing images from synthetic aperture radar's data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP59186317A JPH0646405B2 (ja) 1984-09-07 1984-09-07 合成開口レ−ダ画像再生処理方法

Publications (2)

Publication Number Publication Date
JPS6165366A true JPS6165366A (ja) 1986-04-03
JPH0646405B2 JPH0646405B2 (ja) 1994-06-15

Family

ID=16186213

Family Applications (1)

Application Number Title Priority Date Filing Date
JP59186317A Expired - Lifetime JPH0646405B2 (ja) 1984-09-07 1984-09-07 合成開口レ−ダ画像再生処理方法

Country Status (3)

Country Link
US (1) US4758838A (ja)
JP (1) JPH0646405B2 (ja)
CA (1) CA1257370A (ja)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3688528T2 (de) * 1986-01-22 1994-01-13 Hughes Aircraft Co Optische analoge datenverarbeitungsanordnungen zur behandlung von bipolaren und komplexen daten.
US4829303A (en) * 1988-05-18 1989-05-09 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Data volume reduction for imaging radar polarimetry
IT1231358B (it) * 1989-04-21 1991-12-02 Selenia Ind Elettroniche Dispositivo per migliorare la risoluzione radar
US5053778A (en) * 1989-08-10 1991-10-01 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Generation of topographic terrain models utilizing synthetic aperture radar and surface level data
FR2671879B1 (fr) * 1991-01-22 1993-11-19 Alcatel Espace Dispositif, embarqu2 sur satellite, de mesure du coefficient de retrodiffusion de la mer.
DE4122592C1 (ja) * 1991-07-08 1993-01-28 Deutsche Forschungsanstalt Fuer Luft- Und Raumfahrt Ev, 5300 Bonn, De
US5608404A (en) * 1993-06-23 1997-03-04 The United States Of America As Represented By The United States Department Of Energy Imaging synthetic aperture radar
JP2757845B2 (ja) * 1995-12-26 1998-05-25 日本電気株式会社 合成開口レーダ装置
WO1998002760A1 (en) * 1996-06-28 1998-01-22 Milkovich Systems Engineering Improved fast fourier transform sensor processor using enhanced frequency domain principles
US5943006A (en) * 1996-10-25 1999-08-24 Patrick E. Crane RF image reconstruction and super resolution using fourier transform techniques
WO1999031525A1 (en) * 1997-12-15 1999-06-24 Milkovich Systems Engineering Signal processing architecture which improves sonar and pulse doppler radar performance and tracking capability
US6864828B1 (en) * 2003-02-18 2005-03-08 Lockheed Martin Corporation Method and apparatus for collection and processing of interferometric synthetic aperture radar data
US7307580B2 (en) 2006-01-17 2007-12-11 Raytheon Company Non-statistical method for compressing and decompressing complex SAR data
US7551119B1 (en) * 2008-01-08 2009-06-23 Sandia Corporation Flight path-driven mitigation of wavefront curvature effects in SAR images
US8327838B1 (en) * 2008-10-16 2012-12-11 Lockheed Martin Corporation Solar parabolic trough mirror/receiver alignment
CN103558601B (zh) * 2013-11-01 2016-04-20 北京航空航天大学 一种用于合成孔径雷达过采样数据处理的超分辨方法
CN103645476B (zh) * 2013-12-18 2015-11-04 中国国土资源航空物探遥感中心 一种合成孔径雷达差分干涉图序列的时空同质滤波方法
CN104569945A (zh) * 2014-12-25 2015-04-29 北京无线电计量测试研究所 一种用于合成孔径雷达图像的自适应滤波方法
US9887863B2 (en) * 2015-10-16 2018-02-06 Taiwan Semiconductor Manufacturing Company Ltd. Transceiver group and associated router
CN105353375A (zh) * 2015-10-28 2016-02-24 华中科技大学 一种镜像综合孔径辐射计反演方法
CN105629232B (zh) * 2015-12-24 2018-06-22 中国科学院电子学研究所 用于geosar相位定标的多地面定标站的布局方法
CN105527622B (zh) * 2015-12-24 2018-07-10 中国科学院电子学研究所 用于geosar相位定标的单地面定标站的布局方法
CN105954748B (zh) * 2016-04-26 2018-07-10 中国科学院电子学研究所 基于分块策略的TS-InSAR大气相位滤波方法
CN105974413B (zh) * 2016-06-13 2018-05-04 西安电子科技大学 多基地外辐射源雷达成像系统的自聚焦方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4084158A (en) * 1977-01-03 1978-04-11 Raytheon Company Method of operating synthetic aperture radar
JPS58120771A (ja) * 1982-01-09 1983-07-18 Kawasaki Steel Corp 合金化亜鉛めつき鋼板の製造方法

Also Published As

Publication number Publication date
CA1257370A (en) 1989-07-11
US4758838A (en) 1988-07-19
JPH0646405B2 (ja) 1994-06-15

Similar Documents

Publication Publication Date Title
JPS6165366A (ja) 合成開口レ−ダ画像再生処理方法
Cumming et al. Digital processing of SEASAT SAR data
US4985704A (en) Processing parameter generator for synthetic aperture radar
CN111142105B (zh) 复杂运动目标isar成像方法
CA1222809A (en) Method of reconstructing synthetic aperture radar image
Mason et al. Passive synthetic aperture radar imaging using low-rank matrix recovery methods
KR101839041B1 (ko) 연속 이동 효과 보정 기법을 이용하는 fmcw-sar 시스템 및 이 시스템의 sar 영상 복원 방법
CN107229050B (zh) 一种基于极坐标格式的雷达成像优化方法
CN111175727A (zh) 一种基于条件波数谱密度的宽带信号方位估计的方法
Saeedi et al. Improved navigation-based motion compensation for LFMCW synthetic aperture radar imaging
CN115453530A (zh) 一种基于参数化模型的双基sar滤波反投影两维自聚焦方法
JPS6280767A (ja) 合成開口レ−ダ画像の再生処理方式
JPS61187079A (ja) 合成開口レーダ画像の再生処理方法
Mason et al. Passive synthetic aperture radar imaging based on low-rank matrix recovery
CN112859018A (zh) 一种基于图像几何校正的视频sar成像方法
KR101839042B1 (ko) Fmcw-sar 시스템의 연속 이동 효과 보정 장치 및 방법
JPH0532712B2 (ja)
CN111722224A (zh) 基于Keystone变换的地基合成孔径雷达三维成像方法
JP2511939B2 (ja) 合成開口レ―ダ画像の再成処理方法
Zhang et al. A comparison of PCA based imaging methods for the multireceiver SAS
JPH0432350B2 (ja)
KR102159695B1 (ko) 인터리브드 sar 영상을 복원하는 장치 및 방법
JPH07120405B2 (ja) 合成開口レ−ダの画像再生処理方法
JPS61194570A (ja) 合成開口レーダ画像の焦点化方法
JPS63241484A (ja) デ−タ処理システム