WO2019189765A1 - レーダ画像処理装置及びレーダ画像処理方法 - Google Patents

Abstract

第１の電波受信地点から観測領域が撮像されている第１のレーダ画像における反射信号の抑圧対象範囲である第１の抑圧対象範囲に含まれている複数の画素と、第２の電波受信地点から観測領域が撮像されている第２のレーダ画像における反射信号の抑圧対象範囲である第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点に対する位相と第２の電波受信地点に対する位相との位相差を算出する位相差算出部（７１）と、それぞれの位相差から、第２の抑圧対象範囲に含まれている複数の画素の位相の回転量をそれぞれ算出する回転量算出部（３１）とを設け、差分算出部（３２）が、それぞれの回転量に基づいて、第２の抑圧対象範囲に含まれている複数の画素の位相を回転させ、第１の抑圧対象範囲に含まれている複数の画素と第２の抑圧対象範囲に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出するように、レーダ画像処理装置（１０）を構成した。

Description

レーダ画像処理装置及びレーダ画像処理方法

　この発明は、第１のレーダ画像に含まれている画素と第２のレーダ画像に含まれている位相回転後の画素との差分を算出するレーダ画像処理装置及びレーダ画像処理方法に関するものである。

　レーダ装置により得られるレーダ画像には、高さが高い建物などが散乱体として映っていることがある。
　レーダ装置を実装しているプラットフォームから散乱体の高い位置までの距離は、散乱体の高さの分だけ、プラットフォームから散乱体の低い位置までの距離よりも短くなる。
　プラットフォームから散乱体の高い位置までの距離が、散乱体の低い位置までの距離よりも短くなることで、散乱体の高い位置で反射された信号が、プラットフォーム側に倒れ込む現象であるレイオーバが発生する。
　レイオーバが発生することで、散乱体の高い位置で反射された信号は、倒れ込んだ先の別の反射信号と重なるため、レーダ画像において、複数の反射信号が１つの画素に混在してしまうことがある。

　以下の非特許文献１には、第１のレーダ画像に含まれている画素と第２のレーダ画像に含まれている画素との差分を算出するレーダ画像処理装置が開示されている。
　レーダ画像処理装置は、差分を算出することで、１つの画素に混在している複数の反射信号のうち、第１の電波受信地点に対する位相と、第２の電波受信地点に対する位相との位相差が零である反射信号については抑圧することができる。
　第１の電波受信地点は、第１のレーダ画像を撮像したときのプラットフォームの位置であり、第２の電波受信地点は、第２のレーダ画像を撮像したときのプラットフォームの位置である。

D. L. Bickel, "A null-steering viewpoint of interferometric SAR," IGARSS2000

　従来のレーダ画像処理装置は、１つの画素に混在している複数の反射信号のうち、第１の電波受信地点に対する位相と、第２の電波受信地点に対する位相との位相差が零である反射信号については抑圧することが可能である。
　しかし、１つの画素に混在している複数の反射信号のうち、抑圧可能な反射信号が散乱された位置と同じ高さの位置で散乱されている反射信号は、第１の電波受信地点に対する位相と、第２の電波受信地点に対する位相との位相差が零にならない。
　第１の電波受信地点に対する位相と、第２の電波受信地点に対する位相との位相差が零でない反射信号については、抑圧することができないという課題があった。

　この発明は上記のような課題を解決するためになされたもので、第１の電波受信地点に対する位相と、第２の電波受信地点に対する位相との位相差が零でない反射信号についても抑圧することができるレーダ画像処理装置及びレーダ画像処理方法を得ることを目的とする。

　この発明に係るレーダ画像処理装置は、第１の電波受信地点から観測領域が撮像されている第１のレーダ画像における反射信号の抑圧対象範囲である第１の抑圧対象範囲に含まれている複数の画素と、第２の電波受信地点から観測領域が撮像されている第２のレーダ画像における反射信号の抑圧対象範囲である第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点に対する位相と第２の電波受信地点に対する位相との位相差を算出する位相差算出部と、位相差算出部により算出されたそれぞれの位相差から、第２の抑圧対象範囲に含まれている複数の画素の位相の回転量をそれぞれ算出する回転量算出部とを設け、差分算出部が、回転量算出部により算出されたそれぞれの回転量に基づいて、第２の抑圧対象範囲に含まれている複数の画素の位相を回転させ、第１の抑圧対象範囲に含まれている複数の画素と第２の抑圧対象範囲に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出するようにしたものである。

　この発明によれば、第１の電波受信地点に対する位相と、第２の電波受信地点に対する位相との位相差が零でない反射信号についても抑圧することができる。

実施の形態１によるレーダ画像処理装置１０を示す構成図である。 実施の形態１によるレーダ画像処理装置１０の位相処理部１２を示す構成図である。 実施の形態１によるレーダ画像処理装置１０の画像処理部１３を示す構成図である。 位相処理部１２及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。 位相処理部１２及び画像処理部１３が、ソフトウェア又はファームウェアなどで実現される場合のコンピュータのハードウェア構成図である。 位相処理部１２の処理内容を示すフローチャートである。 傾き面５１、平行面５２及び撮像パラメータを示す説明図である。 スラントレンジ方向の画素のスペーシングΔｓｌと、レーダ画像（第１のレーダ画像、第２のレーダ画像）の範囲Ｓｗと、レーダ画像の中心位置に対応するスラントレンジ方向の位置から観測領域までの距離ｓｌとの関係を示す説明図である。 画像処理部１３の処理内容を示すフローチャートである。 位相回転部３３によって、第２のレーダ画像に含まれている画素の位相が回転されない場合の、１つの画素に混在している反射信号の抑圧を示す説明図である。 位相回転部３３によって、第２のレーダ画像に含まれている画素の位相が回転された場合の、１つの画素に混在している反射信号の抑圧を示す説明図である。 実施の形態２によるレーダ画像処理装置１０の画像処理部１３を示す構成図である。 位相処理部１２及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。 レーダ画像群２に含まれているレーダ画像が２つだけである場合の、１つの画素に混在している複数の反射信号を示す説明図である。 レーダ画像群２に含まれているレーダ画像が２つ以上である場合の、１つの画素に混在している複数の反射信号を示す説明図である。 実施の形態３によるレーダ画像処理装置１０の画像処理部１３を示す構成図である。 位相処理部１２及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。 実施の形態４によるレーダ画像処理装置１０を示す構成図である。 実施の形態４によるレーダ画像処理装置１０の位相処理部１５を示す構成図である。 位相処理部１５及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。 図２１Ａは、第１の抑圧対象範囲を示す説明図、図２１Ｂは、第２の抑圧対象範囲を示す説明図である。 図２２Ａは、第１の抑圧対象範囲を示す説明図、図２２Ｂは、第２の抑圧対象範囲を示す説明図である。

　以下、この発明をより詳細に説明するために、この発明を実施するための形態について、添付の図面に従って説明する。

実施の形態１．
　図１は、実施の形態１によるレーダ画像処理装置１０を示す構成図である。
　図１において、レーダ１は、合成開口レーダ（ＳＡＲ：Ｓｙｎｔｈｅｔｉｃ　Ａｐｅｒｔｕｒｅ　Ｒａｄａｒ）又は実開口レーダなどが該当し、地球等を観測するプラットフォームに実装されている。レーダ１は、レーダ画像を撮像するとともに、レーダ画像の撮像時のパラメータを取得する。プラットフォームとしては、人工衛星又は航空機などが該当する。
　レーダ１は、或る電波受信地点から、観測領域を撮像した後、プラットフォームが、上記の電波受信地点と近い電波受信地点にあるときに、上記の観測領域を再度撮像する。
　リピートパス撮像の場合、プラットフォームが人工衛星であれば、レーダ１が、或る電波受信地点から、観測領域を撮像した後、プラットフォームが地球を周回して、上記の電波受信地点と近い電波受信地点に戻ったときに、レーダ１が、同じ観測領域を再度撮像してレーダ画像を取得する。プラットフォームが航空機であれば、同じ軌道を繰り返し通過するように操縦されており、プラットフォームが、概ね同じ電波受信地点にあるときに、レーダ１が、同じ観測領域を撮像してレーダ画像を取得する。
　シングルパス撮像の場合、複数のレーダ１が同じプラットフォームに実装されており、複数のレーダ１が、或る電波受信地点から、同じ観測領域を撮像してレーダ画像を取得する。このとき、複数のレーダ１は、同じプラットフォームの異なる位置に設置されている。
　また、波長などの撮像パラメータの等しい複数のレーダ１のそれぞれが、異なるプラットフォームに実装されており、複数のレーダ１が、或る電波受信地点から、同じ観測領域を撮像してレーダ画像を取得する。
　したがって、レーダ１は、互いに異なるそれぞれの電波受信地点から、同じ観測領域を２回撮像することで、２つレーダ画像として、第１のレーダ画像と第２のレーダ画像とを取得する。

　以下、第１のレーダ画像を撮像したときのプラットフォームの位置を第１の電波受信地点、第２のレーダ画像を撮像したときのプラットフォームの位置を第２の電波受信地点とする。
　第１のレーダ画像の解像度と第２のレーダ画像の解像度とは、同じ解像度である。したがって、第１のレーダ画像に含まれている複数の画素の画素位置と、第２のレーダ画像に含まれている複数の画素の画素位置とは、同じ（ｐｉｘｅｌ，ｌｉｎｅ）で表される。
　ｐｉｘｅｌは、第１のレーダ画像及び第２のレーダ画像のそれぞれにおけるスラントレンジ方向の画素の位置を示す変数、ｌｉｎｅは、第１のレーダ画像及び第２のレーダ画像のそれぞれにおけるアジマス方向の画素の位置を示す変数である。
　レーダ１は、第１のレーダ画像及び第２のレーダ画像を含んでいるレーダ画像群２をレーダ画像処理装置１０に送信する。
　レーダ１は、第１のレーダ画像に対応する第１の撮像パラメータと、第２のレーダ画像に対応する第２の撮像パラメータとを含んでいる撮像パラメータ群３をレーダ画像処理装置１０に送信する。

　レーダ画像群２は、第１のレーダ画像及び第２のレーダ画像を含んでいる画像群である。
　第１のレーダ画像の撮像及び第２のレーダ画像の撮像に用いられるそれぞれの偏波の種類は限定されないため、第１のレーダ画像及び第２のレーダ画像のそれぞれは、単偏波レーダ画像、二偏波レーダ画像又は四偏波レーダ画像のいずれであってもよい。
　第１のレーダ画像及び第２のレーダ画像のそれぞれは、レーダ１から電波が放射されたのち、観測領域に反射されて、レーダ１により受信された上記電波の強度分布を示すレーダ画像である。
　第１のレーダ画像に含まれている複数の画素及び第２のレーダ画像に含まれている複数の画素は、それぞれ複素数の画素値を有している。
　複素数の画素値は、レーダ１と観測領域に存在している散乱体との距離を示す情報のほか、レーダ１から放射された電波が散乱体に反射された際に発生する位相シフトを示す情報を含んでいる。以降特に断りが無い場合、「画素値」とは複素数の値を持つとする。

　撮像パラメータ群３は、第１の撮像パラメータ及び第２の撮像パラメータを含んでいるパラメータ群である。
　第１の撮像パラメータは、レーダ１によって第１のレーダ画像が撮像されたときのプラットフォームにおける軌道の位置情報と、センサ情報とを含んでいる。
　第２の撮像パラメータは、レーダ１によって第２のレーダ画像が撮像されたときのプラットフォームにおける軌道の位置情報と、センサ情報とを含んでいる。
　軌道の位置情報は、レーダ１によって、第１のレーダ画像又は第２のレーダ画像が撮影されたときのプラットフォームの緯度、経度及び高度を示す情報である。したがって、軌道の位置情報は、第１の電波受信地点又は第２の電波受信地点を示す情報として用いられる。
　センサ情報は、第１のレーダ画像又は第２のレーダ画像が撮像されたときのレーダ１のオフナディア角θを示す情報、レーダ１からの放射電波の波長λを示す情報及びレーダ１から観測領域までの距離の平均値Ｒを示す情報を含んでいる。

　レーダ画像処理装置１０は、レーダ画像取得部１１、位相処理部１２及び画像処理部１３を備えている。
　レーダ画像取得部１１は、レーダ１から送信されたレーダ画像群２及び撮像パラメータ群３のそれぞれを取得する。
　レーダ画像取得部１１は、レーダ画像群２を画像処理部１３に出力し、撮像パラメータ群３を位相処理部１２に出力する。
　位相処理部１２は、レーダ画像取得部１１から出力された撮像パラメータ群３と、グランドレンジ方向に対する２次元の傾き面５１（図７を参照）の傾き角度αとを取得する。
　また、位相処理部１２は、傾き面５１と、傾き面５１と平行な面である平行面５２（図７を参照）との距離ｚ_０を取得する。
　傾き面５１及び平行面５２の詳細については後述する。

　位相処理部１２は、第１の撮像パラメータと、第２の撮像パラメータと、傾き角度αとを用いて、傾き面５１におけるｘ軸（第１の軸）方向の位相変化成分φ（ｘ）を算出する処理を実施する。
　位相処理部１２は、第１の撮像パラメータと、第２の撮像パラメータと、傾き角度αと、距離ｚ_０とを用いて、平行面５２における傾き面５１に対する位相ρ（ｚ_０）を算出する処理を実施する。
　位相処理部１２は、第１のレーダ画像に含まれている複数の画素と、第２のレーダ画像に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点に対する位相と第２の電波受信地点に対する位相との位相差Δφ（ｘ，ｚ_０）を算出する処理を実施する。

　画像処理部１３は、レーダ画像取得部１１から出力されたレーダ画像群２と、位相処理部１２から出力されたそれぞれの位相差Δφ（ｘ，ｚ_０）とを取得する。
　画像処理部１３は、位相処理部１２より出力されたそれぞれの位相差Δφ（ｘ，ｚ_０）から、第２のレーダ画像に含まれている複数の画素の位相の回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］をそれぞれ算出する処理を実施する。
　画像処理部１３は、算出したそれぞれの回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］に基づいて、第２のレーダ画像に含まれている複数の画素の位相を回転させる処理を実施する。
　画像処理部１３は、第１のレーダ画像に含まれている複数の画素と、第２のレーダ画像に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出する処理を実施する。

　図２は、実施の形態１によるレーダ画像処理装置１０の位相処理部１２を示す構成図である。
　図３は、実施の形態１によるレーダ画像処理装置１０の画像処理部１３を示す構成図である。
　図４は、位相処理部１２及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。

　図２において、位相変化成分算出部２１は、例えば、図４に示す位相変化成分算出回路４１で実現される。
　位相変化成分算出部２１は、レーダ画像取得部１１から出力された撮像パラメータ群３と、傾き角度αとを取得する。
　位相変化成分算出部２１は、第１の撮像パラメータと、第２の撮像パラメータと、傾き角度αとを用いて、傾き面５１におけるｘ軸方向の位相変化成分φ（ｘ）を算出する処理を実施する。
　位相変化成分算出部２１は、ｘ軸方向の位相変化成分φ（ｘ）を位相差算出部２３に出力する。

　位相算出部２２は、例えば、図４に示す位相算出回路４２で実現される。
　位相算出部２２は、レーダ画像取得部１１から出力された撮像パラメータ群３と、傾き角度αと、距離ｚ_０とを取得する。
　位相算出部２２は、第１の撮像パラメータと、第２の撮像パラメータと、傾き角度αと、距離ｚ_０とを用いて、平行面５２における傾き面５１に対する位相ρ（ｚ_０）を算出する処理を実施する。
　位相算出部２２は、位相ρ（ｚ_０）を位相差算出部２３に出力する。

　位相差算出部２３は、例えば、図４に示す位相差算出回路４３で実現される。
　位相差算出部２３は、位相変化成分φ（ｘ）と位相ρ（ｚ_０）とから、第１のレーダ画像に含まれている複数の画素と、第２のレーダ画像に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、位相差Δφ（ｘ，ｚ_０）を算出する処理を実施する。
　位相差Δφ（ｘ，ｚ_０）は、それぞれの反射信号において、当該反射信号における第１の電波受信地点に対する位相と、当該反射信号における第２の電波受信地点に対する位相との位相差である。
　位相差算出部２３は、それぞれの位相差Δφ（ｘ，ｚ_０）を画像処理部１３に出力する。

　図３において、回転量算出部３１は、例えば、図４に示す回転量算出回路４４で実現される。
　回転量算出部３１は、位相差算出部２３より出力されたそれぞれの位相差Δφ（ｘ，ｚ_０）から、第２のレーダ画像に含まれている複数の画素の位相の回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］をそれぞれ算出する処理を実施する。
　回転量算出部３１は、それぞれの回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］を位相回転部３３に出力する。

　差分算出部３２は、位相回転部３３及び差分算出処理部３４を備えている。
　位相回転部３３は、例えば、図４に示す位相回転回路４５で実現される。
　位相回転部３３は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、第２のレーダ画像を取得する。
　位相回転部３３は、回転量算出部３１から出力されたそれぞれの回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］に基づいて、第２のレーダ画像に含まれている複数の画素の位相を回転させる処理を実施する。
　位相回転部３３は、位相回転後の複数の画素を含む第２のレーダ画像を差分算出処理部３４に出力する。

　差分算出処理部３４は、例えば、図４に示す差分算出処理回路４６で実現される。
　差分算出処理部３４は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、第１のレーダ画像を取得し、位相回転部３３から出力された第２のレーダ画像を取得する。
　差分算出処理部３４は、第１のレーダ画像に含まれている複数の画素と、第２のレーダ画像に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）を算出する処理を実施する。
　差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）は、不要な散乱体の反射信号が抑圧された抑圧画像の画素である。
　差分算出処理部３４は、それぞれの差分Δｓ（ｐｉｘｅｌ，ｌｉｎｅ）を含む抑圧画像を外部に出力する。

　図２では、位相処理部１２の構成要素である位相変化成分算出部２１、位相算出部２２及び位相差算出部２３のそれぞれが、図４に示すような専用のハードウェアで実現されるものを想定している。
　また、図３では、画像処理部１３の構成要素である回転量算出部３１、位相回転部３３及び差分算出処理部３４のそれぞれが、図４に示すような専用のハードウェアで実現されるものを想定している。
　即ち、位相処理部１２及び画像処理部１３が、位相変化成分算出回路４１、位相算出回路４２、位相差算出回路４３、回転量算出回路４４、位相回転回路４５及び差分算出処理回路４６で実現されるものを想定している。

　ここで、位相変化成分算出回路４１、位相算出回路４２、位相差算出回路４３、回転量算出回路４４、位相回転回路４５及び差分算出処理回路４６のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ　Ｇａｔｅ　Ａｒｒａｙ）、又は、これらを組み合わせたものが該当する。

　位相処理部１２の構成要素及び画像処理部１３の構成要素は、専用のハードウェアで実現されるものに限るものではない。位相処理部１２及び画像処理部１３が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせで実現されるものであってもよい。
　ソフトウェア又はファームウェアは、プログラムとして、コンピュータのメモリに格納される。コンピュータは、プログラムを実行するハードウェアを意味し、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、プロセッサ、あるいは、ＤＳＰ（Ｄｉｇｉｔａｌ　Ｓｉｇｎａｌ　Ｐｒｏｃｅｓｓｏｒ）が該当する。

　図５は、位相処理部１２及び画像処理部１３が、ソフトウェア又はファームウェアなどで実現される場合のコンピュータのハードウェア構成図である。
　位相処理部１２が、ソフトウェア又はファームウェアなどで実現される場合、位相変化成分算出部２１、位相算出部２２及び位相差算出部２３の処理手順をコンピュータに実行させるためのプログラムがメモリ６１に格納される。
　また、画像処理部１３が、ソフトウェア又はファームウェアなどで実現される場合、回転量算出部３１、位相回転部３３及び差分算出処理部３４の処理手順をコンピュータに実行させるためのプログラムがメモリ６１に格納される。
　そして、コンピュータのプロセッサ６２がメモリ６１に格納されているプログラムを実行する。

　また、図４では、位相処理部１２の構成要素及び画像処理部１３の構成要素のそれぞれが専用のハードウェアで実現される例を示し、図５では、位相処理部１２及び画像処理部１３が、ソフトウェア又はファームウェアなどで実現される例を示している。
　しかし、これに限るものではなく、例えば、位相処理部１２における一部の構成要素及び画像処理部１３における一部の構成要素が専用のハードウェアで実現され、残りの構成要素がソフトウェア又はファームウェアなどで実現されるものであってもよい。

　次に、図１に示すレーダ画像処理装置１０の動作について説明する。
　レーダ１は、第１のレーダ画像及び第２のレーダ画像を含んでいるレーダ画像群２と、第１の撮像パラメータ及び第２の撮像パラメータを含んでいる撮像パラメータ群３とをレーダ画像処理装置１０に送信する。
　レーダ画像取得部１１は、レーダ１から送信されたレーダ画像群２及び撮像パラメータ群３のそれぞれを取得する。
　レーダ画像取得部１１は、レーダ画像群２を画像処理部１３に出力し、撮像パラメータ群３を位相処理部１２に出力する。

　レーダ画像（第１のレーダ画像、第２のレーダ画像）に含まれている画素の画素値は、複素数であり、以下の式（１）のように表される。

　式（１）において、Ａｖ（ｐｉｘｅｌ，ｌｉｎｅ）は、画素位置が（ｐｉｘｅｌ，ｌｉｎｅ）の画素の振幅である。
　Ψ（ｐｉｘｅｌ，ｌｉｎｅ）は、画素位置が（ｐｉｘｅｌ，ｌｉｎｅ）の画素の位相（偏角）である。
　ｊは、虚数単位を示す記号である。

　位相処理部１２は、位相差Δφ（ｘ，ｚ_０）を算出する処理を実施する。
　図６は、位相処理部１２の処理内容を示すフローチャートである。
　以下、図６を参照しながら、位相処理部１２の処理内容を具体的に説明する。

　位相変化成分算出部２１は、レーダ画像取得部１１から出力された撮像パラメータ群３と、傾き角度αとを取得する（図６のステップＳＴ１）。
　位相算出部２２は、レーダ画像取得部１１から出力された撮像パラメータ群３と、傾き角度αと、距離ｚ_０とを取得する（図６のステップＳＴ２）。
　傾き角度αは、例えば、ユーザによって事前に設定されるパラメータであり、図７のように表される。
　距離ｚ_０は、例えば、ユーザによって事前に設定されるパラメータであり、図７のように表される。
　傾き角度α及び距離ｚ_０のそれぞれは、例えば、ユーザによる手動の操作で位相算出部２２に与えられるようにしてもよいし、図示せぬ外部機器によって位相算出部２２に与えられるようにしてもよい。

　図７は、傾き面５１、平行面５２及び撮像パラメータを示す説明図である。
　図７において、傾き面５１は、第１のレーダ画像及び第２のレーダ画像のそれぞれに含まれている共通の２次元面である。
　傾き面５１の第１の軸であるｘ軸の方向は、グランドレンジ方向から傾き角度αだけ傾いている方向であり、傾き面５１の第２の軸の方向は、アジマス方向（図７において、紙面奥行き方向）である。
　平行面５２は、傾き面５１と平行な面であり、傾き面５１との距離がｚ_０である。
　傾き面５１が、例えば、水平な地面に垂直に建てられている建物の水平屋根である場合、傾き角度αとして、０度が設定される。
　傾き面５１が、例えば、水平な地面に垂直に建てられている建物の壁面である場合、傾き角度αとして、９０度が設定される。

　Ｐ_１は、第１の電波受信地点を示し、Ｐ_２は、第２の電波受信地点を示している。
　第１の電波受信地点Ｐ_１は、第１のレーダ画像が撮像されたときのプラットフォームの軌道の中心位置であり、第２の電波受信地点Ｐ_２は、第２のレーダ画像が撮像されたときのプラットフォームの軌道の中心位置である。
　Ｂ_１，２は、第１の電波受信地点Ｐ_１と第２の電波受信地点Ｐ_２との距離のうち、レーダ１から放射される電波の方向（以下、「スラントレンジ方向」と称する。）に対して垂直な方向の距離成分である。
　θは、オフナディア角であり、オフナディア角は、プラットフォームの鉛直直下方向と、スラントレンジ方向とのなす角である。
　Ｒは、第１の電波受信地点Ｐ_１及び第２の電波受信地点Ｐ_２のそれぞれと、観測領域との距離の平均値である。
　距離成分Ｂ_１，２、オフナディア角θ及び距離の平均値Ｒは、撮像パラメータに含まれている情報である。
　Ｓｗは、観測対象を撮像している第１のレーダ画像の範囲及び第２のレーダ画像の範囲である。
　第１のレーダ画像の範囲Ｓｗと第２のレーダ画像の範囲Ｓｗとは、同じ範囲である。

　ここでは、第１の電波受信地点Ｐ_１及び第２の電波受信地点Ｐ_２のそれぞれと観測領域との距離が長いため、位相変化成分算出部２１は、オフナディア角θ及び距離の平均値Ｒのそれぞれが変化しないものとして取り扱う。
　即ち、第１の撮像パラメータに含まれているオフナディア角θと、第２の撮像パラメータに含まれているオフナディア角θとが同じ値である。
　また、第１の撮像パラメータに含まれている距離の平均値Ｒと、第２の撮像パラメータに含まれている距離の平均値Ｒとが同じ値である。
　また、第１のレーダ画像に含まれている複数の画素のうち、画素位置が（ｐｉｘｅｌ，ｌｉｎｅ）の画素と、第２のレーダ画像に含まれている複数の画素のうち、画素位置が（ｐｉｘｅｌ，ｌｉｎｅ）の画素とは、同じ画素位置の画素である。

　図８は、スラントレンジ方向の画素のスペーシングΔｓｌと、レンジ画像（第１のレーダ画像、第２のレーダ画像）の範囲Ｓｗと、レーダ画像の中心位置に対応するスラントレンジ方向の位置から観測領域までの距離ｓｌとの関係を示す説明図である。
　図８において、レーダ画像のニアレンジに対応するスラントレンジ方向の位置から、レーダ画像の中心位置に対応するスラントレンジ方向の位置までの距離は、（Ｓｗ／２）・ｓｉｎθである。
　したがって、距離ｓｌは、以下の式（２）のように表される。

　スペーシングΔｓｌ及びレーダ画像の範囲Ｓｗのそれぞれは、撮像パラメータに含まれている情報である。

　また、レーダ画像の中心位置から、傾き面５１におけるｘ軸方向の位置ｘと、距離ｓｌとの関係は、以下の式（３）のように表される。

　式（２）と式（３）より、以下の式（４）が成立する。

　位相変化成分算出部２１は、レーダ画像におけるスラントレンジ方向の画素の位置ｐｉｘｅｌを式（４）に代入することで、位置ｐｉｘｅｌに対応する傾き面５１での位置ｘを算出する。
　式（４）に代入する位置ｐｉｘｅｌにある画素には、複数の散乱体の反射信号が混在している。
　式（４）に代入する位置ｐｉｘｅｌは、例えば、ユーザによる手動の操作で位相変化成分算出部２１に与えられるようにしてもよいし、図示せぬ外部機器によって位相変化成分算出部２１に与えられるようにしてもよい。

　位相変化成分算出部２１は、距離成分Ｂ_１，２と、オフナディア角θと、距離の平均値Ｒと、放射電波の波長λと、傾き角度αと、観測パスパラメータｐとを用いて、傾き面５１において、ｘ軸方向の位置ｘでの位相変化成分φ（ｘ）を算出する（図６のステップＳＴ３）。
　観測パスパラメータｐは、レーダ画像撮像時の観測パスが、リピートパスであるのか、シングルパスであるのかを区別するパラメータであり、リピートパスの場合、ｐ＝２、シングルパスの場合、ｐ＝１となる。観測パスパラメータｐは、例えば、ユーザによる手動の操作によって、位相変化成分算出部２１及び位相算出部２２に与えられるようにしてもよいし、図示せぬ外部機器によって、位相変化成分算出部２１及び位相算出部２２に与えられるようにしてもよい。
　以下の式（５）は、位相変化成分算出部２１が用いる位相変化成分φ（ｘ）の算出式である。

　位相変化成分算出部２１は、ｘ軸方向の位相変化成分φ（ｘ）を位相差算出部２３に出力する。

　位相算出部２２は、距離成分Ｂ_１，２と、オフナディア角θと、距離の平均値Ｒと、放射電波の波長λと、傾き角度αと、距離ｚ_０と、観測パスパラメータｐとを用いて、平行面５２における傾き面５１に対する位相ρ（ｚ_０）を算出する（図６のステップＳＴ４）。
　以下の式（６）は、位相算出部２２が用いる位相ρ（ｚ_０）の算出式である。

　位相算出部２２は、位相ρ（ｚ_０）を位相差算出部２３に出力する。

　位相差算出部２３は、位相変化成分φ（ｘ）と位相ρ（ｚ_０）とを用いて、第１のレーダ画像に含まれている複数の画素と、第２のレーダ画像に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、位相差Δφ（ｘ，ｚ_０）を算出する（図６のステップＳＴ５）。
　位相差Δφ（ｘ，ｚ_０）は、それぞれの反射信号において、当該反射信号における第１の電波受信地点Ｐ_１に対する位相と、当該反射信号における第２の電波受信地点Ｐ_２に対する位相との位相差である。
　以下の式（７）は、位相差算出部２３が用いる位相差Δφ（ｘ，ｚ_０）の算出式である。

　位相差算出部２３は、それぞれの位相差Δφ（ｘ，ｚ_０）を画像処理部１３に出力する。

　画像処理部１３は、抑圧画像を得る処理を実施する。
　図９は、画像処理部１３の処理内容を示すフローチャートである。
　以下、図９を参照しながら、画像処理部１３の処理内容を具体的に説明する。

　回転量算出部３１は、位相差算出部２３から出力されたそれぞれの位相差Δφ（ｘ，ｚ_０）を取得する。
　回転量算出部３１は、それぞれの位相差Δφ（ｘ，ｚ_０）から、第２のレーダ画像に含まれている複数の画素の位相の回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］をそれぞれ算出する（図９のステップＳＴ１１）。
　回転量算出部３１は、それぞれの回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］を位相回転部３３に出力する。

　位相回転部３３は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、第２のレーダ画像を取得する。
　位相回転部３３は、回転量算出部３１から出力されたそれぞれの回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］に基づいて、第２のレーダ画像に含まれている複数の画素の位相を回転させる処理を実施する（図９のステップＳＴ１２）。
　以下の式（８）は、位相回転部３３による位相の回転処理を示す式である。

　式（８）において、Ｓ_２（ｐｉｘｅｌ，ｌｉｎｅ）は、レーダ画像取得部１１から出力された第２のレーダ画像に含まれている画素の画素値、Ｓ_２’（ｐｉｘｅｌ，ｌｉｎｅ）は、位相回転部３３によって、画素の位相が回転された第２のレーダ画像に含まれている画素の画素値である。
　位相回転部３３は、位相回転後の複数の画素を含む第２のレーダ画像を差分算出処理部３４に出力する。

　差分算出処理部３４は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、第１のレーダ画像を取得し、位相回転部３３から出力された位相回転後の複数の画素を含む第２のレーダ画像を取得する。
　差分算出処理部３４は、第１のレーダ画像に含まれている複数の画素と、第２のレーダ画像に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）を算出する（図９のステップＳＴ１３）。
　以下の式（９）は、差分算出処理部３４が用いる差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）の算出式である。

　式（９）において、Ｓ_１（ｐｉｘｅｌ，ｌｉｎｅ）は、第１のレーダ画像に含まれている画素の画素値である。
　差分算出処理部３４は、それぞれの差分Δｓ（ｐｉｘｅｌ，ｌｉｎｅ）を含む抑圧画像を外部に出力する。

　ここで、図１０は、位相回転部３３によって、第２のレーダ画像に含まれている画素の位相が回転されない場合の、１つの画素に混在している反射信号の抑圧を示す説明図である。
　図１０において、“１”が付されている反射信号については、当該反射信号を散乱している散乱体から第１の電波受信地点Ｐ_１までの距離と、当該反射信号を散乱している散乱体から第２の電波受信地点Ｐ_２までの距離とが等しい。したがって、“１”が付されている反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零である。
　よって、“１”が付されている反射信号についての差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）が零になるため、“１”が付されている反射信号は、抑圧される。

　“２”が付されている反射信号については、当該反射信号を散乱している散乱体から第１の電波受信地点Ｐ_１までの距離と、当該反射信号を散乱している散乱体から第２の電波受信地点Ｐ_２までの距離とが等しくない。したがって、“２”が付されている反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零以外である。
　よって、“２”が付されている反射信号についての差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）が零以外になるため、“２”が付されている反射信号は、抑圧されない。

　“３”が付されている反射信号についても、当該反射信号を散乱している散乱体から第１の電波受信地点Ｐ_１までの距離と、当該反射信号を散乱している散乱体から第２の電波受信地点Ｐ_２までの距離とが等しくない。したがって、“３”が付されている反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零以外である。
　よって、“３”が付されている反射信号についての差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）が零以外になるため、“３”が付されている反射信号は、抑圧されない。

　図１１は、位相回転部３３によって、第２のレーダ画像に含まれている画素の位相が回転された場合の、１つの画素に混在している反射信号の抑圧を示す説明図である。
　 “１”が付されている反射信号については、図１０に示すように、当該反射信号を散乱している散乱体から第１の電波受信地点Ｐ_１までの距離と、当該反射信号を散乱している散乱体から第２の電波受信地点Ｐ_２までの距離とが等しい。したがって、“１”が付されている反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零であり、回転量算出部３１によって算出される位相の回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］は零である。
　“１”が付されている反射信号については、位相の回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］が零であるであるため、図１０及び図１１に示すように、位相回転部３３によって、回転されない。したがって、“１”が付されている反射信号については、位相差Δφ（ｘ，ｚ_０）が零のままであるため、差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）が零になり、“１”が付されている反射信号は、抑圧される。

　“２”が付されている反射信号については、図１０に示すように、当該反射信号を散乱している散乱体から第１の電波受信地点Ｐ_１までの距離と、当該反射信号を散乱している散乱体から第２の電波受信地点Ｐ_２までの距離とが等しくない。したがって、“２”が付されている反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零以外である。よって、回転量算出部３１によって算出される位相の回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］は零以外である。
　“２”が付されている位相回転後の反射信号は、図１１に示すように、位相回転部３３によって、回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］だけ回転されても、第１の電波受信地点Ｐ_１までの距離と、第２の電波受信地点Ｐ_２までの距離とが等しくない。したがって、“２”が付されている位相回転後の反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零以外である。
　よって、“２”が付されている位相回転後の反射信号についての差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）が零以外になるため、“２”が付されている位相回転後の反射信号は、抑圧されない。

　“３”が付されている反射信号については、図１０に示すように、当該反射信号を散乱している散乱体から第１の電波受信地点Ｐ_１までの距離と、当該反射信号を散乱している散乱体から第２の電波受信地点Ｐ_２までの距離とが等しくない。したがって、“３”が付されている反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零以外である。よって、回転量算出部３１によって算出される位相の回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］は零以外である。
　“３”が付されている反射信号については、図１１に示すように、位相回転部３３によって、回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］だけ回転されることで、第１の電波受信地点Ｐ_１までの距離と、第２の電波受信地点Ｐ_２までの距離とが等しくなっている。したがって、“３”が付されている位相回転後の反射信号については、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ（ｘ，ｚ_０）が零である。
　よって、“３”が付されている位相回転後の反射信号についての差分ΔＳ（ｐｉｘｅｌ，ｌｉｎｅ）が零になるため、“３”が付されている位相回転後の反射信号は、抑圧される。

　以上の実施の形態１は、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波発射地点に対する位相と第２の電波発射地点に対する位相との位相差を算出する位相差算出部２３と、それぞれの位相差から、第２のレーダ画像に含まれている複数の画素の位相の回転量をそれぞれ算出する回転量算出部３１とを設け、差分算出部３２が、それぞれの回転量に基づいて、第２のレーダ画像に含まれている複数の画素の位相を回転させ、第１のレーダ画像に含まれている複数の画素と第２のレーダ画像に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出するように、レーダ画像処理装置１０を構成した。したがって、レーダ画像処理装置１０は、第１の電波発射地点に対する位相と、第２の電波発射地点に対する位相との位相差が零でない反射信号についても抑圧することができる。

実施の形態２．
　実施の形態１のレーダ画像処理装置１０は、第１のレーダ画像及び第２のレーダ画像を含んでいるレーダ画像群２を取得して、抑圧画像を出力する例を示している。
　実施の形態２では、互いに異なるそれぞれの電波受信地点から同じ観測領域が撮像されている２つ以上のレーダ画像を含んでいるレーダ画像群２を取得して、抑圧画像を出力するレーダ画像処理装置１０について説明する。

　実施の形態２のレーダ画像処理装置１０では、位相処理部１２及び画像処理部１３のそれぞれが、レーダ画像群２に含まれている２つのレーダ画像の組み合わせ毎に、処理を実施する。このとき、それぞれの組み合わせに含まれる一方のレーダ画像を第１のレーダ画像、それぞれの組み合わせに含まれる他方のレーダ画像を第２のレーダ画像とする。
　具体的には、位相変化成分算出部２１、位相算出部２２及び位相差算出部２３は、２つのレーダ画像の全ての組み合わせｉについて、位相差Δφ_ｉ（ｘ，ｚ_０）の算出処理が終了するまで、位相差Δφ_ｉ（ｘ，ｚ_０）の算出処理を繰り返し実施する。ｉは、２つのレーダ画像の組み合わせを示す変数である。
　回転量算出部３１、位相回転部３３及び差分算出処理部３４は、２つのレーダ画像の全ての組み合わせｉについて、差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）の算出処理が終了するまで、差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）の算出処理を繰り返し実施する。

　実施の形態２におけるレーダ画像処理装置１０の構成は、実施の形態１のレーダ画像処理装置１０と同様に、図１である。
　実施の形態２における位相処理部１２の構成は、実施の形態１の位相処理部１２と同様に、図２である。
　ただし、レーダ画像群２は、２つ以上のレーダ画像を含んでおり、撮像パラメータ群３は、２つ以上の撮像パラメータを含んでいる。
　図１２は、実施の形態２によるレーダ画像処理装置１０の画像処理部１３を示す構成図である。
　図１３は、位相処理部１２及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。
　図１２及び図１３において、図３及び図４と同一符号は同一又は相当部分を示すので説明を省略する。
　画像合成部３５は、例えば、図１３に示す画像合成回路４７で実現される。
　画像合成部３５は、抑圧画像の生成に用いる重みパラメータｗ_ｉを取得する。
　画像合成部３５は、重みパラメータｗ_ｉを用いて、差分算出処理部３４により組み合わせｉ毎に算出されたそれぞれの差分のうち、それぞれ画素位置が対応している差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）同士を合成する処理を実施する。
　画像合成部３５は、それぞれの合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）を含む抑圧画像を外部に出力する。

　図２では、位相処理部１２の構成要素である位相変化成分算出部２１、位相算出部２２及び位相差算出部２３のそれぞれが、図１３に示すような専用のハードウェアで実現されるものを想定している。
　また、図１２では、画像処理部１３の構成要素である回転量算出部３１、位相回転部３３、差分算出処理部３４及び画像合成部３５のそれぞれが、図１３に示すような専用のハードウェアで実現されるものを想定している。
　即ち、位相処理部１２及び画像処理部１３が、位相変化成分算出回路４１、位相算出回路４２、位相差算出回路４３、回転量算出回路４４、位相回転回路４５、差分算出処理回路４６及び画像合成回路４７で実現されるものを想定している。
　ここで、位相変化成分算出回路４１、位相算出回路４２、位相差算出回路４３、回転量算出回路４４、位相回転回路４５、差分算出処理回路４６及び画像合成回路４７のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ、ＦＰＧＡ、又は、これらを組み合わせたものが該当する。

　位相処理部１２の構成要素及び画像処理部１３の構成要素は、専用のハードウェアで実現されるものに限るものではない。位相処理部１２及び画像処理部１３が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせで実現されるものであってもよい。
　位相処理部１２が、ソフトウェア又はファームウェアなどで実現される場合、位相変化成分算出部２１、位相算出部２２及び位相差算出部２３の処理手順をコンピュータに実行させるためのプログラムが図５に示すメモリ６１に格納される。
　また、画像処理部１３が、ソフトウェア又はファームウェアなどで実現される場合、回転量算出部３１、位相回転部３３、差分算出処理部３４及び画像合成部３５の処理手順をコンピュータに実行させるためのプログラムがメモリ６１に格納される。
　そして、コンピュータのプロセッサ６２がメモリ６１に格納されているプログラムを実行する。

　次に、レーダ画像処理装置１０の動作について説明する。
　位相処理部１２は、レーダ画像群２に含まれている２つ以上のレーダ画像のうち、２つのレーダ画像の組み合わせｉ毎に、位相差Δφ_ｉ（ｘ，ｚ_０）の算出処理を実施する。
　位相変化成分算出部２１は、レーダ画像取得部１１から出力された撮像パラメータ群３の中から、２つのレーダ画像に対応する２つの撮像パラメータの組み合わせを取得する。
　ここでは、組み合わせｉに含まれる一方のレーダ画像が第１のレーダ画像、組み合わせｉに含まれる他方のレーダ画像が第２のレーダ画像であるとする。
　或る組み合わせに含まれる第１のレーダ画像に係る電波受信地点と、他の組み合わせに含まれる第１のレーダ画像に係る電波受信地点とは、互いに異なる地点である。しかし、ここでは、説明の便宜上、どちらの電波受信地点も、第１の電波受信地点Ｐ_１であるとする。
　また、或る組み合わせに含まれる第２のレーダ画像に係る電波受信地点と、他の組み合わせに含まれる第２のレーダ画像に係る電波受信地点とは、互いに異なる地点である。しかし、ここでは、説明の便宜上、どちらの電波受信地点も、第２の電波受信地点Ｐ_２であるとする。
　第１のレーダ画像に対応する撮像パラメータが第１の撮像パラメータ、第２のレーダ画像に対応する撮像パラメータが第２の撮像パラメータであるとする。
　また、位相変化成分算出部２１は、傾き角度αを取得する。
　位相算出部２２は、第１の撮像パラメータと、第２の撮像パラメータと、傾き角度αと、距離ｚ_０とを取得する。

　位相変化成分算出部２１は、レーダ画像におけるスラントレンジ方向の画素の位置ｐｉｘｅｌを式（４）に代入することで、位置ｐｉｘｅｌに対応する傾き面５１での位置ｘを算出する。
　式（４）に代入する位置ｐｉｘｅｌにある画素は、複数の散乱体の反射信号が混在している画素である。
　位相変化成分算出部２１は、距離成分Ｂ_１，２と、オフナディア角θと、距離の平均値Ｒと、放射電波の波長λと、傾き角度αと、観測パスパラメータｐとを用いて、傾き面５１におけるｘ軸方向の位相変化成分φ_ｉ（ｘ）を算出する。
　以下の式（１０）は、位相変化成分算出部２１が用いる位相変化成分φ_ｉ（ｘ）の算出式である。

　位相変化成分算出部２１は、ｘ軸方向の位相変化成分φ_ｉ（ｘ）を位相差算出部２３に出力する。

　位相算出部２２は、距離成分Ｂ_１，２と、オフナディア角θと、距離の平均値Ｒと、放射電波の波長λと、傾き角度αと、距離ｚ_０と、観測パスパラメータｐとを用いて、平行面５２における傾き面５１に対する位相ρ_ｉ（ｚ_０）を算出する。位相ρ_ｉ（ｚ_０）は、どの組み合わせであっても同じである。
　以下の式（１１）は、位相算出部２２が用いる位相ρ_ｉ（ｚ_０）の算出式である。

　位相算出部２２は、位相ρ_ｉ（ｚ_０）を位相差算出部２３に出力する。

　位相差算出部２３は、組み合わせｉ毎に、位相変化成分φ_ｉ（ｘ）と位相ρ_ｉ（ｚ_０）とを用いて、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点Ｐ_１に対する位相と第２の電波受信地点Ｐ_２に対する位相との位相差Δφ_ｉ（ｘ，ｚ_０）を算出する。
　以下の式（１２）は、位相差算出部２３が用いる位相差Δφ_ｉ（ｘ，ｚ_０）の算出式である。

　位相差算出部２３は、それぞれの位相差Δφ_ｉ（ｘ，ｚ_０）を画像処理部１３に出力する。

　回転量算出部３１は、位相差算出部２３から出力されたそれぞれの位相差Δφ_ｉ（ｘ，ｚ_０）を取得する。
　回転量算出部３１は、組み合わせｉ毎に、それぞれの位相差Δφ_ｉ（ｘ，ｚ_０）から、第２のレーダ画像に含まれている複数の画素の位相の回転量ｅｘｐ［ｊ・Δφ_ｉ（ｘ，ｚ_０）］をそれぞれ算出する。
　回転量算出部３１は、それぞれの回転量ｅｘｐ［ｊ・Δφ_ｉ（ｘ，ｚ_０）］を位相回転部３３に出力する。

　位相回転部３３は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、組み合わせｉに含まれている第２のレーダ画像を取得する。
　位相回転部３３は、回転量算出部３１から出力されたそれぞれの回転量ｅｘｐ［ｊ・Δφ_ｉ（ｘ，ｚ_０）］に基づいて、取得した第２のレーダ画像に含まれている複数の画素の位相を回転させる処理を実施する。
　以下の式（１３）は、位相回転部３３による位相の回転処理を示す式である。

　位相回転部３３は、位相回転後の複数の画素を含む第２のレーダ画像を差分算出処理部３４に出力する。

　差分算出処理部３４は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、組み合わせｉに含まれている第１のレーダ画像を取得し、位相回転部３３から出力された位相回転後の複数の画素を含む第２のレーダ画像を取得する。
　差分算出処理部３４は、取得した第１のレーダ画像に含まれている複数の画素と、取得した第２のレーダ画像に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）を算出する。
　以下の式（１４）は、差分算出処理部３４が用いる差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）の算出式である。

　差分算出処理部３４は、それぞれの差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）を画像合成部３５に出力する。
　回転量算出部３１、位相回転部３３及び差分算出処理部３４は、２つのレーダ画像の全ての組み合わせｉについて、差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）の算出処理が終了するまで、差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）の算出処理を繰り返し実施する。

　画像合成部３５は、抑圧画像の生成に用いる重みパラメータｗ_ｉを取得する。
　重みパラメータｗ_ｉは、例えば、ユーザによる手動の操作で画像合成部３５に与えられるようにしてもよいし、図示せぬ外部機器によって画像合成部３５に与えられるようにしてもよい。
　画像合成部３５は、重みパラメータｗ_ｉを用いて、差分算出処理部３４により組み合わせｉ毎に算出されたそれぞれの差分のうち、それぞれ画素位置が対応している差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）同士を合成する。
　画像合成部３５は、それぞれの合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）を含む抑圧画像を外部に出力する。

　全ての組み合わせの差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）同士を合成する方法として、相加平均を求める方法と、相乗平均を求める方法とがある。
　相加平均を求める方法を用いる場合、画像合成部３５は、以下の式（１５）によって、全ての組み合わせの差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）同士を合成する。

　相乗平均を求める方法を用いる場合、画像合成部３５は、以下の式（１６）によって、全ての組み合わせの差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）同士を合成する。

　式（１５）及び式（１６）において、Ｎは、２つのレーダ画像の組み合わせ数である。

　ここで、図１４は、実施の形態１のレーダ画像処理装置１０のように、レーダ画像群２に含まれているレーダ画像が２つだけである場合の、１つの画素に混在している複数の反射信号を示す説明図である。
　レーダ画像群２に含まれているレーダ画像が２つだけの場合、図１４に示すように、差分算出処理部３４が差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）の算出処理を実施することで、複数のｎｕｌｌ点が形成されることがある。
　図１４の例では、“１”が付されている反射信号、“２”が付されている反射信号及び３”が付されている反射信号の全てにｎｕｌｌ点が形成されている。
　したがって、図１４の例では、“１”が付されている反射信号、“２”が付されている反射信号及び“３”が付されている反射信号の全てが抑圧されてしまう。

　図１５は、実施の形態２のレーダ画像処理装置１０のように、レーダ画像群２に含まれているレーダ画像が２つ以上である場合の、１つの画素に混在している複数の反射信号を示す説明図である。
　図１５では、レーダ画像群２に含まれているレーダ画像がＭ個であり、Ｐ_Ｍは、第Ｍのレーダ画像が撮像されたときのプラットフォームの位置である。
　レーダ画像群２に含まれているレーダ画像が２つ以上であり、画像合成部３５が、それぞれ画素位置が対応している差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）同士を合成することで、形成されるｎｕｌｌ点の数が、レーダ画像が２つの場合よりも減少する。
　図１５の例では、形成されるｎｕｌｌ点の数が１つであり、“２”が付されている反射信号には、ｎｕｌｌ点が形成されていない。

　以上の実施の形態２は、差分算出処理部３４により組み合わせｉ毎に算出されたそれぞれの差分のうち、それぞれ画素位置が対応している差分ΔＳ_ｉ（ｐｉｘｅｌ，ｌｉｎｅ）同士を合成する画像合成部３５を備えるように、レーダ画像処理装置１０を構成した。したがって、レーダ画像処理装置１０は、形成されるｎｕｌｌ点の数を減らして、残す必要がある反射信号の抑圧を防止することができる。

実施の形態３．
　実施の形態２のレーダ画像処理装置１０は、合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）を抑圧画像として出力する例を示している。
　実施の形態３では、画像合成部３５による合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）から、１つの画素に混在している複数の反射信号が抽出された画像を算出するレーダ画像処理装置１０について説明する。

　実施の形態３におけるレーダ画像処理装置１０の構成は、実施の形態１，２のレーダ画像処理装置１０と同様に、図１である。
　実施の形態３における位相処理部１２の構成は、実施の形態１，２の位相処理部１２と同様に、図２である。
　ただし、レーダ画像群２は、２つ以上のレーダ画像を含んでおり、撮像パラメータ群３は、２つ以上の撮像パラメータを含んでいる。
　図１６は、実施の形態３によるレーダ画像処理装置１０の画像処理部１３を示す構成図である。
　図１７は、位相処理部１２及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。
　図１６及び図１７において、図３、図４、図１２及び図１３と同一符号は同一又は相当部分を示すので説明を省略する。
　抽出画像算出部３６は、例えば、図１７に示す抽出画像算出回路４８で実現される。
　抽出画像算出部３６は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、第１のレーダ画像を取得し、画像合成部３５から出力されたそれぞれの合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）を取得する。
　抽出画像算出部３６は、第１のレーダ画像に含まれている複数の画素の画素値と、それぞれの合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）とから、１つの画素に混在している複数の反射信号が抽出された画像を算出する処理を実施する。

　図２では、位相処理部１２の構成要素である位相変化成分算出部２１、位相算出部２２及び位相差算出部２３のそれぞれが、図１７に示すような専用のハードウェアで実現されるものを想定している。
　また、図１６では、画像処理部１３の構成要素である回転量算出部３１、位相回転部３３、差分算出処理部３４、画像合成部３５及び抽出画像算出部３６のそれぞれが、図１７に示すような専用のハードウェアで実現されるものを想定している。
　即ち、位相処理部１２及び画像処理部１３が、位相変化成分算出回路４１、位相算出回路４２、位相差算出回路４３、回転量算出回路４４、位相回転回路４５、差分算出処理回路４６、画像合成回路４７及び抽出画像算出回路４８で実現されるものを想定している。
　位相処理部１２の構成要素及び画像処理部１３の構成要素は、専用のハードウェアで実現されるものに限るものではない。位相処理部１２及び画像処理部１３が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせで実現されるものであってもよい。

　次に、レーダ画像処理装置１０の動作について説明する。
　ただし、レーダ画像処理装置１０が、抽出画像算出部３６を備えている点以外は、実施の形態２のレーダ画像処理装置１０と同様であるため、ここでは、抽出画像算出部３６の動作だけを説明する。

　抽出画像算出部３６は、レーダ画像取得部１１から出力されたレーダ画像群２の中から、第１のレーダ画像を取得し、画像合成部３５から出力されたそれぞれの合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）を取得する。
　抽出画像算出部３６は、第１のレーダ画像に含まれている複数の画素の画素値と、それぞれの合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）とから、複数の反射信号が混在している画素の画素値Ｓ_ｅｘｔ（ｐｉｘｅｌ，ｌｉｎｅ）を算出する。
　以下の式（１７）は、抽出画像算出部３６が用いる画素値Ｓ_ｅｘｔ（ｐｉｘｅｌ，ｌｉｎｅ）の算出式である。

　抽出画像算出部３６は、１つの画素に混在している複数の反射信号が抽出された画像として、画素値Ｓ_ｅｘｔ（ｐｉｘｅｌ，ｌｉｎｅ）を有する画素を含む画像を外部に出力する。

　以上の実施の形態３は、第１のレーダ画像に含まれている複数の画素の画素値と、それぞれの合成後の差分Ｓ_ｓｕｐ（ｐｉｘｅｌ，ｌｉｎｅ）とから、１つの画素に混在している複数の反射信号が抽出された画像を算出する抽出画像算出部３６を備えるように、レーダ画像処理装置１０を構成した。したがって、レーダ画像処理装置１０は、反射信号が抑圧された抑圧画像だけでなく、反射信号が抽出された抽出画像を出力することができる。

実施の形態４．
　実施の形態４では、第１のレーダ画像における反射信号の抑圧対象範囲である第１の抑圧対象範囲と、第２のレーダ画像における反射信号の抑圧対象範囲である第２の抑圧対象範囲との指定を受け付けるレーダ画像処理装置について説明する。

　図１８は、実施の形態４によるレーダ画像処理装置１０を示す構成図である。図１８において、図１と同一符号は同一又は相当部分を示している。
　レーダ画像処理装置１０は、レーダ画像取得部１４、位相処理部１５及び画像処理部１３を備えている。
　レーダ画像取得部１４は、レーダ１から送信されたレーダ画像群２を取得する。レーダ画像群２は、第１のレーダ画像及び第２のレーダ画像を含んでいる。
　レーダ画像取得部１４は、取得したレーダ画像群２を位相処理部１５及び画像処理部１３のそれぞれに出力する。

　位相処理部１５は、範囲指定受付部７０及び位相差算出部７１を備えている。
　位相処理部１５は、第１の抑圧対象範囲として、第１のレーダ画像における反射信号の抑圧対象範囲の指定を受け付け、第２の抑圧対象範囲として、第２のレーダ画像における反射信号の抑圧対象範囲の指定を受け付ける。
　位相処理部１５は、第１の抑圧対象範囲に含まれている複数の画素と、第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点に対する位相と第２の電波受信地点に対する位相との位相差Δφを算出する。
　位相処理部１５は、算出した位相差Δφを、位相差Δφ（ｘ，ｚ_０）として、画像処理部１３に出力する。

　図１８に示す画像処理部１３は、図１に示す画像処理部１３と同様に、画素位置が対応している画素同士の画素値の差分を算出する。
　ただし、図１８に示す画像処理部１３は、図１に示す画像処理部１３と異なり、第２の抑圧対象範囲に含まれている複数の画素の位相の回転量をそれぞれ算出する。
　また、図１８に示す画像処理部１３は、第２の抑圧対象範囲に含まれている複数の画素の位相を回転させ、第１の抑圧対象範囲に含まれている複数の画素と第２の抑圧対象範囲に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出する。

　図１９は、実施の形態４によるレーダ画像処理装置１０の位相処理部１５を示す構成図である。
　図２０は、位相処理部１５及び画像処理部１３におけるそれぞれのハードウェアを示すハードウェア構成図である。図２０において、図４と同一符号は同一又は相当部分を示している。
　図１９において、範囲指定受付部７０は、例えば、図２０に示す範囲指定受付回路８１によって実現される。
　範囲指定受付部７０は、レーダ画像取得部１４から出力されたレーダ画像群２を取得する。
　範囲指定受付部７０は、第１の抑圧対象範囲の指定を受け付け、第２の抑圧対象範囲の指定を受け付ける。
　範囲指定受付部７０は、第１の抑圧対象範囲での第１のレーダ画像及び第２の抑圧対象範囲での第２のレーダ画像のそれぞれを位相差算出部７１に出力する。

　位相差算出部７１は、平面位相変化算出部７２及び位相変化決定部７３を備えている。
　位相差算出部７１は、第１の抑圧対象範囲に含まれている複数の画素と、第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点に対する位相と第２の電波受信地点に対する位相との位相差を算出する。

　平面位相変化算出部７２は、例えば、図２０に示す平面位相変化算出回路８２によって実現される。
　平面位相変化算出部７２は、範囲指定受付部７０から出力された第１の抑圧対象範囲での第１のレーダ画像及び第２の抑圧対象範囲での第２のレーダ画像のそれぞれを取得する。
　平面位相変化算出部７２は、第１の抑圧対象範囲での反射信号と、第２の抑圧対象範囲での反射信号との干渉位相の空間的な変化ψとして、２次元の位相変化成分ψを算出する。
　平面位相変化算出部７２は、第１の抑圧対象範囲での第１のレーダ画像と第２の抑圧対象範囲での第２のレーダ画像との干渉位相∠（Ｓ_１／Ｓ_２）を算出する。
　平面位相変化算出部７２は、算出した２次元の位相変化成分ψ及び干渉位相∠（Ｓ_１／Ｓ_２）のそれぞれを位相変化決定部７３に出力する。

　位相変化決定部７３は、例えば、図２０に示す位相変化決定回路８３によって実現される。
　位相変化決定部７３は、平面位相変化算出部７２により算出された２次元の位相変化成分ψ及び干渉位相∠（Ｓ_１／Ｓ_２）のそれぞれに基づいて、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点に対する位相と第２の電波受信地点に対する位相との位相差Δφを決定する。
　位相変化決定部７３は、決定した位相差Δφを、位相差Δφ（ｘ，ｚ_０）として、画像処理部１３に出力する。

　図１９では、位相処理部１５の構成要素である範囲指定受付部７０、平面位相変化算出部７２及び位相変化決定部７３のそれぞれが、図２０に示すような専用のハードウェアで実現されるものを想定している。
　また、図３では、画像処理部１３の構成要素である回転量算出部３１、位相回転部３３及び差分算出処理部３４のそれぞれが、図２０に示すような専用のハードウェアで実現されるものを想定している。
　即ち、位相処理部１５及び画像処理部１３が、範囲指定受付回路８１、平面位相変化算出回路８２、位相変化決定回路８３、回転量算出回路４４、位相回転回路４５及び差分算出処理回路４６で実現されるものを想定している。

　ここで、範囲指定受付回路８１、平面位相変化算出回路８２、位相変化決定回路８３、回転量算出回路４４、位相回転回路４５及び差分算出処理回路４６のそれぞれは、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ、ＦＰＧＡ、又は、これらを組み合わせたものが該当する。

　位相処理部１５の構成要素及び画像処理部１３の構成要素は、専用のハードウェアで実現されるものに限るものではない。位相処理部１５及び画像処理部１３が、ソフトウェア、ファームウェア、又は、ソフトウェアとファームウェアとの組み合わせで実現されるものであってもよい。
　位相処理部１５が、ソフトウェア又はファームウェアなどで実現される場合、範囲指定受付部７０、平面位相変化算出部７２及び位相変化決定部７３の処理手順をコンピュータに実行させるためのプログラムが図５に示すメモリ６１に格納される。
　また、画像処理部１３が、ソフトウェア又はファームウェアなどで実現される場合、回転量算出部３１、位相回転部３３及び差分算出処理部３４の処理手順をコンピュータに実行させるためのプログラムがメモリ６１に格納される。
　そして、図５に示すコンピュータのプロセッサ６２がメモリ６１に格納されているプログラムを実行する。

　次に、図１９に示すレーダ画像処理装置１０の動作について説明する。
　ただし、レーダ画像取得部１４及び位相処理部１５以外は、図１に示すレーダ画像処理装置１０と同様であるため、ここでは、主にレーダ画像取得部１４及び位相処理部１５の動作を説明する。

　レーダ画像取得部１４は、レーダ１から送信されたレーダ画像群２を取得する。レーダ画像群２は、第１のレーダ画像及び第２のレーダ画像を含んでいる。
　レーダ画像取得部１４は、取得したレーダ画像群２を位相処理部１５及び画像処理部１３のそれぞれに出力する。

　範囲指定受付部７０は、レーダ画像取得部１４から出力されたレーダ画像群２を取得する。
　範囲指定受付部７０は、レーダ画像群２に含まれている第１のレーダ画像における反射信号の抑圧対象範囲である第１の抑圧対象範囲と、レーダ画像群２に含まれている第２のレーダ画像における反射信号の抑圧対象範囲である第２の抑圧対象範囲との指定を受け付ける。
　範囲指定受付部７０は、指定を受け付けた第１の抑圧対象範囲での第１のレーダ画像と、指定を受け付けた第２の抑圧対象範囲での第２のレーダ画像とを平面位相変化算出部７２に出力する。

　第１の抑圧対象範囲が、例えば、図２１Ａに示すような平行四辺形の領域であれば、ユーザは、マウス等のマンマシンインタフェースを用いて、平行四辺形の領域における４つの頂点の座標を指定する。ここでは、第１の抑圧対象範囲が、平行四辺形の領域であるとしているが、これは一例に過ぎず、第１の抑圧対象範囲は、例えば、三角形の領域であってもよいし、矩形の領域であってもよいし、五角形の領域であってもよい。
　範囲指定受付部７０は、マンマシンインタフェースを用いて、ユーザにより指定された領域の頂点の座標データを受けると、座標データが示す頂点を有する領域を、第１の抑圧対象範囲に決定する。
　図２１Ａは、第１の抑圧対象範囲を示す説明図である。

　第２の抑圧対象範囲が、例えば、図２１Ｂに示すような平行四辺形の領域であれば、ユーザは、マウス等のマンマシンインタフェースを用いて、平行四辺形の領域における４つの頂点の座標を指定する。ここでは、第２の抑圧対象範囲が、平行四辺形の領域であるとしているが、これは一例に過ぎず、第２の抑圧対象範囲は、例えば、三角形の領域であってもよいし、矩形の領域であってもよいし、五角形の領域であってもよい。
　範囲指定受付部７０は、マンマシンインタフェースを用いて、ユーザにより指定された領域の頂点の座標データを受けると、座標データが示す頂点を有する領域を、第２の抑圧対象範囲に決定する。
　図２１Ｂは、第２の抑圧対象範囲を示す説明図である。

　図１９に示すレーダ画像処理装置１０では、範囲指定受付部７０が、第１の抑圧対象範囲と第２の抑圧対象範囲との双方の指定を受け付けている。しかし、これは一例に過ぎず、範囲指定受付部７０は、第１の抑圧対象範囲の指定のみを受け付け、第２のレーダ画像において、第１の抑圧対象範囲と画像上の位置が同じ範囲を、第２の抑圧対象範囲に決定するようにしてもよい。

　平面位相変化算出部７２は、範囲指定受付部７０から出力された第１の抑圧対象範囲での第１のレーダ画像及び第２の抑圧対象範囲での第２のレーダ画像のそれぞれを取得する。
　平面位相変化算出部７２は、第１の抑圧対象範囲での反射信号と、第２の抑圧対象範囲での反射信号との干渉位相の空間的な変化ψとして、２次元の位相変化成分ψを算出する。
　平面位相変化算出部７２は、算出した２次元の位相変化成分ψを位相変化決定部７３に出力する。
　以下、平面位相変化算出部７２による２次元の位相変化成分ψの算出処理を具体的に説明する。

　平面位相変化算出部７２は、第１の抑圧対象範囲に含まれている複数の画素の中から、画素位置が（ｐｉｘｅｌ，ｌｉｎｅ）の画素の画素値Ｓ_１（ｐｉｘｅｌ，ｌｉｎｅ）を取得する。
　平面位相変化算出部７２は、第２の抑圧対象範囲に含まれている複数の画素の中から、画素位置が（ｐｉｘｅｌ，ｌｉｎｅ）の画素の画素値Ｓ_２（ｐｉｘｅｌ，ｌｉｎｅ）を取得する。
　平面位相変化算出部７２は、画素値Ｓ_１（ｐｉｘｅｌ，ｌｉｎｅ）と画素値Ｓ_２（ｐｉｘｅｌ，ｌｉｎｅ）との複素共役積（Ｓ_１×Ｓ^＊ _２）として、（Ｓ_１／Ｓ_２）を算出する。＊は、共役複素数を示す記号である。
　平面位相変化算出部７２は、第１及び第２の抑圧対象範囲に含まれている全ての画素について、画素値Ｓ_１（ｐｉｘｅｌ，ｌｉｎｅ）と画素値Ｓ_２（ｐｉｘｅｌ，ｌｉｎｅ）との複素共役積（Ｓ_１×Ｓ^＊ _２）として、（Ｓ_１／Ｓ_２）を算出する。

　平面位相変化算出部７２は、第１及び第２の抑圧対象範囲に含まれている全ての画素についての（Ｓ_１／Ｓ_２）を２次元フーリエ変換する。
　平面位相変化算出部７２は、以下の式（１８）に示すように、２次元フーリエ変換の変換結果である周波数領域の信号から、干渉位相の空間的な変化ψとして、２次元の位相変化成分ψを抽出する。
　２次元の位相変化成分ψは、周波数領域の信号に含まれているレンジ方向の周波数成分ｆ_{ｒａｎｇｅ}の中のピーク値である、ピークを有するレンジ方向の周波数成分ｆ_{ｒａｎｇｅ}に対応する角周波数ｆ_{ｒａｎｇｅ}Δｓｌを含むものである。
　また、２次元の位相変化成分ψは、周波数領域の信号に含まれているアジマス方向の周波数成分ｆ_{ａｚｉｍｕｔｈ}の中のピーク値である、ピークを有するアジマス方向の周波数成分ｆ_{ａｚｉｍｕｔｈ}に対応する角周波数ｆ_{ａｚｉｍｕｔｈ}Δａｚを含むものである。

　ψ＝（ｆ_{ｒａｎｇｅ}Δｓｌ，ｆ_{ａｚｉｍｕｔｈ}Δａｚ）　（１８）
　式（１８）において、Δａｚはアジマス方向のスペーシングである。

　平面位相変化算出部７２は、ピークを有するレンジ方向の周波数成分ｆ_{ｒａｎｇｅ}と、ピークを有するアジマス方向の周波数成分ｆ_{ａｚｉｍｕｔｈ}とから、第１の抑圧対象範囲での第１のレーダ画像と第２の抑圧対象範囲での第２のレーダ画像との干渉位相∠（Ｓ_１／Ｓ_２）を算出する。∠（Ｓ_１／Ｓ_２）は、（Ｓ_１／Ｓ_２）の偏角である。レンジ方向の周波数成分ｆ_{ｒａｎｇｅ}と、アジマス方向の周波数成分ｆ_{ａｚｉｍｕｔｈ}とから、（Ｓ_１／Ｓ_２）の偏角を算出する処理自体は、公知の技術であるため詳細な説明を省略する。
　平面位相変化算出部７２は、第１及び第２の抑圧対象範囲に含まれている全ての画素について、干渉位相∠（Ｓ_１／Ｓ_２）を算出する。
　平面位相変化算出部７２は、第１及び第２の抑圧対象範囲に含まれている全ての画素についての干渉位相∠（Ｓ_１／Ｓ_２）と、２次元の位相変化成分ψとを位相変化決定部７３に出力する。

　位相変化決定部７３は、第１の抑圧対象範囲に含まれている複数の画素と、第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、以下の式（１９）に示すように、干渉位相∠（Ｓ_１／Ｓ_２）から２次元の位相変化成分ψを減算することで、Δψを算出する。

　Δψ＝∠（Ｓ_１／Ｓ_２）－ψ　　　　　　　　　　　（１９）

　位相変化決定部７３は、算出したΔψが０であるか否かを確認し、位相差Δφが０であれば、２次元の位相変化成分ψを、位相差Δφ（ｘ，ｚ_０）として、画像処理部１３の回転量算出部３１に出力する。
　算出したΔψが０であれば、干渉位相∠（Ｓ_１／Ｓ_２）が２次元の位相変化成分ψと一致しているため、回転量算出部３１による回転量ｅｘｐ［ｊ・Δφ（ｘ，ｚ_０）］の算出精度が向上することが想定される。

　位相変化決定部７３は、算出したΔψが０でなければ、例えば、反射信号の抑圧対象範囲の再指定を促すメッセージを図示せぬ表示器に出力する。
　ユーザが表示器に出力されたメッセージを確認して、第１及び第２の抑圧対象範囲を再指定すれば、平面位相変化算出部７２が、ユーザによって再指定された第１及び第２の抑圧対象範囲に基づいて、２次元の位相変化成分ψ及び干渉位相∠（Ｓ_１／Ｓ_２）のそれぞれを再度算出する。
　位相変化決定部７３は、式（１９）に示すように、干渉位相∠（Ｓ_１／Ｓ_２）から２次元の位相変化成分ψを減算することで、Δψを再度算出する。
　位相変化決定部７３は、再度算出したΔψが０であるか否かを確認し、位相差Δφが０であれば、２次元の位相変化成分ψを、位相差Δφ（ｘ，ｚ_０）として、画像処理部１３の回転量算出部３１に出力する。
　位相変化決定部７３は、反射信号の抑圧対象範囲の再指定を促すメッセージを図示せぬ表示器に出力しても、ユーザが第１及び第２の抑圧対象範囲を再指定しなければ、２次元の位相変化成分ψを、位相差Δφ（ｘ，ｚ_０）として、画像処理部１３の回転量算出部３１に出力するようにしてもよい。

　以上の実施の形態４では、第１の抑圧対象範囲に含まれている複数の画素と、第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、第１の電波受信地点に対する位相と第２の電波受信地点に対する位相との位相差を算出する位相差算出部７１と、位相差算出部７１により算出されたそれぞれの位相差から、第２の抑圧対象範囲に含まれている複数の画素の位相の回転量をそれぞれ算出する回転量算出部３１とを設け、差分算出部３２が、回転量算出部３１により算出されたそれぞれの回転量に基づいて、第２の抑圧対象範囲に含まれている複数の画素の位相を回転させ、第１の抑圧対象範囲に含まれている複数の画素と第２の抑圧対象範囲に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出するように、レーダ画像処理装置１０を構成した。したがって、レーダ画像処理装置１０は、第１の電波受信地点に対する位相と、第２の電波受信地点に対する位相との位相差が零でない反射信号についても抑圧することができる。

　図１に示すレーダ画像処理装置１０では、位相処理部１２が、位相差Δφ（ｘ，ｚ_０）を算出するために、撮像パラメータ３、傾き角度α及び距離ｚ_０のそれぞれを取得する必要がある。
　図１８に示すレーダ画像処理装置１０では、位相処理部１５が、第１及び第２の抑圧対象範囲を取得できれば、位相処理部１５が、撮像パラメータ３、傾き角度α及び距離ｚ_０のそれぞれを取得する必要がない。

　図１８に示すレーダ画像処理装置１０では、範囲指定受付部７０が、平行四辺形の領域における４つの頂点の座標から、第１及び第２の抑圧対象範囲を決定している。
　ユーザによる４つの頂点の座標の指定は、抑圧対象範囲の枠を指定していることに相当する。
　抑圧対象範囲の指定は、枠の指定に限るものではなく、第１及び第２の抑圧対象範囲内の３つ以上の点を指定することで、第１及び第２の抑圧対象範囲を指定するものであってもよい。

　第１の抑圧対象範囲が、例えば、図２２Ａに示すような平行四辺形の領域であれば、ユーザは、マウス等のマンマシンインタフェースを用いて、平行四辺形の領域内の３つ以上の点を指定する。
　範囲指定受付部７０は、マンマシンインタフェースを用いて、ユーザにより指定された３つ以上の点の座標データを受けると、３つ以上の点が通る２次元平面を特定し、特定した２次元平面を第１の抑圧対象範囲に決定する。
　範囲指定受付部７０は、第２の抑圧対象範囲についても、第１の抑圧対象範囲と同様に、３つ以上の点が通る２次元平面を特定し、特定した２次元平面を第２の抑圧対象範囲に決定する。
　図２２Ａは、第１の抑圧対象範囲を示す説明図であり、図２２Ｂは、第２の抑圧対象範囲を示す説明図である。

　ここでは、範囲指定受付部７０が、第１の抑圧対象範囲と第２の抑圧対象範囲との双方の指定を受け付けている。しかし、これは一例に過ぎず、範囲指定受付部７０は、第１の抑圧対象範囲の指定のみを受け付け、第２のレーダ画像において、第１の抑圧対象範囲と画像上の位置が同じ範囲を、第２の抑圧対象範囲に決定するようにしてもよい。

　なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。

　この発明は、第１のレーダ画像に含まれている複数の画素と第２のレーダ画像に含まれている位相回転後の複数の画素との差分をそれぞれ算出するレーダ画像処理装置及びレーダ画像処理方法に適している。

　１　レーダ、２　レーダ画像群、３　撮像パラメータ群、１０　レーダ画像処理装置、１１，１４　レーダ画像取得部、１２，１５　位相処理部、１３　画像処理部、２１　位相変化成分算出部、２２　位相算出部、２３　位相差算出部、３１　回転量算出部、３２　差分算出部、３３　位相回転部、３４　差分算出処理部、３５　画像合成部、３６　抽出画像算出部、４１　位相変化成分算出回路、４２　位相算出回路、４３　位相差算出回路、４４　回転量算出回路、４５　位相回転回路、４６　差分算出処理回路、４７　画像合成回路、４８　抽出画像算出回路、５１　傾き面、５２　平行面、６１　メモリ、６２　プロセッサ、７０　範囲指定受付部、７１　位相差算出部、７２　平面位相変化算出部、７３　位相変化決定部、８１　範囲指定受付回路、８２　平面位相変化算出回路、８３　位相変化決定回路。

Claims (5)

1. 　第１の電波受信地点から観測領域が撮像されている第１のレーダ画像における反射信号の抑圧対象範囲である第１の抑圧対象範囲に含まれている複数の画素と、第２の電波受信地点から前記観測領域が撮像されている第２のレーダ画像における反射信号の抑圧対象範囲である第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、前記第１の電波受信地点に対する位相と前記第２の電波受信地点に対する位相との位相差を算出する位相差算出部と、
   　前記位相差算出部により算出されたそれぞれの位相差から、前記第２の抑圧対象範囲に含まれている複数の画素の位相の回転量をそれぞれ算出する回転量算出部と、
   　前記回転量算出部により算出されたそれぞれの回転量に基づいて、前記第２の抑圧対象範囲に含まれている複数の画素の位相を回転させ、前記第１の抑圧対象範囲に含まれている複数の画素と前記第２の抑圧対象範囲に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出する差分算出部と、
   　を備えたレーダ画像処理装置。
2. 　前記位相差算出部は、前記第１の抑圧対象範囲での反射信号と、前記第２の抑圧対象範囲での反射信号との干渉位相の空間的な変化を算出し、前記干渉位相の空間的な変化から、前記第１の電波受信地点に対する位相と前記第２の電波受信地点に対する位相との位相差を算出することを特徴とする請求項１記載のレーダ画像処理装置。
3. 　前記第１の抑圧対象範囲と、前記第２の抑圧対象範囲との指定を受け付ける範囲指定受付部を備え、
   　前記範囲指定受付部は、指定を受け付けた第１の抑圧対象範囲及び指定を受け付けた第２の抑圧対象範囲のそれぞれを前記位相差算出部に出力することを特徴とする請求項１記載のレーダ画像処理装置。
4. 　前記第１の抑圧対象範囲の指定を受け付け、前記第２のレーダ画像において、前記第１の抑圧対象範囲と画像上の位置が同じ範囲を、前記第２の抑圧対象範囲に決定する範囲指定受付部を備え、
   　前記範囲指定受付部は、指定を受け付けた第１の抑圧対象範囲及び決定した第２の抑圧対象範囲のそれぞれを前記位相差算出部に出力することを特徴とする請求項１記載のレーダ画像処理装置。
5. 　位相差算出部が、第１の電波受信地点から観測領域が撮像されている第１のレーダ画像における反射信号の抑圧対象範囲である第１の抑圧対象範囲に含まれている複数の画素と、第２の電波受信地点から前記観測領域が撮像されている第２のレーダ画像における反射信号の抑圧対象範囲である第２の抑圧対象範囲に含まれている複数の画素とのうち、それぞれの画素位置が対応している画素毎に、１つの画素に混在している複数の反射信号のそれぞれにおいて、前記第１の電波受信地点に対する位相と前記第２の電波受信地点に対する位相との位相差を算出し、
   　回転量算出部が、前記位相差算出部により算出されたそれぞれの位相差から、前記第２の抑圧対象範囲に含まれている複数の画素の位相の回転量をそれぞれ算出し、
   　差分算出部が、前記回転量算出部により算出されたそれぞれの回転量に基づいて、前記第２の抑圧対象範囲に含まれている複数の画素の位相を回転させ、前記第１の抑圧対象範囲に含まれている複数の画素と前記第２の抑圧対象範囲に含まれている位相回転後の複数の画素とのうち、それぞれ画素位置が対応している画素同士の画素値の差分を算出する
   　レーダ画像処理方法。

Images