JP6328329B2 - 二重偏波レーダ装置及びレーダ信号処理方法 - Google Patents

Description

この発明は、水平偏波と垂直偏波を用いて、降雨強度やドップラー速度を推定する二重偏波レーダ装置及びレーダ信号処理方法に関するものである。

広く用いられている単一偏波の気象レーダでは、エコー強度と降雨強度の間に成立する関係式を用いて、エコー強度から降雨強度を推定している。しかし、この関係式は、雨滴の粒径分布によって大きく変化するため、単一偏波の気象レーダで正確な降雨強度を推定することは困難である。
一般に雨滴は、球を垂直方向に押し潰した回転楕円体に近い形状をなしている。大きな粒径の雨滴は扁平度が大きく、小さな粒径の雨滴は扁平度が小さくほぼ球形となる。このため、水平偏波と垂直偏波の双方を用いる二重偏波レーダ観測を行えば、雨滴の形状に関する情報として、粒径分布に対応する情報が得られるため、精度の高い降雨強度の計測が可能になる。

二重偏波レーダにおける偏波の送受信において、送信偏波と同一偏波の後方散乱電力と、送信偏波と直交する偏波によって生じる後方散乱電力との比は、交差偏波比（ＬＤＲ）と呼ばれる。
交差偏波比は、氷粒の形状や方向などの定性的な判別の指標となるものであり、ブライトバンドなどにおいて、降水粒子の中で送信偏波と直交する偏波の散乱（Ｚ_ｈｖ，Ｚ_ｖｈ）が比較的強い場合に有効なパラメータとなる。

ここで、交互送信交互受信の運用で交差偏波比を観測する場合、下記に示す３種類の偏波送受信処理を繰り返す必要がある。
（１）水平偏波を送信して水平偏波を受信する送受信処理である「水平偏波送信水平偏波受信」
（２）垂直偏波を送信して垂直偏波を受信する送受信処理である「垂直偏波送信垂直偏波受信」
（３）水平偏波を送信して垂直偏波を受信する送受信処理である「水平偏波送信垂直偏波受信」、あるいは、垂直偏波を送信して水平偏波を受信する送受信処理である「垂直偏波送信水平偏波受信」

３種類の偏波送受信処理を繰り返す場合、偏波の観測時間が、単偏波の送受信処理や、複数の偏波の同時送信同時受信の送受信処理を行う場合と比べて３倍に増加することから、一般的には、ドップラー速度を観測することが可能な範囲であるナイキスト速度が、単偏波の送受信処理等を行う場合と比べて１／３倍に減少することが想定される。
なお、以下の非特許文献１には、「水平偏波送信水平偏波受信」と「垂直偏波送信垂直偏波受信」の２種類の偏波送受信処理を繰り返す場合において、ナイキスト速度が、単偏波の送受信処理等を行う場合と比べて減少しないようにするための方法が開示されている。

Ｍ. ＳＡＣＨＩＤＡＮＡＮＤＡ ＡＮＤ Ｄ. Ｓ. ＺＲＮＩＣ, Ｅｆｆｉｃｉｅｎｔ Ｐｒｏｃｅｓｓｉｎｇ ｏｆ Ａｌｔｅｒｎａｔｅｌｙ Ｐｏｌａｒｉｚｅｄ Ｒａｄａｒ Ｓｉｇｎａｌｓ, １９８９．

従来の二重偏波レーダ装置は以上のように構成されているので、非特許文献１に開示されている方法を用いれば、「水平偏波送信水平偏波受信」と「垂直偏波送信垂直偏波受信」の２種類の偏波送受信処理を繰り返す場合に、ナイキスト速度が、単偏波の送受信処理等を行う場合と比べて減少しないようにすることができる。しかし、非特許文献１に開示されている方法では、３種類の偏波送受信処理を繰り返す場合を想定しておらず、交差偏波比を算出することができないという課題があった。

この発明は上記のような課題を解決するためになされたもので、３種類の偏波送受信処理を繰り返す場合でも、ナイキスト速度が減少することなく、交差偏波比を算出することができる二重偏波レーダ装置及びレーダ信号処理方法を得ることを目的とする。

この発明に係る二重偏波レーダ装置は、水平偏波を送信して水平偏波を受信する第１の送受信処理と、垂直偏波を送信して垂直偏波を受信する第２の送受信処理と、水平偏波を送信して垂直偏波を受信する第３の送受信処理又は垂直偏波を送信して水平偏波を受信する第４の送受信処理とを順番に繰り返し実施する送受信装置と、送受信装置により第１の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第１の反射波強度を算出し、第２の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第２の反射波強度を算出し、第３の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第３の反射波強度又は第４の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第４の反射波強度を算出する反射波強度算出部と、反射波強度算出部により算出された第１の反射波強度と第３の反射波強度との比、あるいは、第２の反射波強度と第４の反射波強度との比である交差偏波比を算出する交差偏波比算出部と、反射波強度算出部により算出された第１の反射波強度と、第２の反射波強度と、第３の反射波強度又は第４の反射波強度とから、送受信装置により偏波が繰り返し送信される周期での各送信偏波における自己相関のラグを算出して、ラグから送受信装置により繰り返し送信される偏波間の位相差を算出し、ラグ、偏波間の位相差及びドップラー速度を観測することが可能な範囲であるナイキスト速度から、偏波を反射している観測対象のドップラー速度を算出する速度算出部と、を備えたものである。

この発明によれば、交差偏波比算出部が、反射波強度算出部により算出された第１の反射波強度と第３の反射波強度との比、あるいは、第２の反射波強度と第４の反射波強度との比である交差偏波比を算出するように構成したので、３種類の偏波送受信処理を繰り返す場合でも、ナイキスト速度が減少することなく、交差偏波比を算出することができる効果がある。

この発明の実施の形態１による二重偏波レーダ装置を示す構成図である。 この発明の実施の形態１による二重偏波レーダ装置の信号処理装置６を示す構成図である。 図２の信号処理装置６のハードウェア構成図である。 信号処理装置６がコンピュータで構成される場合のハードウェア構成図である。 この発明の実施の形態１による二重偏波レーダ装置の処理内容であるレーダ信号処理方法を示すフローチャートである。 送受信装置１により送受信される各偏波の間隔がＴ_ｓである例を示す説明図である。 式（１８）、式（２０）及び式（２２）から分かるＲハット_ａ（Ｔ_ｓ），Ｒハット_ｂ（Ｔ_ｓ），Ｒハット_ｃ（Ｔ_ｓ）の位相の関係を示す説明図である。 送受信装置１により送受信される各偏波の間隔がＴ_ｓである例を示す説明図である。 式（４５）、式（４７）、及び式（４９）及び式（５１）から分かるＲハット_ａ（Ｔ_ｓ），Ｒハット_ｂ（Ｔ_ｓ），Ｒハット_ｃ（Ｔ_ｓ），Ｒハット_ｄ（Ｔ_ｓ）の位相の関係を示す説明図である。 この発明の実施の形態３による二重偏波レーダ装置の信号処理装置６を示す構成図である。 図１０の信号処理装置６のハードウェア構成図である。

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

実施の形態１．
図１はこの発明の実施の形態１による二重偏波レーダ装置を示す構成図である。
また、図２はこの発明の実施の形態１による二重偏波レーダ装置の信号処理装置６を示す構成図であり、図３は図２の信号処理装置６のハードウェア構成図である。
図１〜３において、送受信装置１は送信機２、送受切替器３、アンテナ４及び受信機５から構成されており、下記に示す３種類の偏波送受信処理を繰り返し実施する。
（１）水平偏波を送信して水平偏波を受信する第１の送受信処理である「水平偏波送信水平偏波受信」
（２）垂直偏波を送信して垂直偏波を受信する第２の送受信処理である「垂直偏波送信垂直偏波受信」
（３）水平偏波を送信して垂直偏波を受信する第３の送受信処理である「水平偏波送信垂直偏波受信」、あるいは、垂直偏波を送信して水平偏波を受信する第４の送受信処理である「垂直偏波送信水平偏波受信」

送信機２は水平偏波又は垂直偏波を送受切替器３に出力する装置である。
送受切替器３は送信機２から出力された水平偏波又は垂直偏波をアンテナ４に出力し、アンテナ４により受信された水平偏波又は垂直偏波を受信機５に出力する装置である。
アンテナ４は送受切替器３から出力された水平偏波又は垂直偏波を空間に向けて放射する一方、空間に存在している雨滴等の観測対象に反射されて戻ってきた水平偏波又は垂直偏波を受信する。
受信機５はアンテナ４により受信された水平偏波又は垂直偏波の受信処理を実施し、その水平偏波又は垂直偏波の受信信号をアナログ信号からデジタル信号に変換して、デジタルの受信信号を信号処理装置６に出力する装置である。

信号処理装置６は反射波強度算出部１１、交差偏波比算出部１２及び速度算出部１３から構成されており、送受信装置１から出力された水平偏波と垂直偏波のデジタル受信信号から、水平偏波及び垂直偏波の電力値である反射波強度を算出して、その反射波強度から交差偏波比ＬＤＲを算出するとともに、観測対象のドップラー速度ｖハットを算出する装置である。明細書の文書中では、電子出願の関係上、“ｖ”の文字の上に“＾”の記号を付することができないので、便宜上、“ｖハット”のように表記している。
表示装置７は例えば液晶ディスプレイなどから構成されており、信号処理装置６により算出された交差偏波比ＬＤＲやドップラー速度ｖハットなどを表示する装置である。

反射波強度算出部１１は送受信装置１から出力されたデジタル受信信号を入力する信号入力インタフェースである信号入力処理回路２１と、例えばＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）を実装している半導体処理回路あるいはワンチップマイコンなどから構成されている反射波強度算出処理回路２２と、表示装置７に対するインタフェースである表示処理回路２５とから実現されるものである。
反射波強度算出部１１は送受信装置１により第１の送受信処理で受信された水平偏波のデジタル受信信号から、当該水平偏波の電力値である反射波強度Ｖ_ｈｈ（第１の反射波強度）を算出するとともに、第２の送受信処理で受信された垂直偏波のデジタル受信信号から、当該垂直偏波の電力値である反射波強度Ｖ_ｖｖ（第２の反射波強度）を算出する処理を実施する。
また、反射波強度算出部１１は第３の送受信処理で受信された垂直偏波のデジタル受信信号から、当該垂直偏波の電力値である反射波強度Ｖ_ｖｈ（第３の反射波強度）の算出処理、あるいは、第４の送受信処理で受信された水平偏波のデジタル受信信号から、当該水平偏波の電力値である反射波強度Ｖ_ｈｖ（第４の反射波強度）の算出処理を実施する。

また、反射波強度算出部１１は反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｖｈ（またはＶ_ｈｖ）をそれぞれ積分することで、水平偏波送信水平偏波受信でのレーダ反射因子Ｚ_ｈｈ、垂直偏波送信垂直偏波受信でのレーダ反射因子Ｚ_ｖｖ、水平偏波送信垂直偏波受信でのレーダ反射因子Ｚ_ｖｈ（または垂直偏波送信水平偏波受信でのレーダ反射因子Ｚ_ｈｖ）を算出するとともに、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｖとの比であるレーダ反射因子差Ｚｄｒを算出する処理を実施する。
さらに、反射波強度算出部１１は算出したレーダ反射因子Ｚ_ｈｈ，Ｖ_ｖｖ，Ｚ_ｖｈ（またはＺ_ｈｖ）と、レーダ反射因子差Ｚｄｒとを表示装置７に表示する処理を実施する。

交差偏波比算出部１２は例えばＣＰＵを実装している半導体処理回路あるいはワンチップマイコンなどから構成されている交差偏波比算出処理回路２３と、表示処理回路２５とから実現されるものであり、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈと反射波強度Ｖ_ｖｈとの比、即ち、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈとの比である交差偏波比ＬＤＲ_ｖｈを算出する。あるいは、反射波強度算出部１１により算出された反射波強度Ｖ_ｖｖと反射波強度Ｖ_ｈｖとの比、即ち、レーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖとの比である交差偏波比ＬＤＲ_ｈｖを算出する処理を実施する。
また、交差偏波比算出部１２は算出した交差偏波比ＬＤＲ_ｖｈ、あるいは、交差偏波比ＬＤＲ_ｈｖを表示装置７に表示する処理を実施する。

速度算出部１３は例えばＣＰＵを実装している半導体処理回路あるいはワンチップマイコンなどから構成されている速度算出処理回路２４、表示処理回路２５とから実現されるものであり、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈと、反射波強度Ｖ_ｖｖと、反射波強度Ｖ_ｖｈ又は反射波強度Ｖ_ｈｖとから、パルスペア法を用いて、観測対象のドップラー速度ｖハットを算出する処理を実施する。
また、速度算出部１３は算出したドップラー速度ｖハットを表示装置７に表示する処理を実施する。

図１の例では、信号処理装置６の構成要素である反射波強度算出部１１、交差偏波比算出部１２及び速度算出部１３のそれぞれが専用のハードウェアで構成されているものを想定しているが、信号処理装置６がコンピュータで構成されているものであってもよい。
図４は信号処理装置６がコンピュータで構成される場合のハードウェア構成図である。
信号処理装置６がコンピュータで構成される場合、反射波強度算出部１１、交差偏波比算出部１２及び速度算出部１３の処理内容を記述しているプログラムをコンピュータのメモリ３１に格納し、コンピュータのプロセッサ３２がメモリ３１に格納されているプログラムを実行するようにすればよい。
図５はこの発明の実施の形態１による二重偏波レーダ装置の処理内容であるレーダ信号処理方法を示すフローチャートである。
図４の例では、プロセッサ３２が１個の例を示しているが、コンピュータが複数のプロセッサ３２を実装しているものであってもよい。

次に動作について説明する。
図６は送受信装置１により送受信される各偏波の間隔がＴ_ｓである例を示す説明図である。
図６の例では、第１の送受信処理である「水平偏波送信水平偏波受信」が（３ｉ−２）・Ｔ_ｓの時間に行われ、第２の送受信処理である「垂直偏波送信垂直偏波受信」が（３ｉ−１）・Ｔ_ｓの時間に行われ、第３の送受信処理である「水平偏波送信垂直偏波受信」又は第４の送受信処理である「垂直偏波送信水平偏波受信」が３ｉ・Ｔ_ｓの時間に行われている。
ただし、ｉ＝１，２，３，・・・，Ｉである。
Ｉはサンプルペア数（水平偏波送信水平偏波受信、垂直偏波送信垂直偏波受信、水平偏波送信垂直偏波受信（または垂直偏波送信水平偏波受信）それぞれのヒット数）である。

この実施の形態１では、送受信装置１が、第１の送受信処理である「水平偏波送信水平偏波受信」と、第２の送受信処理である「垂直偏波送信垂直偏波受信」と、第３の送受信処理である「水平偏波送信垂直偏波受信」とを実施するものとする。
図６の例では、実施順序は、「水平偏波送信水平偏波受信」→「垂直偏波送信垂直偏波受信」→「水平偏波送信垂直偏波受信」の順番である。
この実施の形態１では、二重偏波レーダ装置の観測対象が雨滴である例を説明する。

送受信装置１の送信機２は、「水平偏波送信水平偏波受信」→「垂直偏波送信垂直偏波受信」→「水平偏波送信垂直偏波受信」の順番で送受信処理を実施するため、送信周期Ｔ_ｓで、水平偏波→垂直偏波→水平偏波→水平偏波→垂直偏波→水平偏波→・・・→水平偏波→垂直偏波→水平偏波を送信する。
送受切替器３は、送信機２から水平偏波又は垂直偏波を受けると、その水平偏波又は垂直偏波をアンテナ４に出力する。
アンテナ４は、送受切替器３から水平偏波又は垂直偏波を受けると、その水平偏波又は垂直偏波を空間に向けて放射する（図５のステップＳＴ１）。
アンテナ４は、水平偏波又は垂直偏波を空間に向けて放射したのち、観測対象である雨滴に反射されて戻ってきた水平偏波又は垂直偏波を受信する。
図６の例では、送信周期Ｔ_ｓで、水平偏波→垂直偏波→垂直偏波→水平偏波→垂直偏波→垂直偏波→・・・→水平偏波→垂直偏波→垂直偏波を受信する。

送受切替器３は、アンテナ４により受信された水平偏波又は垂直偏波を受信機５に出力する。
受信機５は、アンテナ４により受信された水平偏波又は垂直偏波の受信処理を実施し、その水平偏波又は垂直偏波の受信信号をアナログ信号からデジタル信号に変換して、デジタルの受信信号を信号処理装置６に出力する（ステップＳＴ２）。

信号処理装置６の反射波強度算出部１１は、送受信装置１により第１の送受信処理で受信された水平偏波のデジタル受信信号から、当該水平偏波の電力値である反射波強度Ｖ_ｈｈを算出するとともに、第２の送受信処理で受信された垂直偏波のデジタル受信信号から、当該垂直偏波の電力値である反射波強度Ｖ_ｖｖを算出する。また、第３の送受信処理で受信された垂直偏波のデジタル受信信号から、当該垂直偏波の電力値である反射波強度Ｖ_ｖｈを算出する（ステップＳＴ３）。
以下、反射波強度算出部１１による反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｖｈの算出処理を具体的に説明する。

一般的に、観測対象である雨滴は、空間にランダムに存在しているので、送受信装置１から出力されるデジタル受信信号もランダムな性質を有している。
そのため、反射波強度算出部１１が、反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｖｈを算出する際、複数回の偏波の送信（パルスの送信）によって得られた複数のデジタル受信信号（統計的に独立な複数の受信信号）を平均化処理（積分処理）することで、後述する反射波強度の推定値のゆらぎを抑圧して反射波強度の精度を向上させることが望ましい。ただし、ここでの平均化処理の対象のデジタル受信信号は、複数回の第１の送受信処理で受信された水平偏波のデジタル受信信号同士、複数回の第２の送受信処理で受信された垂直偏波のデジタル受信信号同士であり、また、複数回の第３の送受信処理で受信された垂直偏波のデジタル受信信号同士である。

ここで、アンテナ４の受信信号である後方散乱受信信号の電圧行列と、偏波の送信行列及び後方散乱行列の関係式は、一般的に下記の式（１）〜（３）で表される。

反射波強度算出部１１は、複数回の第１の送受信処理で受信された水平偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、ｎ番目のサンプルにおける反射波強度Ｖ_ｈｈ（ｎ）を下記の式（４）のように算出する。
また、反射波強度算出部１１は、複数回の第２の送受信処理で受信された垂直偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、（ｎ＋１）番目のサンプルにおける反射波強度Ｖ_ｖｖ（ｎ＋１）を下記の式（５）のように算出する。
さらに、反射波強度算出部１１は、複数回の第３の送受信処理で受信された垂直偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、（ｎ＋２）番目のサンプルにおける反射波強度Ｖ_ｖｈ（ｎ＋２）を下記の式（６）のように算出する。なお、これらの反射波強度は、観測領域の方位角と距離の区分（格子）毎に算出される。

ここで、距離ｒ^ｉ（ｎ＋１）は、一つ前の距離ｒ^ｉ（ｎ）にＴ_ｓｖ_ｉだけプラスして移動したものであるため、下記の式（７）のように表される。

反射波強度算出部１１は、反射波強度Ｖ_ｈｈ（ｎ），Ｖ_ｖｖ（ｎ＋１），Ｖ_ｖｈ（ｎ＋２）を算出すると、その反射波強度Ｖ_ｈｈ（ｎ），Ｖ_ｖｖ（ｎ＋１），Ｖ_ｖｈ（ｎ＋２）のそれぞれを積分することで、反射波強度Ｖ_ｈｈ（ｎ）の積分後の反射波強度である水平偏波送信水平偏波受信でのレーダ反射因子Ｚ_ｈｈと、反射波強度Ｖ_ｖｖ（ｎ＋１）の積分後の反射波強度である垂直偏波送信垂直偏波受信でのレーダ反射因子Ｚ_ｖｖと、反射波強度Ｖ_ｖｈ（ｎ＋２）の積分後の反射波強度である水平偏波送信垂直偏波受信でのレーダ反射因子Ｚ_ｖｈとを求める。
この積分範囲は、レーダ装置の精度等を考慮して、ユーザが任意に決定することができるが、例えば、送受信装置１により送受信される各偏波の間隔Ｔ_ｓと、送受信装置１による既知のビーム範囲（レーダを回転させる範囲）とから決定されることが想定される。

また、反射波強度算出部１１は、下記の式（８）に示すように、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｖとの比であるレーダ反射因子差Ｚｄｒを算出する。

さらに、反射波強度算出部１１は、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈを交差偏波比算出部１２に出力するとともに、レーダ反射因子Ｚ_ｈｈ，Ｚ_ｖｖ，Ｚ_ｖｈ及びレーダ反射因子差Ｚｄｒを表示装置７に表示する。

交差偏波比算出部１２は、反射波強度算出部１１からレーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈを受けると、そのレーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈを下記の式（９）に代入することで、交差偏波比ＬＤＲ_ｖｈを算出する（ステップＳＴ４）。

交差偏波比算出部１２は、交差偏波比ＬＤＲ_ｖｈを算出すると、その交差偏波比ＬＤＲ_ｖｈを表示装置７に表示する。

速度算出部１３は、反射波強度算出部１１が反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｖｈを算出すると、その反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｖｈから、パルスペア法を用いて、観測対象のドップラー速度ｖハットを算出する（ステップＳＴ５）。
即ち、速度算出部１３は、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｖｈから、送受信装置１により偏波が繰り返し送信される周期Ｔ_ｓでの各送信偏波における自己相関のラグＲハット（Ｔ_ｓ）を算出して、そのラグＲハット（Ｔ_ｓ）から送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰを算出し、そのラグＲハット（Ｔ_ｓ）、偏波間の位相差φハット_ＤＰ及びドップラー速度ｖハットを観測することが可能な範囲であるナイキスト速度ｖ_ａから、観測対象のドップラー速度ｖハットを算出する。
明細書の文書中では、電子出願の関係上、“Ｒ”、“φ”の文字の上に“＾”の記号を付することができないので、便宜上、“Ｒハット”、“φハット”のように表記している。
以下、速度算出部１３によるドップラー速度ｖハットの算出処理を具体的に説明する。

偏波の交互送信における自己相関の推定法（パルスペア法）では、ラグＲハット（Ｔ_ｓ）は、下記の式（１０）のように表される。一般的には、Ｒハット（Ｔ_ｓ）は“ラグ１”と呼ばれ、Ｒハット（２Ｔ_ｓ）は“ラグ２”と呼ばれる。

式（１０）において、Ｍはサンプルペア数（水平偏波送信水平偏波受信、垂直偏波送信垂直偏波受信、水平偏波送信垂直偏波受信（または垂直偏波送信水平偏波受信）それぞれのヒット数）である。

式（１０）に示すＲハット（Ｔ_ｓ）は、下記の式（１１）〜（１３）に示すように、３つの式に分解することができる。即ち、式（１０）に示すＲハット（Ｔ_ｓ）を構成しているＲハット_ａ（Ｔ_ｓ）、Ｒハット_ｂ（Ｔ_ｓ）、Ｒハット_ｃ（Ｔ_ｓ）を特定することができる。

Ｒハット_ａ（Ｔ_ｓ）は、水平偏波受信時における水平偏波観点の共分散、即ち、第１の送受信処理で水平偏波が受信された時点での当該水平偏波と第２の送受信処理で受信された垂直偏波との共分散（第１の共分散）である。
Ｒハット_ｂ（Ｔ_ｓ）は、垂直偏波観点の共分散、即ち、第２の送受信処理で垂直偏波が受信された時点での当該垂直偏波と第３の送受信処理で受信された垂直偏波との共分散（第２の共分散）である。
Ｒハット_ｃ（Ｔ_ｓ）は、垂直偏波受信時における水平偏波観点の共分散、即ち、第３の送受信処理で垂直偏波が受信された時点での当該垂直偏波と第１の送受信処理で受信された水平偏波との共分散（第３の共分散）である。

これにより、式（１０）に示すＲハット（Ｔ_ｓ）は、下記の式（１４）のように表される。

式（１１）に示すＲハット_ａ（Ｔ_ｓ）は、式（４）に示す反射波強度Ｖ_ｈｈ（ｎ）と、式（５）に示す反射波強度Ｖ_ｖｖ（ｎ＋１）とを用いると、下記の式（１５）のように表される。

ここで、位相である２（ｋ_０＋ｋ_ｈ）ｒ^ｉ（ｎ）は、０〜２πの間に一様に分布するため、式（１５）のｅｘｐ部分の期待値はｉ＝ｌの時以外０になる。

したがって、式（１５）に示すＲハット_ａ（Ｔ_ｓ）は、下記の式（１６）のように表される。

式（１６）において、ｒ^ｉ（ｎ）の平均を取ると、全体の距離の中間になることから、ｒ^ｉ（ｎ）をｒ_０（中間の距離）とおくことができる。また、ｖ_ｉの平均をとると、全体の速度の平均となることから、ｖ_ｉをｖ（平均速度）とおくことができる。

式（１７）において、２（ｋ_ｈ−ｋ_ｖ）ｒ_０は水平偏波と垂直偏波の位相差φ_ＤＰにあたるので、２（ｋ_ｈ−ｋ_ｖ）ｒ_０＝φ_ＤＰであり、−２（ｋ_０＋ｋ_ｖ）Ｔ_ｓｖは垂直偏波のドップラシフト量ψ_ｄにあたるので、−２（ｋ_０＋ｋ_ｖ）Ｔ_ｓｖ＝ψ_ｄである。
また、前提として、水平偏波と垂直偏波のドップラシフト量は変わらないと仮定しているため、−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖ＝ψ_ｄである。
以上より、Ｒハット_ａ（Ｔ_ｓ）は、下記の式（１８）のように表される。

式（１２）に示すＲハット_ｂ（Ｔ_ｓ）は、式（５）に示す反射波強度Ｖ_ｖｖ（ｎ＋１）と、式（６）に示す反射波強度Ｖ_ｖｈ（ｎ＋２）を用いると、下記の式（１９）のように表される。

式（１９）に示すＲハット_ｂ（Ｔ_ｓ）は、
φ_ＤＰ＝２（ｋ_ｈ−ｋ_ｖ）ｒ_０，ψ_ｄ＝−２（ｋ_０＋ｋ_ｖ）Ｔ_ｓｖ，ψ_ｄ＝−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖより、下記の式（２０）のように表される。

式（１３）に示すＲハット_ｃ（Ｔ_ｓ）は、式（４）に示す反射波強度Ｖ_ｈｈ（ｎ）と、式（６）に示す反射波強度Ｖ_ｖｈ（ｎ＋２）を用いると、下記の式（２１）のように表される。

式（２１）に示すＲハット_ｃ（Ｔ_ｓ）は、
φ_ＤＰ＝２（ｋ_ｈ−ｋ_ｖ）ｒ_０，ψ_ｄ＝−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖより、下記の式（２２）のように表される。

ここで、図７は式（１８）、式（２０）及び式（２２）から分かるＲハット_ａ（Ｔ_ｓ），Ｒハット_ｂ（Ｔ_ｓ），Ｒハット_ｃ（Ｔ_ｓ）の位相の関係を示す説明図である。図７では、図面の簡単化のため、（Ｔ_ｓ）の記述を省略して、Ｒハット_ａ，Ｒハット_ｂ，Ｒハット_ｃの
ように表記している。
図７に示すように、Ｒハット_ａと、Ｒハット_ａ＋Ｒハット_ｂまたはＲハット_ａ＋Ｒハット_ｃとの位相差は、φ_ＤＰであり、Ｒハット_ａ＋Ｒハット_ｂまたはＲハット_ａ＋Ｒハット_ｃと、Ｒハット_ｂまたはＲハット_ｃとの位相差は、φ_ＤＰ／２である。
また、Ｒハット_ａ＋Ｒハット_ｂ，Ｒハット_ａ＋Ｒハット_ｃの位相は、ψ_ｄである。

これにより、送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰは、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｂ（Ｔ_ｓ）との偏角の２／３倍、または、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角の２／３倍で表される。
即ち、送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰは、下記の式（２３）のように表される。

ここでは、速度算出部１３が、送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰとして、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｂ（Ｔ_ｓ）との偏角の２／３倍、または、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角の２／３倍を算出する例を示したが、送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰとして、下記の式（２４）に示すように、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｂ（Ｔ_ｓ）との偏角の２／３倍と、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角の２／３倍との平均を算出するようにしてもよい。

速度算出部１３は、偏波間の位相差φハット_ＤＰを算出すると、下記の式（２５）に示すように、ラグＲハット（Ｔ_ｓ）に含まれているＲハット_ａ（Ｔ_ｓ）、偏波間の位相差φハット_ＤＰ及びドップラー速度ｖハットを観測することが可能な範囲であるナイキスト速度ｖ_ａから、観測対象のドップラー速度ｖハットを算出する。

ここで、式（２５）は、単偏波の送受信処理や、複数の偏波の同時送信同時受信の送受信処理を行う場合のドップラー速度ｖハットの算出式と同様であり、ドップラー速度ｖハットを観測することが可能な範囲であるナイキスト速度ｖ_ａが、単偏波の送受信処理や、複数の偏波の同時送信同時受信処理を行う場合と同等のナイキスト速度まで回復している。
したがって、３種類の偏波送受信処理を繰り返す場合でも、単偏波の送受信処理や、複数の偏波の同時送信同時受信処理を行う場合と同等のナイキスト速度で、観測対象のドップラー速度ｖハットを算出することができる。

以上で明らかなように、この実施の形態１によれば、交差偏波比算出部１２が、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈ（ｎ）の積分後の反射波強度である水平偏波送信水平偏波受信でのレーダ反射因子Ｚ_ｈｈと、反射波強度Ｖ_ｖｈ（ｎ＋２）の積分後の反射波強度である水平偏波送信垂直偏波受信でのレーダ反射因子Ｚ_ｖｈとを求め、そのレーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈとの比である交差偏波比ＬＤＲ_ｖｈを算出するように構成したので、３種類の偏波送受信処理を繰り返す場合でも、ナイキスト速度ｖ_ａが減少することなく、交差偏波比ＬＤＲ_ｖｈを算出することができる効果を奏する。

この実施の形態１では、送受信装置１が、第１の送受信処理である「水平偏波送信水平偏波受信」と、第２の送受信処理である「垂直偏波送信垂直偏波受信」と、第３の送受信処理である「水平偏波送信垂直偏波受信」とを実施する例を示したが、第１の送受信処理である「水平偏波送信水平偏波受信」と、第２の送受信処理である「垂直偏波送信垂直偏波受信」と、第４の送受信処理である「垂直偏波送信水平偏波受信」とを実施するものであってもよい。
この場合、反射波強度算出部１１は、複数回の第３の送受信処理で受信された垂直偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、（ｎ＋２）番目のサンプルにおける反射波強度Ｖ_ｖｈ（ｎ＋２）を算出する代わりに、複数回の第４の送受信処理で受信された水平偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、下記の式（２６）に示すように、（ｎ＋２）番目のサンプルにおける反射波強度Ｖ_ｈｖ（ｎ＋２）を算出する。

また、反射波強度算出部１１は、反射波強度Ｖ_ｈｖ（ｎ＋２）を積分することで、反射波強度Ｖ_ｈｖ（ｎ＋２）の積分後の反射波強度である垂直偏波送信水平偏波受信でのレーダ反射因子Ｚ_ｈｖを求める。
この積分範囲は、レーダ装置の精度等を考慮して、ユーザが任意に決定することができるが、例えば、送受信装置１により送受信される各偏波の間隔Ｔ_ｓと、送受信装置１による既知のビーム範囲（レーダを回転させる範囲）とから決定されることが想定される。
さらに、反射波強度算出部１１は、レーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖを交差偏波比算出部１２に出力するとともに、レーダ反射因子Ｚ_ｈｈ，Ｚ_ｖｖ，Ｚ_ｈｖ及びレーダ反射因子差Ｚｄｒを表示装置７に表示する。

交差偏波比算出部１２は、反射波強度算出部１１からレーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖを受けると、そのレーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖを下記の式（２７）に代入することで、交差偏波比ＬＤＲ_ｈｖを算出する。

交差偏波比算出部１２は、交差偏波比ＬＤＲ_ｖｈを算出すると、その交差偏波比ＬＤＲ_ｖｈを表示装置７に表示する。

速度算出部１３は、反射波強度算出部１１が反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｈｖを算出すると、その反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｈｖから、パルスペア法を用いて、観測対象のドップラー速度ｖハットを算出する。
即ち、速度算出部１３は、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈ，Ｖ_ｖｖ，Ｖ_ｈｖから、送受信装置１により偏波が繰り返し送信される周期Ｔ_ｓでの各送信偏波における自己相関のラグＲハット（Ｔ_ｓ）を算出して、そのラグＲハット（Ｔ_ｓ）から送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰを算出し、そのラグＲハット（Ｔ_ｓ）、偏波間の位相差φハット_ＤＰ及びナイキスト速度ｖ_ａから、観測対象のドップラー速度ｖハットを算出する。
反射波強度算出部１１により算出された反射波強度を用いて、ドップラー速度ｖハットを算出する場合、偏波の交互送信における自己相関の推定法（パルスペア法）では、ラグＲハット（Ｔ_ｓ）は、下記の式（２８）のように表される。

式（２８）に示すＲハット（Ｔ_ｓ）は、上記の式（１１）（１２）と、下記の式（２９）とに示すように、３つの式に分解することができる。即ち、式（２８）に示すＲハット（Ｔ_ｓ）を構成しているＲハット_ａ（Ｔ_ｓ）、Ｒハット_ｂ（Ｔ_ｓ）、Ｒハット_ｃ（Ｔ_ｓ）を特定することができる。

Ｒハット_ｃ（Ｔ_ｓ）は、水平偏波受信時における垂直偏波観点の共分散、即ち、第４の送受信処理で水平偏波が受信された時点での当該水平偏波と第１の送受信処理で受信された水平偏波との共分散（第４の共分散）である。

式（２９）に示すＲハット_ｃ（Ｔ_ｓ）は、式（４）に示す反射波強度Ｖ_ｈｈ（ｎ）と、式（２６）に示す反射波強度Ｖ_ｈｖ（ｎ＋２）を用いると、下記の式（３０）のように表される。

式（３０）に示すＲハット_ｃ（Ｔ_ｓ）は、
φ_ＤＰ＝２（ｋ_ｈ−ｋ_ｖ）ｒ_０，ψ_ｄ＝−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖより、下記の式（３１）のように表される。

以降の速度算出部１３の処理内容は、第３の送受信処理である「水平偏波送信垂直偏波受信」を実施する場合と同様である。

実施の形態２．
上記実施の形態１では、送受信装置１が、第３の送受信処理である「水平偏波送信垂直偏波受信」又は第４の送受信処理である「垂直偏波送信水平偏波受信」のいずれかを実施するものを示したが、第３の送受信処理である「水平偏波送信垂直偏波受信」と、第４の送受信処理である「垂直偏波送信水平偏波受信」との双方を実施することで、４種類の偏波送受信処理を繰り返すようにしてもよい。

図８は送受信装置１により送受信される各偏波の間隔がＴ_ｓである例を示す説明図である。
図８の例では、第１の送受信処理である「水平偏波送信水平偏波受信」が（４ｉ−３）・Ｔ_ｓの時間に行われ、第３の送受信処理である「水平偏波送信垂直偏波受信」が（４ｉ−２）・Ｔ_ｓの時間に行われ、第２の送受信処理である「垂直偏波送信垂直偏波受信」が（４ｉ−１）・Ｔ_ｓの時間に行われ、第４の送受信処理である「垂直偏波送信水平偏波受信」が４ｉ・Ｔ_ｓの時間に行われている。
ただし、ｉ＝１，２，３，・・・，Ｉである。
Ｉはサンプルペア数（水平偏波送信水平偏波受信、垂直偏波送信垂直偏波受信、水平偏波送信垂直偏波受信、垂直偏波送信水平偏波受信それぞれのヒット数）である。
即ち、図８の例では、実施順序は、「水平偏波送信水平偏波受信」→「水平偏波送信垂直偏波受信」→「垂直偏波送信垂直偏波受信」→「垂直偏波送信水平偏波受信」の順番である。
この実施の形態２による二重偏波レーダ装置の構成図は、上記実施の形態１による二重偏波レーダ装置と同様に図１である。
次に動作について説明する。

送受信装置１の送信機２は、「水平偏波送信水平偏波受信」→「水平偏波送信垂直偏波受信」→「垂直偏波送信垂直偏波受信」→「垂直偏波送信水平偏波受信」の順番で送受信処理を実施するため、送信周期Ｔ_ｓで、水平偏波→水平偏波→垂直偏波→垂直偏波→水平偏波→水平偏波→垂直偏波→垂直偏波→・・・→水平偏波→水平偏波→垂直偏波→垂直偏波を送信する。
送受切替器３は、送信機２から水平偏波又は垂直偏波を受けると、その水平偏波又は垂直偏波をアンテナ４に出力する。
アンテナ４は、送受切替器３から水平偏波又は垂直偏波を受けると、その水平偏波又は垂直偏波を空間に向けて放射する。
アンテナ４は、水平偏波又は垂直偏波を空間に向けて放射したのち、観測対象である雨滴に反射されて戻ってきた水平偏波又は垂直偏波を受信する。
図８の例では、送信周期Ｔ_ｓで、水平偏波→垂直偏波→垂直偏波→水平偏波→水平偏波→垂直偏波→垂直偏波→水平偏波→・・・→水平偏波→垂直偏波→垂直偏波→水平偏波を受信する。

送受切替器３は、アンテナ４により受信された水平偏波又は垂直偏波を受信機５に出力する。
受信機５は、アンテナ４により受信された水平偏波又は垂直偏波の受信処理を実施し、その水平偏波又は垂直偏波の受信信号をアナログ信号からデジタル信号に変換して、デジタルの受信信号を信号処理装置６に出力する。

信号処理装置６の反射波強度算出部１１は、送受信装置１により第１の送受信処理で受信された水平偏波のデジタル受信信号から、当該水平偏波の電力値である反射波強度Ｖ_ｈｈを算出するとともに、第３の送受信処理で受信された垂直偏波のデジタル受信信号から、当該垂直偏波の電力値である反射波強度Ｖ_ｖｈを算出する。
また、第２の送受信処理で受信された垂直偏波のデジタル受信信号から、当該垂直偏波の電力値である反射波強度Ｖ_ｖvを算出するとともに、第４の送受信処理で受信された水平偏波のデジタル受信信号から、当該水平偏波の電力値である反射波強度Ｖ_ｈｖを算出する。
以下、反射波強度算出部１１による反射波強度Ｖ_ｈｈ，Ｖ_ｖｈ，Ｖ_ｖｖ，Ｖ_ｈｖの算出処理を具体的に説明する。

反射波強度算出部１１は、複数回の第１の送受信処理で受信された水平偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、ｎ番目のサンプルにおける反射波強度Ｖ_ｈｈ（ｎ）を下記の式（３２）のように算出する。
また、反射波強度算出部１１は、複数回の第３の送受信処理で受信された垂直偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、（ｎ＋１）番目のサンプルにおける反射波強度Ｖ_ｖｈ（ｎ＋１）を下記の式（３３）のように算出する。
反射波強度算出部１１は、複数回の第２の送受信処理で受信された垂直偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、（ｎ＋２）番目のサンプルにおける反射波強度Ｖ_ｖｖ（ｎ＋２）を下記の式（３４）のように算出する。
さらに、反射波強度算出部１１は、複数回の第４の送受信処理で受信された水平偏波のデジタル受信信号である後方散乱受信信号の電圧の合計、即ち、（ｎ＋３）番目のサンプルにおける反射波強度Ｖ_ｈｖ（ｎ＋３）を下記の式（３５）のように算出する。なお、これらの反射波強度は、観測領域の方位角と距離の区分（格子）毎に算出される。

反射波強度算出部１１は、反射波強度Ｖ_ｈｈ（ｎ），Ｖ_ｖｈ（ｎ＋１），Ｖ_ｖｖ（ｎ＋２），Ｖ_ｈｖ（ｎ＋３）を算出すると、その反射波強度Ｖ_ｈｈ（ｎ），Ｖ_ｖｈ（ｎ＋１），Ｖ_ｖｖ（ｎ＋２），Ｖ_ｈｖ（ｎ＋３）のそれぞれを積分することで、反射波強度Ｖ_ｈｈ（ｎ）の積分後の反射波強度である水平偏波送信水平偏波受信でのレーダ反射因子Ｚ_ｈｈと、反射波強度Ｖ_ｖｈ（ｎ＋１）の積分後の反射波強度である水平偏波送信垂直偏波受信でのレーダ反射因子Ｚ_ｖｈと、反射波強度Ｖ_ｖｖ（ｎ＋２）の積分後の反射波強度である垂直偏波送信垂直偏波受信でのレーダ反射因子Ｚ_ｖｖと、反射波強度Ｖ_ｈｖ（ｎ＋３））の積分後の反射波強度である垂直偏波送信水平偏波受信でのレーダ反射因子Ｚ_ｈｖとを求める。
この積分範囲は、レーダ装置の精度等を考慮して、ユーザが任意に決定することができるが、例えば、送受信装置１により送受信される各偏波の間隔Ｔ_ｓと、送受信装置１による既知のビーム範囲（レーダを回転させる範囲）とから決定されることが想定される。

また、反射波強度算出部１１は、上記の式（８）に示すように、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｖとの比であるレーダ反射因子差Ｚｄｒを算出する。
さらに、反射波強度算出部１１は、レーダ反射因子Ｚ_ｈｈ，Ｚ_ｖｈ，Ｚ_ｖｖ，Ｚ_ｈｖを交差偏波比算出部１２に出力するとともに、レーダ反射因子Ｚ_ｈｈ，Ｚ_ｖｈ，Ｚ_ｖｖ，Ｚ_ｈｖ及びレーダ反射因子差Ｚｄｒを表示装置７に表示する。

交差偏波比算出部１２は、反射波強度算出部１１からレーダ反射因子Ｚ_ｈｈ，Ｚ_ｖｈ，Ｚ_ｖｖ，Ｚ_ｈｖを受けると、そのレーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈを下記の式（３６）に代入することで、交差偏波比ＬＤＲ_ｖｈを算出する。

また、交差偏波比算出部１２は、そのレーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖを下記の式（３７）に代入することで、交差偏波比ＬＤＲ_ｈｖを算出する。

交差偏波比算出部１２は、交差偏波比ＬＤＲ_ｖｈ，ＬＤＲ_ｈｖを算出すると、その交差偏波比ＬＤＲ_ｖｈ，ＬＤＲ_ｈｖを表示装置７に表示する。

速度算出部１３は、反射波強度算出部１１が反射波強度Ｖ_ｈｈ（ｎ），Ｖ_ｖｈ（ｎ＋１），Ｖ_ｖｖ（ｎ＋２），Ｖ_ｈｖ（ｎ＋３）を算出すると、その反射波強度Ｖ_ｈｈ（ｎ），Ｖ_ｖｈ（ｎ＋１），Ｖ_ｖｖ（ｎ＋２），Ｖ_ｈｖ（ｎ＋３）から、パルスペア法を用いて、観測対象のドップラー速度ｖハットを算出する。
即ち、速度算出部１３は、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈ（ｎ），Ｖ_ｖｈ（ｎ＋１），Ｖ_ｖｖ（ｎ＋２），Ｖ_ｈｖ（ｎ＋３）から、送受信装置１により偏波が繰り返し送信される周期Ｔ_ｓでの各送信偏波における自己相関のラグＲハット（Ｔ_ｓ）を算出して、そのラグＲハット（Ｔ_ｓ）から送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰを算出し、そのラグＲハット（Ｔ_ｓ）、偏波間の位相差φハット_ＤＰ及びナイキスト速度ｖ_ａから、観測対象のドップラー速度ｖハットを算出する。
反射波強度算出部１１により算出された反射波強度を用いて、ドップラー速度ｖハットを算出する場合、偏波の交互送信における自己相関の推定法（パルスペア法）では、ラグＲハット（Ｔ_ｓ）は、下記の式（３８）のように表される。

式（３８）に示すＲハット（Ｔ_ｓ）は、下記の式（３９）〜（４２）に示すように、４つの式に分解することができる。即ち、式（３８）に示すＲハット（Ｔ_ｓ）を構成しているＲハット_ａ（Ｔ_ｓ）、Ｒハット_ｂ（Ｔ_ｓ）、Ｒハット_ｃ（Ｔ_ｓ）、Ｒハット_ｄ（Ｔ_ｓ）を特定することができる。

Ｒハット_ａ（Ｔ_ｓ）は、第１の送受信処理で水平偏波が受信された時点での当該水平偏波と第３の送受信処理で受信された垂直偏波との共分散（第１の共分散）である。
Ｒハット_ｂ（Ｔ_ｓ）は、第３の送受信処理で垂直偏波が受信された時点での当該垂直偏波と第２の送受信処理で受信された垂直偏波との共分散（第２の共分散）である。
Ｒハット_ｃ（Ｔ_ｓ）は、第２の送受信処理で垂直偏波が受信された時点での当該垂直偏波と第４の送受信処理で受信された水平偏波との共分散（第３の共分散）である。
Ｒハット_ｄ（Ｔ_ｓ）は、第４の送受信処理で水平偏波が受信された時点での当該水平偏波と第１の送受信処理で受信された水平偏波との共分散（第４の共分散）である。

これにより、式（３８）に示すＲハット（Ｔ_ｓ）は、下記の式（４３）のように表される。

式（３９）に示すＲハット_ａ（Ｔ_ｓ）は、式（３２）に示す反射波強度Ｖ_ｈｈ（ｎ）と、式（３３）に示す反射波強度Ｖ_ｖｈ（ｎ＋１）とを用いると、下記の式（４４）のように表される。

式（４４）に示すＲハット_ａ（Ｔ_ｓ）は、
φ_ＤＰ＝２（ｋ_ｈ−ｋ_ｖ）ｒ_０，ψ_ｄ＝−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖ＝−２（ｋ_０＋ｋ_ｖ）Ｔ_ｓｖより、下記の式（４５）のように表される。

式（４０）に示すＲハット_ｂ（Ｔ_ｓ）は、式（３３）に示す反射波強度Ｖ_ｖｈ（ｎ＋１）と、式（３４）に示す反射波強度Ｖ_ｖｖ（ｎ＋２）を用いると、下記の式（４６）のように表される。

式（４６）に示すＲハット_ｂ（Ｔ_ｓ）は、
φ_ＤＰ＝２（ｋ_ｈ−ｋ_ｖ）ｒ_０，ψ_ｄ＝−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖより、下記の式（４７）のように表される。

式（４１）に示すＲハット_ｃ（Ｔ_ｓ）は、式（３４）に示す反射波強度Ｖ_ｖｖ（ｎ＋２）と、式（３５）に示す反射波強度Ｖ_ｈｖ（ｎ＋３）とを用いると、下記の式（４８）のように表される。

式（４８）に示すＲハット_ｃ（Ｔ_ｓ）は、
φ_ＤＰ＝２（ｋ_ｈ−ｋ_ｖ）ｒ_０，ψ_ｄ＝−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖ＝−２（ｋ_０＋ｋ_ｖ）Ｔ_ｓｖより、下記の式（４９）のように表される。

式（４２）に示すＲハット_ｄ（Ｔ_ｓ）は、式（３５）に示す反射波強度Ｖ_ｈｖ（ｎ＋３）と、式（３２）に示す反射波強度Ｖ_ｈｈ（ｎ）を用いると、下記の式（５０）のように表される。

式（５０）に示すＲハット_ｄ（Ｔ_ｓ）は、
φ_ＤＰ＝２（ｋ_ｈ−ｋ_ｖ）ｒ_０，ψ_ｄ＝−２（ｋ_０＋ｋ_ｈ）Ｔ_ｓｖ＝−２（ｋ_０＋ｋ_ｖ）Ｔ_ｓｖより、下記の式（５１）のように表される。

ここで、図９は式（４５）、式（４７）、及び式（４９）及び式（５１）から分かるＲハット_ａ（Ｔ_ｓ），Ｒハット_ｂ（Ｔ_ｓ），Ｒハット_ｃ（Ｔ_ｓ），Ｒハット_ｄ（Ｔ_ｓ）の位相の関係を示す説明図である。図９では、図面の簡単化のため、（Ｔ_ｓ）の記述を省略して、Ｒハット_ａ，Ｒハット_ｂ，Ｒハット_ｃ，Ｒハット_ｄのように表記している。
図９に示すように、Ｒハット_ａまたはＲハット_ｂと、Ｒハット_ａ（またはＲハット_ｂ）＋Ｒハット_ｃ（またはＲハット_ｄ）との位相差は、φ_ＤＰ／２であり、Ｒハット_ｃまたはＲハット_ｄと、Ｒハット_ａ（またはＲハット_ｂ）＋Ｒハット_ｃ（またはＲハット_ｄ）との位相差は、φ_ＤＰ／２である。
また、Ｒハット_ａ（またはＲハット_ｂ）＋Ｒハット_ｃ（またはＲハット_ｄ）の位相は、ψ_ｄである。

これにより、送受信装置１により繰り返し送信される偏波間の位相差φ_ＤＰは、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｄ（Ｔ_ｓ）との偏角、Ｒハット_ｂ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角、または、Ｒハット_ｂ（Ｔ_ｓ）とＲハット_ｄ（Ｔ_ｓ）との偏角で表される。
即ち、送受信装置１により繰り返し送信される偏波間の位相差φハット_ＤＰは、下記の式（５２）のように表される。

ここでは、速度算出部１３が、送受信装置１により繰り返し送信される偏波間の位相差φ_ＤＰとして、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｄ（Ｔ_ｓ）との偏角、Ｒハット_ｂ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角、または、Ｒハット_ｂ（Ｔ_ｓ）とＲハット_ｄ（Ｔ_ｓ）との偏角を算出する例を示したが、送受信装置１により繰り返し送信される偏波間の位相差φ_ＤＰとして、下記の式（５３）に示すように、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角と、Ｒハット_ａ（Ｔ_ｓ）とＲハット_ｄ（Ｔ_ｓ）との偏角と、Ｒハット_ｂ（Ｔ_ｓ）とＲハット_ｃ（Ｔ_ｓ）との偏角と、Ｒハット_ｂ（Ｔ_ｓ）とＲハット_ｄ（Ｔ_ｓ）との偏角との平均を算出するようにしてもよい。

速度算出部１３は、偏波間の位相差φハット_ＤＰを算出すると、下記の式（５４）に示すように、ラグＲハット（Ｔ_ｓ）に含まれているＲハット_ａ（Ｔ_ｓ）、偏波間の位相差φハット_ＤＰ及びドップラー速度ｖハットを観測することが可能な範囲であるナイキスト速度ｖ_ａから、観測対象のドップラー速度ｖハットを算出する。

ここで、式（５４）は、単偏波の送受信処理や、複数の偏波の同時送信同時受信処理を行う場合のドップラー速度ｖハットの算出式と同様であり、ドップラー速度ｖハットを観測することが可能な範囲であるナイキスト速度ｖ_ａが、単偏波の送受信処理や、複数の偏波の同時送信同時受信処理を行う場合と同等のナイキスト速度まで回復している。
したがって、４種類の偏波送受信処理を繰り返す場合でも、単偏波の送受信処理や、複数の偏波の同時送信同時受信処理を行う場合と同等のナイキスト速度ｖ_ａで、観測対象のドップラー速度ｖハットを算出することができる。
なお、Ｒハット_ａ（Ｔ_ｓ）の代わりに、Ｒハット_ｂ（Ｔ_ｓ）、Ｒハット_ｃ（Ｔ_ｓ）又はＲハット_ｄ（Ｔ_ｓ）を式（５４）に代入することで、ドップラー速度ｖハットを算出するようにしてもよい。

以上で明らかなように、この実施の形態２によれば、交差偏波比算出部１２が、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈとの比である交差偏波比ＬＤＲ_ｖｈ、または、レーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖとの比である交差偏波比ＬＤＲ_ｈｖを算出するように構成したので、４種類の偏波送受信処理を繰り返す場合でも、ナイキスト速度ｖ_ａが減少することなく、交差偏波比ＬＤＲ_ｖｈまたは交差偏波比ＬＤＲ_ｈｖを算出することができる効果を奏する。

実施の形態３．
図１０はこの発明の実施の形態３による二重偏波レーダ装置の信号処理装置６を示す構成図であり、図１０において、図２と同一符号は同一または相当部分を示すので説明を省略する。
また、図１１は図１０の信号処理装置６のハードウェア構成図であり、図１１において、図３と同一符号は同一または相当部分を示すので説明を省略する。
有効無効判定部１４は例えばＣＰＵを実装している半導体処理回路あるいはワンチップマイコンなどから構成されている有効無効判定処理回路２６から実現されるものであり、反射波強度算出部１１により算出されたレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が事前に設定された閾値Ｔｈ以上であれば、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が有効であると判定し、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が閾値Ｔｈ未満であれば、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が無効であると判定する処理を実施する。
この実施の形態３では、閾値Ｔｈとしては、所定のノイズレベルから規定レベルだけ上の値が設定されるものを想定しているが、ユーザが任意に設定することができる。また、規定レベルは、レーダ装置の性能などから任意に決定される。

交差偏波比算出部１５は例えば交差偏波比算出処理回路２３と、表示処理回路２５とから実現されるものであり、有効無効判定部１４により有効であると判定されれば、図２の交差偏波比算出部１２と同様に、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈと反射波強度Ｖ_ｖｈとの比、即ち、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈとの比である交差偏波比ＬＤＲ_ｖｈを算出する。あるいは、反射波強度算出部１１により算出された反射波強度Ｖ_ｖｖと反射波強度Ｖ_ｈｖとの比、即ち、レーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖとの比である交差偏波比ＬＤＲ_ｈｖを算出する処理を実施する。また、算出した交差偏波比ＬＤＲ_ｖｈ、あるいは、交差偏波比ＬＤＲ_ｈｖを表示装置７に表示する処理を実施する。
また、交差偏波比算出部１５は有効無効判定部１４により無効であると判定されれば、交差偏波比ＬＤＲ_ｖｈ，ＬＤＲ_ｖｈの算出処理を実行しない。

速度算出部１６は速度算出処理回路２４と、表示処理回路２５とから実現されるものであり、有効無効判定部１４により有効であると判定されれば、図２の速度算出部１３と同様に、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈと、反射波強度Ｖ_ｖｖと、反射波強度Ｖ_ｖｈ又は反射波強度Ｖ_ｈｖとから、パルスペア法を用いて、観測対象のドップラー速度ｖハットを算出する処理を実施する。上記実施の形態２と同様に、反射波強度Ｖ_ｈｈと、反射波強度Ｖ_ｖｖと、反射波強度Ｖ_ｖｈと、反射波強度Ｖ_ｈｖとから、観測対象のドップラー速度ｖハットを算出するようにしてもよい。また、速度算出部１６は算出したドップラー速度ｖハットを表示装置７に表示する処理を実施する。
また、速度算出部１６は有効無効判定部１４により無効であると判定されれば、例えば、非特許文献１に記載されているレーダ装置と同様に、反射波強度Ｖ_ｈｈと、反射波強度Ｖ_ｖｖとから、観測対象のドップラー速度ｖハットを算出する処理を実施する。

図１０の例では、信号処理装置６の構成要素である反射波強度算出部１１、有効無効判定部１４、交差偏波比算出部１５及び速度算出部１６のそれぞれが専用のハードウェアで構成されているものを想定しているが、信号処理装置６がコンピュータで構成されているものであってもよい。
信号処理装置６がコンピュータで構成される場合、反射波強度算出部１１、有効無効判定部１４、交差偏波比算出部１５及び速度算出部１６の処理内容を記述しているプログラムを図４に示すコンピュータのメモリ３１に格納し、そのコンピュータのプロセッサ３２がメモリ３１に格納されているプログラムを実行するようにすればよい。

次に動作について説明する。
上記実施の形態１では、交差偏波比算出部１２が、常に、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈとの比である交差偏波比ＬＤＲ_ｖｈ、あるいは、レーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖとの比である交差偏波比ＬＤＲ_ｈｖを算出するものを示したが、交差偏波の電力が低く、信号として不確かである場合には、交差偏波比ＬＤＲ_ｖｈや交差偏波比ＬＤＲ_ｈｖを高精度に算出することができない。
そこで、この実施の形態３では、交差偏波比ＬＤＲ_ｖｈや交差偏波比ＬＤＲ_ｈｖを高精度に算出することができる状況下にある場合だけ、交差偏波比ＬＤＲ_ｖｈまたは交差偏波比ＬＤＲ_ｈｖを算出するようにする。
具体的には、以下の通りである。

有効無効判定部１４は、反射波強度算出部１１がレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）を算出すると、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）と事前に設定された閾値Ｔｈを比較する。
有効無効判定部１４は、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が閾値Ｔｈ以上であれば、交差偏波の電力が高く、信号として確かであるため、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が有効であると判定する。
一方、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が閾値Ｔｈ未満であれば、交差偏波の電力が低く、信号として不確かであるため、そのレーダ反射因子Ｚ_ｖｈ（またはレーダ反射因子Ｚ_ｈｖ）が無効であると判定する。

交差偏波比算出部１５は、有効無効判定部１４により有効であると判定されれば、図２の交差偏波比算出部１２と同様に、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈと反射波強度Ｖ_ｖｈとの比、即ち、レーダ反射因子Ｚ_ｈｈとレーダ反射因子Ｚ_ｖｈとの比である交差偏波比ＬＤＲ_ｖｈを算出する。あるいは、反射波強度算出部１１により算出された反射波強度Ｖ_ｖｖと反射波強度Ｖ_ｈｖとの比、即ち、レーダ反射因子Ｚ_ｖｖとレーダ反射因子Ｚ_ｈｖとの比である交差偏波比ＬＤＲ_ｈｖを算出する。
また、交差偏波比算出部１５は、図２の交差偏波比算出部１２と同様に、その算出した交差偏波比ＬＤＲ_ｖｈ、あるいは、交差偏波比ＬＤＲ_ｈｖを表示装置７に表示する。
交差偏波比算出部１５は、有効無効判定部１４により無効であると判定されれば、交差偏波比ＬＤＲ_ｖｈ，ＬＤＲ_ｖｈの算出処理を実行しない。

速度算出部１６は、有効無効判定部１４により有効であると判定されれば、図２の速度算出部１３と同様に、反射波強度算出部１１により算出された反射波強度Ｖ_ｈｈと、反射波強度Ｖ_ｖｖと、反射波強度Ｖ_ｖｈ又は反射波強度Ｖ_ｈｖとから、パルスペア法を用いて、観測対象のドップラー速度ｖハットを算出し、そのドップラー速度ｖハットを表示装置７に表示する。
速度算出部１６は、有効無効判定部１４により無効であると判定されれば、例えば、非特許文献１に記載されているレーダ装置と同様に、２種類の偏波送受信処理で得られた反射波強度Ｖ_ｈｈと反射波強度Ｖ_ｖｖとから、観測対象のドップラー速度ｖハットを算出し、そのドップラー速度ｖハットを表示装置７に表示する。

以上で明らかなように、この実施の形態３によれば、有効無効判定部１４により有効であると判定された場合に、交差偏波比算出部１５が、交差偏波比ＬＤＲ_ｖｈ、あるいは、交差偏波比ＬＤＲ_ｈｖを算出するように構成したので、高精度に算出された交差偏波比だけを表示装置７に表示することができる効果を奏する。

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

この発明は、広いナイキスト速度でドップラー速度を推定する必要がある二重偏波レーダ装置及びレーダ信号処理方法に適している。

１ 送受信装置、２ 送信機、３ 送受切替器、４ アンテナ、５ 受信機、６ 信号処理装置、７ 表示装置、１１ 反射波強度算出部、１２，１５ 交差偏波比算出部、１３，１６ 速度算出部、１４ 有効無効判定部、２１ 信号入力処理回路、２２ 反射波強度算出処理回路、２３ 交差偏波比算出処理回路、２４ 速度算出処理回路、２５ 表示処理回路、２６ 有効無効判定処理回路、３１ メモリ、３２ プロセッサ。

Claims (8)

1. 水平偏波を送信して水平偏波を受信する第１の送受信処理と、垂直偏波を送信して垂直偏波を受信する第２の送受信処理と、水平偏波を送信して垂直偏波を受信する第３の送受信処理又は垂直偏波を送信して水平偏波を受信する第４の送受信処理とを順番に繰り返し実施する送受信装置と、
   前記送受信装置により前記第１の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第１の反射波強度を算出し、前記第２の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第２の反射波強度を算出し、前記第３の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第３の反射波強度又は前記第４の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第４の反射波強度を算出する反射波強度算出部と、
   前記反射波強度算出部により算出された前記第１の反射波強度と前記第３の反射波強度との比、あるいは、前記第２の反射波強度と前記第４の反射波強度との比である交差偏波比を算出する交差偏波比算出部と、
   前記反射波強度算出部により算出された前記第１の反射波強度と、前記第２の反射波強度と、前記第３の反射波強度又は前記第４の反射波強度とから、前記送受信装置により偏波が繰り返し送信される周期での各送信偏波における自己相関のラグを算出して、前記ラグから前記送受信装置により繰り返し送信される偏波間の位相差を算出し、前記ラグ、前記偏波間の位相差及びドップラー速度を観測することが可能な範囲であるナイキスト速度から、偏波を反射している観測対象のドップラー速度を算出する速度算出部と、
   を備えた二重偏波レーダ装置。
2. 前記速度算出部は、
   前記ラグから、前記送受信装置により前記第１の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第２の送受信処理で受信された垂直偏波との共分散である第１の共分散と、前記第２の送受信処理で垂直偏波が受信された時点での当該垂直偏波と前記第３の送受信処理で受信された垂直偏波又は前記第４の送受信処理で受信された水平偏波との共分散である第２の共分散と、前記第３の送受信処理で垂直偏波が受信された時点での当該垂直偏波又は前記第４の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第１の送受信処理で受信された水平偏波との共分散である第３の共分散とを特定し、
   前記送受信装置により繰り返し送信される偏波間の位相差として、前記第１の共分散と前記第２の共分散との偏角の３分の２倍、または、前記第１の共分散と前記第３の共分散との偏角の３分の２倍を算出することを特徴とする請求項１記載の二重偏波レーダ装置。
3. 前記速度算出部は、
   前記ラグから、前記送受信装置により前記第１の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第２の送受信処理で受信された垂直偏波との共分散である第１の共分散と、前記第２の送受信処理で垂直偏波が受信された時点での当該垂直偏波と前記第３の送受信処理で受信された垂直偏波又は前記第４の送受信処理で受信された水平偏波との共分散である第２の共分散と、前記第３の送受信処理で垂直偏波が受信された時点での当該垂直偏波又は前記第４の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第１の送受信処理で受信された水平偏波との共分散である第３の共分散とを特定し、
   前記送受信装置により繰り返し送信される偏波間の位相差として、前記第１の共分散と前記第２の共分散との偏角の３分の２倍と、前記第１の共分散と前記第３の共分散との偏角の３分の２倍との平均を算出することを特徴とする請求項１記載の二重偏波レーダ装置。
4. 水平偏波を送信して水平偏波を受信する第１の送受信処理と、水平偏波を送信して垂直偏波を受信する第３の送受信処理と、垂直偏波を送信して垂直偏波を受信する第２の送受信処理と、垂直偏波を送信して水平偏波を受信する第４の送受信処理とを順番に繰り返し実施する送受信装置と、
   前記送受信装置により前記第１の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第１の反射波強度を算出し、前記第２の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第２の反射波強度を算出し、前記第３の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第３の反射波強度を算出し、前記第４の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第４の反射波強度を算出する反射波強度算出部と、
   前記反射波強度算出部により算出された前記第１の反射波強度と前記第３の反射波強度との比、あるいは、前記第２の反射波強度と前記第４の反射波強度との比である交差偏波比を算出する交差偏波比算出部と、
   前記反射波強度算出部により算出された前記第１の反射波強度と、前記第２の反射波強度と、前記第３の反射波強度と、前記第４の反射波強度とから、前記送受信装置により偏波が繰り返し送信される周期での各送信偏波における自己相関のラグを算出して、前記ラグから前記送受信装置により繰り返し送信される偏波間の位相差を算出し、前記ラグ、前記偏波間の位相差及びドップラー速度を観測することが可能な範囲であるナイキスト速度から、偏波を反射している観測対象のドップラー速度を算出する速度算出部と、
   を備えた二重偏波レーダ装置。
5. 前記速度算出部は、
   前記ラグから、前記第１の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第３の送受信処理で受信された垂直偏波との共分散である第１の共分散と、前記第３の送受信処理で垂直偏波が受信された時点での当該垂直偏波と前記第２の送受信処理で受信された垂直偏波との共分散である第２の共分散と、前記第２の送受信処理で垂直偏波が受信された時点での当該垂直偏波と前記第４の送受信処理で受信された水平偏波との共分散である第３の共分散と、前記第４の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第１の送受信処理で受信された水平偏波との共分散である第４の共分散とを特定し、
   前記送受信装置により繰り返し送信される偏波間の位相差として、前記第１の共分散と前記第３の共分散との偏角、前記第１の共分散と前記第４の共分散との偏角、前記第２の共分散と前記第３の共分散との偏角、または、前記第２の共分散と前記第４の共分散との偏角を算出することを特徴とする請求項４記載の二重偏波レーダ装置。
6. 前記速度算出部は、
   前記ラグから、前記第１の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第３の送受信処理で受信された垂直偏波との共分散である第１の共分散と、前記第３の送受信処理で垂直偏波が受信された時点での当該垂直偏波と前記第２の送受信処理で受信された垂直偏波との共分散である第２の共分散と、前記第２の送受信処理で垂直偏波が受信された時点での当該垂直偏波と前記第４の送受信処理で受信された水平偏波との共分散である第３の共分散と、前記第４の送受信処理で水平偏波が受信された時点での当該水平偏波と前記第１の送受信処理で受信された水平偏波との共分散である第４の共分散とを特定し、
   前記送受信装置により繰り返し送信される偏波間の位相差として、前記第１の共分散と前記第３の共分散との偏角と、前記第１の共分散と前記第４の共分散との偏角と、前記第２の共分散と前記第３の共分散との偏角と、前記第２の共分散と前記第４の共分散との偏角との平均を算出することを特徴とする請求項４記載の二重偏波レーダ装置。
   
7. 送受信装置が、水平偏波を送信して水平偏波を受信する第１の送受信処理と、垂直偏波を送信して垂直偏波を受信する第２の送受信処理と、水平偏波を送信して垂直偏波を受信する第３の送受信処理又は垂直偏波を送信して水平偏波を受信する第４の送受信処理とを順番に繰り返し実施し、
   反射波強度算出部が、前記送受信装置により前記第１の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第１の反射波強度を算出し、前記第２の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第２の反射波強度を算出し、前記第３の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第３の反射波強度又は前記第４の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第４の反射波強度を算出し、
   交差偏波比算出部が、前記反射波強度算出部により算出された前記第１の反射波強度と前記第３の反射波強度との比、あるいは、前記第２の反射波強度と前記第４の反射波強度との比である交差偏波比を算出し、
   速度算出部が、前記反射波強度算出部により算出された前記第１の反射波強度と、前記第２の反射波強度と、前記第３の反射波強度又は前記第４の反射波強度とから、前記送受信装置により偏波が繰り返し送信される周期での各送信偏波における自己相関のラグを算出して、前記ラグから前記送受信装置により繰り返し送信される偏波間の位相差を算出し、前記ラグ、前記偏波間の位相差及びドップラー速度を観測することが可能な範囲であるナイキスト速度から、偏波を反射している観測対象のドップラー速度を算出する
   レーダ信号処理方法。
8. 送受信装置が、水平偏波を送信して水平偏波を受信する第１の送受信処理と、水平偏波を送信して垂直偏波を受信する第３の送受信処理と、垂直偏波を送信して垂直偏波を受信する第２の送受信処理と、垂直偏波を送信して水平偏波を受信する第４の送受信処理とを順番に繰り返し実施し、
   反射波強度算出部が、前記送受信装置により前記第１の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第１の反射波強度を算出し、前記第２の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第２の反射波強度を算出し、前記第３の送受信処理で受信された垂直偏波から、当該垂直偏波の電力値である第３の反射波強度を算出し、前記第４の送受信処理で受信された水平偏波から、当該水平偏波の電力値である第４の反射波強度を算出し、
   交差偏波比算出部が、前記反射波強度算出部により算出された前記第１の反射波強度と前記第３の反射波強度との比、あるいは、前記第２の反射波強度と前記第４の反射波強度との比である交差偏波比を算出し、
   速度算出部が、前記反射波強度算出部により算出された前記第１の反射波強度と、前記第２の反射波強度と、前記第３の反射波強度と、前記第４の反射波強度とから、前記送受信装置により偏波が繰り返し送信される周期での各送信偏波における自己相関のラグを算出して、前記ラグから前記送受信装置により繰り返し送信される偏波間の位相差を算出し、前記ラグ、前記偏波間の位相差及びドップラー速度を観測することが可能な範囲であるナイキスト速度から、偏波を反射している観測対象のドップラー速度を算出する
   レーダ信号処理方法。

Images