(57)【特許請求の範囲】
【請求項１】
実開口アンテナの受信信号から取得される目標をレンジ及びドップラ周波数の両軸で追跡し、前記目標の画像中心を得てレンジ圧縮及びＡＺ圧縮することにより、レンジ−ドップラ周波数軸で目標を画像化するＩＳＡＲ（Inverse Synthetic Aperture Radar：逆合成開口レーダ）方式のレーダ装置において、
前記実開口アンテナの受信信号からモノパルス出力であるΣ、ΔＡＺ及びΔＥＬそれぞれのビームを生成してサンプリング出力するビーム生成手段と、
前記Σ、ΔＡＺ及びΔＥＬそれぞれのビームのサンプリング出力からレンジ−ドップラ周波数軸におけるレンジ−ドップラ画像を生成する画像生成手段と、
前記Σビームのレンジ−ドップラ画像の中で所定の目標画像を抽出する目標画像抽出手段と、
前記抽出された目標画像について、前記Σビームによるレンジと前記Σ、ΔＡＺ及びΔＥＬそれぞれのビームによるモノパルス測角値とに基づいて前記目標画像のドップラ−レンジセルの位置を算出して目標画像を得る目標画像取得手段と
を具備するレーダ装置。
【請求項２】
前記目標画像取得手段は、前記実開口アンテナの合成開口長におけるＣＰＩ（Coherent Processing Interval：コヒーレント処理間隔）の任意の複数ヒットを２分割し、前記Σビームによる和信号ΣSARによりドップラ−レンジセルを選定し、選定したセルに対応するΣSARと前記Σビームによる差信号ΔSARを用いたモノパルス処理により円錐角θSARを測角し、測角した円錐角θSARの中で最も前記ΔＡＺ及びΔＥＬのビームそれぞれの測角値θAZ、θELに近い角度θAZcalとθELcalにより、前記ドップラ−レンジセルの位置を算出する請求項１記載のレーダ装置。
【請求項３】
前記ビーム生成手段は、前記実開口アンテナの合成開口長におけるＣＰＩ（Coherent Processing Interval：コヒーレント処理間隔）の任意の複数ヒット内で、疑似ランダム順に前記Σビーム、ΔＡＺビーム及びΔＥＬビームを切り替えてサンプリング出力する請求項１または２に記載のレーダ装置。
【請求項４】
前記目標画像取得手段は、前記実開口アンテナの合成開口長におけるＣＰＩ（Coherent Processing Interval：コヒーレント処理間隔）の任意の複数ヒットを２分割し、前記Σ、ΔＡＺ、ΔＥＬの各ビームの出力から第１のモノパルス測角曲線を算出し、前記Σビームによる和信号ΣSARと前記Σビームによる差信号ΔSARから第２のモノパルス測角曲線を算出し、所定のレンジ範囲毎にテーブル化して、観測した目標レンジにより前記第１、第２のモノパルス測角曲線を切り替えて、切り替えられたモノパルス測角曲線に基づいて前記目標画像のドップラ−レンジセルの位置を算出する請求項１記載のレーダ装置。
【請求項５】
実開口アンテナの受信信号から取得される目標をレンジ及びドップラ周波数の両軸で追跡し、前記目標の画像中心を得てレンジ圧縮及びＡＺ圧縮することにより、レンジ−ドップラ周波数軸で目標を画像化するＩＳＡＲ（Inverse Synthetic Aperture Radar：逆合成開口レーダ）方式のレーダ信号処理方法において、
前記実開口アンテナの受信信号からモノパルス出力であるΣ、ΔＡＺ及びΔＥＬそれぞれのビームを生成してサンプリング出力し、
前記Σ、ΔＡＺ及びΔＥＬそれぞれのビームのサンプリング出力からレンジ−ドップラ周波数軸におけるレンジ−ドップラ画像を生成し、
前記Σビームのレンジ−ドップラ画像の中で所定の目標画像を抽出し、
前記抽出された目標画像について、前記Σビームによるレンジと前記Σ、ΔＡＺ及びΔＥＬそれぞれのビームによるモノパルス測角値とに基づいて前記目標画像のドップラ−レンジセルの位置を算出して目標画像を得るレーダ信号処理方法。
【請求項６】
前記目標画像の取得は、前記実開口アンテナの合成開口長におけるＣＰＩ（Coherent Processing Interval：コヒーレント処理間隔）の任意の複数ヒットを２分割し、前記Σビームによる和信号ΣSARによりドップラ−レンジセルを選定し、選定したセルに対応するΣSARと前記Σビームによる差信号ΔSARを用いたモノパルス処理により円錐角θSARを測角し、測角した円錐角θSARの中で最も前記ΔＡＺ及びΔＥＬのビームそれぞれの測角値θAZ、θELに近い角度θAZcalとθELcalにより、前記ドップラ−レンジセルの位置を算出する請求項５記載のレーダ信号処理方法。
【請求項７】
前記ビーム生成は、前記実開口アンテナの合成開口長におけるＣＰＩ（Coherent Processing Interval：コヒーレント処理間隔）の任意の複数ヒット内で、疑似ランダム順に前記Σビーム、ΔＡＺビーム及びΔＥＬビームを切り替えてサンプリング出力する請求項５に記載のレーダ信号処理方法。
【請求項８】
前記目標画像の取得は、前記実開口アンテナの合成開口長におけるＣＰＩ（Coherent Processing Interval：コヒーレント処理間隔）の任意の複数ヒットを２分割し、前記Σ、ΔＡＺ、ΔＥＬの各ビームの出力から第１のモノパルス測角曲線を算出し、前記Σビームによる和信号ΣSARと前記Σビームによる差信号ΔSARから第２のモノパルス測角曲線を算出し、所定のレンジ範囲毎にテーブル化して、観測した目標レンジにより前記第１、第２のモノパルス測角曲線を切り替えて、切り替えられたモノパルス測角曲線に基づいて前記目標画像のドップラ−レンジセルの位置を算出する請求項５記載のレーダ信号処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本実施形態は、レーダ装置及びレーダ信号処理方法に関する。
【背景技術】
【０００２】
レーダ装置にあっては、画像により目標を識別する方法として、目標の重心等をレンジ及びドップラの両軸で追跡して画像中心を取得し、レンジ圧縮及びＡＺ圧縮することにより、レンジ−ドップラ軸で目標を画像化するＩＳＡＲ（Inverse Synthetic Aperture Radar：逆合成開口レーダ）処理が知られている（非特許文献１参照）。従来のＩＳＡＲ処理では、剛体目標が移動している場合、その目標の回転運動により、目標上の各反射点におけるドップラ速度が異なる。このため、ＩＳＡＲ処理によって生成した画像が真の目標形状と異なり、形状寸法の絶対値とも異なるのが通常である。
【０００３】
したがって、ＩＳＡＲ画像により目標を識別する方法では、目標の回転や振動によるドップラ速度の変化により、生成画像の形状や寸法が真の目標とは大きく異なってしまい、識別が困難になる場合があった。
【０００４】
また、ＩＳＡＲ画像を基に、大きさを持つ目標内の特定の点を検出して追跡する場合には、クロスレンジが絶対値と異なるため、誤差が大きく追跡できなくなるという問題があった。
【０００５】
また、ＩＳＡＲ画像を生成しないドップラレーダ等において、大きさを持つ目標を測角する場合には、測角点が広がった複数点のベクトル合成（複素合成）になるため、グリントノイズが大きくなって測角の精度も低下してしまうという問題があった。
【０００６】
また、目標が被搭載機のフライト軸方向にある場合には、剛体目標の回転運動が小さいと、クロスレンジ方向の左右上下等の目標の反射点のドップラが同じになるため、真の目標形状とＳＡＲ画像の形状が大きく異なってしまい、識別できないという問題があった。
【０００７】
また、合成開口長が短い場合には、剛体目標の回転運動が小さくても、画像セルの絶対位置精度が低いという問題があった。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】ＳＡＲ方式（ＩＳＡＲ）、吉田、‘改定レーダ技術’、電子情報通信学会、pp.280-283(1996)
【非特許文献２】ＳＡＲ方式（レンジ圧縮）、大内、“リモートセンシングのための合成開口レーダの基礎”、東京電機大学出版局、pp.131-149(2003)
【非特許文献３】ＳＡＲ方式（ＡＺ圧縮）、大内、“リモートセンシングのための合成開口レーダの基礎”、東京電機大学出版局、pp.171-178(2003)
【非特許文献４】ＳＡＲ処理方式（ポーラフォーマット変換再構成処理）、MEHRDAD SOUMEKH,“Synthetic Aperture Radar Signal Processing”, JOHN WILEY & SONS,INC.,pp.319-325(1999)
【非特許文献５】位相モノパルス（位相比較モノパルス）方式、電子情報通信学会、改訂レーダ技術、pp.262-264（1996）
【非特許文献６】振幅モノパルス（振幅比較モノパルス）方式、電子情報通信学会、改訂レーダ技術、pp.260-262（1996）
【非特許文献７】ファーフィ−ルドレンジ、電子情報通信学会、改訂レーダ技術、pp.144-145（1996）
【特許文献】
【０００９】
【特許文献１】特開２００８−３０４３２１号公報
【発明の概要】
【発明が解決しようとする課題】
【００１０】
以上述べたように、従来のレーダ装置に適用されるＩＳＡＲ処理では、レーダ観測位置と目標の複数反射点の各点との相対速度が異なる場合に、生成画像の形状や寸法が真の目標形状や寸法と大きく異なってしまうという課題があった。この課題に伴い、目標の識別能力や目標内の特定点の位置観測能力が低下する、ドップラレーダ等にあってはグリントノイズが大きくなって測角の精度も低下してしまう、目標の位置を観測可能な被搭載機のフライト軸方向に目標がある場合には、真の目標形状とＳＡＲ画像の形状が大きく異なって識別できない、合成開口長が短い場合には、剛体目標の回転運動が小さくても、画像セルの絶対位置精度が低いという課題があった。
【００１１】
本実施形態は上記課題に鑑みなされたもので、レーダと目標の複数反射点の各点の相対速度が異なる場合でも、真の目標形状や寸法に近い画像を生成することができるレーダ装置とそのレーダ信号処理方法を提供することを目的とする。さらには、識別能力向上、目標内の特定点の位置観測精度の向上、ドップラレーダ等におけるグリントノイズを低減して高精度に測角し目標の位置を観測することができる自機のフライト軸方向に目標がある場合でも目標形状に近いＳＡＲ画像を生成し、また、剛体目標の回転運動が小さく、合成開口長が小さい場合でも絶対位置精度の高いＳＡＲ画像を得ることのできるレーダ装置とレーダ信号処理方法を提供することを目的とする。
【課題を解決するための手段】
【００１２】
上記の課題を解決するために、本実施形態は、実開口アンテナの受信信号から取得される目標をレンジ及びドップラ周波数の両軸で追跡し、前記目標の画像中心を得てレンジ圧縮及びＡＺ圧縮することにより、レンジ−ドップラ周波数軸で目標を画像化するＩＳＡＲ（Inverse Synthetic Aperture Radar：逆合成開口レーダ）方式のレーダ装置において、ビーム生成手段により、前記実開口アンテナの受信信号からモノパルス出力であるΣ、ΔＡＺ及びΔＥＬそれぞれのビームを生成してサンプリング出力し、画像生成手段により、前記Σ、ΔＡＺ及びΔＥＬそれぞれのビームのサンプリング出力からレンジ−ドップラ周波数軸におけるレンジ−ドップラ画像を生成し、目標画像抽出手段により前記Σビームのレンジ−ドップラ画像の中で所定の目標画像を抽出し、前記目標画像取得手段により、前記抽出された目標画像について、前記Σビームによるレンジと前記Σ、ΔＡＺ及びΔＥＬそれぞれのビームによるモノパルス測角値とに基づいて前記目標画像のドップラ−レンジセルの位置を算出して目標画像を得る。
【図面の簡単な説明】
【００１３】
【図１】第１の実施形態に係るレーダ装置の構成を示すブロック図。
【図２】図１に示すレーダ装置において、狭帯域処理器の構成を示すブロック図。
【図３】図１に示すレーダ装置において、合成開口におけるＩＳＡＲ処理の概要を説明するための概念図。
【図４】図１に示すレーダ装置において、合成開口による取得データからＰＲＩ毎のデータ列を抽出する様子を示す図。
【図５】図１に示すレーダ装置において、参照信号として生成される近似曲線を示す図。
【図６】図１に示すレーダ装置において、ＳＡＲ画像を取得する手順を示す図。
【図７】図１に示すレーダ装置において、ポーラフォーマット変換により格子点のデータを生成する様子を示す図。
【図８】図１に示すレーダ装置において、レンジ・測角値に基づく座標変換による画像変換の様子を示す図。
【図９】第２の実施形態に係るレーダ装置の構成を示すブロック図。
【図１０】図９に示すレーダ装置において、合成開口におけるＩＳＡＲ処理の概要を説明するための概念図。
【図１１】図９に示すレーダ装置において、ビームの円錐の曲線特性を示す概念図。
【図１２】第３の実施形態に係るレーダ装置の構成を示すブロック図。
【図１３】図１２に示すレーダ装置におけるサンプル手法を説明するための図。
【図１４】第４の実施形態に係るレーダ装置の構成を示すブロック図。
【図１５】（ａ）はレンジ１〜Ｎのビーム幅の広がり具合を示し、（ｂ）はレンジ１〜Ｎの誤差電圧と目標レンジとの関係を示す図。
【発明を実施するための形態】
【００１４】
以下、実施形態について、図面を参照して説明する。尚、各実施形態の説明において、同一部分には同一符号を付して示し、重複する説明を省略する。
【００１５】
（第１の実施形態）
以下、図１乃至図８を参照して、第１の実施形態に係るレーダ装置について説明する。
【００１６】
図１は第１の実施形態に係るレーダ装置の構成を示すブロック図である。図１に示すレーダ装置において、アンテナ１は複数のアンテナ素子を配列して大開口アレイを形成してなるフェーズドアレイアンテナであり、送受信器２の送受信部２１から特定の周期で繰り返し供給される特定周波数の送信パルス信号（以下、ＰＲＦ（Pulse Repetition Frequency）信号）を指定方向に送出してその反射波を受信する。送受信器２は、送受信部２１において、アンテナ１の複数のアンテナ素子でそれぞれ受信された信号をビーム制御部２２からの指示に従って位相制御を施し合成することで、任意の方向に受信ビームを形成してＰＲＦ受信信号を取得する。ここで、ビーム制御部２２は指定された目標方向の測角値に基づいてΣビーム、ΔＡＺビーム、ΔＥＬビームを形成するように、送受信部２１に対して各ビームに対応する位相制御を施す。これにより、送受信部２１はΣ信号、ΔＡＺ信号、ΔＥＬ信号を生成して信号処理器３へ出力する。
【００１７】
上記信号処理器３に入力されたΣ信号、ΔＡＺ信号、ΔＥＬ信号は、それぞれの系統において、ＡＤ（Analog-Digital）変換部３１１，３１２，３１３で系統別にデジタル信号に変換されて狭帯域処理部３２に送られる。狭帯域処理部３２は、入力したΣ信号、ΔＡＺ信号、ΔＥＬ信号から目標方向を測角演算する。ここで得られた測角値（Σ、ΔＡＺ、ΔＥＬ）はビーム制御部２２に送られる。
【００１８】
図２に上記狭帯域処理部３２の具体的な構成を示す。図２において、デジタル化されたΣ信号、ΔＡＺ信号、ΔＥＬ信号は、それぞれの系統において、データ抽出部３２１１〜３２１３に送られる。データ抽出部３２１１〜３２１３は、サンプリングレートを下げて処理規模を削減するために、入力されたΣ信号、ΔＡＺ信号、ΔＥＬ信号を所定の周波数フィルタに通した後、ＰＲＩ（Pulse repetition Interval）内のデータをレンジセル単位で取得する。各データ抽出部３２１１〜３２１３で得られたΣ、ΔＡＺ、ΔＥＬそれぞれのレンジセル信号は、ＰＲＩ軸ＦＦＴ（Fast Fourier Transformation）処理部３２２１〜３２２３によってＰＲＩ軸の周波数領域信号に変換された後、レンジ（時間）圧縮部３２３１〜３２３３によってレンジ軸上でパルス圧縮されてレンジ圧縮信号となる。
【００１９】
このうち、Σ系統のレンジ圧縮信号については、ＣＦＡＲ（Constant False Alarm Rate）処理部３２４に送られ、そのＣＦＡＲ処理によって複数の極値を持つレンジセルが検出される。ここで検出されたレンジセルに対し、スレショルド検出部３２５は、閾値以上の極値を持つレンジセルを検出してそのレンジ圧縮信号を出力する。また、ΔＡＺ系統、ΔＥＬ系統のレンジ圧縮信号については、それぞれセル検出部３２６１，３２６２に送られ、それぞれΣ系統のスレショルド検出部３２５で検出されたレンジセルに対応するレンジセルの信号が抽出されて、Σ系統の検出セルの信号と共に測角部３２７に送られ、ここで目標方向の測角演算が行われる。
【００２０】
一方、上記ＡＤ変換部３１１でデジタル化されたΣ系統の受信信号は、データ抽出部３３に送られる。このデータ抽出部３３では、目標重心等の代表点のドップラ−レンジを観測してビームを目標方向に指向させ、ドップラ−レンジの情報を基にして、ＡＺ圧縮用の際の参照信号を生成するためのデータを抽出する。ここで抽出されたデータは、ＰＲＩ軸ＦＦＴ処理部３４にてＣＰＩ（Coherent Processing Interval）毎にＰＲＩ軸の周波数領域信号に変換され、レンジ圧縮部３５にてレンジ軸上でパルス圧縮される。このパルス圧縮信号は、ＣＦＡＲ処理部３６に送られ、そのＣＦＡＲ処理によって複数の極値を持つレンジセルが検出される。ここで検出されたレンジセルに対し、スレショルド検出部３７は、閾値以上の極値を持つレンジセルが検出する。レンジ・ドップラ算出部３８は、スレショルド検出部３７で検出されたレンジセルの信号からＣＰＩ毎に目標のレンジとドップラ周波数を算出する。参照信号生成部３９は、ＣＰＩ毎に検出したレンジとドップラ周波数からパルス圧縮のための参照信号を生成する。
【００２１】
また、上記ＡＤ変換部３１１〜３１３でデジタル化されたΣ、ΔＡＺ、ΔＥＬそれぞれの受信信号は、データ抽出部３Ａ１，３Ａ２，３Ａ３にてサンプリング周波数を満足するための所定の周波数フィルタに通される。そして、レンジ圧縮部３Ｂ１，３Ｂ２，３Ｂ３によってレンジ圧縮され、さらにＡＺ圧縮部３Ｃ１，３Ｃ２，３Ｃ３によって参照信号に基づいてクロスレンジ圧縮（ＡＺ圧縮）された後、ポーラフォーマット変換による画像化部３Ｄ１，３Ｄ２，３Ｄ３に送られる。画像化部３Ｄ１，３Ｄ２，３Ｄ３は、レンジ・クロスレンジ圧縮された受信信号に対して２次元ＦＦＴ処理により周波数領域の信号に変換した後、ポーラフォーマット変換による画像処理を行い、２次元逆ＦＦＴ処理により時間領域に信号に戻して画像化する。尚、画像化部３Ｄ１，３Ｄ２，３Ｄ３は、画像品質を向上する場合に用いればよく、省略可能である。
【００２２】
その後、Σ系統では、セル選定部３Ｅにおいて、画像化部３Ｄ１，３Ｄ２，３Ｄ３で生成された画像の中で振幅が所定のスレショルドを超えるセルを選定する。また、ΔＡＺ，ΔＥＬ系統では、セル選定部３Ｅで選定されたセルをセル抽出部３Ｆ１，３Ｆ２に通知し、選定セルに対応するΔＡＺ，ΔＥＬのセルを抽出する。測角部３Ｇは、セル選定部３Ｅ及びセル抽出部３Ｆ１，３Ｆ２で得られたΣ、ΔＡＺ、ΔＥＬのセルの信号から目標方向を測角し、画像化部３Ｈは測角部３Ｇで測角された目標を画像化し出力する。
【００２３】
上記構成によるレーダ装置において、以下、図３乃至図８を参照してそのＩＳＡＲ処理を説明する。
【００２４】
まず、図３を参照してレーダ装置の合成開口の概要を説明する。飛翔体搭載のレーダ装置では、実開口ビームを目標に常に照射するように向けて、合成開口時間（１サイクル）内にＰＲＩ（Pulse Repetition Interval）間隔でパルスを送信し、その反射波を受信してパルス毎にＰＲＩ内のデータをレンジセル単位で取得する。この取得データを用いてＩＳＡＲ処理を実施して目標のＩＳＡＲ画像を得る。尚、図３は飛翔体搭載のレーダ装置の場合を示しているが、ＩＳＡＲ画像を得られれば、レーダ装置は固定の場合でもよい。
【００２５】
図１において、ビーム制御部２２により画像化したい範囲にビームを指向させて送受信し、その受信信号をΣ、ΔＡＺ、ΔＥＬそれぞれの系統振り分けてデジタル信号に変換する。ビーム指向方向は、狭帯域処理部３２により、目標方向を測角した方向とする。この測角については、狭帯域処理を用いなくても、広帯域のモノパルス測角の結果を用いてもよい。
【００２６】
狭帯域処理部３２では、図２に示すように、図１の系統のΣ、ΔＡＺ及びΔＥＬのＡＤ変換出力を入力とし、データ抽出部３２１１〜３２１３により、サンプリングレートを下げて処理規模を削減するために、周波数フィルタで所定の帯域を抽出した後、サンプリングする。そのうち、Σ系統の信号をＰＲＩ軸ＦＦＴ処理部３４でＰＲＩ軸の周波数領域信号に変換して、レンジ圧縮部３５でレンジ方向にパルス圧縮し、ＣＦＡＲ処理部３６で複数の極値を持つレンジセルを検出し、スレショルド検出部３７で閾値以上の極値を持つレンジセルを検出する。この検出したセルに対応するΔＡＺとΔＥＬのセルをセル抽出部３２６１，３２６２で抽出し、測角部３２７でΣ、ΔＡＺ及びΔＥＬ信号を用いて測角する。この測角値を用いて、図１のビーム制御部２２によりビームを目標方向に指向させる。
【００２７】
次に、図１のΣ系のデータ抽出部３３により、必要なサンプルデータを抽出する。データ抽出３３では、目標重心等の代表点のドップラ−レンジを観測し、ビームを目標方向に指向させたときのドップラ−レンジ情報をもとにして、ＡＺ圧縮用の参照信号を生成するためのデータを抽出する。すなわち、図４に示すように、合成開口のために取得したデータＮallの中から、重複を含めてＮヒットのＰＲＩデータをPRI1〜PRINまで抽出する。ここで、CPI(m)（ｍ＝1〜Ｍ）とすると、ＣＰＩ毎にＰＲＩ軸ＦＦＴ処理部３４でＰＲＩ軸の周波数領域信号に変換し、レンジ圧縮部３５でレンジ方向にパルス圧縮し、ＣＦＡＲ処理部３６で複数の極値を持つセルを検出し、スレショルド検出部３７でスレショルド以上の極値を持つセルを検出し、レンジ・ドップラ算出部３８により、目標のレンジ及びドップラ周波数を算出する。次に、参照信号生成部３９により、ＣＰＩ毎に算出したレンジとドップラ周波数により、図５に示す参照信号Refの近似曲線を算出する。これを定式化すると次の通りである。
【数１】
0006363524.tif 000002
【００２８】
一方、アンテナ１のモノパルス出力ΣとΔＡＺとΔＥＬの出力は、それぞれＡＤ変換部３１１，３１２，３１３によりデジタル信号に変換され、データ抽出部３Ａ１，３Ａ２，３Ａ３によりサンプリング出力され、レンジ圧縮部３Ｂ１，３Ｂ２，３Ｂ３により変調信号の波形に対応するレンジ圧縮用の参照信号と相関演算された後、ＡＺ圧縮部３Ｃ１，３Ｃ２，３Ｃ３において（１）式で生成された参照信号との相関演算によりＡＺ圧縮される。これを定式化すると次の通りである。
【００２９】
まずレンジ圧縮について述べる（非特許文献２参照）。レンジ圧縮は、入力信号とレンジ圧縮用信号の相関処理であり、これを周波数領域で行う場合について定式化すると次の通りである。
【数２】
0006363524.tif 000003
【００３０】
また、参照信号sref（線形チャープ信号）は次式で表現できる。
【数３】
0006363524.tif 000004
【００３１】
この参照信号sref(t)のサンプル長を入力信号に合わせて０埋めした信号に置き換える。
【数４】
0006363524.tif 000005
【００３２】
時間軸上の信号にするには、このsを逆フーリエ変換すればよいが、この後クロスレンジ圧縮（ＡＺ圧縮、非特許文献３参照）を行うために、信号sの（ω，ｕ）軸のままとする。次にクロスレンジ圧縮を行うが、その参照信号は（１）式を用い、ω軸については同じ値を入れたものとする。
【数５】
0006363524.tif 000006
【００３３】
（５）式と（６−３）式を乗算して信号csを得る。
【数６】
0006363524.tif 000007
【００３４】
これを用いて、ｕ軸でＦＦＴ処理して周波数領域の信号fcs（ω，ku）を得る。
【数７】
0006363524.tif 000008
【００３５】
ＦＦＴ画像出力は、信号fcsのω軸に関する逆ＦＦＴにより算出することができる。
【数８】
0006363524.tif 000009
【００３６】
次に、より精度の高い画像生成手法として、図６及び図７を参照してポーラフォーマット変換による画像化部３Ｄ１，３Ｄ２，３Ｄ３のＳＡＲ画像処理について述べる（非特許文献４参照）。
【００３７】
図６は第１の実施形態に係るレーダ装置でＳＡＲ画像を取得する手順を示している。まず、サンプリング時間ｔに対する合成開口長のサンプル点のデータｕを入力信号sig(t,u)として取り込んで（Ｓ１１）、ｔ軸に対するフーリエ変換（FFTx(ω,u))を行い（Ｓ１２）、その結果Sin(ω,u)に参照信号のｔ軸に対するフーリエ変換（FFTx(ω,u))の結果を乗算してレンジ圧縮を行い（Ｓ１３）、レンジ圧縮の出力ｓを得る。このレンジ圧縮信号ｓを用いて、ｕ軸に関してフーリエ変換する（Ｓ１４）。
【数９】
0006363524.tif 000010
【００３８】
次にＡＺ圧縮用の参照信号を生成する。
【数１０】
0006363524.tif 000011
【００３９】
また、kx,kyは次式により求めることができる。
【数１１】
0006363524.tif 000012
【００４０】
画像中心については、レンジ、ドップラ、角度を観測して、次式により算出する位置とする。
【数１２】
0006363524.tif 000013
【００４１】
ここで、観測距離、観測ＡＺ角度、観測ＥＬ角度については、例えば、図４の各ＣＰＩにおける観測値の平均値とすればよい。
【００４２】
このＡＺ圧縮参照信号を用いて、次式によりfsmを算出する（Ｓ１５）。
【数１３】
0006363524.tif 000014
【００４３】
fsmを用いて、図７に示すポーラフォーマット変換を行い、kx軸、ky軸で格子点のデータF(kx,ky)を生成する（Ｓ１６）。ポーラフォーマット変換は、取得データを用いて、(kx,ky)軸の格子点の補間後データを内挿手法等を用いて算出する手法であり、細部については非特許文献４の通りであり割愛する。このfsmを用いて、２次元逆ＦＦＴによりポーラフォーマット変換を用いた画像fsarを出力する（Ｓ１７）。
【数１４】
0006363524.tif 000015
【００４４】
このポーラフォーマット変換による画像処理は、常に必要ではなく、画像品質を向上する場合に用いればよい。
【００４５】
次に、Σ系の画像出力を用いて、Ｘ−Ｙ軸（ポーラフォーマット変換を用いない場合はドップラ−レンジ軸に対応）の画像の中で、振幅が所定のスレショルドを超えるセル（Ｐ個）を選定する。このＰ個のセルの各々に対応するΔＡＺ及びΔＥＬ系の画像出力のセルをセル抽出部３Ｆ１，３Ｆ２で抽出し、測角部３Ｇにおいて、Σ、ΔＡＺ及びΔＥＬのセル信号（複素信号）を用いてモノパルス測角する（非特許文献５参照）。この測角値と距離を用いて、画像化部３Ｈにおいて、次式により画像化する。
【数１５】
0006363524.tif 000016
【００４６】
この画像変換の様子を図８に示す。図８において、（ａ）はレンジ−ドップラ画像であり、（ｂ）は（ａ）の画像のレンジ・測角値による座標変換を行った変換画像である。図８から明らかなように、レンジ・測角値による座標変換により、真値の位置に近い画像を生成することができる。
【００４７】
尚、第１の実施形態の手法によれば、被搭載機のフライト方向に目標があり、目標が回転運動を持たず、クロスレンジ方向にアンビギュイティを持ち、ＳＡＲ画像が実際の目標形状と大きく異なる場合でも、実開口の測角はクロスレンジ方向の区別ができるため、実際の形状に近い画像を得ることができる。
【００４８】
（第２の実施形態）
第１の実施形態では、実開口の測角精度が低い場合がある。第２の実施形態は実開口の測角精度を改善する手法について述べる。
【００４９】
図９は第２の実施形態に係るレーダ装置の構成を示すブロック図である。図９において、図１に示した構成と異なる点は、ΔＳＡＲウェイト乗算部３Ｉ、ＡＺ圧縮部３Ｃ４、ポーラフォーマット変換による画像化部３Ｄ４、セル抽出部３Ｆ３を備え、図１に示した測角部３Ｇに代わって、ＳＡＲ測角部３Ｊ、実開口測角部３Ｋ、測角相関部３Ｌを備えるようにしたことにある。
【００５０】
第１の実施形態では、実開口アンテナのモノパルス測角値を用いたが、測角精度を向上するために、図１０に示すような合成開口を用いる。すなわち、目標が飛翔経路軸付近にある場合は、合成開口によるリニアアレイのエンドファイア方向（リニアアレイの軸方向）にビームが指向される。この場合、ブロ−ドサイド方向（リニアアレイの軸に垂直方向）にビームを形成する場合に比べて合成開口ビーム幅は広がるが、合成開口長が大きい場合には、実開口アンテナによるビーム幅よりも十分狭くなる。このエンドファイアアレイのＰＲＩ毎のＮサンプルのデータを使った和信号であるΣ信号と、１〜Ｎ／２とＮ／２＋１〜Ｎの開口２分割した差信号であるΔ信号を用いてモノパルス測角（非特許文献５参照）を行う。これによってＳＡＲ測角値が得られる（特許文献１参照）。
【００５１】
このビ−ムは、エンドファイアアレイの軸を中心に図１０及び図１１に示すような円錐状のビ−ムであり、この円錐状に目標がある場合には、一意に位置を特定できない。この対策のために、第１の実施形態と同じように実開口アンテナのモノパルス測角値（ＡＺ角とＥＬ角）を用いる。具体的には、目標距離ＲとエンドファイアによるＳＡＲ測角値θSARとにより、図１０に示すビームの円錐の曲線を描き、(16)式で決まる実開口アンテナによる目標位置から最も近い円錐状の点を抽出し、第１の実施形態の変換位置と置き換えればよい。
【００５２】
以上の処理を図９に示す系統図で説明する。第１の実施形態と同様の部分の説明は省略する。レンジ圧縮部３Ｂ１でΣ信号のレンジ圧縮した出力のうち、上記のように開口２分割した信号に対してΔ信号用のサイドローブを低減するためのウェイトをΔＳＡＲウェイト乗算部３Ｉにより乗算し、ＡＺ圧縮部３Ｃ４、ポーラフォーマット変換による画像化部３Ｄ１でそれぞれＡＺ圧縮、ポーラフォーマット変換を行って画像化する。さらに、セル選定部３Ｅにて選定されたセルをセル抽出部３Ｆ３に通知し、Σの選定セルに対応するΔSARの信号を抽出する。これにより得られたΣ及びΔSARの信号を用いて、次式の誤差電圧εを演算する。
【数１６】
0006363524.tif 000017
【００５３】
測角部３Ｇにおいて、上記誤差電圧εと予め作成した誤差電圧テーブルとの比較により測角すれば、円錐状の角度を算出することができる。このＳＡＲ測角値を用いて、前述の処理により、高精度な目標位置を算出すればよい。
【００５４】
なお、ＳＡＲ測角について、ΣとΔsar信号による位相モノパルス測角の場合について述べたが、Δsar信号の代わりに、Σ信号ビームからスクイント角分ずらせたΣ２sarビームを用いて、スクイント測角してもよいのは言うまでもない。この場合は、誤差電圧εは、次式となる。
【数１７】
0006363524.tif 000018
【００５５】
この誤差電圧εと予め作成した誤差電圧テーブルを用いて測角することができる。また、振幅のみを用いた振幅モノパルスの場合でも測角できる（非特許文献６参照）。
【数１８】
0006363524.tif 000019
【００５６】
この測角値と距離を用いて、画像化部３Ｈにおいて、（１６）式により３次元位置を画像化する。これにより、第２の実施形態によれば、目標が飛翔経路軸付近にある場合でも、実開口の測角精度を改善して真値に近い画像を取得することができる。
【００５７】
なお、第２の実施形態においては、合成開口長を開口２分割したモノパルス測角手法の主旨で説明しているが、合成開口長内のサンプル列のＦＦＴを用いる手法であり、ドップラ周波数軸におけるモノパルス手法とも言える。これより、剛体目標の回転運動が大きい場合には、そのドップラ周波数成分による位相誤差が生じるため、クロスレンジ方向に画像が歪むのと同様に、測角値も誤差を持つことになる。したがって、本手法により絶対位置精度の向上を図る条件としては剛体目標の回転運動が低い場合か、剛体目標の回転運動を観測し、そのドップラによる位相成分を補正できる場合等となり、モノパルス測角手法を用いる分だけ、クロスレンジの精度向上となる。
【００５８】
（第３の実施形態）
第１の実施形態及び第２の実施形態では、Σの他にΔＡＺとΔＥＬ系統が必要であり、処理規模が増える。第３の実施形態は、Σ、ΔＡＺ、ΔＥＬの各系統の共用を図り、処理規模の増加を抑えるものである。
【００５９】
図１２は第３の実施形態に係るレーダ装置の構成を示すブロック図、図１３は第３の実施形態におけるサンプル手法を説明するための図である。第１及び第２の実施形態と同様の部分の説明は省略する。
【００６０】
第１及び第２の実施形態では、Σ、ΔＡＺ及びΔＥＬの３系統分を同時にサンプルしていたが、第３の実施形態では、送受信器２において、送受信部２１で得られるΣ、ΔＡＺ、ΔＥＬ信号をサンプル切り替え部２３で各信号をサンプリングし、ＡＤ変換を１系統分としてＡＤ変換部３１のみとしている。
【００６１】
すなわち、送受信器２のサンプル切り替え部２３は、図１３（ａ），（ｂ）に示すように、各系統を時分割に切り替える。この場合、３系統分を順にまとめて同時にサンプルする方式とすると、時間差が生じるため、角度誤差が大きくなる。この対策のために、３系統分のサンプル順を擬似ランダム化して、間引きサンプルを実施する。ＡＤ変換部を１ｃｈとするためには、送受信器２において、サンプル切り替え部２３により、擬似ランダム順に送受信部２１からの３系統の順番を切り替える。以降の処理は第１及び第２の実施形態と同様である。
【００６２】
（第４の実施形態）
レーダと目標の距離は遠距離の場合と近距離の場合があるが、第１及び第２の実施形態では、特に近距離の場合に、ファーフィールドレンジ（非特許文献７参照）の条件を満足せず、ビーム幅が広がる場合が想定される。この場合に、遠距離の場合と同じ測角テーブルを用いると、誤差電圧との対応がとれずに、測角値に誤差が生じる。第４の実施形態ではこの対策を施している。図１４は、第４の実施形態に係るレーダ装置の構成を示すブロック図である。
【００６３】
測角曲線としては、一般にファーフィールド条件を満足する遠距離（レンジ１）はビーム幅が狭く、ファーフィールド条件を満足しない近距離（レンジＮに近い方）はビーム幅が広くなる。この様子を図１５に示す。図１５（ａ）はレンジ１〜Ｎのビーム幅の広がり具合を示し、図１５（ｂ）はレンジ１〜Ｎの誤差電圧と目標レンジとの関係を示す図である。図１５から明らかなように、観測レンジに対応したΣとΔビームより、測角用の誤差電圧を予めテーブル化しておき、観測した目標レンジに応じて測角テーブルを切り替えればよい。
【００６４】
図１４に示すブロック図を参照して第４の実施形態を説明する。測角テーブルとしては、Σ系統のセル選定部３Ｅとセル抽出部３Ｆ３でそれぞれ選定・抽出されたセルに基づきＳＡＲ測角値を観測するＳＡＲ可変測角部３Ｍと、ΔＡＺ、ΔＥＬ系統のセル抽出結果に基づいて実開口の可変測角を行う実開口可変測角部３Ｎがあり、各測角部３Ｍ，３ＮはΣ系統のレンジドップラ算出部３８で算出したレンジによりテーブルを切り替える。
【００６５】
尚、第１乃至第４の実施形態において、処理規模を低減するために狭帯域処理部３２により、ビ−ム方向を制御する手法について述べたが、広帯域処理（ＳＡＲ処理）や、他のセンサからの情報等によりビーム方向がわかれば、狭帯域処理を行わなくてもよい。
【００６６】
その他、本実施形態は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。更に、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【符号の説明】
【００６７】
１…アンテナ、１１…送信アンテナ素子、２…送受信器、２１…送受信器、２２…ビーム制御器、２３…サンプル切り替え部、３…信号処理器、３１，３１１，３１２，３１３…ＡＤ変換部、３２…狭帯域処理部、３２１１〜３２１３…データ抽出部、３２２１〜３２２３…ＰＲＩ軸ＦＦＴ処理部、３２３１〜３２３３…レンジ（時間）圧縮部、３２４…ＣＦＡＲ処理部、３２５…スレショルド検出部、３２６１，３２６２…セル検出部、３２７…測角部、３３，３Ａ１，３Ａ２，３Ａ３…データ抽出部、３４…ＰＲＩ軸ＦＦＴ処理部、３５，３Ｂ１，３Ｂ２，３Ｂ３…レンジ圧縮部、３６…ＣＦＡＲ処理部、３７…スレショルド検出部、３８…レンジ・ドップラ算出部、３９…参照信号生成部、３Ｃ１，３Ｃ２，３Ｃ３，３Ｃ４…ＡＺ圧縮部、３Ｄ１，３Ｄ２，３Ｄ３，３Ｄ４…ポーラフォーマット変換による画像化部、３Ｅ…セル選定部、３Ｆ１，３Ｆ２，３Ｆ３…セル抽出部、３Ｇ…測角部、３Ｈ…画像化部、３Ｉ…ΔＳＡＲウェイト算出部、３Ｊ…ＳＡＲ測角部、３Ｋ…実開口測角部、３Ｌ…測角相関部、３Ｍ…ＳＡＲ可変測角部、３Ｎ…実開口可変測角部。
【図１】
0006363524.tif 000020
【図２】
0006363524.tif 000021
【図３】
0006363524.tif 000022
【図４】
0006363524.tif 000023
【図５】
0006363524.tif 000024
【図６】
0006363524.tif 000025
【図７】
0006363524.tif 000026
【図８】
0006363524.tif 000027
【図９】
0006363524.tif 000028
【図１０】
0006363524.tif 000029
【図１１】
0006363524.tif 000030
【図１２】
0006363524.tif 000031
【図１３】
0006363524.tif 000032
【図１４】
0006363524.tif 000033
【図１５】
0006363524.tif 000034