JPS6020162A - 合成開口レ−ダ画像再生処理方式 - Google Patents
合成開口レ−ダ画像再生処理方式Info
- Publication number
- JPS6020162A JPS6020162A JP58127687A JP12768783A JPS6020162A JP S6020162 A JPS6020162 A JP S6020162A JP 58127687 A JP58127687 A JP 58127687A JP 12768783 A JP12768783 A JP 12768783A JP S6020162 A JPS6020162 A JP S6020162A
- Authority
- JP
- Japan
- Prior art keywords
- data
- image
- processing
- picture
- range
- 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
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/9004—SAR image acquisition techniques
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔発明の利用分野〕
本発明は人工衛星あるいは航空機等に搭載される合成開
口レーダ（５Ｙｎｔ１１ｅｔｉＣＡｐｅｒｔｕｒｅ３ａ
ｄａｒ 以下、ｌ’−８ＡＲＪという。）による撮像デ
ータから、人間が理解できる画像を再生するためのディ
ジタ々処理システムにおける高画質画（象再生処理方式
に関する。
口レーダ（５Ｙｎｔ１１ｅｔｉＣＡｐｅｒｔｕｒｅ３ａ
ｄａｒ 以下、ｌ’−８ＡＲＪという。）による撮像デ
ータから、人間が理解できる画像を再生するためのディ
ジタ々処理システムにおける高画質画（象再生処理方式
に関する。
人工衛星あるいは航空機等を用いたリモートセンシング
の分野では、地表を撮像するためのセンサとして、雲を
透過するマイクロ波帯で高解像の画像が得られるＳＡＲ
が注目されている。
の分野では、地表を撮像するためのセンサとして、雲を
透過するマイクロ波帯で高解像の画像が得られるＳＡＲ
が注目されている。
第１図にＳＡＲの全体システムを示した。レーダ・セッ
サ几８．アンテナＡｍを有するＳ Ａ ＩＬは、人工衛
星等に搭載されて飛行経路Ｆｐ上を矢印り方向に移動し
つつ地表の撮像を行う。ＳＡＲからの撮像データは地上
局Ｌ８で受信され、データ・プロセッサＤｐによ逆処理
されて映像フィルムＩＦ’の作成、データ記憶用磁気デ
ープＭＴの作成等が行われる。なお、Ｃは分乃了セルを
、几ａはＳ Ａ Ｒで採取されるデータの地表上のレン
ジ方向を、Ａｚは同アジマス方向を、Ａｎはアンテナゼ
ームをセしてＣｗは刈ｐ幅をそれぞれ示している。
サ几８．アンテナＡｍを有するＳ Ａ ＩＬは、人工衛
星等に搭載されて飛行経路Ｆｐ上を矢印り方向に移動し
つつ地表の撮像を行う。ＳＡＲからの撮像データは地上
局Ｌ８で受信され、データ・プロセッサＤｐによ逆処理
されて映像フィルムＩＦ’の作成、データ記憶用磁気デ
ープＭＴの作成等が行われる。なお、Ｃは分乃了セルを
、几ａはＳ Ａ Ｒで採取されるデータの地表上のレン
ジ方向を、Ａｚは同アジマス方向を、Ａｎはアンテナゼ
ームをセしてＣｗは刈ｐ幅をそれぞれ示している。
以下、ＳＡＲで採取されたデータの処理の概要を述べる
。なお、詳ｒｒ４１Ｉについては、例えば〔安藤”合成
開口レーダと間接Ｕ［測技術″計測と制御。
。なお、詳ｒｒ４１Ｉについては、例えば〔安藤”合成
開口レーダと間接Ｕ［測技術″計測と制御。
ｖｏｚ、ｚ ２ 、 Ａ ２〕を参照されたい。
ＳＡＲの受信画像中においては、原画像上の１点が点像
パターン”（Ｘ＋ｙ）の広が）をもって分布しており、
このま丑ではオＵ用できない。ここでＸは前記レンジ方
向を、ｙは前記アジマス方向を示している。前記受信画
像中で広がっている情報は、まずレンジ方向に圧縮され
、次にアジマス方向に圧縮される。この様子を第２図に
示す。第２図へは地表にマイクロ波反射点が２点だけ存
在したときの受信画像を模式的に示したものであるが、
２方向に圧縮処理を行えば、Ｂのようにもとの地表パタ
ーンをイツることかできる。前記レンジ圧縮処理は画ｌ
象デーク１ラインごとの点１象パターンデータとのコリ
レーション処理によって行う。
パターン”（Ｘ＋ｙ）の広が）をもって分布しており、
このま丑ではオＵ用できない。ここでＸは前記レンジ方
向を、ｙは前記アジマス方向を示している。前記受信画
像中で広がっている情報は、まずレンジ方向に圧縮され
、次にアジマス方向に圧縮される。この様子を第２図に
示す。第２図へは地表にマイクロ波反射点が２点だけ存
在したときの受信画像を模式的に示したものであるが、
２方向に圧縮処理を行えば、Ｂのようにもとの地表パタ
ーンをイツることかできる。前記レンジ圧縮処理は画ｌ
象デーク１ラインごとの点１象パターンデータとのコリ
レーション処理によって行う。
但シ、コリレーション処理全そのまま実行すると、美大
な処理時間がかかるため、高速フーリエ変換（以下、「
ＦＦＴ」という。）複素乗算、高速逆フーリエ菱換（以
下、「■ＦＦＴ」という。）を用いて高速化が図られる
。
な処理時間がかかるため、高速フーリエ変換（以下、「
ＦＦＴ」という。）複素乗算、高速逆フーリエ菱換（以
下、「■ＦＦＴ」という。）を用いて高速化が図られる
。
一方、アジマス圧縮処理に関しては、前記受（ｉ画像上
の情報のアジマス方向の広がシ（点像〕くターン）が、
点像のアジマス方向の位置に応じて変化することが考え
られるため、これまでのＳＡＲ。
の情報のアジマス方向の広がシ（点像〕くターン）が、
点像のアジマス方向の位置に応じて変化することが考え
られるため、これまでのＳＡＲ。
画像再生処理システムでは３つの立場がとられてきた。
第１は、処理画面内で点１象ノくターン（ま一定と仮定
し、前記レンジ圧縮処理と同等の処理を加えるものであ
シ、上記点１象変化が大きい場付には、画面内の一部で
しか焦点の合わない出力画像が得られるという欠点があ
る。第２は、処理画面を細分し、それぞれの分割画像で
点像）くターンを決めるものであシ、分割を細かくすれ
ば焦点は合うが、処理サイズに無関係に必要な固定的演
算Ｊｉｔが１杉響し、全体としてｍ＞Ｎ時間が大きくな
るという欠点がある。第３は、前記レンジ圧縮とは異な
り、点像関数を変化させつつ実空間で相関ｆｉｌ′ＪＴ
：を行うものであり、焦点精度は最良となるが演算時間
はＦＦＴを使用する場合に比べ数十倍になる。
し、前記レンジ圧縮処理と同等の処理を加えるものであ
シ、上記点１象変化が大きい場付には、画面内の一部で
しか焦点の合わない出力画像が得られるという欠点があ
る。第２は、処理画面を細分し、それぞれの分割画像で
点像）くターンを決めるものであシ、分割を細かくすれ
ば焦点は合うが、処理サイズに無関係に必要な固定的演
算Ｊｉｔが１杉響し、全体としてｍ＞Ｎ時間が大きくな
るという欠点がある。第３は、前記レンジ圧縮とは異な
り、点像関数を変化させつつ実空間で相関ｆｉｌ′ＪＴ
：を行うものであり、焦点精度は最良となるが演算時間
はＦＦＴを使用する場合に比べ数十倍になる。
以上、アジマス圧縮処理に関し既存のＳＡＲＡ像再生処
理システムでは、点像のアジマス方向変化を無視しＦＦ
Ｔを使用して高速化を図るか、Ｆ Ｉ’ Ｔを使用せず
に高焦点・高画ＪＪ！ｉ、を倚るかどちらかであった。
理システムでは、点像のアジマス方向変化を無視しＦＦ
Ｔを使用して高速化を図るか、Ｆ Ｉ’ Ｔを使用せず
に高焦点・高画ＪＪ！ｉ、を倚るかどちらかであった。
これらに対し、画像修復の分野では、原画像を座標変換
し、変換画像上では点像（ぶれパターン等）が不変とな
るようにし、実際には点像が画像上の位置によシ変化し
ても、ＦＦＴを使用して高速にしかも高精度に画像修復
を行う方式が提案されている。なお詳細については〔Ａ
。
し、変換画像上では点像（ぶれパターン等）が不変とな
るようにし、実際には点像が画像上の位置によシ変化し
ても、ＦＦＴを使用して高速にしかも高精度に画像修復
を行う方式が提案されている。なお詳細については〔Ａ
。
Ａ、Ｓａｗｃｈｕｋ ： ５ｐａｃｅ−ｖａｒｉａｎｔ
ｉｍａｇｅｒｅｓｔｏｒａｔｉｏｎ ｂｙ ｃｏｏｒ
ｄｉｎａｔｅ ｔｒａｎｓｆｏｒｍａｔｉｏｎ。
ｉｍａｇｅｒｅｓｔｏｒａｔｉｏｎ ｂｙ ｃｏｏｒ
ｄｉｎａｔｅ ｔｒａｎｓｆｏｒｍａｔｉｏｎ。
Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ Ｑｐｔｉｃａｌ ５ｏ
ｃｉｃｔｙ ｏｆＡＪＴｌｅｒｉｃａ 、 Ｆｅｂ、
、 １９７４．：］を参照されたい。しかし、ＳＡＲに
おける下式の点像関数に対しては、変換画像上で点像が
不変となる座標変換は存在しないことが知られており、
Ｆ’Ｆ”Ｔを用いる上記方式は適用できなかった。
ｃｉｃｔｙ ｏｆＡＪＴｌｅｒｉｃａ 、 Ｆｅｂ、
、 １９７４．：］を参照されたい。しかし、ＳＡＲに
おける下式の点像関数に対しては、変換画像上で点像が
不変となる座標変換は存在しないことが知られており、
Ｆ’Ｆ”Ｔを用いる上記方式は適用できなかった。
本発明の目的は、従来の可変点像に対する重速（ｌｉ”
Ｉ” Ｔ ［用）尚焦点ｌＪＭＪ像再生方式における
上述の如き欠点を除去し、ＳＡ几両画像処理システムお
いて、Ｉｉ”ＦＴを使用しかつ画面内いたる所で焦点の
合った高画質のＳＡ几両画像出力可能とするＳＡ几両画
像再生処理方式を提供することにある。
Ｉ” Ｔ ［用）尚焦点ｌＪＭＪ像再生方式における
上述の如き欠点を除去し、ＳＡ几両画像処理システムお
いて、Ｉｉ”ＦＴを使用しかつ画面内いたる所で焦点の
合った高画質のＳＡ几両画像出力可能とするＳＡ几両画
像再生処理方式を提供することにある。
上記目的を達成するため本発明では、衛星あるいは航空
機の軌道データから再生画面内でのドツプラ変化率（式
（１）の係数■ぐ）時間変化データをめ、ドツプラ斐化
率時間変化データから、レンジ圧縮転置画像のアジマス
方向厘ｍＳ＆換デー、！２をめ、アジマス方向座標変換
データによりレンジ圧縮転置画像のアジマス方向内挿演
ｉｔ−処理を行い、座標変換後の画像に関する点１故の
ドツプラ変化率が画面内のアジマス方向位置にかかわら
ず不変としたうえでＦ Ｆ Ｔを使用する点に特徴があ
る。
機の軌道データから再生画面内でのドツプラ変化率（式
（１）の係数■ぐ）時間変化データをめ、ドツプラ斐化
率時間変化データから、レンジ圧縮転置画像のアジマス
方向厘ｍＳ＆換デー、！２をめ、アジマス方向座標変換
データによりレンジ圧縮転置画像のアジマス方向内挿演
ｉｔ−処理を行い、座標変換後の画像に関する点１故の
ドツプラ変化率が画面内のアジマス方向位置にかかわら
ず不変としたうえでＦ Ｆ Ｔを使用する点に特徴があ
る。
実施例の説明に先立ち、本発明の詳細な説明する。第３
図は本光りＪの原理全説明するだめの図である。衛星が
合成囲ロレーダを搭載し図のようにル”６行すると、地
表点ターゲラ）３１と３２では衛星軌道との垂直距Ｐ４
１Ｆ、’１ ＋ ｒ２ 、速度Ｖ１＋ ｖ２の違いによ
り、ドツプラ変化率には、第３図（ｂ）に示すよりに、
下記Δにだけ変化する。
図は本光りＪの原理全説明するだめの図である。衛星が
合成囲ロレーダを搭載し図のようにル”６行すると、地
表点ターゲラ）３１と３２では衛星軌道との垂直距Ｐ４
１Ｆ、’１ ＋ ｒ２ 、速度Ｖ１＋ ｖ２の違いによ
り、ドツプラ変化率には、第３図（ｂ）に示すよりに、
下記Δにだけ変化する。
ここで λ：送信マイクロ波の波長
従って、画面内ｌライン上の全データに対し固冗のドツ
プラ変化率Ｋｏを仮定し、点像（１）式によりアジマス
圧禰処理を行うと再生画面内では、兵のＫとＫｏの走に
応じたふれ劣化が生ずる。
プラ変化率Ｋｏを仮定し、点像（１）式によりアジマス
圧禰処理を行うと再生画面内では、兵のＫとＫｏの走に
応じたふれ劣化が生ずる。
今、画像雁漂系を局所的に伸縮（腺型座標荻換）させる
ことを考える。原画像座標系をｙ１変換画１ぶ座標糸全
ｙ′とし、’／＝ｋＹ’が成立つとする。
ことを考える。原画像座標系をｙ１変換画１ぶ座標糸全
ｙ′とし、’／＝ｋＹ’が成立つとする。
原画像座標系に関する点像の位相ψ（ｙ）がで表わされ
るとすると、変便後の座標系での点像の位相はψ（ｙ′
）は Ｊ（ ψ（ｙ′）＝ψ（ｋｙ’）＝−に２ｙ′２十ｋｂｙ’
・・・・・・（４）で光わされる。したがって、Ｙ ＝
Ｉ（ｙ ’ なる局所的線型座標変換では、ドツプラ
変化率Ｋ ＩｄＪ＜ ” Ｋと変化する。
るとすると、変便後の座標系での点像の位相はψ（ｙ′
）は Ｊ（ ψ（ｙ′）＝ψ（ｋｙ’）＝−に２ｙ′２十ｋｂｙ’
・・・・・・（４）で光わされる。したがって、Ｙ ＝
Ｉ（ｙ ’ なる局所的線型座標変換では、ドツプラ
変化率Ｋ ＩｄＪ＜ ” Ｋと変化する。
以上より、座標ｙごとに変化する真のトップ２変化率Ｋ
（３’）を座標変換ｙ’＝ｙ’（ｙ）によシ固定のド
ツプラ変化率Ｋｏとするためには、下記倣分方程式を前
記座標変換ｙ’＝ｙ’（ｙ）がイｖｔせばよい。
（３’）を座標変換ｙ’＝ｙ’（ｙ）によシ固定のド
ツプラ変化率Ｋｏとするためには、下記倣分方程式を前
記座標変換ｙ’＝ｙ’（ｙ）がイｖｔせばよい。
ドツプラ変化率の画面内変化Ｋ （３’）力補Ａ型であ
る場合すなわちＫ （ｙ） ＝ １Ｃ＋ΔＫ −Ｙの場
合には、とすればよい。本座標変換式の形゛Ｕζを＃＋
Ｔ ３ ｒ−ａ （ｃ）に示す。
る場合すなわちＫ （ｙ） ＝ １Ｃ＋ΔＫ −Ｙの場
合には、とすればよい。本座標変換式の形゛Ｕζを＃＋
Ｔ ３ ｒ−ａ （ｃ）に示す。
以上、本発明の原理を述べたが、以下、第４図によシ本
発明の詳細な説明する。合成開口レーダ原データ５は、
公知のレンジ圧縮処理装置６によシレンジ圧縮され、画
像用ディスク装置１０を用いて転置された佼、座標変換
すなわち画像内挿処理（リサンプル）を後述する画像内
挿処理装置１１において受ける。レンジ圧縮装置６では
、データ５はＦＦＴ装置７によシフーリエ変換され、複
素乗算装置８によシ固定的レンジ参照関数とかけ合わせ
られる。さらにＰＦＴ装置９により逆７−リエ変換され
、レンジ圧縮後のデータとなる。
発明の詳細な説明する。合成開口レーダ原データ５は、
公知のレンジ圧縮処理装置６によシレンジ圧縮され、画
像用ディスク装置１０を用いて転置された佼、座標変換
すなわち画像内挿処理（リサンプル）を後述する画像内
挿処理装置１１において受ける。レンジ圧縮装置６では
、データ５はＦＦＴ装置７によシフーリエ変換され、複
素乗算装置８によシ固定的レンジ参照関数とかけ合わせ
られる。さらにＰＦＴ装置９により逆７−リエ変換され
、レンジ圧縮後のデータとなる。
ＦＦＴ装置９の一実施例を第６図で示す。これは、２″
個の１ラインのデータをパイプライン処理によ９０段の
バタフライ演算器５１ ＋ ５２・・・・・・５ｎを用
いてＦ Ｆ Ｔ処理するものである。上記各＠算器間に
は、データバッファ６１２ ｒ ６２１．・・・・・・
６（ｎ−Ｉｎｎが設けられている。２段目の演算器５２
が第１回目のバタフライ■算で必要とする入力データは
、１段目の演算器５１の第１回目および第２回目の演ｊ
ｉ−＞二結果の両方にまたがっており、従ってデータバ
ッファ６１２は４デ一ク分の容量が必要である。
個の１ラインのデータをパイプライン処理によ９０段の
バタフライ演算器５１ ＋ ５２・・・・・・５ｎを用
いてＦ Ｆ Ｔ処理するものである。上記各＠算器間に
は、データバッファ６１２ ｒ ６２１．・・・・・・
６（ｎ−Ｉｎｎが設けられている。２段目の演算器５２
が第１回目のバタフライ■算で必要とする入力データは
、１段目の演算器５１の第１回目および第２回目の演ｊ
ｉ−＞二結果の両方にまたがっており、従ってデータバ
ッファ６１２は４デ一ク分の容量が必要である。
しかし、上記４デ一タ分のバタフライ演算を終了すれば
、それ以降このデータを入力することはなく、演算器５
１の次のバタフライ演算結果で前記データバッフ７６１
２の内容を破壊しても問題はない。結局、データバッフ
ァ６１２の容量は４デ一タ分で良いことになる。同様に
、データバッファ６ユ、の容量は８デ一タ分であり、最
終段バッファ６（ｎ−１）ｎの料量は２ｎデ一タ分とな
る。
、それ以降このデータを入力することはなく、演算器５
１の次のバタフライ演算結果で前記データバッフ７６１
２の内容を破壊しても問題はない。結局、データバッフ
ァ６１２の容量は４デ一タ分で良いことになる。同様に
、データバッファ６ユ、の容量は８デ一タ分であり、最
終段バッファ６（ｎ−１）ｎの料量は２ｎデ一タ分とな
る。
上記各データバッファは交替バッファであシ、書込み、
読出しを同時に行うことができるものであシ、これによ
シ全演算器が同時に休みなく稼動できる。交替バッファ
の書込みアドレス、読出しアドレスは制御装置７によシ
同期をとって制御される。
読出しを同時に行うことができるものであシ、これによ
シ全演算器が同時に休みなく稼動できる。交替バッファ
の書込みアドレス、読出しアドレスは制御装置７によシ
同期をとって制御される。
上記各演算器による演算内容は、２個の複素データＡ、
Ｂと複素係数ωとを第７図（イ）に示す如く演算するも
のである。第７図（ロ）は上ｉｆ！２演算内存をより具
体的に示したものでおり、２個のｔＮ、ＡデータＡ、Ｂ
の実数都をＡ、、Ｂ、 、虚数部をＡ２１Ｂ２とし、複
素係数ωの実故部をω１、虚数部をω２とするとき、図
のａ −ｊまての処理、すなわち、４回の実数乗算、６
回の火数加ｎを行うことによシ２制の複素データＣ，Ｄ
（その実数都をＣＩ＋ Ｄｔ ｓ 虚Ｍ部を０２ ＋
Ｄｚ とする。）を出力するものである。
Ｂと複素係数ωとを第７図（イ）に示す如く演算するも
のである。第７図（ロ）は上ｉｆ！２演算内存をより具
体的に示したものでおり、２個のｔＮ、ＡデータＡ、Ｂ
の実数都をＡ、、Ｂ、 、虚数部をＡ２１Ｂ２とし、複
素係数ωの実故部をω１、虚数部をω２とするとき、図
のａ −ｊまての処理、すなわち、４回の実数乗算、６
回の火数加ｎを行うことによシ２制の複素データＣ，Ｄ
（その実数都をＣＩ＋ Ｄｔ ｓ 虚Ｍ部を０２ ＋
Ｄｚ とする。）を出力するものである。
演Ｊ７器５１〜５ｎの実施例を第８図に示す。順次入力
し／こ４つの実数データは、係数テーブル８２の係数と
、乗５’１．’ａｓ１により次々に乗算され、相互に加
ｎ器８５にょシ加ｎ（減算を含む）される。乗算器８１
と加算器８５との間には、データレジスタ８４が置かれ
、加算結果の中間バッファとして用いられる。上述の各
装置はマイクロプログラムによる制御装置８３から、動
作およびアドレスの制御を受ける。
し／こ４つの実数データは、係数テーブル８２の係数と
、乗５’１．’ａｓ１により次々に乗算され、相互に加
ｎ器８５にょシ加ｎ（減算を含む）される。乗算器８１
と加算器８５との間には、データレジスタ８４が置かれ
、加算結果の中間バッファとして用いられる。上述の各
装置はマイクロプログラムによる制御装置８３から、動
作およびアドレスの制御を受ける。
煩累乗算装置８は、第８図の（ｆＴＬｎ器において、係
数テーブル８２を取除き、入カデータバッンア８６を付
加しｆｃ第９図の演算器で実現できる。複素データ間の
積は乗算とそれに続く加算にょシ実行できるからである
。これに先立ち、演算装置２は衛星軌道・地球モデルデ
ー月がらドツプラ変化率の画面内変化Ｋ （Ｙ）を式（
匂にょ請求めその出力は、式（５）を屏く機能を有す微
分方程式演算装置３に入力され、座標変換ｙ ’ （ｙ
）ずなわち内挿位置データが装置１１に送られる。座標
敦換後のレンジ圧縮済データは、公知のアジマス圧縮装
置ｆｉ：１２によ）アジマス圧縮を受けるが、この際使
用される点像関数のトップ２変化率は式（５）の定数Ｊ
（Ｏである。なお、画像データのアジマス方向座標変換
のための内挿処理装置１１は、第５図（１））に示すよ
うに、シフトレジスタ２１．乗算器２３、加ｐ−器２４
、内挿荷重ルックアップテーブル２２により実現できる
。第５図（ａ）により、ディジタル画像データの座標変
換のだめのデータ内挿処理を説明する。装置１０からの
画像データ１ラインの値を（・・・ｇａ＋ｇｂ＋・・・
）とすると、データ内挿は、座標変換後のデータ位置（
・・・１，２．・・・）での値（・・・ｇｌ ＋ ｇ２
＋・・・）を４＝Ｊ近のデータからめることであり、
（ｂ）の実施例ではこれを４点の隣接データの荷重和で
めるものとしである。ただし、荷重は内挿する位置によ
シ決まるため、あらかじめめておき、座標変換式により
まる内挿位置（・・・１，２．・・・）による装置２２
のテーブルルックアップ方式とするのが良い。
数テーブル８２を取除き、入カデータバッンア８６を付
加しｆｃ第９図の演算器で実現できる。複素データ間の
積は乗算とそれに続く加算にょシ実行できるからである
。これに先立ち、演算装置２は衛星軌道・地球モデルデ
ー月がらドツプラ変化率の画面内変化Ｋ （Ｙ）を式（
匂にょ請求めその出力は、式（５）を屏く機能を有す微
分方程式演算装置３に入力され、座標変換ｙ ’ （ｙ
）ずなわち内挿位置データが装置１１に送られる。座標
敦換後のレンジ圧縮済データは、公知のアジマス圧縮装
置ｆｉ：１２によ）アジマス圧縮を受けるが、この際使
用される点像関数のトップ２変化率は式（５）の定数Ｊ
（Ｏである。なお、画像データのアジマス方向座標変換
のための内挿処理装置１１は、第５図（１））に示すよ
うに、シフトレジスタ２１．乗算器２３、加ｐ−器２４
、内挿荷重ルックアップテーブル２２により実現できる
。第５図（ａ）により、ディジタル画像データの座標変
換のだめのデータ内挿処理を説明する。装置１０からの
画像データ１ラインの値を（・・・ｇａ＋ｇｂ＋・・・
）とすると、データ内挿は、座標変換後のデータ位置（
・・・１，２．・・・）での値（・・・ｇｌ ＋ ｇ２
＋・・・）を４＝Ｊ近のデータからめることであり、
（ｂ）の実施例ではこれを４点の隣接データの荷重和で
めるものとしである。ただし、荷重は内挿する位置によ
シ決まるため、あらかじめめておき、座標変換式により
まる内挿位置（・・・１，２．・・・）による装置２２
のテーブルルックアップ方式とするのが良い。
ここで公知のアジマス圧縮装置１２は、ＦＦＴ装置１３
、レンジカーバチャ補正装置１４、煩累乗算装置１５、
Ｉｉ”ｌｉ’Ｔ装置１６、画像化装置１７からなる。装
置１４は、レンジ圧縮後の点像が、衛星と地表点との距
離の変化に起因し、画像上でわん曲しているものを補正
するためのもので、レンジ方向の座標変換内挿処理を行
う。装置１７は、複素データを整数データに変換するも
のである。
、レンジカーバチャ補正装置１４、煩累乗算装置１５、
Ｉｉ”ｌｉ’Ｔ装置１６、画像化装置１７からなる。装
置１４は、レンジ圧縮後の点像が、衛星と地表点との距
離の変化に起因し、画像上でわん曲しているものを補正
するためのもので、レンジ方向の座標変換内挿処理を行
う。装置１７は、複素データを整数データに変換するも
のである。
なお、装置１２の再生出力は、原画像と同等の座標系に
変換される。
変換される。
なお、前述のレンジカーバテヤ補正装置１４は、上記の
第５図に示したデータ内挿処理装置と同一のものであシ
、画像化装置１７は複素数のノルム計算であシ、第９図
の核素乗ｎ装置と同一のもので実現される。
第５図に示したデータ内挿処理装置と同一のものであシ
、画像化装置１７は複素数のノルム計算であシ、第９図
の核素乗ｎ装置と同一のもので実現される。
座標変換装置１９で逆換され、出力記憶媒体２０に出力
される。装置１９は、装置１１と同等のものであるが、
逆変換用データｙ（ｙ勺を逆変換算出装置４から入力す
る点が異なる。
される。装置１９は、装置１１と同等のものであるが、
逆変換用データｙ（ｙ勺を逆変換算出装置４から入力す
る点が異なる。
本発明によれば、曾成囲ロレーダの画啄再生処理におい
て、衛星や航を偵の速度、高度の時間変化による１画面
内の点像パタン変化データをめ、原画像の座標変換処理
によｐ変換後の画像上では点像パタンか不変となるよう
にし、その結果高速７−リエ変換（ＦＦＴ）を使用し処
理高速性を達成するとともに画面内いたる所で焦点の合
った高画質のＳＡＲ画像再生処理方式を提供できるとい
う効果がある。
て、衛星や航を偵の速度、高度の時間変化による１画面
内の点像パタン変化データをめ、原画像の座標変換処理
によｐ変換後の画像上では点像パタンか不変となるよう
にし、その結果高速７−リエ変換（ＦＦＴ）を使用し処
理高速性を達成するとともに画面内いたる所で焦点の合
った高画質のＳＡＲ画像再生処理方式を提供できるとい
う効果がある。
第１図は衛星搭載合成υト］ロレーダおよび地上局シス
テム全体の概要を示す図、第２図は愛情画像上で広がっ
た地表の２点の情報を再生画像処理によ９２点に圧縮し
た様子を示す図、第３図は本発明の原理を示す図、第４
図は本光切による衛星地上局の実施例を示す図、第５図
は画像座標変換に必要なデータ内挿処理装置ｆｔのブロ
ック構成を示す図、第６図はＦＦＴ装置の一実施例を示
す図、第７図は複素演算器の４ｔ’ｌ成を示す図、第８
図は俵素乗算器の原形となる演昇器の構成を示す図、第
９図は複素乗算器の構成を示す図である。 罰 １ 図 某 ２ 図 Ａ Ｉ５ 第 ３ 図 Ｃ０−） ’ｌ／’− 片 乙 区 某 ７ 図 （Ａ） 菫 ３ 目 第 ９ （２）
テム全体の概要を示す図、第２図は愛情画像上で広がっ
た地表の２点の情報を再生画像処理によ９２点に圧縮し
た様子を示す図、第３図は本発明の原理を示す図、第４
図は本光切による衛星地上局の実施例を示す図、第５図
は画像座標変換に必要なデータ内挿処理装置ｆｔのブロ
ック構成を示す図、第６図はＦＦＴ装置の一実施例を示
す図、第７図は複素演算器の４ｔ’ｌ成を示す図、第８
図は俵素乗算器の原形となる演昇器の構成を示す図、第
９図は複素乗算器の構成を示す図である。 罰 １ 図 某 ２ 図 Ａ Ｉ５ 第 ３ 図 Ｃ０−） ’ｌ／’− 片 乙 区 某 ７ 図 （Ａ） 菫 ３ 目 第 ９ （２）
Claims (1)
- 合成開口レーダを搭載した飛翔体の軌道および地表のモ
デルから再生画像面内の任煮の位置におけるドツプラ変
化率をめるステップと、該ドツプラ変化率の時間変化を
示すデータからレンジ圧縮された再生画像のアジマス方
向座標変換データをめるステップと、該座標変換データ
にもとづき再生画像データにたいしてアジマス方向に内
挿演算処理をおこなうステップと、内挿演算処理のおこ
なわれた再生画像にたいしてＦＦＴによるアジマス方向
に圧縮演算処理をおこなうステップとからなることを特
徴とする合成−開口レーダ画像再生処理方式。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58127687A JPS6020162A (ja) | 1983-07-15 | 1983-07-15 | 合成開口レ−ダ画像再生処理方式 |
| EP84108152A EP0131907A1 (en) | 1983-07-15 | 1984-07-11 | Method of reconstructing synthetic aperture radar image |
| CA000458614A CA1222809A (en) | 1983-07-15 | 1984-07-11 | Method of reconstructing synthetic aperture radar image |
| US06/630,200 US4616227A (en) | 1983-07-15 | 1984-07-12 | Method of reconstructing synthetic aperture radar image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58127687A JPS6020162A (ja) | 1983-07-15 | 1983-07-15 | 合成開口レ−ダ画像再生処理方式 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6020162A true JPS6020162A (ja) | 1985-02-01 |
| JPH0456268B2 JPH0456268B2 (ja) | 1992-09-07 |
Family
ID=14966213
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58127687A Granted JPS6020162A (ja) | 1983-07-15 | 1983-07-15 | 合成開口レ−ダ画像再生処理方式 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4616227A (ja) |
| EP (1) | EP0131907A1 (ja) |
| JP (1) | JPS6020162A (ja) |
| CA (1) | CA1222809A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7551119B1 (en) * | 2008-01-08 | 2009-06-23 | Sandia Corporation | Flight path-driven mitigation of wavefront curvature effects in SAR images |
Families Citing this family (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4855747A (en) * | 1987-08-17 | 1989-08-08 | Trustees Of The University Of Pennsylvania | Method of target imaging and identification |
| US4801939A (en) * | 1987-09-04 | 1989-01-31 | Unisys Corporation | High-speed data compressor/decompressor for synthetic aperture radar |
| 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 |
| US5043734A (en) * | 1988-12-22 | 1991-08-27 | Hughes Aircraft Company | Discrete autofocus for ultra-high resolution synthetic aperture radar |
| JPH0727021B2 (ja) * | 1989-02-10 | 1995-03-29 | 三菱電機株式会社 | 合成開口レーダ装置 |
| DE4002857C1 (en) * | 1990-02-01 | 1991-09-19 | Dornier Gmbh, 7990 Friedrichshafen, De | Ascertaining mean Doppler frequency of SAR system in real=time - reading complex data in memory module as range lines and reading out as azimuth LINES |
| SE468184B (sv) * | 1991-03-25 | 1992-11-16 | Hans Hellsten | Syntetisk aperturradar |
| DE4122592C1 (ja) * | 1991-07-08 | 1993-01-28 | Deutsche Forschungsanstalt Fuer Luft- Und Raumfahrt Ev, 5300 Bonn, De | |
| US5185609A (en) * | 1991-10-29 | 1993-02-09 | Wilcox Electric, Inc. | Signal monitor utilizing digital signal processing |
| US5281972A (en) * | 1992-09-24 | 1994-01-25 | Hughes Aircraft Company | Beam summing apparatus for RCS measurements of large targets |
| US5343204A (en) * | 1993-07-29 | 1994-08-30 | Unisys Corporation | Auto-focusing correction for rotational acceleration effects on inverse synthetic aperture radar images |
| US5570099A (en) * | 1993-10-15 | 1996-10-29 | Loral Federal Systems Company | TDOA/FDOA technique for locating a transmitter |
| US5469167A (en) * | 1993-10-18 | 1995-11-21 | The United States Of America As Represented By The Secretary Of The Army | Synthetic aperture radar for nonlinear trajectories using range relative doppler processing and invariant mapping |
| SE503650C2 (sv) * | 1994-11-10 | 1996-07-22 | Hans Hellsten | Sätt att vid en radar förbättra radardata |
| JP2757845B2 (ja) * | 1995-12-26 | 1998-05-25 | 日本電気株式会社 | 合成開口レーダ装置 |
| US5745069A (en) * | 1996-09-10 | 1998-04-28 | Ball Corporation | Reduction of radar antenna area |
| DE10346047A1 (de) * | 2003-10-02 | 2005-04-28 | Eads Deutschland Gmbh | Verfahren und Vorrichtung nach dem FMCW-Radarprinzip |
| US7307580B2 (en) | 2006-01-17 | 2007-12-11 | Raytheon Company | Non-statistical method for compressing and decompressing complex SAR data |
| US7920262B2 (en) * | 2008-09-17 | 2011-04-05 | The United States Of America As Represented By The Secretary Of The Army | Systems for measuring backscattered light using rotating mirror |
| US7928901B2 (en) * | 2008-10-16 | 2011-04-19 | The United States Of America As Represented By The Secretary Of The Army | Systems and methods for producing radar images |
| US8362947B2 (en) * | 2011-03-29 | 2013-01-29 | Victor Gorelik | Method for obtaining object-plane field from its two images |
| CN112986996B (zh) * | 2021-02-08 | 2023-08-25 | 中国科学院空天信息创新研究院 | 基于几何熵的多源sar卫星联合三维定位方法 |
| CN112946648B (zh) * | 2021-03-22 | 2023-07-18 | 中国电子科技集团公司第三十八研究所 | 一种大前斜视sar空变运动误差补偿方法 |
| CN113627029B (zh) * | 2021-08-24 | 2022-07-12 | 南京大学 | 一种地月三角平动点附近的轨道设计方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS57179769A (en) * | 1981-04-30 | 1982-11-05 | Toshiba Corp | Correcting system for picture distortion of synthetic aperture-face radar |
| JPS58113877A (ja) * | 1981-12-28 | 1983-07-06 | Nec Corp | 合成開口レ−ダオ−トフオ−カス画像処理装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5822982A (ja) * | 1981-08-03 | 1983-02-10 | Hitachi Ltd | 合成開口レ−ダの画像処理システム |
| DE3279581D1 (en) * | 1981-12-28 | 1989-05-03 | Nec Corp | Synthetic aperture radar image processing system |
-
1983
- 1983-07-15 JP JP58127687A patent/JPS6020162A/ja active Granted
-
1984
- 1984-07-11 EP EP84108152A patent/EP0131907A1/en not_active Withdrawn
- 1984-07-11 CA CA000458614A patent/CA1222809A/en not_active Expired
- 1984-07-12 US US06/630,200 patent/US4616227A/en not_active Expired - Lifetime
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS57179769A (en) * | 1981-04-30 | 1982-11-05 | Toshiba Corp | Correcting system for picture distortion of synthetic aperture-face radar |
| JPS58113877A (ja) * | 1981-12-28 | 1983-07-06 | Nec Corp | 合成開口レ−ダオ−トフオ−カス画像処理装置 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7551119B1 (en) * | 2008-01-08 | 2009-06-23 | Sandia Corporation | Flight path-driven mitigation of wavefront curvature effects in SAR images |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0456268B2 (ja) | 1992-09-07 |
| EP0131907A1 (en) | 1985-01-23 |
| US4616227A (en) | 1986-10-07 |
| CA1222809A (en) | 1987-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPS6020162A (ja) | 合成開口レ−ダ画像再生処理方式 | |
| US10976446B2 (en) | Acceleration of real time computer vision processing on UAVs through GPS attitude estimation | |
| JPS6165366A (ja) | 合成開口レ−ダ画像再生処理方法 | |
| IL97559A (en) | ESAHP ECNEREFFID autofocus for RAS Rader mapping | |
| CN106054188A (zh) | 无人机合成孔径雷达成像的图像偏移自聚焦方法 | |
| WO2024001360A1 (zh) | 绿幕抠图方法、装置及电子设备 | |
| Wang et al. | Application of the spherical harmonic gravity model in high precision inertial navigation systems | |
| Huang et al. | An atomic gravimeter dynamic measurement method based on Kalman filter | |
| SE468184B (sv) | Syntetisk aperturradar | |
| Massonnet et al. | CNES general-purpose SAR correlator | |
| Yang et al. | An LBL positioning algorithm based on an EMD–ML hybrid method | |
| Atemkeng et al. | Baseline-dependent sampling and windowing for radio interferometry: data compression, field-of-interest shaping, and outer field suppression | |
| KR101009967B1 (ko) | 레이더 영상 생성 장치와 그 방법, 및 상기 방법을 구현하는 프로그램이 기록된 기록매체 | |
| Mosavi et al. | Positioning performance analysis using RWLS algorithm based on variance estimation methods | |
| Eineder et al. | A flexible system for the generation of interferometric SAR products | |
| US4283763A (en) | System for graphic representation of vessel's position on Mercator map | |
| JP2016161386A (ja) | 画像レーダ装置 | |
| JPH0532712B2 (ja) | ||
| Lindstrom et al. | Sensitivity of bpa sar image formation to initial position, velocity, and attitude navigation errors | |
| Wu et al. | Sparse-data based 3d surface reconstruction for cartoon and map | |
| Krishnamurthy et al. | Auralization of Unsteady Rotor Noise using a Solution to the Ffowcs Williams-Hawkings Equation | |
| US2888203A (en) | Device for predicting values of a fluctuating system at a predetermined future time | |
| JPH0432350B2 (ja) | ||
| Clark | Geocoding and stereoscopy of synthetic aperture radar images | |
| Golubtsov et al. | Automatic detection and suppression of parasitic chirp signals in the intermediate time–frequency domain |