JP3532924B2 - レーダーによって得られたデータを改善する方法 - Google Patents

レーダーによって得られたデータを改善する方法

Info

Publication number
JP3532924B2
JP3532924B2 JP51598496A JP51598496A JP3532924B2 JP 3532924 B2 JP3532924 B2 JP 3532924B2 JP 51598496 A JP51598496 A JP 51598496A JP 51598496 A JP51598496 A JP 51598496A JP 3532924 B2 JP3532924 B2 JP 3532924B2
Authority
JP
Japan
Prior art keywords
radar
function
frequency
transformed
signal
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
JP51598496A
Other languages
English (en)
Other versions
JPH10508941A (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 JPH10508941A publication Critical patent/JPH10508941A/ja
Application granted granted Critical
Publication of JP3532924B2 publication Critical patent/JP3532924B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • 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
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/288—Coherent receivers
    • G01S7/2883—Coherent receivers using FFT processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Complex Calculations (AREA)

Description

【発明の詳細な説明】 本発明は、レーダーによって得られたデータを改善す
る方法に関するものである。
レーダーに関連して生じ得る問題は、レーダーシグナ
ルの受信がラジオおよびテレビジョンの送信によって妨
害され得るということである。特に、CARABAS（登録商
標）と呼ばれ、スウェーデン特許No.8406007−８（45
6、117）およびヨーロッパ特許86900306.1（202,320）
の発明の要旨をなす広帯域低周波数レーダーにおいて
は、これは極めて重大な問題である。少なくとも局所的
なラジオおよびテレビジョン基地局のシグナル強度は、
その周波数で対応する周波数におけるCARABASのシグナ
ル強度より遥かに強い。さらに、レーダー基地局の出力
をこれらラジオおよびテレビジョンシグナルに対抗する
のに十分な程度まで増大させることは不可能である。さ
らに、これは違法なことである。その代わりに、もし可
能ならばそのような周波数を出力シグナル中において排
除するか、あるいは少なくともこれらの周波数における
受信されたシグナルを無視することが便利である。
レーダーに適用される全周波数インターバルにわたっ
て望まれたシグナル処理を実行することができるように
するため、問題となっているタイプの狭周波数帯干渉
を、フィルタによって除去し、シグナル強度の計算され
た値によって置き換えることが望ましい。本発明は、こ
の問題を解決するものである。
本発明は、添付図面を参照して以下においてより詳細
に説明される。
図１は、「ジャングルロウ」と呼ばれるLabView2にお
いて本発明を実行するときに用いられる第１のサブルー
チンを示した図である。
図２は、「ミューテーション」とよばれるLabView2に
おいて本発明を実行するときに用いられる第２のサブル
ーチンを示した図である。
図３は、「ジェオメトリックマシン」と呼ばれるLabV
iew2において本発明を実行するときに用いられる第３の
サブルーチンを示した図である。
図４は、「ダーウィニスティックリラクゼーション」
と呼ばれるLabView2において本発明を実行するときに用
いられるメインルーチンを示した図である。
本発明は、これまで、上述した広帯域レーダー装置CA
RABASに基づいて発展してきた。そして、以下では、こ
の装置を考慮して主たる説明がなされる。しかしなが
ら、本発明は、それ以外のレーダーに関連して用いられ
ることができ、結果的に、これは請求の範囲に記載され
た構成の範囲内に含まれるものであることが当業者によ
って容易にわかる。
レーダーシグナルが送信されたとき一定の周波数が排
除されたか否かにかかわらず、受信されたレーダーデー
タは多数の帯域に分割される。これら多数の帯域は、互
いに独立しており、また帯域間には、ラジオ送信等から
の干渉が関係するレーダーデータの存在しないギャップ
が生じる。これに関して、SARレーダーにおけるイメー
ジスペクトルと一般的なレーダー未処理データとの間に
は本質的な差異は存在しないが、ここでは、このような
イメージスペクトルによってレーダーデータが意味さ
れ、このレーダデータから実際のイメージが２次元逆フ
ーリエ変換によって得られる。実際、SARイメージは、
送信された周波数およびこれらのデータのドップラー周
波数をそれぞれイメージスペクトルの方位角周波数およ
び動径周波数に関係付けるような処理により、レーダー
未処理データから得られる。ラジオ送信等は、この２次
元スペクトル中において破壊されたデータの同心円とし
て現れる。これら同心円の動径周波数は、干渉が生じる
送信周波数に対応している。
レーダーデータの質の低下に起因してレーダーイメー
ジ中に生じる質の低下を計算するために、レーダーデー
タおよびレーダーイメージの統計に関する仮定がなされ
る。レーダーデータに対して統計はほとんど変動せず、
ランダム位相のスペクトルとみなすことができる。しか
しながら、この場合（レーダーシステムにおける感度揺
らぎの調節の後）、振幅が大きく変動する可能性はほと
んどない。
他方、イメージは非常に不均一な振幅分布によって特
徴付けられており、この振幅分布は、理想的な場合に
は、均一な強度をもつ弱い反射体のバックグラウンドに
対する離散的なより強い反射体から構成されている。こ
の強い反射体はレーダーによって検出されるべき目標物
である。
見失われたポイントにおけるスペクトルの位相値のす
べてにおいて、予備的知識は存在しないので、あらゆる
推測が平等に現実に生じる。本発明においては、したが
って振幅値はこれらのポイントではゼロに設定されてい
る。
レーダーデータは、二つの処理、すなわち、理想的な
レーダーデータに対応する処理と、妨害されたインター
バル中においてこれらのデータを正確に消去するエラー
処理との組み合わせの結果であるとみなされ得る。これ
ら二つの処理は、共にランダム位相を有している。ラン
ダム位相のためにエラー処理と理想的な処理との間の出
力比は、見失われたポイントとすべてのポイントとの間
の比率となる。もし、これらのポイントの半数が見失わ
れるならば、相対的なエラーは−3dBとなる。ここに、
この比率は、REPと略される相対エラー出力と呼ばれ
る。
以下において、ラジオ送信によって妨害された周波数
での受信されたレーダーデータはゼロに設定されるもの
と仮定される。このとき、より強い反射体からのレーダ
反射およびラジオ干渉が混在することは、イメージ全体
にわたって広がり、かつ弱い反射体として誤って解釈さ
れるこれらの反射体からの副ローブに帰着する。しかし
ながら、より強い反射体それ自体に対しては、副ローブ
が合計されて相対的に同等に強い反射体をもたらすこと
は起こりそうもない。
したがって、妨害されたスペクトルインターバル中に
おけるデータを最初に評価する際に、状況は、最も強い
反射体がまさしくそのものであるということ仮定によっ
て改善され得る。この仮定に基づき、妨害されたインタ
ーバル中におけるデータを補間することによって、最も
強い反射体からの副ローブは減少し、その結果、誤った
値を除去するためのスレッショールドは低下し得る。こ
のとき、第２の反復適用が可能となる。この反復適用
は、見失われた周波数におけるスペクトルのより良好な
計算を確立するための強い反射体の広いレンジを含んで
いる。これによって、副ローブレベルがさらに減少す
る。
より弱い反射率を処理する場合、含まれた反射体の個
数は著しく増大する。それ故、誤った目標物が実際の目
標物よりも強くなるという可能性は急激に増大する。も
し、誤った目標物の個数が、反復適用ステップにおいて
現実の目標物の個数を超えるならば、このステップは計
算の質を低下させる。これの解決法は、与えられたステ
ップにおけるスレッショールドを低下されるよりもむし
ろ、与えられた方法でより多くの目標物を含ませること
によって補間を実行することにある。
実験では、２つの級数の累乗が用いられ、そして、最
初のステップでは最も強い反射体だけが、次のステップ
では最も強い二つの反射体が、その次のステップでは４
つの最も強い反射体がというように、反復適用の基礎と
して含まれた。この反復適用は、ピクセルの半数がその
補間中に含まれるようになったときに終了する。この処
理は、指数関数的に進行していくという利点を有してお
り、それによって、非常に膨大な個数のピクセルを処理
するのに適している。それと同時に、SARイメージは非
常に強い散乱体をわずかしか含んでおらず、その一方大
多数の目標物がノイズフロアに対して相対的に小さなコ
ントラストを与えるという事実をうまく考慮する。この
方法は、シミュレートされたデータ上においてうまく機
能することがわかっている。
このように、本発明は、狭周波数帯干渉、すなわち受
信されたレーダー反射の帯域幅よりも小さい帯域幅を有
する干渉によって妨害された周波数帯にわたって補間を
行うことにより、レーダーによって得られたデータを改
善する方法に関するものであり、以下に詳細に説明する
ような構成を備えている。
周知のように、レーダーシグナルはレーダー送信機に
おいて生成され、アンテナシステムを通じて送信され
る。レーダー反射はアンテナシステムを通じてレーダー
受信機において受信され、シグナル処理装置において処
理される。そのとき、次のステップ（Ａ）〜（Ｊ）が実
行される。
（Ａ）モニタリングまたはそれ以外の方法によって、ど
の周波数値が狭周波数帯干渉によって影響を受けたのか
が決定される。
（Ｂ）時間的に連続して接続され、時間離散的な実数ま
たは複素数シグナルｆ（ti）によって表され、Ｎ個のサ
ンプルを含む各レーダー反射が、Ｎ個のサンプルによっ
て与えられた周波数離散的な関数Ｆ（ωｉ）にフーリエ
変換される。なお、ここに、添字ｉは１からＮまでをと
り、ｉ＝１、２、……、Ｎである。
（Ｃ）一組の整数N1、N2、……、NMが、１≦N1＜N2＜…
…＜NM＜Ｎとなるように固定される。
（Ｄ）シグナル処理装置により、狭周波数帯干渉の外側
領域においてFk（ωｉ）＝Ｆ（ωｉ）とおき、狭周波数
帯干渉を有するインターバルの範囲内ではFk（ωｉ）＝
Hk（ωｉ）とおくことによって、Ｆ（ωｉ）がF
k（ωｉ）に変形される。
（Ｅ）以下のステップ（Ｆ）〜（Ｈ）に従って、H1（ω
ｉ）＝０とおかれ、ｋ≧２に対してはHk（ωｉ）が形成
される。
（Ｆ）Fk-1（ωｉ）が時間離散的な関数fk-1（ti）に逆
フーリエ変換される。
（Ｇ）シグナル処理装置により、ｉ＝１、２、……、N
k-1に対してhk（ti）＝fk-1（ti）、ここで|fk-1（ti）
｜は|fk-1（ti）｜のNk-1最大値を構成する、とおき、
関数の残りのＮ−Nk-1個の値に対してはhk（ti）＝０と
おくことによって、fk-1（ti）が関数hk（ti）に変形さ
れる。
（Ｈ）hk（ti）が、周波数離散的な関数Hk（ωｉ）にフ
ーリエ変換される。
（Ｉ）ステップ（Ｄ）〜（Ｈ）がＭ回繰り返される。こ
こで１≦ｋ≦Ｍである。
（Ｊ）FNM（ωｉ）を逆フーリエ変換されることによっ
て、最終的にフィルターをかけられたレーダー反射fNM
（ti）が得られる。
SARレーダーに対して対応するステップは２次元で実
行される。この場合、狭周波数帯干渉は、SARイメージ
のフーリエ変換における妨害されたインターバルとして
生じる。妨害されたインターバルは、２次元スペクトル
における同心円の特性を有しており、すなわち、それら
はインターバルΩ≦（ωi 2＋μj 2）1/2≦Ω’における
妨害されたインターバルとして現れる。ここに、値Ω、
Ω’は各独立した干渉帯域を表している。
ステップ（Ａ）〜（Ｊ）は次の内容を有している。
（Ａ）モニタリングまたはそれ以外の方法によって、ど
の周波数値が狭周波数帯干渉によって影響を受けたのか
が決定される。
（Ｂ）イメージが、Ｎ×Ｎ個のサンプルを含む時間離散
的な実数または複素数シグナルｆ（xi,yj）によって記
述され、シグナルｆ（xi,yj）が、Ｎ×Ｎ個のサンプル
によって与えられた周波数離散的な関数Ｆ（ωi,μｊ）
にフーリエ変換される。なお、ここに、添字ｉおよびｊ
は１からＮまでをとり、ｉ＝１、２、……、Ｎ、ｊ＝
１、２、……、Ｎである。
（Ｃ）一組の整数N1、N2、……、NMが、１≦N1＜N2＜…
…＜NM＜N2となるように固定される。
（Ｄ）シグナル処理装置により、狭周波数帯干渉の外側
領域においてFk（ωi,μｊ）＝Ｆ（ωi,μｊ）とおき、
狭周波数帯干渉を有するインターバルの範囲内ではF
k（ωi,μｊ）＝Hk（ωi,μｊ）とおくことによって、
Ｆ（ωi,μｊ）が関数Fk（ωi,μｊ）に変形される。
（Ｅ）以下のステップ（Ｆ）〜（Ｈ）に従って、H1（ω
i,μｊ）＝０とされ、ｋ≧２に対してはHk（ωi,μｊ）
が形成される。
（Ｆ）Fk-1（ωi,μｊ）が時間離散的な関数fk-1（xi,y
j）に逆フーリエ変換される。
（Ｇ）シグナル処理装置により、ｉ＝１、２、……、N
k-1に対してhk（xi,yj）＝fk-1（xi,yj）、ここで|fk-1
（xi,yj）｜は|fk-1（xi,yj）｜のNk-1最大値を構成す
る、とおき、関数の残りのＮ−Nk-1個の値に対してはhk
（xi,yj）＝０とおくことによって、fk-1（xi,yj）が関
数hk（xi,yj）に変形される。
（Ｈ）hk（xi,yj）が周波数離散的な関数Hk（ωi,
μｊ）にフーリエ変換される。
（Ｉ）ステップ（Ｄ）〜（Ｈ）をＭ回繰り返される。こ
こに、１≦ｋ≦Ｍである。
（Ｊ）FNM（ωi,μｊ）が逆フーリエ変換されることに
よって、最終的にフィルターをかけられたレーダーイメ
ージfNM（xi,yj）が得られる。
上述した本発明の２つの方法のいずれか一方の好まし
い実施例では、ステップ（Ｃ）における数N1、N2、…
…、NMは、Nk＝2k-1、ここでｋは１≦ｋ≦Ｍなる整数で
ある、に従って選択される。
さらには、１次元の場合には、数Ｍを、log2（N/2）
−１に最も近い整数として選択すること、すなわち、NM
≦N/2となるようにすることが便利である。もしＮがＮ
×Ｎによって置き換えられるならば、これと同様のこと
を２次元の場合に適用することができる。
上述のSARの場合において、イメージは、それ自体す
ぐに示唆されるように、Ｎ×Ｎの正方形をなしているも
のと仮定される。しかしながら、その代わりにＮ×N'の
長方形のイメージを用いることも容易に思いつく。長方
形のイメージはより大きな正方形イメージの興味ある一
部とみなされ得る。この長方形イメージにおいて見失わ
れた円はそれイメージの１つの円の一部として見い出さ
れ得る。他方、正方形イメージ中に長方形イメージを埋
め込み、正方形イメージの残りの部分をシグナルに関し
てゼロで満たすことによって長方形イメージから正方形
イメージに移行することが常に可能である。
本発明は、ナショナル インスツルメントからプログ
ラミング言語Ｃに至るグラフィックシェルLabView2によ
って実行され、このプログラムによって生成されたシン
ボルを用いることによって記述される。このプログラム
中に含まれる種々のサブルーチンの名前は、ダーウィン
進化論との外見上の類似によって示唆されている。もち
ろん、現実との関連性は全くない。基礎的なサブルーチ
ンは、この精神から「ジャングルロウ」と呼ばれ、図１
に図示されている。このプログラムは現実の双極子アレ
イに種々の値を保持するか、またはそれらの値をゼロに
設定するルーチンから構成されている。アレイ中のエレ
メントが保持されるか否かを決定する事実は、その値が
スレッショールド値を越えるか否かである。
アレイは、増大する振幅のオーダーにおいてソートさ
れ、そしてランクのオーダーにおいてソートされるべく
反転される。予め決定されたランク数はスレッショール
ド値を設定し、このスレッショールド値以上においてデ
ータ値は不変に保たれる。すなわち、２つの類似アレイ
値が存在しないと仮定することによって、ゼロでない値
の個数はスレッショールド値のランク数に一致する。
次のレベルはプログラム「ミューテーション」であ
る。このプログラムは見失われた周波数データを書き込
むものであり、図２に図示されている。入力値は、周波
数に位相するデータの複素数アレイと、スレッショール
ド値のランク数と、見失われた周波数のアレイである。
データアレイから、イメージアレイは、逆FFT（高速フ
ーリエ変換）の後、実数部分として得られる。「ジャン
グルロウ」は、これに影響を及ぼすことにより、最も強
いイメージアレイ値のみをゼロでないように保持するこ
とを許される。ゼロを虚数部に入力した後、周波数デー
タは順方向FFTによって得られる。値はこのアレイから
拾い上げられ、如何なる周波数データももたないポイン
トとしてリストアップされたアレイ位置に見い出される
べき値に置き換えられ得る。この位置は以下に説明する
ようにプログラム「ダストビン」から得られる。
反復適用においてスレッショールド値のランク数を低
下させるプログラムは、「ジェオメトリックマシン」と
呼ばれ、図３に図示されている。プログラムの入力値は
再び周波数データの複素数アレイからなっている。用い
られるこれらのデータからの情報のみが、アレイの長さ
となる。スレッショールド値のランク数は、シフトレジ
スタから得られる。シフトレジスタは当初１に設定され
ている。プログラムが実行されるとき、シフトレジスタ
内の値は、この値がアレイの長さの半分以下である限
り、数２を掛けられた値によって置き換えられる。シフ
トレジスタ内に生じる値はすべて、スレッショールド値
のランク数のシーケンスを表示する出力アレイに集めら
れる。
最後のステップはメインルーチン「ダーウィニスティ
ックリラクゼーション」であり、図４に図示されてい
る。プログラム「ミューテーション」、「ジェオメトリ
ックマシン」に加えて、このプログラムまたは、プログ
ラム「ダストビン」および「インバージョン」を使用す
る。それらは、現存する形式で簡単に記述される。
この場合、「ダストビン」は、どのアレイ位置におい
て周波数に位相するデータアレイがゼロであるのかをチ
ェックし、それらをアレイとして受け取るサブルーチン
である。「インバージョン」は１次元FFTのみを含んで
おり、適当なパラメータを選択した後、グラフィックデ
ィスプレイに適用されるFFTの実数部分を選び出す。
「ダーウィニスティックリラクゼーション」は、周波
数データの入力アレイによってシフトレジスタを初期化
する。プログラム「ジェオメトリックマシン」によっ
て、「ダーウィニスティックリラクゼーション」は、ス
レッショールド値のランク数を超えるアレイ内の値を読
み取り、プログラム「ミューテーション」により、シフ
トレジスタ内のデータアレイを正確なものとし、そし
て、古いデータアレイを更新されたものによって置き換
える。これはスレッショールド値のランク数を超えるア
レイ内のエレメントのすべてに対して続行される。最後
に、「−１」がランク数として入力され、これはシフト
レジスタのさらにもう１回だけの反復適用を強制し、そ
れによって最終結果が「インバージョンプログラム」に
入力される。
実数のSARデータに対し、種々の方法において、補間
に先立って、どの周波数が妨害されたかを検出すること
ができる。したがって、ラジオ送信等によって妨害され
た周波数の円が、また実際の状況においてゼロに設定さ
れるということを仮定することができる。
干渉がどこで生じているかを決定する簡単な技術は、
レーダー受信機に平行な受信シグナルに単に耳を傾ける
ことである。SARイメージを周波数データのヒストグラ
ムに変換することにより、放送が引き起こす干渉を検出
することも可能である。これらは、受信された振幅を検
出可能なレーダーデータと検出可能な干渉に分解する鋭
いスレッショールド値を生じさせる。
上述したように、LabViewにおいて実行されるフィル
タリングプログラムは、シミュレートされるデータに基
づいてテストされてきた。入力値は、SARイメージのヒ
ストグラムとそれらの１次元の対応物とラジオ送信統計
であった。SARイメージに対して許容され得るREP値は−
10dB以下である。「ダーウィニスティックリラクゼーシ
ョン」によるフィルタリングの結果はREP値の低下を意
味し、もし利用可能な帯域幅の30〜40％がラジオ干渉に
よって妨害されるならば、「ダーウィニスティックリラ
クゼーション」は利用可能なイメージ中において全く同
一の結果を生じさせる。フィルタリングがなければ、帯
域幅の５〜10％だけが見失われ得る。言い換えれば、30
〜40％の妨害時に、REP値はフィルタリングによって低
下せしめられ、フィルタリングなしの５〜10％の妨害時
のREP値と比較できるものとなる。
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平５−249237（ＪＰ，Ａ) 特開 平５−180926（ＪＰ，Ａ) 特開 平４−307802（ＪＰ，Ａ) 米国特許4612545（ＵＳ，Ａ) 米国特許5128679（ＵＳ，Ａ) 特許2815208（ＪＰ，Ｂ２) (58)調査した分野(Int.Cl.7，ＤＢ名) G01S 7/00 - 7/42 G01S 13/00 - 13/95

Claims (6)

    (57)【特許請求の範囲】
  1. 【請求項１】レーダーによって得られたデータを、狭周
    波数帯干渉、すなわち受信されたレーダー反射の帯域幅
    より小さい帯域幅をもつ干渉によって妨害された周波数
    帯にわたって補間することによって改善する方法であっ
    て、レーダー送信機においてレーダーシグナルを生成
    し、前記レーダーシグナルをアンテナシステムを通じて
    送信し、レーダー反射をアンテナシステムを通じてレー
    ダー受信機において受信し、前記レーダー反射をシグナ
    ル処理装置において処理する方法において、 （Ａ）モニタリングまたはそれ以外の方法によって、ど
    の周波数値が前記狭周波数帯干渉によって影響を受けた
    のかを決定し、 （Ｂ）時間的に連続して接続され、時間離散的な実数ま
    たは複素数シグナルｆ（ti）、ここでｉ＝１、２、…
    …、Ｎである、によって表され、Ｎ個のサンプルを含む
    前記各レーダー反射を、Ｎ個のサンプルによって与えら
    れた周波数離散的な関数Ｆ（ωｉ）、ここでｉ＝１、
    ２、……、Ｎである、にフーリエ変換し、 （Ｃ）一組の整数N1、N2、……、NMを、１≦N1＜N2＜…
    …＜NM＜Ｎとなるように固定し、 （Ｄ）前記シグナル処理装置により、前記狭周波数帯干
    渉の外側領域においてFk（ωｉ）＝Ｆ（ωｉ）とおき、
    前記狭周波数帯干渉を有するインターバルの範囲内では
    Fk（ωｉ）＝Hk（ωｉ）とおくことによって、関数Ｆ
    （ωｉ）を関数Fk（ωｉ）に変形し、 （Ｅ）以下のステップ（Ｆ）〜（Ｈ）に従って、H1（ω
    ｉ）＝０とし、ｋ≧２に対してはHk（ωｉ）を形成し、 （Ｆ）関数Fk-1（ωｉ）を時間離散的な関数fk-1（ti）
    に逆フーリエ変換し、 （Ｇ）前記シグナル処理装置により、ｉ＝１、２、…
    …、Nk-1に対してhk（ti）＝fk-1（ti）、ここで|fk-1
    （ti）｜は|fk-1（ti）｜のNk-1最大値を構成する、と
    おき、関数の残りのＮ−Nk-1個の値に対してはhk（ti）
    ＝０とおくことによって、関数fk-1（ti）を関数h
    k（ti）に変形し、 （Ｈ）関数hk（ti）を、周波数離散的な関数Hk（ωｉ）
    にフーリエ変換し、 （Ｉ）前記ステップ（Ｄ）〜（Ｈ）をＭ回繰り返し、こ
    こで１≦ｋ≦Ｍとし、 （Ｊ）FNM（ωｉ）を逆フーリエ変換することによっ
    て、最終的にフィルターをかけられた前記レーダー反射
    fNM（ti）を得ることを特徴とする方法。
  2. 【請求項２】前記ステップ（Ｃ）における数N1、N2、…
    …、NMは、Nk＝2k-1、ここでｋは１≦ｋ≦Ｍなる整数で
    ある、に従って選択されることを特徴とする請求項１に
    記載の方法。
  3. 【請求項３】前記数Ｍは、log2（N/2）−１に最も近い
    整数として選択されることを特徴とする請求項２に記載
    の方法。
  4. 【請求項４】SARレーダーにおけるレーダーイメージ
    を、インターバルΩ≦（ωi 2＋μj 2）1/2≦Ω’におけ
    る妨害されたインターバル中に見い出されるべき周波数
    帯にわたって補間することによって改善する方法であっ
    て、レーダー送信機においてレーダーシグナルを生成
    し、前記レーダーシグナルをアンテナシステムを通じて
    送信し、レーダー反射をアンテナシステムを通じてレー
    ダー受信機において受信し、前記レーダー反射をシグナ
    ル処理装置において処理する方法において、 （Ａ）モニタリングまたはそれ以外の方法によって、ど
    の周波数値が前記狭周波数帯干渉によって影響を受けた
    のかを決定し、 （Ｂ）前記イメージを、N'×N"個のサンプルを含む時間
    離散的な実数または複素数シグナルｆ（xi,yj）、ここ
    でｉ＝１、２、……、N'、ｊ＝１、２、……、N"であ
    る、によって記述し、前記シグナルｆ（xi,yj）を、N'
    ×N"個のサンプルによって与えられた周波数離散的な関
    数Ｆ（ωi,μｊ）、ここでｉ＝１、２、……、N'、ｊ＝
    １、２、……、N"である、にフーリエ変換し、 （Ｃ）一組の整数N1、N2、……、NMを、１≦N1＜N2＜…
    …＜NM＜Ｎ′×N"となるように固定し、 （Ｄ）前記シグナル処理装置により、前記狭周波数帯干
    渉の外側領域においてFk（ωi,μｊ）＝Ｆ（ωi,μｊ）
    とおき、前記狭周波数帯干渉を有するインターバルの範
    囲内ではFk（ωi,μｊ）＝Hk（ωi,μｊ）とおくことに
    よって、関数Ｆ（ωi,μｊ）を関数Fk（ωi,μｊ）に変
    形し、 （Ｅ）以下のステップ（Ｆ）〜（Ｈ）に従って、H1（ω
    i,μｊ）＝０とし、ｋ≧２に対してはHk（ωi,μｊ）を
    形成し、 （Ｆ）関数Fk-1（ωi,μｊ）を時間離散的な関数f
    k-1（xi,yj）に逆フーリエ変換し、 （Ｇ）前記シグナル処理装置により、ｉ＝１、２、…
    …、Nk-1に対してhk（xi,yj）＝fk-1（xi,yj）、ここで
    |fk-1（xi,yj）｜は|fk-1（xi,yj）｜のNk-1最大値を構
    成する、とおき、関数の残りのＮ−Nk-1個の値に対して
    はhk（xi,yj）＝０とおくことによって、関数fk-1（xi,
    yj）を関数hk（xi,yj）に変形し、 （Ｈ）関数hk（xi,yj）を、周波数離散的な関数H
    k（ωi,μｊ）にフーリエ変換し、 （Ｉ）前記ステップ（Ｄ）〜（Ｈ）をＭ回繰り返し、こ
    こで１≦ｋ≦Ｍとし、 （Ｊ）FNM（ωi,μｊ）を逆フーリエ変換することによ
    って、最終的にフィルターをかけられた前記レーダーイ
    メージfNM（xi,yj）を得ることを特徴とする方法。
  5. 【請求項５】前記ステップ（Ｃ）における数N1、N2、…
    …、NMは、Nk＝2k-1、ここでｋは１≦ｋ≦Ｍなる整数で
    ある、に従って選択されることを特徴とする請求項４に
    記載の方法。
  6. 【請求項６】前記数Ｍは、log2（N'×N"/2）−１に最も
    近い整数として選択されることを特徴とする請求項５に
    記載の方法。
JP51598496A 1994-11-10 1995-11-10 レーダーによって得られたデータを改善する方法 Expired - Fee Related JP3532924B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
SE9403872A SE503650C2 (sv) 1994-11-10 1994-11-10 Sätt att vid en radar förbättra radardata
SE9403872-6 1994-11-10
PCT/SE1995/001339 WO1996015462A1 (en) 1994-11-10 1995-11-10 Method of improving data obtained by a radar

Publications (2)

Publication Number Publication Date
JPH10508941A JPH10508941A (ja) 1998-09-02
JP3532924B2 true JP3532924B2 (ja) 2004-05-31

Family

ID=20395925

Family Applications (1)

Application Number Title Priority Date Filing Date
JP51598496A Expired - Fee Related JP3532924B2 (ja) 1994-11-10 1995-11-10 レーダーによって得られたデータを改善する方法

Country Status (7)

Country Link
US (1) US5861835A (ja)
EP (1) EP0805989B1 (ja)
JP (1) JP3532924B2 (ja)
CA (1) CA2204742C (ja)
DE (1) DE69520508T2 (ja)
SE (1) SE503650C2 (ja)
WO (1) WO1996015462A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2299214B (en) * 1995-03-24 1999-09-15 Commw Of Australia Signal processing
WO1998016846A1 (en) 1996-10-17 1998-04-23 Celsiustech Electronics Ab Procedure for the elimination of interference in a radar unit of the fmcw type
DE19822957C1 (de) * 1998-05-22 2000-05-25 Deutsch Zentr Luft & Raumfahrt Verfahren zur Detektion und Unterdrückung von Störsignalen in SAR-Daten und Einrichtung zur Durchführung des Verfahrens
US6094160A (en) * 1999-06-10 2000-07-25 Delco Electronics Corp. Interference rejection method for an automotive radar CW/ICC system
US6608586B1 (en) * 2002-07-25 2003-08-19 Sandia Corporation Method for removing RFI from SAR images
US7397418B1 (en) * 2006-06-05 2008-07-08 Sandia Corporation SAR image formation with azimuth interpolation after azimuth transform
CN103163505B (zh) * 2013-01-31 2014-12-03 西安电子科技大学 基于jade的时变窄带干扰抑制方法
CN106950554B (zh) * 2017-03-21 2019-04-12 中国人民解放军国防科学技术大学 实现振动干扰条件下旋转目标参数估计的方法

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4471357A (en) * 1981-10-26 1984-09-11 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Pipelined digital SAR azimuth correlator using hybrid FFT/transversal filter
US4612545A (en) * 1981-12-03 1986-09-16 Itek Corporation Automatic frequency identifier for radar signals
JPS6020162A (ja) * 1983-07-15 1985-02-01 Hitachi Ltd 合成開口レ−ダ画像再生処理方式
USH741H (en) * 1984-07-12 1990-02-06 The United States Of America As Represented By The Secretary Of The Air Force Discrete complex correlation device for obtaining subpixel accuracy
GB8700973D0 (en) * 1987-01-17 1987-07-08 Scicon Ltd Aperture synthesis
US5128679A (en) * 1989-02-21 1992-07-07 Moffat William V Electronic combat effectiveness monitor
US5061931A (en) * 1990-07-23 1991-10-29 Selenia Industrie Elettroniche Associate S.P.A. Recursive system for image forming by means of a spotlight synthetic aperture radar
US5235338A (en) * 1990-10-31 1993-08-10 Hsiao Stephen S Moving target detection through range cell migration radar
GB9103945D0 (en) * 1991-02-26 1991-06-12 Philips Electronic Associated Linearizing a swept-frequency radar
DE4122592C1 (ja) * 1991-07-08 1993-01-28 Deutsche Forschungsanstalt Fuer Luft- Und Raumfahrt Ev, 5300 Bonn, De
US5179383A (en) * 1991-07-15 1993-01-12 Raney R K Synthetic aperture radar processor to handle large squint with high phase and geometric accuracy
US5138322A (en) * 1991-08-20 1992-08-11 Matrix Engineering, Inc. Method and apparatus for radar measurement of ball in play
US5248976A (en) * 1991-11-27 1993-09-28 Hughes Aircraft Company Multiple discrete autofocus
US5410750A (en) * 1992-02-24 1995-04-25 Raytheon Company Interference suppressor for a radio receiver
US5508605A (en) * 1994-05-24 1996-04-16 Alliedsignal Inc. Method for measuring RF pulse frequency
DE4427657C2 (de) * 1994-08-05 1996-10-24 Deutsche Forsch Luft Raumfahrt Verfahren zur Bilderzeugung mittels einer zweidimensionalen Datenverarbeitung bei einem Radar mit synthetischer Apertur
US5592170A (en) * 1995-04-11 1997-01-07 Jaycor Radar system and method for detecting and discriminating targets from a safe distance

Also Published As

Publication number Publication date
DE69520508T2 (de) 2001-10-11
SE9403872L (sv) 1996-05-11
JPH10508941A (ja) 1998-09-02
EP0805989B1 (en) 2001-03-28
WO1996015462A1 (en) 1996-05-23
DE69520508D1 (de) 2001-05-03
CA2204742C (en) 2007-01-09
US5861835A (en) 1999-01-19
CA2204742A1 (en) 1996-05-23
EP0805989A1 (en) 1997-11-12
SE503650C2 (sv) 1996-07-22

Similar Documents

Publication Publication Date Title
EP2097769B1 (en) System and method for reducing the effect of a radar interference signal
EP1672379B1 (en) System and method for reducing a radar interference signal
US7885688B2 (en) Methods and systems for signal selection
Zaimbashi et al. Invariant target detection in multiband FM-based passive bistatic radar
TWI358913B (en) Active receiver detection and ranging
JP3532924B2 (ja) レーダーによって得られたデータを改善する方法
CN109061580A (zh) 一种调频间断连续波雷达缓和相互干扰方法
EP1521097A2 (en) Pulse compression processor
JP5699404B2 (ja) レーダ受信信号処理装置とその方法
DE2753764A1 (de) Verfahren und vorrichtung zur funkpeilung
CN115629360B (zh) 一种线性调频的短脉冲采样调制转发干扰信号生成方法
JPH0743453A (ja) 短波レーダ装置
CN115825953A (zh) 一种基于随机频率编码信号的前视超分辨成像方法
RU2296345C2 (ru) Способ разрешения целей по дальности радиолокационной станцией и импульсная радиолокационная станция со сжатием импульсов и восстановлением сигналов
EP4040181A1 (en) Radar signal processing device, radar device, radar signal processing method and radar signal processing program
Norman Backscatter ionogram inversion
CN112881984B (zh) 一种雷达信号抗干扰处理方法、装置及存储介质
CN111665475B (zh) 雷达主瓣干扰抑制方法及装置、电子设备、存储介质
RU2216748C2 (ru) Способ распознавания сигналов систем радиосвязи
CN112098984B (zh) 散射信号的目标体检测方法及装置
EP4040182A1 (en) Radar signal processing device, radar device, radar signal processing method and non-transitory computer-readable medium
JPH10206524A (ja) 目標検出装置
JP3886797B2 (ja) トーン検出器
Shinriki et al. Multi-range resolution radar using sideband spectrum energy
CN104535975A (zh) 一种用于合成孔径雷达回波模拟器的图像校准装置及方法

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20031127

TRDD Decision of grant or rejection written
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20040107

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20040203

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040305

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090312

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100312

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110312

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130312

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140312

Year of fee payment: 10

LAPS Cancellation because of no payment of annual fees