JP6066934B2 - 人工衛星またはエアリアルプラットフォーム上で複数のsar画像を生成するシステム、当該システムを備える衛星、人工衛星またはエアリアルプラットフォームで合成開口レーダ(sar)画像を生成する方法 - Google Patents

人工衛星またはエアリアルプラットフォーム上で複数のsar画像を生成するシステム、当該システムを備える衛星、人工衛星またはエアリアルプラットフォームで合成開口レーダ(sar)画像を生成する方法 Download PDF

Info

Publication number
JP6066934B2
JP6066934B2 JP2013557122A JP2013557122A JP6066934B2 JP 6066934 B2 JP6066934 B2 JP 6066934B2 JP 2013557122 A JP2013557122 A JP 2013557122A JP 2013557122 A JP2013557122 A JP 2013557122A JP 6066934 B2 JP6066934 B2 JP 6066934B2
Authority
JP
Japan
Prior art keywords
range
azimuth
frequency
compression
image
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.)
Expired - Fee Related
Application number
JP2013557122A
Other languages
English (en)
Other versions
JP2014510913A (ja
Inventor
ウォーカー ウィシャート アレクサンダー
ウォーカー ウィシャート アレクサンダー
Original Assignee
エアバス ディフェンス アンド スペイス リミテッド
エアバス ディフェンス アンド スペイス リミテッド
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 エアバス ディフェンス アンド スペイス リミテッド, エアバス ディフェンス アンド スペイス リミテッド filed Critical エアバス ディフェンス アンド スペイス リミテッド
Publication of JP2014510913A publication Critical patent/JP2014510913A/ja
Application granted granted Critical
Publication of JP6066934B2 publication Critical patent/JP6066934B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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
    • 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/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric 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)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Description

本願発明は、合成開口レーダ(SAR)データの処理に関する。
目標物の高解像度画像を生成するために合成開口レーダを用いることが出来る。放射パルスによって目標物が照射され、パルス反射が検出され処理され、画像データが形成される。衛星などの宇宙探査機にSARが用いられ得る。デジタル形式の未加工のエコーデータは、典型的にはレベル0のプロダクトと呼ばれる。レベル0のプロダクトから形成されたSAR画像は、典型的にはレベル1のプロダクトと呼ばれる。
典型的には、衛星でSARにより得られるレベル0の未加工データは地上でレベル1のプロダクトへと処理され、よって衛星は、レベル0のデータを、地上への送信が可能となるまで格納しておかなければならない。このことは、衛星が地上局と通信を行う機会が頻繁にない場合において課題となる。例えば、惑星ミッションにおいて、衛星はより長い時間、受信したデータを格納しておかなければならなくなる。さらに、地球低軌道(LEO)衛星などの衛星は、地上局と通信を行っていない間、データを格納しておかなければならない。よって、得られる画像の数は、衛星の記憶装置の容量によって制限されてしまう。
例えば干渉画像または複素コヒーレンス関数を生成すべく、衛星への2以上のSARアンテナの搭載が所望されることがある。SARアンテナの数を増やすことにより、生成されるレベル0の未加工データの量も増え、課題はさらに深刻なものとなる。
衛生で画像処理を実行することの利点が提案されている。例えば、低解像度画像が許容される場合、レベル0のデータを衛星で低解像度画像へと処理し、これにより、レベル0のデータを格納する場合と比較してより多くの画像のデータを衛星が格納出来るようにすることが出来る。さらに、いくつかの適用例においては、目標領域のリアルタイム、またはほぼリアルタイムのモニタリングを行うことが所望され、地上で処理を行うべく画像データを送信するにはあまりにも長い時間がかかり得る。
従来技術のSAR画像処理は地上のソフトウェアによって行われている。宇宙で船上にて行われる処理に適した現在利用可能なCPUは、地上で用いられる計算処理量が多いソフトウェアを効果的に実行するには処理速度が遅すぎる。これは、宇宙で船上にて行われる処理に適したフォールトトレラントな汎用CPUの処理速度は、地上のシステムと比較し数桁遅いからである。
本願発明は上記した課題およびその他の課題を解決する。
本願発明によると、人工衛星またはエアリアルプラットフォームに用いる装置は、SARデータのレンジラインまたは方位角ラインに対応し複数の目標点の情報を含む信号を、複数の周波数チャネルへ周波数多重分離する周波数多重分離器と、複数の周波数チャネルのそれぞれに対し画像形成を実行する圧縮手段とを含み、複数の周波数チャネルのそれぞれの信号は複数の目標点の情報を含む。
データが複数のチャネルへ周波数多重分離されるので、圧縮が実行されるチャネルの帯域幅はより狭くなり、圧縮をハードウェアで行うことが出来る。周波数多重分離器および圧縮手段の両方をハードウェアで実装できる。
このため、または他の理由により、本願発明は、人工衛星またはエアリアルプラットフォームでレンジ圧縮および方位角圧縮を行えるようにする。例えば、本願発明は、衛星または他の宇宙探査機でレンジ圧縮または方位角圧縮を行えるようにする。また本願発明は、無人機(UAV)などの空中車両でレンジ圧縮または方位角圧縮を行えるようにする。
本願発明によると、マルチレートデジタル信号処理技術を新規なやり方で組み合わせ、SAR画像生成の実施の複雑さを最小限とする。
圧縮手段は1以上の圧縮フィルタを有してよい。
周波数多重分離器の出力は、複数の圧縮フィルタの入力に提供されてもよい。周波数多重分離器は周波数チャネルを分離し、各周波数チャネルがそれぞれの圧縮フィルタに提供されてもよい。
他の実施形態において、1以上の圧縮フィルタは、周波数多重分離器へ出力信号を提供する圧縮フィルタを含み、当該圧縮フィルタh'(n)は、周波数チャネルのそれぞれに対してイメージ応答特性を有するフィルタリング機能を実行する。イメージ応答特性を有する圧縮フィルタの周波数応答特性は、H'(ejω)=H(ejωK)であり、ここでH(ejω)はプロトタイプフィルタh(n)の周波数応答特性であり、
Figure 0006066934
であり、ここでKは周波数チャネルの数である。
周波数多重分離器は、圧縮された形式の周波数チャネルを分離する。
周波数多重分離および圧縮の順序を入れ替え、圧縮機能を多重分離機能の前に実行することにより、装置の計算効率がより一層向上する。
周波数多重分離器はブロックチャネル形成器を含んでよい。ブロックチャネル形成器は、ブロック微細処理段およびブロック粗処理段を有し、ブロック微細処理段は、レンジラインまたは方位角ラインに対応する信号に基づき、複数の偶数インデックスチャネルを含み、奇数インデックスチャネルに対応するスペクトルの領域がゼロにされている1以上の第1周波数分割多重(FDM)信号、および複数の奇数インデックスチャネルを含み、偶数インデックスチャネルに対応するスペクトルの領域がゼロにされている1以上の第2FDM信号を生成する。1以上の第1FDM信号および1以上の第2FDM信号は、別々のブロック粗処理段に提供され、ブロック粗処理段のそれぞれは対応するFDM信号に基づき複数の狭帯域幅チャネルを生成してよい。ブロック粗処理段はそれぞれ、対応するFDM信号のスペクトルのゼロにされた領域に遷移域がある応答特性を有するフィルタリング機能を実行する。
本願発明によると、人工衛星またはエアリアルプラットフォーム上でSAR画像を生成する構造は、受信したSARエコーをデジタルサンプルに変換するアナログ−デジタル変換手段と、レンジラインに対応するサンプルのアレイに対してレンジ圧縮を実行するレンジ処理手段、および方位角ラインに対応するサンプルのアレイに対して方位角圧縮を実行する方位角処理手段を有するプロセッサと、サブ画像を組み合わせる画像生成手段とを含み、レンジ処理手段および方位角処理手段のうち少なくとも一方が上述した装置を含み、プロセッサは、レンジラインまたは方位角ラインが対応する信号が多重分離された周波数チャネルに対応する複数のサブ画像を生成する。
レンジ処理手段および方位角処理手段の両方が、周波数多重分離器および圧縮手段を有する上記装置を含んでよい。レンジ処理手段の周波数多重分離器は、レンジラインに対応する信号をKrange個の周波数チャネルに多重分離し、方位角処理手段の周波数多重分離器は、方位角ラインに対応する信号をKazimuth個の周波数チャネルに多重分離し、プロセッサは、Krangeazimuth個のサブ画像を生成する。
画像生成手段は、サブ画像を平均化し、SAR画像を生成してよい。未加工データから生成された、マルチルックされ平均化された画像を格納するのに必要なメモリのボリュームは、典型的には、未加工の高解像度データを格納するのに必要なメモリのボリュームよりもはるかに小さい。サブ画像を平均化することにより、地上へ送信されるデータを格納するのに必要なメモリははるかに小さくて済み、より多くの画像のデータが地上への送信と送信との合間に格納することが出来る。加えて、衛星で画像生成を実行することにより、画像化された目標物の、リアルタイムの、または少なくとも迅速な分析を実行することが出来る。目標物をモニタリングし、衛星で生成された画像に変化が生じたことに応答して警告が発せられてもよい。
受信器は2以上のSARアンテナからエコーを受信し、画像生成手段は、サブ画像に基づきインターフェログラムまたは複素コヒーレンスマップを生成してよい。
本願発明によると、衛星は上述した装置または構造を含む。衛星は処理されたSARデータをダウンリンク送信する手段も含んでよい。
本願発明によると、人工衛星またはエアリアルプラットフォームでSAR画像を生成する方法は、SARアンテナから送信されたチャープのエコーを受信する工程と、エコーをデジタルサンプルにデジタル化する工程と、レンジラインに対応するサンプルのアレイに対してレンジ圧縮を実行する工程と、方位角ラインに対応するサンプルのアレイに対して方位角圧縮を実行する工程とを含み、レンジ圧縮を実行する工程および方位角圧縮を実行する工程のうち少なくとも一方は、複数の目標点の情報を含むアレイにそれぞれが対応する信号を、複数の目標点の情報をそれぞれが含む複数の周波数チャネルに周波数多重分離する工程、およびレンジ圧縮および方位角圧縮が複数の周波数チャネルに対応する複数のサブ画像を生成するよう各チャネルに対して画像形成を実行する工程を有し、方法はさらに、複数のサブ画像を組み合わせる工程を含む。
周波数多重分離する工程および画像形成を実行する工程は、ハードウェアで実行されてもよい。
信号を周波数多重分離する工程は、各周波数チャネルに対して画像形成を実行する工程の前に実行されてもよい。
他の実施形態において、各チャネルに対して画像形成を実行する工程は、信号を分離する工程の前に実行され、周波数多重分離する工程は、圧縮された形式の周波数チャネルを周波数多重分離する工程を含む。
レンジ圧縮を実行する工程は、レンジラインを複数のチャネルに周波数多重分離し、複数のチャネルを圧縮フィルタを用いて個別にレンジ圧縮し、複数の中間行列を形成する工程を有し、方位角圧縮を実行する工程は、各中間行列の各方位角ラインを複数の周波数チャネルへ周波数多重分離し、各周波数チャネルを圧縮フィルタを用いて個別に方位角圧縮し、複数のサブ画像を形成する工程を有する。
他の実施形態において、レンジ圧縮実行する工程は、レンジラインを複数のチャネルに周波数多重分離し、複数のチャネルを圧縮フィルタを用いて個別にレンジ圧縮し、複数の中間行列を形成する工程を有し、方位角圧縮を実行する工程は、各方位角ラインのサンプルを平均化する、または他のアンフォーカスSAR処理を行い、複数のサブ画像を形成する工程を有する。
複数のサブ画像を組み合わせる工程は、複数のサブ画像を平均化し、圧縮画像を生成する工程を有してよい。
SARアンテナから送信されたチャープのエコーを受信する工程は、第1アンテナおよび第2アンテナからエコーを受信する工程を有し、複数のサブ画像を組み合わせる工程は、第1アンテナおよび第2アンテナからの対応するサブ画像を関連付け、複数のサブインターフェログラムを形成する工程と、複数のサブインターフェログラムを平均化し、圧縮されたインターフェログラムを生成する工程とを有する。
添付の図面を参照し本願発明の実施形態を例示として説明する。
図1は、SARアンテナ構造を備える宇宙探査機のペイロードを示す概略図である。 図2は、SARアンテナ構造のデジタルプロセッサのコンポーネントを示す概略図である。 図3は、図2に示すデジタルプロセッサの処理部を示す概略図である。 図4は、図3の処理部で処理される信号のサブチャネルスタック構造を示す概略図である。 図5は、図3の処理部のブロックチャネル形成器を示す機能図である。 図6は、図2の処理部の代替的な実施形態を示す。 図7は、SARデータを受信し処理する方法を示す。 図8は、図3の処理部でレンジラインまたは方位角ラインに対応する信号を処理する方法を示す。 図9は、図6の処理部でレンジラインまたは方位角ラインに対応する信号を処理する方法を示す。
図1は、合成開口レーダ(SAR)アンテナ構造2、メモリ3、および地上の制御局(図示せず)から命令を受信し、データを送信する送受信器4を備える衛星ペイロード1を示す。衛星ペイロードは、SARアンテナ構造2、メモリ3、および送受信器4を制御する制御器5も備える。制御器5はさらに、地上の制御局から命令を受信してもよい。SAR構造2は、パルスの送信およびエコーの受信の両方に用いることが出来る。図1は、受信チェーンにおけるコンポーネントを示す。図1に示すように、SAR構造は、アンテナ6、低ノイズ増幅器(LNA)7、受信したSARエコーをベースバンド周波数にダウンコンバートするダウンコンバータ8、アナログ−デジタルコンバータ(ADC)9、およびデジタル信号プロセッサ10を備える。
アンテナ6がレーダパルスを送信するのにも用いられた場合、構造2は、デジタル−アナログコンバータ、アップコンバータ、および電力増幅器を含む送信チェーンも備える。なお、ADCおよびDACコンバータの十分に高速な動作のためには、ダウンコンバータ8およびアップコンバータは用いられなくてもよい。
いくつかの実施形態によると、衛星ペイロード1は2以上のSARアンテナを含む。この場合、別個のノイズ増幅器、ダウンコンバータ、およびADCを含む別個の受信チェーンがアンテナ毎に存在し得る。アンテナがレーダパルスを送信するのにも用いられる場合、別個の送信チェーンがアンテナ毎に存在し得る。
アンテナは広い範囲にパルスビームを照射する。各パルスは、経時的に変化する周波数を有し得る。時間tにおけるパルスまたは「チャープ」の瞬時RF周波数は、f+Kchirptである。ここでfは中心周波数の、または搬送波周波数であり、Kchirpは、パルスまたはチャープの周波数の変化率である。搬送波からの最大瞬時周波数オフセットは、Kchirpchirpであり、ここでTchirpは、パルス持続時間である。パルスの帯域幅Bchirpは、最大瞬時周波数オフセットにほぼ等しいものとして見なすことが出来る。
chirp≒Kchirpchirp 数式1
一例として、各チャープは、およそ100MHzの帯域幅Bchirpを有する。
衛星は、アロングトラックまたは「方位角」方向に移動し、パルスは、アクロストラックまたは「レンジ」方向に伝搬する。各パルスは目標領域またはスワスを照射し、送信されたパルス毎に、アンテナ構造2はパルスにより照射されたスワス内の散乱体からの反射を受信し処理する。アンテナ6が受信したエコーは、チャープの帯域幅に関連する周波数でサンプリングされる。特定のパルスの受信されたエコーのサンプルは、「レンジライン」と呼ばれるサンプルの1次元のアレイを形成する。いくつかの実施形態において、各レンジラインのサンプルの典型的な数は、20,000である。しかし、レンジラインのサンプルの数は、適用例毎に大幅に異なり得る。
レンジラインの各サンプルは、アンテナと目標物との間の往復距離と比例するインデックスに関連付けられる。パルスのサンプリングされたエコーを含むレンジラインが、レンジ方向に対応する1つの軸と、方位角方向に対応する1つの軸とを有するデータ行列を形成する各パルスに関して生成される。同じインデックスを有するサンプルによって形成されるデータアレイは、「方位角ライン」と呼ばれる。いくつかの実施形態において、方位角ラインのサンプルの典型的な数は、10,000である。しかし、方位角ラインのサンプルの数は、適用例毎に大幅に異なり得る。
信号プロセッサ10は、「レンジドップラー」アルゴリズムとして知られるアルゴリズムに基づきチャープのエコーを処理する。レベル0のSARデータは、信号プロセッサ10によりレベル1の画像データへと処理される。処理は、エコー行列内の点の、SAR画像内の対応する点へのマッピングを含む。より詳細には、SARレーダは、アンテナと目標点との間の距離が相対的な移動に合わせて変化するという事実を利用する。目標物の特定の点からの信号エネルギーがレンジおよび方位角に拡散させられ、SARデータ処理によってこの分散されたエネルギーが出力画像の単一の画素へと集約される。レンジにおいて、目標物の特定の点からの信号が送信されたパルスの持続時間によって拡散させられ、方位角において信号は、目標点がアンテナビームによって照射される時間によって拡散させられる。
図2は、信号を複素形式に変換するプリプロセッサ11、レンジ方向のデータ行列を画像処理するレンジプロセッサ12、レンジ方向の処理に続き行列内のデータを調整するレンジマイグレーション部13、方位角方向のデータ行列を画像処理する方位角プロセッサ14、および画像を生成する画像生成部15を備える信号プロセッサ10を示す。
データの画像処理には、エコー信号の、元の送信された信号との比較が含まれる。このことは、エコー信号と送信されたパルスとの畳み込みによって実行される。理想的には、エコー行列内の各点または画素がSAR画像内の対応する画素にマッピングされる2次元畳み込みを実行して画像を生成する。しかしこれは処理量が多過ぎるので、実際にはレンジ圧縮および方位角圧縮は典型的には、新たな中間行列を蓄積するレンジプロセッサ12によるレンジラインの1式の1次元畳み込み、およびその後の、方位角プロセッサ14による方位角ラインの1式の1次元畳み込みへと分けられ、画像に対応する最終的な行列が形成される。
レンジラインはデジタル化されてプリプロセッサ11を介してレンジプロセッサ12へ提供される。プリプロセッサ11はデジタルアンチエイリアシングフィルタ(DAAF)を備える。各レンジラインの未加工データが実デジタル信号としてDAAFに入力され、DAAFは、受信した実デジタル信号に対して動作し、当該信号を複素ベースバンド形式に変換し、レンジラインをレンジプロセッサ12に渡す。またDAAFは、サンプリングレートを2の倍数でデシメートする。レンジプロセッサ12および方位角処理部内の信号のサンプリング周波数は以降、fで示す。従って、DAAFへの入力信号のサンプリング周波数は2fとなる。
レンジプロセッサ12は各レンジラインに対して動作し、未加工データをレンジ圧縮する。以下により詳細に説明するように、レンジプロセッサ12は各レンジラインの1次元畳み込みを実行し、メモリ3内に1以上の中間行列を蓄積する。
レンジ圧縮された出力は歪んでいることが多く、中間行列内の全てのレンジ行が埋められると、レンジマイグレーション部13はメモリ内の中間行列を読み出し、方位角処理を行うべく正しい列にサンプルを整列させる。当業者には理解いただけるように、レンジマイグレーション機能は、メモリ内のサンプルに対して動作する適切なアルゴリズムを用いて実行される。その後、方位角プロセッサ14は、調整された1以上の中間行列をメモリ3から読み出し、以下により詳細に説明するように方位角ラインの1次元畳み込みを実行し、最終的な行列をメモリ3内に格納する。以下により詳細に説明するように、画像生成部15は、最終的な行列に対しさらに処理を実行してもよい。
畳み込みは異なる技術を用いて実施することも可能である。例えば、直接的な畳み込みを用いてもよい。さらに、高速フーリエ変換(FFT)技術を用いた高速畳み込み、または機能的に同等のいわゆるSPECAN(スペクトル解析)方法も用いることが出来る。本願発明の実施形態によると、デジタルプロセッサのいくつかまたは全てのコンポーネントがハードウェアで実装される。当業者であれば、DAAFおよびレンジマイグレーション機能をハードウェアで容易に実装できることを理解されるであろう。当業者には理解いただけるように、DAAFは理想的にはハードウェアに適している、高速かつ非常に単純な機能である。レンジマイグレーションアルゴリズムは、適用例に応じてハードウェアおよびソフトウェアのいずれかで実装され得る。
図5を参照してより詳細に説明するように、本願発明によると、レンジ圧縮部12および方位角処理部14の一方または両方が、ハードウェアで実装される。加えて、いくつかの実施形態において、画像生成部15もハードウェアを用いて実装される。他の実施形態において、画像生成部15はソフトウェア機能である。制御器5はソフトウェアで実装される制御機能を有する小型の汎用CPUであってよく、全てではないがデジタルプロセッサのコンポーネントの1以上が、ソフトウェア制御機能の一部として設けられる。
図3は、レンジ処理部12または方位角処理部14のコンポーネントを示す。レンジ処理部12および方位角処理部14は互いに対応するコンポーネントを有しているが、それらコンポーネントは、レンジラインまたは方位角ラインに対して動作するかによって異なる構成を有する。以下、図3の説明においてはレンジ処理部としての処理部を説明し、当該処理部がレンジラインに対しどのように動作するかを説明する。しかし、対応する処理段が、異なる畳み込み関数を用いてではあるが、方位角処理部において各方位角ラインに対し実行され得る。
本願発明の実施形態によると、処理部12、14は、データをK個のサブバンドへ周波数多重分離し、各サブバンドのサブ画像を作成し、複数の当該サブ画像を組み合わせる、より効率的な処理を提供する。サブバンドは本明細書において、「周波数チャネル」とも呼ぶ。全帯域幅信号から多重分離された各周波数チャネルは、全帯域幅信号よりも低い解像度ではあるが、目標領域の全ての散乱体の情報を含む。
図3を参照すると、レンジプロセッサ構造はハードウェアで実装され、周波数チャネル形成を行い、デシメート機能を有するブロックプロセッサ16、および複数のレンジ圧縮フィルタ17a〜17kを備える。帯域幅をより小さなチャネルに分けることにより、ハードウェアで容易に実装可能な直接的な畳み込みフィルタリングを各サブチャネルに用いることが出来る。
複素信号が信号プロセッサのDAAF11からブロックプロセッサ16へ提供される。ブロックプロセッサ16において、複数の異なる周波数帯域に対応する複数のサブバンドへ複素信号が多重分離され、各サブバンドのサンプリング周波数が低下させられる。図4は、レンジラインの場合にはチャープ帯域幅Bchirpであるスペクトル全帯域幅信号を示す。スペクトル全帯域幅信号は、K個の等しく離間するサブバンドチャネルからなる、ω=2πk/K(k=0,...K−lであり、Kは偶数である)を中心とした均等なスタック構造を有する周波数多重信号である。言い換えると、レンジラインに関しては、信号が全チャープ帯域幅の複数の周波数サブバンドへ多重分離される。例えばKは100に等しい。しかし、Kは任意の適切な値であってもよい。図5を参照してより詳細に説明するブロックチャネル形成器16は、効率的に計算を行う、構成要素である対応するK個のサブバンドチャネルチャネルに入力を分割するFFTベースのフィルタリングプロセッサである。各サブバンドのサンプリングレートは、ブロックチャネル形成器においてMの倍数で減少させられる。ここでMはチャネル数以下の数である。なお、M=Kの場合、最大のデシメートが行われ、サブバンドに対し臨界サンプリングが行われる。
その後、各サブバンドチャネルはfs/Mのサンプリングレートで別個の圧縮フィルタ17a〜17kへ提供され、圧縮フィルタ17a〜17kは各サブバンドの直接的な畳み込みを実行する。言い換えると、信号の圧縮が各サブバンドに対し別々に実行される。
なお、レンジラインに対応する全帯域幅信号は、複数の目標点の情報を含み、各周波数サブバンドは、当該複数の目標点の情報をそれぞれ含む。言い換えると、各周波数サブバンドは、同じ目標点の情報を含むが、その解像度は全帯域幅信号よりも低い。
全帯域幅ではなく複数のサブバンドに対しレンジ圧縮まはた方位角圧縮を実行することにより、フィルタ17a〜17kの処理速度が遅くてもよいこととなる。全帯域幅圧縮フィルタの長さをNとすると、このフィルタによる入力の直接的な畳み込みにおける計算速度は1秒あたりNfの複素乗算である。 計算速度no channelisaztion=Nf 数式2
いくつかの実施形態において、圧縮フィルタ17a〜17kは、それぞれが自身のサブチャネルに対して動作する同一のインパルス応答フィルタである。各フィルタのフィルタインパルス応答は、送信されたチャープ波形の複素共役であり、フィルタNの長さは、チャープ期間Tchirpのチャープ波形サンプルの数である。結果的に、送信されたチャープ波形が最高レートfでサンプリングされた場合、
N=Tchirp/(1/f) 数式3
となる。
レンジ処理における全帯域幅圧縮フィルタは、帯域幅Bchirpに亘り、上述した数式1を参照して言及したようにBchirp≒Kchirpchirpで表される。
臨界サンプリングされたサブバンドはBchirp/Kの帯域幅を有する。結果的に、変化率KchirpでチャープがサブバンドTsub−band,chirpをスイープするのに必要な時間は、Tchirp/K秒である。さらに、サブバンドはより小さな帯域幅を有するので、サンプリングレートは図3に示すようにf/M(M≦K)まで低下させられ得る。結果的に、サブバンドチャネルのサンプリング時間工程はM/f秒である。当業者には理解いただけるように、M=Kの場合、サンプリングレートは帯域幅に対応させられ、サブバンドに対し臨界サンプリングが行われる。反対に、M<K、例えばK/2の場合、サブバンドはオーバーサンプリングされる。
サブバンドフィルタ圧縮フィルタの長さをNsubとすると、サブバンドフィルタの長さは、サンプリング時間工程において分割されたサブバンドをチャープがスイープするのに必要な時間に等しくなる。
Figure 0006066934
数式4
これは数式3を用いると、
Figure 0006066934
数式5
となる。
サブバンドが臨界サンプリングされると、M=Kであるので、Nsub=N/Kとなる。サブバンドが例えばM=K/2など2の倍数でオーバーサンプリングされた場合、フィルタはより長く、例えばNsub=2N/Kとなる。
一連のサブバンド圧縮フィルタの計算速度は、フィルタ数のKに各サブバンドの速度f/Mおよび各フィルタの長さNsubを乗算したものである。
計算速度channellsed=K(f/M)NsubK(f/M)N/(MK)=Nf/M 数式6
臨界サンプリングされた場合、計算速度は、
計算速度channellsed−Nf/K 数式7
と表すことが出来る。
結果的に、サブバンドに対して動作する圧縮フィルタの畳み込みの計算負荷は、M=Kの場合、全帯域幅信号の直接的な畳み込みの計算負荷の1/Kの倍数である。
全てのレンジラインが処理されると、各チャネルの完全な中間行列がメモリ3内で生成されている。レンジマイグレーション部13において行列が調整されると、中間行列の各方位角ラインが方位角プロセッサ14へ供給され個別に処理される。レンジ圧縮器により生成された中間行列は順次処理されてもよい。上述したように、方位角プロセッサ14も図3に示すように実装されてもよい。各中間行列の各方位角ラインは、K'個のチャネルにチャネル形成されてもよく、K'個のチャネルのそれぞれに対して個別に実行される方位角圧縮は、各中間行列のK'個の新たな行列を作成する。方位角処理部において、ブロックチャネル形成器16は典型的にはデータを異なる数の複数のチャネルに分割し、圧縮フィルタは異なるインパルス応答特性に対応する複数の異なる係数を有する。さらに、方位角処理部において、周波数帯域幅はドップラー帯域幅であり、各パルスの帯域幅ではない。SARアンテナは目標点と相対的に移動するので、目標点までの距離は変化し、このことにより、方位角の関数として受信信号の位相変化が引き起こされる。方位角ラインの隣接するサンプル間の位相変化は、パルス繰り返し数および目標物と相対的にアンテナがどれくらいの速さで移動するかに依存する。位相変化によって生じる方位角の帯域幅は、ドップラー帯域幅として知られている。従って、方位角ラインに関し、信号は受信したエコーのドップラー帯域幅の複数の周波数サブバンドへ多重分離される。なお、方位角ラインに対応する全帯域幅信号は、複数の目標点の情報を含み、信号が多重分離される各周波数サブバンドは、当該複数の目標点の情報をそれぞれ含む。言い換えると、各周波数サブバンドは、同じ目標点の情報を含むが、その解像度は全帯域幅信号よりも低い。
全方位角帯域幅圧縮フィルタもチャープ波の形態をとり得るが、レンジに関してと同じように、異なる中心周波数およびチャープレートを有する。結果的に、周波数チャネルに対する方位角圧縮における計算量の節約も、全帯域幅信号の直接的な畳み込みとは対照的に、1/K'と比例する。ここでK'は方位角帯域幅が分割されるチャネル数である。しかし典型的には、方位角帯域幅はレンジ帯域幅よりもはるかに小さく、方位角圧縮の際に用いられるチャネルの数は、典型的には、レンジ処理の際に用いられるチャネルの数よりもはるかに少ない。
レンジ圧縮および方位角圧縮が完了すると、元の未加工データ行列はメモリ3内で複数のサブ画像に変換されている。全てのサブ画像は同じ目標領域に関するものである。レンジ処理部および方位角処理部の両方がデータを複数の周波数チャネルに分割する実施形態において、サブ画像の数はKK'である。ここでKはレンジ処理の際に用いられるチャネルの数であり、K'は方位角処理の際に用いられるチャネルの数である。
なお、サブ画像はそれ自体が完全に形成されたSAR画像であり、未加工データがその全帯域幅で処理された場合に得られるであろう画像よりも解像度が低いだけである。
再び図2を参照すると、画像生成器15はサブ画像に基づきSAR画像を形成する。いくつかの実施形態において、画像生成器15は複数のサブ画像を組み合わせて全帯域幅の出力信号を生成する。しかし他の実施形態において、画像生成器15は、平均化により複数のサブ画像を組み合わせて、サブ画像と比較し向上した信号対ノイズ比を有する低解像度の画像を生成する。画像は単純に画像を重ね合わせて、対応する画素を追加することにより平均化することが出来る。KK'=100の場合、格納する必要があるデータの量は、画像を組み合わせて完全な解像度の画像を得る代わりに、サブ画像を平均化することにより100の倍数で減少させられる。いくつかの実施形態において、複数の画像からなるサブセットを組み合わせ、複数の異なるサブセットを平均化し、より高い解像度の画像を形成する。平均化されたレベル1のプロダクトを得ることが出来るレベル0のデータを格納するのに必要な記憶ボリュームと比較し、平均化されたレベル1のプロダクトを格納するのに必要な記憶ボリュームははるかに小さい。
いくつかの実施形態において、画像生成器はさらに、レベル1のプロダクトの一部としてインターフェログラムを生成する。宇宙探査機が2以上のSARアンテナを有する場合、複数の異なるアンテナからのサブ画像を画像生成器15内で関連付け、複数のサブインターフェログラムを生成することが出来る。信号対ノイズ比が向上したインターフェログラムは、サブインターフェログラムを平均化することにより生成することが出来る。
さらにいくつかの実施形態において、さらにデータ量を減らすべく画像生成器は、標準的な画像圧縮技術を平均化されたレベル1のプロダクトに対し適用してもよい。一例として、jpeg圧縮技術を用いてさらにデータを圧縮することが出来る。しかし、適切な圧縮技術のいずれを用いてもよい。平均化されたSAR画像または平均化されたインターフェログラムはさらに、格納および/またはダウンリンク送信の前に、標準的な画像圧縮技術を用いて圧縮することが出来る。
図5は、図3のブロックプロセッサ16の適切なアーキテクチャを示す。要約すると、ブロックプロセッサ16は、ブロック微細処理段18a、18b、およびブロック粗処理段19a、19bを有する。ブロックプロセッサおよびブロック処理段は、それらが入力信号の全ての周波数チャネルに対し同時に単一の計算関数を適用するという意味で、「ブロック」プロセッサおよび「ブロック」処理段である。ブロック微細処理段18a、18b、およびブロック粗処理段19a、19bは協働してレンジラインまたは方位角ラインに対応する信号を複数の周波数チャネルに多重分離する。図5に示す実施形態において、ブロック微細処理段は2つのイメージハーフバンドフィルタ18a、18bによって実行される。ハーフバンドフィルタ18a、18bは、レンジラインまたは方位角ラインに対応する信号に基づき、複数の偶数インデックス周波数チャネルを有する周波数分割多重(FDM)信号、および複数の奇数インデックス周波数チャネルを有するFDM信号を生成する。奇数インデックスチャネルに対応するスペクトルの領域は偶数インデックスチャネルを有する信号においてゼロにされており、偶数インデックスチャネルに対応するスペクトルの領域は、奇数インデックスチャネルを有する信号においてゼロにされている。ブロックチャネル形成器のハーフバンドフィルタは、1つの高域フィルタおよび1つの低域フィルタを有し得る。ハーフバンドフィルタのインパルス応答は、N−1個のゼロが埋められ、画像のコピーのオリジナルである基礎となる詰まったプロトタイプの周波数応答のN個の偶数の画像およびN個の奇数の画像を生成する。
そしてブロック粗処理段19a、19bは微細処理段で出力されたFDM信号から偶数のチャネルおよび奇数のチャネルを分離する。より詳細には、各ブロック処理段はWeighted OverLap Add(WOLA)ユニット20a、20bおよびリンクされたK/2ポイントFFTユニット21a、21bを備える。微細処理段からの各FDM信号は、別々のWOLAユニット20a、20b、および対応するリンクされたK/2ポイントFFTユニット21a、21bに提供される。各WOLAユニットおよび関連付けられたK/2ポイントFFTユニットは、それぞれのFDM信号から複数の狭帯域幅チャネルを生成する。各WOLAは、関連付けられたFDM信号のスペクトルのゼロにされた領域に遷移域がある応答特性を有するフィルタリング機能を実行する。いくつかの実施形態において、WOLAユニットおよびリンクされたK/2ポイントFFTユニットはそれぞれ、多相のDiscrete Fourier Transform(DFT)チャネル形成器に置き換えられてもよい。Mの倍数のデシメートは、ブロックチャネル形成器16内で実行される。概念としては、ブロックチャネル形成は、入力の最大のサンプリングレートで起こり、各サブバンド出力は、Mの倍数でオーバーサンプリングされる。そしてサブバンドチャネル出力はエイリアシングを生じることなくMの倍数でデシメートされる。実際には、デシメートの係数は、ブロックチャネル形成器の構成によって吸収され、出力の所望されないサンプルは単純に計算されず、計算量が節約される。
ブロック信号生成器が特定用途集積回路(ASIC)またはフィールドプログラム可能ゲートアレイ(FPGA)に組み込まれてもよい。ブロックプロセッサの適切なアーキテクチャ、およびその修正および変形は、EP0695054およびEP0831611により詳細に説明されている。しかし、ブロックチャネル形成器の適切なアーキテクチャの何れも用いることが出来る。
図6は、レンジ処理部12または方位角処理部14の他の実施形態を示す。図3および6において、同様の参照符号は同様のコンポーネントを指す。しかし図6においては、チャネル毎に別々のフィルタを設ける代わりに、イメージ周波数応答特性を有する圧縮フィルタ17が設けられて信号に対して動作し、その後、信号はブロックプロセッサ16に提供される。圧縮フィルタ17は、図3のプロトタイプフィルタ17a〜17kに基づく。圧縮フィルタの応答は、K−1個のゼロサンプルを各入力サンプル間に埋めることにより補間され、周波数応答H'(ejω)=H(ejωK)が形成される。ここでH(ejω)は、プロトタイプフィルタh(n)の周波数応答であり、
Figure 0006066934
である。
よってイメージフィルタは、全帯域幅に対しては動作せず、各チャネルに対し低解像度圧縮を実行し、ブロック処理で全帯域の信号に対して動作する。続くブロックチャネル形成器16は、既に圧縮された形式のサブバンドを周波数多重分離する。言い換えると、圧縮は周波数チャネルが分離される前に、「概念的な」周波数チャネルに対し実行されるものと見なすことが出来る。イメージ圧縮フィルタはそれでも各サブチャネルに関し、臨界サンプリングされた場合、Nsub個のゼロ以外の係数を有するので、図6の構造の計算のワークロードは、fsub=fN/Kとなり、数式7に示されるように臨界サンプリングされた場合の図3の構造と同じである。
図6のデジタルプロセッサは、サブバンドがオーバーサンプリングされた(例えば2だけ)場合に計算負荷を低減し、続くSAR処理の補助をする。例えば、複数の異なるSARアンテナが用いられた場合、インターフェログラムまたは複素コヒーレンスマップを形成する前に、対応するサブ画像を位置合わせ(co−register)し、補間するのにオーバーサンプリングが必要となり得る。この場合、第1の実施形態における一連の圧縮フィルタのワークロードは、第2の実施形態の単一の圧縮フィルタのワークロードよりも大きくなり、1/Kに対し1/Mとなる。例えばM=K/2の場合、第2の実施形態におけるブロック圧縮は、4倍の効率性となる(1/Kに対し4/K)。
図7は、SAR信号を受信し処理する方法を示す。S7.1において、複数の異なるチャープのエコーがアンテナ6により受信される。そしてS7.2において、エコーがサンプリングされる。S7.3において、信号は、デジタル化されると直ぐに、DAAFフィルタに渡され、複素形式へと処理される。レンジラインが受信され複素形式に変換されると、S7.4において、レンジ処理部12へ渡されレンジ圧縮される。S7.4においてレンジラインが処理されると、中間行列の行が埋められ、複数の完全な中間行列が形成される。なお、S7.1、S7.2、S7.3、およびS7.4の工程は、次の工程が始まる前に終了しなくてもよい。レンジプロセッサは、レンジラインのエコーが受信される際に、同じ画像に関する他のレンジラインを処理してもよい。
ステップS7.5において、方位角処理を行うべく行列内のデータが調整され正しい列へ移動させられ、中間行列の方位角ラインは、ステップS7.6で方位角圧縮され、複数のサブ画像行列が形成される。ステップS7.7において、サブ画像は平均化され、最終的な圧縮された低解像度画像が形成される。地上局への送信が可能となる適切な時間が来るまで圧縮画像はさらに処理され、および/または格納され得る。
図8および図9は、工程S7.4またはS7.6をより詳細に説明する。なおいくつかの実施形態において、図3または6に示すタイプの、信号を複数のチャネルに分割する処理部をレンジ圧縮のみ、または方位角圧縮のみに用い、レンジ圧縮および方位角圧縮のうち一方が全帯域幅に対し実行される、または他の代替的な技術が用いられる。例えば、アンフォーカスSARの方位角圧縮には代替的な技術を用いるのが適切であり得る。しかし他の実施形態においては、データアレイに対応する信号が、レンジ圧縮および方位角圧縮の両方の間に複数の周波数チャネルに分割される。図8および図9を参照すると、図3に示す処理構成でS8.1においてレンジラインまたは方位角ラインに対応するデータ信号が、複数のチャネルへ多重分離され、S8.2において、複数の圧縮フィルタを用いて個々のチャネルに対し圧縮が実行される。または、図6に示す処理構成で、S9.1において信号のチャネルが圧縮され、S9.2において、圧縮されたチャネルが多重分離される。
本願発明の特定の例を説明してきたが、本願発明の範囲は添付の請求項によってのみ定められ、特定の例によって限定されない。よって、当業者には理解いただけるように、本願発明は他のやり方によって実施することが出来る。
例えば上述したように、図3および6を参照して説明したようにレンジ圧縮および方位角圧縮のうち一方のみを実施してもよい。図3および6を参照して説明した処理技術および処理部はレンジラインに関してのみ用いてもよく、方位角圧縮は、各方位角ラインのサンプルを平均化することにより、または、アンフォーカスSARシステムに用いるのに適切な代替的な技術により行ってもよい。いくつかの実施形態において、レンジ圧縮は、図3または6を参照して説明したようにハードウェアで実行され、方位角圧縮はソフトウェアで実行される。方位角処理に関する帯域幅は一般的にはるかに低いので、方位角処理の処理速度は場合によっては、宇宙空間に用いられるCPUのソフトウェアが対応出来る程度の速さであってもよい。結果的に、方位角処理が単純化される適用例において、方位角処理はソフトウェアで実行されてもよい。
なお、信号プロセッサ10のコンポーネントはハードウェアまたはソフトウェアとして実装されるものとして説明してきたが、当該コンポーネントによって実現される機能のいくつかは、ハードウェアとソフトウェアとの組み合わせによって実装することも可能である。
さらに、本願発明は衛星および他の宇宙基地以外にも適用が可能である。本願発明は、適切な適用例のいずれ、特にプロセッサマス、電力、および特にTM帯域幅が重要である適用例において用いることが出来る。例えば、本願発明は無人機(UAV)または他のエアリアルプラットフォームに用いることが出来る。UAVの場合には、リアルタイムで画像化する態様が特に関連性があるであろう。
加えて、図2を参照して、レンジラインがプリプロセッサ11からレンジ処理部12へ直接提供されるものとして説明したが、レンジ処理部12はADC9からのサンプルに対しリアルタイムで動作しなくてもよい。代わりに、ADCの出力は、後に処理を行うべくメモリに格納されてもよい。このようなアプローチは、比較的データ取得/格納段階が短く、時間をおいて処理段が行われるミッションにおいて、ダウンリンク送信の前のデータボリュームを減らすのに有用である。このようなミッションの一例は、惑星の接近通過のミッションである。
また、ブロックチャネル形成器および圧縮フィルタの特定の例を説明してきたが、適切なアーキテクチャのいずれであっても、信号を複数のチャネルに周波数多重分離し、チャネルに対し個別に画像形成を実行するのに用いることが出来る。さらに、圧縮フィルタおよびチャネル形成器の数および順序は、適用例に応じて異なり得る。図6を参照して単一のイメージ応答フィルタを説明したが、チャネルが分離される前にチャネルに対し別々に画像形成を実行するべく複数のイメージ応答フィルタが設けられてもよい。

Claims (14)

  1. 受信した複数の合成開口レーダ(SAR)エコーを複数のデジタルサンプルに変換するアナログ−デジタル変換手段と、
    複数のレンジラインに対応する複数のサンプルの複数のアレイに対してレンジ圧縮を実行するレンジ処理手段、および複数の方位角ラインに対応する複数のサンプルの複数のアレイに対して方位角圧縮を実行する方位角処理手段を有するプロセッサと
    を備え、
    前記レンジ処理手段が、前記SARエコーのデータの1つのレンジラインに対応する領域に含まれる複数の目標点の情報を含む信号を、 range の周波数チャネルへ周波数多重分離する周波数多重分離器と、前記 range の周波数チャネルのそれぞれに対して圧縮を実行する圧縮手段とを有し、各周波数チャネルの信号は前記複数の目標点の情報を含
    前記方位角処理手段が、前記SARエコーのデータの1つの方位角ラインに対応する領域に含まれる複数の目標点の情報を含む信号を、K azimuth 個の周波数チャネルへ周波数多重分離し、K range azimuth 個の周波数チャネルを生成する周波数多重分離器と、前記K range azimuth 個の周波数チャネルのそれぞれに対して圧縮を実行する圧縮手段とを有し、各周波数チャネルの信号は前記複数の目標点の情報を含み、
    前記プロセッサは、前記K range azimuth 個の周波数チャネルに対応するK range azimuth 個のサブ画像を生成し、
    前記K range azimuth 個のサブ画像を組み合わせる画像生成手段をさらに備える、人工衛星またはエアリアルプラットフォーム上で複数のSAR画像を生成するシステム。
  2. 前記レンジ処理手段の前記周波数多重分離器および前記圧縮手段、並びに、前記方位角処理手段の前記周波数多重分離器および前記圧縮手段、の少なくとも一方がハードウェアで実装される、請求項1に記載のシステム
  3. 前記レンジ処理手段の前記圧縮手段および前記方位角処理手段の前記圧縮手段の少なくとも一方が、1以上の圧縮フィルタを有する、請求項1または2に記載のシステム
  4. 前記1以上の圧縮フィルタは、前記周波数多重分離器へ出力信号を提供し、
    応答を提供する前記1以上の圧縮フィルタの周波数応答特性は、H'(ejω)=H(ejωK)であり、ここでH(ejω)はプロトタイプフィルタh(n)の周波数応答特性であり、
    Figure 0006066934
    であり、ここでKは周波数チャネルの数である、請求項3に記載のシステム
  5. 前記レンジ処理手段の前記周波数多重分離器および前記方位角処理手段の前記周波数多重分離器の少なくとも一方はブロック微細処理段およびブロック粗処理段を有し、
    前記ブロック微細処理段は、前記レンジラインまたは方位角に対応する信号に基づき、複数の偶数インデックスチャネルを含み、複数の奇数インデックスチャネルに対応するスペクトルの領域がゼロにされている1以上の第1周波数分割多重(FDM)信号と、前記複数の奇数インデックスチャネルを含み、前記複数の偶数インデックスチャネルに対応するスペクトルの領域がゼロにされている1以上の第2FDM信号とを生成し、
    前記1以上の第1FDM信号および前記1以上の第2FDM信号は、別々の前記ブロック粗処理段に提供され、
    前記ブロック粗処理段のそれぞれは、
    対応するFDM信号に基づき複数の狭帯域幅チャネルを生成し、
    前記対応するFDM信号のスペクトルのゼロにされた領域に遷移域がある応答特性を有するフィルタリング機能を実行する、請求項1から4のいずれか1項に記載のシステム
  6. 前記画像生成手段は、前記 range azimuth のサブ画像を平均化し、SAR画像を生成する、請求項1から5のいずれか1項に記載のシステム。
  7. 2以上のSARアンテナから複数のエコーを受信する受信器をさらに備え、
    前記画像生成手段は、前記 range azimuth のサブ画像からインターフェログラムを生成する、請求項1から6のいずれか1項に記載のシステム。
  8. 請求項1から7のいずれか1項に記載のシステムを備える衛星。
  9. 人工衛星またはエアリアルプラットフォームで合成開口レーダ(SAR)画像を生成する方法であり、
    SARアンテナから送信されたチャープの複数のエコーを受信する工程と、
    前記複数のエコーを複数のデジタルサンプルにデジタル化する工程と、
    複数のレンジラインに対応する複数のサンプルの複数アレイに対してレンジ圧縮を実行する工程と、
    複数の方位角ラインに対応する複数のサンプルの複数アレイに対して方位角圧縮を実行する工程と
    を備え、
    前記レンジ圧縮を実行する工程は前記SARエコーのデータの1つのレンジラインに対応する領域に含まれる複数の目標点の情報を含む信号を、K range の周波数チャネルに周波数多重分離する工程と、前記K range 個の周波数チャネルのそれぞれに対して圧縮を実行する圧縮工程とを有し、各周波数チャネルは前記複数の目標点の情報を含み、
    前記方位角圧縮を実行する工程は、前記SARエコーの前記データの1つの方位角ラインに対応する領域に含まれる複数の目標点の情報を含む信号を、K azimuth 個の周波数チャネルへ周波数多重分離し、K range azimuth 個の周波数チャネルを生成する周波数多重分離する工程と、前記K range azimuth 個の周波数チャネルのそれぞれに対して圧縮を実行する圧縮工程とを有し、各周波数チャネルは前記複数の目標点の情報を含み、
    前記方法はさらに、
    前記K range azimuth 個の周波数チャネルに対応するK range azimuth 個のサブ画像を生成する工程と、
    前記K range azimuth 個のサブ画像を組み合わせる工程と、を備える、方法。
  10. 前記複数の周波数チャネルのそれぞれに対して圧縮を実行する工程が、前記複数の信号を前記周波数多重分離する工程の前に実行され、
    前記周波数多重分離する工程は、圧縮された形式の前記複数の周波数チャネルを周波数多重分離する工程を含む、請求項に記載の方法。
  11. 前記レンジ圧縮を実行する工程は、前記 range の周波数チャネルを圧縮フィルタを用いて個別にレンジ圧縮し、複数の中間行列を形成する工程を有し、
    前記方位角圧縮を実行する工程は、各中間行列の各方位角ラインを前記K azimuth の周波数チャネルへ周波数多重分離し、各周波数チャネルを圧縮フィルタを用いて個別に方位角圧縮し、前記 range azimuth のサブ画像を形成する工程を有する、請求項9または10に記載の方法。
  12. 前記方位角圧縮を実行する工程は、各方位角ラインの前記サンプルを平均化し、前記 range azimuth のサブ画像を形成する工程を有する、請求項9から11のいずれか1項に記載の方法。
  13. 前記 range azimuth のサブ画像を組み合わせる工程は、前記 range azimuth のサブ画像を平均化し、圧縮画像を生成する工程を有する、請求項9から12のいずれか1項に記載の方法。
  14. 前記SARアンテナから送信されたチャープの複数のエコーを受信する工程は、第1アンテナおよび第2アンテナから複数のエコーを受信する工程を有し、前記K range azimuth 個のサブ画像を組み合わせる工程は、前記第1アンテナおよび前記第2アンテナからの対応するサブ画像を関連付け、複数のサブインターフェログラムを形成する工程と、前記複数のサブインターフェログラムを平均化し、圧縮されたインターフェログラムを生成する工程とを有する、請求項9から12のいずれか1項に記載の方法。
JP2013557122A 2011-03-10 2012-03-09 人工衛星またはエアリアルプラットフォーム上で複数のsar画像を生成するシステム、当該システムを備える衛星、人工衛星またはエアリアルプラットフォームで合成開口レーダ(sar)画像を生成する方法 Expired - Fee Related JP6066934B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP11275041 2011-03-10
EP11275041.9 2011-03-10
PCT/EP2012/054164 WO2012120137A1 (en) 2011-03-10 2012-03-09 Sar data processing

Publications (2)

Publication Number Publication Date
JP2014510913A JP2014510913A (ja) 2014-05-01
JP6066934B2 true JP6066934B2 (ja) 2017-01-25

Family

ID=44247017

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013557122A Expired - Fee Related JP6066934B2 (ja) 2011-03-10 2012-03-09 人工衛星またはエアリアルプラットフォーム上で複数のsar画像を生成するシステム、当該システムを備える衛星、人工衛星またはエアリアルプラットフォームで合成開口レーダ(sar)画像を生成する方法

Country Status (5)

Country Link
US (1) US9684071B2 (ja)
EP (1) EP2684072A1 (ja)
JP (1) JP6066934B2 (ja)
CA (1) CA2828923A1 (ja)
WO (1) WO2012120137A1 (ja)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201204792D0 (en) * 2012-03-19 2012-05-02 Qinetiq Ltd Detection techniques
US9442193B2 (en) * 2012-07-12 2016-09-13 Mitsubishi Electric Corporation Radar system and data processing device
RU2526850C2 (ru) * 2012-11-28 2014-08-27 Открытое акционерное общество "Научно-исследовательский институт точных приборов" Способ получения радиолокационного изображения участка земной поверхности и радиолокационная станция с синтезированной апертурой антенны (варианты)
US9103918B2 (en) * 2012-11-29 2015-08-11 Raytheon Company Enhanced radar range resolution
RU2511216C1 (ru) * 2012-12-17 2014-04-10 Открытое акционерное общество "Научно-исследовательский институт приборостроения имени В.В. Тихомирова" Способ формирования изображения поверхности в радиолокационной станции с синтезированием апертуры антенны
US10230925B2 (en) 2014-06-13 2019-03-12 Urthecast Corp. Systems and methods for processing and providing terrestrial and/or space-based earth observation video
WO2016067419A1 (ja) * 2014-10-30 2016-05-06 三菱電機株式会社 合成開口レーダ装置
ES2691496T3 (es) * 2014-11-14 2018-11-27 Airbus Defence and Space GmbH Reducción de datos de recepción de un radar, en particular de un radar de apertura sintética
CA2980920C (en) 2015-03-25 2023-09-26 King Abdulaziz City Of Science And Technology Apparatus and methods for synthetic aperture radar with digital beamforming
EP3311449B1 (en) 2015-06-16 2019-12-11 King Abdulaziz City for Science and Technology Efficient planar phased array antenna assembly
WO2017091747A1 (en) 2015-11-25 2017-06-01 Urthecast Corp. Synthetic aperture radar imaging apparatus and methods
KR101746505B1 (ko) * 2015-11-25 2017-06-13 국방과학연구소 데이터 처리장치, 데이터 처리시스템 및 데이터 처리방법
DE112018000997T5 (de) * 2017-02-24 2019-11-07 Japan Aerospace Exploration Agency Flugkörper und Programm
US10551491B2 (en) * 2017-03-10 2020-02-04 Planet Labs Inc. Systems and methods for performing synthetic aperture radar imaging using superchirps
DE102017205649B3 (de) * 2017-04-03 2018-03-22 Deutsches Zentrum für Luft- und Raumfahrt e.V. Verfahren und Vorrichtung zur rechnergestützten Verarbeitung von SAR-Rohdaten
CA3064735C (en) 2017-05-23 2022-06-21 Urthecast Corp. Synthetic aperture radar imaging apparatus and methods
EP3646054A4 (en) 2017-05-23 2020-10-28 King Abdulaziz City for Science and Technology APERTURE SYNTHESIS RADAR IMAGING APPARATUS FOR MOBILE TARGETS
US11460535B2 (en) * 2017-08-24 2022-10-04 Iceye Oy System and method for transmitting information from synthetic aperture radar satellite to client receiver
US10670711B2 (en) 2017-09-29 2020-06-02 Planet Labs Inc. Systems for synthetic aperture radar transmit and receive antennas
EP3698167A4 (en) 2017-11-22 2021-11-17 Urthecast Corp. SYNTHETIC OPENING RADAR FORMING APPARATUS AND ASSOCIATED PROCESSES
US10690767B2 (en) * 2017-12-11 2020-06-23 Raytheon Company System and method for synthetic aperture radar image formation
CN108562878B (zh) * 2018-03-29 2021-08-31 南京长峰航天电子科技有限公司 一种雷达信号模拟器宽带捷变频测频方法
RU2709623C1 (ru) * 2019-04-30 2019-12-19 Общество с ограниченной ответственностью "Когнитив Роботикс" Способ получения трёхмерного изображения в радаре бокового обзора с синтезированием апертуры антенны
RU2709484C1 (ru) * 2019-04-30 2019-12-18 Общество с ограниченной ответственностью "Когнитив Роботикс" Способ получения трёхмерного изображения в радаре переднебокового обзора с синтезированием апертуры антенны
US11360209B2 (en) 2019-05-03 2022-06-14 Eagle Technology, Llc Multiple beam antenna for wide swath satellite based SAR
US11378681B2 (en) * 2019-10-30 2022-07-05 Nxp Usa, Inc. Systems and methods for automotive synthetic aperture radar
CN110907933B (zh) * 2019-11-26 2022-12-27 西安空间无线电技术研究所 一种基于分布式的综合孔径相关处理系统及方法
EP3913398B1 (en) * 2020-05-20 2022-11-16 Airbus Defence and Space GmbH Synthetic aperture radar data reduction for satellites
WO2022162868A1 (ja) * 2021-01-29 2022-08-04 三菱電機株式会社 物体透視装置、制御回路、記憶媒体、および物体透視方法
JP2022146506A (ja) * 2021-03-22 2022-10-05 三菱重工業株式会社 信号処理装置、信号処理方法及び信号処理プログラム
WO2023233563A1 (ja) * 2022-06-01 2023-12-07 三菱電機株式会社 飛翔体に搭載されるセンサシステム、及びセンサシステム用のデータ処理システム
CN115291218B (zh) * 2022-10-10 2022-12-09 中国电子科技集团公司第十四研究所 一种同源共视多波段干涉sar试验系统

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5646623A (en) 1978-05-15 1997-07-08 Walters; Glenn A. Coherent, frequency multiplexed radar
JPH0682386B2 (ja) 1985-02-22 1994-10-19 株式会社日立製作所 合成開口レーダ画像の焦点化方法
WO1986006177A1 (en) * 1985-04-08 1986-10-23 Hughes Aircraft Company High speed synthetic aperture radar processing system
IT1181998B (it) * 1985-12-24 1987-09-30 Selenia Ind Elettroniche Elaboratore di segnale per radar ad apertura sintetica particolarmente adatto al calcolo parallelo
DE4037725A1 (de) * 1990-11-27 1992-06-11 Deutsche Forsch Luft Raumfahrt Verfahren zur digitalen generierung von sar-bildern und einrichtung zu dessen durchfuehrung
GB9414664D0 (en) 1994-07-20 1994-09-07 British Aerospace Digital signal processing apparatus
GB2317537B (en) 1996-09-19 2000-05-17 Matra Marconi Space Digital signal processing apparatus for frequency demultiplexing or multiplexing
JP2005024395A (ja) 2003-07-02 2005-01-27 Toshiba Corp レーダ装置およびレーダシステム

Also Published As

Publication number Publication date
US9684071B2 (en) 2017-06-20
JP2014510913A (ja) 2014-05-01
EP2684072A1 (en) 2014-01-15
CA2828923A1 (en) 2012-09-13
US20140009326A1 (en) 2014-01-09
WO2012120137A1 (en) 2012-09-13

Similar Documents

Publication Publication Date Title
JP6066934B2 (ja) 人工衛星またはエアリアルプラットフォーム上で複数のsar画像を生成するシステム、当該システムを備える衛星、人工衛星またはエアリアルプラットフォームで合成開口レーダ(sar)画像を生成する方法
JP6632342B2 (ja) 合成開口レーダなどのレーダの受信データの圧縮
EP2743727B1 (en) Method for implementing high-resolution wide-swath spaceborne SAR system
JP5411506B2 (ja) 高分解能合成開口レーダ装置及び高分解能合成開口レーダ装置用アンテナ
US10228459B2 (en) Radar system and radar signal processing device
JPH0727021B2 (ja) 合成開口レーダ装置
CN106932778B (zh) 方位多通道fmcw sar滑动聚束成像方法
US10690767B2 (en) System and method for synthetic aperture radar image formation
CN107544068A (zh) 一种基于频域bp的图像域宽带合成方法
CN103969644A (zh) 一种多通道调频连续波sar成像方法
CN111929681B (zh) 一种基于轻小型无人机载sar的实时成像处理系统
JP5219897B2 (ja) 画像レーダ装置
JP5106323B2 (ja) レーダ画像再生装置
CN105974409B (zh) 基于多频子带并发的星载滑动聚束mimo-sar成像方法
US20170010345A1 (en) Direct Sampling of Received Signals in Radar
WO2013067594A1 (en) An apparatus and a method for obtaining information about at least one target
CN109001733B (zh) 一种geo星机双基合成孔径雷达频谱模糊抑制方法
GB2086169A (en) Method of processing data generated by a synthetic aperture radar system
CN113189547B (zh) 一种基于合成带宽频率变标的sar成像的方法及系统
CN111624600B (zh) 一种基于一比特量化的sar动目标检测方法及系统
US20170269193A1 (en) Signal processing device
He et al. Large-Scene Imaging Scheme for Geosynchronous Spaceborne-Airborne Bistatic SAR
CN110687505A (zh) 一种适用于非零带宽信号的频带采样微波被动观测系统
DE102004045273B4 (de) Verfahren zur Bilderzeugung bei einem Synthetischen Apertur Radar
Luo et al. Imaging for MIMO sliding spotlight SAR using stepped frequency chirps

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150203

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150929

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160607

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20160810

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20161107

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: 20161129

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20161220

R150 Certificate of patent or registration of utility model

Ref document number: 6066934

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees