JP6415288B2 - レーダ装置 - Google Patents

Description

この発明は、例えば、海面の波の速度や、上空の風の速度を推定するレーダ装置に関するものである。

レーダ装置である海洋レーダは、海の流れを面的に効率よく観測できる装置として知られており、短波帯の電波を海に向かって照射すると、視線方向に伝搬する波の成分の中で、電波の１／２となる波長と海面波の波長の成分が一致して共鳴することで、電波が後方に強く散乱する現象が生じる。この現象はブラッグ共鳴散乱と呼ばれている。
海面に反射された電波の反射波を受信し、その受信データを周波数解析することで、レンジビン（距離分解能）毎にドップラー（周波数）スペクトルを抽出し、最も強いピークとなる１次散乱のドップラーシフトを測定すれば、海面の波のレーダ視線方向の速度（以下、「流速」と称する）を推定することができる。

しかし、或るレンジビンに、船等の観測対象以外の物体が存在する場合、海流からの反射波よりも強い反射波（クラッタ、外部ノイズ等の呼び方が存在するが、ここでは不要波で統一する）を受信するため、観測対象以外の物体からの不要波のピーク周波数を誤って抽出して、正しい流速を推定できないことがある。
以下の特許文献１には、船等の観測対象以外の物体が存在する場合でも、正確な流速を推定することができるようにしているレーダ装置が開示されている。

特許文献１に開示されているレーダ装置では海流の空間分布特性を利用しており、このレーダ装置は、流れが同時観測されたビーム上の空間に対しては滑らかな分布条件が成立するため、レンジ方向の観測範囲を複数のブロックに分割して、ブロック毎に、流れの回帰曲線とのズレから不要波を検出し、２つのスペクトルピーク、前後の正常な流速の平均値や回帰曲線から求めた流速を使って補正することで、正常な流速を推定するようにしている。
この不要波の検出は、流れの回帰曲線からのズレが、標準偏差に対して経験的に決められた定数を掛けた値より大きいか否かにより判定している。しかしながら、多数の不要波がレンジ方向に分布する場合には、標準偏差に対して掛ける適正な定数を決定することが困難であり、流速の推定精度が劣化してしまうことがある。

以下の特許文献２には、多数の不要波がレンジ方向に分布する場合でも、正確な流速を推定することができるようにするために、同一ビームの複数のレンジビンの流速をメディアンフィルタに供給して海流の速度成分の代表値を抽出し、この代表値から所定の閾値を超える流速を不要波として除去したのち、流速を再検出することで流速を推定しているレーダ装置が開示されている。
実際に、メディアンフィルタを用いて不要波を除去するには、入り込む不要波の数の最大値Ｍを予測し、（２Ｍ＋１）のレンジビンを必要とする。

特開平１１−８３９９２号公報（請求項４） 特開２００２−３２３５５８号公報（段落番号［００１１］）

従来のレーダ装置は以上のように構成されているので、特許文献２の場合、多数の不要波がレンジ方向に分布する場合でも、正確な流速を推定することができるが、入り込む不要波の数を予測する必要があり、その予測数よりも多数の不要波が入り込む状況下では、誤って不要波を抽出してしまって、流速の推定精度が劣化してしまうことがある課題があった。

この発明は上記のような課題を解決するためになされたもので、入り込む不要波の数を予測することなく、観測対象の速度を高精度に推定することができるレーダ装置を得ることを目的とする。

この発明に係るレーダ装置は、観測対象に反射されたビームの反射波を周波数解析することで、レンジビン毎にドップラースペクトルを算出するスペクトル算出手段と、スペクトル算出手段により算出されたドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算し、加算処理後のドップラースペクトルのスペクトル値から観測対象の速度を推定する速度推定手段とを備え、その速度推定手段が、スペクトル算出手段により算出されたドップラースペクトルを構成する各々の周波数ビンを注目セルとして、当該注目セルのレンジビンより１レンジビン前のドップラースペクトルの中に、レンジ方向の１次散乱ピークのばらつき分に相当する幅を有する探索領域を設定し、その探索領域内の複数の周波数ビンの中で、スペクトル値が最も大きい周波数ビンを、加算対象の周波数ビンとして選択して、加算対象の周波数ビンにおけるスペクトル値の累積加算値であるスコア値を、注目セルとしている周波数ビンのスペクトル値に加算し、その加算結果を注目セルとしている周波数ビンのスコア値とするスペクトル値加算手段を備えるようにしたものである。

この発明によれば、スペクトル算出手段により算出されたドップラースペクトルを構成する各々の周波数ビンを注目セルとして、当該注目セルのレンジビンより１レンジビン前のドップラースペクトルの中に、レンジ方向の１次散乱ピークのばらつき分に相当する幅を有する探索領域を設定し、その探索領域内の複数の周波数ビンの中で、スペクトル値が最も大きい周波数ビンを、加算対象の周波数ビンとして選択して、加算対象の周波数ビンにおけるスペクトル値の累積加算値であるスコア値を、注目セルとしている周波数ビンのスペクトル値に加算し、その加算結果を注目セルとしている周波数ビンのスコア値とするスペクトル値加算手段を備えるように構成したので、入り込む不要波の数を予測することなく、観測対象の速度を高精度に推定することができる効果がある。

この発明の実施の形態１によるレーダ装置を示す構成図である。 この発明の実施の形態１によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。 レーダ装置によって観測されるドップラースペクトルを示す説明図である。 ＴＢＤ処理部１１により設定される探索領域と加算対象の周波数ビンの選択処理を示す説明図である。 バックトラック処理部１３による最大の周波数ビンの特定処理を示す説明図である。 平滑処理部１４による流速の平滑化処理を示す説明図である。 この発明の実施の形態２によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。 オフセット量算出部２２による相関ピークとなる周波数ビンの特定処理を示す説明図である。 ＴＢＤ処理部２３により設定される探索領域を示す説明図である。 この発明の実施の形態３によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。 探索領域決定部３１による探索領域の決定処理を示す説明図である。 この発明の実施の形態４によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。 この発明の実施の形態５によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。 異常値の検出処理を示す説明図である。 補正後のオフセット量を示す説明図である。 この発明の実施の形態６によるレーダ装置を示す構成図である。 この発明の実施の形態６によるレーダ装置の異常ピーク判定部５３を示す構成図である。 波形推定処理部５１による流速の予測値を示す説明図である。 流速偏差抽出部５２により算出される流速偏差（観測値と予測値の差分値）を示す説明図である。 異常ピークによる流速と津波による流速が混在している場合の推定波形を示す説明図である。 異常ピークによる流速と津波による流速が混在している場合の流速偏差を示す説明図である。 異常ピークによる流速が観測値に含まれている場合の異常ピーク候補の検出例を示す説明図である。 津波による流速が観測値に含まれている場合の異常ピーク候補の検出例を示す説明図である。 不要波の振幅値が１次散乱の加算値よりも大きいことによる不要波の検出例を示す説明図である。 １次散乱と不要波の振幅の確率密度分布を示す説明図である。 この発明の実施の形態７によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。 １次散乱の振幅特性を示す説明図である。 この発明の実施の形態７のＴＢＤ処理部７２による不要波の判定処理を示す説明図である。 この発明の実施の形態８によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。 この発明の実施の形態８のＴＢＤ処理部８１による不要波の判定処理を示す説明図である。 平均化による１次散乱の振幅比の確率密度分布を示す説明図である。

実施の形態１．
図１はこの発明の実施の形態１によるレーダ装置を示す構成図である。
この実施の形態１では、海面の波の速度である流速を推定する海洋レーダを例にして説明するが、海洋レーダに限るものではなく、例えば、上空の風のように空間的な広がりを持ち、高度方向に相関を持つ場合にも適用可能である。また、レンジ方向だけでなくビーム方向へ相関持つ場合にも適用が可能である。

図１において、ビーム送受信系１はレーダビームを海面（観測対象）に放射する一方、海面に反射されたレーダビームの反射波を受信する。
ビーム送受信系１は、例えば、フェーズドアレーアンテナと、レーダビームを生成して、そのレーダビームをフェーズドアレーアンテナに出力することで、そのレーダビームを海面に放射する送信機と、海面に反射されたレーダビームの反射波がフェーズドアレーアンテナに入射されると、その反射波を受信して、その受信信号に対する所定の信号受信処理（例えば、受信信号をデジタル信号に変換するＡ／Ｄ変換処理、ベースバンド信号に変換する周波数変換処理など）を実施した後の受信信号であるデジタル受信データを出力する受信機とから構成されている。

信号処理部２は例えばＣＰＵを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、ビーム送受信系１からデジタル受信データを受けると、そのデジタル受信データに対するＦＦＴ（Ｆａｓｔ Ｆｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ：高速フーリエ変換）処理等を実施することで、そのデジタル受信データの周波数を解析して、レンジビン（距離分解能）毎にドップラースペクトルを算出する処理を実施する。なお、信号処理部２はスペクトル算出手段を構成している。

２次元ピーク検出処理部３は例えばＣＰＵを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、信号処理部２により算出されたドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算し、加算処理後のドップラースペクトルのスペクトル値から海面の波のレーダ視線方向の速度である流速を推定する処理を実施する。なお、２次元ピーク検出処理部３は速度推定手段を構成している。
出力データ保存部４は例えばＲＡＭやハードディスクなどの記憶媒体から構成されており、２次元ピーク検出処理部３によりビーム毎に推定された流速を保存する。

図１の例では、レーダ装置の構成要素である信号処理部２、２次元ピーク検出処理部３及び出力データ保存部４のそれぞれが専用のハードウェアで構成されているものを想定しているが、信号処理部２、２次元ピーク検出処理部３及び出力データ保存部４がコンピュータで構成されていてもよい。
信号処理部２、２次元ピーク検出処理部３及び出力データ保存部４をコンピュータで構成する場合、出力データ保存部４をコンピュータのメモリ上に構成するとともに、信号処理部２及び２次元ピーク検出処理部３の処理内容を記述しているプログラムを当該コンピュータのメモリに格納し、当該コンピュータのＣＰＵが当該メモリに格納されているプログラムを実行するようにすればよい。

図２はこの発明の実施の形態１によるレーダ装置の２次元ピーク検出処理部３を示す構成図である。
図２において、ＴＢＤ（Ｔｒａｃｋ Ｂｅｆｏｒｅ Ｄｅｔｅｃｔ）処理部１１は信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セルとして、当該注目セルのレンジビンより１レンジビン前のドップラースペクトルの中に、レンジ方向の１次散乱ピークのばらつき分に相当する幅を有する探索領域を設定し、その探索領域内の複数の周波数ビンの中から、加算対象の周波数ビン（例えば、スペクトル値である振幅値が最大の周波数ビン）を選択して、加算対象の周波数ビンにおけるスペクトル値の累積加算値であるスコア値を注目セルのスペクトル値に加算し、その加算結果を注目セルのスコア値とする処理を実施する。なお、ＴＢＤ処理部１１はスペクトル値加算手段を構成している。以下、スコア値を「スコア」と表記する。
ＴＢＤ履歴結果保存部１２はＴＢＤ処理部１１により選択された加算対象の周波数ビンを示す周波数ビン番号を保存するとともに、加算処理後の注目セルのスペクトル値を保存する。

バックトラック処理部１３はＴＢＤ処理部１１による加算処理後のドップラースペクトルにおける各周波数ビンのスペクトル値を比較して、海面の波の流速に対応する周波数ビン（例えば、スペクトル値である振幅値が最大の周波数ビン）を特定し、その周波数ビンと予め設定された基準周波数から、海面の波の流速を算出する処理を実施する。
即ち、バックトラック処理部１３は加算処理後のドップラースペクトルのうち、加算処理の最終処理レンジビンにおけるドップラースペクトルの中で、例えば、スペクトル値である振幅値が最大の周波数ビンを特定したのち、ＴＢＤ履歴結果保存部１２に保存されている周波数ビン番号を参照して、最大の周波数ビンに対応する加算対象の周波数ビン（ＴＢＤ処理部１１により選択された加算対象の周波数ビン）を最初の処理ビンまで辿りながら、各レンジビンにおける最大の周波数ビンを特定し、レンジビン毎に、最大の周波数ビンと予め設定された基準周波数から、海面の波の流速を算出する処理を実施する。

平滑処理部１４はバックトラック処理部１３によりレンジビン毎に算出された流速を平滑化し、平滑化後の流速を出力データ保存部４に出力する処理を実施する。
なお、バックトラック処理部１３及び平滑処理部１４から速度算出手段が構成されている。

次に動作について説明する。
ビーム送受信系１は、例えば、海洋レーダ等で一般的に用いられるＦＭＣＷ（Ｆｒｅｑｕｅｎｃｙ Ｍｏｄｕｌａｔｉｏｎ Ｃｏｎｔｉｎｕｏｕｓ Ｗａｖｅ）方式や、ＦＭＩＣＷ（Ｆｒｅｑｕｅｎｃｙ Ｍｏｄｕｌａｔｉｏｎ Ｉｎｔｅｒｒｕｐｔｅｄ Ｃｏｎｔｉｎｕｏｕｓ Ｗａｖｅ）方式のレーダビームを海面に照射する。
ビーム送受信系１は、レーダビームを海面に照射したのち、海面に反射されたレーダビームの反射波を受信し、その受信信号に対する所定の信号受信処理を実施した後の受信信号であるデジタル受信データを信号処理部２に出力する。

信号処理部２は、ビーム送受信系１からデジタル受信データを受けると、そのデジタル受信データに対するＦＦＴ処理等を実施することで、そのデジタル受信データの周波数を解析して、レンジビン毎にドップラースペクトルを算出する。
ここで、図３はレーダ装置によって観測されるドップラースペクトルを示す説明図である。
図３に示すように、レーダ装置のアンテナから放射されるレーダビームの方位（もしくは角度）方向であるビーム方向に照射されたレーダビーム毎に、距離方向（以下、「レンジ方向」と称する）について、レンジビン単位のドップラースペクトルが得られる。このとき、ドップラースペクトルが周波数分解能毎に一定に区切られているセルを周波数ビンと称する。

ドップラースペクトルは、図３に示すように、レーダ装置に近づいてくる波の成分と、レーダ装置から遠ざかる波の成分とを有しており、ドップラースペクトルのピーク周波数として、プラス周波数領域のｆ_ｐと、マイナス周波数領域のｆ_ｍとがあるため、２つの流れが存在している。
図３において、ｆ_０は予めレーダ装置で設定される基準周波数であり、プラス周波数領域の基準周波数ではｆ_０、マイナス周波数領域の基準周波数では−ｆ_０となる。
２つの流れが同一の流れであれば、下記の式（１）が成立し、周波数ｆ_ｐ，ｆ_ｍから基準周波数ｆ_０，−ｆ_０をそれぞれ差し引けば、下記の式（２）のようになるため、ドップラーシフト分を計測することができ、流速の推定が可能になる。

２次元ピーク検出処理部３のＴＢＤ処理部１１は、信号処理部２からレンジビン毎のドップラースペクトルを受けると、ドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算する処理を実施する。
図４はＴＢＤ処理部１１により設定される探索領域と加算対象の周波数ビンの選択処理を示す説明図である。
以下、図４を参照しながら、ＴＢＤ処理部１１の処理内容を具体的に説明する。ここでは、ドップラースペクトルのスペクトル値として、ドップラースペクトルの振幅値を繰り返し加算する例を説明するが、ドップラースペクトルの振幅値に限るものではなく、例えば、ドップラースペクトルの尤度を繰り返し加算するものであってもよい。また、繰り返し加算する際に推移確率の概念を入れて加算するものであってもよい。

この実施の形態１では、レンジビンの数がＮ個あるものとし、ｉ番目のレンジビンをレンジビンｉで表記する。ｉ＝１，２，・・・，Ｎであり、１レンジ側をニアレンジ側、Ｎレンジ側をファーレンジ側と称する。
また、１つのレンジビンにおけるドップラースペクトルは、Ｍ個の周波数ビンからなり、ｊ番目の周波数ビンを周波数ビンｊで表記する。ｊ＝１，２，・・・，Ｍである。
ＴＢＤ処理（振幅値の繰り返し加算処理）を実施する方向は、ニアレンジ側からファーレンジ側への方向でもよいし、ファーレンジ側からニアレンジ側への方向でもよいが、この実施の形態１では、ニアレンジ側からファーレンジ側への方向にＴＢＤ処理を実施するものとする。

ＴＢＤ処理部１１は、ニアレンジ側のレンジビンからＴＢＤ処理を実施する際、各レンジビンにおける全ての周波数ビンを順番に注目セルとして振幅値の加算処理を実施するが、ここでは、レンジビンｉにおけるドップラースペクトルのｊ番目の注目セル（ｉ，ｊ）に着目して説明する。
このとき、注目セル（ｉ，ｊ）の振幅値をｘ_ｊ（ｉ）で表記し、加算処理後のスコア（振幅値の累積加算値）をＳ_ｊ（ｉ）で表記する。

ＴＢＤ処理部１１は、レンジビンｉにおけるドップラースペクトルのｊ番目の周波数ビンを注目セル（ｉ，ｊ）に設定すると、図４に示すように、レンジビンｉより１レンジビン前のレンジビン（ｉ−１）のドップラースペクトルの中に、レンジ方向の１次散乱ピークのばらつき分に相当する幅を有する探索領域を設定する。
図４の例では、５個の周波数ビンを含む探索領域（探索領域の範囲をΔの記号で表している）を設定しており、その探索領域の中心の周波数ビンは、注目セル（ｉ，ｊ）の周波数ビンと同じｊ番目の周波数ビンである。
なお、探索領域の範囲Δを広げすぎると、不要波を誤って抽出する可能性がある。一方、同時に観測される同一ビーム上では、１次散乱ピークはレンジ方向に連続的に分布し、２次散乱を含むブラッグ散乱スペクトルはレンジ方向に相関を持つことになる。そのため、探索領域の範囲Δは、レンジ方向の１次散乱ピークのばらつき分を想定して設定すればよく、必要以上に探索領域を広げる必要がない。

ＴＢＤ処理部１１は、レンジビン（ｉ−１）のドップラースペクトルの中に探索領域を設定すると、その探索領域内の複数の周波数ビン（図４の例では、（ｉ−１，ｊ−２）、（ｉ−１，ｊ−１）、・・・（ｉ−１，ｊ＋２）の周波数ビン）の中から、加算対象の周波数ビンを選択する。例えば、振幅値が最大の周波数ビンを選択する。
ＴＢＤ処理部１１は、加算対象の周波数ビンを選択すると、下記の式（３）に示すように、加算対象の周波数ビンにおける振幅値の累積加算値であるスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算することで、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出する。

式（３）において、λは忘却係数であり、０〜１の値が設定される。ｊ’は探索領域内の周波数ビンを示す変数であり、Ｓ_ｊ’（ｉ）は探索領域内の周波数ビンのスコア（振幅値の累積加算値）である。

ＴＢＤ処理部１１は、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出すると、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存するとともに、ＴＢＤ履歴結果として、注目セル（ｉ，ｊ）に対する加算対象の周波数ビンを示す周波数ビン番号ＳＰ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存する。

バックトラック処理部１３は、ＴＢＤ処理部１１のＴＢＤ処理が完了すると、最終処理レンジビンであるＮ番目のレンジビンのドップラースペクトルの中で、スコアが最大の周波数ビンを特定する。
ここで、図５はバックトラック処理部１３による最大の周波数ビンの特定処理を示す説明図である。
図５の例では、４番目の周波数ビンを最大の周波数ビンとして特定している。

バックトラック処理部１３は、最終処理レンジビンであるＮ番目のレンジビンのドップラースペクトルの中で、スコアが最大の周波数ビンを特定すると、図５に示すように、ＴＢＤ履歴結果保存部１２に保存されている周波数ビン番号ＳＰ_ｊ（ｉ）を参照して、ＴＢＤ処理の実施方向と逆方向にＴＢＤ履歴結果を辿ることで（図５の例では、ｉ＝Ｎ−１，Ｎ−２，・・・，１の順にＴＢＤ履歴結果を辿っている）、各レンジビンにおける最大の周波数ビンを特定する。
図５の例では、（Ｎ−１）番目のレンジビンでは、５番目の周波数ビンを最大の周波数ビンとして特定し、（Ｎ−２）番目のレンジビンでは、３番目の周波数ビンを最大の周波数ビンとして特定し、１番目のレンジビンでは、３番目の周波数ビンを最大の周波数ビンとして特定している。

バックトラック処理部１３は、各レンジビンにおける最大の周波数ビンを特定すると、レンジビン毎に、最大の周波数ビンと予め設定された基準周波数ｆ_０（あるいは、−ｆ_０）から、海面の波の流速を算出する。
平滑処理部１４は、バックトラック処理部１３がレンジビン毎に流速を算出すると、その流速のばらつきをレンジ方向で平滑化し、平滑化後の流速を出力データ保存部４に出力する。
平滑処理部１４による流速の平滑化処理には、例えば、ローパスフィルタ等の平滑化フィルタが用いられる。
図６は平滑処理部１４による流速の平滑化処理を示す説明図である。
これにより、２次元ピーク検出処理部３からばらつきが抑えられた流速が出力データ保存部４に出力される。

以上で明らかなように、この実施の形態１によれば、ＴＢＤ処理部１１が、信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セルとして、当該注目セルのレンジビンより１レンジビン前のドップラースペクトルの中に、レンジ方向の１次散乱ピークのばらつき分に相当する幅を有する探索領域を設定し、その探索領域内の複数の周波数ビンの中から、加算対象の周波数ビン（例えば、スペクトル値である振幅値が最大の周波数ビン）を選択して、加算対象の周波数ビンのスコアを注目セルの振幅値に加算し、その加算結果を注目セルのスコアとするように構成したので、入り込む不要波の数を予測することなく、海面の波の流速を高精度に推定することができる効果を奏する。

また、この実施の形態１によれば、バックトラック処理部１３が、加算処理後のドップラースペクトルのうち、加算処理の最終処理レンジビンにおけるドップラースペクトルの中で、例えば、振幅値が最大の周波数ビンを特定したのち、ＴＢＤ履歴結果保存部１２に保存されている周波数ビン番号を参照して、最大の周波数ビンに対応する加算対象の周波数ビンを最初の処理ビンまで辿りながら、各レンジビンにおける最大の周波数ビンを特定するように構成したので、不要波の数に依存せずに、高精度に１次散乱ピーク（各レンジビンにおける最大の周波数ビン）を抽出することができる効果を奏する。

実施の形態２．
上記実施の形態１では、探索領域の中心の周波数ビンが、注目セル（ｉ，ｊ）の周波数ビンと同じｊ番目の周波数ビンである例を示したが、探索領域を広げることなく、津波のような大きな流速変化に追従できるようにするために、注目セル（ｉ，ｊ）の探索領域にオフセットを与えるようにしてもよい。

図７はこの発明の実施の形態２によるレーダ装置の２次元ピーク検出処理部３を示す構成図であり、図７において、図２と同一符号は同一または相当部分を示すので説明を省略する。
相互相関処理部２１は信号処理部２により算出された各レンジビンのドップラースペクトルからレンジビン間の相互相関値を算出する処理を実施する。
オフセット量算出部２２は相互相関処理部２１により算出されたレンジビン間の相互相関値から相関ピークとなる周波数ビンを特定し、その周波数ビンを示す周波数ビン番号をオフセット量としてＴＢＤ処理部２３に出力する処理を実施する。

ＴＢＤ処理部２３は中心の周波数ビンが注目セル（ｉ，ｊ）の周波数ビンと同じである探索領域（図４に示す探索領域）をオフセット量算出部２２から出力されたオフセット量だけずらすことで（図９を参照）、探索領域に対応する周波数ビンの範囲を決定する処理を実施する。
また、ＴＢＤ処理部２３は信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セル（ｉ，ｊ）として、オフセット量だけずらしている探索領域内の複数の周波数ビンの中から、加算対象の周波数ビン（例えば、振幅値が最大の周波数ビン）を選択して、加算対象の周波数ビンのスコアを注目セル（ｉ，ｊ）の振幅値に加算し、その加算結果を注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）とする処理を実施する。なお、相互相関処理部２１、オフセット量算出部２２及びＴＢＤ処理部２３からスペクトル値加算手段が構成されている。

次に動作について説明する。
ただし、相互相関処理部２１、オフセット量算出部２２及びＴＢＤ処理部２３以外は、上記実施の形態１と同様であるため、ここでは、相互相関処理部２１、オフセット量算出部２２及びＴＢＤ処理部２３の処理内容だけを説明する。

２次元ピーク検出処理部３の相互相関処理部２１は、信号処理部２から各レンジビンのドップラースペクトルを受けると、下記の式（５）に示すように、各レンジビンのドップラースペクトルからレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）を算出する。

式（５）において、ｘ_ｊ（ｉ）はｉ番目のレンジビンのドップラースペクトルにおけるｊ番目の周波数ビンの振幅値である。
また、ＦＦＴはＦＦＴ処理を示す記号、ＩＦＦＴはＩｎｖｅｒｓｅＦＦＴ（ＦＦＴの逆変換）処理を示す記号、ｃｏｎｊは複素共役を示す記号である。
なお、相互相関処理は、畳み込み演算によって算出することができるので、式（５）以外の公知の算出方法で、レンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）を算出するようにしてもよい。

オフセット量算出部２２は、相互相関処理部２１がレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）を算出すると、その相互相関値ＣＯＲ_ｊ（ｉ）から相関ピークとなる周波数ビンを特定し、その周波数ビンを示す周波数ビン番号をオフセット量としてＴＢＤ処理部２３に出力する。
ここで、図８はオフセット量算出部２２による相関ピークとなる周波数ビンの特定処理を示す説明図である。
図８の例では、複数個の相関ピークが得られているが、相関ピークのサーチ範囲を所定の範囲で設定しておけば、サーチ範囲内の最大の相互相関値を検出することで、相関ピークを抽出することができる。

具体的には、ｉ番目のレンジビンについて、相互相関値ＣＯＲ_ｊ（ｉ）が下記の式（６）を満足するｊ番目の周波数ビン（相関ピークとなる周波数ビン）を選択し、その周波数ビンを示す周波数ビン番号をオフセット量ＯＳ（ｉ）としてＴＢＤ処理部２３に出力する。

式（６）において、ＯＳ（ｉ）はｉ番目のレンジビンの注目セル（ｉ，ｊ）に対する探索領域のオフセット量、Ｑ＝｛ｊ｜１＜＝ｊ＜＝Ｓ，Ｍ−Ｓ＋１＜＝ｊ＜＝Ｍ｝である。
また、Ｓは相関ピークのサーチ範囲、Ｍは周波数ビン数である。

ＴＢＤ処理部２３は、オフセット量算出部２２からオフセット量ＯＳ（ｉ）を受けると、図４に示す探索領域（中心の周波数ビンが注目セル（ｉ，ｊ）の周波数ビンと同じである探索領域）を、そのオフセット量ＯＳ（ｉ）だけずらすことで、探索領域に対応する周波数ビンの範囲を決定する。
ここで、図９はＴＢＤ処理部２３により設定される探索領域を示す説明図である。
図９の例では、オフセット量が周波数ビンの７個分であり、探索領域を図中右方向に、周波数ビンの７個分だけずらしている。
これにより、探索領域はオフセット量だけずれるが、探索領域の範囲Δは、上記実施の形態１と同様に、レンジ方向の１次散乱ピークのばらつき分に相当する幅であり、図９の例では、探索領域が５個の周波数ビンを含んでいる。

ＴＢＤ処理部２３は、レンジビン（ｉ−１）のドップラースペクトルの中に探索領域を設定すると、図２のＴＢＤ処理部１１と同様に、その探索領域内の複数の周波数ビンの中から、加算対象の周波数ビンを選択する。例えば、振幅値が最大の周波数ビンを選択する。
ＴＢＤ処理部２３は、加算対象の周波数ビンを選択すると、図２のＴＢＤ処理部１１と同様に、上記の式（３）を用いて、加算対象の周波数ビンのスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算することで、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出する。
ＴＢＤ処理部２３は、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出すると、図２のＴＢＤ処理部１１と同様に、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存するとともに、ＴＢＤ履歴結果として、注目セル（ｉ，ｊ）に対する加算対象の周波数ビンを示す周波数ビン番号ＳＰ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存する。
バックトラック処理部１３及び平滑処理部１４の処理内容は、上記実施の形態１と同様である。

以上で明らかなように、この実施の形態２によれば、信号処理部２により算出された各レンジビンのドップラースペクトルからレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）を算出する相互相関処理部２１と、相互相関処理部２１により算出されたレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）から相関ピークとなる周波数ビンを特定し、その周波数ビンを示す周波数ビン番号をオフセット量ＯＳ（ｉ）として出力するオフセット量算出部２２とを設け、ＴＢＤ処理部２３が、オフセット量算出部２２から出力されたオフセット量ＯＳ（ｉ）だけ探索領域をずらして、その探索領域に対応する周波数ビンの範囲を決定するように構成したので、探索領域を広げることなく、津波のような大きな流速変化に追従できるようになり、その結果、不要波の数に依存せずに、流速変化が大きい津波等の流速を高精度に推定することができる効果を奏する。

実施の形態３．
上記実施の形態１，２では、予め、レンジ方向の１次散乱ピークのばらつき分に相当する幅を有する探索領域を設定するものを示したが、その探索領域の周波数位置及び範囲Δを自動的に算出するようにしてもよい。

図１０はこの発明の実施の形態３によるレーダ装置の２次元ピーク検出処理部３を示す構成図であり、図１０において、図７と同一符号は同一または相当部分を示すので説明を省略する。
探索領域決定部３１は相互相関処理部２１により算出されたレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）から、その相互相関値ＣＯＲ_ｊ（ｉ）の相関値幅（最大値ＣＯＲ_ｍａｘと最小値ＣＯＲ_ｍｉｎの差分）を算出して、予め設定された係数ｒを当該相関値幅に乗算し、相互相関値ＣＯＲ_ｊ（ｉ）が当該乗算結果より大きくなる周波数ビンの範囲を特定し、その周波数ビンの範囲を探索領域に決定する処理を実施する。

ＴＢＤ処理部３２は信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セルとして、当該注目セルのレンジビンより１レンジビン前のドップラースペクトルの中に、探索領域決定部３１により決定された探索領域を設定し、その探索領域内の複数の周波数ビンの中から、加算対象の周波数ビン（例えば、振幅値が最大の周波数ビン）を選択して、加算対象の周波数ビンのスコアを注目セル（ｉ，ｊ）の振幅値に加算し、その加算結果を注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）とする処理を実施する。
なお、相互相関処理部２１、探索領域決定部３１及びＴＢＤ処理部３２からスペクトル値加算手段が構成されている。

次に動作について説明する。
ただし、探索領域決定部３１及びＴＢＤ処理部３２以外は、上記実施の形態２と同様であるため、ここでは、探索領域決定部３１及びＴＢＤ処理部３２の処理内容だけを説明する。
図１１は探索領域決定部３１による探索領域の決定処理を示す説明図である。

探索領域決定部３１は、相互相関処理部２１がレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）を算出すると、ｉ番目のレンジビンにおける注目セル（ｉ，ｊ）に対する探索領域を決定する際、その相互相関値ＣＯＲ_ｊ（ｉ）の相関値幅（最大値ＣＯＲ_ｍａｘと最小値ＣＯＲ_ｍｉｎの差分）を算出し、予め設定された係数ｒを当該相関値幅に乗算する。
図１１では、この乗算結果をｒ（ＣＯＲ_ｍａｘ−ＣＯＲ_ｍｉｎ）で表記している。

そして、探索領域決定部３１は、相関ピークを中心にして、相互相関値ＣＯＲ_ｊ（ｉ）が乗算結果ｒ（ＣＯＲ_ｍａｘ−ＣＯＲ_ｍｉｎ）より大きくなる周波数ビンの範囲を特定し、その周波数ビンの範囲を探索領域に決定する。
具体的には、相互相関値ＣＯＲ_ｊ（ｉ）が下記の式（７）を満足する周波数ビンの範囲を探索領域に決定する。

式（７）において、集合Ｐ＝｛ｊ｜１＜＝ｊ＜＝Ｍ｝、集合Ｑ＝｛ｊ｜１＜＝ｊ＜＝Ｓ，Ｍ−Ｓ＋１＜＝ｊ＜＝Ｍ｝である。Ｓは相関ピークのサーチ範囲、Ｍは周波数ビン数である。

ＴＢＤ処理部３２は、探索領域決定部３１が探索領域を決定すると、注目セル（ｉ，ｊ）のレンジビンより１レンジビン前のドップラースペクトルの中に、その探索領域を設定する。
ＴＢＤ処理部３２は、レンジビン（ｉ−１）のドップラースペクトルの中に探索領域を設定すると、図２のＴＢＤ処理部１１と同様に、その探索領域内の複数の周波数ビンの中から、加算対象の周波数ビンを選択する。例えば、振幅値が最大の周波数ビンを選択する。
ＴＢＤ処理部３２は、加算対象の周波数ビンを選択すると、図２のＴＢＤ処理部１１と同様に、上記の式（３）を用いて、加算対象の周波数ビンのスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算することで、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出する。
ＴＢＤ処理部３２は、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出すると、図２のＴＢＤ処理部１１と同様に、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存するとともに、ＴＢＤ履歴結果として、注目セル（ｉ，ｊ）に対する加算対象の周波数ビンを示す周波数ビン番号ＳＰ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存する。
バックトラック処理部１３及び平滑処理部１４の処理内容は、上記実施の形態１と同様である。

以上で明らかなように、この実施の形態３によれば、探索領域決定部３１が、相互相関処理部２１により算出されたレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）から、その相互相関値ＣＯＲ_ｊ（ｉ）の相関値幅（最大値ＣＯＲ_ｍａｘと最小値ＣＯＲ_ｍｉｎの差分）を算出して、予め設定された係数ｒを当該相関値幅に乗算し、相互相関値ＣＯＲ_ｊ（ｉ）が当該乗算結果ｒ（ＣＯＲ_ｍａｘ−ＣＯＲ_ｍｉｎ）より大きくなる周波数ビンの範囲を特定し、その周波数ビンの範囲を探索領域に決定するように構成したので、上記実施の形態１と同様の効果を奏する他に、事前に探索領域を設定することなく、適正な探索領域を自動的に設定することができる効果を奏する。

実施の形態４．
上記実施の形態３では、探索領域決定部３１を実装している２次元ピーク検出処理部３を示したが、２次元ピーク検出処理部３がオフセット量算出部２２及び探索領域決定部３１を実装しているものであってもよい。
図１２はこの発明の実施の形態４によるレーダ装置の２次元ピーク検出処理部３を示す構成図であり、図１２において、図７及び図１０と同一符号は同一または相当部分を示すので説明を省略する。
ＴＢＤ処理部３３は探索領域決定部３１により決定された探索領域（図１１に示す探索領域）をオフセット量算出部２２から出力されたオフセット量だけずらすことで、探索領域に対応する周波数ビンの範囲を決定する処理を実施する。
また、ＴＢＤ処理部３３は信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セル（ｉ，ｊ）として、オフセット量だけずらしている探索領域内の複数の周波数ビンの中から、加算対象の周波数ビン（例えば、振幅値が最大の周波数ビン）を選択して、加算対象の周波数ビンのスコアを注目セル（ｉ，ｊ）の振幅値に加算し、その加算結果を注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）とする処理を実施する。なお、相互相関処理部２１、オフセット量算出部２２、探索領域決定部３１及びＴＢＤ処理部３３からスペクトル値加算手段が構成されている。

この実施の形態４では、探索領域決定部３１を実装しているので、上記実施の形態３と同様に、事前に探索領域を設定することなく、適正な探索領域を自動的に設定することができる。
また、この実施の形態４では、図１０の２次元ピーク検出処理部３に対して、オフセット量算出部２２を付加しているので、上記実施の形態２と同様に、探索領域を広げることなく、津波のような大きな流速変化に追従することができる。

実施の形態５．
図１３はこの発明の実施の形態５によるレーダ装置の２次元ピーク検出処理部３を示す構成図であり、図１３において、図１２と同一符号は同一または相当部分を示すので説明を省略する。
異常時処理部４１は相互相関処理部２１により算出されたレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）における相関ボトム（最小値ＣＯＲ_ｍｉｎ）に対する相関ピーク（最大値ＣＯＲ_ｍａｘ）の比（相関ピーク（最大値ＣＯＲ_ｍａｘ）と相関ボトム（最小値ＣＯＲ_ｍｉｎ）の比）を算出し、その比が予め設定された閾値ＣＯＲ_ｔｈ（第１の閾値）より小さい場合、異常時の処理として、オフセット量算出部２２により算出されたオフセット量を補正する処理を実施する。
また、異常時処理部４１は相関対象である（ｉ＋１）番目のレンジビンのドップラースペクトルをオフセット量算出部２２により算出されたオフセット量だけずらし、そのオフセット量だけずらした（ｉ＋１）番目のレンジビンのドップラースペクトルと、相関元であるｉ番目のレンジビンのドップラースペクトルとの相関係数ＣＣ（ｉ）を算出し、その相関係数ＣＣ（ｉ）が予め設定された閾値ＣＣ_ｔｈ（第２の閾値）より小さい場合、異常時の処理として、オフセット量算出部２２により算出されたオフセット量を補正する処理を実施する。

さらに、異常時処理部４１は相関ピーク（最大値ＣＯＲ_ｍａｘ）と相関ボトム（最小値ＣＯＲ_ｍｉｎ）が閾値ＣＯＲ_ｔｈより大きく、かつ、相関係数ＣＣ（ｉ）が閾値ＣＣ_ｔｈより大きい場合、各レンジビンのオフセット量の絶対値と予め設定された閾値ＯＳ_ｔｈ１（第３の閾値）を比較し、当該レンジビンのオフセット量の絶対値が閾値ＯＳ_ｔｈ１より大きく、かつ、当該レンジビンと隣接しているレンジビンのオフセット量の絶対値が予め設定された閾値ＯＳ_ｔｈ２（第４の閾値）より小さい場合、当該レンジビンのオフセット量を隣接しているレンジビンのオフセット量を用いて補正する処理を実施する。

ＴＢＤ処理部４２は探索領域決定部３１により決定された探索領域（図１１に示す探索領域）をオフセット量算出部２２から出力されたオフセット量だけずらすことで、探索領域に対応する周波数ビンの範囲を決定する処理を実施する。ただし、異常時処理部４１によりオフセット量が補正された場合、補正後のオフセット量だけずらすようにする。
また、ＴＢＤ処理部４２は信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セル（ｉ，ｊ）として、オフセット量だけずらしている探索領域内の複数の周波数ビンの中から、加算対象の周波数ビン（例えば、振幅値が最大の周波数ビン）を選択して、加算対象の周波数ビンのスコアを注目セル（ｉ，ｊ）の振幅値に加算し、その加算結果を注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）とする処理を実施する。なお、相互相関処理部２１、オフセット量算出部２２、探索領域決定部３１、異常時処理部４１及びＴＢＤ処理部４２からスペクトル値加算手段が構成されている。

次に動作について説明する。
ただし、異常時処理部４１及びＴＢＤ処理部４２以外は、上記実施の形態４と同様であるため、ここでは、異常時処理部４１及びＴＢＤ処理部４２の処理内容だけを説明する。
レンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）における相関ピーク（最大値ＣＯＲ_ｍａｘ）から、１次散乱のピークのオフセット量を算出することができるが、ドップラースペクトルが多峰性となる場合や、１次散乱が他の散乱よりも弱くなる場合には、相関ピークが不正常になって、オフセット量が異常値になることがある。
そのため、上記実施の形態２，４では、オフセット量が異常値になることがあり、その結果として、流速の推定結果に誤差が含まれることがある。
この実施の形態５では、オフセット量が異常値となる場合には、そのオフセット量を補正して、正しい流速の推定結果を得ることができるようにしている。
具体的には、以下の通りである。

異常時処理部４１は、相互相関処理部２１がレンジビン間の相互相関値ＣＯＲ_ｊ（ｉ）を算出すると、その相互相関値ＣＯＲ_ｊ（ｉ）における相関ボトム（最小値ＣＯＲ_ｍｉｎ）に対する相関ピーク（最大値ＣＯＲ_ｍａｘ）の比（ＣＯＲ_ｍａｘ／ＣＯＲ_ｍｉｎ）を算出する。
異常時処理部４１は、相関ボトムに対する相関ピークの比（ＣＯＲ_ｍａｘ／ＣＯＲ_ｍｉｎ）を算出すると、下記の式（８）に示すように、相関ボトムに対する相関ピークの比（ＣＯＲ_ｍａｘ／ＣＯＲ_ｍｉｎ）と予め設定された閾値ＣＯＲ_ｔｈを比較し、相関ボトムに対する相関ピークの比（ＣＯＲ_ｍａｘ／ＣＯＲ_ｍｉｎ）が閾値ＣＯＲ_ｔｈより小さい場合、オフセット量算出部２２により算出されたオフセット量は異常値であると判定する。

また、異常時処理部４１は、相関ボトムに対する相関ピークの比（ＣＯＲ_ｍａｘ／ＣＯＲ_ｍｉｎ）が閾値ＣＯＲ_ｔｈより大きい場合、相関対象である（ｉ＋１）番目のレンジビンのドップラースペクトルをオフセット量算出部２２から出力されたオフセット量だけずらし、下記の式（９）に示すように、そのオフセット量だけずらした（ｉ＋１）番目のレンジビンのドップラースペクトルと、相関元であるｉ番目のレンジビンのドップラースペクトルとの相関係数ＣＣ（ｉ）を算出する。

式（９）において、ｘ_ｊ（ｉ）は相関元の注目セル（ｉ，ｊ）の振幅値、ｘ_{ｊ−ｏｆｆｓｅｔ}（ｉ＋１）はオフセット補正後の相関対象のセル（ｉ＋１，ｊ−ｏｆｆｓｅｔ）の振幅値、ｏｆｆｓｅｔはオフセット量算出部２２により算出されたオフセット量である。
また、ｘ（ｉ）バー（明細書の文章中では、電子出願の関係上、文字の上に“−”の記号を付することができないので、ｘ（ｉ）バーのように表記している）はｉ番目のレンジビンにおけるドップラースペクトルの振幅平均値、ｘ（ｉ＋１）バーは（ｉ＋１）番目のレンジビンにおけるドップラースペクトルの振幅平均値である。

異常時処理部４１は、相関係数ＣＣ（ｉ）を算出すると、下記の式（１０）に示すように、その相関係数ＣＣ（ｉ）と予め設定された閾値ＣＣ_ｔｈを比較し、その相関係数ＣＣ（ｉ）が閾値ＣＣ_ｔｈより小さい場合、オフセット量算出部２２により算出されたオフセット量は異常値であると判定する。

異常時処理部４１は、相関ボトムに対する相関ピークの比（ＣＯＲ_ｍａｘ／ＣＯＲ_ｍｉｎ）が閾値ＣＯＲ_ｔｈより小さいために、オフセット量算出部２２により算出されたオフセット量が異常値であると判定した場合、あるいは、相関係数ＣＣ（ｉ）が閾値ＣＣ_ｔｈより小さいために、オフセット量算出部２２により算出されたオフセット量が異常値であると判定した場合、そのオフセット量を零に補正（あるいは、事前に設定された値に補正）し、補正後のオフセット量をＴＢＤ処理部４２に出力する。

異常時処理部４１は、レンジビン毎に、上記の異常判定を行っても、オフセット量算出部２２により算出されたオフセット量が異常値であることが認められない場合（相関ボトムに対する相関ピークの比（ＣＯＲ_ｍａｘ／ＣＯＲ_ｍｉｎ）が閾値ＣＯＲ_ｔｈ大きく、かつ、相関係数ＣＣ（ｉ）が閾値ＣＣ_ｔｈより大きい場合）、全相関レンジビン処理後のオフセット量から異常値の判定を行う。
即ち、異常時処理部４１は、下記の式（１１）が成立するか否かを判定し、式（１１）が成立する場合、ｉ番目のオフセット量ＯＳ（ｉ）を異常値候補として検出する。

式（１１）において、ＯＳ_ｔｈ１は異常値候補検出用の閾値（第３の閾値）である。

異常時処理部４１は、異常値候補を検出すると、その異常値候補のレンジビンと隣接しているレンジビンのオフセット量が正常値であるか否かを判定し、隣接しているレンジビンのオフセット量が正常値であれば、その異常値候補は異常値であると判定する。
隣接しているレンジビンのオフセット量は、下記の式（１２）が成立すれば、正常値であると判定する。

式（１２）において、ＯＳ_ｔｈ２は正常値判定用の閾値（第４の閾値）である。

異常時処理部４１は、異常値候補が異常値であると判定すると、その異常値候補と隣接しているレンジビンのオフセット量を用いて補正する。
例えば、異常値候補のレンジビンと隣接しているレンジビンのオフセット量の平均値を算出し、その異常値候補のオフセット量を当該平均値に置き換える補正を行う。
なお、異常時処理部４１は、各レンジビンのオフセット量が正常値であれば、そのオフセット量をＴＢＤ処理部４２に出力し、各レンジビンのオフセット量が異常値であって、そのオフセット量を補正していれば、補正後のオフセット量をＴＢＤ処理部４２に出力する。

ここで、図１４は異常値の検出処理を示す説明図であり、図１５は補正後のオフセット量を示す説明図である。
異常値の検出は、図１４に示すように、異常値候補検出用の閾値ＯＳ_ｔｈ１を超えるオフセット量を異常値候補として検出し、隣接している複数個のレンジビンのオフセット量が正常値であれば、異常値候補を異常値と判定している。
異常値であるオフセット量の補正は、図１５に示すように、異常値と判定されたオフセット量と隣接している複数個のレンジビンのオフセット量が正常値であれば、隣接しているレンジビンの複数個のオフセット量を用いて行う。
図１５の例では、異常値の前後２個のオフセット量の平均値を異常値のオフセット量に置き換える補正を行っている。
ここでは、異常値の前後２個のオフセット量を用いているが、これは一例に過ぎず、例えば、異常値の前後１個のオフセット量を用いてもよい。また、異常値の前後２個のオフセット量の平均値ではなく、前後の正常値に置き換える補正を行ってもよい。

ＴＢＤ処理部４２は、異常時処理部４１から各レンジビンのオフセット量を受けると、探索領域決定部３１により決定された探索領域（図１１に示す探索領域）を、そのオフセット量だけずらすことで、探索領域に対応する周波数ビンの範囲を決定する。ただし、例えば、そのオフセット量が零に補正されていれば、探索領域はずらさない。
ＴＢＤ処理部４２は、レンジビン（ｉ−１）のドップラースペクトルの中に探索領域を設定すると、図２のＴＢＤ処理部１１と同様に、その探索領域内の複数の周波数ビンの中から、加算対象の周波数ビンを選択する。例えば、振幅値が最大の周波数ビンを選択する。
ＴＢＤ処理部４２は、加算対象の周波数ビンを選択すると、図２のＴＢＤ処理部１１と同様に、上記の式（３）を用いて、加算対象の周波数ビンのスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算することで、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出する。
ＴＢＤ処理部４２は、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）を算出すると、図２のＴＢＤ処理部１１と同様に、注目セル（ｉ，ｊ）のスコアＳ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存するとともに、ＴＢＤ履歴結果として、注目セル（ｉ，ｊ）に対する加算対象の周波数ビンを示す周波数ビン番号ＳＰ_ｊ（ｉ）をＴＢＤ履歴結果保存部１２に保存する。
バックトラック処理部１３及び平滑処理部１４の処理内容は、上記実施の形態１と同様である。

以上で明らかなように、この実施の形態５によれば、オフセット量算出部２２により算出されたオフセット量が異常値であると判定すれば、そのオフセット量を補正するように構成したので、ドップラースペクトルが多峰性となる場合や、１次散乱が他の散乱よりも弱くなる場合でも、正しい流速の推定結果を得ることができる効果を奏する。

実施の形態６．
図１６はこの発明の実施の形態６によるレーダ装置を示す構成図であり、図１６において、図１と同一符号は同一または相当部分を示すので説明を省略する。
波形推定処理部５１は例えばＣＰＵを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、２次元ピーク検出処理部３により推定された流速である観測値から、その流速のトレンドを把握することで、その流速の予測値を算出する処理を実施する。
流速偏差抽出部５２は例えばＣＰＵを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、２次元ピーク検出処理部３により推定された流速である観測値と波形推定処理部５１により算出された予測値との差分値を算出する処理を実施する。
なお、波形推定処理部５１及び流速偏差抽出部５２から速度偏差算出手段が構成されている。

異常ピーク判定部５３は例えばＣＰＵを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、各レーダビームのレンジビン毎に、流速偏差抽出部５２により算出された差分値と予め設定された第５の閾値を比較して、その差分値が第５の閾値より大きいレンジビンを異常ピーク候補として検出し、その異常ピーク候補を有するレーダビームと隣接しているレーダビームが有する異常ピーク候補の数が予め設定された第６の閾値より少なければ、第５の閾値より大きい差分値の算出元の観測値は異常値であると判定する処理を実施する。なお、異常ピーク判定部５３は正否判定手段を構成している。

図１６の例では、レーダ装置の構成要素である信号処理部２、２次元ピーク検出処理部３、出力データ保存部４、波形推定処理部５１、流速偏差抽出部５２及び異常ピーク判定部５３のそれぞれが専用のハードウェアで構成されているものを想定しているが、信号処理部２、２次元ピーク検出処理部３、出力データ保存部４、波形推定処理部５１、流速偏差抽出部５２及び異常ピーク判定部５３がコンピュータで構成されていてもよい。
信号処理部２、２次元ピーク検出処理部３、出力データ保存部４、波形推定処理部５１、流速偏差抽出部５２及び異常ピーク判定部５３をコンピュータで構成する場合、出力データ保存部４をコンピュータのメモリ上に構成するとともに、信号処理部２、２次元ピーク検出処理部３、波形推定処理部５１、流速偏差抽出部５２及び異常ピーク判定部５３の処理内容を記述しているプログラムを当該コンピュータのメモリに格納し、当該コンピュータのＣＰＵが当該メモリに格納されているプログラムを実行するようにすればよい。

図１７はこの発明の実施の形態６によるレーダ装置の異常ピーク判定部５３を示す構成図である。
図１７において、異常ピーク候補検出部６１は各レーダビームのレンジビン毎に、流速偏差抽出部５２により算出された差分値と予め設定された第５の閾値を比較して、その差分値が第５の閾値より大きいレンジビンを異常ピーク候補として検出する処理を実施する。
異常ピーク判定処理部６２は異常ピーク候補検出部６１により検出された異常ピーク候補を有するレーダビームと隣接しているレーダビームが有する異常ピーク候補の数が予め設定された第６の閾値より少なければ、第５の閾値より大きい差分値の算出元の観測値は異常値であると判定し、その判定結果を出力データ保存部４に格納する処理を実施する。

次に動作について説明する。
ただし、波形推定処理部５１、流速偏差抽出部５２及び異常ピーク判定部５３以外は、上記実施の形態１〜５と同様であるため、ここでは、波形推定処理部５１、流速偏差抽出部５２及び異常ピーク判定部５３の処理内容だけを説明する。

ドップラースペクトルが多峰性となる場合や、１次散乱が他の散乱よりも弱くなる場合、特に２つ山のスペクトルがある距離の範囲で連続して存在する（以下、「異常ピーク」と称する）場合には、１次散乱ピークを精度良く抽出することができず、不正常な流速の推定となる。
このため、上記実施の形態２，４，５のように、津波のような大きな流速を検出することが可能なレーダ装置では、異常ピークによる流速と津波による流速を識別することができない。
即ち、上記実施の形態２，４，５では、同一ビームのドップラースペクトルから各レンジの流速を推定しているので、異常ピークによる流速であるのか、津波による流速であるのかを識別することができないが、この実施の形態６では、複数のビームのドップラースペクトルから異常ピークを判定することで、異常ピークによる流速であるのか、津波による流速であるのかを識別することができるようにしている。

波形推定処理部５１は、２次元ピーク検出処理部３により推定された流速である観測値
から、その流速のトレンドを把握することで、その流速の予測値を算出する。
例えば、ローパスフィルタやカルマンフィルタ等を用いて、流速の観測値をフィルタリングすることで、その流速の予測値を算出する。
ローパスフィルタを用いる場合には、ローパスフィルタの通過帯域を設定する必要がある。カルマンフィルタを用いる場合、観測値から流速と流速変化率を推定することができるため、次時刻の流速の予測値の算出が可能になる。
図１８は波形推定処理部５１による流速の予測値を示す説明図である。
図１８において、ビームとレンジ毎に区切られているセルはレンジセルと呼ばれ、レンジセル毎に予測値が算出される。

流速偏差抽出部５２は、波形推定処理部５１が流速の予測値を算出すると、レンジセル毎に、その予測値と２次元ピーク検出処理部３から出力された流速の観測値との差分値を算出し、その差分値である流速偏差を異常ピーク判定部５３に出力する。
ここで、図１９は流速偏差抽出部５２により算出される流速偏差（観測値と予測値の差分値）を示す説明図である。
異常ピークによる流速と津波による流速は、平時（通常時）の流速よりも大きな流速となる。
図２０は異常ピークによる流速と津波による流速が混在している場合の推定波形を示す説明図である。
また、図２１は異常ピークによる流速と津波による流速が混在している場合の流速偏差を示す説明図である。
異常ピークによる流速と津波による流速が混在している観測値から予測値が算出され、その予測値と観測値から算出された図２１に示すような流速偏差だけでは、異常ピークによる成分であるのか、津波による成分であるのかは不明である。

異常ピーク判定部５３の異常ピーク候補検出部６１は、流速偏差抽出部５２がレンジセル毎に流速偏差を算出すると、レンジセル毎に、当該レンジセルの流速偏差と予め設定された第５の閾値を比較し、その流速偏差が第５の閾値より大きければ、当該レンジセルを異常ピーク候補として検出する。
ここで、図２２は異常ピークによる流速が観測値に含まれている場合の異常ピーク候補の検出例を示す説明図であり、図２３は津波による流速が観測値に含まれている場合の異常ピーク候補の検出例を示す説明図である。
図２２及び図２３において、レンジ方向とビーム方向で区切られているセルがレンジセルであり、黒色のセルは異常ピーク候補が検出されたセルであることを示している。

異常ピークが発生している状況下では、図２２に示すように、注目ビームに隣接しているビームには、異常ピーク候補の数が少ない。
一方、津波が発生している状況下では、図２３に示すように、流速偏差が閾値を超えるレンジセルが空間的に広がるため、注目ビームに隣接しているビームには、多くの数の異常ピーク候補が検出される。
したがって、津波が発生している状況下では、ビーム方向の相関が大きく、異常ピークが発生している状況下では、ビーム方向の相関が小さくなる。

異常ピーク判定部５３の異常ピーク判定処理部６２は、異常ピーク候補検出部６１が異常ピーク候補を検出すると、異常ピーク候補を有する注目のレーダビームと隣接しているレーダビームが有する異常ピーク候補の数を計数し、隣接しているレーダビームが有する異常ピーク候補の数が予め設定された第６の閾値より少なければ、注目のレーダビームが有する異常ピーク候補に係る流速偏差の算出元の観測値は異常値であると判定する。
例えば、異常ピークによる流速が観測値に含まれているため、図２２のように異常ピーク候補が検出されている場合、仮に、ｋ番目のレーダビームを注目ビームとすると、注目ビームｋは、３個の異常ピーク候補を有している。
しかし、注目ビームｋの隣接ビームｋ−１，ｋ＋１では、注目ビームｋが有する３個の異常ピーク候補と同一のレンジビンにおいて、異常ピーク候補が１個も検出されていない。
したがって、隣接ビームｋ−１，ｋ＋１が有している異常ピーク候補の数は０個であり、例えば、第６の閾値が“１”であるとすれば、注目ビームｋが有する異常ピーク候補に係る流速偏差の算出元の観測値は異常値であると判定される。

また、津波による流速が観測値に含まれているため、図２３のように異常ピーク候補が検出されている場合、仮に、ｋ番目のレーダビームを注目ビームとすると、注目ビームｋは、３個の異常ピーク候補を有している。
このとき、注目ビームｋの隣接ビームｋ−１，ｋ＋１は、注目ビームｋが有する３個の異常ピーク候補と同一のレンジビンにおいて、６個の異常ピーク候補が検出されている。
したがって、隣接ビームｋ−１，ｋ＋１が有している異常ピーク候補の数は６個であり、例えば、第６の閾値が“１”であるとすれば、注目ビームｋが有する異常ピーク候補に係る流速偏差の算出元の観測値は異常値でないと判定される。

この実施の形態６では、隣接ビームｋ−１，ｋ＋１が有している異常ピーク候補の数を第６の閾値と比較している例を示しているが、隣接ビームの対象を増やし、例えば、隣接ビームｋ−２，ｋ−１，ｋ＋１，ｋ＋２が有している異常ピーク候補の数を第６の閾値と比較するようにしてもよい。
また、隣接ビームｋ−１，ｋ＋１を除外し、例えば、隣接ビームｋ−２，ｋ＋２が有している異常ピーク候補の数を第６の閾値と比較するようにしてもよい。

以上で明らかなように、この実施の形態６によれば、異常ピーク判定部５３が、流速偏差抽出部５２により算出されたレンジセル毎の流速偏差と第５の閾値を比較して、その流速偏差が第５の閾値より大きいレンジセルを異常ピーク候補として検出し、その異常ピーク候補を有するレーダビームと隣接しているレーダビームが有する異常ピーク候補の数が第６の閾値より少なければ、第５の閾値より大きい流速偏差の算出元の観測値は異常値であると判定するように構成したので、異常ピークによる流速であるのか、津波による流速であるのかを識別することができる効果を奏する。

実施の形態７．
海域によっては、船等による不要波の振幅値が１次散乱の振幅値よりも非常に大きくなることがある。このような場合、上記実施の形態１のように、振幅値が最大の周波数ビンを加算対象の周波数ビンとして、注目セルのスペクトル値に加算するようにすると、不要波の周波数ビンを誤って抽出してしまって、海面の波の流速推定に誤差が生じることがある。

図２４は不要波の振幅値が１次散乱の加算値よりも大きいことによる不要波の検出例を示す説明図である。
図２４に示すように、不要波の周波数ビンが２レンジビンに連続して存在する場合、１次散乱を加算した最終レンジビンのセルのスコアよりも、不要波の周波数ビンを加算した最終レンジビンのセルのスコアの方が大きいと、バックトラックで抽出される周波数ビンとして、不要波のセルを加算した経路を抽出することになる。
バックトラックで抽出される周波数ビンとして、不要波のセルを加算した経路を抽出する状況は、特に１次散乱を複数回加算したスコアよりも不要波のセルの振幅値の方が大きい場合に生じる。

上記のような状況を回避するには、注目セルの振幅値が不要波である否かを判定する必要がある。
図２５は１次散乱と不要波の振幅の確率密度分布を示す説明図であるが、上記のような状況を回避する方法として、図２５に示すような確率密度分布において、例えば、分布の９９．７％となる値を閾値に設定し、その閾値を超える振幅値を有する信号を不要波として判定する方法が考えられる。
ただし、この方法を用いる場合、海域によって１次散乱の振幅が変わるため、海域毎に１次散乱の振幅の確率密度分布を抽出する必要がある。

上記実施の形態１〜６のレーダ装置では、注目セルの振幅値が不要波であるか否かを判定する処理を実施しないため、非常に大きな振幅を有する不要波を加算することがあり、その結果、バックトラックで不要波のセルを加算した経路を誤って抽出してしまって、不正常な流速推定となることがある。
そこで、この実施の形態７では、１レンジビン前の探索範囲内で最大のスコアとなるセルの振幅値に対する注目セルの振幅値の比を算出し、その比が予め設定された不要波判定用の閾値より大きければ、注目セルの振幅値が不要波であると判定するようにする。また、注目セルの振幅値が不要波であると判定すると、１レンジビン前の探索範囲内で最大のスコアを注目セルの振幅値に加算せずに、注目セルの加算処理後のスコアとして、注目セルの１レンジビン前のスコアを用いるようにすることで、正常な流速推定を行えるようにする。

図２６はこの発明の実施の形態７によるレーダ装置の２次元ピーク検出処理部３を示す構成図であり、図２６において、図２と同一符号は同一または相当部分を示すので説明を省略する。
不要波判定用閾値設定部７１は隣接レンジビンの１次散乱の振幅値の比の確率密度分布から、不要波判定用の閾値（第７の閾値）を設定する処理を実施する。不要波判定用閾値設定部７１は閾値設定手段を構成している。
なお、１次散乱はレンジ方向に連続的に変化しており、ある海域の１次散乱の振幅値の比の確率密度分布は他の海域でも同等であるため、海域毎に確率密度分布を抽出する必要がない。

ＴＢＤ処理部７２は信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セルとして、加算対象の周波数ビンのスペクトル値である振幅値（１レンジビン前の探索範囲内で最大のスコアとなるセルの振幅値）に対する当該注目セルの振幅値の比（注目セルの振幅値と加算対象の周波数ビンの振幅値との比）を算出し、その比が不要波判定用閾値設定部７１により設定された不要波判定用の閾値より大きければ、当該注目セルの振幅値が不要波であると判定する処理を実施する。
また、ＴＢＤ処理部７２は注目セルの振幅値が不要波でないと判定すれば、図２のＴＢＤ処理部１１と同様に、ドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算し、加算後の最終処理レンジビンにおける最大のスコアとなる周波数ビン番号をバックトラック処理部１３へ出力し、また、１レンジビン前から選択した最大のスコアとなる周波数ビン番号をＴＢＤ履歴結果保存部１２へ出力する。
一方、注目セルの振幅値が不要波であると判定すれば、加算対象の周波数ビンの振幅値を当該注目セルの振幅値に加算せずに、注目セルの加算処理後のスコアとして、注目セルの１レンジビン前のスコアを用いるようにする。なお、ＴＢＤ処理部７２はスペクトル値加算手段を構成している。

次に動作について説明する。
ただし、不要波判定用閾値設定部７１及びＴＢＤ処理部７２以外は、上記実施の形態１〜６と同様であるため、ここでは、不要波判定用閾値設定部７１及びＴＢＤ処理部７２の処理内容だけを説明する。また、ＴＢＤ処理部７２がドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算する処理については、図２のＴＢＤ処理部１１と同様であるため説明を省略する。

不要波判定用閾値設定部７１は、隣接レンジビンの１次散乱の振幅値の比の確率密度分布から、不要波判定用の閾値η_ｔｈを設定する。
即ち、不要波判定用閾値設定部７１は、事前に海域データのレンジビン毎のドップラースペクトルを取得し、レンジビン毎のドップラースペクトルから１次散乱の周波数ビンの振幅値を抽出する。
そして、不要波判定用閾値設定部７１は、（ｉ−1）番目のレンジビンの１次散乱の振幅値に対するｉ番目のレンジビンの１次散乱の振幅値の比を示す確率密度分布を求め、その確率密度分布において、例えば、分布の９９．７％となる値を不要波判定用の閾値η_ｔｈに設定する。

ここで、図２７は１次散乱の振幅特性を示す説明図である。
１次散乱の振幅特性は、図２７に示すように、レンジ方向に連続的に変化しており、隣接レンジの信号は同等の振幅を有している。このため、振幅値の比の確率密度分布とすることで、ある海域の１次散乱の振幅比の分布は他の海域でも同等になり、海域毎に確率密度分布を抽出する必要がない。また、ここでは振幅比として説明したが、振幅差として計算するようにしてもよい。

ＴＢＤ処理部７２は、１レンジビン前の探索範囲内で最大のスコアを注目セルの振幅値に加算する前に、その注目セルの振幅値が不要波であるか否かを判定し、その注目セルの振幅値が不要波であれば、加算処理を行わず、その注目セルの加算処理後のスコアとして、その注目セルの１レンジビン前のスコアを用いるようにする。
図２８はこの発明の実施の形態７のＴＢＤ処理部７２による不要波の判定処理を示す説明図である。
以下、ＴＢＤ処理部７２による不要波の判定処理を具体的に説明する。

ＴＢＤ処理部７２による不要波の判定処理では、注目セル（ｉ，ｊ）のスコアＳｊ（ｉ）の算出をする前に、下記の式（１３）に示すように、注目セル（ｉ，ｊ）の１レンジビン前の最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）に対する注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）の比η_ｊ（ｉ）を算出する。

ただし、ｊ’_ｍａｘはＴＢＤ履歴結果保存部１２により格納されている注目セル（ｉ，ｊ）の１レンジビン前の最大のスコアとなるセルの周波数ビンを示している。ここでは、下記の式（１４）のように定義する。

ＴＢＤ処理部７２は、振幅値の比η_ｊ（ｉ）を算出すると、その比η_ｊ（ｉ）と不要波判定用閾値設定部７１により設定された不要波判定用の閾値η_ｔｈを比較し、下記の式（１５）に示すように、その比η_ｊ（ｉ）が不要波判定用の閾値η_ｔｈより大きければ、注目セル（ｉ，ｊ）の振幅値が不要波であると判定する。

ＴＢＤ処理部７２は、注目セル（ｉ，ｊ）の振幅値が不要波であると判定すると、１レンジビン前の最大のスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算せず、下記の式（１６）に示すように、注目セル（ｉ，ｊ）の加算処理後のスコアＳ_ｊ（ｉ）として、注目セル（ｉ，ｊ）の１レンジビン前のスコアＳ_ｊ（ｉ−１）を用いるようにする。

以上で明らかなように、この実施の形態７によれば、ＴＢＤ処理部７２が、注目セル（ｉ，ｊ）の１レンジビン前の最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）に対する注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）の比η_ｊ（ｉ）を算出して、その比η_ｊ（ｉ）が不要波判定用の閾値η_ｔｈより大きければ、注目セル（ｉ，ｊ）の振幅値が不要波であると判定し、注目セル（ｉ，ｊ）の振幅値が不要波であると判定すれば、１レンジビン前の最大のスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算せずに、注目セル（ｉ，ｊ）の加算処理後のスコアＳ_ｊ（ｉ）として、注目セル（ｉ，ｊ）の１レンジビン前のスコアＳ_ｊ（ｉ−１）を用いるように構成したので、船等による不要波の振幅値が１次散乱の振幅値よりも非常に大きい場合でも、海面の波の流速を高精度に推定することができる効果を奏する。
また、１次散乱の振幅はレンジ方向に連続的に変化しているので、ある海域の１次散乱の振幅比の確率密度分布から不要波判定用の閾値η_ｔｈを設定することで、他の海域でも、その閾値η_ｔｈを用いることができる。

実施の形態８．
上記実施の形態７では、注目セル（ｉ，ｊ）の１レンジビン前の最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）に対する注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）の比η_ｊ（ｉ）を算出して、その比η_ｊ（ｉ）と不要波判定用の閾値η_ｔｈを比較するものを示したが、注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）は遠距離であるほど、ランダム的なノイズによって揺らいでいるため、１次散乱の振幅比の確率密度分布が広くなる。このため、振幅値の小さな不要波を判定することができずに、その不要波の周波数ビンを誤って抽出してしまって、海面の波の流速推定に誤差が生じることがある。

そこで、この実施の形態８では、注目セル（ｉ，ｊ）の１レンジビン前から事前に設定されたＬレンジビン前までの最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）〜ｘ_{ｊ’ｍａｘ}（ｉ−Ｌ）の平均値を算出して、その平均値に対する注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）の比η_ｊ（ｉ）を算出し、その比η_ｊ（ｉ）が不要波判定用の閾値η_ｔｈより大きければ、１レンジビン前の最大のスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算せずに、注目セル（ｉ，ｊ）の加算処理後のスコアＳ_ｊ（ｉ）として、注目セル（ｉ，ｊ）の１レンジビン前のスコアＳ_ｊ（ｉ−１）を用いるようにすることで、正常な流速推定を行えるようにする。

図２９はこの発明の実施の形態８によるレーダ装置の２次元ピーク検出処理部３を示す構成図であり、図２９において、図２６と同一符号は同一または相当部分を示すので説明を省略する。
ＴＢＤ処理部８１は信号処理部２により算出されたドップラースペクトルを構成する各周波数ビンを注目セル（ｉ，ｊ）として、注目セル（ｉ，ｊ）の１レンジビン前から事前に設定されたＬレンジビン前までの最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）〜ｘ_{ｊ’ｍａｘ}（ｉ−Ｌ）の平均値を算出して、その平均値に対する注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）の比η_ｊ（ｉ）（注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）と当該平均値の比）を算出し、その比η_ｊ（ｉ）が不要波判定用閾値設定部７１により設定された不要波判定用の閾値η_ｔｈより大きければ、注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）が不要波であると判定する処理を実施する。

また、ＴＢＤ処理部８１は注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）が不要波でないと判定すれば、図２のＴＢＤ処理部１１と同様に、ドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算し、加算後の最終処理レンジビンにおける最大のスコアとなる周波数ビン番号をバックトラック処理部１３へ出力し、また、１レンジビン前から選択した最大のスコアとなる周波数ビン番号をＴＢＤ履歴結果保存部１２へ出力する。
一方、注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）が不要波であると判定すれば、１レンジビン前の最大のスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算せずに、注目セル（ｉ，ｊ）の加算処理後のスコアＳ_ｊ（ｉ）として、注目セル（ｉ，ｊ）の１レンジビン前のスコアＳ_ｊ（ｉ−１）を用いるようにする。なお、ＴＢＤ処理部８１はスペクトル値加算手段を構成している。

次に動作について説明する。
ただし、ＴＢＤ処理部８１以外は、上記実施の形態７と同様であるため、ここでは、ＴＢＤ処理部８１の処理内容だけを説明する。また、ＴＢＤ処理部８１がドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算する処理については、図２のＴＢＤ処理部１１と同様であるため説明を省略する。

図３０はこの発明の実施の形態８のＴＢＤ処理部８１による不要波の判定処理を示す説明図である。
ＴＢＤ処理部８１による不要波の判定処理では、注目セル（ｉ，ｊ）のスコアＳｊ（ｉ）の算出をする前に、注目セル（ｉ，ｊ）の１レンジビン前から事前に設定されたＬレンジビン前までの最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）〜ｘ_{ｊ’ｍａｘ}（ｉ−Ｌ）の平均値ｘ_{ｊ’ａｖｅ}を算出して、その平均値ｘ_{ｊ’ａｖｅ}に対する注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）の比η_ｊ（ｉ）を算出する。

