JP3794361B2 - レーダ信号処理装置及びレーダ信号処理方法 - Google Patents
レーダ信号処理装置及びレーダ信号処理方法 Download PDFInfo
- Publication number
- JP3794361B2 JP3794361B2 JP2002249336A JP2002249336A JP3794361B2 JP 3794361 B2 JP3794361 B2 JP 3794361B2 JP 2002249336 A JP2002249336 A JP 2002249336A JP 2002249336 A JP2002249336 A JP 2002249336A JP 3794361 B2 JP3794361 B2 JP 3794361B2
- Authority
- JP
- Japan
- Prior art keywords
- power
- echo
- signal processing
- abnormal echo
- power spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
- G01S13/951—Radar or analogous systems specially adapted for specific applications for meteorological use ground based
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Description
【０００１】
【発明の属する技術分野】
この発明は、大気を計測するレーダのレーダ信号処理装置及びレーダ信号処理方法に関するものである。
【０００２】
【従来の技術】
風向・風速の情報は気象予報に必要な情報の一つである。風向・風速を計測する最も一般的な方法は、地上に風向・風速計を設置することであるが、その場合は、地上の風向・風速計では地表付近の風しか測定できない。気象予報をより正確にするためには、上空の風向・風速も知ることが必要であるが、従来はゾンデ等による観測でしか上空の風を計測することができなかった。ゾンデ観測では、ゾンデを上げた時刻のデータしか得ることができないため、観測の時間分解能が数時間以上と低いことが欠点であった。
【０００３】
それに対して、近年ウィンドプロファイラと呼ばれる大気レーダによって上空の風向風速を計測する技術が確立されつつある。ウィンドプロファイラでは、1分〜数分毎に上空の風向・風速を計測することが可能となる。このような高い時間分解能での上空の風情報は、気象予報の精度向上に有効となると期待されている。
【０００４】
ここで、ウィンドプロファイラで大気を計測する原理を説明する。ウィンドプロファイラはドップラレーダの一種であり、一般的に図２０のような構成を持つ。この図において、１０１は空中線、１０２は送受信装置、１０３はレーダ信号処理装置、１０４は品質管理処理装置、１０５は風速ベクトル算出装置、１０６は表示・記録装置である。
【０００５】
次に動作を説明する。送受信装置１０２で生成された電波は、空中線１０１によって空中に放射される。空中に放射された電波は、大気の屈折率の粗密が散乱体となって反射される。反射された電波は空中線によって受信され、送受信装置１０２に入力される。散乱体が上空の風と共に流れていると、受信された電波はドップラ効果によって周波数が変化する。この周波数の変化を一般のドップラレーダと同様に検出することにより、上空の風速を測定する。具体的には、受信された電波は送受信装置１０２によって、増幅、周波数変換、検波されることにより、ビデオ信号に変換され、レーダ信号処理装置１０３へと出力される。レーダ信号処理装置１０３では、ビデオ信号に周波数解析処理を行うことにより、ドップラ周波数を算出する。さらにドップラ周波数からドップラ速度を算出する。
【０００６】
図２１にレーダ信号処理装置１０３の内部構成を表すブロック図を示す。この図において、１はフーリエ変換部、２は周波数領域電力算出部、４はインコヒーレント積分部、５は信号検出部である。フーリエ変換部１では、各距離毎に得られる受信信号を、時間領域信号から周波数領域信号へと変換する。変換には、例えば高速フーリエ変換（Fast Fourier Transform; FFT)を用いれば良いが、その他の手法を用いても構わない。周波数領域信号に変換された信号は周波数領域電力算出部２に入力され、各周波数成分毎に電力を算出することにより、パワースペクトルが得られる。インコヒーレント積分部４では、パワースペクトルを積分することにより、パワースペクトルのゆらぎを低減させる。信号検出部５では、インコヒーレント積分後のパワースペクトルから大気エコーのスペクトルピークを検出し、その周波数からドップラ速度、すなわち視線方向風速などの大気の物理量を算出する。
【０００７】
品質管理処理装置１０４では、レーダ信号処理装置１０３から出力されたドップラ速度データを入力し、大気の時間的・空間的連続性を仮定することにより、時間軸及び空間軸方向で不連続なデータを精度の低いデータと判断し、欠損データとすることにより除去する。そして、除去されずに残ったデータを精度の高いデータとして出力する。
【０００８】
風速ベクトル算出装置１０５は、異なるビーム方向で観測されたドップラ速度データを入力し、それらを合成することにより、３次元風速ベクトルを算出する。図２２は風速ベクトル算出の原理を模式的に示したものである。図２２（ａ）の左側は上空の風速計測を行う際の典型的なビーム方向を示したものである。ここでは、天頂方向及び天頂角１０度程度で東西南北にビームを切り替えるとしている。図２２（ｂ）は図２２（ａ）の状況を東西方向に平行な垂直断面で切って表示したものである。西から東へと風が吹いている場合、西側のビームで観測すれば近づく方向のドップラ速度が、東側のビームで観測すれば遠ざかる方向のドップラ速度が観測される。西側のビームと東側のビームで観測されたドップラ速度の差異を利用することにより、東西方向の風速を求めることができる。同様にして、北側のビームと南側のビームを用いれば、南北方向の風速を得ることができる。さらに天頂方向のドップラ速度は風速の鉛直成分に相当する。その結果、３次元風速ベクトルが得られることになる。
【０００９】
表示・記録装置１０６では、風速ベクトル算出手段から出力された風速ベクトルデータを表示又は記録する。
【００１０】
従来の品質管理処理の方式として、「ウィンドプロファイラによる風観測の定常化に関する業務実験 − 調査結果と提言 −、測候時報 ６５．３、１９９８」で提案されている二次曲面近似方式がある。図２３は、この二次曲面近似方式を実現するための品質管理処理装置の構成を示すブロック図である。この図において、４０１は近似曲面定義範囲設定部、４０２は近似二次曲面算出部、４０３は近似精度判定部、４０４はデータ除去部である。
【００１１】
次に図２３の動作を説明する。近似曲面定義範囲設定部４０１は観測データに対して二次曲面をフィッティングする際の、近似二次曲面を定義する時間範囲と距離範囲を設定する。近似二次曲面算出部４０２では、入力データに対して、近似曲面定義範囲設定部４０１で設定された時間・距離範囲において、二次曲面をフィッティングする。二次曲面のフィッティングを行うには、例えば最小二乗法が用いられる。
【００１２】
近似精度判定部４０３では、近似二次曲面算出部４０２で算出された近似二次曲面の精度が高いかどうかを判定する。例えば前述の文献では、各時間・距離の格子点の観測値と近似二次曲面の差を算出し、差の最大値が予め設定されたしきい値よりも大きいか、又は近似曲面定義範囲内での差の二乗平均の平方根が予め設定されたしきい値よりも大きいかのいずれかが成立する場合に、近似二次曲面の精度が低いと判定される。
【００１３】
近似二次曲面の精度が低い場合、ドップラ速度データと近似二次曲面はデータ除去部４０４に入力される。データ除去部４０４では、近似二次曲面とドップラ速度データの差が大きい格子点において、ドップラ速度データを欠損とすることにより除去する。データ除去部４０４でデータ除去処理されたドップラ速度データは、近似二次曲面算出部４０２に入力され、二次曲面のフィッティングが再度行われる。このような処理を、近似二次曲面の精度が高いと近似精度判定部４０３で判定されるまで繰り返す。
【００１４】
近似二次曲面の精度が高いと判定されると、近似精度判定部４０３は、現時点でのドップラ速度データを品質管理処理済みのデータとして出力する。その後、近似曲面定義範囲設定部４０１は、近似曲面の定義範囲を次の高度又は次の時間へと更新して、同じ処理を繰り返す。
【００１５】
【発明が解決しようとする課題】
以上のような二次曲面近似方式による品質管理処理は、データ品質の劣化した格子点が少ない場合には良好に動作する。例えば、レーダの上空に一羽の鳥が存在するだけであれば、時空間座標の１点だけで鳥の移動速度に対応するドップラ速度が観測され、その他の点では風に対応するドップラ速度、すなわち視線方向風速が観測される。したがって、先の図２３の構成による品質管理装置を用いることにより、鳥エコーだけを除去することができる。
【００１６】
しかし、多数の鳥が時間的に連続して飛来するような場合、品質の劣化した格子点の数が多くなる。すると、時空間座標において、大気エコーよりも渡り鳥エコーが卓越して検出されるようになる。この場合、図２３のような従来の品質管理処理装置による異常エコーの除去は困難となる。このような状況は、渡り鳥が一斉に移動する季節になると頻繁に発生する。
【００１７】
レーダのビームは細いため、鳥がビーム内を通過する時間は限られている。したがって、時間的には断続的に鳥エコーが受信されることになる。鳥エコーが途切れた時間には大気エコーのみを受信することが可能である。しかし、大気に比べて鳥のレーダ断面積は大きいため、時間方向にデータを積分する積分処理を行うと、レーダ断面積の大きい鳥エコーのみが卓越し、大気エコーは埋もれてしまう。これにより、信号処理後のデータでは、鳥エコーは時間的に連続して現れることになる。
【００１８】
本発明は上記のような問題を解決するためになされたものであり、レーダによる大気観測を行う際に、多数の鳥が上空を通過することにより生じる異常エコーのみを除去し、大気エコーのデータを残すことにより、風速計測などの大気観測を正確に行うことができるようにすることを目的とする。
【００１９】
【課題を解決するための手段】
本発明の請求項１に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記異常エコー除去部は、上記パワースペクトルの電力値としてパワースペクトルを周波数軸上で積分した値を用い、この積分値がしきい値を越える場合に異常エコーであると判別し、異常エコーに相当しないパワースペクトルのみ出力することを特徴とするものである。
【００２０】
また、本発明の請求項２に係るレーダ信号処理装置は、上記パワースペクトルの電力値は、パワースペクトルを地形エコー成分が卓越する周波数範囲を除く周波数軸上で積分した値とすることを特徴とする請求項１に記載のものである。
【００２１】
また、本発明の請求項３に係るレーダ信号処理装置は、上記パワースペクトルの電力値は、パワースペクトルを周波数０の成分を除く周波数軸上で積分した値とすることを特徴とする請求項１に記載のものである。
【００２２】
また、本発明の請求項４に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記異常エコー除去部は、上記パワースペクトルの電力値としてパワースペクトルのピーク値を用い、このピーク値がしきい値を越える場合に異常エコーであると判別し、異常エコーに相当しないパワースペクトルのみ出力することを特徴とするものである。
【００２３】
また、本発明の請求項５に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、一定の時間範囲の受信信号の全電力を算出する時間領域電力算出部と、この全電力がしきい値を越える場合に当該時間範囲の受信信号を異常エコーであると判別し異常エコーに相当しない受信信号を出力する異常エコー除去部と、上記異常エコーに相当しない受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備えたことを特徴とするものである。
【００２４】
また、本発明の請求項６に係るレーダ信号処理装置は、日付又は時刻の条件に応じ、異常エコー判定を行うしきい値を設定するしきい値設定部を備え、この設定されたしきい値を用い異常エコー除去判定を行うことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２５】
また、本発明の請求項７に係るレーダ信号処理装置は、日付又は時刻の条件に応じ、異常エコー除去処理を行うか否かを切り替えることができるようにしたことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２６】
また、本発明の請求項８に係るレーダ信号処理装置は、上記日付又は時刻の条件として渡り鳥の飛来の多い季節でかつ夜間を条件としたことを特徴とする請求項６又は７に記載のものである。
【００２７】
また、本発明の請求項９に係るレーダ信号処理装置は、降雨強度に応じ、異常エコー判定を行うしきい値を設定するしきい値設定部を備え、この設定されたしきい値を用いて異常エコー除去判定を行うことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２８】
また、本発明の請求項１０に係るレーダ信号処理装置は、降雨強度に応じ、異常エコー除去処理を行うか否かを切り替えることができるようにしたことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２９】
また、本発明の請求項１１に係るレーダ信号処理装置は、上記信号検出部は鉛直方向の風速を算出し、この風速値に応じて異常エコー判定のしきい値を変化させるようにしたことを特徴とする請求項１乃至１０のいずれか一項に記載のものである。
【００３０】
また、本発明の請求項１２に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記信号検出部が算出した鉛直方向の風速に基づき降雨の存在を判定し、降雨が存在する場合には異常エコー除去処理を行わないことを特徴とするものである。
【００３１】
また、本発明の請求項１３に係るレーダ信号処理装置は、上記信号検出部が算出した鉛直方向の風速に基づき降雨の存在を判定し、降雨が存在する場合には異常エコー除去処理を行わないことを特徴とする請求項１乃至１１のいずれか一項に記載のものである。
【００３２】
また、本発明の請求項１４に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルの電力を周波数軸上で積分することにより信号の全電力を算出する第３のステップと、各距離毎に全電力が予め設定したしきい値を超える場合に信号中に異常エコーが混入していると判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３３】
また、本発明の請求項１５に係るレーダ信号処理方法は、上記第３のステップは、各距離において、地形エコー成分である固定クラッタ領域を除く周波数範囲において、パワースペクトルの電力を周波数軸上で積分することにより、固定クラッタ電力を除く信号の全電力を算出するようにしたものである。
【００３４】
また、本発明の請求項１６に係るレーダ信号処理方法は、上記第３のステップは、周波数０の点のみを固定クラッタ成分とみなし、それ以外の周波数点のデータを積分することにより、固定クラッタ電力を除く信号の全電力を算出するものである。
【００３５】
また、本発明の請求項１７に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルのピーク電力を算出する第３のステップと、各距離毎にピーク電力が予め設定したしきい値を超える場合に信号中に異常エコーが混入していると判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３６】
また、本発明の請求項１８に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルの電力を周波数軸上で積分することにより信号の全電力を算出する第３のステップと、各距離で得られた全電力のうちの最大値が予め設定したしきい値を超える場合に、その時刻の全高度のパワースペクトルデータを使用しないと判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３７】
また、本発明の請求項１９に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルのピーク電力値を算出する第３のステップと、各距離で得られたピーク電力値のうちの最大値が予め設定したしきい値を超える場合に、その時刻の全高度のパワースペクトルデータを使用しないと判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３８】
また、本発明の請求項２０に係るレーダ信号処理方法は、異常エコーの判定に用いるしきい値として、予め設定した固定値を用いるものである。
【００３９】
また、本発明の請求項２１に係るレーダ信号処理方法は、複数時刻で得られたパワースペクトルの電力値を比較して、相対的に高い電力値を持つ時刻以外の電力値を抽出し、抽出した電力値に予め定めた係数を乗算することにより異常エコーの判定に用いるしきい値を得るものである。
【００４０】
また、本発明の請求項２２に係るレーダ信号処理方法は、各距離において、一定時間内の複数の時間で得られたパワースペクトルから算出された全電力又はピーク電力のうち、最低の電力値を大気エコーの電力として抽出し、この抽出された最低電力値に予め設定した係数を乗算して得られる値をしきい値として設定するようにしたものである。
【００４１】
また、本発明の請求項２３に係るレーダ信号処理方法は、異常エコーを含まないデータが異常エコー判定しきい値を超える確率がある一定値以下となるように、最低電力値に乗算する係数を設定するものである。
【００４２】
また、本発明の請求項２４に係るレーダ信号処理方法は、異常エコーとして除去されたデータ数によらず、一定の積分時間でインコヒーレント積分処理を行うものである。
【００４３】
また、本発明の請求項２５に係るレーダ信号処理方法は、異常エコーとして除去されたデータ数によらず、インコヒーレント積分数が一定となるようにインコヒーレント積分処理を行うものである。
【００４４】
【発明の実施の形態】
実施の形態１．
図１はこの発明の実施の形態１のレーダ信号処理装置の構成を示すブロック図である。図において、１はフーリエ変換部、２は周波数領域電力算出部、３は異常エコー除去部、４はインコヒーレント積分部、５は信号検出部である。
【００４５】
次に動作について説明する。フーリエ変換部１では、各距離毎に得られる受信信号を、時間領域信号から周波数領域信号へと変換する。変換には、例えば高速フーリエ変換（Fast Fourier Transform; FFT)を用いれば良いが、その他の手法を用いても構わない。周波数領域信号に変換された信号は周波数領域電力算出部２に入力され、各周波数成分毎に電力を算出することにより、パワースペクトルが得られる。
【００４６】
異常エコー除去部３では、周波数領域電力算出部２で得られたパワースペクトルの電力値が、設定されたしきい値よりも大きい場合に、異常エコーが混入していると判定する。異常エコーが混入していると判定されたパワースペクトルは除去する。インコヒーレント積分部４では、異常エコー除去部３で除去されずに残ったパワースペクトルを積分することにより、パワースペクトルのゆらぎを低減させる。信号検出部５では、インコヒーレント積分後のパワースペクトルから大気エコーのスペクトルピークを検出し、その周波数からドップラ速度、すなわち視線方向風速を算出する。
【００４７】
なお、異常エコー除去部３におけるパワースペクトルの除去の方法としては、除去するパワースペクトルを欠損とする、すなわち「除去」のフラグを立てることが考えられる。インコヒーレント積分部４では、「除去」のフラグのついていないパワースペクトルのみを用いて積分処理を行うことになる。あるいは、異常エコー除去部３におけるパワースペクトルの除去の他の方法として、除去するパワースペクトルの電力値を全て０とすることも考えられる。この場合、インコヒーレント積分部４では、異常エコー除去部３でパワースペクトルが除去されたか否かを考慮せずに積分処理を行っても、異常エコーは結果的に除去できることになる。
【００４８】
異常エコー判定に用いるパワースペクトルの電力値としては、例えばパワースペクトルを周波数軸上で積分した全電力を用いることが考えられる。図２はそのような場合の異常エコー除去部３の内部構成を示すブロック図である。この図において、３０１は電力積分部、３０２は電力判定部である。電力積分部３０１では、パワースペクトルを周波数軸上で積分することにより、パワースペクトルの全電力を算出する。電力判定部３０２では、電力積分部３０１で算出された全電力を、予め設定したしきい値と比較し、全電力の方が大きければパワースペクトルに異常エコー成分が含まれると判定し、このパワースペクトルを除去する。逆に全電力の方が小さければ、パワースペクトルをそのままインコヒーレント積分部４へと出力する。
【００４９】
ただし、近距離のデータでは、地形エコー成分の電力が大きい場合がある。地形エコー成分は周波数０近辺のみに存在し、スペクトル幅も大気エコーに比べて明らかに狭いため、両者を識別することは比較的容易である。信号検出部５にて地形エコーの影響を取り除くことも通常良く行われることから、異常エコー除去部３においては地形エコー以外の異常エコー、すなわち鳥エコーのみを除去するようにしても良い。この場合、電力積分部３０１では、周波数０の近傍を除く周波数範囲でパワースペクトルを積分するようにすれば良い。
【００５０】
地形エコー成分を取り除く最も簡易な方法としては、周波数０成分のみを除いてパワースペクトルの積分値を求めることである。図３はそのような動作をする異常エコー除去部の構成を示す図である。この図において、３０３は周波数０成分減算部である。その他の符号は前述のものと同じである。周波数０成分減算部３０３では、電力積分部３０１において全周波数範囲でパワースペクトルを積分することにより得られた電力値から、周波数０成分を減算する。これにより、周波数０の成分を除く電力積分値を得ることができる。
【００５１】
以上では、パワースペクトルを積分することにより得られる電力を用いて異常エコー判定を行う例を示したが、図４では、ピーク電力値を用いて判定する場合の異常エコー除去部３の内部構成を示すブロック図である。この図において、３０４はピーク電力抽出部である。ピーク電力抽出部では、パワースペクトル中の最大の電力値を抽出し、この値を信号電力の代表値と見なして、３０２電力判定部にて異常エコー判定を行うようにしている。
【００５２】
図５は本発明の実施の形態のレーダ信号処理方法の流れを示す図である。この図では、パワースペクトルを周波数軸上で積分する場合、すなわち図２あるいは図３の異常エコー除去部の動作に対応したレーダ信号処理方法の流れを示している。
【００５３】
ステップｓｔ００１とステップｓｔ０１２で囲まれた区間は、インコヒーレント積分を行うループを表すものである。ステップ００２では、フーリエ変換を行う際の時間点数分の受信信号を入力する。ステップｓｔ００３からステップｓｔ００７で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。ｓｔ００４では、各距離において、時系列信号である受信信号にフーリエ変換を施すことにより、周波数領域の信号へと変換する。ｓｔ００５では、各周波数成分毎に電力値を算出することにより、パワースペクトルを得る。ｓｔ００６では、パワースペクトルを周波数軸上で積分することにより、全電力を算出する。周波数軸上の積分では、図２の異常エコー除去部３の動作のように、周波数０の近傍を除く周波数領域で積分するようにしても良いし、あるいは図３の異常エコー除去部３のように、周波数０の成分のみを除いた積分を行うようにしても良い。
【００５４】
ステップｓｔ００８からステップｓｔ０１１で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。ステップｓｔ００９では、各距離毎に、全電力と予め設定されたしきい値と比較し、全電力がしきい値を下回る場合にのみ、ステップｓｔ０１０にて処理中のパワースペクトルを積算する。インコヒーレント積分ループは、設定された回数Ｎinだけｓｔ００２からｓｔ０１１の処理を行うと終了するようにする。この場合、ｓｔ０１０で積算されるパワースペクトルの数、すなわちインコヒーレント積分数は最大でＮinとなり、異常エコーの影響を受ける場合には、これよりも小さいインコヒーレント積分数となる。インコヒーレント積分ループを、設定された回数Ｎinだけ行うようにした場合、インコヒーレント積分後のパワースペクトルの得られる時間間隔が、実際に実行したインコヒーレント積分数によらず一定となるため、品質管理処理装置によるデータ品質管理処理などの信号処理より後段のデータ処理を行う上で、入力データの扱いが容易になる。ただし、インコヒーレント積分数が少なくなった場合には、それだけ信号検出の性能が劣化することになる。
【００５５】
あるいは、インコヒーレント積分数が所定数以上になった場合にステップｓｔ０１２でインコヒーレント積分ループを終了するようにしても良い。この場合には、常にインコヒーレント積分数を一定以上にすることができるため、一定以上の信号検出性能を常に確保することができる。
【００５６】
なお、図５では、ステップｓｔ００１からステップｓｔ００７と、ステップｓｔ００８からステップｓｔ０１１の２つのレンジループを持っているが、これを一つにまとめるようにしても良い。あるいは処理が等価であればレンジループをさらに細かく分けるようにしても良い。このように、処理が等価であれば、本発明はループの区切り方を限定するものではない。
【００５７】
図５に示した本発明のレーダ信号処理方法の流れを示す図では、レーダのビーム方向の切り替えについては特に表現されていない。実際のレーダでは、例えば図６に示すように、１つの時間範囲でのフーリエ変換が終わると次のビーム方向に切り替え、必要な観測方向を全てカバーするというようにインコヒーレント積分を行う時間範囲の途中でビームを切り替えるような手順が取られることが多い。このような場合にも、個々のビーム方向のデータのみに着目すれば、図５の処理の流れは変わらない。よって、本発明は、このようなビーム切り替えの手順を特に限定するものではない。
【００５８】
異常エコー判定で用いるしきい値については、大気のレーダ断面積と鳥のレーダ断面積を考慮して固定的に設定すれば良い。一般に鳥のレーダ断面積は大気のレーダ断面積と比較して十分大きい。そこで、レーダ方程式から予測される大気エコーの受信電力と鳥エコーの受信電力の間の電力値をしきい値として採用すれば良い。
【００５９】
図７は、パワースペクトルのピーク電力を用いて異常エコー判定を行う場合、すなわち図４の異常エコー除去部の動作に対応したレーダ信号処理方法の流れを示す図である。この図では、ステップ００６ｂでパワースペクトルのピーク電力を抽出するようにし、そのピーク電力についてしきい値による判定をステップ００９ｂで行うようにしている以外は、図５に示したレーダ信号処理方法と同じである。
【００６０】
以上のように、本発明の実施の形態では、エコー強度によって異常エコーと大気エコーを識別し、大気エコーのみを用いて信号検出を行うようにしているため、渡り鳥が上空を飛来するような場合にも正確に大気観測を行うことが可能になる。特に、インコヒーレント積分前のデータに対して異常エコーの除去を行うようにしているため、従来のレーダ信号処理装置又は品質管理処理装置に比べて、異常エコー除去の処理を行う時間単位が短いため、断続的に異常エコーが発生している状況において、正常エコーが異常エコーに埋もれる前に異常エコーを除去し、正常エコーのみを残すことが可能となる。
【００６１】
また、パワースペクトルを算出するために算出した周波数領域上の電力値を用いるため、後述の実施の形態７と比較して時間領域での電力値算出が不要となるため、実施の形態７の場合よりも信号処理演算量は少なくて済む。
【００６２】
なお、本実施の形態では、風速計測用レーダへの適用を想定したレーダ信号処理装置及びレーダ信号処理方法について記述しているが、本レーダ信号処理装置及びレーダ信号処理方法は、電波を送受信するレーダだけに限定されるものではなく、音波を送受信するソーダや、光波を送受信するレーザレーダに対しても全く同様に適用することができる。また、風速計測を目的としたレーダだけでなく、降雨レーダを始めとする大気観測用レーダにも適用できる。このことは、以後に記述する別の実施の形態においても同様である。
【００６３】
実施の形態２．
以上の実施の形態１では、異常エコー判定及び異常エコー除去を距離毎に行うようにしたものであるが、次に、異常エコー判定及び異常エコー除去を全レンジで一括して行うような場合の実施の形態を示す。本発明の実施の形態において、レーダ信号処理装置の構成は、実施の形態１と同様に図１で表されるものであるが、レーダ信号処理方法が異なるものとなる。
【００６４】
図８は本発明の実施の形態２のレーダ信号処理方法の流れを示すものである。ステップｓｔ１０１とステップｓｔ１１３で囲まれた区間は、インコヒーレント積分を行うループを表すものである。ステップ１０２では、フーリエ変換を行う際の時間点数分の受信信号を入力する。ステップｓｔ１０３からステップｓｔ１０７で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。ｓｔ１０４では、各距離において、時系列信号である受信信号にフーリエ変換を施すことにより、周波数領域の信号へと変換する。ｓｔ１０５では、各周波数成分毎に電力値を算出することにより、パワースペクトルを算出する。ｓｔ１０６では、パワースペクトルを周波数軸上で積分することにより、全電力を算出する。周波数軸上の積分では、図２の異常エコー除去部３の動作のように、周波数０の近傍を除く周波数領域で積分するようにしても良いし、あるいは図３の異常エコー除去部３のように、周波数０の成分のみを除いた積分を行うようにしても良い。ステップｓｔ１０８では、ステップｓｔ１０６でレンジ毎に算出された全電力のうちの最大値を抽出する。ステップｓｔ１０９では、ステップｓｔ１０８でえられた最大全電力と予め設定したしきい値とを比較し、最大全電力の方が小さい場合にのみ、ステップｓｔ１１０からステップｓｔ１１２で表されるレンジループの処理を行う。レンジループ内のステップｓｔ１１１では、ステップｓｔ１０５で得られたパワースペクトルを積分後パワースペクトルに加算することにより、インコヒーレント積分を実行する。
【００６５】
このように図８のレーダ信号処理方法では、ステップｓｔ１０９のしきい値判定で異常エコーが受信信号に混入していると判定された場合には、全距離でパワースペクトルの積算が行われないため、全距離で一括して異常エコーが除去されることになる。
【００６６】
図８では、異常エコー判定を全電力に対して行うのに対して、図９のようにピーク電力に対して異常エコー判定を行うことも考えられる。すなわち図４の異常エコー除去部の動作に対応したレーダ信号処理方法の流れを示す図である。この図では、ステップ１０６ｂでパワースペクトルのピーク電力を抽出するようにし、そのピーク電力の最大値をステップ１０８ｂで抽出し、抽出された最大ピーク電力についてしきい値による判定をステップ１０９ｂで行うようにしている以外は、図８に示したレーダ信号処理方法と同じである。
【００６７】
本発明の実施の形態のレーダ信号処理方法によれば、全距離で一括して異常エコー判定を行うようにしているため、距離毎に異常エコー判定を行うよりも演算時間を短縮することができる。特に、以前から用いられているレーダ信号処理装置に異常エコー除去機能を加える場合のように、信号処理ハードウェアの性能に制約があるような場合に有効である。
【００６８】
実施の形態３．
以上の実施の形態では、異常エコー判定に用いるしきい値は固定したものを用いるとしていたが、本実施の形態では、しきい値を観測データから決定するような実施形態について述べる。
【００６９】
図１０は本発明の実施の形態におけるレーダ信号処理装置の構成を示すブロック図である。この図において、７はしきい値設定部である。その他の符号は前述のものと同じである。次に動作を説明する。しきい値設定部７は周波数領域電力算出部２で得られた各距離のパワースペクトルデータを入力し、複数時間分だけ蓄積する。各レンジ毎に、蓄積されたパワースペクトルの中から異常エコーが含まれないと判断されるデータを抽出し、その電力値をもとに異常エコー判定のしきい値を設定する。具体的には、電力値から判断して異常エコーが含まれないと判断される時刻のデータを選択し、その電力値に予め定義した方法でマージンを加えた電力値を、異常エコー判定しきい値として異常エコー除去部３へ出力する。このしきい値を用いて、異常エコー除去部３は周波数領域電力算出部２から出力されたパワースペクトルを評価し、異常エコー成分が含まれると判断されたデータを除去する。その他の動作は前述の図１のレーダ信号処理装置と同じである。
【００７０】
次に図１０のレーダ信号処理装置の動作を、処理アルゴリズムという観点から説明する。本実施の形態においても、しきい値を設定する方法以外の処理アルゴリズムについては、実施の形態１又は実施の形態２のレーダ信号処理方法と同じである。ここでは、レーダ信号処理方法のアルゴリズムのうちの、異常エコー判定に用いるしきい値を設定する部分のみを特に説明する。
【００７１】
図１１は本発明の実施の形態のレーダ信号処理方法における異常エコー判定しきい値の設定部分の流れを示す図である。この図において、ステップｓｔ２０１からステップｓｔ２０４で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。
【００７２】
例えば実施の形態１で示した図５では、ステップｓｔ００３からステップｓｔ００７の処理において、レンジ毎のフーリエ変換及びパワースペクトル算出を複数の時間ブロックで行うことにより、各レンジにおいて複数時刻でのパワースペクトル及び全電力値が得られている。
【００７３】
図１１のｓｔ２０２では、複数時刻の全電力値から、最小の全電力値を抽出する。最小電力値となっている時間は、鳥による異常エコーが混入している可能性が最も低いと考えられる。そこで、最小全電力値が得られている時刻の観測データが異常エコーを含まないデータのサンプルであると見なして、ステップｓｔ２０３ではその全電力値にある係数を乗算した電力値を異常エコー判定のしきい値を設定する。乗算する係数の値は経験的に固定値を設定しても良い。あるいは大気エコーの電力値の統計的性質、すなわち確率密度分布を考慮して、異常エコーを含まない大気エコーがしきい値を超える確率、即ち誤って正常エコーが除去される確率が一定となるように設定しても良い。図１１における全電力としては、全周波数成分から計算した全電力を用いても良いし、あるいは、ドップラー周波数が０近傍である地形エコー成分を除いた全電力を用いても良い。
【００７４】
以上の方法によってしきい値設定を行えば、例えば大気エコーの受信強度が弱い場合には、異常エコーの影響を受けやすくなるが、異常エコー判定しきい値が低くなるため、より多くの異常エコー成分を除去することが可能となる。逆に大気エコーの受信強度が比較的強い場合には、しきい値が大きくなるため、大気エコーが強いために誤って異常エコーと判定される確率が低くなるため、多くの大気エコー成分を利用することが可能となる。
【００７５】
図１１では全電力値を利用して異常エコー判定しきい値を設定していたが、図１２のようにパワースペクトルのピーク値を用いて異常エコー判定しきい値を設定しても良い。図１２では、図１１のステップｓｔ２０２の代わりにステップｓｔ２０２ｂを用いることにより、複数時刻のピーク電力値から、最小のピーク電力値を抽出する。また、ステップｓｔ２０３の代わりにステップｓｔ２０３ｂを用いることにより、ピーク電力値にマージンを加算した電力レベルを異常エコー判定のしきい値を設定する。全電力の代わりにピーク電力を用いること以外は、図１１と図１２は同様の処理により異常エコー検出しきい値を設定する。
【００７６】
図１１又は図１２のしきい値設定処理は、レーダ信号処理の流れの中の適当なところにおけば良い。例えば図５に示したレーダ信号処理方法を用いる場合は、ｓｔ００７の直後に入れることが考えられる。あるいは、ｓｔ０１２の直後にしきい値設定処理を入れ、設定されたしきい値は次の時刻の信号処理に用いるようにしても良い。本発明では、しきい値処理を入れる段階を特に制限するものではない。
【００７７】
本実施の形態によれば、大気エコーのエコー強度に応じたしきい値を用いて異常エコー判定を行うようにしているため、誤って大気エコーを除去する確率を低くしたまま、異常エコーを正確に除去することが可能となり、観測精度を高く保つことができる。
【００７８】
実施の形態４．
前述の実施の形態３では、異常エコー判定しきい値を観測されたデータを用いて変更していたが、渡り鳥飛来の頻度は季節と時刻に依存することから、日時によって判定しきい値を変更するようにしても良い。具体的には、渡り鳥が移動する期間は毎年一定しており、しかも外敵から目立たない夜間に移動することが多い。そのため、渡り鳥が移動する期間、例えば春と秋の夜間に、異常エコー判定しきい値を低くし、異常エコーの残留が少なくなるようにすれば良い。その他の日時では、鳥が観測範囲内を通過するとしても、その頻度はまれであるため、従来から行われている信号処理後の品質管理処理にて除去できるため、異常エコー判定しきい値を高めにして、正常エコーを誤って除去する確率が低くなるようにすることが考えられる。
【００７９】
そのようなレーダ信号処理装置の構成例のブロック図を図１３に示す。この図において、８は時計である。その他の符号は前述のものと同じである。時計８から出力される日時のデータはしきい値設定部７へ渡される。しきい値設定部７では、日時毎に予め設定されたしきい値を異常エコー除去部３へ出力する。異常エコー除去部３はしきい値設定部７から入力したしきい値を用いて、周波数領域電力算出部２から入力したパワースペクトルに対して異常エコーか否かの判定を行い、異常エコーであると判定されれば、そのパワースペクトルデータを除去する。その他の動作については、図１０に示したレーダ信号処理装置と同じである。
【００８０】
図１３のレーダ信号処理装置の他の例を図１４に示す。この図において、９は処理切替スイッチである。その他の符号は前述のものと同じである。図１３では、日時によってしきい値を変化させていたが、図１４では異常エコー除去部３による処理を、日時によってオン／オフするようにしている。例えば、渡りの季節の夜間のみ処理切替スイッチをオンにして、異常エコー除去部３による処理を行い、それ以外の日時では処理切替スイッチをオフにして、異常エコー除去部３による処理を行わないようにする。
【００８１】
以上のように、日時によってしきい値を変更することにより、特に異常エコーが良く除去できるようにしたため、あるいは異常エコー除去処理をオン／オフすることにより渡り鳥による異常エコーの影響が大きく現れるときだけ異常エコーを除去するようにした。よって、渡り鳥の影響が小さい日時において、正常エコーが誤って除去されることがなくなり、インコヒーレント積分数が少なくならないため、正常エコーの検出確率が減らず、データ取得率を維持することが可能となる。
【００８２】
実施の形態５．
前述の実施の形態４では、日時によってレーダ信号処理方法を変更するようにしていたが、渡り鳥の渡りは天候にも左右される。本発明の実施の形態５では、天候によってレーダ信号処理方法を変化させるようにしたものである。図１５は本発明の実施の形態のレーダ信号処理装置の構成を示すブロック図である。この図の符号は全て前述のものと同じである。図１５では、しきい値設定部７が雨量データを入力し、降雨強度によって変化するしきい値を異常エコー除去部３へと出力する。雨量の強いときには、渡り鳥が飛来する可能性が低いため、しきい値を大きく設定することにより、正常エコーを誤って除去しないようにする。逆に降雨のないときには、渡り鳥が飛来する可能性があるため、しきい値を小さく設定することにより、異常エコーが除去できるようにする。
【００８３】
図１６に示すブロック図は図１５を変形したもので、雨量データによって、異常エコー除去部３による処理をオン／オフしている。すなわち、降雨のあるときには異常エコー除去部３の処理をオフにし、降雨のないときは逆にオンにするようにしている。
【００８４】
雨量データとしては、一般的な雨量計で観測されるものを用いてもよいし、あるいは風速計測用のレーダとは別に設置された降雨レーダで計測された雨量強度を用いても良いが、本発明では特に雨量データの観測方法を限定するものではない。
【００８５】
また、渡り鳥飛来の特性が雨量以外の気象データに依存する場合には、そのようなデータも合わせてしきい値設定部７に入力するようにしても良い。
【００８６】
本発明の実施の形態によれば、降雨の存在するときには渡り鳥が飛来しないと判断して信号処理を制御するため、降雨時に誤って正常エコーを除去することがなくなる。
【００８７】
実施の形態６．
前述の実施の形態５では、雨量データを用いて、降雨時の異常エコー除去処理を制御していたが、本発明の実施の形態では、レーダ装置で観測されたデータから降雨の有無を判定するような実施形態を説明する。
【００８８】
図１７及び図１８は本発明の実施の形態のレーダ信号処理装置の構成を示す図である。この図に含まれる符号は全て前述のものと同じである。次に動作を説明する。信号検出部から出力されるドップラ速度データのうち、ビームを天頂方向に向けて観測されたドップラ速度を処理切替スイッチ９へ入力する。降雨が存在しない場合、一般的に鉛直方向の風は１〜２ｍ／ｓ程度である。それに対し、雨滴の落下速度はそれよりも大きくなる。また、渡り鳥はほとんどの場合に水平飛行をすると考えられる。したがって、天頂ビームのドップラ速度が近づく方向に大きい場合には、降雨が存在し、渡り鳥は上空を通過していないと判断することができる。そこで、図１７のレーダ信号処理装置では、ドップラ速度が近づく方向に大きい場合に、異常エコー判定しきい値を大きくし、あるいは図１８のレーダ信号処理装置では、ドップラ速度が近づく方向に大きい場合に、異常エコー判定部３の処理をオフにする。
【００８９】
なお、鉛直方向の風速を得るのには、必ずしも天頂方向へビームを向けて観測しなくても良い。例えば、天頂角θで東方向と西方向を観測して得られるドップラ速度をVe とVw とすると、風速の鉛直成分ｗは下記式（１）で算出することができる。
ｗ＝(Ve + Vw)/２cosθ （１）
式（１）のような計算は信号検出部５で行うとしても良いし、あるいは図２０の風速ベクトル算出装置で算出される風速ベクトルの鉛直成分を利用するようにしても良い。
【００９０】
本発明の実施の形態によれば、別の装置で観測された雨量データを入力しなくても、降雨の有無を判定し、その結果で異常エコー除去部の動作を制御するようにしているため、他の装置からのデータ入力が不要となり、簡易な構成で適切な異常エコー除去処理が可能となる。
【００９１】
実施の形態７．
以上に述べた実施の形態では、時系列領域の受信信号を周波数領域に変換した後のデータに対して、異常エコーの判定を行っていた。ただし、異常エコーの有無の判定は電力値で行っているため、周波数領域に変換する前の時系列データに対して行うことも可能である。図１９は本発明の実施の形態のレーダ信号処理装置の構成を示すブロック図である。この図において、７０１は時間領域電力算出部である。
【００９２】
次に動作を説明する。時間領域電力算出部７０１では、時系列データである受信データを設定した時間範囲で区切り、その時間範囲における全電力を算出する。異常エコー除去部１では、時間領域電力算出部７０１から出力された電力を予め設定したしきい値と比較し、全電力の方がしきい値よりも大きければ、その時間範囲の時系列データを受信信号から除去する。除去する場合には、データを０としても良いし、あるいは欠損データであるというマークをつけておいても良い。それ以外の処理、すなわちフーリエ変換部２以降の処理については、実施の形態１などの前述の実施の形態の動作と同じである。
【００９３】
なお、時間領域電力算出部７０１の全電力算出における時間範囲の設定は、フーリエ変換部２の処理においてフーリエ変換を実行する際の時間窓と同じとしても良いし、あるいは全く別の時間範囲設定としても良い。同じ時間範囲設定とすれば、本実施の形態のように時間領域で全電力を算出しても、あるいは前述の実施の形態のように周波数領域で全電力を算出するようにしても、両者は同等となる。逆に時間範囲設定をフーリエ変換の時間窓と異なるように設定することも、本実施の形態では可能となるため、例えばフーリエ変換処理の時間長よりも短い時間で異常エコーが断続して生じるような場合にも、全電力を算出する時間範囲設定を短く設定することにより、異常エコーの存在する時間のデータだけを的確に除去することが可能であり、正常データの除去を最低限に抑えることが可能となる。
【００９４】
なお、全電力の代わりに、区切った時間範囲内の平均電力を時間領域全電力算出部７０１で算出しても良い。何故なら、全電力を区切った時間範囲の時間点数で割ったものであり、異常エコー判定のしきい値を時間点数で割ったものを用いれば、除去処理の原理としては実質的な差異がないからである。
【００９５】
【発明の効果】
本発明の請求項１によれば、受信信号の電力値として、周波数軸上で得られているパワースペクトルを積分した電力値を用い、インコヒーレント積分前の電力値をしきい値と比較することにより、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【００９６】
また、本発明の請求項２によれば、地形エコーが卓越する低周波領域を除いて全電力値を算出するようにしたため、強い地形エコーが混入している場合にも、渡り鳥による異常エコーの混入の有無の判定を正確に行うことが可能となる。
【００９７】
また、本発明の請求項３によれば、全ドップラ周波数範囲でパワースペクトルを積分して得られる全電力値から、周波数０の点の電力値のみを除くことにより、地形エコーを除いた全電力値を算出するようにしたため、簡易な処理で上記請求項３と同様の効果を得ることが可能となる。
【００９８】
また、本発明の請求項４によれば、受信信号の電力値として、周波数軸上で得られているパワースペクトルのピーク値を抽出したものを用いることにより、上記請求項１と同様の効果を得ることができる。
【００９９】
また、本発明の請求項５によれば、時系列データから算出した電力値をしきい値と比較することにより、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１００】
また、本発明の請求項６によれば、季節や時刻によって異常エコー判定のしきい値を調整することにより、渡り鳥の影響が小さい日時において、正常エコーが誤って除去されることがなくなり、インコヒーレント積分数が少なくならないため、正常エコーの検出確率を維持することができ、その結果としてデータ取得率を維持することが可能となる。
【０１０１】
また、本発明の請求項７によれば、季節や時刻によって異常エコー判定をオン／オフすることにより、渡り鳥の影響が小さい日時において、正常エコーが誤って除去されることがなくなり、インコヒーレント積分数が少なくならないため、正常エコーの検出確率を維持することができ、その結果として減らず、データ取得率を維持することが可能となる。
【０１０２】
また、本発明の請求項８によれば、渡り鳥の特性を考慮して、渡り鳥の飛来する季節の夜間に異常エコー処理をオンにし、それ以外の場合に異常エコー処理をオフとすることにより、上記請求項７と同様の効果を得ることが可能となる。
【０１０３】
また、本発明の請求項９によれば、降雨強度に応じて異常エコー判定のしきい値を調整するようにしたため、降雨時に誤って正常エコーを除去することがなくなるという効果がある。
【０１０４】
また、本発明の請求項１０によれば、降雨の存在するときには渡り鳥が飛来しないと判断してレーダ信号処理を制御するため、降雨時に誤って正常エコーを除去することがなくなるという効果がある。
【０１０５】
また、本発明の請求項１１によれば、レーダ装置で観測された雨滴の落下速度に応じて異常エコー判定のしきい値を調整するようにしたため、他の装置からのデータを入力することなく、降雨時に誤って正常エコーを除去することのないレーダ信号処理が可能となる。
【０１０６】
また、本発明の請求項１２によれば、レーダ装置の観測データから降雨の有無を判定し、その結果で異常エコー除去部の動作を制御するようにしているため、他の装置からのデータを入力することなく、降雨時に誤って正常エコーを除去することのないレーダ信号処理が可能となる。
【０１０７】
また、本発明の請求項１３によれば、レーダ装置の観測データから降雨の有無を判定し、その結果で異常エコー除去部の動作を制御するようにしているため、他の装置からのデータを入力することなく、降雨時に誤って正常エコーを除去することのないレーダ信号処理が可能となる。
【０１０８】
また、本発明の請求項１４によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１０９】
また、本発明の請求項１５によれば、強い地形エコーが混入している場合にも、渡り鳥による異常エコーの混入の有無の判定を正確に行い、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１０】
また、本発明の請求項１６によれば、全ドップラ周波数範囲でパワースペクトルを積分して得られる全電力値から、周波数０の点の電力値のみを除き、地形エコーを除いた全電力値を算出することにより、強い地形エコーが混入している場合にも、簡易な処理で渡り鳥による異常エコーの混入の有無の判定を正確に行い、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１１】
また、本発明の請求項１７によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１２】
また、本発明の請求項１８によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となり、かつ全距離で一括して異常エコー判定を行うようにしているため、距離毎に異常エコー判定を行うよりも演算時間を短縮することができる。
【０１１３】
また、本発明の請求項１９によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となり、かつ全距離で一括して異常エコー判定を行うようにしているため、距離毎に異常エコー判定を行うよりも演算時間を短縮することができる。
【０１１４】
また、本発明の請求項２０によれば、異常エコーの判定に用いるしきい値として、予め設定した固定値を用いることにより、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１５】
また、本発明の請求項２１によれば、大気エコーのエコー強度に応じたしきい値を用いて異常エコー判定を行うようにしているため、誤って大気エコーを除去する確率を低くしたまま、異常エコーを正確に除去することが可能となり、観測精度を高く保つことができる。
【０１１６】
また、本発明の請求項２２によれば、大気エコーのエコー強度に応じたしきい値を用いて異常エコー判定を行うようにしているため、誤って大気エコーを除去する確率を低くしたまま、異常エコーを正確に除去することが可能となり、観測精度を高く保つことができる。
【０１１７】
また、本発明の請求項２３によれば、異常エコーを含まないデータが異常エコー判定しきい値を超える確率がある一定値以下となるように、最低電力値に乗算する係数を設定することにより、正常エコーの検出確率を一定に保ったまま、異常エコーを除去により観測精度を高く保つことができる。
【０１１８】
また、本発明の請求項２４によれば、異常エコーとして除去されたデータ数によらず、一定の積分時間でインコヒーレント積分処理を行うようにしたため、品質管理処理装置によるデータ品質管理処理などのレーダ信号処理より後段のデータ処理を行う上で、入力データの扱いが容易になる。
【０１１９】
また、本発明の請求項２５によれば、異常エコーとして除去されたデータ数によらず、インコヒーレント積分数が一定となるようにインコヒーレント積分処理を行うようにしたため、一定以上の信号検出性能を常に確保することができる。
【図面の簡単な説明】
【図１】 本発明の実施の形態１のレーダ信号処理装置の構成を示すブロック図。
【図２】 本発明の実施の形態１のレーダ信号処理装置の異常エコー除去部の構成を示すブロック図。
【図３】 本発明の実施の形態１のレーダ信号処理装置の異常エコー除去部の別の構成を示すブロック図。
【図４】 本発明の実施の形態１のレーダ信号処理装置の異常エコー除去部の別の構成を示すブロック図。
【図５】 本発明の実施の形態１のレーダ信号処理方法の流れを示す図。
【図６】 風速観測における典型的な観測手順を示す図。
【図７】 本発明の実施の形態１の別のレーダ信号処理方法の流れを示す図。
【図８】 本発明の実施の形態２のレーダ信号処理方法の流れを示す図。
【図９】 本発明の実施の形態２の別のレーダ信号処理方法の流れを示す図。
【図１０】 本発明の実施の形態３のレーダ信号処理装置の構成を示すブロック図。
【図１１】 本発明の実施の形態３におけるレーダ信号処理方法におけるしきい値設定アルゴリズムの流れを示す図。
【図１２】 本発明の実施の形態３におけるレーダ信号処理方法における別のしきい値設定アルゴリズムの流れを示す図。
【図１３】 本発明の実施の形態４のレーダ信号処理装置の構成を示すブロック図。
【図１４】 本発明の実施の形態４の別のレーダ信号処理装置の構成を示すブロック図。
【図１５】 本発明の実施の形態５のレーダ信号処理装置の構成を示すブロック図。
【図１６】 本発明の実施の形態５の別のレーダ信号処理装置の構成を示すブロック図。
【図１７】 本発明の実施の形態６のレーダ信号処理装置の構成を示すブロック図。
【図１８】 本発明の実施の形態６の別のレーダ信号処理装置の構成を示すブロック図。
【図１９】 本発明の実施の形態７のレーダ信号処理装置の構成を示すブロック図。
【図２０】 一般的なウィンドプロファイラの構成を示すブロック図
【図２１】 従来のウィンドプロファイラのレーダ信号処理装置の構成を示すブロック図。
【図２２】 ウィンドプロファイラにおける風速ベクトル算出の原理を示す図。
【図２３】 従来のウィンドプロファイラの品質管理処理装置の構成例を示すブロック図。
【符号の説明】
１ フーリエ変換部、 ２ 周波数領域電力算出部、 ３ 異常エコー除去部、 ４ インコヒーレント積分部、 ５ 信号検出部、 ７ しきい値設定部、８ 時計、 ９ 処理切替スイッチ、 １０１ 空中線、 １０２ 送受信装置、 １０３ レーダ信号処理装置、 １０４ 品質管理処理装置、 １０５ 風速ベクトル算出装置、 １０６ 表示記録装置、 ３０１ 電力積分部、 ３０２ 電力判定部、 ３０３ 周波数０成分減算部、 ３０４ ピーク電力抽出部、 ４０１ 近似曲面定義範囲部、 ４０２ 近似二次曲面算出部、 ４０３近似精度判定部、 ４０４ データ除去部、 ７０１ 時間領域電力算出部。
【発明の属する技術分野】
この発明は、大気を計測するレーダのレーダ信号処理装置及びレーダ信号処理方法に関するものである。
【０００２】
【従来の技術】
風向・風速の情報は気象予報に必要な情報の一つである。風向・風速を計測する最も一般的な方法は、地上に風向・風速計を設置することであるが、その場合は、地上の風向・風速計では地表付近の風しか測定できない。気象予報をより正確にするためには、上空の風向・風速も知ることが必要であるが、従来はゾンデ等による観測でしか上空の風を計測することができなかった。ゾンデ観測では、ゾンデを上げた時刻のデータしか得ることができないため、観測の時間分解能が数時間以上と低いことが欠点であった。
【０００３】
それに対して、近年ウィンドプロファイラと呼ばれる大気レーダによって上空の風向風速を計測する技術が確立されつつある。ウィンドプロファイラでは、1分〜数分毎に上空の風向・風速を計測することが可能となる。このような高い時間分解能での上空の風情報は、気象予報の精度向上に有効となると期待されている。
【０００４】
ここで、ウィンドプロファイラで大気を計測する原理を説明する。ウィンドプロファイラはドップラレーダの一種であり、一般的に図２０のような構成を持つ。この図において、１０１は空中線、１０２は送受信装置、１０３はレーダ信号処理装置、１０４は品質管理処理装置、１０５は風速ベクトル算出装置、１０６は表示・記録装置である。
【０００５】
次に動作を説明する。送受信装置１０２で生成された電波は、空中線１０１によって空中に放射される。空中に放射された電波は、大気の屈折率の粗密が散乱体となって反射される。反射された電波は空中線によって受信され、送受信装置１０２に入力される。散乱体が上空の風と共に流れていると、受信された電波はドップラ効果によって周波数が変化する。この周波数の変化を一般のドップラレーダと同様に検出することにより、上空の風速を測定する。具体的には、受信された電波は送受信装置１０２によって、増幅、周波数変換、検波されることにより、ビデオ信号に変換され、レーダ信号処理装置１０３へと出力される。レーダ信号処理装置１０３では、ビデオ信号に周波数解析処理を行うことにより、ドップラ周波数を算出する。さらにドップラ周波数からドップラ速度を算出する。
【０００６】
図２１にレーダ信号処理装置１０３の内部構成を表すブロック図を示す。この図において、１はフーリエ変換部、２は周波数領域電力算出部、４はインコヒーレント積分部、５は信号検出部である。フーリエ変換部１では、各距離毎に得られる受信信号を、時間領域信号から周波数領域信号へと変換する。変換には、例えば高速フーリエ変換（Fast Fourier Transform; FFT)を用いれば良いが、その他の手法を用いても構わない。周波数領域信号に変換された信号は周波数領域電力算出部２に入力され、各周波数成分毎に電力を算出することにより、パワースペクトルが得られる。インコヒーレント積分部４では、パワースペクトルを積分することにより、パワースペクトルのゆらぎを低減させる。信号検出部５では、インコヒーレント積分後のパワースペクトルから大気エコーのスペクトルピークを検出し、その周波数からドップラ速度、すなわち視線方向風速などの大気の物理量を算出する。
【０００７】
品質管理処理装置１０４では、レーダ信号処理装置１０３から出力されたドップラ速度データを入力し、大気の時間的・空間的連続性を仮定することにより、時間軸及び空間軸方向で不連続なデータを精度の低いデータと判断し、欠損データとすることにより除去する。そして、除去されずに残ったデータを精度の高いデータとして出力する。
【０００８】
風速ベクトル算出装置１０５は、異なるビーム方向で観測されたドップラ速度データを入力し、それらを合成することにより、３次元風速ベクトルを算出する。図２２は風速ベクトル算出の原理を模式的に示したものである。図２２（ａ）の左側は上空の風速計測を行う際の典型的なビーム方向を示したものである。ここでは、天頂方向及び天頂角１０度程度で東西南北にビームを切り替えるとしている。図２２（ｂ）は図２２（ａ）の状況を東西方向に平行な垂直断面で切って表示したものである。西から東へと風が吹いている場合、西側のビームで観測すれば近づく方向のドップラ速度が、東側のビームで観測すれば遠ざかる方向のドップラ速度が観測される。西側のビームと東側のビームで観測されたドップラ速度の差異を利用することにより、東西方向の風速を求めることができる。同様にして、北側のビームと南側のビームを用いれば、南北方向の風速を得ることができる。さらに天頂方向のドップラ速度は風速の鉛直成分に相当する。その結果、３次元風速ベクトルが得られることになる。
【０００９】
表示・記録装置１０６では、風速ベクトル算出手段から出力された風速ベクトルデータを表示又は記録する。
【００１０】
従来の品質管理処理の方式として、「ウィンドプロファイラによる風観測の定常化に関する業務実験 − 調査結果と提言 −、測候時報 ６５．３、１９９８」で提案されている二次曲面近似方式がある。図２３は、この二次曲面近似方式を実現するための品質管理処理装置の構成を示すブロック図である。この図において、４０１は近似曲面定義範囲設定部、４０２は近似二次曲面算出部、４０３は近似精度判定部、４０４はデータ除去部である。
【００１１】
次に図２３の動作を説明する。近似曲面定義範囲設定部４０１は観測データに対して二次曲面をフィッティングする際の、近似二次曲面を定義する時間範囲と距離範囲を設定する。近似二次曲面算出部４０２では、入力データに対して、近似曲面定義範囲設定部４０１で設定された時間・距離範囲において、二次曲面をフィッティングする。二次曲面のフィッティングを行うには、例えば最小二乗法が用いられる。
【００１２】
近似精度判定部４０３では、近似二次曲面算出部４０２で算出された近似二次曲面の精度が高いかどうかを判定する。例えば前述の文献では、各時間・距離の格子点の観測値と近似二次曲面の差を算出し、差の最大値が予め設定されたしきい値よりも大きいか、又は近似曲面定義範囲内での差の二乗平均の平方根が予め設定されたしきい値よりも大きいかのいずれかが成立する場合に、近似二次曲面の精度が低いと判定される。
【００１３】
近似二次曲面の精度が低い場合、ドップラ速度データと近似二次曲面はデータ除去部４０４に入力される。データ除去部４０４では、近似二次曲面とドップラ速度データの差が大きい格子点において、ドップラ速度データを欠損とすることにより除去する。データ除去部４０４でデータ除去処理されたドップラ速度データは、近似二次曲面算出部４０２に入力され、二次曲面のフィッティングが再度行われる。このような処理を、近似二次曲面の精度が高いと近似精度判定部４０３で判定されるまで繰り返す。
【００１４】
近似二次曲面の精度が高いと判定されると、近似精度判定部４０３は、現時点でのドップラ速度データを品質管理処理済みのデータとして出力する。その後、近似曲面定義範囲設定部４０１は、近似曲面の定義範囲を次の高度又は次の時間へと更新して、同じ処理を繰り返す。
【００１５】
【発明が解決しようとする課題】
以上のような二次曲面近似方式による品質管理処理は、データ品質の劣化した格子点が少ない場合には良好に動作する。例えば、レーダの上空に一羽の鳥が存在するだけであれば、時空間座標の１点だけで鳥の移動速度に対応するドップラ速度が観測され、その他の点では風に対応するドップラ速度、すなわち視線方向風速が観測される。したがって、先の図２３の構成による品質管理装置を用いることにより、鳥エコーだけを除去することができる。
【００１６】
しかし、多数の鳥が時間的に連続して飛来するような場合、品質の劣化した格子点の数が多くなる。すると、時空間座標において、大気エコーよりも渡り鳥エコーが卓越して検出されるようになる。この場合、図２３のような従来の品質管理処理装置による異常エコーの除去は困難となる。このような状況は、渡り鳥が一斉に移動する季節になると頻繁に発生する。
【００１７】
レーダのビームは細いため、鳥がビーム内を通過する時間は限られている。したがって、時間的には断続的に鳥エコーが受信されることになる。鳥エコーが途切れた時間には大気エコーのみを受信することが可能である。しかし、大気に比べて鳥のレーダ断面積は大きいため、時間方向にデータを積分する積分処理を行うと、レーダ断面積の大きい鳥エコーのみが卓越し、大気エコーは埋もれてしまう。これにより、信号処理後のデータでは、鳥エコーは時間的に連続して現れることになる。
【００１８】
本発明は上記のような問題を解決するためになされたものであり、レーダによる大気観測を行う際に、多数の鳥が上空を通過することにより生じる異常エコーのみを除去し、大気エコーのデータを残すことにより、風速計測などの大気観測を正確に行うことができるようにすることを目的とする。
【００１９】
【課題を解決するための手段】
本発明の請求項１に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記異常エコー除去部は、上記パワースペクトルの電力値としてパワースペクトルを周波数軸上で積分した値を用い、この積分値がしきい値を越える場合に異常エコーであると判別し、異常エコーに相当しないパワースペクトルのみ出力することを特徴とするものである。
【００２０】
また、本発明の請求項２に係るレーダ信号処理装置は、上記パワースペクトルの電力値は、パワースペクトルを地形エコー成分が卓越する周波数範囲を除く周波数軸上で積分した値とすることを特徴とする請求項１に記載のものである。
【００２１】
また、本発明の請求項３に係るレーダ信号処理装置は、上記パワースペクトルの電力値は、パワースペクトルを周波数０の成分を除く周波数軸上で積分した値とすることを特徴とする請求項１に記載のものである。
【００２２】
また、本発明の請求項４に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記異常エコー除去部は、上記パワースペクトルの電力値としてパワースペクトルのピーク値を用い、このピーク値がしきい値を越える場合に異常エコーであると判別し、異常エコーに相当しないパワースペクトルのみ出力することを特徴とするものである。
【００２３】
また、本発明の請求項５に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、一定の時間範囲の受信信号の全電力を算出する時間領域電力算出部と、この全電力がしきい値を越える場合に当該時間範囲の受信信号を異常エコーであると判別し異常エコーに相当しない受信信号を出力する異常エコー除去部と、上記異常エコーに相当しない受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備えたことを特徴とするものである。
【００２４】
また、本発明の請求項６に係るレーダ信号処理装置は、日付又は時刻の条件に応じ、異常エコー判定を行うしきい値を設定するしきい値設定部を備え、この設定されたしきい値を用い異常エコー除去判定を行うことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２５】
また、本発明の請求項７に係るレーダ信号処理装置は、日付又は時刻の条件に応じ、異常エコー除去処理を行うか否かを切り替えることができるようにしたことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２６】
また、本発明の請求項８に係るレーダ信号処理装置は、上記日付又は時刻の条件として渡り鳥の飛来の多い季節でかつ夜間を条件としたことを特徴とする請求項６又は７に記載のものである。
【００２７】
また、本発明の請求項９に係るレーダ信号処理装置は、降雨強度に応じ、異常エコー判定を行うしきい値を設定するしきい値設定部を備え、この設定されたしきい値を用いて異常エコー除去判定を行うことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２８】
また、本発明の請求項１０に係るレーダ信号処理装置は、降雨強度に応じ、異常エコー除去処理を行うか否かを切り替えることができるようにしたことを特徴とする請求項１乃至５のいずれか一項に記載のものである。
【００２９】
また、本発明の請求項１１に係るレーダ信号処理装置は、上記信号検出部は鉛直方向の風速を算出し、この風速値に応じて異常エコー判定のしきい値を変化させるようにしたことを特徴とする請求項１乃至１０のいずれか一項に記載のものである。
【００３０】
また、本発明の請求項１２に係るレーダ信号処理装置は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記信号検出部が算出した鉛直方向の風速に基づき降雨の存在を判定し、降雨が存在する場合には異常エコー除去処理を行わないことを特徴とするものである。
【００３１】
また、本発明の請求項１３に係るレーダ信号処理装置は、上記信号検出部が算出した鉛直方向の風速に基づき降雨の存在を判定し、降雨が存在する場合には異常エコー除去処理を行わないことを特徴とする請求項１乃至１１のいずれか一項に記載のものである。
【００３２】
また、本発明の請求項１４に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルの電力を周波数軸上で積分することにより信号の全電力を算出する第３のステップと、各距離毎に全電力が予め設定したしきい値を超える場合に信号中に異常エコーが混入していると判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３３】
また、本発明の請求項１５に係るレーダ信号処理方法は、上記第３のステップは、各距離において、地形エコー成分である固定クラッタ領域を除く周波数範囲において、パワースペクトルの電力を周波数軸上で積分することにより、固定クラッタ電力を除く信号の全電力を算出するようにしたものである。
【００３４】
また、本発明の請求項１６に係るレーダ信号処理方法は、上記第３のステップは、周波数０の点のみを固定クラッタ成分とみなし、それ以外の周波数点のデータを積分することにより、固定クラッタ電力を除く信号の全電力を算出するものである。
【００３５】
また、本発明の請求項１７に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルのピーク電力を算出する第３のステップと、各距離毎にピーク電力が予め設定したしきい値を超える場合に信号中に異常エコーが混入していると判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３６】
また、本発明の請求項１８に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルの電力を周波数軸上で積分することにより信号の全電力を算出する第３のステップと、各距離で得られた全電力のうちの最大値が予め設定したしきい値を超える場合に、その時刻の全高度のパワースペクトルデータを使用しないと判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３７】
また、本発明の請求項１９に係るレーダ信号処理方法は、受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルのピーク電力値を算出する第３のステップと、各距離で得られたピーク電力値のうちの最大値が予め設定したしきい値を超える場合に、その時刻の全高度のパワースペクトルデータを使用しないと判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなるものである。
【００３８】
また、本発明の請求項２０に係るレーダ信号処理方法は、異常エコーの判定に用いるしきい値として、予め設定した固定値を用いるものである。
【００３９】
また、本発明の請求項２１に係るレーダ信号処理方法は、複数時刻で得られたパワースペクトルの電力値を比較して、相対的に高い電力値を持つ時刻以外の電力値を抽出し、抽出した電力値に予め定めた係数を乗算することにより異常エコーの判定に用いるしきい値を得るものである。
【００４０】
また、本発明の請求項２２に係るレーダ信号処理方法は、各距離において、一定時間内の複数の時間で得られたパワースペクトルから算出された全電力又はピーク電力のうち、最低の電力値を大気エコーの電力として抽出し、この抽出された最低電力値に予め設定した係数を乗算して得られる値をしきい値として設定するようにしたものである。
【００４１】
また、本発明の請求項２３に係るレーダ信号処理方法は、異常エコーを含まないデータが異常エコー判定しきい値を超える確率がある一定値以下となるように、最低電力値に乗算する係数を設定するものである。
【００４２】
また、本発明の請求項２４に係るレーダ信号処理方法は、異常エコーとして除去されたデータ数によらず、一定の積分時間でインコヒーレント積分処理を行うものである。
【００４３】
また、本発明の請求項２５に係るレーダ信号処理方法は、異常エコーとして除去されたデータ数によらず、インコヒーレント積分数が一定となるようにインコヒーレント積分処理を行うものである。
【００４４】
【発明の実施の形態】
実施の形態１．
図１はこの発明の実施の形態１のレーダ信号処理装置の構成を示すブロック図である。図において、１はフーリエ変換部、２は周波数領域電力算出部、３は異常エコー除去部、４はインコヒーレント積分部、５は信号検出部である。
【００４５】
次に動作について説明する。フーリエ変換部１では、各距離毎に得られる受信信号を、時間領域信号から周波数領域信号へと変換する。変換には、例えば高速フーリエ変換（Fast Fourier Transform; FFT)を用いれば良いが、その他の手法を用いても構わない。周波数領域信号に変換された信号は周波数領域電力算出部２に入力され、各周波数成分毎に電力を算出することにより、パワースペクトルが得られる。
【００４６】
異常エコー除去部３では、周波数領域電力算出部２で得られたパワースペクトルの電力値が、設定されたしきい値よりも大きい場合に、異常エコーが混入していると判定する。異常エコーが混入していると判定されたパワースペクトルは除去する。インコヒーレント積分部４では、異常エコー除去部３で除去されずに残ったパワースペクトルを積分することにより、パワースペクトルのゆらぎを低減させる。信号検出部５では、インコヒーレント積分後のパワースペクトルから大気エコーのスペクトルピークを検出し、その周波数からドップラ速度、すなわち視線方向風速を算出する。
【００４７】
なお、異常エコー除去部３におけるパワースペクトルの除去の方法としては、除去するパワースペクトルを欠損とする、すなわち「除去」のフラグを立てることが考えられる。インコヒーレント積分部４では、「除去」のフラグのついていないパワースペクトルのみを用いて積分処理を行うことになる。あるいは、異常エコー除去部３におけるパワースペクトルの除去の他の方法として、除去するパワースペクトルの電力値を全て０とすることも考えられる。この場合、インコヒーレント積分部４では、異常エコー除去部３でパワースペクトルが除去されたか否かを考慮せずに積分処理を行っても、異常エコーは結果的に除去できることになる。
【００４８】
異常エコー判定に用いるパワースペクトルの電力値としては、例えばパワースペクトルを周波数軸上で積分した全電力を用いることが考えられる。図２はそのような場合の異常エコー除去部３の内部構成を示すブロック図である。この図において、３０１は電力積分部、３０２は電力判定部である。電力積分部３０１では、パワースペクトルを周波数軸上で積分することにより、パワースペクトルの全電力を算出する。電力判定部３０２では、電力積分部３０１で算出された全電力を、予め設定したしきい値と比較し、全電力の方が大きければパワースペクトルに異常エコー成分が含まれると判定し、このパワースペクトルを除去する。逆に全電力の方が小さければ、パワースペクトルをそのままインコヒーレント積分部４へと出力する。
【００４９】
ただし、近距離のデータでは、地形エコー成分の電力が大きい場合がある。地形エコー成分は周波数０近辺のみに存在し、スペクトル幅も大気エコーに比べて明らかに狭いため、両者を識別することは比較的容易である。信号検出部５にて地形エコーの影響を取り除くことも通常良く行われることから、異常エコー除去部３においては地形エコー以外の異常エコー、すなわち鳥エコーのみを除去するようにしても良い。この場合、電力積分部３０１では、周波数０の近傍を除く周波数範囲でパワースペクトルを積分するようにすれば良い。
【００５０】
地形エコー成分を取り除く最も簡易な方法としては、周波数０成分のみを除いてパワースペクトルの積分値を求めることである。図３はそのような動作をする異常エコー除去部の構成を示す図である。この図において、３０３は周波数０成分減算部である。その他の符号は前述のものと同じである。周波数０成分減算部３０３では、電力積分部３０１において全周波数範囲でパワースペクトルを積分することにより得られた電力値から、周波数０成分を減算する。これにより、周波数０の成分を除く電力積分値を得ることができる。
【００５１】
以上では、パワースペクトルを積分することにより得られる電力を用いて異常エコー判定を行う例を示したが、図４では、ピーク電力値を用いて判定する場合の異常エコー除去部３の内部構成を示すブロック図である。この図において、３０４はピーク電力抽出部である。ピーク電力抽出部では、パワースペクトル中の最大の電力値を抽出し、この値を信号電力の代表値と見なして、３０２電力判定部にて異常エコー判定を行うようにしている。
【００５２】
図５は本発明の実施の形態のレーダ信号処理方法の流れを示す図である。この図では、パワースペクトルを周波数軸上で積分する場合、すなわち図２あるいは図３の異常エコー除去部の動作に対応したレーダ信号処理方法の流れを示している。
【００５３】
ステップｓｔ００１とステップｓｔ０１２で囲まれた区間は、インコヒーレント積分を行うループを表すものである。ステップ００２では、フーリエ変換を行う際の時間点数分の受信信号を入力する。ステップｓｔ００３からステップｓｔ００７で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。ｓｔ００４では、各距離において、時系列信号である受信信号にフーリエ変換を施すことにより、周波数領域の信号へと変換する。ｓｔ００５では、各周波数成分毎に電力値を算出することにより、パワースペクトルを得る。ｓｔ００６では、パワースペクトルを周波数軸上で積分することにより、全電力を算出する。周波数軸上の積分では、図２の異常エコー除去部３の動作のように、周波数０の近傍を除く周波数領域で積分するようにしても良いし、あるいは図３の異常エコー除去部３のように、周波数０の成分のみを除いた積分を行うようにしても良い。
【００５４】
ステップｓｔ００８からステップｓｔ０１１で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。ステップｓｔ００９では、各距離毎に、全電力と予め設定されたしきい値と比較し、全電力がしきい値を下回る場合にのみ、ステップｓｔ０１０にて処理中のパワースペクトルを積算する。インコヒーレント積分ループは、設定された回数Ｎinだけｓｔ００２からｓｔ０１１の処理を行うと終了するようにする。この場合、ｓｔ０１０で積算されるパワースペクトルの数、すなわちインコヒーレント積分数は最大でＮinとなり、異常エコーの影響を受ける場合には、これよりも小さいインコヒーレント積分数となる。インコヒーレント積分ループを、設定された回数Ｎinだけ行うようにした場合、インコヒーレント積分後のパワースペクトルの得られる時間間隔が、実際に実行したインコヒーレント積分数によらず一定となるため、品質管理処理装置によるデータ品質管理処理などの信号処理より後段のデータ処理を行う上で、入力データの扱いが容易になる。ただし、インコヒーレント積分数が少なくなった場合には、それだけ信号検出の性能が劣化することになる。
【００５５】
あるいは、インコヒーレント積分数が所定数以上になった場合にステップｓｔ０１２でインコヒーレント積分ループを終了するようにしても良い。この場合には、常にインコヒーレント積分数を一定以上にすることができるため、一定以上の信号検出性能を常に確保することができる。
【００５６】
なお、図５では、ステップｓｔ００１からステップｓｔ００７と、ステップｓｔ００８からステップｓｔ０１１の２つのレンジループを持っているが、これを一つにまとめるようにしても良い。あるいは処理が等価であればレンジループをさらに細かく分けるようにしても良い。このように、処理が等価であれば、本発明はループの区切り方を限定するものではない。
【００５７】
図５に示した本発明のレーダ信号処理方法の流れを示す図では、レーダのビーム方向の切り替えについては特に表現されていない。実際のレーダでは、例えば図６に示すように、１つの時間範囲でのフーリエ変換が終わると次のビーム方向に切り替え、必要な観測方向を全てカバーするというようにインコヒーレント積分を行う時間範囲の途中でビームを切り替えるような手順が取られることが多い。このような場合にも、個々のビーム方向のデータのみに着目すれば、図５の処理の流れは変わらない。よって、本発明は、このようなビーム切り替えの手順を特に限定するものではない。
【００５８】
異常エコー判定で用いるしきい値については、大気のレーダ断面積と鳥のレーダ断面積を考慮して固定的に設定すれば良い。一般に鳥のレーダ断面積は大気のレーダ断面積と比較して十分大きい。そこで、レーダ方程式から予測される大気エコーの受信電力と鳥エコーの受信電力の間の電力値をしきい値として採用すれば良い。
【００５９】
図７は、パワースペクトルのピーク電力を用いて異常エコー判定を行う場合、すなわち図４の異常エコー除去部の動作に対応したレーダ信号処理方法の流れを示す図である。この図では、ステップ００６ｂでパワースペクトルのピーク電力を抽出するようにし、そのピーク電力についてしきい値による判定をステップ００９ｂで行うようにしている以外は、図５に示したレーダ信号処理方法と同じである。
【００６０】
以上のように、本発明の実施の形態では、エコー強度によって異常エコーと大気エコーを識別し、大気エコーのみを用いて信号検出を行うようにしているため、渡り鳥が上空を飛来するような場合にも正確に大気観測を行うことが可能になる。特に、インコヒーレント積分前のデータに対して異常エコーの除去を行うようにしているため、従来のレーダ信号処理装置又は品質管理処理装置に比べて、異常エコー除去の処理を行う時間単位が短いため、断続的に異常エコーが発生している状況において、正常エコーが異常エコーに埋もれる前に異常エコーを除去し、正常エコーのみを残すことが可能となる。
【００６１】
また、パワースペクトルを算出するために算出した周波数領域上の電力値を用いるため、後述の実施の形態７と比較して時間領域での電力値算出が不要となるため、実施の形態７の場合よりも信号処理演算量は少なくて済む。
【００６２】
なお、本実施の形態では、風速計測用レーダへの適用を想定したレーダ信号処理装置及びレーダ信号処理方法について記述しているが、本レーダ信号処理装置及びレーダ信号処理方法は、電波を送受信するレーダだけに限定されるものではなく、音波を送受信するソーダや、光波を送受信するレーザレーダに対しても全く同様に適用することができる。また、風速計測を目的としたレーダだけでなく、降雨レーダを始めとする大気観測用レーダにも適用できる。このことは、以後に記述する別の実施の形態においても同様である。
【００６３】
実施の形態２．
以上の実施の形態１では、異常エコー判定及び異常エコー除去を距離毎に行うようにしたものであるが、次に、異常エコー判定及び異常エコー除去を全レンジで一括して行うような場合の実施の形態を示す。本発明の実施の形態において、レーダ信号処理装置の構成は、実施の形態１と同様に図１で表されるものであるが、レーダ信号処理方法が異なるものとなる。
【００６４】
図８は本発明の実施の形態２のレーダ信号処理方法の流れを示すものである。ステップｓｔ１０１とステップｓｔ１１３で囲まれた区間は、インコヒーレント積分を行うループを表すものである。ステップ１０２では、フーリエ変換を行う際の時間点数分の受信信号を入力する。ステップｓｔ１０３からステップｓｔ１０７で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。ｓｔ１０４では、各距離において、時系列信号である受信信号にフーリエ変換を施すことにより、周波数領域の信号へと変換する。ｓｔ１０５では、各周波数成分毎に電力値を算出することにより、パワースペクトルを算出する。ｓｔ１０６では、パワースペクトルを周波数軸上で積分することにより、全電力を算出する。周波数軸上の積分では、図２の異常エコー除去部３の動作のように、周波数０の近傍を除く周波数領域で積分するようにしても良いし、あるいは図３の異常エコー除去部３のように、周波数０の成分のみを除いた積分を行うようにしても良い。ステップｓｔ１０８では、ステップｓｔ１０６でレンジ毎に算出された全電力のうちの最大値を抽出する。ステップｓｔ１０９では、ステップｓｔ１０８でえられた最大全電力と予め設定したしきい値とを比較し、最大全電力の方が小さい場合にのみ、ステップｓｔ１１０からステップｓｔ１１２で表されるレンジループの処理を行う。レンジループ内のステップｓｔ１１１では、ステップｓｔ１０５で得られたパワースペクトルを積分後パワースペクトルに加算することにより、インコヒーレント積分を実行する。
【００６５】
このように図８のレーダ信号処理方法では、ステップｓｔ１０９のしきい値判定で異常エコーが受信信号に混入していると判定された場合には、全距離でパワースペクトルの積算が行われないため、全距離で一括して異常エコーが除去されることになる。
【００６６】
図８では、異常エコー判定を全電力に対して行うのに対して、図９のようにピーク電力に対して異常エコー判定を行うことも考えられる。すなわち図４の異常エコー除去部の動作に対応したレーダ信号処理方法の流れを示す図である。この図では、ステップ１０６ｂでパワースペクトルのピーク電力を抽出するようにし、そのピーク電力の最大値をステップ１０８ｂで抽出し、抽出された最大ピーク電力についてしきい値による判定をステップ１０９ｂで行うようにしている以外は、図８に示したレーダ信号処理方法と同じである。
【００６７】
本発明の実施の形態のレーダ信号処理方法によれば、全距離で一括して異常エコー判定を行うようにしているため、距離毎に異常エコー判定を行うよりも演算時間を短縮することができる。特に、以前から用いられているレーダ信号処理装置に異常エコー除去機能を加える場合のように、信号処理ハードウェアの性能に制約があるような場合に有効である。
【００６８】
実施の形態３．
以上の実施の形態では、異常エコー判定に用いるしきい値は固定したものを用いるとしていたが、本実施の形態では、しきい値を観測データから決定するような実施形態について述べる。
【００６９】
図１０は本発明の実施の形態におけるレーダ信号処理装置の構成を示すブロック図である。この図において、７はしきい値設定部である。その他の符号は前述のものと同じである。次に動作を説明する。しきい値設定部７は周波数領域電力算出部２で得られた各距離のパワースペクトルデータを入力し、複数時間分だけ蓄積する。各レンジ毎に、蓄積されたパワースペクトルの中から異常エコーが含まれないと判断されるデータを抽出し、その電力値をもとに異常エコー判定のしきい値を設定する。具体的には、電力値から判断して異常エコーが含まれないと判断される時刻のデータを選択し、その電力値に予め定義した方法でマージンを加えた電力値を、異常エコー判定しきい値として異常エコー除去部３へ出力する。このしきい値を用いて、異常エコー除去部３は周波数領域電力算出部２から出力されたパワースペクトルを評価し、異常エコー成分が含まれると判断されたデータを除去する。その他の動作は前述の図１のレーダ信号処理装置と同じである。
【００７０】
次に図１０のレーダ信号処理装置の動作を、処理アルゴリズムという観点から説明する。本実施の形態においても、しきい値を設定する方法以外の処理アルゴリズムについては、実施の形態１又は実施の形態２のレーダ信号処理方法と同じである。ここでは、レーダ信号処理方法のアルゴリズムのうちの、異常エコー判定に用いるしきい値を設定する部分のみを特に説明する。
【００７１】
図１１は本発明の実施の形態のレーダ信号処理方法における異常エコー判定しきい値の設定部分の流れを示す図である。この図において、ステップｓｔ２０１からステップｓｔ２０４で囲まれた区間は、レンジループ、すなわち各距離での処理を行うループを表すものである。
【００７２】
例えば実施の形態１で示した図５では、ステップｓｔ００３からステップｓｔ００７の処理において、レンジ毎のフーリエ変換及びパワースペクトル算出を複数の時間ブロックで行うことにより、各レンジにおいて複数時刻でのパワースペクトル及び全電力値が得られている。
【００７３】
図１１のｓｔ２０２では、複数時刻の全電力値から、最小の全電力値を抽出する。最小電力値となっている時間は、鳥による異常エコーが混入している可能性が最も低いと考えられる。そこで、最小全電力値が得られている時刻の観測データが異常エコーを含まないデータのサンプルであると見なして、ステップｓｔ２０３ではその全電力値にある係数を乗算した電力値を異常エコー判定のしきい値を設定する。乗算する係数の値は経験的に固定値を設定しても良い。あるいは大気エコーの電力値の統計的性質、すなわち確率密度分布を考慮して、異常エコーを含まない大気エコーがしきい値を超える確率、即ち誤って正常エコーが除去される確率が一定となるように設定しても良い。図１１における全電力としては、全周波数成分から計算した全電力を用いても良いし、あるいは、ドップラー周波数が０近傍である地形エコー成分を除いた全電力を用いても良い。
【００７４】
以上の方法によってしきい値設定を行えば、例えば大気エコーの受信強度が弱い場合には、異常エコーの影響を受けやすくなるが、異常エコー判定しきい値が低くなるため、より多くの異常エコー成分を除去することが可能となる。逆に大気エコーの受信強度が比較的強い場合には、しきい値が大きくなるため、大気エコーが強いために誤って異常エコーと判定される確率が低くなるため、多くの大気エコー成分を利用することが可能となる。
【００７５】
図１１では全電力値を利用して異常エコー判定しきい値を設定していたが、図１２のようにパワースペクトルのピーク値を用いて異常エコー判定しきい値を設定しても良い。図１２では、図１１のステップｓｔ２０２の代わりにステップｓｔ２０２ｂを用いることにより、複数時刻のピーク電力値から、最小のピーク電力値を抽出する。また、ステップｓｔ２０３の代わりにステップｓｔ２０３ｂを用いることにより、ピーク電力値にマージンを加算した電力レベルを異常エコー判定のしきい値を設定する。全電力の代わりにピーク電力を用いること以外は、図１１と図１２は同様の処理により異常エコー検出しきい値を設定する。
【００７６】
図１１又は図１２のしきい値設定処理は、レーダ信号処理の流れの中の適当なところにおけば良い。例えば図５に示したレーダ信号処理方法を用いる場合は、ｓｔ００７の直後に入れることが考えられる。あるいは、ｓｔ０１２の直後にしきい値設定処理を入れ、設定されたしきい値は次の時刻の信号処理に用いるようにしても良い。本発明では、しきい値処理を入れる段階を特に制限するものではない。
【００７７】
本実施の形態によれば、大気エコーのエコー強度に応じたしきい値を用いて異常エコー判定を行うようにしているため、誤って大気エコーを除去する確率を低くしたまま、異常エコーを正確に除去することが可能となり、観測精度を高く保つことができる。
【００７８】
実施の形態４．
前述の実施の形態３では、異常エコー判定しきい値を観測されたデータを用いて変更していたが、渡り鳥飛来の頻度は季節と時刻に依存することから、日時によって判定しきい値を変更するようにしても良い。具体的には、渡り鳥が移動する期間は毎年一定しており、しかも外敵から目立たない夜間に移動することが多い。そのため、渡り鳥が移動する期間、例えば春と秋の夜間に、異常エコー判定しきい値を低くし、異常エコーの残留が少なくなるようにすれば良い。その他の日時では、鳥が観測範囲内を通過するとしても、その頻度はまれであるため、従来から行われている信号処理後の品質管理処理にて除去できるため、異常エコー判定しきい値を高めにして、正常エコーを誤って除去する確率が低くなるようにすることが考えられる。
【００７９】
そのようなレーダ信号処理装置の構成例のブロック図を図１３に示す。この図において、８は時計である。その他の符号は前述のものと同じである。時計８から出力される日時のデータはしきい値設定部７へ渡される。しきい値設定部７では、日時毎に予め設定されたしきい値を異常エコー除去部３へ出力する。異常エコー除去部３はしきい値設定部７から入力したしきい値を用いて、周波数領域電力算出部２から入力したパワースペクトルに対して異常エコーか否かの判定を行い、異常エコーであると判定されれば、そのパワースペクトルデータを除去する。その他の動作については、図１０に示したレーダ信号処理装置と同じである。
【００８０】
図１３のレーダ信号処理装置の他の例を図１４に示す。この図において、９は処理切替スイッチである。その他の符号は前述のものと同じである。図１３では、日時によってしきい値を変化させていたが、図１４では異常エコー除去部３による処理を、日時によってオン／オフするようにしている。例えば、渡りの季節の夜間のみ処理切替スイッチをオンにして、異常エコー除去部３による処理を行い、それ以外の日時では処理切替スイッチをオフにして、異常エコー除去部３による処理を行わないようにする。
【００８１】
以上のように、日時によってしきい値を変更することにより、特に異常エコーが良く除去できるようにしたため、あるいは異常エコー除去処理をオン／オフすることにより渡り鳥による異常エコーの影響が大きく現れるときだけ異常エコーを除去するようにした。よって、渡り鳥の影響が小さい日時において、正常エコーが誤って除去されることがなくなり、インコヒーレント積分数が少なくならないため、正常エコーの検出確率が減らず、データ取得率を維持することが可能となる。
【００８２】
実施の形態５．
前述の実施の形態４では、日時によってレーダ信号処理方法を変更するようにしていたが、渡り鳥の渡りは天候にも左右される。本発明の実施の形態５では、天候によってレーダ信号処理方法を変化させるようにしたものである。図１５は本発明の実施の形態のレーダ信号処理装置の構成を示すブロック図である。この図の符号は全て前述のものと同じである。図１５では、しきい値設定部７が雨量データを入力し、降雨強度によって変化するしきい値を異常エコー除去部３へと出力する。雨量の強いときには、渡り鳥が飛来する可能性が低いため、しきい値を大きく設定することにより、正常エコーを誤って除去しないようにする。逆に降雨のないときには、渡り鳥が飛来する可能性があるため、しきい値を小さく設定することにより、異常エコーが除去できるようにする。
【００８３】
図１６に示すブロック図は図１５を変形したもので、雨量データによって、異常エコー除去部３による処理をオン／オフしている。すなわち、降雨のあるときには異常エコー除去部３の処理をオフにし、降雨のないときは逆にオンにするようにしている。
【００８４】
雨量データとしては、一般的な雨量計で観測されるものを用いてもよいし、あるいは風速計測用のレーダとは別に設置された降雨レーダで計測された雨量強度を用いても良いが、本発明では特に雨量データの観測方法を限定するものではない。
【００８５】
また、渡り鳥飛来の特性が雨量以外の気象データに依存する場合には、そのようなデータも合わせてしきい値設定部７に入力するようにしても良い。
【００８６】
本発明の実施の形態によれば、降雨の存在するときには渡り鳥が飛来しないと判断して信号処理を制御するため、降雨時に誤って正常エコーを除去することがなくなる。
【００８７】
実施の形態６．
前述の実施の形態５では、雨量データを用いて、降雨時の異常エコー除去処理を制御していたが、本発明の実施の形態では、レーダ装置で観測されたデータから降雨の有無を判定するような実施形態を説明する。
【００８８】
図１７及び図１８は本発明の実施の形態のレーダ信号処理装置の構成を示す図である。この図に含まれる符号は全て前述のものと同じである。次に動作を説明する。信号検出部から出力されるドップラ速度データのうち、ビームを天頂方向に向けて観測されたドップラ速度を処理切替スイッチ９へ入力する。降雨が存在しない場合、一般的に鉛直方向の風は１〜２ｍ／ｓ程度である。それに対し、雨滴の落下速度はそれよりも大きくなる。また、渡り鳥はほとんどの場合に水平飛行をすると考えられる。したがって、天頂ビームのドップラ速度が近づく方向に大きい場合には、降雨が存在し、渡り鳥は上空を通過していないと判断することができる。そこで、図１７のレーダ信号処理装置では、ドップラ速度が近づく方向に大きい場合に、異常エコー判定しきい値を大きくし、あるいは図１８のレーダ信号処理装置では、ドップラ速度が近づく方向に大きい場合に、異常エコー判定部３の処理をオフにする。
【００８９】
なお、鉛直方向の風速を得るのには、必ずしも天頂方向へビームを向けて観測しなくても良い。例えば、天頂角θで東方向と西方向を観測して得られるドップラ速度をVe とVw とすると、風速の鉛直成分ｗは下記式（１）で算出することができる。
ｗ＝(Ve + Vw)/２cosθ （１）
式（１）のような計算は信号検出部５で行うとしても良いし、あるいは図２０の風速ベクトル算出装置で算出される風速ベクトルの鉛直成分を利用するようにしても良い。
【００９０】
本発明の実施の形態によれば、別の装置で観測された雨量データを入力しなくても、降雨の有無を判定し、その結果で異常エコー除去部の動作を制御するようにしているため、他の装置からのデータ入力が不要となり、簡易な構成で適切な異常エコー除去処理が可能となる。
【００９１】
実施の形態７．
以上に述べた実施の形態では、時系列領域の受信信号を周波数領域に変換した後のデータに対して、異常エコーの判定を行っていた。ただし、異常エコーの有無の判定は電力値で行っているため、周波数領域に変換する前の時系列データに対して行うことも可能である。図１９は本発明の実施の形態のレーダ信号処理装置の構成を示すブロック図である。この図において、７０１は時間領域電力算出部である。
【００９２】
次に動作を説明する。時間領域電力算出部７０１では、時系列データである受信データを設定した時間範囲で区切り、その時間範囲における全電力を算出する。異常エコー除去部１では、時間領域電力算出部７０１から出力された電力を予め設定したしきい値と比較し、全電力の方がしきい値よりも大きければ、その時間範囲の時系列データを受信信号から除去する。除去する場合には、データを０としても良いし、あるいは欠損データであるというマークをつけておいても良い。それ以外の処理、すなわちフーリエ変換部２以降の処理については、実施の形態１などの前述の実施の形態の動作と同じである。
【００９３】
なお、時間領域電力算出部７０１の全電力算出における時間範囲の設定は、フーリエ変換部２の処理においてフーリエ変換を実行する際の時間窓と同じとしても良いし、あるいは全く別の時間範囲設定としても良い。同じ時間範囲設定とすれば、本実施の形態のように時間領域で全電力を算出しても、あるいは前述の実施の形態のように周波数領域で全電力を算出するようにしても、両者は同等となる。逆に時間範囲設定をフーリエ変換の時間窓と異なるように設定することも、本実施の形態では可能となるため、例えばフーリエ変換処理の時間長よりも短い時間で異常エコーが断続して生じるような場合にも、全電力を算出する時間範囲設定を短く設定することにより、異常エコーの存在する時間のデータだけを的確に除去することが可能であり、正常データの除去を最低限に抑えることが可能となる。
【００９４】
なお、全電力の代わりに、区切った時間範囲内の平均電力を時間領域全電力算出部７０１で算出しても良い。何故なら、全電力を区切った時間範囲の時間点数で割ったものであり、異常エコー判定のしきい値を時間点数で割ったものを用いれば、除去処理の原理としては実質的な差異がないからである。
【００９５】
【発明の効果】
本発明の請求項１によれば、受信信号の電力値として、周波数軸上で得られているパワースペクトルを積分した電力値を用い、インコヒーレント積分前の電力値をしきい値と比較することにより、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【００９６】
また、本発明の請求項２によれば、地形エコーが卓越する低周波領域を除いて全電力値を算出するようにしたため、強い地形エコーが混入している場合にも、渡り鳥による異常エコーの混入の有無の判定を正確に行うことが可能となる。
【００９７】
また、本発明の請求項３によれば、全ドップラ周波数範囲でパワースペクトルを積分して得られる全電力値から、周波数０の点の電力値のみを除くことにより、地形エコーを除いた全電力値を算出するようにしたため、簡易な処理で上記請求項３と同様の効果を得ることが可能となる。
【００９８】
また、本発明の請求項４によれば、受信信号の電力値として、周波数軸上で得られているパワースペクトルのピーク値を抽出したものを用いることにより、上記請求項１と同様の効果を得ることができる。
【００９９】
また、本発明の請求項５によれば、時系列データから算出した電力値をしきい値と比較することにより、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１００】
また、本発明の請求項６によれば、季節や時刻によって異常エコー判定のしきい値を調整することにより、渡り鳥の影響が小さい日時において、正常エコーが誤って除去されることがなくなり、インコヒーレント積分数が少なくならないため、正常エコーの検出確率を維持することができ、その結果としてデータ取得率を維持することが可能となる。
【０１０１】
また、本発明の請求項７によれば、季節や時刻によって異常エコー判定をオン／オフすることにより、渡り鳥の影響が小さい日時において、正常エコーが誤って除去されることがなくなり、インコヒーレント積分数が少なくならないため、正常エコーの検出確率を維持することができ、その結果として減らず、データ取得率を維持することが可能となる。
【０１０２】
また、本発明の請求項８によれば、渡り鳥の特性を考慮して、渡り鳥の飛来する季節の夜間に異常エコー処理をオンにし、それ以外の場合に異常エコー処理をオフとすることにより、上記請求項７と同様の効果を得ることが可能となる。
【０１０３】
また、本発明の請求項９によれば、降雨強度に応じて異常エコー判定のしきい値を調整するようにしたため、降雨時に誤って正常エコーを除去することがなくなるという効果がある。
【０１０４】
また、本発明の請求項１０によれば、降雨の存在するときには渡り鳥が飛来しないと判断してレーダ信号処理を制御するため、降雨時に誤って正常エコーを除去することがなくなるという効果がある。
【０１０５】
また、本発明の請求項１１によれば、レーダ装置で観測された雨滴の落下速度に応じて異常エコー判定のしきい値を調整するようにしたため、他の装置からのデータを入力することなく、降雨時に誤って正常エコーを除去することのないレーダ信号処理が可能となる。
【０１０６】
また、本発明の請求項１２によれば、レーダ装置の観測データから降雨の有無を判定し、その結果で異常エコー除去部の動作を制御するようにしているため、他の装置からのデータを入力することなく、降雨時に誤って正常エコーを除去することのないレーダ信号処理が可能となる。
【０１０７】
また、本発明の請求項１３によれば、レーダ装置の観測データから降雨の有無を判定し、その結果で異常エコー除去部の動作を制御するようにしているため、他の装置からのデータを入力することなく、降雨時に誤って正常エコーを除去することのないレーダ信号処理が可能となる。
【０１０８】
また、本発明の請求項１４によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１０９】
また、本発明の請求項１５によれば、強い地形エコーが混入している場合にも、渡り鳥による異常エコーの混入の有無の判定を正確に行い、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１０】
また、本発明の請求項１６によれば、全ドップラ周波数範囲でパワースペクトルを積分して得られる全電力値から、周波数０の点の電力値のみを除き、地形エコーを除いた全電力値を算出することにより、強い地形エコーが混入している場合にも、簡易な処理で渡り鳥による異常エコーの混入の有無の判定を正確に行い、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１１】
また、本発明の請求項１７によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１２】
また、本発明の請求項１８によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となり、かつ全距離で一括して異常エコー判定を行うようにしているため、距離毎に異常エコー判定を行うよりも演算時間を短縮することができる。
【０１１３】
また、本発明の請求項１９によれば、異常エコーと大気エコーを識別し、異常エコーを除去するようにしたため、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となり、かつ全距離で一括して異常エコー判定を行うようにしているため、距離毎に異常エコー判定を行うよりも演算時間を短縮することができる。
【０１１４】
また、本発明の請求項２０によれば、異常エコーの判定に用いるしきい値として、予め設定した固定値を用いることにより、断続的に生じる渡り鳥による異常エコーのみを除去し、大気エコーである正常エコーのみを用いて大気の物理量を求めるため、正確な測定が可能となる。
【０１１５】
また、本発明の請求項２１によれば、大気エコーのエコー強度に応じたしきい値を用いて異常エコー判定を行うようにしているため、誤って大気エコーを除去する確率を低くしたまま、異常エコーを正確に除去することが可能となり、観測精度を高く保つことができる。
【０１１６】
また、本発明の請求項２２によれば、大気エコーのエコー強度に応じたしきい値を用いて異常エコー判定を行うようにしているため、誤って大気エコーを除去する確率を低くしたまま、異常エコーを正確に除去することが可能となり、観測精度を高く保つことができる。
【０１１７】
また、本発明の請求項２３によれば、異常エコーを含まないデータが異常エコー判定しきい値を超える確率がある一定値以下となるように、最低電力値に乗算する係数を設定することにより、正常エコーの検出確率を一定に保ったまま、異常エコーを除去により観測精度を高く保つことができる。
【０１１８】
また、本発明の請求項２４によれば、異常エコーとして除去されたデータ数によらず、一定の積分時間でインコヒーレント積分処理を行うようにしたため、品質管理処理装置によるデータ品質管理処理などのレーダ信号処理より後段のデータ処理を行う上で、入力データの扱いが容易になる。
【０１１９】
また、本発明の請求項２５によれば、異常エコーとして除去されたデータ数によらず、インコヒーレント積分数が一定となるようにインコヒーレント積分処理を行うようにしたため、一定以上の信号検出性能を常に確保することができる。
【図面の簡単な説明】
【図１】 本発明の実施の形態１のレーダ信号処理装置の構成を示すブロック図。
【図２】 本発明の実施の形態１のレーダ信号処理装置の異常エコー除去部の構成を示すブロック図。
【図３】 本発明の実施の形態１のレーダ信号処理装置の異常エコー除去部の別の構成を示すブロック図。
【図４】 本発明の実施の形態１のレーダ信号処理装置の異常エコー除去部の別の構成を示すブロック図。
【図５】 本発明の実施の形態１のレーダ信号処理方法の流れを示す図。
【図６】 風速観測における典型的な観測手順を示す図。
【図７】 本発明の実施の形態１の別のレーダ信号処理方法の流れを示す図。
【図８】 本発明の実施の形態２のレーダ信号処理方法の流れを示す図。
【図９】 本発明の実施の形態２の別のレーダ信号処理方法の流れを示す図。
【図１０】 本発明の実施の形態３のレーダ信号処理装置の構成を示すブロック図。
【図１１】 本発明の実施の形態３におけるレーダ信号処理方法におけるしきい値設定アルゴリズムの流れを示す図。
【図１２】 本発明の実施の形態３におけるレーダ信号処理方法における別のしきい値設定アルゴリズムの流れを示す図。
【図１３】 本発明の実施の形態４のレーダ信号処理装置の構成を示すブロック図。
【図１４】 本発明の実施の形態４の別のレーダ信号処理装置の構成を示すブロック図。
【図１５】 本発明の実施の形態５のレーダ信号処理装置の構成を示すブロック図。
【図１６】 本発明の実施の形態５の別のレーダ信号処理装置の構成を示すブロック図。
【図１７】 本発明の実施の形態６のレーダ信号処理装置の構成を示すブロック図。
【図１８】 本発明の実施の形態６の別のレーダ信号処理装置の構成を示すブロック図。
【図１９】 本発明の実施の形態７のレーダ信号処理装置の構成を示すブロック図。
【図２０】 一般的なウィンドプロファイラの構成を示すブロック図
【図２１】 従来のウィンドプロファイラのレーダ信号処理装置の構成を示すブロック図。
【図２２】 ウィンドプロファイラにおける風速ベクトル算出の原理を示す図。
【図２３】 従来のウィンドプロファイラの品質管理処理装置の構成例を示すブロック図。
【符号の説明】
１ フーリエ変換部、 ２ 周波数領域電力算出部、 ３ 異常エコー除去部、 ４ インコヒーレント積分部、 ５ 信号検出部、 ７ しきい値設定部、８ 時計、 ９ 処理切替スイッチ、 １０１ 空中線、 １０２ 送受信装置、 １０３ レーダ信号処理装置、 １０４ 品質管理処理装置、 １０５ 風速ベクトル算出装置、 １０６ 表示記録装置、 ３０１ 電力積分部、 ３０２ 電力判定部、 ３０３ 周波数０成分減算部、 ３０４ ピーク電力抽出部、 ４０１ 近似曲面定義範囲部、 ４０２ 近似二次曲面算出部、 ４０３近似精度判定部、 ４０４ データ除去部、 ７０１ 時間領域電力算出部。
Claims (25)
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記異常エコー除去部は、上記パワースペクトルの電力値としてパワースペクトルを周波数軸上で積分した値を用い、この積分値がしきい値を越える場合に異常エコーであると判別し、異常エコーに相当しないパワースペクトルのみ出力することを特徴とするレーダ信号処理装置。
- 上記パワースペクトルの電力値は、パワースペクトルを地形エコー成分が卓越する周波数範囲を除く周波数軸上で積分した値とすることを特徴とする請求項１に記載のレーダ信号処理装置。
- 上記パワースペクトルの電力値は、パワースペクトルを周波数０の成分を除く周波数軸上で積分した値とすることを特徴とする請求項１に記載のレーダ信号処理装置。
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記異常エコー除去部は、上記パワースペクトルの電力値としてパワースペクトルのピーク値を用い、このピーク値がしきい値を越える場合に異常エコーであると判別し、異常エコーに相当しないパワースペクトルのみ出力することを特徴とするレーダ信号処理装置。
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、一定の時間範囲の受信信号の全電力を算出する時間領域電力算出部と、この全電力がしきい値を越える場合に当該時間範囲の受信信号を異常エコーであると判別し異常エコーに相当しない受信信号を出力する異常エコー除去部と、上記異常エコーに相当しない受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備えたことを特徴とするレーダ信号処理装置。
- 日付又は時刻の条件に応じ、異常エコー判定を行うしきい値を設定するしきい値設定部を備え、この設定されたしきい値を用い異常エコー除去判定を行うことを特徴とする請求項１乃至５のいずれか一項に記載のレーダ信号処理装置。
- 日付又は時刻の条件に応じ、異常エコー除去処理を行うか否かを切り替えることができるようにしたことを特徴とする請求項１乃至５のいずれか一項に記載のレーダ信号処理装置。
- 上記日付又は時刻の条件として渡り鳥の飛来の多い季節でかつ夜間を条件としたことを特徴とする請求項６又は７に記載のレーダ信号処理装置。
- 降雨強度に応じ、異常エコー判定を行うしきい値を設定するしきい値設定部を備え、この設定されたしきい値を用いて異常エコー除去判定を行うことを特徴とする請求項１乃至５のいずれか一項に記載のレーダ信号処理装置。
- 降雨強度に応じ、異常エコー除去処理を行うか否かを切り替えるこ とができるようにしたことを特徴とする請求項１乃至５のいずれか一項に記載のレーダ信号処理装置。
- 上記信号検出部は鉛直方向の風速を算出し、この風速値に応じて異常エコー判定のしきい値を変化させるようにしたことを特徴とする請求項１乃至１０のいずれか一項に記載のレーダ信号処理装置。
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理装置において、受信信号を周波数信号にフーリエ変換するフーリエ変換部と、この周波数信号から周波数毎の電力であるパワースペクトルを算出する周波数領域電力算出部と、このパワースペクトルの電力値に基づき上記異常エコーを判別し異常エコーに相当しないパワースペクトルのみ出力する異常エコー除去部と、この異常エコーに相当しないパワースペクトルのみをインコヒーレント積分し平均化するインコヒーレント処理部と、このインコヒーレント積分したパワースペクトルから大気の物理量を算出する信号検出部とを備え、上記信号検出部が算出した鉛直方向の風速に基づき降雨の存在を判定し、降雨が存在する場合には異常エコー除去処理を行わないことを特徴とするレーダ信号処理装置。
- 上記信号検出部が算出した鉛直方向の風速に基づき降雨の存在を判定し、降雨が存在する場合には異常エコー除去処理を行わないことを特徴とする請求項１乃至１１のいずれか一項に記載のレーダ信号処理装置。
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルの電力を周波数軸上で積分することにより信号の全電力を算出する第３のステップと、各距離毎に全電力が予め設定したしきい値を超える場合に信号中に異常エコーが混入していると判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなることを特徴とするレーダ信号処理方法。
- 上記第３のステップは、各距離において、地形エコー成分である固定クラッタ領域を除く周波数範囲において、パワースペクトルの電力を周波数軸上で積分することにより、固定クラッタ電力を除く信号の全電力を算出するようにしたことを特徴とする請求項１４に記載のレーダ信号処理方法。
- 上記第３のステップは、周波数０の点のみを固定クラッタ成分とみなし、それ以外の周波数点のデータを積分することにより、固定クラッタ電力を除く信号の全電力を算出することを特徴とする請求項１４に記載のレーダ信号処理方法。
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルのピーク電力を算出する第３のステップと、各距離毎にピーク電力が予め設定したしきい値を超える場合に信号中に異常エコーが混入していると判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなることを特徴とするレーダ信号処理方法。
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルの電力を周波数軸上で積分することにより信号の全電力を算出する第３のステップと、各距離で得られた全電力のうちの最大値が予め設定したしきい値を超える場合に、その時刻の全高度のパワースペクトルデータを使用しないと判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなることを特徴とするレーダ信号処理方法。
- 受信信号から異常エコーを除去し大気の物理量を測定するレーダ装置に適用するレーダ信号処理方法において、各距離で得られる受信信号からある時間範囲の信号を切り出し、フーリエ変換を施すことにより周波数領域のデータに変換する第１のステップと、上記第１のステップで算出された周波数領域のデータについて周波数毎に電力値を算出してパワースペクトルを求める第２のステップと、各距離においてパワースペクトルのピーク電力値を算出する第３のステップと、各距離で得られたピーク電力値のうちの最大値が予め設定したしきい値を超える場合に、その時刻の全高度のパワースペクトルデータを使用しないと判定する第４のステップと、上記第１のステップから第４のステップの処理を複数の時間で行うことにより得られる複数の時間でのパワースペクトルのうち、異常エコーが混入していないと判定された場合のパワースペクトルのみを用いてインコヒーレント積分を行う第５のステップと、上記第５のステップで得られたインコヒーレント積分後のパワースペクトルから大気エコーを検出し、大気の物理量を算出する第６のステップからなることを特徴とするレーダ信号処理方法。
- 異常エコーの判定に用いるしきい値として、予め設定した固定値を用いることを特徴とする請求項１４乃至１９のいずれか一項に記載のレーダ信号処理方法。
- 複数時刻で得られたパワースペクトルの電力値を比較して、相対的に高い電力値を持つ時刻以外の電力値を抽出し、抽出した電力値に予め定めた係数を乗算することにより異常エコーの判定に用いるしきい値を得ることを特徴とする請求項１４乃至１９のいずれか一項に記載のレーダ信号処理方法。
- 各距離において、一定時間内の複数の時間で得られたパワースペクトルから算出された全電力又はピーク電力のうち、最低の電力値を大気エコーの電力として抽出し、この抽出された最低電力値に予め設定した係数を乗算して得られる値をしきい値として設定するようにしたことを特徴とする請求項２１に記載のレーダ信号処理方法。
- 異常エコーを含まないデータが異常エコー判定しきい値を超える確率がある一定値以下となるように、最低電力値に乗算する係数を設定することを特徴とする請求項２２に記載のレーダ信号処理方法。
- 異常エコーとして除去されたデータ数によらず、一定の積分時間でインコヒーレント積分処理を行うことを特徴とする請求項１４乃至２３のいずれか一項に記載のレーダ信号処理方法。
- 異常エコーとして除去されたデータ数によらず、インコヒーレント積分数が一定となるようにインコヒーレント積分処理を行うことを特徴とする請求項１４乃至２３のいずれか一項に記載のレーダ信号処理方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002249336A JP3794361B2 (ja) | 2002-08-28 | 2002-08-28 | レーダ信号処理装置及びレーダ信号処理方法 |
| US10/396,517 US6894638B2 (en) | 2002-08-28 | 2003-03-26 | Radar signal processing unit and radar signal processing method for abnormal signal extraction |
| FR0304891A FR2844057B1 (fr) | 2002-08-28 | 2003-04-18 | Unite de traitement de signaux de radar et procede de traite ment de signaux de radar |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002249336A JP3794361B2 (ja) | 2002-08-28 | 2002-08-28 | レーダ信号処理装置及びレーダ信号処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004085472A JP2004085472A (ja) | 2004-03-18 |
| JP3794361B2 true JP3794361B2 (ja) | 2006-07-05 |
Family
ID=31712250
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002249336A Expired - Lifetime JP3794361B2 (ja) | 2002-08-28 | 2002-08-28 | レーダ信号処理装置及びレーダ信号処理方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6894638B2 (ja) |
| JP (1) | JP3794361B2 (ja) |
| FR (1) | FR2844057B1 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8830113B2 (en) | 2010-11-18 | 2014-09-09 | Furuno Electric Company Limited | Method and program for detecting object of interest, and radar device |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0411097D0 (en) * | 2004-05-19 | 2004-06-23 | Qinetiq Ltd | Laser radar device and method |
| JP3858098B2 (ja) * | 2004-05-26 | 2006-12-13 | 防衛庁技術研究本部長 | 目標検出装置 |
| WO2006000835A1 (en) * | 2004-06-24 | 2006-01-05 | Bae Systems Integrated System Technologies Limited | Improvements relating to velocity extraction |
| US20070171122A1 (en) * | 2006-01-25 | 2007-07-26 | Fujitsu Ten Limited | Radar apparatus and interference detection method |
| JP5111778B2 (ja) * | 2006-04-28 | 2013-01-09 | 株式会社東芝 | レーダ装置とその信号処理方法 |
| JP2008014837A (ja) * | 2006-07-06 | 2008-01-24 | Toshiba Corp | レーダ装置とその信号処理方法 |
| JP4903094B2 (ja) * | 2007-07-19 | 2012-03-21 | 三菱電機株式会社 | レーダ装置 |
| DE102009000472A1 (de) * | 2009-01-29 | 2010-08-05 | Robert Bosch Gmbh | Verfahren zur Detektion von Niederschlag mit einem Radarortungsgerät für Kraftfahrzeuge |
| JP5325361B2 (ja) * | 2009-02-06 | 2013-10-23 | 古野電気株式会社 | レーダ装置 |
| US8264908B2 (en) * | 2009-03-09 | 2012-09-11 | Second Wind Systems, Inc. | Method of detecting and compensating for precipitation in sodar systems |
| JP2011002425A (ja) * | 2009-06-22 | 2011-01-06 | Furuno Electric Co Ltd | レーダ装置 |
| JP6150474B2 (ja) * | 2012-08-09 | 2017-06-21 | 三菱電機株式会社 | ウィンドプロファイラにおける信号処理装置および信号処理方法 |
| JP5875492B2 (ja) * | 2012-09-13 | 2016-03-02 | 三菱電機株式会社 | 風計測装置 |
| JP6167644B2 (ja) * | 2013-04-26 | 2017-07-26 | 住友電気工業株式会社 | 気象観測装置、気象観測プログラム及び気象観測方法 |
| EP3056924B1 (en) * | 2013-10-11 | 2020-01-01 | Mitsubishi Electric Corporation | Wind measurement lidar device |
| CN103983975B (zh) * | 2014-06-05 | 2017-01-11 | 中国气象局气象探测中心 | 基于两种雷达的大气运动垂直速度检测方法及系统 |
| JP6435040B2 (ja) * | 2015-04-16 | 2018-12-05 | 古野電気株式会社 | レーダシステム |
| JP6053866B2 (ja) * | 2015-05-20 | 2016-12-27 | 株式会社日立パワーソリューションズ | 風力発電システムおよび鳥飛来予測装置 |
| US10048370B2 (en) * | 2015-07-07 | 2018-08-14 | Honeywell International Inc. | Estimating weather and ground reflectivity with doppler spectral information |
| JP6562583B2 (ja) * | 2016-05-31 | 2019-08-21 | メトロウェザー株式会社 | ドップラーシフト解析装置 |
| KR102176230B1 (ko) * | 2017-07-14 | 2020-11-09 | 주식회사 엘지화학 | 고분자 막의 분석 방법 |
| CN109254291B (zh) * | 2018-10-22 | 2021-03-23 | 中国气象科学研究院 | 云雷达的多普勒功率谱数据融合方法及装置 |
| CN112526547B (zh) * | 2020-11-18 | 2024-02-23 | 南京泰爱信科技有限公司 | 一种基于测风激光雷达的大气边界层分类方法及装置 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62100676A (ja) | 1985-10-28 | 1987-05-11 | Toshiba Corp | レ−ダ情報処理装置 |
| JPS62161070A (ja) | 1986-01-09 | 1987-07-17 | Toshiba Corp | レ−ダ情報処理装置 |
| JPS62161071A (ja) | 1986-01-09 | 1987-07-17 | Toshiba Corp | レ−ダ情報処理装置 |
| US5262782A (en) * | 1991-04-09 | 1993-11-16 | Unisys Corporation | Microburst precursor detector |
| US5157403A (en) * | 1991-10-04 | 1992-10-20 | General Electric Co. | Noise correction in radars for diffuse targets |
| US5175551A (en) * | 1991-12-18 | 1992-12-29 | Unisys Corporation | Downdraft velocity estimator for a microburst precursor detection system |
| US5563604A (en) * | 1995-02-08 | 1996-10-08 | Alliedsignal Inc. | Weather radar using spectral gaussian envelope discrimination for clutter rejection |
| US5568151A (en) * | 1995-05-17 | 1996-10-22 | Merritt; David A. | Statistical averaging method for wind profiler doppler spectra |
| US5592171A (en) * | 1995-08-17 | 1997-01-07 | The United States Of America As Represented By The Secretary Of Commerce | Wind profiling radar |
| US6563452B1 (en) * | 1998-07-06 | 2003-05-13 | Honeywell International Inc. | Apparatus and method for determining wind profiles and for predicting clear air turbulence |
| US6307500B1 (en) * | 1999-08-13 | 2001-10-23 | University Corporation For Atmospheric Research | Method of moment estimation and feature extraction for devices which measure spectra as a function of range or time |
| US6480142B1 (en) * | 2001-05-17 | 2002-11-12 | William L. Rubin | Method and apparatus for measuring velocity and turbulence of atmospheric flows |
-
2002
- 2002-08-28 JP JP2002249336A patent/JP3794361B2/ja not_active Expired - Lifetime
-
2003
- 2003-03-26 US US10/396,517 patent/US6894638B2/en not_active Expired - Lifetime
- 2003-04-18 FR FR0304891A patent/FR2844057B1/fr not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8830113B2 (en) | 2010-11-18 | 2014-09-09 | Furuno Electric Company Limited | Method and program for detecting object of interest, and radar device |
Also Published As
| Publication number | Publication date |
|---|---|
| FR2844057A1 (fr) | 2004-03-05 |
| US20040041725A1 (en) | 2004-03-04 |
| JP2004085472A (ja) | 2004-03-18 |
| US6894638B2 (en) | 2005-05-17 |
| FR2844057B1 (fr) | 2004-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3794361B2 (ja) | レーダ信号処理装置及びレーダ信号処理方法 | |
| US8818029B1 (en) | Weather forecasting systems and methods | |
| KR101431707B1 (ko) | 엑스밴드 이중편파 레이더 관측자료를 이용한 통합형 강우량 산정 방법 | |
| US5568151A (en) | Statistical averaging method for wind profiler doppler spectra | |
| Goudenhoofdt et al. | Generation and verification of rainfall estimates from 10-yr volumetric weather radar measurements | |
| CN104508515B (zh) | 用于天气检测的毫米波雷达系统和天气检测方法 | |
| EP3951439B1 (en) | Advection fog forecasting system and forecasting method | |
| Zhuge et al. | Summertime convective initiation nowcasting over southeastern China based on Advanced Himawari Imager observations | |
| Chang et al. | An application of reflectivity climatology in constructing radar hybrid scans over complex terrain | |
| CN110749871B (zh) | 双偏振天气雷达的参量估算方法 | |
| EP2990820B1 (en) | Clutter suppressing device and radar apparatus provided with the same | |
| JP2023047319A (ja) | 気象レーダーのためのノイズ除去装置および方法 | |
| KR101263121B1 (ko) | 위성 적외영상 자료를 이용한 현업용 기상레이더 반사도 합성자료의 채프에코 제거 방법 | |
| Koistinen et al. | Operational measurement of precipitation in cold climates | |
| Sugier et al. | Detection and removal of clutter and anaprop in radar data using a statistical scheme based on echo fluctuation | |
| Cifelli et al. | Dual-polarization radar rainfall estimation | |
| US5689444A (en) | Statistical quality control of wind profiler data | |
| Hartung et al. | An intercomparison of UW cloud-top cooling rates with WSR-88D radar data | |
| Ruan et al. | Analysis of the microphysical properties of a stratiform rain event using an L-Band profiler radar | |
| Li et al. | A surface wind model–based method to estimate rain-induced radar path attenuation over ocean | |
| CN113987094B (zh) | 一种基于气象雷达的gis地图预警方法 | |
| CN115932859A (zh) | 基于毫米波云雷达电网覆冰水粒子相态识别方法及系统 | |
| Capozzi et al. | X-band weather radar monitoring of precipitation fields in Naples urban areas: data quality, comparison and analysis | |
| JP3702347B2 (ja) | ウィンドプロファイラにおける信号処理方法及びそのプログラムと装置 | |
| JP2007003396A (ja) | レーダ信号処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20031211 |
|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20040712 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050802 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20051129 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060124 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060322 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060404 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 3794361 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100421 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100421 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110421 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110421 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120421 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120421 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130421 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130421 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140421 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |