JP6246024B2 - レーダ信号処理装置 - Google Patents
レーダ信号処理装置 Download PDFInfo
- Publication number
- JP6246024B2 JP6246024B2 JP2014038708A JP2014038708A JP6246024B2 JP 6246024 B2 JP6246024 B2 JP 6246024B2 JP 2014038708 A JP2014038708 A JP 2014038708A JP 2014038708 A JP2014038708 A JP 2014038708A JP 6246024 B2 JP6246024 B2 JP 6246024B2
- Authority
- JP
- Japan
- Prior art keywords
- component
- signal
- clutter
- azimuth ambiguity
- region
- 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.)
- Active
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Description
この発明は、合成開口レーダ（ＳＡＲ：Ｓｙｎｔｈｅｔｉｃ Ａｐｅｒｔｕｒｅ Ｒａｄａｒ）の画像であるＳＡＲ画像に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧して、移動目標を検出するレーダ信号処理装置に関するものである。
以下の非特許文献１に開示されているレーダ信号処理装置では、複数の送受信アンテナを軌道方向に具備することで複数のチャンネルを有し、送受信の位相中心が軌道方向で異なる状態で観測された複数のＳＡＲ画像（観測領域が同一で観測時刻が異なる複数のＳＡＲ画像）を取得する。
そして、このレーダ信号処理装置では、複数のＳＡＲ画像間の位置合わせ（レジストレーション）を実施した後に、複数のＳＡＲ画像間の差分をとることで、クラッタ成分を抑圧して移動目標を検出するＭＴＩ（Ｍｏｖｉｎｇ Ｔａｒｇｅｔ Ｉｎｄｉｃａｔｏｒ）処理を実施する。
そして、このレーダ信号処理装置では、複数のＳＡＲ画像間の位置合わせ（レジストレーション）を実施した後に、複数のＳＡＲ画像間の差分をとることで、クラッタ成分を抑圧して移動目標を検出するＭＴＩ（Ｍｏｖｉｎｇ Ｔａｒｇｅｔ Ｉｎｄｉｃａｔｏｒ）処理を実施する。
ただし、複数のＳＡＲ画像には、クラッタ成分の他に、アンテナパターンのサイドローブに起因するアジマスアンビギュイティ成分が発生している場合がある。
このとき、チャンネル間の送受信の位相中心が重複する条件であるＤＰＣＡ（Ｄｉｐｌａｃｅｄ Ｐｈａｓｅ Ｃｅｎｔｅｒ Ａｎｔｅｎｎａ）条件を満足している場合には、複数のＳＡＲ画像間のレジストレーションをピクセル単位で行うことができるため、アジマスアンビギュイティ成分もクラッタ成分と同様に抑圧して移動目標を検出することが可能である。
しかし、ＤＰＣＡ条件を満足していない場合、複数のＳＡＲ画像間のレジストレーションをサブピクセル単位で行うことになる。サブピクセル単位でのリサンプリングを伴うレジストレーションを行う場合、アジマスアンビギュイティ成分はクラッタ成分と異なる振る舞いをするため、レジストレーション後のＳＡＲ画像内のアジマスアンビギュイティ成分に位相回転が発生する。その結果、複数のＳＡＲ画像間の差分をとってもアジマスアンビギュイティ成分が残余する問題が生じる。
このとき、チャンネル間の送受信の位相中心が重複する条件であるＤＰＣＡ（Ｄｉｐｌａｃｅｄ Ｐｈａｓｅ Ｃｅｎｔｅｒ Ａｎｔｅｎｎａ）条件を満足している場合には、複数のＳＡＲ画像間のレジストレーションをピクセル単位で行うことができるため、アジマスアンビギュイティ成分もクラッタ成分と同様に抑圧して移動目標を検出することが可能である。
しかし、ＤＰＣＡ条件を満足していない場合、複数のＳＡＲ画像間のレジストレーションをサブピクセル単位で行うことになる。サブピクセル単位でのリサンプリングを伴うレジストレーションを行う場合、アジマスアンビギュイティ成分はクラッタ成分と異なる振る舞いをするため、レジストレーション後のＳＡＲ画像内のアジマスアンビギュイティ成分に位相回転が発生する。その結果、複数のＳＡＲ画像間の差分をとってもアジマスアンビギュイティ成分が残余する問題が生じる。
ここで、上記の問題を確認するため、非特許文献１に開示されているレーダ信号処理装置の目標検知方式を簡単に説明する。
図９は非特許文献１に開示されているレーダ信号処理装置を示す構成図であり、図１０は図９のレーダ信号処理装置の処理内容を示すフローチャートである。
このレーダ信号処理装置は、１つの送信アンテナからパルスが繰り返し放射され、２つの受信アンテナが目標に反射された上記パルスの反射波を受信することで、２枚のＳＡＲ画像（１）（２）を取得し、２枚のＳＡＲ画像（１）（２）から移動目標を検出するものである。
このとき、２つの受信アンテナ間の距離は、プラットフォーム速度Ｖｐに対して、２ｍＶｐ／ＰＲＦ＋２ＬＶｐ／ＰＲＦで表されるものとする。ＰＦＲ（Ｐｕｌｓｅ Ｒｅｐｅｔｉｔｉｏｎ Ｆｒｅｑｕｅｎｃｙ）はパルス繰返し周波数、ｍ整数、Ｌは実数、０＜Ｌ＜１である。
Ｌ＝０の場合には、図１１に示すように、送受のアンテナ位相中心が重複するＤＰＣＡ条件を満足するが、ここでは、Ｌ≠０であるために、図１２に示すように、ＤＰＣＡ条件を満足しないものとする。
図９は非特許文献１に開示されているレーダ信号処理装置を示す構成図であり、図１０は図９のレーダ信号処理装置の処理内容を示すフローチャートである。
このレーダ信号処理装置は、１つの送信アンテナからパルスが繰り返し放射され、２つの受信アンテナが目標に反射された上記パルスの反射波を受信することで、２枚のＳＡＲ画像（１）（２）を取得し、２枚のＳＡＲ画像（１）（２）から移動目標を検出するものである。
このとき、２つの受信アンテナ間の距離は、プラットフォーム速度Ｖｐに対して、２ｍＶｐ／ＰＲＦ＋２ＬＶｐ／ＰＲＦで表されるものとする。ＰＦＲ（Ｐｕｌｓｅ Ｒｅｐｅｔｉｔｉｏｎ Ｆｒｅｑｕｅｎｃｙ）はパルス繰返し周波数、ｍ整数、Ｌは実数、０＜Ｌ＜１である。
Ｌ＝０の場合には、図１１に示すように、送受のアンテナ位相中心が重複するＤＰＣＡ条件を満足するが、ここでは、Ｌ≠０であるために、図１２に示すように、ＤＰＣＡ条件を満足しないものとする。
まず、レジストレーション部１０１は、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを行う（図１０のステップＳＴ１０１）。
このとき、レジストレーション部１０１は、アジマス方向の（ｍ＋Ｌ）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
このとき、レジストレーション部１０１は、アジマス方向の（ｍ＋Ｌ）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
式（１）において、ｓ（ｐ，ｑ）は、図１３に示すようなクラッタ成分１１１であり、ａ＋（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１２、ａ−（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１３である。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
一方、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）は以下のように記述される。
式（２）において、ｎ２（ｐ，ｑ）はＳＡＲ画像（２）内の雑音成分、ｔ２（ｐ，ｑ）はＳＡＲ画像（２）内の移動目標成分である。
上記のように、ＳＡＲ画像（１）（２）に生じているアジマスアンビギュイティ成分に位相回転が発生するのは、クラッタ成分に対してアジマスアンビギュイティ成分の周波数帯が異なるにもかかわらず、クラッタ成分と同様のアジマスリサンプリングが施されるためである。
ただし、実際にはＰＲＦ未満のドップラー周波数帯域制限を帯域幅ＢＤで行っているため、図１３の下図のような信号成分になる。しかし、この場合も、アンテナのサイドローブパターンに起因するアジマスアンビギュイティ成分は残余している。
ただし、実際にはＰＲＦ未満のドップラー周波数帯域制限を帯域幅ＢＤで行っているため、図１３の下図のような信号成分になる。しかし、この場合も、アンテナのサイドローブパターンに起因するアジマスアンビギュイティ成分は残余している。
差分画像生成部１０２は、ＳＡＲ画像（１）と、レジストレーション後のＳＡＲ画像（２）との差分を取ることで、ＳＡＲ画像（１）（２）に含まれているクラッタ成分を抑圧し（ステップＳＴ１０２）、両画像の差分画像を出力する（ステップＳＴ１０３）。
ＳＡＲ画像（１）と、レジストレーション後のＳＡＲ画像（２）との差分画像の信号成分は、下記の式（３）のように表される。
ＳＡＲ画像（１）と、レジストレーション後のＳＡＲ画像（２）との差分画像の信号成分は、下記の式（３）のように表される。
Ｌ＝０でＤＰＣＡ条件を満足していれば、アジマスアンビギュイティ成分１１２，１１３が残ることなく、差分画像から移動目標成分を確認することが可能である。
しかし、Ｌ≠０でＤＰＣＡ条件を満足していなければ、図１４に示すように、アジマスアンビギュイティ成分が残余するため、移動目標の検出が困難になる。
しかし、Ｌ≠０でＤＰＣＡ条件を満足していなければ、図１４に示すように、アジマスアンビギュイティ成分が残余するため、移動目標の検出が困難になる。
C.E. Livingstone，I. Sikaneta，C.H. Gierull，S. Chiu，A. Beaudoin，J. ampbell，. Beaudoin，S. Gong，and T.A. Knight，"An airborne synthetic aperture radar（SAR）experiment to support RADARSAT-2ground moving target indication（GMTI）"，Can. J. Remote Sensing，Vol. 28，No. 6，pp. 794-813，2002
従来のレーダ信号処理装置は以上のように構成されているので、目標に反射されたパルスの反射波を受信する複数のアンテナの位相中心が重複しているというＤＰＣＡ条件を満足している場合、ＳＡＲ画像（１）と、レジストレーション後のＳＡＲ画像（２）との差分をとれば、アジマスアンビギュイティ成分を除去することができる。しかし、ＤＰＣＡ条件を満足していない場合、ＳＡＲ画像（１）と、レジストレーション後のＳＡＲ画像（２）との差分をとっても、アジマスアンビギュイティ成分が残余するため、移動目標を検出することができなくなることがある課題があった。
この発明は上記のような課題を解決するためになされたもので、ＤＰＣＡ条件を満足していない場合でも、クラッタ成分とアジマスアンビギュイティ成分の両方を抑圧して、移動目標を検出することができるレーダ信号処理装置を得ることを目的とする。
この発明に係るレーダ信号処理装置は、観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別手段と、信号成分分別手段により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と信号成分分別手段により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分をウィナーフィルタを用いてを推定し、正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを設け、
ウィナーフィルタは、信号成分分別手段により正の領域に分別された信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、正の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるとともに、信号成分分別手段により負の領域の信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、負の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるようにしたものである。
ウィナーフィルタは、信号成分分別手段により正の領域に分別された信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、正の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるとともに、信号成分分別手段により負の領域の信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、負の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるようにしたものである。
この発明によれば、位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別手段と、信号成分分別手段により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と信号成分分別手段により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分をウィナーフィルタを用いてを推定して、正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを設け、ウィナーフィルタは、信号成分分別手段により正の領域に分別された信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、正の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるとともに、信号成分分別手段により負の領域の信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、負の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるように構成したので、複数のアンテナの位相中心が重複していない場合でも、クラッタ成分とアジマスアンビギュイティ成分の両方を抑圧して、移動目標を検出することができる効果がある。
実施の形態１．
図１はこの発明の実施の形態１によるレーダ信号処理装置を示す構成図である。
図１において、レジストレーション部１は観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像であるＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。なお、レジストレーション部１は位置合わせ手段を構成している。
図１はこの発明の実施の形態１によるレーダ信号処理装置を示す構成図である。
図１において、レジストレーション部１は観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像であるＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。なお、レジストレーション部１は位置合わせ手段を構成している。
信号成分推定部２は信号成分分別部３と不要成分抑圧部４から構成されており、クラッタ成分及びアジマスアンビギュイティ成分を抑圧したＳＡＲ画像（１）とＳＡＲ画像（２）を出力する処理を実施する。
信号成分分別部３はレジストレーション部１により位置合わせが行われたＳＡＲ画像（１）（２）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する処理を実施する。なお、信号成分分別部３は信号成分分別手段を構成している。
信号成分分別部３はレジストレーション部１により位置合わせが行われたＳＡＲ画像（１）（２）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する処理を実施する。なお、信号成分分別部３は信号成分分別手段を構成している。
不要成分抑圧部４は信号成分分別部３により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、信号成分分別部３により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する処理を実施する。なお、不要成分抑圧部４は不要成分抑圧手段を構成している。
目標検出部５は不要成分抑圧部４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する処理を実施する。なお、目標検出部５は目標検出手段を構成している。
目標検出部５は不要成分抑圧部４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する処理を実施する。なお、目標検出部５は目標検出手段を構成している。
図１の例では、レーダ信号処理装置の構成要素であるレジストレーション部１、信号成分推定部２及び目標検出部５のそれぞれが専用のハードウェア（例えば、ＣＰＵを実装している半導体集積回路、あるいは、ワンチップマイコンなど）で構成されているものを想定しているが、レーダ信号処理装置がコンピュータで構成されていてもよい。
レーダ信号処理装置をコンピュータで構成する場合、レジストレーション部１、信号成分推定部２及び目標検出部５の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのＣＰＵが当該メモリに格納されているプログラムを実行するようにすればよい。
図２はこの発明の実施の形態１によるレーダ信号処理装置の処理内容を示すフローチャートである。
レーダ信号処理装置をコンピュータで構成する場合、レジストレーション部１、信号成分推定部２及び目標検出部５の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのＣＰＵが当該メモリに格納されているプログラムを実行するようにすればよい。
図２はこの発明の実施の形態１によるレーダ信号処理装置の処理内容を示すフローチャートである。
次に動作について説明する。
この実施の形態１では、１つの送信アンテナからパルスが繰り返し放射され、２つの受信アンテナが目標に反射された上記パルスの反射波を受信することで、２枚のＳＡＲ画像（１）（２）を取得し、２枚のＳＡＲ画像（１）（２）から移動目標を検出するものとする。
このとき、２つの受信アンテナ間の距離は、プラットフォーム速度Ｖｐに対して、２ｍＶｐ／ＰＲＦ＋２ＬＶｐ／ＰＲＦで表されるものとする。ＰＦＲ（Ｐｕｌｓｅ Ｒｅｐｅｔｉｔｉｏｎ Ｆｒｅｑｕｅｎｃｙ）はパルス繰返し周波数、ｍ整数、Ｌは実数、−１＜Ｌ＜０又は０＜Ｌ＜１である。
Ｌ＝０の場合には、図１１に示すように、送受のアンテナ位相中心が重複するＤＰＣＡ条件を満足するが、この実施の形態１では、Ｌ≠０であるために、図１２に示すように、ＤＰＣＡ条件を満足しないものとする。
この実施の形態１では、１つの送信アンテナからパルスが繰り返し放射され、２つの受信アンテナが目標に反射された上記パルスの反射波を受信することで、２枚のＳＡＲ画像（１）（２）を取得し、２枚のＳＡＲ画像（１）（２）から移動目標を検出するものとする。
このとき、２つの受信アンテナ間の距離は、プラットフォーム速度Ｖｐに対して、２ｍＶｐ／ＰＲＦ＋２ＬＶｐ／ＰＲＦで表されるものとする。ＰＦＲ（Ｐｕｌｓｅ Ｒｅｐｅｔｉｔｉｏｎ Ｆｒｅｑｕｅｎｃｙ）はパルス繰返し周波数、ｍ整数、Ｌは実数、−１＜Ｌ＜０又は０＜Ｌ＜１である。
Ｌ＝０の場合には、図１１に示すように、送受のアンテナ位相中心が重複するＤＰＣＡ条件を満足するが、この実施の形態１では、Ｌ≠０であるために、図１２に示すように、ＤＰＣＡ条件を満足しないものとする。
レジストレーション部１は、観測領域が同一で観測時刻が異なるＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する（図２のステップＳＴ１）。
このとき、レジストレーション部１は、アジマス方向の（ｍ＋Ｌ）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
このとき、レジストレーション部１は、アジマス方向の（ｍ＋Ｌ）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
式（４）において、ｓ（ｐ，ｑ）は、図１３に示すようなクラッタ成分１１１であり、ａ＋（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１２、ａ−（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１３である。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
一方、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）は以下のように記述される。
式（５）において、ｎ２（ｐ，ｑ）はＳＡＲ画像（２）内の雑音成分、ｔ２（ｐ，ｑ）はＳＡＲ画像（２）内の移動目標成分である。
信号成分推定部２は、レジストレーション部１により位置合わせが行われたＳＡＲ画像（１）とＳＡＲ画像（２）を受けると、そのＳＡＲ画像（１）とＳＡＲ画像（２）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定して、そのクラッタ成分及びアジマスアンビギュイティ成分を抑圧し、クラッタ成分及びアジマスアンビギュイティ成分を抑圧したＳＡＲ画像（１）とＳＡＲ画像（２）を出力する。
以下、信号成分推定部２の処理内容を具体的に説明する。
以下、信号成分推定部２の処理内容を具体的に説明する。
ＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）には、２種類のアジマスアンビギュイティ成分が混在している。
この問題を簡単化するため、信号成分推定部２の信号成分分別部３は、ＳＡＲ画像（１）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別し、ＳＡＲ画像（２）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する（ステップＳＴ２）。
この問題を簡単化するため、信号成分推定部２の信号成分分別部３は、ＳＡＲ画像（１）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別し、ＳＡＲ画像（２）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する（ステップＳＴ２）。
ＳＡＲ画像（１）及びレジストレーション後のＳＡＲ画像（２）における２Ｐ個のドップラー周波数の信号成分は、図１３の下図のように表される。
この信号成分をそれぞれ２つにコピーして、一方の信号成分をドップラー周波数領域の負の領域に０系列を与え、他方の信号成分をドップラー周波数領域の正の領域に０系列を与えると、前者のスペクトルは図３のようになり、後者のスペクトルは図４のようになる。
上記のスペクトルを時間領域に戻すと、ドップラー周波数領域の負の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１＋（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（２）の受信信号ｒ２＋（ｐ，ｑ）は、下記の式（６）（７）のように表され、ドップラー周波数領域の正の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１−（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（２）の受信信号ｒ２−（ｐ，ｑ）は、下記の式（８）（９）のように表される。
この信号成分をそれぞれ２つにコピーして、一方の信号成分をドップラー周波数領域の負の領域に０系列を与え、他方の信号成分をドップラー周波数領域の正の領域に０系列を与えると、前者のスペクトルは図３のようになり、後者のスペクトルは図４のようになる。
上記のスペクトルを時間領域に戻すと、ドップラー周波数領域の負の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１＋（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（２）の受信信号ｒ２＋（ｐ，ｑ）は、下記の式（６）（７）のように表され、ドップラー周波数領域の正の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１−（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（２）の受信信号ｒ２−（ｐ，ｑ）は、下記の式（８）（９）のように表される。
式（６）〜（９）において、クラッタ成分、移動目標成分及び雑音成分に付記されている下字の＋は、受信信号と同様に、ドップラー周波数領域の負の領域に０系列が与えられて得られる信号であることを表し、下字の−は、ドップラー周波数領域の正の領域に０系列が与えられて得られる信号であることを表している。このような変換を行うことで、推定問題を簡単化することができる。
不要成分抑圧部４は、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ）を用いて、なるべく移動目標成分が含まれないように、クラッタ成分ｓ＋（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定するとともに、受信信号ｒ１−（ｐ，ｑ），ｒ２−（ｐ，ｑ）を用いて、なるべく移動目標成分が含まれないように、クラッタ成分ｓ−（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定する（ステップＳＴ３）。
以下、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ）を用いて、クラッタ成分ｓ＋（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定する処理を説明する。
以下、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ）を用いて、クラッタ成分ｓ＋（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定する処理を説明する。
まず、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ）の定式化を整理する。
受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ）の定式を行列とベクトルで表すと、以下のように与えられる。
受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ）の定式を行列とベクトルで表すと、以下のように与えられる。
式（１０）のように定式化された信号の推定方法の１つとして、最小二乗法が知られている（例えば、非特許文献２）。
［非特許文献２］
M. K. Steven，“Fundamentals of Statistical Signal Processing : Estimation Theory”，Prentice Hall PTR Prentice-Hall，Inc.，1993.
最小二乗法による具体的な推定方法を式で表すと、下記の式（１６）のようになる。
［非特許文献２］
M. K. Steven，“Fundamentals of Statistical Signal Processing : Estimation Theory”，Prentice Hall PTR Prentice-Hall，Inc.，1993.
最小二乗法による具体的な推定方法を式で表すと、下記の式（１６）のようになる。
また、式（１６）において、−１の添え字は逆行列を表し、Ｈの添え字は共役転置を表している。
しかし、最小二乗法で推定する場合、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）以外の雑音成分ｎ１＋（ｐ，ｑ），ｎ２＋（ｐ，ｑ）や移動目標成分ｔ１＋（ｐ，ｑ）の存在を仮定していない。
このため、例えば、受信信号ｒ１＋（ｐ，ｑ）に含まれているクラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を抑圧するために、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を用いて、ｒ１＋（ｐ，ｑ）−ｓ＋（ｐ，ｑ）−ａ＋（ｐ，ｑ）を計算しても０となり、移動目標成分ｔ１＋（ｐ，ｑ）が残余しない。
したがって、最小二乗法で推定する方式は適用することができず、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）以外の信号の存在を仮定する推定方法を用いる必要がある。
そこで、この実施の形態１では、ウィナーフィルタを用いて、信号の推定を行うものとする。
このため、例えば、受信信号ｒ１＋（ｐ，ｑ）に含まれているクラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を抑圧するために、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を用いて、ｒ１＋（ｐ，ｑ）−ｓ＋（ｐ，ｑ）−ａ＋（ｐ，ｑ）を計算しても０となり、移動目標成分ｔ１＋（ｐ，ｑ）が残余しない。
したがって、最小二乗法で推定する方式は適用することができず、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）以外の信号の存在を仮定する推定方法を用いる必要がある。
そこで、この実施の形態１では、ウィナーフィルタを用いて、信号の推定を行うものとする。
ウィナーフィルタを用いる推定方法を式で表すと、下記の式（１７）のようになる。
ただし、Ｅ（ａ＋／ｓ＋）はアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値であり、Ｅ（ｎ＋／ｓ＋）は雑音成分ｎ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値である。
この推定方式は、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）以外の信号の存在を仮定し、それらの信号の存在が、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）の推定の際に含まれる量を統計的に最小化した方式である。
したがって、ｒ１＋（ｐ，ｑ）−ｓ＋（ｐ，ｑ）−ａ＋（ｐ，ｑ）を計算した場合には、移動目標成分ｔ１＋（ｐ，ｑ）が残余する。
したがって、ｒ１＋（ｐ，ｑ）−ｓ＋（ｐ，ｑ）−ａ＋（ｐ，ｑ）を計算した場合には、移動目標成分ｔ１＋（ｐ，ｑ）が残余する。
受信信号ｒ１−（ｐ，ｑ），ｒ２−（ｐ，ｑ）についても同様に、行列とベクトルで表すと、以下のように与えられる。
また、ウィナーフィルタを用いる推定方法は、下記の式（２６）のようになる。
不要成分抑圧部４は、ウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定すると、それらの推定結果を用いて、ＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、ＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する（ステップＳＴ４）。
クラッタ成分及びアジマスアンビギュイティ成分の抑圧方法は、いくつかのバリエーションが考えられる。
最も単純な抑圧方法としては、以下の計算を行うものがある。
クラッタ成分及びアジマスアンビギュイティ成分の抑圧方法は、いくつかのバリエーションが考えられる。
最も単純な抑圧方法としては、以下の計算を行うものがある。
式（２９）の計算結果又は式（３０）の計算結果を用いれば、移動目標を検出することが可能であるが、下記の式（３１）に示すように、２チャンネルの成分を合成すれば、より移動目標の検出性能を高めることができる。
なお、不要成分抑圧部４により推定されたクラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）には推定誤差が含まれている。
特に、クラッタ成分の推定誤差は、クラッタ成分の電力が大きいことから、移動目標の検出性能劣化に大きく寄与すると考えられる。
そこで、不要成分抑圧部４により推定されたクラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）を用いずにクラッタを除去し、残ったアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を抑圧する下記の処理が考えられる。
特に、クラッタ成分の推定誤差は、クラッタ成分の電力が大きいことから、移動目標の検出性能劣化に大きく寄与すると考えられる。
そこで、不要成分抑圧部４により推定されたクラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）を用いずにクラッタを除去し、残ったアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を抑圧する下記の処理が考えられる。
上記の不要成分抑圧部４によるクラッタ成分及びアジマスアンビギュイティ成分の抑圧は、振幅の次元で行っているが、下記に示すように、電力の次元で行ってから、振幅の次元に戻すことで、移動目標の検出性能を高めることも可能である。
図５は不要成分抑圧部４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像を示す説明図である。
図５の例では、４個の移動目標成分が表されており、クラッタ成分及びアジマスアンビギュイティ成分が抑圧されている。
目標検出部５は、不要成分抑圧部４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像（以下、「抑圧後画像」と称する）から移動目標を検出する（ステップＳＴ５）。
移動目標の検出処理としては、閾値判定処理などが考えられるが、抑圧後画像を電力成分に変換してから閾値判定処理を行うようにしてもよい。また、雑音を低減するために、移動平均フィルタなどを施してから閾値判定処理を行うようにしてもよい。
図５の例では、４個の移動目標成分が表されており、クラッタ成分及びアジマスアンビギュイティ成分が抑圧されている。
目標検出部５は、不要成分抑圧部４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像（以下、「抑圧後画像」と称する）から移動目標を検出する（ステップＳＴ５）。
移動目標の検出処理としては、閾値判定処理などが考えられるが、抑圧後画像を電力成分に変換してから閾値判定処理を行うようにしてもよい。また、雑音を低減するために、移動平均フィルタなどを施してから閾値判定処理を行うようにしてもよい。
また、ＳＡＲ画像（１）とＳＡＲ画像（２）を干渉させて、移動目標の速度推定を行うＡＴＩ（Ａｌｏｎｇ Ｔｒａｃｋ Ｉｎｔｅｒｆｅｒｏｍｅｔｒｙ）にも用いることが可能である。
具体的には、以下の偏角による位相マップを計算する。
式（３７）において、＊は共役を示す記号である。
具体的には、以下の偏角による位相マップを計算する。
以上で明らかなように、この実施の形態１によれば、レジストレーション部１により位置合わせが行われたＳＡＲ画像（１）（２）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別部３と、信号成分分別部３により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、信号成分分別部３により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧部４とを設けるように構成したので、複数のアンテナの位相中心が重複していない場合でも、クラッタ成分とアジマスアンビギュイティ成分の両方を抑圧して、移動目標を検出することができる効果を奏する。
この実施の形態１では、上記の式（１８）（１９）に示すように、アジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値Ｅ（ａ＋／ｓ＋）と、雑音成分ｎ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値Ｅ（ｎ＋／ｓ＋）とを用いて、ウィナーフィルタを生成するとともに、式（２７）（２８）に示すように、アジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比の期待値Ｅ（ａ−／ｓ−）と、雑音成分ｎ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比の期待値Ｅ（ｎ−／ｓ−）とを用いて、ウィナーフィルタを生成し、それらのウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、その推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、不要成分抑圧部４は、期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比（ａ＋／ｓ＋）と、その推定したクラッタ成分ｓ＋（ｐ，ｑ）と雑音成分ｎ＋（ｐ，ｑ）の電力比（ｎ＋／ｓ＋）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定するようにする。
同様に、期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
同様に、期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部４は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
なお、クラッタ成分及びアジマスアンビギュイティ成分の収束判定としては、例えば、前回の推定結果と今回推定結果を比較して、それらの推定結果の差分が予め設定された閾値以内であれば、収束していると判定する方法が考えられる。
なお、クラッタ成分及びアジマスアンビギュイティ成分の収束判定としては、例えば、前回の推定結果と今回推定結果を比較して、それらの推定結果の差分が予め設定された閾値以内であれば、収束していると判定する方法が考えられる。
この実施の形態１では、不要成分抑圧部４が、期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）を用いて、ウィナーフィルタを生成するものを示したが、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、クラッタ成分の信号電力を用いて、ウィナーフィルタを生成するようにしてもよい。
即ち、上記の式（１８）（１９）における期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）、，式（２７）（２８）における期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）の代わりに、クラッタ成分の信号電力に反比例する値を代入して、ウィナーフィルタを生成するようにしてもよい。
この場合も、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、上記の式（１８）（１９）における期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）、，式（２７）（２８）における期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）の代わりに、クラッタ成分の信号電力に反比例する値を代入して、ウィナーフィルタを生成するようにしてもよい。
この場合も、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、不要成分抑圧部４は、クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比（ａ＋／ｓ＋）と、その推定したクラッタ成分ｓ＋（ｐ，ｑ）と雑音成分ｎ＋（ｐ，ｑ）の電力比（ｎ＋／ｓ＋）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定するようにする。
同様に、クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部４は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
同様に、クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部４は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
この実施の形態１では、不要成分抑圧部４が、ウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、正則化最小二乗法を行うフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するようにしてもよい。正則化最小二乗法を行うフィルタは、下記の式（３８）のように表される。
式（３８）において、Ｉは単位行列、λは任意の実数である。
ここでは、λが任意の実数であるものを示したが、λの代わりに、雑音成分ｎ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値Ｅ（ｎ＋／ｓ＋）（あるいは、雑音成分ｎ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比の期待値Ｅ（ｎ−／ｓ−））を式（３８）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
また、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、λの代わりに、クラッタ成分の信号電力を式（３８）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
また、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、λの代わりに、クラッタ成分の信号電力を式（３８）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
この実施の形態１では、不要成分抑圧部４が、ウィナーフィルタ又は正則化最小二乗法を行うフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、これに限るものではなく、推定したい成分の他の信号成分の存在を仮定するフィルタや推定方式を適用することが可能である。
また、この実施の形態１では、２本の受信アンテナを用いる例を示したが、複数本の送信アンテナを用いる場合にも対応可能である。
また、この実施の形態１では、２本の受信アンテナを用いる例を示したが、複数本の送信アンテナを用いる場合にも対応可能である。
実施の形態２．
この実施の形態２では、１つの送信アンテナからパルスが繰り返し放射され、３つの受信アンテナが目標に反射された上記パルスの反射波を受信することで、３枚のＳＡＲ画像（１）（２）（３）を取得し、３枚のＳＡＲ画像（１）（２）（３）から移動目標を検出するレーダ信号処理装置について説明する。
このとき、プラットフォーム速度Ｖｐに対して、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（２）に対応する受信アンテナとの距離は、２ｍ２Ｖｐ／ＰＲＦ＋２Ｌ２Ｖｐ／ＰＲＦで表されるものとする。また、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（３）に対応する受信アンテナとの距離は、２ｍ３Ｖｐ／ＰＲＦ＋２Ｌ３Ｖｐ／ＰＲＦで表されるものとする。
ＰＦＲはパルス繰返し周波数、ｍ２，ｍ３は整数、Ｌ２，Ｌ３は実数、０＜Ｌ２、Ｌ３＜１である。
この実施の形態２では、１つの送信アンテナからパルスが繰り返し放射され、３つの受信アンテナが目標に反射された上記パルスの反射波を受信することで、３枚のＳＡＲ画像（１）（２）（３）を取得し、３枚のＳＡＲ画像（１）（２）（３）から移動目標を検出するレーダ信号処理装置について説明する。
このとき、プラットフォーム速度Ｖｐに対して、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（２）に対応する受信アンテナとの距離は、２ｍ２Ｖｐ／ＰＲＦ＋２Ｌ２Ｖｐ／ＰＲＦで表されるものとする。また、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（３）に対応する受信アンテナとの距離は、２ｍ３Ｖｐ／ＰＲＦ＋２Ｌ３Ｖｐ／ＰＲＦで表されるものとする。
ＰＦＲはパルス繰返し周波数、ｍ２，ｍ３は整数、Ｌ２，Ｌ３は実数、０＜Ｌ２、Ｌ３＜１である。
図６はこの発明の実施の形態２によるレーダ信号処理装置を示す構成図である。
図６において、レジストレーション部１１ａは観測領域が同一で観測時刻が異なるＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
レジストレーション部１１ｂは観測領域が同一で観測時刻が異なるＳＡＲ画像（１）とＳＡＲ画像（３）が与えられると、ＳＡＲ画像（３）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。なお、レジストレーション部１１ａ，１１ｂは位置合わせ手段を構成している。
図６において、レジストレーション部１１ａは観測領域が同一で観測時刻が異なるＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
レジストレーション部１１ｂは観測領域が同一で観測時刻が異なるＳＡＲ画像（１）とＳＡＲ画像（３）が与えられると、ＳＡＲ画像（３）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。なお、レジストレーション部１１ａ，１１ｂは位置合わせ手段を構成している。
信号成分推定部１２は信号成分分別部１３と不要成分抑圧部１４から構成されており、クラッタ成分及びアジマスアンビギュイティ成分を抑圧したＳＡＲ画像（１）（２）（３）を出力する処理を実施する。
信号成分分別部１３はレジストレーション部１１ａ，１１ｂにより位置合わせが行われたＳＡＲ画像（１）（２）（３）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する処理を実施する。なお、信号成分分別部１３は信号成分分別手段を構成している。
信号成分分別部１３はレジストレーション部１１ａ，１１ｂにより位置合わせが行われたＳＡＲ画像（１）（２）（３）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する処理を実施する。なお、信号成分分別部１３は信号成分分別手段を構成している。
不要成分抑圧部１４は信号成分分別部１３により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、信号成分分別部１３により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する処理を実施する。なお、不要成分抑圧部１４は不要成分抑圧手段を構成している。
目標検出部１５は不要成分抑圧部１４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する処理を実施する。なお、目標検出部１５は目標検出手段を構成している。
目標検出部１５は不要成分抑圧部１４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する処理を実施する。なお、目標検出部１５は目標検出手段を構成している。
図６の例では、レーダ信号処理装置の構成要素であるレジストレーション部１１ａ，１１ｂ、信号成分推定部１２及び目標検出部１５のそれぞれが専用のハードウェア（例えば、ＣＰＵを実装している半導体集積回路、あるいは、ワンチップマイコンなど）で構成されているものを想定しているが、レーダ信号処理装置がコンピュータで構成されていてもよい。
レーダ信号処理装置をコンピュータで構成する場合、レジストレーション部１１ａ，１１ｂ、信号成分推定部１２及び目標検出部１５の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのＣＰＵが当該メモリに格納されているプログラムを実行するようにすればよい。
レーダ信号処理装置をコンピュータで構成する場合、レジストレーション部１１ａ，１１ｂ、信号成分推定部１２及び目標検出部１５の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのＣＰＵが当該メモリに格納されているプログラムを実行するようにすればよい。
次に動作について説明する。
この実施の形態２では、観測領域が同一で観測時刻が異なる３枚のＳＡＲ画像（１）（２）（３）が与えられると、ＳＡＲ画像（１）を基準にして、以下の処理を実施する。
レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
また、レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（３）が与えられると、ＳＡＲ画像（３）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
このとき、レジストレーション部１１ａは、アジマス方向の（ｍ２＋Ｌ２）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
また、レジストレーション部１１ｂは、アジマス方向の（ｍ３＋Ｌ３）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
この実施の形態２では、観測領域が同一で観測時刻が異なる３枚のＳＡＲ画像（１）（２）（３）が与えられると、ＳＡＲ画像（１）を基準にして、以下の処理を実施する。
レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
また、レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（３）が与えられると、ＳＡＲ画像（３）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
このとき、レジストレーション部１１ａは、アジマス方向の（ｍ２＋Ｌ２）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
また、レジストレーション部１１ｂは、アジマス方向の（ｍ３＋Ｌ３）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
式（３９）において、ｓ（ｐ，ｑ）は、図１３に示すようなクラッタ成分１１１であり、ａ＋（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１２、ａ−（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１３である。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
一方、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）は以下のように記述される。
式（４０）において、ｎ２（ｐ，ｑ）はＳＡＲ画像（２）内の雑音成分、ｔ２（ｐ，ｑ）はＳＡＲ画像（２）内の移動目標成分である。
また、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）は以下のように記述される。
式（４１）において、ｎ３（ｐ，ｑ）はＳＡＲ画像（３）内の雑音成分、ｔ３（ｐ，ｑ）はＳＡＲ画像（３）内の移動目標成分である。
信号成分推定部１２は、レジストレーション部１１ａ，１１ｂにより位置合わせが行われたＳＡＲ画像（１）（２）（３）を受けると、そのＳＡＲ画像（１）（２）（３）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定して、そのクラッタ成分及びアジマスアンビギュイティ成分を抑圧し、クラッタ成分及びアジマスアンビギュイティ成分を抑圧したＳＡＲ画像（１）（２）（３）を出力する。
以下、信号成分推定部１２の処理内容を具体的に説明する。
以下、信号成分推定部１２の処理内容を具体的に説明する。
ＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）、レジストレーション後のＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）には、２種類のアジマスアンビギュイティ成分が混在している。
この問題を簡単化するため、信号成分推定部１２の信号成分分別部１３は、ＳＡＲ画像（１）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別し、ＳＡＲ画像（２）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別し、ＳＡＲ画像（３）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する。
この問題を簡単化するため、信号成分推定部１２の信号成分分別部１３は、ＳＡＲ画像（１）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別し、ＳＡＲ画像（２）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別し、ＳＡＲ画像（３）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する。
例えば、ＳＡＲ画像（１）及びレジストレーション後のＳＡＲ画像（２）における２Ｐ個のドップラー周波数の信号成分は、図１３の下図のように表される。
この信号成分をそれぞれ２つにコピーして、一方の信号成分をドップラー周波数領域の負の領域に０系列を与え、他方の信号成分をドップラー周波数領域の正の領域に０系列を与えると、前者のスペクトルは図３のようになり、後者のスペクトルは図４のようになる。
上記のスペクトルを時間領域に戻すと、ドップラー周波数領域の負の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１＋（ｐ，ｑ）、レジストレーション後のＳＡＲ画像（２）の受信信号ｒ２＋（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（３）の受信信号ｒ３＋（ｐ，ｑ）は、下記の式（４２）〜（４４）のように表され、ドップラー周波数領域の正の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１−（ｐ，ｑ）、レジストレーション後のＳＡＲ画像（２）の受信信号ｒ２−（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（３）の受信信号ｒ３−（ｐ，ｑ）は、下記の式（４５）〜（４７）のように表される。
この信号成分をそれぞれ２つにコピーして、一方の信号成分をドップラー周波数領域の負の領域に０系列を与え、他方の信号成分をドップラー周波数領域の正の領域に０系列を与えると、前者のスペクトルは図３のようになり、後者のスペクトルは図４のようになる。
上記のスペクトルを時間領域に戻すと、ドップラー周波数領域の負の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１＋（ｐ，ｑ）、レジストレーション後のＳＡＲ画像（２）の受信信号ｒ２＋（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（３）の受信信号ｒ３＋（ｐ，ｑ）は、下記の式（４２）〜（４４）のように表され、ドップラー周波数領域の正の領域に０系列が与えたられたＳＡＲ画像（１）の受信信号ｒ１−（ｐ，ｑ）、レジストレーション後のＳＡＲ画像（２）の受信信号ｒ２−（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（３）の受信信号ｒ３−（ｐ，ｑ）は、下記の式（４５）〜（４７）のように表される。
式（４２）〜（４７）において、クラッタ成分、移動目標成分及び雑音成分に付記されている下字の＋は、受信信号と同様に、ドップラー周波数領域の負の領域に０系列が与えられて得られる信号であることを表し、下字の−は、ドップラー周波数領域の正の領域に０系列が与えられて得られる信号であることを表している。このような変換を行うことで、推定問題を簡単化することができる。
不要成分抑圧部１４は、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ），ｒ３＋（ｐ，ｑ）を用いて、なるべく移動目標成分が含まれないように、クラッタ成分ｓ＋（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定するとともに、受信信号ｒ１−（ｐ，ｑ），ｒ２−（ｐ，ｑ），ｒ３−（ｐ，ｑ）を用いて、なるべく移動目標成分が含まれないように、クラッタ成分ｓ−（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定する。
以下、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ），ｒ３＋（ｐ，ｑ）を用いて、クラッタ成分ｓ＋（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定する処理を説明する。
以下、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ），ｒ３＋（ｐ，ｑ）を用いて、クラッタ成分ｓ＋（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定する処理を説明する。
まず、受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ），ｒ３＋（ｐ，ｑ）の定式化を整理する。
受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ），ｒ３＋（ｐ，ｑ）の定式を行列とベクトルで表すと、以下のように与えられる。
受信信号ｒ１＋（ｐ，ｑ），ｒ２＋（ｐ，ｑ），ｒ３＋（ｐ，ｑ）の定式を行列とベクトルで表すと、以下のように与えられる。
この実施の形態２では、上記実施の形態１と同様の理由で、ウィナーフィルタを用いて、信号の推定を行う。
ウィナーフィルタを用いる推定方法を式で表すと、下記の式（５４）のようになる。
ただし、Ｅ（ａ＋／ｓ＋）はアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値であり、Ｅ（ｎ＋／ｓ＋）は雑音成分ｎ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値である。
ウィナーフィルタを用いる推定方法を式で表すと、下記の式（５４）のようになる。
この推定方式は、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）以外の信号の存在を仮定し、それらの信号の存在が、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）の推定の際に含まれる量を統計的に最小化した方式である。
したがって、ｒ１＋（ｐ，ｑ）−ｓ＋（ｐ，ｑ）−ａ＋（ｐ，ｑ）を計算した場合には、移動目標成分ｔ１＋（ｐ，ｑ）が残余する。
したがって、ｒ１＋（ｐ，ｑ）−ｓ＋（ｐ，ｑ）−ａ＋（ｐ，ｑ）を計算した場合には、移動目標成分ｔ１＋（ｐ，ｑ）が残余する。
受信信号ｒ１−（ｐ，ｑ），ｒ２−（ｐ，ｑ），ｒ３−（ｐ，ｑ）についても同様に、行列とベクトルで表すと、以下のように与えられる。
また、ウィナーフィルタを用いる推定方法は、下記の式（６３）のようになる。
不要成分抑圧部１４は、ウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定すると、それらの推定結果を用いて、ＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、ＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する。
また、ＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する。
最も単純な抑圧方法としては、以下の計算を行うものがある。
また、ＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する。
最も単純な抑圧方法としては、以下の計算を行うものがある。
式（６６）の計算結果、式（６７）の計算結果又は式（６８）の計算結果を用いれば、移動目標を検出することが可能であるが、下記の式（６９）に示すように、３チャンネルの成分を合成すれば、より移動目標の検出性能を高めることができる。
なお、不要成分抑圧部１４により推定されたクラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）には推定誤差が含まれている。
特に、クラッタ成分の推定誤差は、クラッタ成分の電力が大きいことから、移動目標の検出性能劣化に大きく寄与すると考えられる。
そこで、不要成分抑圧部１４により推定されたクラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）を用いずにクラッタを除去し、残ったアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を抑圧する下記の処理が考えられる。
特に、クラッタ成分の推定誤差は、クラッタ成分の電力が大きいことから、移動目標の検出性能劣化に大きく寄与すると考えられる。
そこで、不要成分抑圧部１４により推定されたクラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）を用いずにクラッタを除去し、残ったアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を抑圧する下記の処理が考えられる。
上記の不要成分抑圧部１４によるクラッタ成分及びアジマスアンビギュイティ成分の抑圧は、振幅の次元で行っているが、下記に示すように、電力の次元で行ってから、振幅の次元に戻すことで、移動目標の検出性能を高めることも可能である。
目標検出部１５は、不要成分抑圧部１４によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像である抑圧後画像から移動目標を検出する。
移動目標の検出処理としては、閾値判定処理などが考えられるが、抑圧後画像を電力成分に変換してから閾値判定処理を行うようにしてもよい。また、雑音を低減するために、移動平均フィルタなどを施してから閾値判定処理を行うようにしてもよい。
移動目標の検出処理としては、閾値判定処理などが考えられるが、抑圧後画像を電力成分に変換してから閾値判定処理を行うようにしてもよい。また、雑音を低減するために、移動平均フィルタなどを施してから閾値判定処理を行うようにしてもよい。
また、ＳＡＲ画像（１）とＳＡＲ画像（２）、ＳＡＲ画像（１）とＳＡＲ画像（３）、ＳＡＲ画像（２）とＳＡＲ画像（３）を干渉させて、移動目標の速度推定を行うＡＴＩにも用いることが可能である。
具体的には、以下の偏角による位相マップを計算する。
式（７６）〜（７８）において、＊は共役を示す記号である。
具体的には、以下の偏角による位相マップを計算する。
式（７６）〜（７８）において、＊は共役を示す記号である。
以上で明らかなように、この実施の形態２によれば、レジストレーション部１１ａ，１１ｂにより位置合わせが行われたＳＡＲ画像（１）（２）（３）におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別部１３と、信号成分分別部１３により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、信号成分分別部１３により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、そのクラッタ成分及びアジマスアンビギュイティ成分を推定することで、負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧部１４とを設けるように構成したので、複数のアンテナの位相中心が重複していない場合でも、クラッタ成分とアジマスアンビギュイティ成分の両方を抑圧して、移動目標を検出することができる効果を奏する。
また、上記実施の形態１よりも、チャンネル数が多いので、位相のアンラップを解くことで、速度推定範囲を拡大することができる効果を奏する。
また、上記実施の形態１よりも、チャンネル数が多いので、位相のアンラップを解くことで、速度推定範囲を拡大することができる効果を奏する。
この実施の形態２では、上記の式（５５）（５６）に示すように、アジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値Ｅ（ａ＋／ｓ＋）と、雑音成分ｎ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値Ｅ（ｎ＋／ｓ＋）とを用いて、ウィナーフィルタを生成するとともに、式（６４）（６５）に示すように、アジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比の期待値Ｅ（ａ−／ｓ−）と、雑音成分ｎ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比の期待値Ｅ（ｎ−／ｓ−）とを用いて、ウィナーフィルタを生成し、それらのウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、その推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、不要成分抑圧部１４は、期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比（ａ＋／ｓ＋）と、その推定したクラッタ成分ｓ＋（ｐ，ｑ）と雑音成分ｎ＋（ｐ，ｑ）の電力比（ｎ＋／ｓ＋）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定するようにする。
同様に、期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
同様に、期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部１４は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
なお、クラッタ成分及びアジマスアンビギュイティ成分の収束判定としては、例えば、前回の推定結果と今回推定結果を比較して、それらの推定結果の差分が予め設定された閾値以内であれば、収束していると判定する方法が考えられる。
なお、クラッタ成分及びアジマスアンビギュイティ成分の収束判定としては、例えば、前回の推定結果と今回推定結果を比較して、それらの推定結果の差分が予め設定された閾値以内であれば、収束していると判定する方法が考えられる。
この実施の形態２では、不要成分抑圧部１４が、期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）を用いて、ウィナーフィルタを生成するものを示したが、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、クラッタ成分の信号電力を用いて、ウィナーフィルタを生成するようにしてもよい。
即ち、上記の式（５５）（５６）における期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）、，式（６４）（６５）における期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）の代わりに、クラッタ成分の信号電力に反比例する値を代入して、ウィナーフィルタを生成するようにしてもよい。
この場合も、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、上記の式（５５）（５６）における期待値Ｅ（ａ＋／ｓ＋），Ｅ（ｎ＋／ｓ＋）、，式（６４）（６５）における期待値Ｅ（ａ−／ｓ−），Ｅ（ｎ−／ｓ−）の代わりに、クラッタ成分の信号電力に反比例する値を代入して、ウィナーフィルタを生成するようにしてもよい。
この場合も、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、不要成分抑圧部１４は、クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比（ａ＋／ｓ＋）と、その推定したクラッタ成分ｓ＋（ｐ，ｑ）と雑音成分ｎ＋（ｐ，ｑ）の電力比（ｎ＋／ｓ＋）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）を推定するようにする。
同様に、クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部１４は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
同様に、クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比（ａ−／ｓ−）と、その推定したクラッタ成分ｓ−（ｐ，ｑ）と雑音成分ｎ−（ｐ，ｑ）の電力比（ｎ−／ｓ−）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部１４は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
この実施の形態２では、不要成分抑圧部１４が、ウィナーフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、正則化最小二乗法を行うフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するようにしてもよい。正則化最小二乗法を行うフィルタは、下記の式（７９）のように表される。
式（７９）において、Ｉは単位行列、λは任意の実数である。
式（７９）において、Ｉは単位行列、λは任意の実数である。
ここでは、λが任意の実数であるものを示したが、λの代わりに、雑音成分ｎ＋（ｐ，ｑ）とクラッタ成分ｓ＋（ｐ，ｑ）の電力比の期待値Ｅ（ｎ＋／ｓ＋）（あるいは、雑音成分ｎ−（ｐ，ｑ）とクラッタ成分ｓ−（ｐ，ｑ）の電力比の期待値Ｅ（ｎ−／ｓ−））を式（７９）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
また、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、λの代わりに、クラッタ成分の信号電力を式（７９）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
また、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、λの代わりに、クラッタ成分の信号電力を式（７９）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
この実施の形態２では、不要成分抑圧部１４が、ウィナーフィルタ又は正則化最小二乗法を行うフィルタを用いて、クラッタ成分ｓ＋（ｐ，ｑ），ｓ−（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、これに限るものではなく、推定したい成分の他の信号成分の存在を仮定するフィルタや推定方式を適用することが可能である。
また、この実施の形態２では、３本の受信アンテナを用いる例を示したが、４本以上の受信アンテナを用いる場合にも適用可能である。また、複数本の送信アンテナを用いる場合や、複数本の送信アンテナ及び複数本の受信アンテナを用いる場合にも対応可能である。
また、この実施の形態２では、３本の受信アンテナを用いる例を示したが、４本以上の受信アンテナを用いる場合にも適用可能である。また、複数本の送信アンテナを用いる場合や、複数本の送信アンテナ及び複数本の受信アンテナを用いる場合にも対応可能である。
実施の形態３．
この実施の形態３では、上記実施の形態２と同様に、３本の受信アンテナを用いる例を説明するが、図７に示すように、ドップラー周波数領域におけるアジマスアンビギュイティ成分の影響が大きく、ドップラー周波数領域で信号成分の分別を行っても、２種類のアジマスアンビギュイティ成分が混在してしまうことである。
この実施の形態３では、図７に示すように、２種類のアジマスアンビギュイティ成分が混在してしまう場合でも、クラッタ成分とアジマスアンビギュイティ成分の両方を抑圧して、移動目標を検出することができるレーダ信号処理装置を説明する。
この実施の形態３では、上記実施の形態２と同様に、３本の受信アンテナを用いる例を説明するが、図７に示すように、ドップラー周波数領域におけるアジマスアンビギュイティ成分の影響が大きく、ドップラー周波数領域で信号成分の分別を行っても、２種類のアジマスアンビギュイティ成分が混在してしまうことである。
この実施の形態３では、図７に示すように、２種類のアジマスアンビギュイティ成分が混在してしまう場合でも、クラッタ成分とアジマスアンビギュイティ成分の両方を抑圧して、移動目標を検出することができるレーダ信号処理装置を説明する。
図８はこの発明の実施の形態３によるレーダ信号処理装置を示す構成図であり、図において、図６と同一符号は同一または相当部分を示すので説明を省略する。
不要成分抑圧部２０はレジストレーション部１１ａ，１１ｂにより位置合わせが行われたＳＡＲ画像（１）（２）（３）におけるドップラー周波数領域の正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定することで、正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、そのドップラー周波数領域の負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定することで、負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する処理を実施する。なお、不要成分抑圧部２０は不要成分抑圧手段を構成している。
不要成分抑圧部２０はレジストレーション部１１ａ，１１ｂにより位置合わせが行われたＳＡＲ画像（１）（２）（３）におけるドップラー周波数領域の正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定することで、正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、そのドップラー周波数領域の負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定することで、負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する処理を実施する。なお、不要成分抑圧部２０は不要成分抑圧手段を構成している。
この実施の形態３では、１つの送信アンテナからパルスが繰り返し放射され、３つの受信アンテナが目標に反射された上記パルスの反射波を受信することで、３枚のＳＡＲ画像（１）（２）（３）を取得し、３枚のＳＡＲ画像（１）（２）（３）から移動目標を検出するものとする。
このとき、プラットフォーム速度Ｖｐに対して、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（２）に対応する受信アンテナとの距離は、２ｍ２Ｖｐ／ＰＲＦ＋２Ｌ２Ｖｐ／ＰＲＦで表されるものとする。また、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（３）に対応する受信アンテナとの距離は、２ｍ３Ｖｐ／ＰＲＦ＋２Ｌ３Ｖｐ／ＰＲＦで表されるものとする。
ＰＦＲはパルス繰返し周波数、ｍ２，ｍ３は整数、Ｌ２，Ｌ３は実数、０＜Ｌ２、Ｌ３＜１である。
このとき、プラットフォーム速度Ｖｐに対して、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（２）に対応する受信アンテナとの距離は、２ｍ２Ｖｐ／ＰＲＦ＋２Ｌ２Ｖｐ／ＰＲＦで表されるものとする。また、ＳＡＲ画像（１）に対応する受信アンテナと、ＳＡＲ画像（３）に対応する受信アンテナとの距離は、２ｍ３Ｖｐ／ＰＲＦ＋２Ｌ３Ｖｐ／ＰＲＦで表されるものとする。
ＰＦＲはパルス繰返し周波数、ｍ２，ｍ３は整数、Ｌ２，Ｌ３は実数、０＜Ｌ２、Ｌ３＜１である。
次に動作について説明する。
この実施の形態３では、観測領域が同一で観測時刻が異なる３枚のＳＡＲ画像（１）（２）（３）が与えられると、ＳＡＲ画像（１）を基準にして、以下の処理を実施する。
レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
また、レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（３）が与えられると、ＳＡＲ画像（３）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
このとき、レジストレーション部１１ａは、アジマス方向の（ｍ２＋Ｌ２）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
また、レジストレーション部１１ｂは、アジマス方向の（ｍ３＋Ｌ３）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
この実施の形態３では、観測領域が同一で観測時刻が異なる３枚のＳＡＲ画像（１）（２）（３）が与えられると、ＳＡＲ画像（１）を基準にして、以下の処理を実施する。
レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（２）が与えられると、ＳＡＲ画像（２）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
また、レジストレーション部１１ａは、ＳＡＲ画像（１）とＳＡＲ画像（３）が与えられると、ＳＡＲ画像（３）の位置をＳＡＲ画像（１）の位置に合わせるレジストレーションを実施する。
このとき、レジストレーション部１１ａは、アジマス方向の（ｍ２＋Ｌ２）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
また、レジストレーション部１１ｂは、アジマス方向の（ｍ３＋Ｌ３）ピクセルの画像ずれを補正するため、アジマスリサンプリングを行う。
この結果、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）は以下のように記述される。
式（８０）において、ｓ（ｐ，ｑ）は、図１３に示すようなクラッタ成分１１１であり、ａ＋（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１２、ａ−（ｐ，ｑ）は図１３に示すような１次のアジマスアンビギュイティ成分１１３である。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
１次のアジマスアンビギュイティ成分は、アンテナのサイドローブパターンによって発生したものであるとして、２次以降の成分は無視できるほど小さいものとする。
ｎ１（ｐ，ｑ）はＳＡＲ画像（１）内の雑音成分、ｔ１（ｐ，ｑ）はＳＡＲ画像（１）内の移動目標成分である。
なお、レンジビンはＰ個、アジマスビンはＱ個存在するものとしている。
一方、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）は以下のように記述される。
式（８１）において、ｎ２（ｐ，ｑ）はＳＡＲ画像（２）内の雑音成分、ｔ２（ｐ，ｑ）はＳＡＲ画像（２）内の移動目標成分である。
また、ｐ番目のレンジビン、ｑ番目のアジマスビンにおけるＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）は以下のように記述される。
式（８２）において、ｎ３（ｐ，ｑ）はＳＡＲ画像（３）内の雑音成分、ｔ３（ｐ，ｑ）はＳＡＲ画像（３）内の移動目標成分である。
信号成分推定部１２の不要成分抑圧部２０は、レジストレーション部１１ａ，１１ｂにより位置合わせが行われたＳＡＲ画像（１）（２）（３）を受けると、そのＳＡＲ画像（１）（２）（３）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定して、そのクラッタ成分及びアジマスアンビギュイティ成分を抑圧し、クラッタ成分及びアジマスアンビギュイティ成分を抑圧したＳＡＲ画像（１）（２）（３）を出力する。
ＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）、レジストレーション後のＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）には、図７に示すような２種類のアジマスアンビギュイティ成分が混在している。
ＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）、レジストレーション後のＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）及びレジストレーション後のＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）には、図７に示すような２種類のアジマスアンビギュイティ成分が混在している。
不要成分抑圧部２０は、受信信号ｒ１（ｐ，ｑ），ｒ２（ｐ，ｑ），ｒ３（ｐ，ｑ）を用いて、なるべく移動目標成分が含まれないように、クラッタ成分ｓ（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定する。
以下、受信信号ｒ１（ｐ，ｑ），ｒ２（ｐ，ｑ），ｒ３（ｐ，ｑ）を用いて、クラッタ成分ｓ（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定する処理を説明する。
以下、受信信号ｒ１（ｐ，ｑ），ｒ２（ｐ，ｑ），ｒ３（ｐ，ｑ）を用いて、クラッタ成分ｓ（ｐ，ｑ）及びアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定する処理を説明する。
まず、受信信号ｒ１（ｐ，ｑ），ｒ２（ｐ，ｑ），ｒ３（ｐ，ｑ）の定式化を整理する。
受信信号ｒ１（ｐ，ｑ），ｒ２（ｐ，ｑ），ｒ３（ｐ，ｑ）の定式を行列とベクトルで表すと、以下のように与えられる。
受信信号ｒ１（ｐ，ｑ），ｒ２（ｐ，ｑ），ｒ３（ｐ，ｑ）の定式を行列とベクトルで表すと、以下のように与えられる。
この実施の形態３では、上記実施の形態１と同様の理由で、ウィナーフィルタを用いて、信号の推定を行う。
ウィナーフィルタを用いる推定方法を式で表すと、下記の式（８９）のようになる。
ただし、Ｅ（ａ＋／ｓ）はアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値であり、Ｅ（ａ−／ｓ）はアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値である。
また、Ｅ（ｎ／ｓ）は雑音成分ｎ（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値である。
ウィナーフィルタを用いる推定方法を式で表すと、下記の式（８９）のようになる。
また、Ｅ（ｎ／ｓ）は雑音成分ｎ（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値である。
この推定方式は、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）以外の信号の存在を仮定し、それらの信号の存在が、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）の推定の際に含まれる量を統計的に最小化した方式である。
したがって、ｒ１（ｐ，ｑ）−ｓ（ｐ，ｑ）−ａ＋（ｐ，ｑ）−ａ−（ｐ，ｑ）を計算した場合には、移動目標成分ｔ１（ｐ，ｑ）が残余する。
したがって、ｒ１（ｐ，ｑ）−ｓ（ｐ，ｑ）−ａ＋（ｐ，ｑ）−ａ−（ｐ，ｑ）を計算した場合には、移動目標成分ｔ１（ｐ，ｑ）が残余する。
不要成分抑圧部２０は、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定すると、それらの推定結果を用いて、ＳＡＲ画像（１）の受信信号ｒ１（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、ＳＡＲ画像（２）の受信信号ｒ２（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する。
また、ＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する。
最も単純な抑圧方法としては、以下の計算を行うものがある。
また、ＳＡＲ画像（３）の受信信号ｒ３（ｐ，ｑ）に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する。
最も単純な抑圧方法としては、以下の計算を行うものがある。
式（９２）の計算結果、式（９３）の計算結果又は式（９４）の計算結果を用いれば、移動目標を検出することが可能であるが、下記の式（９５）に示すように、３チャンネルの成分を合成すれば、より移動目標の検出性能を高めることができる。
なお、不要成分抑圧部２０により推定されたクラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）には推定誤差が含まれている。
特に、クラッタ成分の推定誤差は、クラッタ成分の電力が大きいことから、移動目標の検出性能劣化に大きく寄与すると考えられる。
そこで、不要成分抑圧部２０により推定されたクラッタ成分ｓ（ｐ，ｑ）を用いずにクラッタを除去し、残ったアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を抑圧する下記の処理が考えられる。
特に、クラッタ成分の推定誤差は、クラッタ成分の電力が大きいことから、移動目標の検出性能劣化に大きく寄与すると考えられる。
そこで、不要成分抑圧部２０により推定されたクラッタ成分ｓ（ｐ，ｑ）を用いずにクラッタを除去し、残ったアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を抑圧する下記の処理が考えられる。
上記の不要成分抑圧部２０によるクラッタ成分及びアジマスアンビギュイティ成分の抑圧は、振幅の次元で行っているが、下記に示すように、電力の次元で行ってから、振幅の次元に戻すことで、移動目標の検出性能を高めることも可能である。
目標検出部１５は、不要成分抑圧部２０によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像である抑圧後画像から移動目標を検出する。
移動目標の検出処理としては、閾値判定処理などが考えられるが、抑圧後画像を電力成分に変換してから閾値判定処理を行うようにしてもよい。また、雑音を低減するために、移動平均フィルタなどを施してから閾値判定処理を行うようにしてもよい。
移動目標の検出処理としては、閾値判定処理などが考えられるが、抑圧後画像を電力成分に変換してから閾値判定処理を行うようにしてもよい。また、雑音を低減するために、移動平均フィルタなどを施してから閾値判定処理を行うようにしてもよい。
また、ＳＡＲ画像（１）とＳＡＲ画像（２）、ＳＡＲ画像（１）とＳＡＲ画像（３）、ＳＡＲ画像（２）とＳＡＲ画像（３）を干渉させて、移動目標の速度推定を行うＡＴＩにも用いることが可能である。
具体的には、以下の偏角による位相マップを計算する。
式（１０２）〜（１０４）において、＊は共役を示す記号である。
具体的には、以下の偏角による位相マップを計算する。
以上で明らかなように、この実施の形態３によれば、図７に示すような２種類のアジマスアンビギュイティ成分が混在している場合でも、クラッタ成分とアジマスアンビギュイティ成分の両方を抑圧して、移動目標を検出することができる効果を奏する。
この実施の形態３では、上記の式（９０）（９１）に示すように、アジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値Ｅ（ａ＋／ｓ）と、アジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値Ｅ（ａ−／ｓ）と、雑音成分ｎ（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値Ｅ（ｎ／ｓ）とを用いて、ウィナーフィルタを生成し、そのウィナーフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、その推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、不要成分抑圧部２０は、期待値Ｅ（ａ＋／ｓ），Ｅ（ａ−／ｓ），Ｅ（ｎ／ｓ）を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比（ａ＋／ｓ）と、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比（ａ−／ｓ）と、その推定したクラッタ成分ｓ（ｐ，ｑ）と雑音成分ｎ（ｐ，ｑ）の電力比（ｎ／ｓ）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部２０は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
なお、クラッタ成分及びアジマスアンビギュイティ成分の収束判定としては、例えば、前回の推定結果と今回推定結果を比較して、それらの推定結果の差分が予め設定された閾値以内であれば、収束していると判定する方法が考えられる。
なお、クラッタ成分及びアジマスアンビギュイティ成分の収束判定としては、例えば、前回の推定結果と今回推定結果を比較して、それらの推定結果の差分が予め設定された閾値以内であれば、収束していると判定する方法が考えられる。
この実施の形態３では、不要成分抑圧部２０が、期待値Ｅ（ａ＋／ｓ），Ｅ（ａ−／ｓ），Ｅ（ｎ／ｓ）を用いて、ウィナーフィルタを生成するものを示したが、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、クラッタ成分の信号電力を用いて、ウィナーフィルタを生成するようにしてもよい。
即ち、上記の式（９０）（９１）における期待値Ｅ（ａ＋／ｓ），Ｅ（ａ−／ｓ），Ｅ（ｎ／ｓ）の代わりに、クラッタ成分の信号電力に反比例する値を代入して、ウィナーフィルタを生成するようにしてもよい。
この場合も、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、上記の式（９０）（９１）における期待値Ｅ（ａ＋／ｓ），Ｅ（ａ−／ｓ），Ｅ（ｎ／ｓ）の代わりに、クラッタ成分の信号電力に反比例する値を代入して、ウィナーフィルタを生成するようにしてもよい。
この場合も、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返し実施して、その推定精度を高めるようにしてもよい。
即ち、不要成分抑圧部２０は、クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定したのち、その推定したアジマスアンビギュイティ成分ａ＋（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比（ａ＋／ｓ）と、その推定したアジマスアンビギュイティ成分ａ−（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比（ａ−／ｓ）と、その推定したクラッタ成分ｓ（ｐ，ｑ）と雑音成分ｎ（ｐ，ｑ）の電力比（ｎ／ｓ）とを用いて、ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するようにする。
不要成分抑圧部２０は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
不要成分抑圧部２０は、アジマスアンビギュイティ成分等の推定処理とウィナーフィルタの生成処理を再帰的に繰り返す場合、例えば、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、その推定処理と生成処理を繰り返すようにしてもよいし、予め設定された回数だけ、その推定処理と生成処理を繰り返すようにしてもよい。
この実施の形態３では、不要成分抑圧部２０が、ウィナーフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、正則化最小二乗法を行うフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するようにしてもよい。 正則化最小二乗法を行うフィルタは、下記の式（１０５）のように表される。
式（１０５）において、Ｉは単位行列、λは任意の実数である。
ここでは、λが任意の実数であるものを示したが、λの代わりに、雑音成分ｎ（ｐ，ｑ）とクラッタ成分ｓ（ｐ，ｑ）の電力比の期待値Ｅ（ｎ／ｓ）を式（１０５）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
また、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、λの代わりに、クラッタ成分の信号電力を式（１０５）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
また、信号成分におけるピクセル単位の信号電力を計算し、クラッタ成分の電力が支配的である場合（クラッタ成分の信号電力の割合が予め設定された閾値以上である場合）、λの代わりに、クラッタ成分の信号電力を式（１０５）に代入して、正則化最小二乗法を行うフィルタを生成するようにしてもよい。
この実施の形態３では、不要成分抑圧部２０が、ウィナーフィルタ又は正則化最小二乗法を行うフィルタを用いて、クラッタ成分ｓ（ｐ，ｑ）とアジマスアンビギュイティ成分ａ＋（ｐ，ｑ），ａ−（ｐ，ｑ）を推定するものを示したが、これに限るものではなく、推定したい成分の他の信号成分の存在を仮定するフィルタや推定方式を適用することが可能である。
また、この実施の形態３では、３本の受信アンテナを用いる例を示したが、４本以上の受信アンテナを用いる場合にも適用可能である。また、複数本の送信アンテナを用いる場合や、複数本の送信アンテナ及び複数本の受信アンテナを用いる場合にも対応可能である。
また、この実施の形態３では、３本の受信アンテナを用いる例を示したが、４本以上の受信アンテナを用いる場合にも適用可能である。また、複数本の送信アンテナを用いる場合や、複数本の送信アンテナ及び複数本の受信アンテナを用いる場合にも対応可能である。
なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。
１ レジストレーション部（位置合わせ手段）、２ 信号成分推定部、３ 信号成分分別部（信号成分分別手段）、４ 不要成分抑圧部（不要成分抑圧手段）、５ 目標検出部（目標検出手段）、１１ａ，１１ｂ レジストレーション部（位置合わせ手段）、１２ 信号成分推定部、１３ 信号成分分別部（信号成分分別手段）、１４，２０ 不要成分抑圧部（不要成分抑圧手段）、１５ 目標検出部（目標検出手段）、１０１ レジストレーション部、１０２ 差分画像生成部、１１１ クラッタ成分、１１２，１１３ アジマスアンビギュイティ成分。
Claims (11)
- 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別手段と、
前記信号成分分別手段により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記信号成分分別手段により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分をウィナーフィルタを用いて推定し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを備え、
前記ウィナーフィルタは、前記信号成分分別手段により正の領域に分別された信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、前記正の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるとともに、前記信号成分分別手段により負の領域の信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、前記負の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されることを特徴とするレーダ信号処理装置。 - 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分とをウィナーフィルタを用いて推定し、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを備え、
前記ウィナーフィルタは、前記正の領域に分別された信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、前記正の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されるとともに、前記負の領域の信号成分に含まれているアジマスアンビギュイティ成分とクラッタ成分の電力比の期待値と、前記負の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて生成されることを特徴とするレーダ信号処理装置。 - 前記不要成分抑圧手段は、前記電力比の期待値を用いて生成したウィナーフィルタを用いて、前記正の領域及び前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定したのち、その推定したアジマスアンビギュイティ成分とクラッタ成分の電力比と、その推定したクラッタ成分と前記雑音成分の電力比とを用いて、前記ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、前記正の領域及び前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を再推定することを特徴とする請求項１又は請求項２記載のレーダ信号処理装置。
- 前記不要成分抑圧手段は、前記正の領域及び前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定すると、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束しているか否かを判定し、その推定したクラッタ成分及びアジマスアンビギュイティ成分が収束するまで、前記ウィナーフィルタを再生成しながら、その再生成したウィナーフィルタを用いて、前記正の領域及び前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定する処理を繰り返し実施することを特徴とする請求項３記載のレーダ信号処理装置。
- 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別手段と、
前記信号成分分別手段により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記信号成分分別手段により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを備え、
前記不要成分抑圧手段は、前記信号成分分別手段により正の領域に分別された信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上であり、前記信号成分分別手段により負の領域に分別された信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上である場合、前記正の領域に分別された信号成分に含まれているクラッタ成分の信号電力と前記負の領域に分別された信号成分に含まれているクラッタ成分の信号電力を用いて、ウィナーフィルタを生成し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分の推定は前記ウィナーフィルタを用いて実施することを特徴とするレーダ信号処理装置。 - 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定し、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを備え、
前記不要成分抑圧手段は、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の正の領域の信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上であり、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の負の領域の信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上である場合、前記正の領域の信号成分に含まれているクラッタ成分の信号電力と前記負の領域の信号成分に含まれているクラッタ成分の信号電力を用いて、ウィナーフィルタを生成し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分の推定を前記ウィナーフィルタを用いて実施することを特徴とするレーダ信号処理装置。 - 前記不要成分抑圧手段は、前記クラッタ成分の信号電力を用いて生成したウィナーフィルタを用いて、前記正の領域及び前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定したのち、その推定したアジマスアンビギュイティ成分とクラッタ成分の電力比を用いて、前記ウィナーフィルタを再生成し、その再生成したウィナーフィルタを用いて、前記正の領域及び前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を再推定することを特徴とする請求項５又は請求項６記載のレーダ信号処理装置。
- 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別手段と、
前記信号成分分別手段により正の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値と、前記信号成分分別手段により負の領域に分別された信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて、正則化最小二乗法を行うフィルタを生成し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を前記フィルタを用いて推定し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段と
を備えたレーダ信号処理装置。 - 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の正の領域の信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値と、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の負の領域の信号成分に含まれている雑音成分とクラッタ成分の電力比の期待値とを用いて、正則化最小二乗法を行うフィルタを生成し、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を前記フィルタを用いて推定することで前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を前記フィルタを用いて推定することで前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段と
を備えたレーダ信号処理装置。 - 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の信号成分を正の領域と負の領域に分別する信号成分分別手段と、
前記信号成分分別手段により正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記信号成分分別手段により負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを備え、
前記不要成分抑圧手段は、前記信号成分分別手段により正の領域に分別された信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上であり、前記信号成分分別手段により負の領域に分別された信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上である場合、前記正の領域に分別された信号成分に含まれているクラッタ成分の信号電力と前記信号成分分別手段により正の領域に分別された信号成分に含まれているクラッタ成分の信号電力を用いて、正則化最小二乗法を行うフィルタを生成し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分の推定を前記フィルタを用いて実施することを特徴とするレーダ信号処理装置。 - 観測領域が同一で観測時刻が異なる複数の合成開口レーダ画像の位置合わせを行う位置合わせ手段と、
前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分以外の成分を仮定して、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を推定し、前記正の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧するとともに、前記負の領域の信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分を抑圧する不要成分抑圧手段と、
前記不要成分抑圧手段によりクラッタ成分及びアジマスアンビギュイティ成分が抑圧された信号成分を有する画像から移動目標を検出する目標検出手段とを備え、
前記不要成分抑圧手段は、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の正の領域の信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上であり、前記位置合わせ手段により位置合わせが行われた複数の合成開口レーダ画像におけるドップラー周波数領域の負の領域の信号成分におけるピクセル単位の信号電力のうち、前記クラッタ成分の信号電力の割合が予め設定された閾値以上である場合、前記正の領域の信号成分に含まれているクラッタ成分の信号電力と前記負の領域の信号成分に含まれているクラッタ成分の信号電力を用いて、正則化最小二乗法を行うフィルタを生成し、前記正の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分と前記負の領域に分別された信号成分に含まれているクラッタ成分及びアジマスアンビギュイティ成分の推定を前記フィルタを用いて実施することを特徴とするレーダ信号処理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2014038708A JP6246024B2 (ja) | 2014-02-28 | 2014-02-28 | レーダ信号処理装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2014038708A JP6246024B2 (ja) | 2014-02-28 | 2014-02-28 | レーダ信号処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2015161662A JP2015161662A (ja) | 2015-09-07 |
| JP6246024B2 true JP6246024B2 (ja) | 2017-12-13 |
Family
ID=54184830
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2014038708A Active JP6246024B2 (ja) | 2014-02-28 | 2014-02-28 | レーダ信号処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6246024B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017022123A1 (ja) * | 2015-08-06 | 2017-02-09 | 三菱電機株式会社 | 画像処理装置および画像処理方法 |
| JP6177446B1 (ja) * | 2016-01-12 | 2017-08-09 | 三菱電機株式会社 | 信号処理装置、信号処理方法及びレーダ装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2972266B1 (fr) * | 2011-03-03 | 2013-04-26 | Thales Sa | Procede de detection et de caracterisation d'une cible mobile sur une image radar |
| JP5933295B2 (ja) * | 2012-03-05 | 2016-06-08 | 三菱電機株式会社 | 合成開口レーダ装置およびその移動目標検出方法 |
-
2014
- 2014-02-28 JP JP2014038708A patent/JP6246024B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2015161662A (ja) | 2015-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103869313B (zh) | 一种多通道合成孔径雷达距离向模糊的抑制方法及装置 | |
| CN101609151B (zh) | 一种基于单通道sar图像序列分解的动目标检测方法 | |
| CN110412559A (zh) | 分布式无人机mimo雷达的非相参融合目标检测方法 | |
| US20140285371A1 (en) | Bistatic inverse synthetic aperture radar imaging | |
| US20170299694A1 (en) | Antenna specification estimation device and radar device | |
| JP2019086469A (ja) | レーダ装置及び到来方向推定装置 | |
| EP3751309B1 (en) | Radar image processing device and radar image processing method | |
| CN110806577B (zh) | 合成孔径雷达的聚焦成像方法及装置、设备、存储介质 | |
| Liu et al. | A novel channel phase bias estimation method for spaceborne along-track multi-channel HRWS SAR in time-domain | |
| JP5811931B2 (ja) | 位相モノパルスレーダ装置 | |
| Shu et al. | Robust radial velocity estimation of moving targets based on adaptive data reconstruction and subspace projection algorithm | |
| US8760340B2 (en) | Processing radar return signals to detect targets | |
| CN101509975B (zh) | 一种基于双子孔径不同重叠度的运动目标检测方法 | |
| CN114966657A (zh) | 使用稀疏阵列处理的基于雷达的检测 | |
| Xu et al. | Joint Doppler and DOA estimation using (Ultra-) Wideband FMCW signals | |
| JP2010175457A (ja) | レーダ装置 | |
| JP5047002B2 (ja) | 波数推定装置 | |
| JP6246024B2 (ja) | レーダ信号処理装置 | |
| CN114265058B (zh) | Mimo雷达目标测角方法、装置、电子设备及存储介质 | |
| JP5116590B2 (ja) | 波数推定装置 | |
| JP2013044642A (ja) | パッシブレーダ装置 | |
| JP6287674B2 (ja) | 遅延時間推定装置及び測高装置 | |
| CN111880154B (zh) | 基于对称波数谱对消的复图像域运动目标检测方法 | |
| WO2023186386A1 (en) | Detecting and suppressing ambiguities in synthetic aperture radar data and images | |
| Nguyen et al. | Analysis of 2D CA-CFAR for DDMA FMCW MIMO Radar |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20161109 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170809 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170822 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20171003 |
|
| 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: 20171017 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20171114 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6246024 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |