【０００１】
【発明の属する技術分野】
この発明は、ＣＡＰＰＩ観測等により観測された観測データを地図上に投影させるための等高度データに変換するデータ処理方法、及びこれを用いた気象レーダ装置に関する。
【０００２】
【従来の技術】
気象レーダは、空中の雲の様子や、降雨状態を検出する装置である。また、ＣＡＰＰＩ（Ｃｏｎｓｔａｎｔ Ａｌｔｉｔｕｄｅ Ｐｌａｎｅ Ｐｏｓｉｔｉｏｎ Ｉｎｄｉｃａｔｏｒ）観測は、空間の気象状況を立体的に捉えることができる気象レーダの運用形態の一つである。このＣＡＰＰＩ観測により得られた観測データは、さらに等高度変換というデータ処理が行われ、地上からの一定高度面における気象情報に加工される。
【０００３】
図５はＣＡＰＰＩ観測と等高度変換される観測データとの関係を示す概念図である。例えば、山頂等に設置された気象レーダ装置９は、所定ビーム幅のアンテナビーム１０をその仰角を変えて方位方向に走査する。例えば、全方位に走査することにより気象レーダ装置９の設置点を中心とした全方位の気象状況を観測することができる。そして、このようなアンテナビーム１０の走査により得られた観測データ（例えば、極座標系の観測データ）を地図上に投影させるための等高度データに変換する。この場合、等高度面Ｚの等高度データは、図５に示すようにアンテナビーム１０の各仰角毎に得られた観測データのうち等高度面Ｚに接する部分（網掛け部分）の観測データがそれぞれ等高度面Ｚ上の等高度データとして算出され、変換される。
【０００４】
従って、等高度面Ｚの高度が高高度になるとアンテナビーム１２を高仰角の位置で走査しなければならず、それに伴い処理ループ回数も増加する。そこで、従来の気象レーダ装置においては、処理時間短縮のため、あらかじめオフラインにて等高度変換テーブルを作成し、このような等高度変換テーブルを用いて等高度データを算出するデータ処理方法が利用されていた。
【０００５】
【発明が解決しようとする課題】
しかし、上述した従来のデータ処理方法は、あらかじめ作成したオフラインの等高度変換テーブルを用いて等高度データを算出するので、以下のような問題点があった。
（１）オフラインテーブルを確保するために余分なリソース（メモリ、ディスク領域など）を準備する必要がある。そのリソースは、変換高度数が多いほど、また方位分解能が細かいほど大量に必要である。
（２）各種の観測方式を適宜組み合わせて運用するような場合、テーブルのパターンをそれに伴い多数準備しておく必要がある。
（３）また、運用中にアンテナビーム１２の観測仰角の間隔を変更するような自由度の高い運用形態を行うことが出来ない。この場合、指示された実際の仰角値が等高度変換テーブルの仰角値と異なる場合でも、あらかじめ作成された仰角値のままで演算を行うため、誤差が生ずる。
【０００６】
この発明は、上記のような課題を解消するためになされたもので、ＣＡＰＰＩ観測等により観測された観測データを短い処理時間で等高度データに変換することができる一方、アンテナビームの仰角間隔の変更等にも拘らず、観測データの正確な等高度変換処理を行うことができる新規なデータ処理方法、およびこれを用いた気象レーダ装置を提供することを目的としている。
【０００７】
【課題を解決するための手段】
請求項１の発明に係るデータ処理方法は、空中線を仰角方向及び方位方向に走査し、前記空中線の走査を介して観測された極座標系の観測データを等高度座標系のデータに変換するデータ処理方法において、等高度面を求めるステップと、前記等高度面上の所定位置と前記空中線の設置点との間の水平距離を求めるステップと、前記等高度面までの海抜高度と前記水平距離とに基づいて前記所定位置に対応する前記空中線の仮想仰角を求めるステップと、前記空中線による現実の観測仰角のうち前記仮想仰角を挟む２つの観測仰角を求めるステップと、これら２つの観測仰角と前記水平距離までの前記２つの観測仰角における仰角方向距離とから前記等高度面を挟む前記２つの観測仰角に対応した高度を求めるステップと、前記ステップにより求めた観測仰角、仰角方向距離及び高度からなる２組の極座標系の観測データから近似により前記等高度面上の所定位置における極座標系の観測データを求めるステップと、前記近似により求めた前記所定位置の極座標系の観測データを前記等高度面における等高度座標系のデータに変換するステップとを備えたものである。
【０００８】
請求項２の発明に係るデータ処理方法は、前記観測データが極座標系の観測データであるものである。
【０００９】
請求項３の発明に係るデータ処理方法は、前記等高度面上の所定位置における観測データの近似が直線補完方式又は最大値方式を用いたものである。
【００１０】
請求項４の発明に係る気象レーダ装置は、送信部と、この送信部から出力された送信信号を送信電波として大気中に放射し、その反射電波を受信する空中線と、この空中線に観測方式に応じた方位及び仰角方向を指示する空中線制御部と、前記空中線により受信された前記反射電波をビデオ信号に変換する受信部と、この受信部から出力されたビデオ信号と前記空中線の方位及び仰角方向とに基づいて極座標系の観測データを出力する信号処理部と、この信号処理部により得られた前記極座標系の観測データを等高度座標系のデータに変換する等高度変換部とを備え、この等高度変換部は、等高度面の所定位置に対応する前記空中線の仮想仰角を求め、この仮想仰角を挟む２つの観測仰角に対応した２組の観測データから等高度座標系のデータを算出するものである。
【００１１】
【発明の実施の形態】
実施の形態１．
以下、この発明の一実施形態について図１乃至図４を用いて説明する。図１はこの発明の実施形態１による気象レーダ装置を示すブロック構成図である。図１において、１は送信信号を出力する送信部、２は送信部１から出力された送信信号を送信電波として大気中に放射し、その反射電波を受信する空中線、３は空中線２に観測方式に応じた方位及び仰角方向を指示する命令信号を出力し、空中線から送信される送信電波の方位及び仰角方向を制御する空中線制御部、４は空中線２により受信された反射電波を受信処理し、ビデオ信号に変換する受信部、５は受信部４により変換されたビデオ信号を受信強度値（ｄＢｍ値）あるいは反射強度値（ｄＢＺ値）（以下、単に受信強度データという。）に変換する信号処理部、６は信号処理部５により変換された受信強度データを等高度座標系の気象データに変換する等高度変換部である。
【００１２】
なお、７は現実の観測データの中から等高度面Ｚにおける仮想仰角φｖを挟む観測仰角φ１，φ２を求め、この現実の観測仰角φ１，φ２から等高度座標系の気象データを求めるための２組の観測データを算出する仰角算出部、８は仰角算出部７により算出された２組の観測データに基づいて等高度座標系の気象データを求める変換処理部８であり、図１に示すように、実施の形態１による気象レーダ装置には、等高度変換処理のためのオフラインのデータテーブル（以下、オフラインテーブルという。）が予め設けられていない。従って、気象レーダ装置の運用パターンが変更された場合には、その都度変更された運用パターンによる現実の観測が行われ、実際に観測された観測データに基づいて等高度座標系の気象データが算出される。
【００１３】
また、信号処理部５から出力された受信強度値（ｄＢｍ値）あるいは反射強度値（ｄＢＺ値）等の受信強度データは、信号処理された時点では、極座標系（ｒ，θ，φ）の観測データであり、この極座標系の観測データが等高度変換部６の処理により等高度座標系（Ｒ，θ，Ｚ）のデータに変換される。
【００１４】
次に、図１に示す気象レーダ装置の全体動作について簡単に説明する。まず、空中線２が空中線制御部３により制御され、図５に示すように空中線２上に形成されたアンテナビーム（以下、観測ビームという。）が走査される。送信部１から出力された送信信号はこの観測ビームの方向に向かって送信電波として放射される。空中線２から放射された送信電波は大気中の雨粒、塵等によって反射され、その反射電波が空中線２により受信される。空中線２により受信された反射電波は受信部４においてビデオ信号に変換され、信号処理部５において極座標系の受信強度データに変換される。信号処理部５から出力された極座標系の受信強度データは、等高度変換部６に入力され、等高度座標系の気象データに変換される。なお、実施の形態１による気象レーダ装置はＣＡＰＰＩ観測を行うものとし、空中線制御部３はＣＡＰＰＩ観測方式に応じた方位及び仰角方向を指示する命令信号を空中線２に対して出力する。これにより、気象レーダ装置の設置点を中心とした３次元の観測データを得ることができる。
【００１５】
次に、等高度変換部６の変換処理について図２乃至図４を用いて説明する。図２は等高度変換部６の処理内容を示すフローチャート図である。なお、信号処理部５から出力された観測データの座標系を（ｒ，θ，φ）、求める等高度面Ｚの座標系（Ｒ，θ，Ｚ）と定義する。ここで、ｒは観測ビームにおけるレーダからの距離、θは観測ビームの方位角、φは観測ビームの仰角、Ｒは気象レーダ装置の設置点であるレーダ中心からの水平距離、Ｚは等高度面の高度である。図３はこれらの関係を示す座標系の模式図である。
【００１６】
図２に示すように、等高度変換部６は、まず、求めたい高度Ｈの等高度面Ｚを決め、この等高度面Ｚについて以下の処理を行う（Ｓ０１）。すなわち、等高度面Ｚが決まると、次にレーダ中心から等高度面Ｚ上の求めたい位置までの水平距離Ｒを決め（Ｓ０２）、式（１）によりこの水平距離Ｒの位置における仮想仰角φｖを算出する（Ｓ０３）。
【００１７】
【数１】
【００１８】
仰角算出部７は、式（１）により仮想仰角φｖが求まると、実際に観測した観測仰角の中からこの仮想仰角φｖを挟む２つの現実の観測仰角（φ１，φ２）を抽出し、これら２つの観測仰角φ１，φ２の水平距離Ｒの位置までの仰角距離ｒ１，ｒ２、及びそれに対応した高度Ｈ１，Ｈ２を式（２）、式（３）によりそれぞれ算出する（Ｓ０４）。なお、一般に空中線の仰角φ、当該空中線から放射される観測ビームの仰角方向距離ｒ、及び海抜高度Ｈとの間には、式（３）に示す関係があることが幾何学的に求められている。Ｈ０は空中線２の地上海抜高度である。
【００１９】
【数２】
【００２０】
【数３】
【００２１】
これら２組の観測データ（φ１，ｒ１，Ｈ１）及び（φ２，ｒ２，Ｈ２）は、等高度座標系の気象データを求めるため、変換処理部８の記憶手段（図示省略）に一時的に記憶される（Ｓ０５）。そして、変換処理部８は、これら２組の観測データ（φ１，ｒ１，Ｈ１）及び（φ２，ｒ２，Ｈ２）から近似方式（近い方の点を採用）、最大値方式（強度値の大きい方を採用）、距離による直線補完法等の近似式を用いて等高度面Ｚ上の水平距離Ｒの位置における等高度データを算出する（Ｓ０５）。なお、図４は仮想仰角φｖに基づき抽出された２つの観測仰角φ１，φ２と等高度面Ｚとの関係を示す模式図である。
【００２２】
このようにして、等高度面Ｚ上のある方位の水平距離Ｒの点における等高度データ、すなわち等高度座標系のデータが算出され、極座標系の観測データから等高度座標系のデータへの変換処理が行われる（Ｓ０６）。なお、水平距離Ｒを変更して以上の処理を繰り返すことにより、等高度面Ｚ上のある一方位についての等高度データを算出することができ（Ｓ０７）、また、以上の処理を方位を変更して繰り返すことにより、等高度面Ｚ上の全方位についての等高度データを算出することができる（Ｓ０８）。
【００２３】
ここで、式（１）の導出を示す。式（２）を式（３）に代入すると、式（４）、さらに式（５）となる。ゆえに、これから等高度面Ｚ上の所定位置における空中線２の仮想仰角φｖを求める式（３）が導出される。
【００２４】
【数４】
【００２５】
【数５】
【００２６】
以上のように、この実施の形態１によるデータ処理方法によれば、等高度面上の所定位置における空中線２の仮想仰角φｖを求め、この仮想仰角φｖを挟む２つの観測仰角φ１，φ２を抽出し、これら２つの観測仰角φ１，φ２に対応した観測仰角、仰角方向距離及び高度からなる２組の観測データから近似により前記等高度面上の所定位置における観測データを求め、前記近似により求めた前記所定位置の観測データを前記等高度面における等高度座標系のデータに変換するようにしたので、オフラインテーブルを確保するための余分なリソース（メモリ、ディスク領域など）を準備することなくＣＡＰＰＩ観測等により観測された極座標系の観測データを短時間で等高度データに変換することができる一方、アンテナビームの仰角間隔の変更等にも拘らず、観測データの正確な等高度変換処理を行うことができる。
【００２７】
実施の形態２．
次にこの発明の実施の形態２について説明する。実施形態１による気象レーダ装置では、極座標系の観測データ（（ｒ，θ，φ）−＞（Ｒ，θ，Ｚ））を等高度座標系のデータに変換する処理を行ったが、直交座標系（（Ｘ，Ｙ，φ）−＞（Ｘ，Ｙ，Ｚ））の観測データから等高度座標系のデータに変換することも可能である。すなわち、仰角別の観測データをあらかじめ直交座標に変換した後に、等高度面に変換することにより、等高度座標系のデータを得ることができる。
【００２８】
この場合、メッシュ情報そのものにはレーダ中心からの距離ｒ及び地上面水平距離Ｒが保持されていないため、処理の都度、空中線２の設置点、すなわちレーダ中心座標位置からｒ及びＲを求める必要がある。地上面水平距離Ｒは、レーダ中心座標をＸ０、Ｙ０、求めたい座標位置をＸ、Ｙとすると、式（６）により求めることができる。
【００２９】
【数６】
【００３０】
そして、極座標系の観測データの場合と同様の演算を行うことにより、等高度座標系のデータに変換することができる。このように、直交座標系のデータの場合には、求めたい座標位置毎にＲが異なること・各高度面で共通の値であること、観測毎に変動する要素ではないことから、オフラインにてあらかじめ計算しておく、あるいは一度計算した値を保持するようにしておくなどの効率化を図ることは有意義である。
【００３１】
【発明の効果】
以上のように、この発明によれば、ＣＡＰＰＩ観測等により観測された観測データを仮想仰角を用いたデータ処理によって等高度データに変換するので、オフラインテーブルを確保するための余分なリソース（メモリ、ディスク領域など）を準備することなく短時間で等高度データに変換することができる。また、アンテナビームの仰角間隔の変更等にも拘らず、観測データの正確な等高度変換処理を行うことができる。
【図面の簡単な説明】
【図１】実施の形態１による気象レーダ装置を示すブロック構成図である。
【図２】実施の形態１によるデータ処理方法であって、図１に示す等高度変換部６の処理内容を示すフローチャート図である。
【図３】実施の形態１による座標系の模式図である。
【図４】仮想仰角φｖに基づき抽出された２つの観測仰角φ１，φ２と等高度面Ｚとの関係を示す模式図である。
【図５】一般的なＣＡＰＰＩ観測と等高度変換される観測データとの関係を示す概念図である。
【符号の説明】
１ 空中線部、２ 空中線制御部、３ 送信部、４ 受信部、
５ 信号処理部、６ 等高度変換部、７ 仰角算出部、８ 変換処理部。[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a data processing method for converting observation data observed by CAPPI observation or the like into iso-altitude data for projecting on a map, and a weather radar device using the same.
[0002]
[Prior art]
A weather radar is a device that detects the state of clouds in the air and the state of rainfall. Further, CAPPI (Constant Altitude Plane Position Indicator) observation is one of the operation modes of a weather radar capable of three-dimensionally capturing the weather conditions in space. Observation data obtained by the CAPPI observation is further subjected to data processing called iso-altitude conversion, and processed into weather information at a constant altitude plane from the ground.
[0003]
FIG. 5 is a conceptual diagram showing the relationship between CAPPI observation and observation data subjected to iso-altitude conversion. For example, the weather radar device 9 installed on a mountaintop or the like scans the antenna beam 10 having a predetermined beam width in the azimuth direction while changing its elevation angle. For example, by scanning in all directions, it is possible to observe the weather conditions in all directions centering on the installation point of the weather radar device 9. Then, the observation data (for example, observation data in the polar coordinate system) obtained by scanning the antenna beam 10 is converted into iso-altitude data for projecting on a map. In this case, the iso-altitude data of the iso-altitude plane Z is, as shown in FIG. 5, the observation data of the portion (shaded area) in contact with the iso-altitude plane Z among the observation data obtained for each elevation angle of the antenna beam 10. Each is calculated and converted as iso-altitude data on the iso-altitude plane Z.
[0004]
Therefore, when the height of the iso-altitude plane Z becomes higher, the antenna beam 12 must be scanned at a position with a high elevation angle, and the number of processing loops increases accordingly. Therefore, in a conventional weather radar device, a data processing method of creating an iso-altitude conversion table in advance and calculating iso-altitude data using such an iso-altitude conversion table is used in order to reduce processing time. I was
[0005]
[Problems to be solved by the invention]
However, the above-described conventional data processing method calculates iso-altitude data using an offline iso-altitude conversion table created in advance, and thus has the following problems.
(1) Extra resources (memory, disk area, etc.) need to be prepared to secure an offline table. The resources are required in large amounts as the number of converted altitudes is large and the azimuth resolution is fine.
(2) In the case where various observation methods are used in combination as appropriate, it is necessary to prepare a large number of table patterns accordingly.
(3) In addition, an operation mode with a high degree of freedom, such as changing the interval between the observation elevation angles of the antenna beam 12 during operation, cannot be performed. In this case, even when the specified actual elevation angle value is different from the elevation angle value in the iso-altitude conversion table, the calculation is performed with the elevation angle value created in advance, so that an error occurs.
[0006]
The present invention has been made in order to solve the above-described problems, and can convert observation data observed by CAPPI observation or the like into iso-altitude data in a short processing time, while maintaining the elevation angle interval of an antenna beam. It is an object of the present invention to provide a new data processing method capable of performing accurate iso-altitude conversion processing of observation data despite a change or the like, and a weather radar apparatus using the same.
[0007]
[Means for Solving the Problems]
The data processing method according to claim 1 is a data processing for scanning an aerial in elevation and azimuth directions, and converting observation data of a polar coordinate system observed through the aerial scanning into data of an iso-altitude coordinate system. Determining a contour level, determining a horizontal distance between a predetermined position on the contour plane and the installation point of the aerial, and calculating a sea level altitude to the contour plane and the horizontal distance. Calculating a virtual elevation angle of the antenna corresponding to the predetermined position based on the predetermined position; obtaining two observation elevation angles sandwiching the virtual elevation angle among actual observation elevation angles by the antenna; and calculating the two observation elevation angles and the horizontal distance. Obtaining an altitude corresponding to the two observation elevation angles sandwiching the iso-altitude plane from the elevation angle distances at the two observation elevation angles up to; Obtaining the observation data of the polar coordinate system at the predetermined position on the iso-altitude plane by approximation from the observation data of the two sets of polar coordinate systems consisting of the obtained observation elevation angle, elevation direction distance, and altitude; and the predetermined position obtained by the approximation. Converting the observation data in the polar coordinate system into data in the iso-altitude coordinate system on the iso-altitude plane.
[0008]
In a data processing method according to a second aspect of the present invention, the observation data is observation data in a polar coordinate system.
[0009]
In a data processing method according to a third aspect of the present invention, the approximation of the observation data at the predetermined position on the iso-altitude plane uses a linear interpolation method or a maximum value method.
[0010]
A meteorological radar device according to a fourth aspect of the present invention provides a transmitting unit, an antenna for transmitting a transmission signal output from the transmitting unit to the atmosphere as a transmission radio wave, and receiving the reflected radio wave, and an observation system for the antenna. An antenna control unit for instructing the corresponding azimuth and elevation direction, a reception unit for converting the reflected radio wave received by the antenna into a video signal, a video signal output from the reception unit and the azimuth and elevation direction of the antenna A signal processing unit that outputs the observation data of the polar coordinate system based on and an iso-altitude conversion unit that converts the observation data of the polar coordinate system obtained by the signal processing unit to data of the iso-altitude coordinate system, The iso-altitude conversion unit obtains a virtual elevation angle of the antenna corresponding to a predetermined position on the iso-altitude plane, and obtains data of an iso-altitude coordinate system from two sets of observation data corresponding to two observation elevation angles sandwiching the virtual elevation angle. And it calculates the.
[0011]
BEST MODE FOR CARRYING OUT THE INVENTION
Embodiment 1 FIG.
Hereinafter, an embodiment of the present invention will be described with reference to FIGS. FIG. 1 is a block diagram showing a weather radar device according to Embodiment 1 of the present invention. In FIG. 1, reference numeral 1 denotes a transmitting unit that outputs a transmission signal, 2 denotes an antenna that emits the transmission signal output from the transmitting unit 1 into the atmosphere as a transmission radio wave, and an antenna that receives the reflected radio wave, and 3 denotes an observation system to the antenna 2. An antenna control unit that outputs a command signal instructing the azimuth and the elevation angle direction according to, and controls the azimuth and the elevation direction of the transmission radio wave transmitted from the antenna, receives the reflected radio wave received by the antenna 2, and performs a reception process. The receiver 5 converts the video signal into a video signal, and the signal processing converts the video signal converted by the receiver 4 into a reception intensity value (dBm value) or a reflection intensity value (dBZ value) (hereinafter, simply referred to as reception intensity data). Units 6 are iso-altitude conversion units for converting the received intensity data converted by the signal processing unit 5 into weather data in an iso-altitude coordinate system.
[0012]
Reference numeral 7 denotes an observation elevation angle φ1, φ2 sandwiching the virtual elevation angle φv on the iso-altitude plane Z from the actual observation data, and 2 for obtaining weather data in the iso-altitude coordinate system from the actual observation elevation angles φ1, φ2. An elevation angle calculation unit 8 for calculating a set of observation data is a conversion processing unit 8 for obtaining weather data in an iso-altitude coordinate system based on the two sets of observation data calculated by the elevation angle calculation unit 7, as shown in FIG. In addition, the weather radar device according to the first embodiment is not provided with an offline data table (hereinafter, referred to as an offline table) for iso-altitude conversion processing in advance. Therefore, when the operation pattern of the weather radar device is changed, actual observations are performed based on the changed operation pattern each time, and weather data in the iso-altitude coordinate system is calculated based on the actually observed observation data. Is done.
[0013]
The reception intensity data such as the reception intensity value (dBm value) or the reflection intensity value (dBZ value) output from the signal processing unit 5 is obtained by observing the polar coordinate system (r, θ, φ) at the time of signal processing. The observation data of the polar coordinate system is converted into data of the iso-altitude coordinate system (R, θ, Z) by the processing of the iso-altitude conversion unit 6.
[0014]
Next, the overall operation of the weather radar device shown in FIG. 1 will be briefly described. First, the antenna 2 is controlled by the antenna control unit 3, and an antenna beam (hereinafter, referred to as an observation beam) formed on the antenna 2 is scanned as shown in FIG. The transmission signal output from the transmission unit 1 is radiated as a transmission radio wave in the direction of the observation beam. The transmission radio wave radiated from the antenna 2 is reflected by raindrops, dust and the like in the atmosphere, and the reflected radio wave is received by the antenna 2. The reflected radio wave received by the antenna 2 is converted into a video signal in the receiving unit 4 and converted into reception intensity data in a polar coordinate system in the signal processing unit 5. The reception intensity data of the polar coordinate system output from the signal processing unit 5 is input to the iso-altitude conversion unit 6 and converted into weather data of the iso-altitude coordinate system. The weather radar device according to the first embodiment performs CAPPI observation, and the antenna control unit 3 outputs to the antenna 2 a command signal indicating an azimuth and an elevation direction according to the CAPPI observation method. Thereby, three-dimensional observation data centering on the installation point of the weather radar device can be obtained.
[0015]
Next, the conversion process of the equal height conversion unit 6 will be described with reference to FIGS. FIG. 2 is a flowchart showing the processing contents of the equal height conversion unit 6. The coordinate system of the observation data output from the signal processing unit 5 is defined as (r, θ, φ) and the coordinate system (R, θ, Z) of the iso-altitude plane Z to be obtained. Here, r is the distance of the observation beam from the radar, θ is the azimuth of the observation beam, φ is the elevation angle of the observation beam, R is the horizontal distance from the radar center where the weather radar device is installed, and Z is the iso-altitude plane. Of altitude. FIG. 3 is a schematic diagram of a coordinate system showing these relationships.
[0016]
As shown in FIG. 2, the iso-altitude conversion section 6 first determines the iso-altitude plane Z of the altitude H to be obtained, and performs the following processing on the iso-altitude plane Z (S01). That is, when the iso-altitude plane Z is determined, the horizontal distance R from the radar center to the desired position on the iso-altitude plane Z is determined (S02), and the virtual elevation angle φv at the position of the horizontal distance R is calculated by Expression (1). Is calculated (S03).
[0017]
(Equation 1)
[0018]
When the virtual elevation angle φv is obtained from Expression (1), the elevation angle calculation unit 7 extracts two actual observation elevation angles (φ1, φ2) sandwiching the virtual elevation angle φv from the observation elevation angles actually observed. The elevation angles r1 and r2 to the position of the horizontal distance R between the two observation elevation angles φ1 and φ2, and the altitudes H1 and H2 corresponding to the elevation angles r1 and r2 are calculated by equations (2) and (3), respectively (S04). In general, it is geometrically determined that there is a relationship shown in Expression (3) between the elevation angle φ of the antenna, the elevation distance r of the observation beam emitted from the antenna, and the altitude H above sea level. I have. H 0 is the altitude above Antenna 2 above Shanghai.
[0019]
(Equation 2)
[0020]
[Equation 3]
[0021]
These two sets of observation data (φ1, r1, H1) and (φ2, r2, H2) are temporarily stored in storage means (not shown) of the conversion processing unit 8 in order to obtain meteorological data in an iso-altitude coordinate system. Is performed (S05). Then, the conversion processing unit 8 uses the two sets of observation data (φ1, r1, H1) and (φ2, r2, H2) to calculate the approximation method (using the closest point) and the maximum value method (the one with the larger intensity value). Is used, and the iso-altitude data at the position of the horizontal distance R on the iso-altitude plane Z is calculated using an approximate expression such as a straight line interpolation method based on distance (S05). FIG. 4 is a schematic diagram showing the relationship between two observation elevation angles φ1 and φ2 extracted based on the virtual elevation angle φv and the iso-altitude plane Z.
[0022]
In this manner, iso-altitude data at a point at a horizontal distance R in a certain direction on the iso-altitude plane Z, that is, data in the iso-altitude coordinate system is calculated, and conversion from observation data in the polar coordinate system to data in the iso-altitude coordinate system is performed. Processing is performed (S06). By changing the horizontal distance R and repeating the above processing, it is possible to calculate the iso-altitude data for a certain position on the iso-altitude plane Z (S07). By repeating the above, iso-altitude data for all directions on the iso-altitude plane Z can be calculated (S08).
[0023]
Here, the derivation of Expression (1) will be described. When equation (2) is substituted into equation (3), equation (4) and further equation (5) are obtained. Therefore, the equation (3) for obtaining the virtual elevation angle φv of the antenna 2 at a predetermined position on the iso-altitude plane Z is derived from this.
[0024]
(Equation 4)
[0025]
(Equation 5)
[0026]
As described above, according to the data processing method according to the first embodiment, the virtual elevation angle φv of the antenna 2 at a predetermined position on the iso-altitude plane is obtained, and two observation elevation angles φ1 and φ2 sandwiching the virtual elevation angle φv are extracted. Observation data at a predetermined position on the iso-altitude plane is obtained by approximation from two sets of observation data consisting of the observation elevation angles corresponding to these two observation elevation angles φ1 and φ2, the elevation direction distance, and the altitude, and obtained by the approximation. Since the observation data at the predetermined position is converted into data in the iso-altitude coordinate system on the iso-altitude plane, CAPPI observation can be performed without preparing extra resources (memory, disk area, etc.) for securing an offline table. Can be converted to iso-altitude data in a short period of time, while changing the elevation angle interval of the antenna beam. It can be performed irrespective the exact equal high conversion of the observed data to.
[0027]
Embodiment 2 FIG.
Next, a second embodiment of the present invention will be described. In the meteorological radar device according to the first embodiment, the process of converting observation data ((r, θ, φ) → (R, θ, Z)) in the polar coordinate system into data in the iso-altitude coordinate system is performed. It is also possible to convert observation data of the system ((X, Y, φ) → (X, Y, Z)) into data of an iso-altitude coordinate system. That is, the data of the iso-altitude coordinate system can be obtained by converting the observation data for each elevation angle into orthogonal coordinates in advance, and then converting the data to the iso-altitude plane.
[0028]
In this case, since the mesh information itself does not hold the distance r from the radar center and the ground surface horizontal distance R, it is necessary to obtain r and R from the installation point of the antenna 2, that is, the radar center coordinate position each time the processing is performed. is there. Assuming that the radar center coordinates are X 0 and Y 0 and the desired coordinate positions are X and Y, the ground surface horizontal distance R can be obtained by Expression (6).
[0029]
(Equation 6)
[0030]
Then, by performing the same calculation as that of the observation data in the polar coordinate system, the data can be converted into data in the iso-altitude coordinate system. As described above, in the case of the data of the orthogonal coordinate system, R is different for each coordinate position to be obtained. It is significant to improve the efficiency by calculating in advance or by holding the value once calculated.
[0031]
【The invention's effect】
As described above, according to the present invention, observation data observed by CAPPI observation or the like is converted into iso-altitude data by data processing using a virtual elevation angle, so that extra resources (memory, memory, ) Can be converted to iso-altitude data in a short time without preparing a disk area. In addition, accurate iso-altitude conversion processing of observation data can be performed irrespective of a change in an elevation angle interval between antenna beams.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a weather radar device according to a first embodiment.
FIG. 2 is a flowchart showing a data processing method according to the first embodiment, showing processing contents of an iso-altitude conversion unit 6 shown in FIG. 1;
FIG. 3 is a schematic diagram of a coordinate system according to the first embodiment.
FIG. 4 is a schematic diagram showing a relationship between two observation elevation angles φ1 and φ2 extracted based on a virtual elevation angle φv and an iso-altitude plane Z;
FIG. 5 is a conceptual diagram showing a relationship between general CAPPI observation and observation data subjected to iso-altitude conversion.
[Explanation of symbols]
1 antenna unit, 2 antenna control unit, 3 transmission unit, 4 reception unit,
5 signal processing unit, 6 iso-altitude conversion unit, 7 elevation angle calculation unit, 8 conversion processing unit.