ＴＢＤ処理部８１は、振幅値の比η_ｊ（ｉ）を算出すると、図２６のＴＢＤ処理部７２と同様に、その比η_ｊ（ｉ）と不要波判定用閾値設定部７１により設定された不要波判定用の閾値η_ｔｈを比較し、その比η_ｊ（ｉ）が不要波判定用の閾値η_ｔｈより大きければ、注目セル（ｉ，ｊ）の振幅値が不要波であると判定する。
ＴＢＤ処理部８１は、注目セル（ｉ，ｊ）の振幅値が不要波であると判定すると、図２６のＴＢＤ処理部７２と同様に、１レンジビン前の最大のスコアを注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）に加算せず、注目セル（ｉ，ｊ）の加算処理後のスコアＳ_ｊ（ｉ）として、注目セル（ｉ，ｊ）の１レンジビン前のスコアＳ_ｊ（ｉ−１）を用いるようにする。

ここで、図３１は平均化による１次散乱の振幅比の確率密度分布を示す説明図である。
注目セル（ｉ，ｊ）の１レンジビン前からＬレンジビン前までの最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）〜ｘ_{ｊ’ｍａｘ}（ｉ−Ｌ）を平均化することで、１次散乱の確率密度分布と不要波の確率密度分布が分離することができるため、不要波の検出を減らすことができる。
１次散乱の確率密度分布を平均化していない場合、図３１（ａ）のように、１次散乱の確率密度分布と不要波の確率密度分布が一部重なっているため、両者を分離することができないが、１次散乱の確率密度分布を平均化している場合、図３１（ｂ）のように、１次散乱の確率密度分布と不要波の確率密度分布が重なっていないため、両者を分離することができる。

以上で明らかなように、この実施の形態８によれば、ＴＢＤ処理部８１が、注目セル（ｉ，ｊ）の１レンジビン前から事前に設定されたＬレンジビン前までの最大のスコアとなるセルの振幅値ｘ_{ｊ’ｍａｘ}（ｉ−１）〜ｘ_{ｊ’ｍａｘ}（ｉ−Ｌ）の平均値を算出して、その平均値に対する注目セル（ｉ，ｊ）の振幅値ｘ_ｊ（ｉ）の比η_ｊ（ｉ）を算出するように構成したので、振幅値の小さな不要波が含まれている場合でも、海面の波の流速を高精度に推定することができる効果を奏する。

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

１ ビーム送受信系、２ 信号処理部（スペクトル算出手段）、３ ２次元ピーク検出処理部（速度推定手段）、４ 出力データ保存部、１１ ＴＢＤ処理部（スペクトル値加算手段）、１２ ＴＢＤ履歴結果保存部、１３ バックトラック処理部（速度算出手段）、１４ 平滑処理部（速度算出手段）、２１ 相互相関処理部（スペクトル値加算手段）、２２ オフセット量算出部（スペクトル値加算手段）、２３ ＴＢＤ処理部（スペクトル値加算手段）、３１ 探索領域決定部（スペクトル値加算手段）、３２，３３ ＴＢＤ処理部（スペクトル値加算手段）、４１ 異常時処理部（スペクトル値加算手段）、４２ ＴＢＤ処理部（スペクトル値加算手段）、５１ 波形推定処理部（速度偏差算出手段）、５２ 流速偏差抽出部（速度偏差算出手段）、５３ 異常ピーク判定部（正否判定手段）、６１ 異常ピーク候補検出部、６２ 異常ピーク判定処理部、７１ 不要波判定用閾値設定部（閾値設定手段）、７２ ＴＢＤ処理部（スペクトル値加算手段）、８１ ＴＢＤ処理部（スペクトル値加算手段）。

Claims (13)

1. 観測対象に反射されたビームの反射波を周波数解析することで、レンジビン毎にドップラースペクトルを算出するスペクトル算出手段と、
   前記スペクトル算出手段により算出されたドップラースペクトルのスペクトル値をレンジ方向に繰り返し加算し、加算処理後のドップラースペクトルのスペクトル値から前記観測対象の速度を推定する速度推定手段とを備えたレーダ装置において、
   前記速度推定手段は、前記スペクトル算出手段により算出されたドップラースペクトルを構成する各々の周波数ビンを注目セルとして、当該注目セルのレンジビンより１レンジビン前のドップラースペクトルの中に、レンジ方向の１次散乱ピークのばらつき分に相当する幅を有する探索領域を設定し、前記探索領域内の複数の周波数ビンの中で、スペクトル値が最も大きい周波数ビンを、加算対象の周波数ビンとして選択して、前記加算対象の周波数ビンにおけるスペクトル値の累積加算値であるスコア値を、前記注目セルとしている周波数ビンのスペクトル値に加算し、その加算結果を前記注目セルとしている周波数ビンのスコア値とするスペクトル値加算手段を備えたことを特徴とするレーダ装置。
2. 前記速度推定手段は、前記スペクトル値加算手段による加算処理後のドップラースペクトルにおける各々の周波数ビンのスコア値を比較して、前記観測対象の速度に対応する周波数ビンを特定し、当該周波数ビンから前記観測対象の速度を算出する速度算出手段を備えたことを特徴とする請求項１記載のレーダ装置。
3. 前記速度算出手段は、前記注目セルと前記加算対象の周波数ビンとの対応関係に基づいて、レンジビン毎に前記観測対象の速度に対応する周波数ビンを特定して、当該周波数ビンから前記観測対象の速度を算出し、レンジビン毎に算出した速度を平滑化することを特徴とする請求項２記載のレーダ装置。
4. 前記スペクトル値加算手段は、前記スペクトル算出手段により算出された各レンジビンのドップラースペクトルからレンジビン間の相互相関値を算出して、前記相互相関値の相関値幅を算出し、前記相関値幅に対して予め設定された係数が乗算された値より前記相互相関値が大きくなる周波数ビンの範囲を特定して、前記周波数ビンの範囲を前記探索領域に決定することを特徴とする請求項１から請求項３のうちのいずれか１項記載のレーダ装置。
5. 前記スペクトル値加算手段は、前記スペクトル算出手段により算出された各レンジビンのドップラースペクトルからレンジビン間の相互相関値を算出し、前記相互相関値に基づいて前記探索領域のオフセット量を算出し、前記オフセット量だけ前記探索領域をずらすことを特徴とする請求項１から請求項４のうちのいずれか１項記載のレーダ装置。
6. 前記スペクトル値加算手段は、前記レンジビン間の相互相関値における相関ピークと相関ボトムの比が予め設定された第１の閾値より小さい場合、異常時の処理として、前記オフセット量を補正することを特徴とする請求項５記載のレーダ装置。
7. 前記スペクトル値加算手段は、相関対象のレンジビンのドップラースペクトルを前記オフセット量だけずらし、前記オフセット量だけずらした相関対象のレンジビンのドップラースペクトルと、相関元のレンジビンのドップラースペクトルとの相関係数が予め設定された第２の閾値より小さい場合、異常時の処理として、前記オフセット量を補正することを特徴とする請求項５記載のレーダ装置。
8. 前記スペクトル値加算手段は、各レンジビンのオフセット量の絶対値と予め設定された第３の閾値を比較し、当該レンジビンのオフセット量の絶対値が前記第３の閾値より大きく、かつ、当該レンジビンと隣接しているレンジビンのオフセット量の絶対値が予め設定された第４の閾値より小さい場合、当該レンジビンのオフセット量を前記隣接しているレンジビンのオフセット量を用いて補正することを特徴とする請求項５から請求項７のうちのいずれか１項記載のレーダ装置。
9. 前記速度推定手段により推定された速度である観測値から、前記速度の予測値を算出して、前記観測値と前記予測値との差分値を算出する速度偏差算出手段と、
   前記速度偏差算出手段により算出された差分値から前記観測値の正否を判定する正否判定手段とを備えたことを特徴とする請求項１から請求項８のうちのいずれか１項記載のレーダ装置。
10. 前記正否判定手段は、各ビームのレンジビン毎に、前記速度偏差算出手段により算出された差分値と予め設定された第５の閾値を比較して、前記差分値が前記第５の閾値より大きいレンジビンを異常ピーク候補として検出し、前記異常ピーク候補を有するビームと隣接しているビームが有する異常ピーク候補の数が予め設定された第６の閾値より少なければ、前記第５の閾値より大きい前記差分値の算出元の観測値は異常値であると判定することを特徴とする請求項９記載のレーダ装置。
11. 前記スペクトル値加算手段は、前記注目セルとしている周波数ビンのスペクトル値と前記加算対象の周波数ビンのスペクトル値との比を算出して、前記比と予め設定された第７の閾値を比較し、前記比が前記第７の閾値より大きければ、前記加算対象の周波数ビンのスコア値を前記注目セルとしている周波数ビンのスペクトル値に加算せずに、前記注目セルとしている周波数ビンのスコア値として、前記注目セルとしている周波数ビンの１レンジビン前の加算処理後のスコア値を用いることを特徴とする請求項１から請求項１０のうちのいずれか１項記載のレーダ装置。
12. 前記スペクトル値加算手段は、前記注目セルとしている周波数ビンの１レンジビン前から事前に設定されたレンジビン前までの前記加算対象の周波数ビンのスペクトル値の平均値を算出し、前記注目セルとしている周波数ビンのスペクトル値と前記平均値の比を算出して、前記比と予め設定された第７の閾値を比較し、前記比が前記第７の閾値より大きければ、前記加算対象の周波数ビンのスコア値を前記注目セルとしている周波数ビンのスペクトル値に加算せずに、前記注目セルとしている周波数ビンのスコア値として、前記注目セルとしている周波数ビンの１レンジビン前の加算処理後のスコア値を用いることを特徴とする請求項１から請求項１０のうちのいずれか１項記載のレーダ装置。
13. 隣接レンジビンの１次散乱の振幅値の比の確率密度分布から、前記第７の閾値を設定する閾値設定手段を備えたことを特徴とする請求項１１または請求項１２記載のレーダ装置。

Images