JP3626089B2 - ウィンドプロファイラにおける信号処理装置および信号処理方法 - Google Patents
ウィンドプロファイラにおける信号処理装置および信号処理方法 Download PDFInfo
- Publication number
- JP3626089B2 JP3626089B2 JP2000365808A JP2000365808A JP3626089B2 JP 3626089 B2 JP3626089 B2 JP 3626089B2 JP 2000365808 A JP2000365808 A JP 2000365808A JP 2000365808 A JP2000365808 A JP 2000365808A JP 3626089 B2 JP3626089 B2 JP 3626089B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- altitude
- time
- integration time
- incoherent integration
- 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 - Fee Related
Links
Images
Classifications
-
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Description
【０００１】
【発明の属する技術分野】
この発明は、上空の風速分布を測定するウィンドプロファイラにおける信号処理装置および信号処理方法に関するものである。
【０００２】
【従来の技術】
風向・風速の情報は気象予報に必要な情報の一つである。風向・風速を計測する最も一般的な方法としては、地上に風向・風速計を設置することである。しかしながら、その場合には地表付近の風しか測定できない。気象予報をより正確にするためには上空の風向・風速のデータも知ることが必要である。そのため、従来、ゾンデ等を用いて上空の風向・風速の観測を行っていた。このゾンデを用いた観測では、ゾンデを上げた時刻における風向・風速のデータしか得ることができないため、観測の時間分解能が数時間以上と低いことが欠点であった。
【０００３】
それに対して、近年、ウィンドプロファイラと呼ばれる大気レーダを用いて上空の風向・風速を計測する技術が確立されつつある。ウィンドプロファイラでは、１分〜数分毎に上空の風向・風速を計測することが可能となる。このような高い時間分解能による上空の風向・風速のデータは、気象予報の精度向上に有効となると期待されている。このようなウィンドプロファイラについては、例えばＣａｒｔｅｒ ｅｔ ａｌ．， Ｄｅｖｅｌｏｐｍｅｎｔ ｉｎ ＵＨＦ ｌｏｗｅｒ ｔｒｏｐｏｓｈｅｒｉｃ ｗｉｎｄ ｐｒｏｆｉｌｉｎｇ ａｔ ＮＯＡＡ’ｓ Ａｅｒｏｎｏｍｙ Ｌａｂｏｒａｔｏｒｙ， Ｒａｄｉｏ Ｓｃｉｅｎｃｅ， ｖｏｌ．３０， ｎｏ．４， ｐｐ．９７〜１００１， １９９５．に紹介されたものなどがある。
【０００４】
ここで、ウィンドプロファイラで大気を計測する原理を説明する。
ウィンドプロファイラはドップラレーダの一種であり、一般的に図１３に示すように、空中線１０１、送受信装置１０２、信号処理装置１０３、風速ベクトル算出装置１０４、表示・記録装置１０５によって構成されている。このように構成されたウィンドプロファイラでは、空中線１０１より空中に放射された電磁波の反射波が空中線１０１にて受信され、その受信された電磁波が送受信装置１０２によって増幅、周波数変換、検波されてビデオ信号に変換される。このビデオ信号に信号処理装置１０３で周波数解析処理を行うことにより、ドップラ周波数を算出し、そのドップラ周波数からドップラ速度を算出して風速ベクトル算出装置１０４に送り、そこで算出された風速ベクトルを表示・記録装置１０５に表示もしくは記録する。
【０００５】
なお、上記従来のウィンドプロファイラでは信号処理装置１０３として、例えば図１４に示すような構成のものが一般に用いられている。この図１４において、１はＦＦＴ処理部、２はインコヒーレント積分部、３はドップラ速度算出部、４は品質管理・時間平均処理部である。
【０００６】
次にその動作について説明する。
送受信装置１０２より出力された受信信号のデータは、ＦＦＴ処理部１に入力される。ＦＦＴ処理部１では受信信号のデータに対してフーリエ変換を施す。このフーリエ変換後の受信信号のデータはインコヒーレント積分部２に入力される。ここでは、受信信号のデータのフーリエ変換から電力値を求めることによって得られるパワースペクトルを算出した後、複数の時刻に得られたパワースペクトルを積算する処理、すなわちインコヒーレント積分の処理を行う。受信信号のデータのパワースペクトルには、大気乱流エコーのスペクトル、すなわち信号スペクトルと雑音スペクトルが含まれる。このインコヒーレント積分処理により、雑音成分のゆらぎおよび信号成分のゆらぎの幅が小さくなる。そのため、インコヒーレント積分を行うことで信号検出の精度をあげることができる。ただし、インコヒーレント積分は複数の時刻のデータを用いるため、時間分解能を低下させることになる。
【０００７】
ドップラ速度算出部３では、インコヒーレント積分後の安定化されたパワースペクトルから信号スペクトルのピークを検出し、その中心ドップラ周波数から平均ドップラ周波数を算出し、さらに、そのドップラ周波数をドップラ速度に変換することにより、ターゲットである大気乱流の視線方向速度を算出する。
【０００８】
ただし、ドップラレーダで直接計測される風速は、実際の風速の視線方向への射影成分のみとなる。そこで、風速ベクトル算出装置１０４では、レーダ上空のある一定領域の風速分布が一様であると仮定して、レーダの観測方向を複数方向に変えて測定を行うことにより、３次元の風速ベクトルの合成を行う。
【０００９】
なお、このようなウィンドプロファイラでは、インコヒーレント積分時間が高度によらず一定である。しかし、ウィンドプロファイラのターゲットである大気乱流のレーダ反射率は、高度の上昇とともに急激に減少するため、実際には高度によって最適な積分時間が異なる。すなわち、ＳＮ比の高い低高度よりも、ＳＮ比が低下する高高度の方が長い積分時間を必要とする。
【００１０】
また、高高度の風速は地表面の影響をあまり受けないため、風速の時間変化が低高度の場合よりも小さい。このことからも、高高度ではインコヒーレント積分時間を長くすることが有利となる。
【００１１】
このように、ウィンドプロファイラのインコヒーレント積分時間は、高度によって最適な値が異なるものである。
【発明が解決しようとする課題】
【００１２】
従来のウィンドプロファイラにおける信号処理装置は以上のように構成されているので、インコヒーレント積分時間が高度によらず一定であり、高高度ではインコヒーレント積分時間が不足するため、十分な信号検出性能が得られず、データ取得率が低下するという課題があった。
【００１３】
この発明は上記のような課題を解決するためになされたもので、高度に応じてインコヒーレント積分時間を設定することにより、できるだけ広い高度範囲でデータ取得率を向上させた、ウィンドプロファイラにおける信号処理装置および信号処理方法を得ることを目的とする。
【００１４】
【課題を解決するための手段】
この発明に係るウィンドプロファイラにおける信号処理装置は、受信信号のデータにフーリエ変換を施すフーリエ変換処理部と、高度毎に最適なインコヒーレント積分時間を設定する積分時間設定部と、フーリエ変換されたデータからパワースペクトルを算出し、設定されたインコヒーレント積分時間だけそのパワースペクトルを時間積分するインコヒーレント積分部と、パワースペクトルから信号スペクトルのピークを検出し、そのドップラ周波数から大気の視線方向速度を算出するドップラ速度算出部と、低品質のドップラ速度のデータを取り除くとともに時間平均を行う品質管理・時間平均処理部とを備え、風速の測定を行うようにしたものである。
【００１５】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分時間を、低高度においては短く、高高度においては長く設定するようにしたものである。
【００１６】
この発明に係るウィンドプロファイラにおける信号処理装置は、積分時間設定部が、ドップラ速度算出部の出力するドップラ速度データよりデータ取得率を算出するデータ取得率算出部からのデータ取得率に応じて、インコヒーレント積分時間を設定するようにしたものである。
【００１７】
この発明に係るウィンドプロファイラにおける信号処理装置は、積分時間設定部が、品質管理・時間平均処理部の出力するドップラ速度データよりデータ取得率を算出するデータ取得率算出部からのデータ取得率に応じて、インコヒーレント積分時間を設定するようにしたものである。
【００１８】
この発明に係るウィンドプロファイラにおける信号処理装置は、各高度のインコヒーレント積分時間を、積分時間設定部が、データ取得率算出部にて算出されたその高度のデータ取得率をもとに設定するようにしたものである。
【００１９】
この発明に係るウィンドプロファイラにおける信号処理装置は、あらかじめ分割された高度範囲毎のインコヒーレント積分時間を、積分時間設定部が、データ取得率算出部にて算出されたあらかじめ分割された高度範囲毎のデータ取得率をもとに設定するようにしたものである。
【００２０】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分時間を長く設定した高度については、データ補間処理部で時間方向にデータを補間して、高度によらず一定の時間分解能のドップラ速度を品質管理・時間平均処理部に出力するようにしたものである。
【００２１】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分時間に反比例するデータ数を用いて、品質管理・時間平均処理部による時間平均処理を行い、データの時間分解能が高度によらず一定となるドップラ速度を出力するようにしたものである。
【００２２】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分の時間方向の積分範囲の更新幅をインコヒーレント積分時間によらず一定とし、インコヒーレント積分部からの積分後のデータの時間分解能を高度によらず一定となるようにしたものである。
【００２３】
この発明に係るウィンドプロファイラにおける信号処理方法は、受信信号のデータを読込んでフーリエ変換を施し、高度毎にインコヒーレント積分時間を設定し、フーリエ変換されたデータからパワースペクトルを算出して、設定されたインコヒーレント積分時間だけインコヒーレント積分を行い、そのインコヒーレント積分後のパワースペクトルから信号スペクトルを検出して、ドップラ周波数からドップラ速度を算出し、処理対象となっている時刻に関する全高度のデータの処理が未終了であれば、処理対象の高度を次の高度へ更新した後、インコヒーレント積分時間の設定に処理を戻し、終了していれば、処理対象の時刻を次の時刻に進めるとともに処理対象となる高度を初期化し、処理をデータの読込みに戻すことによって、風速の測定を行うようにしたものである。
【００２４】
この発明に係るウィンドプロファイラにおける信号処理方法は、受信信号のデータを読込んでフーリエ変換を施し、インコヒーレント積分時間の初期値を設定し、フーリエ変換されたデータからパワースペクトルを算出して、設定されたインコヒーレント積分時間だけインコヒーレント積分を行い、そのインコヒーレント積分後のパワースペクトルから信号スペクトルを検出して、ドップラ周波数からドップラ速度を算出し、そのドップラ速度から処理対象となっている高度のデータ取得率の算出を行って、そのデータ取得率があらかじめ定められたしきい値未満であり、かつインコヒーレント積分時間がその上限未満であれば、インコヒーレント積分時間の設定を増加させて、インコヒーレント積分の処理に戻し、それ以外の場合には、処理対象となっている時刻に関する全高度のデータの処理が終了していることを確認し、未終了であれば、処理対象の高度を次の高度へ更新した後、インコヒーレント積分時間の設定に処理を戻し、終了していれば、処理対象の時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、処理をデータの読込みに戻すことによって、風速の測定を行うようにしたものである。
【００２５】
【発明の実施の形態】
以下、この発明の実施の一形態を説明する。
実施の形態１．
図１はこの発明の実施の形態１によるウィンドプロファイラにおける信号処理装置を示すブロック図である。なお、ウィンドプロファイラの全体構成は図１３に示した従来の場合と同様であり、この図１に示した信号処理装置１０３は、図１３に示されるウィンドプロファイラの一部分を形成している。
【００２６】
図１３において、１０１は電磁波の送受信を行う空中線であり、１０２は空中線１０１で送受信される電磁波の入出力処理を行う送受信装置である。１０３は図１にその構成を示すこの発明の実施の形態１による信号処理装置である。１０４は信号処理装置１０３の出力より３次元の風速ベクトルの合成を行う風速ベクトル算出装置であり、１０５はこの風速ベクトル算出装置１０４の計算結果の表示あるいは記録を行う表示・記録装置である。
【００２７】
また、図１において、１は送受信装置１０２で受信された受信信号のデータにフーリエ変換を施すフーリエ変換処理部であり、ここでは、例えば高速フーリエ変換（ＦＦＴ）を施すＦＦＴ処理部が用いられている。２はＦＦＴ処理部１にてフーリエ変換されて出力された受信信号のデータからパワースペクトルを算出するとともに、複数の時刻で得られたパワースペクトルをあつめて積分するインコヒーレント積分部である。３はインコヒーレント積分部２から出力されたパワースペクトルから信号スペクトルを検出し、その中心周波数からターゲットのドップラ速度、すなわちターゲットの視線方向速度を算出するドップラ速度算出部である。４はドップラ速度算出部３で算出されたドップラ速度データのうちから品質の低下したデータを除去した後、その時間平均を行って風速ベクトル算出装置１０４に出力する品質管理・時間平均処理部である。５は高度に応じてインコヒーレント積分時間を設定し、インコヒーレント積分部２に積分時間を指示する積分時間設定部である。
【００２８】
次に動作について説明する。
図１３に示すように構成されたウィンドプロファイラでは、送受信装置１０２で生成した電磁波を空中線１０１より空中に放射する。空中に放射された電磁波は大気の屈折率の粗密が散乱体となって反射される。反射された電磁波は空中線１０１によって受信され、送受信装置１０２に入力される。散乱体が上空の風と共に流れていると、受信された電磁波はドップラ効果によって周波数が変化する。この周波数の変化を一般のドップラレーダと同様に検出することにより、上空の風速を測定する。具体的には、送受信装置１０２によって、受信された電磁波の増幅、周波数変換、検波を行うことによりビデオ信号に変換し、それを信号処理装置１０３へ出力する。信号処理装置１０３では、そのビデオ信号に周波数解析処理を行うことにより、ドップラ周波数を算出する。さらにドップラ周波数からドップラ速度を算出する。
【００２９】
以下に、実施の形態１による信号処理装置１０３の詳細な処理動作について説明する。
図１３に示す送受信装置１０２から出力された受信信号のデータは、図１のＦＦＴ処理部１に入力されてフーリエ変換が施される。ＦＦＴ処理部１にてフーリエ変換された受信信号のデータはインコヒーレント積分部２に入力される。インコヒーレント積分部２では、周波数毎に振幅の絶対値の二乗を計算することによってパワースペクトルを算出する。インコヒーレント積分部２はさらに、複数の時刻のデータから算出されたパワースペクトルを集めてそれを積分する。なお、積分処理において積算する時間方向のデータ数は、積分時間設定部５において高度毎に設定された積分時間によって定める。
【００３０】
インコヒーレント積分後のパワースペクトルはインコヒーレント積分部２よりドップラ速度算出部３に入力される。ドップラ速度算出部３では、インコヒーレント積分後の安定化されたパワースペクトルから信号スペクトルのピークを検出して、その中心ドップラ周波数から平均ドップラ周波数を算出し、さらに、そのドップラ周波数をドップラ速度に変換することにより、ターゲットである大気乱流の視線方向速度を算出する。
【００３１】
ただし、ドップラレーダで直接計測される風速は、実際の風速の視線方向への射影成分のみとなる。そこで、レーダ上空のある一定領域の風速分布が一様であると仮定して、レーダの観測方向を複数方向に変えて測定を行い、風速ベクトル算出装置１０４でそれを３次元の風速ベクトルに合成する。
【００３２】
この３次元風速ベクトルの算出原理を図２を用いて説明する。図２（ａ）はウィンドプロファイラにおけるビーム方向の例である。この例では天頂、および東、西、南、北の計５つの方向にビームを向けて観測している。天頂方向以外のビームの天頂角は、例えば１０度程度に設定される。なお、３次元風速ベクトルは東西方向成分、南北方向成分、鉛直方向成分に分けることができる。鉛直方向成分は天頂方向に指向したビーム（以下天頂ビームという）で得られるドップラ速度そのものとなる。東西方向成分は東方向に指向したビーム（以下東ビームといい、他の方角についても同様とする）のドップラ速度と西ビームのドップラ速度を用いて求めることができる。
【００３３】
図２（ｂ）は、東西方向に水平に風が吹いている場合に、東ビームと西ビームでどのようなドップラ速度が観測されるかを示したものである。図示のように、西から東に風が吹いている場合、西ビームでは近づく方向のドップラ速度が観測され、東ビームでは遠ざかる方向のドップラ速度が観測される。このように、異なるビーム方向で同一の風を観測すると、異なるドップラ速度が得られる。従って、ビーム間のドップラ速度の差異から水平風成分を算出することが可能となる。同様にして、北ビームと南ビームを用いれば、風速の南北方向成分を算出することができる。
【００３４】
なお、図３はウィンドプロファイラで観測される受信信号のデータのパワースペクトルを模式的に示した説明図である。ＳＮ比が高い場合には、信号スペクトルのピークを容易に検出することができるが、ＳＮ比が小さくなると、信号スペクトルのピーク値が低くなり、雑音のゆらぎと識別しずらくなる。インコヒーレント積分部２による積分は、雑音および信号の統計的なゆらぎを小さくする効果がある。そのため、積分時間設定部５によってインコヒーレント積分時間を長くすれば、図３に示した雑音電力密度の標準偏差が小さくなり、ＳＮ比が小さく、信号スペクトルのピーク値が低い場合にも、信号スペクトルのピークを検出することが可能となる。
【００３５】
ドップラ速度算出部３はこのようにしてパワースペクトルから大気乱流エコーのピークを検出し、その中心ドップラ周波数から、大気乱流エコーのドップラ速度、すなわち視線方向風速を求める。ドップラ速度算出部３から出力された視線方向風速は品質管理・時間平均処理部４において低品質なデータを除去した後、時間平均処理が施されて風速ベクトル算出装置１０４へと出力される。風速ベクトル算出装置１０４では、異なるビーム方向で観測したデータを合成することによって３次元風速ベクトルを算出する。この風速ベクトル算出装置１０４の計算結果は、表示・記録装置１０５に送られて表示あるいは記録される。
【００３６】
積分時間設定部５では、高度毎に最適となるインコヒーレント積分時間を設定する。低高度の範囲ではレーダとターゲットの距離が近いこと、および大気濃度が濃いため、大気乱流による電磁波の反射率が高い。そのため、比較的高いＳＮ比の受信信号が観測される。よって、インコヒーレント積分時間が短くても、パワースペクトルから信号スペクトルを容易に検出することができる。また、低高度の風は地表面の影響を受け易く、風向・風速の時間変化が短時間で生じる。そのため、３次元風速ベクトルの測定は高い時間分解能で行うことが望ましい。インコヒーレント積分後に時間平均処理を行う場合においても、地表面の影響を受けてある時刻のみに特異な風が生じた場合においても、高い時間分解能で風速ベクトルを求めておけば、特異な風速を品質管理処理によって除去し、卓越する傾向の風速ベクトルデータのみを用いて平均風速を求めることができる。
【００３７】
一方、高高度の範囲ではレーダとターゲットの距離が大きいこと、および大気濃度が薄くなって大気乱流の電磁波反射率が急激に減少することにより、受信信号のＳＮ比は低くなることが多い。従って、信号スペクトルを正確に検出するためには、インコヒーレント積分数を多くしてパワースペクトルのゆらぎを小さくする必要がある。また、高高度では、風が地表面の影響をあまり受けないことから、風速ベクトルの時間変化は小さい。そのことから、風速ベクトル算出の時間分解能は低高度の場合より低くてもあまり問題とならない。
【００３８】
以上のことから、インコヒーレント積分時間を、低高度で短く、高高度で長くすることが、一般的には適切な信号処理方法となると考えられる。典型的な例を図４に示す。この図では、インコヒーレント積分時間を、低高度である高度範囲１では１分、中程度の高度である高度範囲２では５分、高高度である高度範囲３では１０分とすることを示している。
【００３９】
以下、このような信号処理装置１０３の動作について、図５に示すフローチャートを用いて詳細に説明する。
まずデータ読込みステップＳＴ１において、送受信装置１０２からの受信信号のデータを信号処理装置１０３に入力する。信号処理装置１０３ではその入力された受信信号のデータに対して、フーリエ変換ステップ（以下ＦＦＴステップという）ＳＴ２において高速フーリエ変換を施す。次に積分時間設定ステップＳＴ３に進み、高度毎のインコヒーレント積分時間を設定する。次にインコヒーレント積分ステップＳＴ４において、ＦＦＴステップＳＴ２で得られた受信信号のデータのフーリエ変換からパワースペクトルを算出し、積分時間設定ステップＳＴ３で設定されたインコヒーレント積分時間だけ、インコヒーレント積分を実施する。
【００４０】
次にドップラ速度算出ステップＳＴ５において、インコヒーレント積分ステップＳＴ４で算出されたインコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出する。その後、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを、全高度処理終了判定ステップＳＴ６において判定する。判定の結果、全高度についての処理が済んでいなければ、高度更新ステップＳＴ７に進んで、処理対象となる高度を次の高度に更新した後、処理を積分時間設定ステップＳＴ３に戻す。一方、全高度についての処理が済んでいれば、時刻更新ステップＳＴ８において処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化して、処理をデータ読込みステップＳＴ１に戻す。
【００４１】
なお、品質管理・時間平均処理部４では品質管理処理および時間平均処理を実施する。この品質管理処理は品質が劣化したと判断されるデータを除去するための処理であり、ある時間・空間範囲でのデータの連続性を仮定することにより、不連続なデータを除去して品質の高いデータのみを残している。
【００４２】
この品質管理・時間平均処理部４で不連続なデータを除去する際、従来の品質管理処理では、データが全ての高度で同じ時間分解能で得られるということを前提にしているか、もしくは全ての高度で同じ時間分解能であることが望ましい処理方式となっていることがある。そのような場合、品質管理・時間平均処理部４において、入力されたドップラ速度データに補間処理を行うことも有効であると考えられる。
【００４３】
このようなデータ補間処理が行われる、この発明の実施の形態１におけるウィンドプロファイラの信号処理装置の構成例を図６に示す。図において、６はドップラ速度算出部３より入力されたドップラ速度データに補間処理を施し、品質管理・時間平均処理部４に出力するデータ補間処理部である。
【００４４】
このデータ補間処理部６では、全高度範囲のうち、最も時間分解能の高い高度における時間分解能と同一の時間分解能となるように、他の高度において時間方向にデータ補間を行う。例えば、インコヒーレント積分時間を長くしたことによって、品質管理処理を行う時間範囲に時間方向に一つのデータしか存在しない場合には、全補間点で同一の風速値を仮定すればよい。
【００４５】
以上のように、この実施の形態１によれば、高度毎にインコヒーレント積分時間の設定を行っているので、高度毎のターゲットの性質に応じて、適切な時間分解能で、かつ検出精度の高い信号処理を実現することができるという効果が得られる。
【００４６】
また、インコヒーレント積分時間を、低高度では短く、高高度では長く設定しているので、地表面に近い低高度では高い時間分解能で風速ベクトルを測定するとともに、ＳＮ比が低下する高高度でも検出精度を高めることが可能になるという効果も得られる。
【００４７】
さらに、データ補間処理を組合せることにより、高度毎にインコヒーレント積分時間が異なる場合にも、出力されるドップラ速度の時間分解能は高度によらず一定となるという効果も得られる。
【００４８】
なお、全高度において同じ時間間隔でドップラ速度データ、あるいは風速データが得られれば、数値気象予測の初期値として風速データを利用することをはじめとする風速データの２次利用において、データを利用しやすい形で出力することが可能になるという効果が得られる。
【００４９】
実施の形態２．
ここで、大気の状態は天候によって変化するため、最適なインコヒーレント積分時間は時間とともに変化する。上記実施の形態１では、高度毎にインコヒーレント積分時間をあらかじめ定めるものについて説明したが、大気の状態によって適応的に高度毎にインコヒーレント積分時間を設定するようにしてもよい。この実施の形態２は、そのような大気の状態によって適応的に高度毎にインコヒーレント積分時間を設定する、ウィンドプロファイラにおける信号処理装置および信号処理方法に関するものである。
【００５０】
図７はこの発明の実施の形態２による、そのようなウィンドプロファイラの信号処理装置を示すブロック図である。図において、１０３は信号処理装置、１はＦＦＴ処理部、２はインコヒーレント積分部、３はドップラ速度算出部、４は品質管理・時間平均処理部、５は積分時間設定部であり、これらは図１に同一符号を付したものと同等のものであるため、詳細な説明は省略する。また、７はドップラ速度算出部３の出力するドップラ速度データを入力してデータ取得率を算出し、それを積分時間設定部５に入力するデータ取得率算出部である。なお、積分時間設定部５は、データ取得率算出部７で算出されたデータ取得率をもとにインコヒーレント積分時間を設定している点で、図１に同一符号を付したものとは異なっている。
【００５１】
次に動作について説明する。
ここで、この図７に示す構成の信号処理装置１０３において、ＦＦＴ処理部１による受信信号のデータのフーリエ変換、インコヒーレント積分部２によるインコヒーレント積分、ドップラ速度算出部３による信号スペクトルの検出およびドップラ速度算出、品質管理・時間平均処理部４によるドップラ速度の時間平均の各処理については、図１に示した実施の形態１の場合と同様の動作となる。なお、この実施の形態２において特徴的であるのは、データ取得率算出部７の算出したデータ取得率をもとに、積分時間設定部５がインコヒーレント積分時間を設定することにある。
【００５２】
ドップラ速度算出部３の算出した各高度、各ビームでのドップラ速度は、品質管理・時間平均処理部４とともに、データ取得率算出部７にも入力される。このドップラ速度算出部３で得られたドップラ速度は、ＳＮ比の低下などにより信号スペクトルが検出されなかった場合には、その高度およびビームのドップラ速度は欠損となる。データ取得率算出部７では、このドップラ速度算出部３からのドップラ速度データより、全データに対する欠損でないデータの比率をデータ取得率として算出する。なお、このデータ取得率は各高度、各ビーム毎に、複数の時刻のデータから算出することができる。
【００５３】
また、一つの時刻のデータからデータ取得率を算出するために、観測高度範囲をいくつかの高度範囲に分割し、その高度範囲内のデータからデータ取得率を算出するようにしてもよい。その場合、データ取得率算出部７であらかじめ分割された高度範囲毎のデータ取得率を算出し、その高度範囲毎のデータ取得率をもとに、積分時間設定部５にてあらかじめ分割された高度範囲毎のインコヒーレント積分時間の設定を行う。
【００５４】
データ取得率算出部７で算出されたデータ取得率は、積分時間設定部５に入力される。積分時間設定部５では、あらかじめ設定したしきい値よりも取得率が低くなる高度においてインコヒーレント積分時間を増加させ、標準値よりも大きくなるようにする。データ取得率に対するしきい値は全高度で同じでもよいし、あるいは高度によって異なるしきい値を用いてもよい。
【００５５】
以下、このような信号処理装置１０３の動作について、図８に示すフローチャートを用いて詳細に説明する。
まずデータ読込みステップＳＴ１０にて、送受信装置１０２の出力する受信信号のデータを信号処理装置１０３に入力し、ＦＦＴステップＳＴ１１においてその入力された受信信号のデータに対して高速フーリエ変換を施す。次に積分時間初期値設定ステップＳＴ１２に進み、インコヒーレント積分時間の初期値の設定を行う。次にインコヒーレント積分ステップＳＴ１３にて、ＦＦＴステップＳＴ１１で得られたフーリエ変換された受信信号のデータからパワースペクトルを算出し、積分時間初期値設定ステップＳＴ１２において設定されたインコヒーレント積分時間だけ、インコヒーレント積分を実施する。
【００５６】
次にドップラ速度算出ステップＳＴ１４において、上記インコヒーレント積分ステップＳＴ１３で算出されたインコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出する。次にデータ取得率算出ステップＳＴ１５に進み、上記ドップラ速度算出ステップＳＴ１４にて算出されたドップラ速度から、現在処理対象となっている高度のデータ取得率を算出する。その後、データ取得率判定ステップＳＴ１６において、上記データ取得率算出ステップＳＴ１５にて算出されたデータ取得率をあらかじめ定めたしきい値と比較する。
【００５７】
比較の結果、ドップラ速度のデータ取得率がしきい値を下回る場合には積分時間判定ステップＳＴ１７へ進み、しきい値以上のデータ取得率でドップラ速度が得られている場合には、後述する全高度処理終了判定ステップＳＴ１９に進む。積分時間判定ステップＳＴ１７では、現時点のインコヒーレント積分時間があらかじめ定められたインコヒーレント積分時間（上限）を超えているか否かを識別し、現時点のインコヒーレント積分時間が上限を超えている場合には後述する全高度処理終了判定ステップＳＴ１９に進み、それ以外の場合には積分時間増加ステップＳＴ１８に進む。この積分時間増加ステップＳＴ１８ではインコヒーレント積分時間の設定を増加させた後、インコヒーレント積分ステップＳＴ１３に処理を戻す。
【００５８】
また、上記全高度処理終了判定ステップＳＴ１９では、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かについて判定する。判定の結果、全高度についての処理が済んでいなければ、高度更新ステップＳＴ２０に進んで、処理対象となる高度を次の高度へと更新し、処理を積分時間初期値設定ステップＳＴ１２に戻す。一方、全高度についての処理が済んでいれば、時刻更新ステップＳＴ２１において処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化して、処理をデータ読込みステップＳＴ１０に戻す。
【００５９】
以上の処理ステップについて、データ数や積分時間をもう少し具体的に記述することにより補足する。
データ読込みステップＳＴ１０において入力される受信信号のデータは、インコヒーレント積分時間の上限と等しい時間分のデータ量を有している。ここで、このインコヒーレント積分時間の上限をＴｍａｘとする。また、ＦＦＴステップＳＴ１１にて、各高度において、時間方向にＦＦＴ点数で区切られたデータに対してフーリエ変換を施す。次に積分時間初期値設定ステップＳＴ１２にてインコヒーレント積分時間の最小値となるインコヒーレント積分時間の初期値を設定する。
【００６０】
次にインコヒーレント積分ステップＳＴ１３では、フーリエ変換されたデータをインコヒーレント積分時間毎にまとめ、インコヒーレント積分を行う。今、インコヒーレント積分時間がＴｉであるとすると、１つの高度に対してＴｍａｘ／Ｔｉ個のインコヒーレント積分後のパワースペクトルが得られる。それぞれのパワースペクトルから信号スペクトルピークを検出し、ドップラ速度算出ステップＳＴ１４にてそのドップラ速度を算出する。この時点で最大Ｔｍａｘ／Ｔｉ個のドップラ速度が算出される。ただし、ＳＮ比が低い場合などには、信号スペクトルが検出できない場合があるため、実際に得られるドップラ速度はＴｍａｘ／Ｔｉ個よりも小さい。
【００６１】
そこで、データ取得率算出ステップＳＴ１５において、現在注目している高度のデータ取得率を算出する。次にデータ取得率判定ステップＳＴ１６でそのデータ取得率をしきい値と比較し、データ取得率があらかじめ設定したしきい値よりも大きければ、全高度処理終了判定ステップＳＴ１９に進む。逆にデータ取得率がしきい値よりも小さければ、積分時間判定ステップＳＴ１７において、現時点のインコヒーレント積分時間が上限を超えていないことを確認する。上限を超えていないことが確認されれば、積分時間増加ステップＳＴ１８にてインコヒーレント積分時間を増加させて、再度インコヒーレント積分ステップＳＴ１３に戻ってインコヒーレント積分からの処理を繰り返す。以降、このインコヒーレント積分ステップＳＴ１３から積分時間増加ステップＳＴ１８までの処理を、処理中の高度のデータ取得率がしきい値を超えるか、あるいはインコヒーレント積分時間が上限を超えるまで繰り返す。
【００６２】
以上、１高度毎にデータ取得率を算出する場合についての処理フローを説明したが、高度範囲毎にデータ取得率を算出する場合には、ＦＦＴステップＳＴ１１から積分時間増加ステップＳＴ１８までの処理を１つの高度範囲毎に実施し、高度更新ステップＳＴ２０で高度範囲の単位で高度の更新を行うことになる。この場合、高度範囲毎にインコヒーレント積分時間の修正処理を行っているため、１高度毎にインコヒーレント積分時間を設定するのに比べて処理が簡略化され、信号処理の負荷を小さくすることができる。
【００６３】
なお、上記実施の形態２の説明では、データ取得率の算出を、ドップラ速度算出部３から出力されたドップラ速度データを用いて行っているものについて示したが、品質管理・時間平均処理部４から出力された時間平均後のドップラ速度をデータ取得率算出部７に入力して、データ取得率を算出するようにしてもよい。図９のブロック図にはそのような構成によるウィンドプロファイラにおける信号処理装置を示す。この場合、データ取得率算出部７へのデータ入力が、ドップラ速度算出部３からではなく、品質管理・時間平均処理部４から行われているだけで、その他は図７と同一である。
【００６４】
また、この実施の形態２においても、上記実施の形態１の場合と同様に、ドップラ速度算出部３と品質管理・時間平均処理部４との間にデータ補間処理部６を設け、ドップラ速度算出部３より入力されたドップラ速度データに補間処理を施して、品質管理・時間平均処理部４に高度によらず一定の時間分解能のドップラ速度を出力するようにしてもよい。
【００６５】
以上のように、この実施の形態２によれば、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することが可能になるため、信号処理の性能を常に最適な状態に保つことができ、また、上空風の特性が観測地によって異なっても、この実施の形態２では自動的にインコヒーレント積分時間を設定しているため、設置時の信号処理設定が容易になるなどの効果が得られる。
【００６６】
実施の形態３．
なお、上記各実施の形態では、高度毎にインコヒーレント積分時間が異なるため、最終的に得られる風速データの時間間隔は高度によってまちまちになることがあるが、気象予測のためなどにデータを利用する際は、全高度において同じ時間間隔で風速データが得られる方が便利なことが多いと考えられる。そのため、最終的に全ての高度で同じ時間間隔で風速データが出力されるように、インコヒーレント積分後のデータの時間分解能の差を、時間平均を行う際の時間平均数の差によって吸収することも考えられる。この実施の形態３は、そのようなウィンドプロファイラにおける信号処理装置に関するものである。
【００６７】
図１０はこの発明の実施の形態３による、そのようなウィンドプロファイラの信号処理装置を示すブロック図であり、実施の形態２の各部に相当する部分には図７と同一符号を付してその説明を省略する。図において、８は積分時間設定部５から出力される各高度におけるインコヒーレント積分時間を受けて、そのインコヒーレント積分時間に反比例するデータ数を各高度の時間平均数として設定し、品質管理・時間平均処理部４に各高度の時間平均数を指示する高度毎時間平均数設定部である。
【００６８】
次に動作について説明する。
ここで、この図１０に示した信号処理装置１０３における、ＦＦＴ処理部１による受信信号のデータのフーリエ変換、インコヒーレント積分部２によるインコヒーレント積分、ドップラ速度算出部３による信号スペクトルの検出およびドップラ速度算出、データ取得率算出部７によるデータ取得率の算出、および積分時間設定部５による高度毎のインコヒーレント積分時間の設定の各処理については、図７に示した実施の形態２と同一の動作となる。この実施の形態３で特徴的なのは、データを平均する時間が高度によらず一定となるように、インコヒーレント積分時間に反比例する数のデータを用いて時間平均処理を行うようにしたことにある。
【００６９】
データ取得率算出部７ではドップラ速度算出部３からのドップラ速度データよりデータ取得率を算出し、積分時間設定部５に入力する。積分時間設定部５はこのデータ取得率算出部７の算出したデータ取得率をもとに、各高度範囲のインコヒーレント積分時間を設定する。この積分時間設定部５から出力される各高度におけるインコヒーレント積分時間はインコヒーレント積分部２とともに、高度毎時間平均数設定部８にも入力される。高度毎時間平均数設定部８では入力されたインコヒーレント積分時間に反比例するデータ数を各高度の時間平均数に設定する。品質管理・時間平均処理部４はこの時間平均数を用いて時間平均処理を行う。その結果、品質管理・時間平均処理部４から出力されるドップラ速度データは高度によらず同じ時間間隔で出力される。
【００７０】
なお、上記実施の形態３の説明では、ドップラ速度算出部３から出力されたドップラ速度データを用いてデータ取得率を算出しているものを示したが、品質管理・時間平均処理部４から出力された時間平均後のドップラ速度をデータ取得率算出部７に入力して、データ取得率を算出するようにしてもよい。図１１のブロック図にはそのような構成によるウィンドプロファイラにおける信号処理装置を示す。この場合、データ取得率算出部７へのデータ入力を、ドップラ速度算出部３に代えて品質管理・時間平均処理部４から行っているだけで、その他は図１０と同一である。
【００７１】
また、この実施の形態３においても、上記各実施の形態の場合と同様に、ドップラ速度算出部３と品質管理・時間平均処理部４との間にデータ補間処理部６を設け、ドップラ速度算出部３より入力されたドップラ速度データに補間処理を施して、品質管理・時間平均処理部４に高度によらず一定の時間分解能のドップラ速度を出力するようにしてもよい。
【００７２】
以上のように、この実施の形態３によれば、高度毎時間平均数設定部８より品質管理・時間平均処理部４に、各高度の時間平均数を指示しているので、全高度において同じ時間間隔でドップラ速度データあるいは風速データを得ることができるため、数値気象予測の初期値として風速データを利用することをはじめとする風速データの２次利用において、データを利用しやすい形で出力することが可能になるという効果が得られる。
【００７３】
実施の形態４．
また、上記実施の形態３では、ドップラ速度データが高度によらず同じ時間間隔で出力されるものについて説明したが、この発明はそれにのみ限定されるものではない。以下、図１２を用いてそのようなこの発明の実施の形態４について説明する。
【００７４】
図１２はこの発明の実施の形態４におけるインコヒーレント積分の時間範囲の更新方法を示す説明図である。この図１２では、低高度でインコヒーレント積分時間が短く、高高度でインコヒーレント積分が長く設定されている状況を想定している。高高度ではインコヒーレント積分の時間範囲を更新する際に、前後の時間帯のインコヒーレント積分において、インコヒーレント積分時間範囲が重なるようにし、低高度と同じ時間間隔で更新されるようにしている。これにより、高度によってインコヒーレント積分時間が異なっていても、出力されるドップラ速度の時間間隔は高度に依存しなくなる。
【００７５】
以上のように、この実施の形態４によれば、インコヒーレント積分における時間方向の積分範囲の更新幅を、インコヒーレント積分時間によらず一定としているので、インコヒーレント積分部２から出力されるインコヒーレント積分後のデータの時間分解能を、高度によらず一定とすることが可能になるという効果が得られる。
【００７６】
【発明の効果】
以上のように、この発明によれば、入力された受信信号のデータにＦＦＴ処理部でフーリエ変換を施して、積分時間設定部で高度毎に最適なインコヒーレント積分時間を設定し、そのフーリエ変換されたデータから、インコヒーレント積分部においてパワースペクトルを算出するとともに、設定されたインコヒーレント積分時間だけパワースペクトルを時間積分し、ドップラ速度算出部で、そのパワースペクトルから信号スペクトルのピークを検出して、そのドップラ周波数から大気の視線方向速度を算出し、品質管理・時間平均処理部にて、低品質のドップラ速度のデータを取り除くとともに時間平均を行うことによって、風速の測定を行うように構成したので、高度毎にインコヒーレント積分時間が設定され、高度毎のターゲットの性質に応じた適切な時間分解能で、高い検出性能の信号処理を実現することができる、ウィンドプロファイラにおける信号処理装置が得られるという効果がある。
【００７７】
この発明によれば、低高度におけるインコヒーレント積分時間を短くし、高高度におけるインコヒーレント積分時間を長く設定するように構成したので、地表面に近い低高度では高い時間分解能で風速ベクトルを測定することが可能になるとともに、ＳＮ比が低下する高高度でも高いデータ取得率を維持することができるという効果がある。
【００７８】
この発明によれば、データ取得率算出部がドップラ速度算出部の出力するドップラ速度データより算出したデータ取得率に応じて、積分時間設定部がインコヒーレント積分時間を設定するように構成したので、データ取得率の低い高度においてはインコヒーレント積分時間を増加させることにより、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することができるため、信号処理の性能を常に最適な状態に保つことが可能となり、また、上空風の特性が観測地によって異なっても、自動的にインコヒーレント積分時間を設定するため、設置時の信号処理設定が容易になるなどの効果がある。
【００７９】
この発明によれば、データ取得率算出部が品質管理・時間平均処理部の出力するドップラ速度データより算出したデータ取得率に応じて、積分時間設定部がインコヒーレント積分時間を設定するように構成したので、データ取得率の低い高度においてはインコヒーレント積分時間を増加させることにより、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することができるため、信号処理の性能を常に最適な状態に保つことが可能となり、また、上空風の特性が観測地によって異っても、自動的にインコヒーレント積分時間を設定するため、設置時の信号処理設定が容易になるなどの効果がある。
【００８０】
この発明によれば、データ取得率算出部の算出した高度毎のデータ取得率をもとに、積分時間設定部が各高度のインコヒーレント積分時間を設定するように構成したので、その高度のインコヒーレント積分時間をより的確に設定することができるという効果がある。
【００８１】
この発明によれば、データ取得率算出部からのあらかじめ分割された高度範囲毎のデータ取得率をもとに、積分時間設定部があらかじめ分割された高度範囲毎のインコヒーレント積分時間を設定するように構成したので、各高度毎にインコヒーレント積分時間を設定するよりも処理が簡略化され、インコヒーレント積分時間決定に要する処理負荷を低減することができるという効果がある。
【００８２】
この発明によれば、インコヒーレント積分時間が長く設定された高度では、データ補間処理部で時間方向にデータを補間し、ドップラ速度を高度によらず一定の時間分解能で品質管理・時間平均処理部に出力するように構成したので、全高度において同じ時間間隔で、ドップラ速度データあるいは風速データを得ることとができ、数値気象予測の初期値として風速データを利用する場合など、風速データを２次利用する際に利用しやすい形でデータ出力することが可能になるという効果がある。
【００８３】
この発明によれば、品質管理・時間平均処理部による時間平均処理を、インコヒーレント積分時間に反比例するデータ数を用いて行い、データの時間分解能が高度によらず一定となるドップラ速度を出力するように構成したので、全高度において同じ時間間隔で、ドップラ速度データあるいは風速データを得ることができ、数値気象予測の初期値として風速データを利用する場合など、風速データを２次利用する際に、利用しやすい形でデータを出力することが可能となるという効果がある。
【００８４】
この発明によれば、インコヒーレント積分における時間方向の積分範囲の更新幅を、インコヒーレント積分時間によらず一定とし、インコヒーレント積分部が出力する積分後のデータの時間分解能を、高度によらず一定となるように構成したので、全高度において同じ時間間隔で、ドップラ速度データあるいは風速データを得ることができ、数値気象予測の初期値として風速データを利用する場合などに、風速データを２次利用しやすい形で出力することが可能になるという効果がある。
【００８５】
この発明によれば、受信信号のデータを読込むデータ読込みステップと、受信信号のデータにフーリエ変換を施すＦＦＴステップと、高度毎のインコヒーレント積分時間を設定する積分時間設定ステップと、受信信号のデータのフーリエ変換からパワースペクトルを算出し、設定されたインコヒーレント積分時間だけインコヒーレント積分を実施するインコヒーレント積分ステップと、インコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出するドップラ速度算出ステップと、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを判定する全高度処理終了判定ステップと、処理が未終了の場合に処理対象となる高度を次の高度へと更新し、積分時間設定ステップに戻る高度更新ステップと、処理が終了している場合に処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、データ読込みステップに戻る時刻更新ステップとからなる信号処理方法により、高度毎にインコヒーレント積分時間を設定するように構成したので、高度毎のターゲットの性質に応じて適切な時間分解能で、かつ高い検出性能の信号処理を実現できるウィンドプロファイラにおける信号処理方法が得られるという効果がある。
【００８６】
この発明によれば、受信信号のデータを読込むデータ読込みステップと、受信信号のデータにフーリエ変換を施すＦＦＴステップと、インコヒーレント積分時間の初期値を設定する積分時間初期値設定ステップと、受信信号のデータのフーリエ変換からパワースペクトルを算出し、設定されたインコヒーレント積分時間だけインコヒーレント積分を実施するインコヒーレント積分ステップと、インコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出するドップラ速度算出ステップと、算出されたドップラ速度から、現在処理対象となっている高度のデータ取得率を算出するデータ取得率算出ステップと、算出されたデータ取得率をあらかじめ定めたしきい値と比較するデータ取得率判定ステップと、データ取得率がしきい値未満のとき、現時点のインコヒーレント積分時間とあらかじめ定められたインコヒーレント積分時間を比較する積分時間判定ステップと、現時点のインコヒーレント積分時間があらかじめ定められたインコヒーレント積分時間未満の場合に、インコヒーレント積分時間の設定を増加させて、インコヒーレント積分ステップに戻る積分時間増加ステップと、データ取得率がしきい値以上の場合、あるいは現時点のインコヒーレント積分時間があらかじめ定められたインコヒーレント積分時間以上の場合に、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを判定する全高度処理終了判定ステップと、処理が未終了の場合に処理対象となる高度を次の高度へと更新し、積分時間初期値設定ステップに戻る高度更新ステップと、処理が終了している場合に処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、データ読込みステップに戻る時刻更新ステップとからなる信号処理方法により、高度毎にインコヒーレント積分時間を設定するように構成したので、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することができるため、信号処理の性能を常に最適な状態に保つことが可能となり、また、上空風の特性が観測地によって異なっても、自動的にインコヒーレント積分時間を設定するため、設置時の信号処理設定が容易となるウィンドプロファイラにおける信号処理方法が得られるという効果がある。
【図面の簡単な説明】
【図１】この発明の実施の形態１のウィンドプロファイラにおける信号処理装置を示すブロック図である。
【図２】実施の形態１における風速ベクトル算出の原理を示す説明図である。
【図３】実施の形態１における大気乱流エコーの検出の原理を示す説明図である。
【図４】実施の形態１における信号処理装置で出力される速度データが高度によって異なる時間分解能を持つことを示す説明図である。
【図５】実施の形態１のウィンドプロファイラにおける信号処理方法を記すフローチャートである。
【図６】実施の形態１における信号処理装置の他の構成例を示すブロック図である。
【図７】この発明の実施の形態２のウィンドプロファイラにおける信号処理装置を示すブロック図。
【図８】実施の形態２のウィンドプロファイラにおける信号処理方法を示すフローチャートである。
【図９】実施の形態２における信号処理装置の他の構成例を示すブロック図である。
【図１０】この発明の実施の形態３のウィンドプロファイラにおける信号処理装置を示すブロック図である。
【図１１】実施の形態３における信号処理装置の他の構成例を示すブロック図である。
【図１２】この発明の実施の形態４におけるインコヒーレント積分の時間範囲の更新方法を示す説明図である。
【図１３】この発明および従来の信号処理装置および信号処理方法が適用される、一般的なウィンドプロファイラを示すブロック図である。
【図１４】従来のウィンドプロファイラにおける信号処理装置を示すブロック図である。
【符号の説明】
１ ＦＦＴ処理部、２ インコヒーレント積分部、３ ドップラ速度算出部、４ 品質管理・時間平均処理部、５ 積分時間設定部、６ データ補間処理部、７ データ取得率算出部、８ 高度毎時間平均数設定部、１０１ 空中線、１０２ 送受信装置、１０３ 信号処理装置、１０４ 風速ベクトル算出装置、１０５ 表示・記録装置。
【発明の属する技術分野】
この発明は、上空の風速分布を測定するウィンドプロファイラにおける信号処理装置および信号処理方法に関するものである。
【０００２】
【従来の技術】
風向・風速の情報は気象予報に必要な情報の一つである。風向・風速を計測する最も一般的な方法としては、地上に風向・風速計を設置することである。しかしながら、その場合には地表付近の風しか測定できない。気象予報をより正確にするためには上空の風向・風速のデータも知ることが必要である。そのため、従来、ゾンデ等を用いて上空の風向・風速の観測を行っていた。このゾンデを用いた観測では、ゾンデを上げた時刻における風向・風速のデータしか得ることができないため、観測の時間分解能が数時間以上と低いことが欠点であった。
【０００３】
それに対して、近年、ウィンドプロファイラと呼ばれる大気レーダを用いて上空の風向・風速を計測する技術が確立されつつある。ウィンドプロファイラでは、１分〜数分毎に上空の風向・風速を計測することが可能となる。このような高い時間分解能による上空の風向・風速のデータは、気象予報の精度向上に有効となると期待されている。このようなウィンドプロファイラについては、例えばＣａｒｔｅｒ ｅｔ ａｌ．， Ｄｅｖｅｌｏｐｍｅｎｔ ｉｎ ＵＨＦ ｌｏｗｅｒ ｔｒｏｐｏｓｈｅｒｉｃ ｗｉｎｄ ｐｒｏｆｉｌｉｎｇ ａｔ ＮＯＡＡ’ｓ Ａｅｒｏｎｏｍｙ Ｌａｂｏｒａｔｏｒｙ， Ｒａｄｉｏ Ｓｃｉｅｎｃｅ， ｖｏｌ．３０， ｎｏ．４， ｐｐ．９７〜１００１， １９９５．に紹介されたものなどがある。
【０００４】
ここで、ウィンドプロファイラで大気を計測する原理を説明する。
ウィンドプロファイラはドップラレーダの一種であり、一般的に図１３に示すように、空中線１０１、送受信装置１０２、信号処理装置１０３、風速ベクトル算出装置１０４、表示・記録装置１０５によって構成されている。このように構成されたウィンドプロファイラでは、空中線１０１より空中に放射された電磁波の反射波が空中線１０１にて受信され、その受信された電磁波が送受信装置１０２によって増幅、周波数変換、検波されてビデオ信号に変換される。このビデオ信号に信号処理装置１０３で周波数解析処理を行うことにより、ドップラ周波数を算出し、そのドップラ周波数からドップラ速度を算出して風速ベクトル算出装置１０４に送り、そこで算出された風速ベクトルを表示・記録装置１０５に表示もしくは記録する。
【０００５】
なお、上記従来のウィンドプロファイラでは信号処理装置１０３として、例えば図１４に示すような構成のものが一般に用いられている。この図１４において、１はＦＦＴ処理部、２はインコヒーレント積分部、３はドップラ速度算出部、４は品質管理・時間平均処理部である。
【０００６】
次にその動作について説明する。
送受信装置１０２より出力された受信信号のデータは、ＦＦＴ処理部１に入力される。ＦＦＴ処理部１では受信信号のデータに対してフーリエ変換を施す。このフーリエ変換後の受信信号のデータはインコヒーレント積分部２に入力される。ここでは、受信信号のデータのフーリエ変換から電力値を求めることによって得られるパワースペクトルを算出した後、複数の時刻に得られたパワースペクトルを積算する処理、すなわちインコヒーレント積分の処理を行う。受信信号のデータのパワースペクトルには、大気乱流エコーのスペクトル、すなわち信号スペクトルと雑音スペクトルが含まれる。このインコヒーレント積分処理により、雑音成分のゆらぎおよび信号成分のゆらぎの幅が小さくなる。そのため、インコヒーレント積分を行うことで信号検出の精度をあげることができる。ただし、インコヒーレント積分は複数の時刻のデータを用いるため、時間分解能を低下させることになる。
【０００７】
ドップラ速度算出部３では、インコヒーレント積分後の安定化されたパワースペクトルから信号スペクトルのピークを検出し、その中心ドップラ周波数から平均ドップラ周波数を算出し、さらに、そのドップラ周波数をドップラ速度に変換することにより、ターゲットである大気乱流の視線方向速度を算出する。
【０００８】
ただし、ドップラレーダで直接計測される風速は、実際の風速の視線方向への射影成分のみとなる。そこで、風速ベクトル算出装置１０４では、レーダ上空のある一定領域の風速分布が一様であると仮定して、レーダの観測方向を複数方向に変えて測定を行うことにより、３次元の風速ベクトルの合成を行う。
【０００９】
なお、このようなウィンドプロファイラでは、インコヒーレント積分時間が高度によらず一定である。しかし、ウィンドプロファイラのターゲットである大気乱流のレーダ反射率は、高度の上昇とともに急激に減少するため、実際には高度によって最適な積分時間が異なる。すなわち、ＳＮ比の高い低高度よりも、ＳＮ比が低下する高高度の方が長い積分時間を必要とする。
【００１０】
また、高高度の風速は地表面の影響をあまり受けないため、風速の時間変化が低高度の場合よりも小さい。このことからも、高高度ではインコヒーレント積分時間を長くすることが有利となる。
【００１１】
このように、ウィンドプロファイラのインコヒーレント積分時間は、高度によって最適な値が異なるものである。
【発明が解決しようとする課題】
【００１２】
従来のウィンドプロファイラにおける信号処理装置は以上のように構成されているので、インコヒーレント積分時間が高度によらず一定であり、高高度ではインコヒーレント積分時間が不足するため、十分な信号検出性能が得られず、データ取得率が低下するという課題があった。
【００１３】
この発明は上記のような課題を解決するためになされたもので、高度に応じてインコヒーレント積分時間を設定することにより、できるだけ広い高度範囲でデータ取得率を向上させた、ウィンドプロファイラにおける信号処理装置および信号処理方法を得ることを目的とする。
【００１４】
【課題を解決するための手段】
この発明に係るウィンドプロファイラにおける信号処理装置は、受信信号のデータにフーリエ変換を施すフーリエ変換処理部と、高度毎に最適なインコヒーレント積分時間を設定する積分時間設定部と、フーリエ変換されたデータからパワースペクトルを算出し、設定されたインコヒーレント積分時間だけそのパワースペクトルを時間積分するインコヒーレント積分部と、パワースペクトルから信号スペクトルのピークを検出し、そのドップラ周波数から大気の視線方向速度を算出するドップラ速度算出部と、低品質のドップラ速度のデータを取り除くとともに時間平均を行う品質管理・時間平均処理部とを備え、風速の測定を行うようにしたものである。
【００１５】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分時間を、低高度においては短く、高高度においては長く設定するようにしたものである。
【００１６】
この発明に係るウィンドプロファイラにおける信号処理装置は、積分時間設定部が、ドップラ速度算出部の出力するドップラ速度データよりデータ取得率を算出するデータ取得率算出部からのデータ取得率に応じて、インコヒーレント積分時間を設定するようにしたものである。
【００１７】
この発明に係るウィンドプロファイラにおける信号処理装置は、積分時間設定部が、品質管理・時間平均処理部の出力するドップラ速度データよりデータ取得率を算出するデータ取得率算出部からのデータ取得率に応じて、インコヒーレント積分時間を設定するようにしたものである。
【００１８】
この発明に係るウィンドプロファイラにおける信号処理装置は、各高度のインコヒーレント積分時間を、積分時間設定部が、データ取得率算出部にて算出されたその高度のデータ取得率をもとに設定するようにしたものである。
【００１９】
この発明に係るウィンドプロファイラにおける信号処理装置は、あらかじめ分割された高度範囲毎のインコヒーレント積分時間を、積分時間設定部が、データ取得率算出部にて算出されたあらかじめ分割された高度範囲毎のデータ取得率をもとに設定するようにしたものである。
【００２０】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分時間を長く設定した高度については、データ補間処理部で時間方向にデータを補間して、高度によらず一定の時間分解能のドップラ速度を品質管理・時間平均処理部に出力するようにしたものである。
【００２１】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分時間に反比例するデータ数を用いて、品質管理・時間平均処理部による時間平均処理を行い、データの時間分解能が高度によらず一定となるドップラ速度を出力するようにしたものである。
【００２２】
この発明に係るウィンドプロファイラにおける信号処理装置は、インコヒーレント積分の時間方向の積分範囲の更新幅をインコヒーレント積分時間によらず一定とし、インコヒーレント積分部からの積分後のデータの時間分解能を高度によらず一定となるようにしたものである。
【００２３】
この発明に係るウィンドプロファイラにおける信号処理方法は、受信信号のデータを読込んでフーリエ変換を施し、高度毎にインコヒーレント積分時間を設定し、フーリエ変換されたデータからパワースペクトルを算出して、設定されたインコヒーレント積分時間だけインコヒーレント積分を行い、そのインコヒーレント積分後のパワースペクトルから信号スペクトルを検出して、ドップラ周波数からドップラ速度を算出し、処理対象となっている時刻に関する全高度のデータの処理が未終了であれば、処理対象の高度を次の高度へ更新した後、インコヒーレント積分時間の設定に処理を戻し、終了していれば、処理対象の時刻を次の時刻に進めるとともに処理対象となる高度を初期化し、処理をデータの読込みに戻すことによって、風速の測定を行うようにしたものである。
【００２４】
この発明に係るウィンドプロファイラにおける信号処理方法は、受信信号のデータを読込んでフーリエ変換を施し、インコヒーレント積分時間の初期値を設定し、フーリエ変換されたデータからパワースペクトルを算出して、設定されたインコヒーレント積分時間だけインコヒーレント積分を行い、そのインコヒーレント積分後のパワースペクトルから信号スペクトルを検出して、ドップラ周波数からドップラ速度を算出し、そのドップラ速度から処理対象となっている高度のデータ取得率の算出を行って、そのデータ取得率があらかじめ定められたしきい値未満であり、かつインコヒーレント積分時間がその上限未満であれば、インコヒーレント積分時間の設定を増加させて、インコヒーレント積分の処理に戻し、それ以外の場合には、処理対象となっている時刻に関する全高度のデータの処理が終了していることを確認し、未終了であれば、処理対象の高度を次の高度へ更新した後、インコヒーレント積分時間の設定に処理を戻し、終了していれば、処理対象の時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、処理をデータの読込みに戻すことによって、風速の測定を行うようにしたものである。
【００２５】
【発明の実施の形態】
以下、この発明の実施の一形態を説明する。
実施の形態１．
図１はこの発明の実施の形態１によるウィンドプロファイラにおける信号処理装置を示すブロック図である。なお、ウィンドプロファイラの全体構成は図１３に示した従来の場合と同様であり、この図１に示した信号処理装置１０３は、図１３に示されるウィンドプロファイラの一部分を形成している。
【００２６】
図１３において、１０１は電磁波の送受信を行う空中線であり、１０２は空中線１０１で送受信される電磁波の入出力処理を行う送受信装置である。１０３は図１にその構成を示すこの発明の実施の形態１による信号処理装置である。１０４は信号処理装置１０３の出力より３次元の風速ベクトルの合成を行う風速ベクトル算出装置であり、１０５はこの風速ベクトル算出装置１０４の計算結果の表示あるいは記録を行う表示・記録装置である。
【００２７】
また、図１において、１は送受信装置１０２で受信された受信信号のデータにフーリエ変換を施すフーリエ変換処理部であり、ここでは、例えば高速フーリエ変換（ＦＦＴ）を施すＦＦＴ処理部が用いられている。２はＦＦＴ処理部１にてフーリエ変換されて出力された受信信号のデータからパワースペクトルを算出するとともに、複数の時刻で得られたパワースペクトルをあつめて積分するインコヒーレント積分部である。３はインコヒーレント積分部２から出力されたパワースペクトルから信号スペクトルを検出し、その中心周波数からターゲットのドップラ速度、すなわちターゲットの視線方向速度を算出するドップラ速度算出部である。４はドップラ速度算出部３で算出されたドップラ速度データのうちから品質の低下したデータを除去した後、その時間平均を行って風速ベクトル算出装置１０４に出力する品質管理・時間平均処理部である。５は高度に応じてインコヒーレント積分時間を設定し、インコヒーレント積分部２に積分時間を指示する積分時間設定部である。
【００２８】
次に動作について説明する。
図１３に示すように構成されたウィンドプロファイラでは、送受信装置１０２で生成した電磁波を空中線１０１より空中に放射する。空中に放射された電磁波は大気の屈折率の粗密が散乱体となって反射される。反射された電磁波は空中線１０１によって受信され、送受信装置１０２に入力される。散乱体が上空の風と共に流れていると、受信された電磁波はドップラ効果によって周波数が変化する。この周波数の変化を一般のドップラレーダと同様に検出することにより、上空の風速を測定する。具体的には、送受信装置１０２によって、受信された電磁波の増幅、周波数変換、検波を行うことによりビデオ信号に変換し、それを信号処理装置１０３へ出力する。信号処理装置１０３では、そのビデオ信号に周波数解析処理を行うことにより、ドップラ周波数を算出する。さらにドップラ周波数からドップラ速度を算出する。
【００２９】
以下に、実施の形態１による信号処理装置１０３の詳細な処理動作について説明する。
図１３に示す送受信装置１０２から出力された受信信号のデータは、図１のＦＦＴ処理部１に入力されてフーリエ変換が施される。ＦＦＴ処理部１にてフーリエ変換された受信信号のデータはインコヒーレント積分部２に入力される。インコヒーレント積分部２では、周波数毎に振幅の絶対値の二乗を計算することによってパワースペクトルを算出する。インコヒーレント積分部２はさらに、複数の時刻のデータから算出されたパワースペクトルを集めてそれを積分する。なお、積分処理において積算する時間方向のデータ数は、積分時間設定部５において高度毎に設定された積分時間によって定める。
【００３０】
インコヒーレント積分後のパワースペクトルはインコヒーレント積分部２よりドップラ速度算出部３に入力される。ドップラ速度算出部３では、インコヒーレント積分後の安定化されたパワースペクトルから信号スペクトルのピークを検出して、その中心ドップラ周波数から平均ドップラ周波数を算出し、さらに、そのドップラ周波数をドップラ速度に変換することにより、ターゲットである大気乱流の視線方向速度を算出する。
【００３１】
ただし、ドップラレーダで直接計測される風速は、実際の風速の視線方向への射影成分のみとなる。そこで、レーダ上空のある一定領域の風速分布が一様であると仮定して、レーダの観測方向を複数方向に変えて測定を行い、風速ベクトル算出装置１０４でそれを３次元の風速ベクトルに合成する。
【００３２】
この３次元風速ベクトルの算出原理を図２を用いて説明する。図２（ａ）はウィンドプロファイラにおけるビーム方向の例である。この例では天頂、および東、西、南、北の計５つの方向にビームを向けて観測している。天頂方向以外のビームの天頂角は、例えば１０度程度に設定される。なお、３次元風速ベクトルは東西方向成分、南北方向成分、鉛直方向成分に分けることができる。鉛直方向成分は天頂方向に指向したビーム（以下天頂ビームという）で得られるドップラ速度そのものとなる。東西方向成分は東方向に指向したビーム（以下東ビームといい、他の方角についても同様とする）のドップラ速度と西ビームのドップラ速度を用いて求めることができる。
【００３３】
図２（ｂ）は、東西方向に水平に風が吹いている場合に、東ビームと西ビームでどのようなドップラ速度が観測されるかを示したものである。図示のように、西から東に風が吹いている場合、西ビームでは近づく方向のドップラ速度が観測され、東ビームでは遠ざかる方向のドップラ速度が観測される。このように、異なるビーム方向で同一の風を観測すると、異なるドップラ速度が得られる。従って、ビーム間のドップラ速度の差異から水平風成分を算出することが可能となる。同様にして、北ビームと南ビームを用いれば、風速の南北方向成分を算出することができる。
【００３４】
なお、図３はウィンドプロファイラで観測される受信信号のデータのパワースペクトルを模式的に示した説明図である。ＳＮ比が高い場合には、信号スペクトルのピークを容易に検出することができるが、ＳＮ比が小さくなると、信号スペクトルのピーク値が低くなり、雑音のゆらぎと識別しずらくなる。インコヒーレント積分部２による積分は、雑音および信号の統計的なゆらぎを小さくする効果がある。そのため、積分時間設定部５によってインコヒーレント積分時間を長くすれば、図３に示した雑音電力密度の標準偏差が小さくなり、ＳＮ比が小さく、信号スペクトルのピーク値が低い場合にも、信号スペクトルのピークを検出することが可能となる。
【００３５】
ドップラ速度算出部３はこのようにしてパワースペクトルから大気乱流エコーのピークを検出し、その中心ドップラ周波数から、大気乱流エコーのドップラ速度、すなわち視線方向風速を求める。ドップラ速度算出部３から出力された視線方向風速は品質管理・時間平均処理部４において低品質なデータを除去した後、時間平均処理が施されて風速ベクトル算出装置１０４へと出力される。風速ベクトル算出装置１０４では、異なるビーム方向で観測したデータを合成することによって３次元風速ベクトルを算出する。この風速ベクトル算出装置１０４の計算結果は、表示・記録装置１０５に送られて表示あるいは記録される。
【００３６】
積分時間設定部５では、高度毎に最適となるインコヒーレント積分時間を設定する。低高度の範囲ではレーダとターゲットの距離が近いこと、および大気濃度が濃いため、大気乱流による電磁波の反射率が高い。そのため、比較的高いＳＮ比の受信信号が観測される。よって、インコヒーレント積分時間が短くても、パワースペクトルから信号スペクトルを容易に検出することができる。また、低高度の風は地表面の影響を受け易く、風向・風速の時間変化が短時間で生じる。そのため、３次元風速ベクトルの測定は高い時間分解能で行うことが望ましい。インコヒーレント積分後に時間平均処理を行う場合においても、地表面の影響を受けてある時刻のみに特異な風が生じた場合においても、高い時間分解能で風速ベクトルを求めておけば、特異な風速を品質管理処理によって除去し、卓越する傾向の風速ベクトルデータのみを用いて平均風速を求めることができる。
【００３７】
一方、高高度の範囲ではレーダとターゲットの距離が大きいこと、および大気濃度が薄くなって大気乱流の電磁波反射率が急激に減少することにより、受信信号のＳＮ比は低くなることが多い。従って、信号スペクトルを正確に検出するためには、インコヒーレント積分数を多くしてパワースペクトルのゆらぎを小さくする必要がある。また、高高度では、風が地表面の影響をあまり受けないことから、風速ベクトルの時間変化は小さい。そのことから、風速ベクトル算出の時間分解能は低高度の場合より低くてもあまり問題とならない。
【００３８】
以上のことから、インコヒーレント積分時間を、低高度で短く、高高度で長くすることが、一般的には適切な信号処理方法となると考えられる。典型的な例を図４に示す。この図では、インコヒーレント積分時間を、低高度である高度範囲１では１分、中程度の高度である高度範囲２では５分、高高度である高度範囲３では１０分とすることを示している。
【００３９】
以下、このような信号処理装置１０３の動作について、図５に示すフローチャートを用いて詳細に説明する。
まずデータ読込みステップＳＴ１において、送受信装置１０２からの受信信号のデータを信号処理装置１０３に入力する。信号処理装置１０３ではその入力された受信信号のデータに対して、フーリエ変換ステップ（以下ＦＦＴステップという）ＳＴ２において高速フーリエ変換を施す。次に積分時間設定ステップＳＴ３に進み、高度毎のインコヒーレント積分時間を設定する。次にインコヒーレント積分ステップＳＴ４において、ＦＦＴステップＳＴ２で得られた受信信号のデータのフーリエ変換からパワースペクトルを算出し、積分時間設定ステップＳＴ３で設定されたインコヒーレント積分時間だけ、インコヒーレント積分を実施する。
【００４０】
次にドップラ速度算出ステップＳＴ５において、インコヒーレント積分ステップＳＴ４で算出されたインコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出する。その後、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを、全高度処理終了判定ステップＳＴ６において判定する。判定の結果、全高度についての処理が済んでいなければ、高度更新ステップＳＴ７に進んで、処理対象となる高度を次の高度に更新した後、処理を積分時間設定ステップＳＴ３に戻す。一方、全高度についての処理が済んでいれば、時刻更新ステップＳＴ８において処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化して、処理をデータ読込みステップＳＴ１に戻す。
【００４１】
なお、品質管理・時間平均処理部４では品質管理処理および時間平均処理を実施する。この品質管理処理は品質が劣化したと判断されるデータを除去するための処理であり、ある時間・空間範囲でのデータの連続性を仮定することにより、不連続なデータを除去して品質の高いデータのみを残している。
【００４２】
この品質管理・時間平均処理部４で不連続なデータを除去する際、従来の品質管理処理では、データが全ての高度で同じ時間分解能で得られるということを前提にしているか、もしくは全ての高度で同じ時間分解能であることが望ましい処理方式となっていることがある。そのような場合、品質管理・時間平均処理部４において、入力されたドップラ速度データに補間処理を行うことも有効であると考えられる。
【００４３】
このようなデータ補間処理が行われる、この発明の実施の形態１におけるウィンドプロファイラの信号処理装置の構成例を図６に示す。図において、６はドップラ速度算出部３より入力されたドップラ速度データに補間処理を施し、品質管理・時間平均処理部４に出力するデータ補間処理部である。
【００４４】
このデータ補間処理部６では、全高度範囲のうち、最も時間分解能の高い高度における時間分解能と同一の時間分解能となるように、他の高度において時間方向にデータ補間を行う。例えば、インコヒーレント積分時間を長くしたことによって、品質管理処理を行う時間範囲に時間方向に一つのデータしか存在しない場合には、全補間点で同一の風速値を仮定すればよい。
【００４５】
以上のように、この実施の形態１によれば、高度毎にインコヒーレント積分時間の設定を行っているので、高度毎のターゲットの性質に応じて、適切な時間分解能で、かつ検出精度の高い信号処理を実現することができるという効果が得られる。
【００４６】
また、インコヒーレント積分時間を、低高度では短く、高高度では長く設定しているので、地表面に近い低高度では高い時間分解能で風速ベクトルを測定するとともに、ＳＮ比が低下する高高度でも検出精度を高めることが可能になるという効果も得られる。
【００４７】
さらに、データ補間処理を組合せることにより、高度毎にインコヒーレント積分時間が異なる場合にも、出力されるドップラ速度の時間分解能は高度によらず一定となるという効果も得られる。
【００４８】
なお、全高度において同じ時間間隔でドップラ速度データ、あるいは風速データが得られれば、数値気象予測の初期値として風速データを利用することをはじめとする風速データの２次利用において、データを利用しやすい形で出力することが可能になるという効果が得られる。
【００４９】
実施の形態２．
ここで、大気の状態は天候によって変化するため、最適なインコヒーレント積分時間は時間とともに変化する。上記実施の形態１では、高度毎にインコヒーレント積分時間をあらかじめ定めるものについて説明したが、大気の状態によって適応的に高度毎にインコヒーレント積分時間を設定するようにしてもよい。この実施の形態２は、そのような大気の状態によって適応的に高度毎にインコヒーレント積分時間を設定する、ウィンドプロファイラにおける信号処理装置および信号処理方法に関するものである。
【００５０】
図７はこの発明の実施の形態２による、そのようなウィンドプロファイラの信号処理装置を示すブロック図である。図において、１０３は信号処理装置、１はＦＦＴ処理部、２はインコヒーレント積分部、３はドップラ速度算出部、４は品質管理・時間平均処理部、５は積分時間設定部であり、これらは図１に同一符号を付したものと同等のものであるため、詳細な説明は省略する。また、７はドップラ速度算出部３の出力するドップラ速度データを入力してデータ取得率を算出し、それを積分時間設定部５に入力するデータ取得率算出部である。なお、積分時間設定部５は、データ取得率算出部７で算出されたデータ取得率をもとにインコヒーレント積分時間を設定している点で、図１に同一符号を付したものとは異なっている。
【００５１】
次に動作について説明する。
ここで、この図７に示す構成の信号処理装置１０３において、ＦＦＴ処理部１による受信信号のデータのフーリエ変換、インコヒーレント積分部２によるインコヒーレント積分、ドップラ速度算出部３による信号スペクトルの検出およびドップラ速度算出、品質管理・時間平均処理部４によるドップラ速度の時間平均の各処理については、図１に示した実施の形態１の場合と同様の動作となる。なお、この実施の形態２において特徴的であるのは、データ取得率算出部７の算出したデータ取得率をもとに、積分時間設定部５がインコヒーレント積分時間を設定することにある。
【００５２】
ドップラ速度算出部３の算出した各高度、各ビームでのドップラ速度は、品質管理・時間平均処理部４とともに、データ取得率算出部７にも入力される。このドップラ速度算出部３で得られたドップラ速度は、ＳＮ比の低下などにより信号スペクトルが検出されなかった場合には、その高度およびビームのドップラ速度は欠損となる。データ取得率算出部７では、このドップラ速度算出部３からのドップラ速度データより、全データに対する欠損でないデータの比率をデータ取得率として算出する。なお、このデータ取得率は各高度、各ビーム毎に、複数の時刻のデータから算出することができる。
【００５３】
また、一つの時刻のデータからデータ取得率を算出するために、観測高度範囲をいくつかの高度範囲に分割し、その高度範囲内のデータからデータ取得率を算出するようにしてもよい。その場合、データ取得率算出部７であらかじめ分割された高度範囲毎のデータ取得率を算出し、その高度範囲毎のデータ取得率をもとに、積分時間設定部５にてあらかじめ分割された高度範囲毎のインコヒーレント積分時間の設定を行う。
【００５４】
データ取得率算出部７で算出されたデータ取得率は、積分時間設定部５に入力される。積分時間設定部５では、あらかじめ設定したしきい値よりも取得率が低くなる高度においてインコヒーレント積分時間を増加させ、標準値よりも大きくなるようにする。データ取得率に対するしきい値は全高度で同じでもよいし、あるいは高度によって異なるしきい値を用いてもよい。
【００５５】
以下、このような信号処理装置１０３の動作について、図８に示すフローチャートを用いて詳細に説明する。
まずデータ読込みステップＳＴ１０にて、送受信装置１０２の出力する受信信号のデータを信号処理装置１０３に入力し、ＦＦＴステップＳＴ１１においてその入力された受信信号のデータに対して高速フーリエ変換を施す。次に積分時間初期値設定ステップＳＴ１２に進み、インコヒーレント積分時間の初期値の設定を行う。次にインコヒーレント積分ステップＳＴ１３にて、ＦＦＴステップＳＴ１１で得られたフーリエ変換された受信信号のデータからパワースペクトルを算出し、積分時間初期値設定ステップＳＴ１２において設定されたインコヒーレント積分時間だけ、インコヒーレント積分を実施する。
【００５６】
次にドップラ速度算出ステップＳＴ１４において、上記インコヒーレント積分ステップＳＴ１３で算出されたインコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出する。次にデータ取得率算出ステップＳＴ１５に進み、上記ドップラ速度算出ステップＳＴ１４にて算出されたドップラ速度から、現在処理対象となっている高度のデータ取得率を算出する。その後、データ取得率判定ステップＳＴ１６において、上記データ取得率算出ステップＳＴ１５にて算出されたデータ取得率をあらかじめ定めたしきい値と比較する。
【００５７】
比較の結果、ドップラ速度のデータ取得率がしきい値を下回る場合には積分時間判定ステップＳＴ１７へ進み、しきい値以上のデータ取得率でドップラ速度が得られている場合には、後述する全高度処理終了判定ステップＳＴ１９に進む。積分時間判定ステップＳＴ１７では、現時点のインコヒーレント積分時間があらかじめ定められたインコヒーレント積分時間（上限）を超えているか否かを識別し、現時点のインコヒーレント積分時間が上限を超えている場合には後述する全高度処理終了判定ステップＳＴ１９に進み、それ以外の場合には積分時間増加ステップＳＴ１８に進む。この積分時間増加ステップＳＴ１８ではインコヒーレント積分時間の設定を増加させた後、インコヒーレント積分ステップＳＴ１３に処理を戻す。
【００５８】
また、上記全高度処理終了判定ステップＳＴ１９では、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かについて判定する。判定の結果、全高度についての処理が済んでいなければ、高度更新ステップＳＴ２０に進んで、処理対象となる高度を次の高度へと更新し、処理を積分時間初期値設定ステップＳＴ１２に戻す。一方、全高度についての処理が済んでいれば、時刻更新ステップＳＴ２１において処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化して、処理をデータ読込みステップＳＴ１０に戻す。
【００５９】
以上の処理ステップについて、データ数や積分時間をもう少し具体的に記述することにより補足する。
データ読込みステップＳＴ１０において入力される受信信号のデータは、インコヒーレント積分時間の上限と等しい時間分のデータ量を有している。ここで、このインコヒーレント積分時間の上限をＴｍａｘとする。また、ＦＦＴステップＳＴ１１にて、各高度において、時間方向にＦＦＴ点数で区切られたデータに対してフーリエ変換を施す。次に積分時間初期値設定ステップＳＴ１２にてインコヒーレント積分時間の最小値となるインコヒーレント積分時間の初期値を設定する。
【００６０】
次にインコヒーレント積分ステップＳＴ１３では、フーリエ変換されたデータをインコヒーレント積分時間毎にまとめ、インコヒーレント積分を行う。今、インコヒーレント積分時間がＴｉであるとすると、１つの高度に対してＴｍａｘ／Ｔｉ個のインコヒーレント積分後のパワースペクトルが得られる。それぞれのパワースペクトルから信号スペクトルピークを検出し、ドップラ速度算出ステップＳＴ１４にてそのドップラ速度を算出する。この時点で最大Ｔｍａｘ／Ｔｉ個のドップラ速度が算出される。ただし、ＳＮ比が低い場合などには、信号スペクトルが検出できない場合があるため、実際に得られるドップラ速度はＴｍａｘ／Ｔｉ個よりも小さい。
【００６１】
そこで、データ取得率算出ステップＳＴ１５において、現在注目している高度のデータ取得率を算出する。次にデータ取得率判定ステップＳＴ１６でそのデータ取得率をしきい値と比較し、データ取得率があらかじめ設定したしきい値よりも大きければ、全高度処理終了判定ステップＳＴ１９に進む。逆にデータ取得率がしきい値よりも小さければ、積分時間判定ステップＳＴ１７において、現時点のインコヒーレント積分時間が上限を超えていないことを確認する。上限を超えていないことが確認されれば、積分時間増加ステップＳＴ１８にてインコヒーレント積分時間を増加させて、再度インコヒーレント積分ステップＳＴ１３に戻ってインコヒーレント積分からの処理を繰り返す。以降、このインコヒーレント積分ステップＳＴ１３から積分時間増加ステップＳＴ１８までの処理を、処理中の高度のデータ取得率がしきい値を超えるか、あるいはインコヒーレント積分時間が上限を超えるまで繰り返す。
【００６２】
以上、１高度毎にデータ取得率を算出する場合についての処理フローを説明したが、高度範囲毎にデータ取得率を算出する場合には、ＦＦＴステップＳＴ１１から積分時間増加ステップＳＴ１８までの処理を１つの高度範囲毎に実施し、高度更新ステップＳＴ２０で高度範囲の単位で高度の更新を行うことになる。この場合、高度範囲毎にインコヒーレント積分時間の修正処理を行っているため、１高度毎にインコヒーレント積分時間を設定するのに比べて処理が簡略化され、信号処理の負荷を小さくすることができる。
【００６３】
なお、上記実施の形態２の説明では、データ取得率の算出を、ドップラ速度算出部３から出力されたドップラ速度データを用いて行っているものについて示したが、品質管理・時間平均処理部４から出力された時間平均後のドップラ速度をデータ取得率算出部７に入力して、データ取得率を算出するようにしてもよい。図９のブロック図にはそのような構成によるウィンドプロファイラにおける信号処理装置を示す。この場合、データ取得率算出部７へのデータ入力が、ドップラ速度算出部３からではなく、品質管理・時間平均処理部４から行われているだけで、その他は図７と同一である。
【００６４】
また、この実施の形態２においても、上記実施の形態１の場合と同様に、ドップラ速度算出部３と品質管理・時間平均処理部４との間にデータ補間処理部６を設け、ドップラ速度算出部３より入力されたドップラ速度データに補間処理を施して、品質管理・時間平均処理部４に高度によらず一定の時間分解能のドップラ速度を出力するようにしてもよい。
【００６５】
以上のように、この実施の形態２によれば、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することが可能になるため、信号処理の性能を常に最適な状態に保つことができ、また、上空風の特性が観測地によって異なっても、この実施の形態２では自動的にインコヒーレント積分時間を設定しているため、設置時の信号処理設定が容易になるなどの効果が得られる。
【００６６】
実施の形態３．
なお、上記各実施の形態では、高度毎にインコヒーレント積分時間が異なるため、最終的に得られる風速データの時間間隔は高度によってまちまちになることがあるが、気象予測のためなどにデータを利用する際は、全高度において同じ時間間隔で風速データが得られる方が便利なことが多いと考えられる。そのため、最終的に全ての高度で同じ時間間隔で風速データが出力されるように、インコヒーレント積分後のデータの時間分解能の差を、時間平均を行う際の時間平均数の差によって吸収することも考えられる。この実施の形態３は、そのようなウィンドプロファイラにおける信号処理装置に関するものである。
【００６７】
図１０はこの発明の実施の形態３による、そのようなウィンドプロファイラの信号処理装置を示すブロック図であり、実施の形態２の各部に相当する部分には図７と同一符号を付してその説明を省略する。図において、８は積分時間設定部５から出力される各高度におけるインコヒーレント積分時間を受けて、そのインコヒーレント積分時間に反比例するデータ数を各高度の時間平均数として設定し、品質管理・時間平均処理部４に各高度の時間平均数を指示する高度毎時間平均数設定部である。
【００６８】
次に動作について説明する。
ここで、この図１０に示した信号処理装置１０３における、ＦＦＴ処理部１による受信信号のデータのフーリエ変換、インコヒーレント積分部２によるインコヒーレント積分、ドップラ速度算出部３による信号スペクトルの検出およびドップラ速度算出、データ取得率算出部７によるデータ取得率の算出、および積分時間設定部５による高度毎のインコヒーレント積分時間の設定の各処理については、図７に示した実施の形態２と同一の動作となる。この実施の形態３で特徴的なのは、データを平均する時間が高度によらず一定となるように、インコヒーレント積分時間に反比例する数のデータを用いて時間平均処理を行うようにしたことにある。
【００６９】
データ取得率算出部７ではドップラ速度算出部３からのドップラ速度データよりデータ取得率を算出し、積分時間設定部５に入力する。積分時間設定部５はこのデータ取得率算出部７の算出したデータ取得率をもとに、各高度範囲のインコヒーレント積分時間を設定する。この積分時間設定部５から出力される各高度におけるインコヒーレント積分時間はインコヒーレント積分部２とともに、高度毎時間平均数設定部８にも入力される。高度毎時間平均数設定部８では入力されたインコヒーレント積分時間に反比例するデータ数を各高度の時間平均数に設定する。品質管理・時間平均処理部４はこの時間平均数を用いて時間平均処理を行う。その結果、品質管理・時間平均処理部４から出力されるドップラ速度データは高度によらず同じ時間間隔で出力される。
【００７０】
なお、上記実施の形態３の説明では、ドップラ速度算出部３から出力されたドップラ速度データを用いてデータ取得率を算出しているものを示したが、品質管理・時間平均処理部４から出力された時間平均後のドップラ速度をデータ取得率算出部７に入力して、データ取得率を算出するようにしてもよい。図１１のブロック図にはそのような構成によるウィンドプロファイラにおける信号処理装置を示す。この場合、データ取得率算出部７へのデータ入力を、ドップラ速度算出部３に代えて品質管理・時間平均処理部４から行っているだけで、その他は図１０と同一である。
【００７１】
また、この実施の形態３においても、上記各実施の形態の場合と同様に、ドップラ速度算出部３と品質管理・時間平均処理部４との間にデータ補間処理部６を設け、ドップラ速度算出部３より入力されたドップラ速度データに補間処理を施して、品質管理・時間平均処理部４に高度によらず一定の時間分解能のドップラ速度を出力するようにしてもよい。
【００７２】
以上のように、この実施の形態３によれば、高度毎時間平均数設定部８より品質管理・時間平均処理部４に、各高度の時間平均数を指示しているので、全高度において同じ時間間隔でドップラ速度データあるいは風速データを得ることができるため、数値気象予測の初期値として風速データを利用することをはじめとする風速データの２次利用において、データを利用しやすい形で出力することが可能になるという効果が得られる。
【００７３】
実施の形態４．
また、上記実施の形態３では、ドップラ速度データが高度によらず同じ時間間隔で出力されるものについて説明したが、この発明はそれにのみ限定されるものではない。以下、図１２を用いてそのようなこの発明の実施の形態４について説明する。
【００７４】
図１２はこの発明の実施の形態４におけるインコヒーレント積分の時間範囲の更新方法を示す説明図である。この図１２では、低高度でインコヒーレント積分時間が短く、高高度でインコヒーレント積分が長く設定されている状況を想定している。高高度ではインコヒーレント積分の時間範囲を更新する際に、前後の時間帯のインコヒーレント積分において、インコヒーレント積分時間範囲が重なるようにし、低高度と同じ時間間隔で更新されるようにしている。これにより、高度によってインコヒーレント積分時間が異なっていても、出力されるドップラ速度の時間間隔は高度に依存しなくなる。
【００７５】
以上のように、この実施の形態４によれば、インコヒーレント積分における時間方向の積分範囲の更新幅を、インコヒーレント積分時間によらず一定としているので、インコヒーレント積分部２から出力されるインコヒーレント積分後のデータの時間分解能を、高度によらず一定とすることが可能になるという効果が得られる。
【００７６】
【発明の効果】
以上のように、この発明によれば、入力された受信信号のデータにＦＦＴ処理部でフーリエ変換を施して、積分時間設定部で高度毎に最適なインコヒーレント積分時間を設定し、そのフーリエ変換されたデータから、インコヒーレント積分部においてパワースペクトルを算出するとともに、設定されたインコヒーレント積分時間だけパワースペクトルを時間積分し、ドップラ速度算出部で、そのパワースペクトルから信号スペクトルのピークを検出して、そのドップラ周波数から大気の視線方向速度を算出し、品質管理・時間平均処理部にて、低品質のドップラ速度のデータを取り除くとともに時間平均を行うことによって、風速の測定を行うように構成したので、高度毎にインコヒーレント積分時間が設定され、高度毎のターゲットの性質に応じた適切な時間分解能で、高い検出性能の信号処理を実現することができる、ウィンドプロファイラにおける信号処理装置が得られるという効果がある。
【００７７】
この発明によれば、低高度におけるインコヒーレント積分時間を短くし、高高度におけるインコヒーレント積分時間を長く設定するように構成したので、地表面に近い低高度では高い時間分解能で風速ベクトルを測定することが可能になるとともに、ＳＮ比が低下する高高度でも高いデータ取得率を維持することができるという効果がある。
【００７８】
この発明によれば、データ取得率算出部がドップラ速度算出部の出力するドップラ速度データより算出したデータ取得率に応じて、積分時間設定部がインコヒーレント積分時間を設定するように構成したので、データ取得率の低い高度においてはインコヒーレント積分時間を増加させることにより、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することができるため、信号処理の性能を常に最適な状態に保つことが可能となり、また、上空風の特性が観測地によって異なっても、自動的にインコヒーレント積分時間を設定するため、設置時の信号処理設定が容易になるなどの効果がある。
【００７９】
この発明によれば、データ取得率算出部が品質管理・時間平均処理部の出力するドップラ速度データより算出したデータ取得率に応じて、積分時間設定部がインコヒーレント積分時間を設定するように構成したので、データ取得率の低い高度においてはインコヒーレント積分時間を増加させることにより、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することができるため、信号処理の性能を常に最適な状態に保つことが可能となり、また、上空風の特性が観測地によって異っても、自動的にインコヒーレント積分時間を設定するため、設置時の信号処理設定が容易になるなどの効果がある。
【００８０】
この発明によれば、データ取得率算出部の算出した高度毎のデータ取得率をもとに、積分時間設定部が各高度のインコヒーレント積分時間を設定するように構成したので、その高度のインコヒーレント積分時間をより的確に設定することができるという効果がある。
【００８１】
この発明によれば、データ取得率算出部からのあらかじめ分割された高度範囲毎のデータ取得率をもとに、積分時間設定部があらかじめ分割された高度範囲毎のインコヒーレント積分時間を設定するように構成したので、各高度毎にインコヒーレント積分時間を設定するよりも処理が簡略化され、インコヒーレント積分時間決定に要する処理負荷を低減することができるという効果がある。
【００８２】
この発明によれば、インコヒーレント積分時間が長く設定された高度では、データ補間処理部で時間方向にデータを補間し、ドップラ速度を高度によらず一定の時間分解能で品質管理・時間平均処理部に出力するように構成したので、全高度において同じ時間間隔で、ドップラ速度データあるいは風速データを得ることとができ、数値気象予測の初期値として風速データを利用する場合など、風速データを２次利用する際に利用しやすい形でデータ出力することが可能になるという効果がある。
【００８３】
この発明によれば、品質管理・時間平均処理部による時間平均処理を、インコヒーレント積分時間に反比例するデータ数を用いて行い、データの時間分解能が高度によらず一定となるドップラ速度を出力するように構成したので、全高度において同じ時間間隔で、ドップラ速度データあるいは風速データを得ることができ、数値気象予測の初期値として風速データを利用する場合など、風速データを２次利用する際に、利用しやすい形でデータを出力することが可能となるという効果がある。
【００８４】
この発明によれば、インコヒーレント積分における時間方向の積分範囲の更新幅を、インコヒーレント積分時間によらず一定とし、インコヒーレント積分部が出力する積分後のデータの時間分解能を、高度によらず一定となるように構成したので、全高度において同じ時間間隔で、ドップラ速度データあるいは風速データを得ることができ、数値気象予測の初期値として風速データを利用する場合などに、風速データを２次利用しやすい形で出力することが可能になるという効果がある。
【００８５】
この発明によれば、受信信号のデータを読込むデータ読込みステップと、受信信号のデータにフーリエ変換を施すＦＦＴステップと、高度毎のインコヒーレント積分時間を設定する積分時間設定ステップと、受信信号のデータのフーリエ変換からパワースペクトルを算出し、設定されたインコヒーレント積分時間だけインコヒーレント積分を実施するインコヒーレント積分ステップと、インコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出するドップラ速度算出ステップと、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを判定する全高度処理終了判定ステップと、処理が未終了の場合に処理対象となる高度を次の高度へと更新し、積分時間設定ステップに戻る高度更新ステップと、処理が終了している場合に処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、データ読込みステップに戻る時刻更新ステップとからなる信号処理方法により、高度毎にインコヒーレント積分時間を設定するように構成したので、高度毎のターゲットの性質に応じて適切な時間分解能で、かつ高い検出性能の信号処理を実現できるウィンドプロファイラにおける信号処理方法が得られるという効果がある。
【００８６】
この発明によれば、受信信号のデータを読込むデータ読込みステップと、受信信号のデータにフーリエ変換を施すＦＦＴステップと、インコヒーレント積分時間の初期値を設定する積分時間初期値設定ステップと、受信信号のデータのフーリエ変換からパワースペクトルを算出し、設定されたインコヒーレント積分時間だけインコヒーレント積分を実施するインコヒーレント積分ステップと、インコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出するドップラ速度算出ステップと、算出されたドップラ速度から、現在処理対象となっている高度のデータ取得率を算出するデータ取得率算出ステップと、算出されたデータ取得率をあらかじめ定めたしきい値と比較するデータ取得率判定ステップと、データ取得率がしきい値未満のとき、現時点のインコヒーレント積分時間とあらかじめ定められたインコヒーレント積分時間を比較する積分時間判定ステップと、現時点のインコヒーレント積分時間があらかじめ定められたインコヒーレント積分時間未満の場合に、インコヒーレント積分時間の設定を増加させて、インコヒーレント積分ステップに戻る積分時間増加ステップと、データ取得率がしきい値以上の場合、あるいは現時点のインコヒーレント積分時間があらかじめ定められたインコヒーレント積分時間以上の場合に、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを判定する全高度処理終了判定ステップと、処理が未終了の場合に処理対象となる高度を次の高度へと更新し、積分時間初期値設定ステップに戻る高度更新ステップと、処理が終了している場合に処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、データ読込みステップに戻る時刻更新ステップとからなる信号処理方法により、高度毎にインコヒーレント積分時間を設定するように構成したので、気象状況の変化によってデータ取得率が変化した場合にも、適切なインコヒーレント積分時間を適応的に設定することができるため、信号処理の性能を常に最適な状態に保つことが可能となり、また、上空風の特性が観測地によって異なっても、自動的にインコヒーレント積分時間を設定するため、設置時の信号処理設定が容易となるウィンドプロファイラにおける信号処理方法が得られるという効果がある。
【図面の簡単な説明】
【図１】この発明の実施の形態１のウィンドプロファイラにおける信号処理装置を示すブロック図である。
【図２】実施の形態１における風速ベクトル算出の原理を示す説明図である。
【図３】実施の形態１における大気乱流エコーの検出の原理を示す説明図である。
【図４】実施の形態１における信号処理装置で出力される速度データが高度によって異なる時間分解能を持つことを示す説明図である。
【図５】実施の形態１のウィンドプロファイラにおける信号処理方法を記すフローチャートである。
【図６】実施の形態１における信号処理装置の他の構成例を示すブロック図である。
【図７】この発明の実施の形態２のウィンドプロファイラにおける信号処理装置を示すブロック図。
【図８】実施の形態２のウィンドプロファイラにおける信号処理方法を示すフローチャートである。
【図９】実施の形態２における信号処理装置の他の構成例を示すブロック図である。
【図１０】この発明の実施の形態３のウィンドプロファイラにおける信号処理装置を示すブロック図である。
【図１１】実施の形態３における信号処理装置の他の構成例を示すブロック図である。
【図１２】この発明の実施の形態４におけるインコヒーレント積分の時間範囲の更新方法を示す説明図である。
【図１３】この発明および従来の信号処理装置および信号処理方法が適用される、一般的なウィンドプロファイラを示すブロック図である。
【図１４】従来のウィンドプロファイラにおける信号処理装置を示すブロック図である。
【符号の説明】
１ ＦＦＴ処理部、２ インコヒーレント積分部、３ ドップラ速度算出部、４ 品質管理・時間平均処理部、５ 積分時間設定部、６ データ補間処理部、７ データ取得率算出部、８ 高度毎時間平均数設定部、１０１ 空中線、１０２ 送受信装置、１０３ 信号処理装置、１０４ 風速ベクトル算出装置、１０５ 表示・記録装置。
Claims (11)
- 空間に電磁波を放射し、大気で反射された前記電磁波を受信して、そのドップラ周波数から風速を測定するウィンドプロファイラにおける信号処理装置において、
受信した前記電磁波による受信信号のデータにフーリエ変換を施すフーリエ変換処理部と、
前記フーリエ変換処理部でフーリエ変換された受信信号のデータからパワースペクトルを算出するとともに、複数の時刻で得られるパワースペクトルを、入力されたインコヒーレント積分時間だけ時間積分するインコヒーレント積分部と、
前記パワースペクトルを時間積分する際の最適なインコヒーレント積分時間を高度に応じて設定し、それを前記インコヒーレント積分部に入力する積分時間設定部と、
前記インコヒーレント積分部で出力されたパワースペクトルから信号スペクトルのピークを検出し、その中心ドップラ周波数から大気の視線方向速度を算出するドップラ速度算出部と、
前記ドップラ速度算出部から出力されたドップラ速度のうちの、低品質のデータを取り除くとともに時間平均を行う品質管理・時間平均処理部とを備えたことを特徴とするウィンドプロファイラにおける信号処理装置。 - 積分時間設定部が、低高度においては短いインコヒーレント積分時間を設定し、高高度においては長いインコヒーレント積分時間を設定するものであることを特徴とする請求項１記載のウィンドプロファイラにおける信号処理装置。
- ドップラ速度算出部から出力されるドップラ速度データを入力して、データ取得率を算出するデータ取得率算出部を有し、
積分時間設定部が、前記データ取得率算出部で算出されたデータ取得率に応じてインコヒーレント積分時間を設定するものであることを特徴とする請求項１記載のウィンドプロファイラにおける信号処理装置。 - 品質管理・時間平均処理部から出力されるドップラ速度データを入力して、データ取得率を算出するデータ取得率算出部を有し、
積分時間設定部が、前記データ取得率算出部で算出されたデータ取得率に応じてインコヒーレント積分時間を設定するものであることを特徴とする請求項１記載のウィンドプロファイラにおける信号処理装置。 - データ取得率算出部が各高度毎のデータ取得率を算出するものであり、
積分時間設定部が、前記データ取得率算出部にて算出されたデータ取得率をもとに、その高度のインコヒーレント積分時間の設定を行うものであることを特徴とする請求項３または請求項４記載のウィンドプロファイラにおける信号処理装置。 - データ取得率算出部があらかじめ分割された高度範囲毎のデータ取得率を算出するものであり、
積分時間設定部が、前記データ取得率算出部にて算出されたデータ取得率をもとに、あらかじめ分割された高度範囲毎のインコヒーレント積分時間の設定を行うものであることを特徴とする請求項３または請求項４記載のウィンドプロファイラにおける信号処理装置。 - ドップラ速度算出部から出力されたデータを入力し、インコヒーレント積分時間を長く設定した高度において時間方向にデータを補間することにより、高度によらず一定の時間分解能でドップラ速度を出力するデータ補間処理部を有し、
品質管理・時間平均処理部が、前記データ補間処理部から出力されるドップラ速度を入力するものであることを特徴とする請求項１、請求項３、あるいは請求項４のうちのいずれか１項記載のウィンドプロファイラにおける信号処理装置。 - 積分時間設定部から出力される各高度におけるインコヒーレント積分時間が入力され、そのインコヒーレント積分時間に反比例するデータ数を各高度の時間平均数として設定する高度毎時間平均数設定部を有し、
品質管理・時間平均処理部が、前記高度毎時間平均数設定部から出力される各高度の時間平均数を用いて時間平均の処理を行うことを特徴とする請求項１、請求項３、あるいは請求項４のうちのいずれか１項記載のウィンドプロファイラにおける信号処理装置。 - インコヒーレント積分部が、インコヒーレント積分を行う際に、時間方向の積分範囲の更新幅をインコヒーレント積分時間によらず一定にすることを特徴とする請求項１、請求項３、あるいは請求項４のうちのいずれか１項記載のウィンドプロファイラにおける信号処理装置。
- 空間に電磁波を放射し、大気で反射された前記電磁波を受信して、そのドップラ周波数から風速を測定するウィンドプロファイラにおける信号処理方法において、
受信信号のデータを入力するデータ読込みステップと、
前記データ読込みステップにて入力された受信信号のデータにフーリエ変換を施すフーリエ変換ステップと、
高度毎にインコヒーレント積分時間を設定する積分時間設定ステップと、
前記フーリエ変換ステップにてフーリエ変換された受信信号のデータからパワースペクトルを算出し、前記積分時間設定ステップにて設定されたインコヒーレント積分時間だけインコヒーレント積分を実施するインコヒーレント積分ステップと、
前記インコヒーレント積分ステップによるインコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出するドップラ速度算出ステップと、
処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを判定する全高度処理終了判定ステップと、
前記全高度処理終了判定ステップにて全高度データの処理が終了していないと判定された場合に、処理対象となる高度を次の高度へと更新し、前記積分時間設定ステップに処理を戻す高度更新ステップと、
前記全高度処理終了判定ステップにて全高度データの処理が終了したと判定された場合に、処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、前記データ読込みステップに処理を戻す時刻更新ステップを有することを特徴とするウィンドプロファイラにおける信号処理方法。 - 空間に電磁波を放射し、大気で反射された前記電磁波を受信して、そのドップラ周波数から風速を測定するウィンドプロファイラにおける信号処理方法において、
受信信号のデータを入力するデータ読込みステップと、
前記データ読込みステップにて入力された受信信号のデータにフーリエ変換を施すフーリエ変換ステップと、
インコヒーレント積分時間の初期値の設定を行う積分時間初期値設定ステップと、
前記フーリエ変換ステップにてフーリエ変換された受信信号のデータからパワースペクトルを算出し、前記積分時間設定ステップにて設定されたインコヒーレント積分時間だけインコヒーレント積分を実施するインコヒーレント積分ステップと、
前記インコヒーレント積分ステップによるインコヒーレント積分後のパワースペクトルから信号スペクトルを検出し、そのドップラ周波数からドップラ速度を算出するドップラ速度算出ステップと、
前記ドップラ速度算出ステップで算出されたドップラ速度から、現在処理対象となっている高度のデータ取得率を算出するデータ取得率算出ステップと、
前記データ取得率算出ステップで算出されたデータ取得率をあらかじめ定められたしきい値と比較するデータ取得率判定ステップと、
前記データ取得率判定ステップにてデータ取得率がしきい値未満と判定された場合に、現時点のインコヒーレント積分時間とあらかじめ定められたインコヒーレント積分時間とを比較する積分時間判定ステップと、
前記積分時間判定ステップにてインコヒーレント積分時間があらかじめ定められたインコヒーレント積分時間値未満と判定された場合に、インコヒーレント積分時間の設定を増加させて、処理をインコヒーレント積分ステップに戻す積分時間増加ステップと、
前記データ取得率判定ステップにてデータ取得率がしきい値以上と判定された場合、もしくは前記積分時間判定ステップにてインコヒーレント積分時間が上限値以上と判定された場合に、処理対象となっている時刻のデータについて、全高度のデータが処理されたか否かを判定する全高度処理終了判定ステップと、
前記全高度処理終了判定ステップにて全高度データの処理が終了していないと判定された場合に、処理対象となる高度を次の高度へと更新し、前記積分時間設定ステップに処理を戻す高度更新ステップと、
前記全高度処理終了判定ステップにて全高度データの処理が終了したと判定された場合に、処理対象となる時刻を次の時刻に進めるとともに、処理対象となる高度を初期化し、前記データ読込みステップに処理を戻す時刻更新ステップを有することを特徴とするウィンドプロファイラにおける信号処理方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000365808A JP3626089B2 (ja) | 2000-11-30 | 2000-11-30 | ウィンドプロファイラにおける信号処理装置および信号処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000365808A JP3626089B2 (ja) | 2000-11-30 | 2000-11-30 | ウィンドプロファイラにおける信号処理装置および信号処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002168948A JP2002168948A (ja) | 2002-06-14 |
| JP3626089B2 true JP3626089B2 (ja) | 2005-03-02 |
Family
ID=18836517
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000365808A Expired - Fee Related JP3626089B2 (ja) | 2000-11-30 | 2000-11-30 | ウィンドプロファイラにおける信号処理装置および信号処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3626089B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10234556B2 (en) | 2012-08-09 | 2019-03-19 | Mitsubishi Electric Corporation | Signal processing device and signal processing method in wind profiler |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4819343B2 (ja) * | 2004-11-09 | 2011-11-24 | Ｎｅｃエンジニアリング株式会社 | ウィンドプロファイラ・システム |
| JP4806949B2 (ja) * | 2005-03-31 | 2011-11-02 | 三菱電機株式会社 | レーザレーダ装置 |
| JP5111778B2 (ja) * | 2006-04-28 | 2013-01-09 | 株式会社東芝 | レーダ装置とその信号処理方法 |
| JP5046793B2 (ja) * | 2007-08-24 | 2012-10-10 | 三菱電機株式会社 | 風計測装置 |
| JP5605983B2 (ja) * | 2008-09-01 | 2014-10-15 | 三菱電機株式会社 | レーダ装置 |
| JP5781018B2 (ja) * | 2012-06-08 | 2015-09-16 | 三菱電機株式会社 | 風計測装置 |
| JPWO2014024508A1 (ja) | 2012-08-08 | 2016-07-25 | 三菱電機株式会社 | レーダ装置 |
| EP3296765B1 (en) * | 2015-05-12 | 2019-11-13 | Mitsubishi Electric Corporation | Laser radar device and wind speed observation method |
| CN109916852B (zh) * | 2019-01-31 | 2020-11-03 | 上海禾赛科技股份有限公司 | 一种激光气体遥测仪信号采集方法和系统 |
| CN109828276B (zh) * | 2019-03-06 | 2023-04-21 | 南京微麦科斯电子科技有限责任公司 | 一种基于干涉算法的风廓线雷达 |
-
2000
- 2000-11-30 JP JP2000365808A patent/JP3626089B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10234556B2 (en) | 2012-08-09 | 2019-03-19 | Mitsubishi Electric Corporation | Signal processing device and signal processing method in wind profiler |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002168948A (ja) | 2002-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110609287B (zh) | 一种双频雷达散射计及同时测量海面风场和流场的方法 | |
| CN102176000B (zh) | 一种船用雷达海杂波的抑制方法 | |
| JP3335544B2 (ja) | レーダ装置及びそのレーダ信号処理方法 | |
| JP3626089B2 (ja) | ウィンドプロファイラにおける信号処理装置および信号処理方法 | |
| Vicen-Bueno et al. | Real-time ocean wind vector retrieval from marine radar image sequences acquired at grazing angle | |
| US7132974B1 (en) | Methods and systems for estimating three dimensional distribution of turbulence intensity using radar measurements | |
| US9097805B2 (en) | Systems and methods for calibrating dual polarization radar systems | |
| JP3794361B2 (ja) | レーダ信号処理装置及びレーダ信号処理方法 | |
| US8456350B2 (en) | Methods and apparatus for sea state measurement via radar sea clutter eccentricity | |
| CN116106855B (zh) | 摔倒检测方法及摔倒检测装置 | |
| Haines et al. | Implementing quality control of high-frequency radar estimates and application to Gulf Stream surface currents | |
| WO2018168165A1 (ja) | 気象予測装置、気象予測方法、およびプログラム | |
| AU2010222857B2 (en) | Method of detecting and compensating for precipitation in sodar systems | |
| JP3613112B2 (ja) | レーダ装置及びその制御方法 | |
| Sujithlal et al. | Identification of the tropopause using the Jenks natural breaks classification from 205 MHz stratosphere–troposphere wind profiler radar | |
| CN113945988B (zh) | 基于雷达测量的覆冰检测方法、覆冰检测系统及设备 | |
| CN114384489A (zh) | 一种浴室跌倒检测雷达的水雾杂波滤除方法 | |
| Wang et al. | An improved asymmetric hurricane parametric model based on cross-polarization SAR observations | |
| CN105527624B (zh) | 一种雷达回波动态估算噪声的方法和气象雷达系统 | |
| JP3702347B2 (ja) | ウィンドプロファイラにおける信号処理方法及びそのプログラムと装置 | |
| CN117491998A (zh) | 一种步进频率合成孔径成像方法及系统 | |
| Furumoto et al. | Continuous humidity monitoring in a tropical region with the equatorial atmosphere radar | |
| US4143547A (en) | Method and apparatus for remote measurement of variation in wind velocity | |
| Anandan et al. | First results of experimental tests of the newly developed NARL phased-array Doppler Sodar | |
| CN116736242A (zh) | 一种降水测量雷达稳定性分析的数据处理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20041001 |
|
| 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: 20041102 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20041201 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071210 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081210 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091210 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091210 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101210 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111210 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111210 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121210 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121210 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131210 Year of fee payment: 9 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |