JP2010014700A - 画像レーダ装置 - Google Patents

Abstract

【課題】目標の相対運動と３次元形状を推定することのできる画像レーダ装置を得る。
【解決手段】目標に対して送信電波を照射するための送信機１及び送受信アンテナ３と、前記送信電波が前記目標で反射された反射電波を受信する、互いに異なる位置に配置された複数の受信用補助アンテナ４，５および受信機６〜８と、前記複数の受信用アンテナおよび受信機で受信された前記受信電波に基づき、前記目標のレーダ動画像を再生するレーダ画像再生部９と、前記レーダ画像再生部で連続的に再生された複数のレーダ動画像を用いて目標の運動及び形状を推定する目標運動・形状推定手段１１とを備える。
【選択図】図１

Description

この発明は、送信電波と反射電波からなる受信電波とに基づき遠方の目標を観測して、目標の形状、反射強度分布、種類などを特定する画像レーダ装置に関するものである。

従来の画像レーダ装置においては、異なる時刻に送受信された電波から得られる２枚のレーダ画像（以下、単に「画像」ともいう）を比較して、レーダ画像上の輝点の画像上での速度分布を収集し、適当な方法で画像レーダ装置と目標との相対運動を推定したうえで、相対運動情報を用いて目標上の電波反射点の３次元空間における位置を算出することにより、目標の３次元形状を得ている（例えば、特許文献１参照）。

特開２００１−３３７１６３号公報（第５−８頁、図１−図５）

レーダ画像は、目標の三次元形状を、レーダの視線方向と目標の相対運動で決まる二次元の投影面に投影したものとして得られる。レーダの観測においては、目標は、電波を反射する複数の孤立反射点の集合と考えることが出来る。レーダ画像においては、これらの孤立反射点の位置を二次元平面に投影した位置に、輝点が現れる。

ところで、レーダ画像を用いて目標の三次元形状を推定するためには、レーダ画像の投影面の向きを知る必要がある。前述の通り、レーダ画像の投影面の向きを決定する要素の一つに、レーダと目標の相対運動が含まれるが、三次元形状を推定する従来の画像レーダ装置においては、レーダと目標との相対運動の推定方法が明確に設定されておらず、画像レーダ装置以外の何らかの手段で画像レーダ装置と目標との相対運動を別途推定する必要があった。特に、ＩＳＡＲ（Inverse Synthetic Aperture Radar）方式によるレーダ画像生成の際は、目標運動は一般に未知であり、目標の相対運動を事前に推定するのは困難であった。

この発明は、上記課題を解決するためになされたものであり、目標の相対運動と３次元形状を推定することのできる画像レーダ装置を得ることを目的とする。

この発明に係る画像レーダ装置は、目標に対して送信電波を照射する送信手段と、前記送信電波が前記目標で反射された反射電波を受信する、互いに異なる位置に配置された複数の受信用アンテナおよび受信機と、前記複数の受信用アンテナおよび受信機で受信された前記受信電波に基づき、前記目標のレーダ動画像を再生するレーダ画像再生部と、前記レーダ画像再生部で連続的に再生された複数のレーダ動画像を用いて目標の運動及び形状を推定する目標運動・形状推定手段とを備えたものである。

この発明によれば、複数のレーダ動画像から観測された視線方向の位置およびドップラー周波数に加えて、観測信号の位相差の情報を利用できるように構成したので、目標の運動と三次元形状を同時に推定することが可能となる。

この発明の実施の形態１に係る画像レーダ装置の機能構成を示すブロック図である。 この発明の実施の形態１に係る画像レーダ装置の動作原理を説明するもので、観測のジオメトリを示す図である。 この発明の実施の形態１に係る画像レーダ装置の動作原理を説明するもので、各アンテナで受信された信号を用いて得られるＩＳＡＲ動画像の概念図である。 この発明の実施の形態１に係る画像レーダ装置の動作原理を説明するもので、観測された位相差の値に直線を当てはめて位相差の変化率を求める概念図である。 この発明の実施の形態２に係る画像レーダ装置の機能構成を示すブロック図である。 この発明の実施の形態３に係る画像レーダ装置の機能構成を示すブロック図である。 この発明の実施の形態４に係る画像レーダ装置の機能構成を示すブロック図である。 この発明の実施の形態５に係る画像レーダ装置の機能構成を示すブロック図である。 この発明の実施の形態５に係る画像レーダ装置の動作を説明するための説明図である。 図９のレンジ方向加速度利用型目標形状・角速度算出部１７ａの内部構成を示すブロック図ある。 この発明の実施の形態６に係る画像レーダ装置の機能構成を示すブロック図である。 この発明の実施の形態６に係る画像レーダ装置の動作を説明するための説明図である。 図１１の加速度利用型目標形状・角速度算出部１７ｂの内部構成を示すブロック図ある。 この発明の実施の形態７に係る画像レーダ装置の機能構成を示すブロック図である。 図１４の時変角速度考慮型目標形状・角速度算出部１７ｃの内部構成を示すブロック図ある。

実施の形態１．
図１は、この発明の実施の形態１に係る画像レーダ装置の機能構成を示すブロック図である。また、図２乃至図４はこの発明の実施の形態１に係る画像レーダ装置の動作原理を説明するための説明図である。以下では、まず、図２乃至図４を用いて、本実施の形態１に係る画像レーダ装置の動作原理を説明する。なお、本明細書においては、画像レーダの方式としてＩＳＡＲ方式を対象として記載するが、ＳＡＲ（Synthetic Aperture Radar）方式や、その他の画像レーダ方式を用いたものについても、同様の原理が成立することは明らかである。

図２は、観測のジオメトリを示す図である。まず、目標は複数の孤立反射点から成る一つの剛体と考える。図２において、１９はこの目標を構成する孤立反射点のうちの一つＰ_ｋ（ｋ＝１，２，・・・，Ｋ；Ｋは孤立反射点の数）の位置を示す。また、孤立反射点Ｐ_ｋの位置ベクトルを

とする。

次に、ここでは、目標が等角速度運動をしていることを前提とし、その角速度ベクトルを

とする。

また、図２において、２０は目標の回転を表す角速度ベクトル

を表す。一般に、船舶や航空機などの目標が等角速度運動を続けることはないが、短い時間の間、近似的に等角速度運動をしているとみなすことは可能である。また、図２において、２１は孤立反射点Ｐ_ｋの速度を表すベクトル

である。目標上の孤立反射点は全て同じ角速度で運動するため、各孤立反射点の速度

は、次式で表される。

また、図２に示すとおり、目標回転軸上の任意の点Ｏを原点とする直交座標系ｘ、ｙ、ｚ座標を定義する。ｘ、ｙ、ｚ座標の軸の向きは任意であるが、例えばｘ、ｙ平面が水平面を表し、ｚが鉛直方向上向きを表すと考えて良い。

なお、本明細書において、

のように上線のついた変数は空間内のベクトルを表し、

は変数ｘの時間微分を表す。また、

のように二重上線のついた変数は行列を表し、

のように

記号のついた変数は長さが１の単位ベクトルを表す。

図２において、３、４、５はそれぞれ送受信アンテナと受信用補助アンテナの位置を表している。また、送受信アンテナ３、受信用補助アンテナ４、受信用補助アンテナ５の位置ベクトルを、それぞれ、

とする。さらに、送受信アンテナ３からそれぞれ受信用補助アンテナ４と受信用補助アンテナ５へ向かうベクトルを基線ベクトルと呼び、次式で定義する。

また、図２において、２２は、送受信アンテナ３から原点に向かう方向（レーダの視線方向と一致するものとする）を表す単位ベクトルであり、次式で定義される。

なお、以下では、送受信アンテナ３から目標までの距離

が、基線ベクトルの長さと目標のサイズに対して十分に長いものとする。

観測において、送受信アンテナ３は、繰り返し周期（ＰＲＩ：Pulse Repetition Interval）ｔ_ｐｒｉで広帯域パルスの送受信を繰り返す。また、送受信アンテナ３と同時に受信用補助アンテナ４および受信用補助アンテナ５においても信号を受信する。すると、各アンテナで受信された信号を用いてそれぞれＩＳＡＲの動画像を再生することができる。

図３は、各アンテナで受信された信号を用いて得られるＩＳＡＲ動画像の概念図である。図３において、２３、２４、２５は送受信アンテナ３、受信用補助アンテナ４、受信用補助アンテナ５の受信信号を用いてそれぞれ再生されたＩＳＡＲ動画像である。ＩＳＡＲ動画像の１枚目のフレームは、１パルス目からＨパルス目のデータを用いて生成され、２フレーム目は１＋Ｊパルス目からＨ＋Ｊパルス目を用いて生成されるものとする。３フレーム目以降も同様で、結局、ｎフレーム目は１＋（ｎ−１）Ｊパルス目からＨ＋（ｎ−１）Ｊパルス目のデータを用いて生成される。このとき、動画像のフレーム間の時間間隔ｔ_ｆｒは次式で表される。

また、目標の運動および３次元形状を推定するのに用いる動画像のフレーム数をＮとすると、全体の観測に必要な送信パルス数はＨ＋（Ｎ−１）Ｊパルスであり、観測時間は｛Ｈ＋（Ｎ−１）Ｊ｝ｔ_ｐｒｉである。

上記の観測およびＩＳＡＲ処理の結果、動画像２３、２４、２５のｎ番目のフレームにおいて、孤立反射点Ｐ_ｋに対応する輝点が現れるものとする。すると、このフレームにおける輝点の座標の読み取り値から孤立反射点Ｐ_ｋのレーダの視線方向の位置ｒ_ｋｎとドップラー周波数ｆ_ｄｋｎ、および輝点のピーク複素振幅ｓ_０ｋｎ、ｓ_１ｋｎ、ｓ_２ｋｎが観測される（図３参照）。ここで、厳密には各アンテナから孤立反射点Ｐ_ｋまでの距離や、各アンテナで観測されるドップラー周波数は異なるが、式（５）の条件より、ｒ_ｋｎとｆ_ｄｋｎは送受信アンテナ３における観測量で代表できるものとする。

以下では、これらの観測量と、求めたい未知パラメータ

および

との関係を導く。まず、ｒ_ｋｎとｆ_ｄｋｎは次式の関係を満たす。

ただし、ＩＳＡＲ画像上におけるレンジ方向の原点の位置は、空間内に定義した原点Ｏであるとして定義した。なお、孤立反射点の位置

は当然時変であるが、以下ではｎ番目のフレームの観測時（より具体的には、

パルス目の送受信時；ただし，

は小数点以下の切り捨てを表す）における位置を

と考える。

次に、ピーク複素振幅ｓ_０ｋｎ、ｓ_１ｋｎ、ｓ_２ｋｎについては、これらの観測値の間の位相差と

が次式の関係を満たす。

なお、式（９）と式（１０）は、アレイアンテナを用いて遠方からの到来波を受信する際に、その到来角度によって決まるアレイ素子間の位相差の式に他ならない。目標の回転に伴い、

が変化するため、φ_１ｋｎ、φ_２ｋｎは変化するが、短い時間の間ではその変化はほぼ線形とみなすことができ、その変化率は次式で表される。

ここで、式（１１）、式（１２）の導出には以下の関係および近似を用いた。

なお、

のとき、送受信アンテナ３から各孤立反射点への方向はほぼ平行であるため、次式の近似が成り立つ。

また、位相差の変化率

は、ｎ番目のフレームの前後のＮフレームにおける位相差の値

を用い、式（１５）、（１６）で表される最小二乗法による直線あてはめによって推定できる。

なお、ここで、

は小数点以下の切り捨て、

は小数点以下の切り上げをそれぞれ表す。

図４は、観測された位相差の値に直線を当てはめて式（１５）、（１６）に示す位相差の変化率を求める概念図である。

原理的には、上記で導いた式（１）、（７）、（８）、（１１）、（１２）の関係を用いて、目標の運動および三次元形状を推定することが可能である。以下では、図１を用いて、この発明の実施の形態１に係る画像レーダ装置の機能構成および動作について説明する。

図１に示す画像レーダ装置は、送信機１、送受切換器２、送受信アンテナ３、受信用補助アンテナ４及び５、受信機６〜８、レーダ画像生成部９、画像履歴蓄積部１０、目標運動・形状推定部１１及び表示部１８を備えている。ここで、目標運動・形状推定部１１は、反射点抽出部１２、フレーム間反射点対応部１３、位相差算出部１４、位相差変異算出部１５、反射点速度算出部１６、目標形状・角速度算出部１７を有し、目標形状・角速度算出部１７で算出された結果は表示部１８で表示される。

送信機１は、周知の発振器および変調器などを含み、送受切換器２および受信機６〜８とともに、送受信制御部（図示せず）の制御下で駆動される。送受信アンテナ３は、目標に向けて送信パルスを出射するとともに、送信パルスが目標で反射されて戻ってくる反射電波を受信電波として受信する。このとき、受信用補助アンテナ４と５も同時に該反射電波を受信電波として受信する。

受信機６〜８は、それぞれ送受信アンテナ３および受信用補助アンテナ４、５で受信された電波を増幅、検波して得られた受信信号をレーダ画像再生部９に送る。そして、前述の通り、本実施の形態１に係る画像レーダ装置においては、以上のようなパルスの送受信を繰返し周期ｔ_ｐｒｉで繰り返す。

レーダ画像再生部９では、前述の通り、それぞれ送受信アンテナ３、受信用補助アンテナ４、受信用補助アンテナ５の受信信号を用いてＩＳＡＲ動画像２３、２４、２５を再生する。ここで、ｎフレーム目を再生するためには１＋（ｎ−１）Ｊパルス目からＨ＋（ｎ−１）Ｊパルス目のデータを用いる。

画像履歴蓄積部１０では、レーダ画像再生部９で再生された動画像２３、２４、２５をメモリに蓄積する。目標運動・形状推定部１１は、画像履歴蓄積部１０に一旦格納された動画像２３、２４、２５の情報を利用して、目標の運動と形状を推定して、推定結果の形状を表示部１８に出力する。

以下では、目標運動・形状推定部１１の構成と動作を説明する。まず、反射点抽出部１２では、動画像の各フレームにおいて、輝点（輝度の強い点）を検出し、その位置と複素振幅を計測する。輝点の検出には一般に良く知られたＣＦＡＲ（Constant False Alarm Rate）処理などを用いればよい。以下では、反射点抽出部１２において、動画像の各フレームにおいてＫ点ずつの点目標が検出されたと仮定して説明を進める。反射点抽出部１２は、検出処理の結果として、Ｋ点の点目標について、レーダの視線方向の位置ｒ_ｋｎとドップラー周波数ｆ_ｄｋｎ、および輝点のピーク複素振幅ｓ_０ｋｎ、ｓ_１ｋｎ、ｓ_２ｋｎを出力する（ｋ＝１，２，・・・，Ｋ）。

フレーム間反射点対応部１３では、各フレームで検出された輝点同士を対応付ける。ここでは、レーダの目標追尾などに利用される追尾方式などを利用して対応付けを行う。反射点抽出部１２の出力結果において、ｎ＋１フレーム目におけるｋ番目の輝点の信号と、ｎフレーム目におけるｋ番目の輝点の信号が、同じ孤立反射点からの信号に対応しているとは限らない。フレーム間反射点対応部１３における対応付け処理（追尾処理）によって、フレーム間での輝点の対応付けがなされ、その結果、同じ孤立反射点Ｐ_ｋの信号が全てのフレームにおいて、ｋ番目の信号としてラベル付けされる。

なお、反射点抽出部１２において、全てのフレームにおいて点目標が必ずＫ点検出されるとは限らない。例えば、ｎ番目のフレームにおいてＫ点検出されていて、ｎ＋１番目のフレームにおいてＫ−１点しか検出されなければ、ｎ番目で検出された輝点のうち一つは、ｎ＋１番目のフレームで検出された輝点と対応付けられないことになる。このような場合は、フレーム間反射点対応部１３において対応付けできない旨を判定して記録する。

次いで、位相差算出部１４では、送受信アンテナ３と、受信用補助アンテナ４，５で得られた信号の位相差を、各フレーム、各検出信号に対して次式によって算出する。

次に、位相差変化率算出部１５は、各フレーム、各検出信号に対して次式に示す位相差の変化率を算出する。

なお、位相差変化率算出部１５における上記の最小二乗法による変化率の算出の代わりに、各輝点の反射強度（反射電力や振幅など）で重み付けを行い、重み付け最小二乗法によって変化率を算出しても良い。重み付けをすることによって、ＳＮＲ（Signal to Noise Ratio）の高い信号の情報をより重視することが出来、位相差変化率の推定精度を向上することが可能である。

次に、反射点速度算出部１６は、各孤立反射点の速度

を求める。反射点速度算出部１６の処理を説明するために、まず、式（８）、（１１）、（１２）を整理してまとめると次式を得る。

ただし、上付きのＴは行列の転置を表す。ここで、行列

を構成する要素は、レーダの視線方向ベクトル

と基線ベクトル

であり、これらは既知であると考えることができる。また、左辺を構成する要素はレンジｒ_ｋｎ、ドップラー周波数ｆ_ｄｋｎと位相差の変化率

であるが、レンジｒ_ｋｎ、ドップラー周波数ｆ_ｄｋｎは、反射点抽出部１２の出力（ただし、フレーム間反射点対応部１３によってラベル付けされたもの）として得られており、位相差の変化率

は位相差変化率算出部１５の出力として得られている。さらに、原点までの距離ｒ_０と波長λも含まれるが、これらは既知の値として扱うことができる。

以上より、反射点速度算出部１６は、各孤立反射点の速度

を次式によって求める。

目標形状・角速度算出部１７は、角速度

および各孤立反射点の位置

を算出する。以下に、式（１）、（７）を再度示す。

これらの方程式の左辺において、

は式（１８）で得られており、ｒ_ｋｎは観測値である。また、これらの方程式はＫ個の孤立反射点についてそれぞれ成り立つので、得られる方程式の数は３Ｋ＋Ｋ本である。一方、未知のパラメータは

と

であり、未知変数の数は３＋３Ｋ個である。したがって、連立方程式が解を持つための必要条件は、次式で表される。

すなわち、目標上に孤立反射点が３点以上存在して、それが観測できることが条件である。一般に、レーダの観測対象となる目標を構成する孤立反射点の数は３つ以上であるから、この条件は概して成立すると考えてよい。

以上を踏まえ、目標形状・角速度算出部１７は、式（１９）、（２０）の連立方程式を解くことによって、角速度

および各孤立反射点の位置

を算出する。なお、Ｋ＞３の場合は最小二乗法によって解を求める。また、ここでも、位相差変化率算出部１５の場合と同様に、各輝点の強度によって重み付けをした重み付け最小二乗法を用いて解を求めても良い。重み付けをすることによって、ＳＮＲの高い信号の情報を重視するために、角速度の推定精度および目標形状の推定精度を向上することが可能である。

最後に、目標形状・角速度算出部１７において算出された孤立反射点の位置

は表示部１８に送られ、表示部１８は推定された目標の三次元形状を表示する。

ところで、反射点速度算出部１６の処理において、

の逆行列が存在する条件は、次式で表される。

すなわち、２つの基線ベクトルが互いに平行でなく、かつ、レーダの視線ベクトルが２つの基線ベクトルの張る平面に直交する成分を持っていることが条件となる。したがって、受信用補助アンテナ４、５は上記の条件を満たすように配置する必要がある。なお、受信用補助アンテナ４、５は２つの基線ベクトルが直交するように配置するのが望ましい。

また、目標形状・角速度算出部１７において、連立方程式（１９）、（２０）が線形独立である条件は、次式で表される。

つまり、目標の回転軸がレーダの視線方向と一致しないこと、および回転軸がレーダの視線方向と直交しないことを示している。この条件が成立するためには、目標の運動が既知である必要があるため、条件によっては目標の三次元形状が推定できない場合があることを示唆する。この問題に対する対策は実施の形態３に述べる。

以上のように、この発明の実施の形態１によれば、送受信アンテナ３に加えて受信用補助アンテナ４、５を具備し、各アンテナに接続された受信機６〜８で受信した信号を用いてそれぞれＩＳＡＲ動画像を生成することにより、これらのＩＳＡＲ動画像から観測された視線方向の位置およびドップラー周波数に加えて、観測信号の位相差の情報を利用できるように構成したので、レーダの観測信号を用いることで、目標の運動と三次元形状を同時に推定することが可能となる効果を奏する。

なお、以上においては、受信用補助アンテナの数は２つに限定して説明していたが、２つ以上ある場合も、同様の原理によって、実施が可能である。

実施の形態２．
図５は、この発明の実施の形態２に係る画像レーダ装置の機能構成を示すブロック図である。図５において、新たな符号として、２６は相関処理部、１１ａは相関動画像利用目標運動・形状推定部、２７は相関画像履歴蓄積部、２８は相関画像反射点抽出部であり、その他は図１と同様である。

ＩＳＡＲ動画像２３、２４、２５の各フレームの画素数をレンジ軸方向にＩ個、ドップラー軸方向にＪ個とする。動画像２３、２４、２５の（ｉ，ｊ）における画素の観測値をそれぞれｓ_０ｎ（ｉ，ｊ）、ｓ_１ｎ（ｉ，ｊ）、ｓ_２ｎ（ｉ，ｊ）とすると、相関処理部２６では、動画像２３と動画像２４の相関および動画像２３と動画像２５の相関を次式によって算出する。

ここで、算出されたｃ_１０ｎ（ｉ，ｊ）とｃ_２０ｎ（ｉ，ｊ）を、以下では相関動画像と呼ぶ。式（２１）、（２２）では、簡単のため、画素単位で複素共役を掛け合わせる処理のみを記してあるが、相関処理部２６では、相関を計算する画像のうち一方をサブピクセル単位でずらしながら、画素毎の相関値の二乗和が最大となるような計算をして、２つの画像の位置あわせを実施するように構成しても良い。このようにすることによって、送受信アンテナ３および受信用補助アンテナ４，５で観測された動画像２３、２４、２５の位置あわせを実施することが出来る。相関処理部２６によって算出されたこれらの相関動画像は、相関画像履歴蓄積部２７に送られてメモリに一旦格納される。

相関画像反射点抽出部２８は、相関動画像の中から、振幅あるいは電力値の大きい輝点を検出する。ここでの処理は、反射点抽出部１２と同様にＣＦＡＲ処理などによって実現できる。相関画像反射点抽出部２８は検出処理の結果として、Ｋ点の点目標について、レーダの視線方向の位置ｒ_ｋｎとドップラー周波数ｆ_ｄｋｎ、および輝点のピーク位置における位相差φ_１ｋｎ，φ_２ｋｎを出力する（ｋ＝１，２，・・・，Ｋ）。

以下の動作は実施の形態１と同様であるが、相関画像反射点抽出部２８においてすでに位相差が算出されているので、実施の形態１で必要であった位相差算出部１４は不要となる。

以上のように、この発明の実施の形態２によれば、送受信アンテナ３と受信用補助アンテナ４，５に接続された受信機で受信した信号を用いてそれぞれＩＳＡＲ動画像を生成し、最初にこれらのＩＳＡＲ動画像のフレーム毎に相関動画像を算出する際にサブピクセル単位で画像間の位置あわせを実施するため、位相差の推定精度が向上し、ひいては目標運動及び形状の推定精度を向上できる効果を奏する。

実施の形態３．
図６は、この発明の実施の形態３に係る画像レーダ装置の機能構成を示すブロック図である。図６において、新たな符号として、２９は形状推定可能性判定部であり、その他は図１と同様である。

式（２３）に示したとおり、目標の回転軸がレーダの視線方向と一致する場合、または、回転軸がレーダの視線方向と直交する場合、目標形状・角速度推定部１７においては、連立方程式の線形独立性がくずれるため、孤立反射点の位置ベクトル

及び角速度

を算出することができない。そこで、形状推定可能性判定部２９では、反射点速度算出部１６で算出された孤立反射点の速度ベクトル

を用いて、まず、回転軸の向きを推定する。

図２のジオメトリからも明らかなように、孤立反射点の速度ベクトル

は回転軸と直交する関係にある。したがって、形状推定可能性判定部２９においては、反射点速度算出部１６で算出された孤立反射点の速度ベクトル

の全てに直交するベクトルを算出することで、回転軸の向きを算出する。ただし、一般的には、速度ベクトル

の推定量には誤差分が含まれているため、速度ベクトル

の全てに直交するベクトルは存在しない場合がある。

このような場合にも対応するため、例えば、形状推定可能性判定部２９では、以下のような方法で回転軸を推定する。まず、次式によって、速度ベクトル

の共分散行列

を算出する。

次に、この共分散行列

の固有値分解を実施し、最小固有値に対応する固有ベクトル

を得る。この最小固有値に対応する固有ベクトルの向きを回転軸とする。なお、反射点速度算出部１６で算出された孤立反射点の速度ベクトル

が誤差を持たなければ、最小固有値の値は０となる。

次いで、最小固有値に対応する固有ベクトル

と視線方向ベクトル

を比較する。この２つのベクトルが、平行でもなく、直交もしていなければ、目標形状・角速度算出部１７において、孤立反射点の位置ベクトル

及び角速度

を算出する。一方、この２つのベクトルが平行であるか、あるいは直交している場合については、表示部１８にその旨を出力し、表示部１８においては、現在の位置関係では目標形状を正しく推定できないことをオペレータに伝える。画像レーダ装置が例えばヘリコプタなどの移動体に装備されている場合、オペレータは、この情報に基づいて画像レーダ装置を搭載した移動体の位置を移動させ、再度観測をやり直すことにより、目標の運動及び三次元形状を推定する。

以上のように、この発明の実施の形態３によれば、形状推定可能性判定部２９を備えるように構成したので、目標の孤立反射点の速度を推定した後に、形状推定の可能性を判定し、可能であれば形状推定を実施し、不可能であれば、プラットフォームを移動することで、形状推定を可能とする効果を奏する。

実施の形態４．
図７は、この発明の実施の形態４に係る画像レーダ装置の機能構成を示すブロック図である。図７において、新たな符号として、３０は形状合成部であり、その他は図１と同様である。

実施の形態１乃至３においては、ＩＳＡＲ動画像のうちＮ個のフレームの情報を用いて目標運動及び形状を推定していた。しかし、ＩＳＡＲ動画像は連続的に観測し続けることが可能であるため、連続するＮ個のフレームを用いて、目標形状を次々と繰り返し推定することが可能である。形状合成部３０では、これらの連続して推定された形状を合成して、より精度のよい形状推定結果を得る。

形状合成部３０は、例えば、連続して得られた推定結果を単に重畳するように構成することが出来る。特に、目標形状が複雑で、一方向からの観測だけでは遮蔽の影響により全体の形状が推定できないような場合、長時間の観測により複数の方向から観測された信号を元に推定された結果を重畳して表示することができ、より正確な三次元形状を得ることが可能となる。

または、形状合成部３０は、個々の推定結果の平均処理をするような形で、推定精度を向上するように構成しても良い。

以上のように、この発明の実施の形態４によれば、形状合成部３０を備えるように構成したので、複数回の形状推定結果を合成することができる。そのため、形状推定精度が向上したり、複数の角度から観測された形状を重畳することで、遮蔽の影響を低減する効果を奏する。

実施の形態５．
図８は、この発明の実施の形態５に係る画像レーダ装置の機能構成を示すブロック図である。図８において、新たな符号として、１１ｂはドップラー周波数変化率利用型目標運動・形状推定部、１７ａはレンジ方向加速度利用型目標形状・角速度算出部、３１はドップラー周波数変化率算出部、３２は反射点レンジ方向加速度算出部であり、その他は図１と同様である。

また、図９と図１０は、この発明の実施の形態５に係る画像レーダ装置の動作を説明するための説明図である。

実施の形態１乃至４の目標運動・形状推定部１１および相関動画像利用目標運動・形状推定部１１ａの処理において、ドップラー周波数についてはｎ番目のフレームにおけるドップラー周波数ｆ_ｄｋｎの情報のみを利用していた。しかし、一般に孤立反射点のドップラー周波数の値は、図９に示すように時間と共に変化する。したがって、ドップラー周波数の変化率の情報も用いたほうが目標形状の推定精度を向上できることが期待される。そこで、ドップラー周波数変化率利用型目標運動・形状推定部１１ｂでは、目標形状および角速度の推定に、ドップラー周波数の変化率

の情報も用いる。

以下では、まず、レンジ方向加速度利用型目標形状・角速度算出部１７ａの動作を説明するため、ドップラー周波数の変化率

が得られたと仮定した場合の形状推定方式について説明し、ドップラー周波数の変化率

の推定方式については後述する。

ドップラー周波数の変化率

を観測すると、孤立反射点Ｐ_ｋのレンジ方向の速度の変化率

が次式によって求められる。

一方で、レンジ方向の速度の変化率

と孤立反射点の位置

は、次式の関係を満たす。

なお、ここでは、実施の形態１乃至４と同様に目標が等角速度運動をしていることを前提とし、その角速度ベクトルの変化率はゼロとした。

このように、ドップラー周波数の変化率

から求められたレンジ方向の速度の変化率

は、孤立反射点の位置

と式（２４）の関係を満たすので、この関係を目標形状の推定に利用することが出来る。具体的には、式（１９）、（２０）および式（２４）の連立方程式を解くことによって、各孤立反射点の位置

を算出できる。なお、式（１９）、（２０）および（２４）の連立方程式をベクトル表現すると次式で表される。

ここで、上付きのＴはベクトルおよび行列の転置を表す。また、

であり、

としたとき、行列

は次式によって定義される。

なお、式（２６）において、各孤立反射点の速度については、ｎ番目のフレーム観測時における速度であることを明示するため、

と表記している。

以上を踏まえ、レンジ方向加速度利用型目標形状・角速度算出部１７ａは、式（２６）の連立方程式を非線形の最小二乗法で解くことによって、角速度

および各孤立反射点の位置

を算出する。

式（２６）の解法については如何なる方法を用いてもかまわないが、角速度

について、その向きを表す単位ベクトル

と角速度の大きさを表す

を分離して推定することによって、式（２６）をより効率良く解くことが可能となる。

図１０は、この考え方に基づいてレンジ方向加速度利用型目標形状・角速度算出部１７ａを構成した場合の構成図である。図１０において、新たな符号として、１７ａ−１は回転軸推定部、１７ａ−２はレンジ方向加速度利用型角速度・反射点座標算出部であり、その他は図８と同様である。

回転軸推定部１７ａ−１においては、実施の形態３に示した形状推定可能性判定部２９と同様に、反射点速度算出部１６で算出された孤立反射点の速度ベクトル

を用いて、回転軸の向きを推定する。すなわち、まず、次式によって、速度ベクトル

の共分散行列

を算出する。

次に、この共分散行列

の固有値分解を実施し、最小固有値に対応する固有ベクトル

を得る。この最小固有値に対応する固有ベクトルの向きを回転軸とみなし、角速度

の向きを表す単位ベクトル

の推定値を

で与える。

次いで、レンジ方向加速度利用型角速度・反射点座標算出部１７ａ−２においては、

の値を一つ仮定して、行列

を与えた上で、式（２６）を次の線形最小二乗法によって解く。

このとき、孤立反射点の座標

の推定値は仮定した

の値の関数であるから、

と表記している。また、式（２９）の最小二乗法による解における二乗誤差は次式で表される。

次いで、レンジ方向加速度利用型角速度・反射点座標算出部１７ａ−２においては、

の値をとりうる値の範囲で様々に仮定して、式（２９）を解く。その上で、二乗誤差

を最小にする

と、そのときの孤立反射点の座標の推定値

を最終的な推定値として出力する。このように構成することにより、解の探索は回転の速さ

についてのみ実施すればよく、その他は線形最小二乗法によって求めることが可能である。

以上においては、ドップラー周波数の変化率

が得られたものと仮定して議論を進めたが、ドップラー周波数の変化率

は、実施の形態１における位相差変化率算出部１５の動作と同様に、ｎ番目のフレームの前後のＮフレームにおけるドップラー周波数の値

から推定することが可能である。すなわち、ドップラー周波数変化率算出部３１は、次式によってドップラー周波数の変化率

を推定する。これはすなわち、図９に示すようにｎ番目のフレームの前後のＮフレームにおけるドップラー周波数の観測値に当てはめた直線の傾きに相当する。

ここで、

はドップラー周波数の変化率の推定値を表す。なお、ドップラー周波数変化率算出部３１における上記の最小二乗法による変化率の算出の代わりに、各輝点の反射強度（反射電力や振幅など）で重み付けを行い、重み付け最小二乗法によって変化率を算出しても良い。重み付けをすることによって、ＳＮＲ（Signal to Noise Ratio）の高い信号の情報をより重視することが出来、ドップラー周波数変化率の推定精度を向上することが可能である。

ドップラー周波数変化率算出部３１においてドップラー周波数の変化率の推定値

が算出された後は、反射点レンジ方向加速度算出部３２において式（２３）によりレンジ方向の加速度が算出される。

以上のように、この発明の実施の形態５によれば、ドップラー周波数変化率算出部３１および反射点レンジ方向加速度算出部３２を備え、さらに、レンジ方向加速度利用型目標形状・角速度算出部１７ａを備えるように構成したので、孤立反射点のレンジ方向の位置および孤立反射点の速度に加えて、レンジ方向の加速度の情報も用いることができるようになり、目標形状の推定精度を向上する効果を奏する。

また、レンジ方向加速度利用型目標形状・角速度算出部１７ａが、回転軸推定部１７ａ−１と、レンジ方向加速度利用型角速度・反射点座標算出部１７ａ−２を備えるように構成したので、目標の角速度と孤立反射点の座標の推定に用いる非線形最小二乗問題を効率良く解くことができる効果を奏する。

実施の形態６．
図１１は、この発明の実施の形態６に係る画像レーダ装置の機能構成を示すブロック図である。図１１において、新たな符号として、１１ｃは位相差２次微分利用型目標運動・形状推定部、１５ａは位相差変化率・２次微分算出部、１７ｂは加速度利用型目標形状・角速度算出部、３２ａは反射点加速度算出部であり、その他は図８と同様である。

また、図１２と図１３は、この発明の実施の形態６に係る画像レーダ装置の動作を説明するための説明図である。

実施の形態５のレンジ方向加速度利用型目標形状・角速度算出部１７ａでは、目標形状および角速度の推定に、ドップラー周波数の変化率の推定値から求められるレンジ方向の加速度の情報を利用している。しかし、一般に孤立反射点の加速度はレンジ方向のみならず、レンジに直交する成分も有している。したがって、加速度のレンジ方向に直交する成分を用いることができれば、目標形状の推定精度を向上できると期待される。

本実施の形態６に係る画像レーダ装置においては、位相差の２次微分成分の推定値から加速度のレンジ方向に直交する成分を推定する。まず、位相差変化率・２次微分算出部１５ａにおいては、図１２に示すように位相差の観測値に最小二乗法によって２次の曲線を当てはめることによって、位相差の変化率と２次微分成分を算出する。すなわち、ｎフレーム目における位相差、位相差の変化率、位相差の２次微分を次式によって推定する。

ここで、

はそれぞれ、送受信アンテナ３と受信用補助アンテナ４で得られた信号のｎ番目のフレームにおける位相差、位相差の変化率、位相差の２次微分の推定値を表す。また、

はそれぞれ、送受信アンテナ３と受信用補助アンテナ５で得られた信号のｎ番目のフレームにおける位相差、位相差の変化率、位相差の２次微分の推定値を表す。位相差変化率・２次微分算出部１５ａは、以上によって算出された

を出力する。

次いで、反射点加速度算出部３２ａは、ドップラー周波数変化率算出部３１によって算出されたドップラー周波数の変化率の推定値

と位相差変化率・２次微分算出部１５ａにおいて算出された位相差、位相差の変化率、位相差の２次微分の推定値

を利用して、孤立反射点の加速度を算出する。

ここで、位相差の２次微分と加速度の関係は次式で表される。

ただし、式(３４)、式(３５)の導出には以下の関係および近似を用いた。

観測対象の目標のサイズに対して十分に遠方から観測していることを想定しているため、

としている。なお、

のとき、送受信アンテナ３から各孤立反射点への方向はほぼ平行であるため、次式の近似が成り立つ。

反射点加速度算出部３２ａは、式（３４）、（３５）および（２３）より、反射点の加速度を算出する。まず、式（３４）、（３５）および（２３）を整理してベクトル表現で表すと次式で表される。

したがって、加速度

は次式によって算出できる。

反射点加速度算出部３２ａは、式（３９）によって算出された加速度

を出力する。

次いで、加速度利用型目標形状・角速度算出部１７ｂは、実施の形態５におけるレンジ方向加速度利用型目標形状・角速度算出部１７ａと同様の操作により、角速度

および各孤立反射点の位置

を算出する。

まず、加速度

と孤立反射点の位置

の関係を次式に示す。

なお、ここでは、実施の形態５と同様に目標が等角速度運動をしていることを前提とし、その角速度ベクトルの変化率はゼロとした。

このように、加速度

は、孤立反射点の位置

と式（４０）の関係を満たすので、この関係を目標形状の推定に利用することが出来る。具体的には、式（１９）、（２０）および（４０）の連立方程式を解くことによって、各孤立反射点の位置

を算出できる。なお、式（１９）、（２０）および（４０）の連立方程式をベクトル表現すると次式で表される。

ここで、行列

は式（２７）によって定義されるものである。

以上を踏まえ、加速度利用型目標形状・角速度算出部１７ｂは、式（４１）の連立方程式を非線形の最小二乗法で解くことによって、角速度

および各孤立反射点の位置

を算出する。なお、実施の形態５で示したレンジ方向加速度利用型目標形状・角速度算出部１７ａと同様に、角速度

について、その向きを表す単位ベクトル

と回転の速さを表す

を分離して推定することによって、式（４１）をより効率良く解くことが可能である。

図１３は、この考え方に基づいて加速度利用型目標形状・角速度算出部１７ｂを構成した場合の構成図である。図１３において、新たな符号として、１７ｂ−１は回転軸推定部、１７ｂ−２は加速度利用型角速度・反射点座標算出部であり、その他は図１１と同様である。なお、回転軸推定部１７ｂ−１と加速度利用型角速度・反射点座標算出部１７ｂ−２の動作は、それぞれ、実施の形態５における回転軸推定部１７ａ−１とレンジ方向加速度利用型角速度・反射点座標算出部１７ａ−２と同様であり、説明を要さない。

以上のように、この発明の実施の形態６によれば、位相差変化率・２次微分算出部１５ａおよび反射点加速度算出部３２ａを備え、さらに、加速度利用型目標形状・角速度算出部１７ｂを備えるように構成したので、孤立反射点のレンジ方向の位置および孤立反射点の速度に加えて加速度の情報も用いることができるようになり、目標形状の推定精度をさらに向上する効果を奏する。

加えて、加速度利用型目標形状・角速度算出部１７ｂが、回転軸推定部１７ｂ−１と、加速度利用型角速度・反射点座標算出部１７ｂ−２を備えるように構成したので、目標の角速度と孤立反射点の座標の推定に用いる非線形最小二乗問題を効率良く解くことができる効果を奏する。

実施の形態７．
図１４は、この発明の実施の形態７に係る画像レーダ装置の機能構成を示すブロック図である。図１４において、新たな符号として、１１ｄは時変角速度考慮型目標運動・形状推定部、１７ｃは時変角速度考慮型目標形状・角速度算出部であり、その他は図８と同様である。

また、図１５は、この発明の実施の形態７に係る画像レーダ装置の動作を説明するための説明図である。

実施の形態１乃至６に係る画像レーダ装置においては、目標の運動が等角速度運動であることを前提としていた。目標の運動の変動に対して十分に短い時間の間であれば、この前提は成立するが、フレーム数を増やして動画像の観測時間を増加すると等角速度運動の前提が成立しない可能性がある。そこで、本実施の形態７に係る画像レーダ装置においては、角速度が時間の関数であることを想定して目標運動および形状を推定する。ただし、ここでは角速度の時間変化は次式で示す１次式で表されることを前提とする。

ここで、

は角速度の変化率であり、

はｎ番目のフレーム観測時における角速度を表す。

角速度の時間変化を考慮する場合、加速度

と孤立反射点の位置

の関係は、次式で表される。

このように、加速度

は、孤立反射点の位置

と式（４３）の関係を満たすので、この関係を目標形状の推定に利用することが出来る。具体的には、式（１９）、（２０）および（４３）の連立方程式を解くことによって、各孤立反射点の位置

を算出できる。なお、式（１９）、（２０）および（４３）の連立方程式をベクトル表現すると次式で表される。

ここで、行列

は式（２７）によって定義されるものである（ただし、

で置き換える）。また、

であり、

としたとき、行列

は次式によって定義される。

以上を踏まえ、時変角速度考慮型目標形状・角速度算出部１７ｃは、式（４４）の連立方程式を非線形の最小二乗法で解くことによって、n番目のフレーム観測時における角速度

、角速度の変化率

、および各孤立反射点の位置

を算出する。

なお、実施の形態５で示したレンジ方向加速度利用型目標形状・角速度算出部１７ａと同様に、角速度

について、その向きを表す単位ベクトル

と角速度の大きさを表す

を分離して推定する。また、角速度の変化率

が、次の関係を満たすことに着目する。

ここで、表記の都合上、

と定義した。式（４６）より、角速度の変化率

については、向きを表す単位ベクトル

の変化率

と、角速度の大きさの変化率

を分離して推定できることが分かる。

図１５は、上記の考え方に基づいて時変角速度考慮型目標形状・角速度算出部１７ｃを構成した場合の構成図である。図１５において、新たな符号として、１７ｃ−１は複数回転軸推定部、１７ｃ−２は平均回転軸・回転軸変化率推定部、１７ｃ−３は時変角速度考慮型角速度・反射点座標算出部であり、その他は図１４と同様である。

まず、複数回転軸推定部１７ｃ−１は、実施の形態５に示した回転軸推定部１７ａ−１と同様に、反射点速度算出部１６で算出された孤立反射点の速度ベクトル

を用いて、回転軸の向きを複数回推定する。実施の形態５に示した回転軸推定部１７ａ−１においては、ｎ番目のフレームを中心とする近傍Ｎフレームの観測値を用いて推定した速度ベクトル

を用いて、ｎ番目のフレームの観測時における回転軸を推定していたが、本実施の形態における複数回転軸推定部１７ｃ−１は、ｎ番目のフレーム以外にもｎ＋ｉ

番目のフレームを中心としたＮフレームの観測値を用いて、ｎ＋ｉ番目のフレームの観測時における回転軸もそれぞれ求めるようにする。以下では、このようにして算出された回転軸の推定値を

と表す。

次いで、平均回転軸・回転軸変化率推定部１７ｃ−２は、

の値を用いて、ｎフレーム目における回転軸の向きとその変化率を次式によって推定する。

ここで、

はそれぞれ、ｎ番目のフレームにおける回転軸の向きとその変化率の推定値を表す。平均回転軸・回転軸変化率推定部１７ｃ−２は、上記のようにして推定された回転軸の向きとその変化率の推定値

を出力する。

次に、時変角速度考慮型角速度・反射点座標算出部１７ｃ−３においては、角速度の大きさ

の値とその変化率

の組を一つ仮定して、行列

を与えた上で、式（４４）を次の線形最小二乗法によって解く。

このとき、孤立反射点の座標

の推定値は仮定した角速度の大きさ

の値とその変化率

の値の関数であるから、

と表記している。また、式（４８）の最小二乗法による解における二乗誤差は次式で表される。

次いで、時変角速度考慮型角速度・反射点座標算出部１７ｃ−３においては、角速度の大きさ

の値とその変化率

の組を、とりうる値の範囲で様々に仮定して、式（４８）を解く。その上で、二乗誤差

を最小にする

の組と、そのときの孤立反射点の座標の推定値

を最終的な推定値として出力する。このように構成することにより、解の探索は

の組についてのみ実施すればよく、その他は線形最小二乗法によって求めることが可能である。

以上のように、この発明の実施の形態７によれば、時変角速度考慮型目標形状・角速度算出部１７ｃを備えるように構成したので、目標角速度が時変である場合についても、目標運動と目標形状を推定できる効果を奏する。このことは、観測時間を延長してより多数のフレーム数の動画を利用して目標形状を推定できることを示唆しており、目標形状の推定精度をさらに向上する効果をも奏する。

加えて、時変角速度考慮型目標形状・角速度算出部１７ｃが、複数回転軸推定部１７ｃ−１、平均回転軸・回転軸変化率推定部１７ｃ−２、時変角速度考慮型角速度・反射点座標算出部１７ｃ−３を備えるように構成したので、目標の運動と孤立反射点の座標の推定に用いる非線形最小二乗問題を効率良く解くことができる効果を奏する。

１ 送信機、２ 送受切換器、３ 送受信アンテナ、４，５ 受信用補助アンテナ、６〜８ 受信機、９ レーダ画像生成部、１０ 画像履歴蓄積部、１１ 目標運動・形状推定部、１１ａ 相関動画像利用目標運動・形状推定部、１１ｂ ドップラー周波数変化率利用型目標運動・形状推定部、１１ｃ 位相差２次微分利用型目標運動・形状推定部、１１ｄ 時変角速度考慮型目標運動・形状推定部、１２ 反射点抽出部、１３ フレーム間反射点対応部、１４ 位相差算出部、１５ 位相差変異算出部、１５ａ 位相差変化率・２次微分算出部、１６ 反射点速度算出部、１７ 目標形状・角速度算出部、１７ａ レンジ方向加速度利用型目標形状・角速度算出部、１７ａ−１ 回転軸推定部、１７ａ−２ レンジ方向加速度利用型角速度・反射点座標算出部、１７ｂ 加速度利用型目標形状・角速度算出部、１７ｂ−１ 回転軸推定部、１７ｂ−２ 加速度利用型角速度・反射点座標算出部、１７ｃ 時変角速度考慮型目標形状・角速度算出部、１７ｃ−１ 複数回転軸推定部、１７ｃ−２ 平均回転軸・回転軸変化率推定部、１７ｃ−３ 時変角速度考慮型角速度・反射点座標算出部、２６ 相関処理部、２７ 相関画像履歴蓄積部、２８ 相関画像反射点抽出部、２９ 形状推定可能性判定部、３０ 形状合成部、３１ ドップラー周波数変化率算出部、３２ 反射点レンジ方向加速度算出部、３２ａ 反射点加速度算出部。

Claims (20)

1. 目標に対して送信電波を照射する送信手段と、
   前記送信電波が前記目標で反射された反射電波を受信する、互いに異なる位置に配置された複数の受信用アンテナおよび受信機と、
   前記複数の受信用アンテナおよび受信機で受信された前記受信電波に基づき、前記目標のレーダ動画像を再生するレーダ画像再生部と、
   前記レーダ画像再生部で連続的に再生された複数のレーダ動画像を用いて目標の運動及び形状を推定する目標運動・形状推定手段と
   を備えた画像レーダ装置。
2. 請求項１に記載の画像レーダ装置において、
   前記目標運動・形状推定手段は、
   前記レーダ画像再生部において再生されたレーダ動画像の各フレームから、輝度の強い点の視線方向の位置とドップラー周波数および複素振幅の情報を抽出する反射点抽出手段と、
   各フレームにおいて前記反射点抽出手段で抽出された点をフレーム間で対応付けるフレーム間反射点対応手段と、
   異なる受信機で受信された信号に基づく動画像で観測された対応する点同士の前記複素振幅の位相差を算出する位相差算出手段と、
   前記位相差算出手段において算出された位相差の変化率を算出する位相差変化率算出手段と、
   前記反射点抽出手段で抽出された点のドップラー周波数および前記位相差変化率算出手段において算出された位相差の変化率を利用して、対応する反射点の速度を算出する反射点速度算出手段と、
   前記反射点抽出手段で抽出された点の視線方向の位置および前記反射点速度算出手段において算出された反射点の速度を利用して、目標の角速度および反射点の三次元の座標を算出する目標形状・角速度算出手段と
   を有することを特徴とする画像レーダ装置。
3. 請求項１に記載の画像レーダ装置において、
   前記レーダ画像再生部で連続的に再生された異なる受信機で受信された信号に基づき生成されたレーダ動画像間の画素ごとの相関を算出して相関動画像を得る相関処理部をさらに備え、
   前記目標形状・角速度算出手段は、前記レーダ画像再生部で連続的に再生された複数のレーダ動画像に対して、前記相関処理部で算出された相関動画像を用いて目標の運動及び形状を推定する相関動画像利用目標運動・形状推定手段でなる
   ことを特徴とする画像レーダ装置。
4. 請求項３に記載の画像レーダ装置において、
   前記相関動画像利用目標運動・形状推定手段は、
   相関動画像の各フレームから、相関値の大きい点の視線方向の位置とドップラー周波数および複素振幅の情報を抽出する反射点抽出手段と、
   各フレームにおいて前記反射点抽出手段で抽出された点をフレーム間で対応付けるフレーム間反射点対応手段と、
   前記相関動画像における位相値として得られている位相差の変化率を算出する位相差変化率算出手段と、
   前記反射点抽出手段で抽出された点のドップラー周波数および、前記位相差変化率算出手段において算出された位相差の変化率を利用して、対応する反射点の速度を算出する反射点速度算出手段と、
   前記反射点抽出手段で抽出された点の視線方向の位置および、前記反射点速度算出手段において算出された反射点の速度を利用して、目標の角速度および反射点の三次元の座標を算出する目標形状・角速度算出手段と
   を有することを特徴とする画像レーダ装置。
5. 請求項１から４までのいずれか１項に記載の画像レーダ装置において、
   前記受信用アンテナは、送受信を兼ねる送受信用アンテナを１つ、受信用補助アンテナを２つ有し、送受信用アンテナと受信用補助アンテナのうち一つをつなぐ２本の基線ベクトルが互いに直交するようにして、３つのアンテナが配置されている
   ことを特徴とする画像レーダ装置。
6. 請求項１から５までのいずれか１項に記載の画像レーダ装置において、
   前記反射点速度算出手段において算出された各反射点の速度から、目標運動の回転軸の向きを推定することにより、現在の観測のジオメトリの条件で目標の三次元形状を推定できるか否かを判定する形状推定可能性判定手段をさらに備えた
   ことを特徴とする画像レーダ装置。
7. 請求項６に記載の画像レーダ装置において、
   前記形状推定可能性判定手段は、前記反射点速度算出手段において算出された各反射点の速度ベクトルの共分散行列を算出し、前記共分散行列を固有値展開した後、その最小固有値に対応する固有ベクトルの向きを目標運動の回転軸の向きと推定する
   ことを特徴とする画像レーダ装置。
8. 請求項１から７までのいずれか１項に記載の画像レーダ装置において、
   前記目標運動・形状推定部から出力される複数回の観測の実施に基づく複数の目標形状の推定結果を合成する形状合成手段をさらに備えた
   ことを特徴とする画像レーダ装置。
9. 請求項８に記載の画像レーダ装置において、
   前記形状合成手段は、遮蔽による形状の欠損を補うために、異なる方向からの観測結果から推定された形状を重畳する
   ことを特徴とする画像レーダ装置。
10. 請求項８に記載の画像レーダ装置において、
    前記形状合成手段は、形状の推定精度を向上するために、複数の観測結果の平均処理を実施する
    ことを特徴とする画像レーダ装置。
11. 請求項１から１０までのいずれか１項に記載の画像レーダ装置において、
    前記位相差変化率算出手段は、複数フレームにおいて算出された位相差の値に対して最小二乗法によって直線あてはめを行うことによって、位相差変化率を算出する
    ことを特徴とする画像レーダ装置。
12. 請求項１から１０までのいずれか１項に記載の画像レーダ装置において、
    前記位相差変化率算出手段は、複数フレームにおいて算出された位相差の値に対して、前記反射点抽出手段で抽出された、対応する点の複素振幅から算出された反射強度の情報を用いて重み付けを行うことにより、重み付け最小二乗法によって直線あてはめを行うことによって位相差変化率を算出する
    ことを特徴とする画像レーダ装置。
13. 請求項１から１２までのいずれか１項に記載の画像レーダ装置において、
    前記目標形状・角速度算出部は、前記反射点抽出手段で抽出された３個以上の点のドップラー周波数および前記位相差変化率算出手段において算出された位相差の変化率を利用して、最小二乗法によって目標の角速度および反射点の三次元の座標を算出する
    ことを特徴とする画像レーダ装置。
14. 請求項１から１２までのいずれか１項に記載の画像レーダ装置において、
    前記目標形状・角速度算出部は、前記反射点抽出手段で抽出された３個以上の点のドップラー周波数および前記位相差変化率算出手段において算出された位相差の変化率を利用して、さらに、前記反射点抽出手段で抽出された、対応する点の複素振幅から算出された反射強度の情報を用いて重み付けを行うことにより、重み付け最小二乗法によって目標の角速度および反射点の三次元の座標を算出する
    ことを特徴とする画像レーダ装置。
15. 請求項１に記載の画像レーダ装置において、
    前記目標運動・形状推定手段は、目標形状および角速度の推定に、ドップラー周波数の変化率の情報を用いるドップラー周波数変化率利用型目標運動・形状推定手段でなり、
    前記ドップラー周波数変化率利用型目標運動・形状推定手段は、
    前記レーダ画像再生部において再生されたレーダ動画像の各フレームから、輝度の強い点の視線方向の位置とドップラー周波数および複素振幅の情報を抽出する反射点抽出手段と、
    各フレームにおいて前記反射点抽出手段で抽出された点をフレーム間で対応付けるフレーム間反射点対応手段と、
    異なる受信機で受信された信号に基づく動画像で観測された対応する点同士の前記複素振幅の位相差を算出する位相差算出手段と、
    前記位相差算出手段において算出された位相差の変化率を算出する位相差変化率算出手段と、
    前記反射点抽出手段で抽出された点のドップラー周波数および前記位相差変化率算出手段において算出された位相差の変化率を利用して、対応する反射点の速度を算出する反射点速度算出手段と、
    前記反射点抽出手段で抽出された点のドップラー周波数の変化率を算出するドップラー周波数変化率算出手段と、
    前記ドップラー周波数変化率算出手段で算出されたドップラー周波数の変化率を、レンジ方向の加速度に換算する反射点レンジ方向加速度算出手段と、
    前記反射点抽出手段で抽出された点の視線方向の位置および前記反射点速度算出手段において算出された反射点の速度および前記反射点レンジ方向加速度算出手段において算出された反射点のレンジ方向の加速度を利用して、目標の角速度および反射点の三次元の座標を算出するレンジ方向加速度利用型目標形状・角速度算出手段と
    を有することを特徴とする画像レーダ装置。
16. 請求項１５に記載の画像レーダ装置において、
    前記レンジ方向加速度利用型目標形状・角速度算出手段は、
    前記反射点速度算出部において算出された反射点の速度の集合の共分散行列を固有値分解して、最小固有値に対応する固有ベクトルを算出し、同固有ベクトルの向きを目標回転運動の回転軸方向の推定値として出力する回転軸推定手段と、
    前記回転軸推定手段によって算出された回転軸方向の推定値と、前記反射点抽出手段で抽出された点の視線方向の位置および前記反射点速度算出手段において算出された反射点の速度および前記反射点レンジ方向加速度算出手段において算出された反射点のレンジ方向の加速度とを利用して、角速度の大きさをパラメ−タとしつつ線形最小二乗法によって反射点の三次元の座標の推定値を算出することを繰り返し、推定誤差を最小とする角速度の大きさと、反射点の三次元の座標の推定値を出力するレンジ方向加速度利用型角速度・反射点座標算出部と
    を有することを特徴とする画像レーダ装置。
17. 請求項１に記載の画像レーダ装置において、
    前記目標運動・形状推定手段は、目標形状および角速度の推定に、位相差の２次微分成分の情報を用いる位相差２次微分利用型目標運動・形状推定手段でなり、
    前記位相差２次微分利用型目標運動・形状推定手段は、
    前記レーダ画像再生部において再生されたレーダ動画像の各フレームから、輝度の強い点の視線方向の位置とドップラー周波数および複素振幅の情報を抽出する反射点抽出手段と、
    各フレームにおいて前記反射点抽出手段で抽出された点をフレーム間で対応付けるフレーム間反射点対応手段と、
    異なる受信機で受信された信号に基づく動画像で観測された対応する点同士の前記複素振幅の位相差を算出する位相差算出手段と、
    前記位相差算出手段において算出された位相差の変化率および２次微分成分を算出する位相差変化率・２次微分算出手段と、
    前記反射点抽出手段で抽出された点のドップラー周波数および前記位相差変化率・２次微分算出手段において算出された位相差の変化率を利用して、対応する反射点の速度を算出する反射点速度算出手段と、
    前記反射点抽出手段で抽出された点のドップラー周波数の変化率を算出するドップラー周波数変化率算出手段と、
    前記ドップラー周波数変化率算出手段で算出されたドップラー周波数の変化率と、前記位相差変化率・２次微分算出手段において算出された位相差の２次微分成分を利用して反射点の加速度を推定する反射点加速度算出手段と、
    前記反射点抽出手段で抽出された点の視線方向の位置および前記反射点速度算出手段において算出された反射点の速度および前記反射点加速度算出手段において算出された反射点の加速度を利用して、目標の角速度および反射点の三次元の座標を算出する加速度利用型目標形状・角速度算出手段と
    を有することを特徴とする画像レーダ装置。
18. 請求項１７に記載の画像レーダ装置において、
    前記加速度利用型目標形状・角速度算出手段は、
    前記反射点速度算出部において算出された反射点の速度の集合の共分散行列を固有値分解して、最小固有値に対応する固有ベクトルを算出し、同固有ベクトルの向きを目標回転運動の回転軸方向の推定値として出力する回転軸推定手段と、
    前記回転軸推定手段によって算出された回転軸方向の推定値と、前記反射点抽出手段で抽出された点の視線方向の位置および前記反射点速度算出手段において算出された反射点の速度および前記反射点加速度算出手段において算出された反射点の加速度とを利用して、角速度の大きさをパラメ−タとしつつ線形最小二乗法によって反射点の三次元の座標の推定値を算出することを繰り返し、推定誤差を最小とする角速度の大きさと、反射点の三次元の座標の推定値を出力する加速度利用型角速度・反射点座標算出部と
    を有することを特徴とする画像レーダ装置。
19. 請求項１に記載の画像レーダ装置において、
    前記目標運動・形状推定手段は、目標形状および角速度の推定に、角速度が時間の関数であることを想定する時変角速度考慮型目標運動・形状推定手段でなり、
    前記時変角速度考慮型目標運動・形状推定手段は、
    前記レーダ画像再生部において再生されたレーダ動画像の各フレームから、輝度の強い点の視線方向の位置とドップラー周波数および複素振幅の情報を抽出する反射点抽出手段と、
    各フレームにおいて前記反射点抽出手段で抽出された点をフレーム間で対応付けるフレーム間反射点対応手段と、
    異なる受信機で受信された信号に基づく動画像で観測された対応する点同士の前記複素振幅の位相差を算出する位相差算出手段と、
    前記位相差算出手段において算出された位相差の変化率および２次微分成分を算出する位相差変化率・２次微分算出手段と、
    前記反射点抽出手段で抽出された点のドップラー周波数および前記位相差変化率・２次微分算出手段において算出された位相差の変化率を利用して、対応する反射点の速度を算出する反射点速度算出手段と、
    前記反射点抽出手段で抽出された点のドップラー周波数の変化率を算出するドップラー周波数変化率算出手段と、
    前記ドップラー周波数変化率算出手段で算出されたドップラー周波数の変化率と、前記位相差変化率・２次微分算出手段において算出された位相差の２次微分成分を利用して反射点の加速度を推定する反射点加速度算出手段と、
    前記反射点抽出手段で抽出された点の視線方向の位置および前記反射点速度算出手段において算出された反射点の速度および前記反射点加速度算出手段において算出された反射点の加速度を利用して、目標の角速度および反射点の三次元の座標を算出する時変角速度考慮型目標形状・角速度算出手段と
    を有することを特徴とする画像レーダ装置。
20. 請求項１９に記載の画像レーダ装置において、
    前記時変角速度考慮型目標形状・角速度算出手段は、
    前記反射点速度算出部において算出された反射点の速度の集合の共分散行列を固有値分解して、最小固有値に対応する固有ベクトルを算出し、同固有ベクトルの向きを目標回転運動の回転軸方向の推定値として出力することを複数回繰り返す複数回転軸推定手段と、
    前記複数回転軸推定手段によって算出された複数の回転軸方向の推定値から、中心時刻における回転軸方向および回転軸変化率を推定する平均回転軸・回転軸変化率推定手段と、
    前記反射点抽出手段で抽出された点の視線方向の位置および前記反射点速度算出手段において算出された反射点の速度および前記反射点加速度算出手段において算出された反射点の加速度と前記平均回転軸・回転軸変化率推定手段において算出された回転軸方向および回転軸変化率を利用して、角速度の大きさおよび角速度の大きさの変化率をパラメ−タとしつつ線形最小二乗法によって反射点の三次元の座標の推定値を算出することを繰り返し、このときの推定誤差を最小とする角速度の大きさと角速度の大きさの変化率と、そのときの反射点の三次元の座標の推定値を出力する時変角速度考慮型角速度・反射点座標算出部と
    を有することを特徴とする画像レーダ装置。

Images