JP2011033528A - Method and system of automatically tracking weather phenomenon - Google Patents
Method and system of automatically tracking weather phenomenon Download PDFInfo
- Publication number
- JP2011033528A JP2011033528A JP2009181539A JP2009181539A JP2011033528A JP 2011033528 A JP2011033528 A JP 2011033528A JP 2009181539 A JP2009181539 A JP 2009181539A JP 2009181539 A JP2009181539 A JP 2009181539A JP 2011033528 A JP2011033528 A JP 2011033528A
- Authority
- JP
- Japan
- Prior art keywords
- echo
- observation
- rhi
- weather
- azimuth angle
- 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.)
- Pending
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)
Abstract
Description
この発明は、気象レーダが雷雲等の気象現象を観測する際に、該気象現象をRHI観測により追尾するための自動追尾方法及び自動追尾システムに関する。 The present invention relates to an automatic tracking method and an automatic tracking system for tracking a weather phenomenon such as a thundercloud by RHI observation when the weather radar observes the weather phenomenon.
特許文献1〜3には、気象レーダがPPI(Plane Position Indicator)観測を行うことにより雷雲等の気象現象を追尾することが提案されている。
ところで、近年、大都市におけるヒートアイランド現象によって、夏季に局地的な積乱雲等の雷雲(気象現象)が頻発し、この結果、前記雷雲に起因した突発的な豪雨による市街地での浸水等の局地的な自然災害が増加している。前記積乱雲は、局地的且つ鉛直方向に成長する気象現象であるため、仰角を固定し且つ方位角を可変にした状態でレーダビームを方位角方向に走査する特許文献1〜3のPPI観測では、前記雷雲を追尾することは難しい。そこで、仰角を可変し且つ方位角を固定した状態でレーダビームを仰角方向に走査するRHI(Range Height Indicator)観測により前記雷雲を追尾することが考えられるが、気象レーダは、予め設定された観測スケジュールに従って、固定の方位角でRHI観測を行うため、前記積乱雲を追尾することは難しい。一方、PPI観測の結果に基づいて、ユーザが前記観測スケジュール中の方位角を変更し、前記気象レーダが前記変更した方位角にてRHI観測を行うことも考えられるが、前記観測スケジュールの変更作業中に、変更後の方位角とは異なる方位角の方向に前記雷雲が移動すれば、変更後の前記方位角でRHI観測を行っても前記雷雲を追尾することはできない。
By the way, in recent years, thunderclouds (meteorological phenomena) such as local cumulonimbus clouds frequently occur in summer due to the heat island phenomenon in large cities, and as a result, local flooding in urban areas due to sudden heavy rain caused by the thunderclouds. Natural disasters are increasing. Since the cumulonimbus is a local and vertical weather phenomenon, the PPI observation in
この発明は、上述した問題を解決するためになされたものであり、積乱雲を含む雷雲等の気象現象をRHI観測により確実且つ自動的に追尾することができる自動追尾方法及び自動追尾システムを提供することを目的とする。 The present invention has been made to solve the above-described problems, and provides an automatic tracking method and an automatic tracking system that can reliably and automatically track meteorological phenomena such as thunderclouds including cumulonimbus clouds by RHI observation. For the purpose.
この発明に係る気象現象の自動追尾方法は、
気象レーダが観測した気象現象を示すエコー強度データを用いてエコー頂高度及び鉛直積算雨水量を求める第1工程と、
求めた前記エコー頂高度及び前記鉛直積算雨水量に基づいて、前記エコー強度データ中の所定のエコーを、追尾すべき注目エコーとして決定する第2工程と、
前記注目エコーに応じた気象現象に対するRHI観測を前記気象レーダが行う際の方位角を算出する第3工程とを有することを特徴とする。
The automatic tracking method of the meteorological phenomenon according to the present invention,
A first step of determining an echo top height and a vertical accumulated rainwater amount using echo intensity data indicating a meteorological phenomenon observed by a weather radar;
A second step of determining a predetermined echo in the echo intensity data as a target echo to be tracked based on the obtained echo peak altitude and the vertically integrated rainwater amount;
And a third step of calculating an azimuth angle when the weather radar performs RHI observation for a weather phenomenon corresponding to the attention echo.
この発明に係る気象現象の自動追尾システムは、
観測シーケンスに従ってPPI観測を行うことにより気象現象を示すエコー強度データを取得する気象レーダと、
前記エコー強度データを用いてエコー頂高度及び鉛直積算雨水量を求め、求めた前記エコー頂高度及び前記鉛直積算雨水量に基づいて、前記エコー強度データ中の所定のエコーを、追尾すべき注目エコーとして決定し、決定した該注目エコーに応じた気象現象に対するRHI観測を前記気象レーダが行う際の方位角を算出し、算出した前記方位角を示す方位書き替えコマンドを発行するコマンド発行手段とを有し、
前記気象レーダは、前記コマンド発行手段が前記方位書き替えコマンドを発行した際に、前記注目エコーに応じた気象現象に対するRHI観測の方位角を、前記観測シーケンスに予め設定された方位角から、前記方位書き替えコマンドの方位角に変更して、前記気象現象に対するRHI観測を行うことを特徴とする。
An automatic tracking system for meteorological phenomena according to the present invention includes:
A weather radar for acquiring echo intensity data indicating a weather phenomenon by performing PPI observation according to an observation sequence;
Using the echo intensity data, the echo top height and the vertical accumulated rainwater amount are obtained, and based on the obtained echo top altitude and the vertical accumulated rainwater amount, a predetermined echo in the echo intensity data to be tracked to be tracked Command issuing means for calculating an azimuth angle when the weather radar performs RHI observation for the meteorological phenomenon corresponding to the determined echo of interest, and issuing an azimuth rewriting command indicating the calculated azimuth angle Have
When the command issuing means issues the direction rewriting command, the weather radar determines an azimuth angle of RHI observation for a meteorological phenomenon corresponding to the echo of interest from an azimuth angle set in advance in the observation sequence. The azimuth angle of the azimuth rewriting command is changed to perform RHI observation for the meteorological phenomenon.
これらの発明によれば、気象レーダが観測した気象現象を示すエコー強度データより得られたエコー頂高度及び鉛直積算雨水量を用いて、所定のエコーを注目エコーとして自動的に決定し、さらに、決定した前記注目エコーに応じた気象現象に対してRHI観測を行う際の方位角を自動的に算出するので、前記注目エコーに応じた積乱雲を含む雷雲等の気象現象を確実且つ自動的に追尾することができる。 According to these inventions, a predetermined echo is automatically determined as an attention echo using an echo peak height and a vertical accumulated rainwater amount obtained from echo intensity data indicating a weather phenomenon observed by a weather radar, Since the azimuth angle at the time of RHI observation is automatically calculated for the meteorological phenomenon corresponding to the determined echo of interest, the meteorological phenomenon such as thunderclouds including the cumulonimbus corresponding to the echo of interest is reliably and automatically tracked. can do.
この発明に係る気象現象の自動追尾方法及び自動追尾システムの好適な実施形態について、図1〜図11を参照しながら説明する。 A preferred embodiment of a weather phenomenon automatic tracking method and automatic tracking system according to the present invention will be described with reference to FIGS.
本実施形態に係る自動追尾システム10は、図1に示すように、アンテナ12、アンテナ駆動部14、送受信機16、信号処理機18、ホストコンピュータ(以下、ホストPCともいう。)20及びデータベース24より構成される気象ドップラレーダ23と、コマンド発行手段としてのオンライン化解析用コンピュータ(以下、解析用PCともいう。)22とを有する。
As shown in FIG. 1, the
アンテナ12は、送受信機16からアンテナ駆動部14を介して出力された送信信号をレーダビームとして外部に放射し、一方で、積乱雲等の雷雲(以下、気象現象ともいう。)で反射した前記レーダビーム(反射波)を受信した際に、該反射波を受信信号としてアンテナ駆動部14を介して送受信機16に出力する。この場合、気象ドップラレーダ23は、アンテナ駆動部14によるアンテナ12の回転によって、PPI観測(仰角を固定した状態で方位角方向にアンテナ12を回転することによりレーダビームの送受信を行う観測)、CAPPI観測(複数の仰角に対してPPI観測をそれぞれ行う観測)、及び、RHI観測(方位角を固定した状態で仰角方向にアンテナ12を回転することによりレーダビームの送受信を行う観測)を行うことが可能である。
The
信号処理機18は、PPI観測、CAPPI観測及びRHI観測による受信信号に基づいて、前記反射波の強度を示すエコー強度データを生成してホストPC20に出力する。なお、前記反射波は、雷雲内の水滴に対するレーダビームの反射により発生するので、前記エコー強度データは、前記水滴(水分)の量に応じたデータを示す。
The
ホストPC20は、信号処理機18からのエコー強度データを用いて、図3Aに示すように、同一高度でのエコー強度データ(画像データ)をディスプレイに表示し、一方で、データベース24に記憶されている観測シーケンスに基づいて、アンテナ駆動部14、送受信機16及び信号処理機18を制御する。
The host PC 20 uses the echo intensity data from the
前記観測シーケンスは、PPI観測又はCAPPI観測での仰角やRHI観測での方位角が予め設定されたアンテナ12の走査スケジュール(例えば、CAPPI観測を行った後に所定の方位角でRHI観測を行うスケジュールを1サイクルとし、このサイクルを繰り返し行う走査スケジュール)をいい、ユーザがホストPC20を操作して前記走査スケジュールの設定変更をすることができないようにされている。従って、ホストPC20は、前記観測シーケンスに設定されたPPI観測、CAPPI観測及びRHI観測の走査スケジュールに従ってアンテナ駆動部14、送受信機16及び信号処理機18を制御することによりエコー強度データ(の画像データ)を取得する。
The observation sequence includes a scanning schedule of the
なお、1サイクルのPPI観測又はCAPPI観測及びRHI観測の走査スケジュール(例えば、10分間〜20分間のPPI観測又はCAPPI観測及びRHI観測を1サイクルとする走査スケジュール)を実行することでエコー強度データが得られるので、図3A〜図4Dの順に示すように、1サイクルの観測が終了する毎に、ホストPC20(又は解析用PC22)のディスプレイの画面30の表示内容は順次更新される。すなわち、画面30の表示内容が切り替わるまでは(例えば、画面30が図3Aの表示を行ってから図3Bの表示を行うまでの時間内では)、気象ドップラレーダ23によりPPI観測又はCAPPI観測及びRHI観測が行われていることになる。また、図3A〜図4Dの画面30において、気象ドップラレーダ23の観測可能範囲を示す円32の中心34がアンテナ12の位置を示し、円32内には時間経過に伴って移動する雷雲等の気象現象を示すエコー36〜60が表示されている。
By executing a scanning schedule for one cycle of PPI observation or CAPPI observation and RHI observation (for example, a scanning schedule for 10 to 20 minutes of PPI observation or CAPPI observation and RHI observation as one cycle), echo intensity data is obtained. Thus, as shown in the order of FIGS. 3A to 4D, the display content on the
解析用PC22は、積乱雲等の雷雲のような局地的且つ鉛直方向に成長する気象現象に対して、前記観測シーケンスに基づきアンテナ12を回転しても該気象現象を追尾することが難しい場合に、ホストPC20からのエコー強度データを用いて、RHI観測の方位角を、前記観測シーケンスに予め設定されている方位角から前記気象現象の方位角に変更するためのRHI方位角書き替えコマンドを発行してデータベース24に格納する。ホストPC20は、データベース24に前記RHI方位角書き替えコマンドが格納されている場合にのみ、前記観測シーケンスに設定されているRHI観測の方位角を、前記RHI方位角書き替えコマンドに設定されている方位角に変更し、変更後の前記方位角(観測シーケンス)に基づいてアンテナ駆動部14、送受信機16及び信号処理機18を制御する。
The
次に、解析用PC22におけるRHI方位角書き替えコマンドの発行処理について、図2〜図9Bを参照しながら説明する。
Next, the RHI azimuth rewriting command issuance process in the
図2は、RHI方位角書き替えコマンドの発行処理を示すフローチャートであり、このフローチャートの処理は、上記の1サイクルにつき1回実行される。 FIG. 2 is a flowchart showing the RHI azimuth rewriting command issuance process, and the process of this flowchart is executed once per cycle.
ステップS1において、ホストPC20(図1参照)から解析用PC22に1次データとしてのPPI観測によるエコー強度データが入力される。この場合、気象ドップラレーダ23がCAPPI観測を行っていれば、それぞれの仰角に対するPPI観測を行う毎にホストPC20から解析用PC22にエコー強度データが入力されるので、解析用PC20は、全てのPPI観測のエコー強度データを取得するまで各エコー強度データを蓄積する(ステップS2)。
In step S1, echo intensity data by PPI observation as primary data is input from the host PC 20 (see FIG. 1) to the
次に、ステップS3(第1工程)において、解析用PC22は、蓄積された各エコー強度データに基づいて、同一高度のエコー強度(図3A参照)、エコー頂高度(図5A参照)及び鉛直積算雨水量(図5B参照)を示すプロダクト(画像データ)を生成する。図3A、図5A及び図5Bは、解析用PC22(又はホストPC20)のディスプレイの画面30に表示されるエコー強度、エコー頂高度及び鉛直積算雨水量の画像データを示す模式的説明図である。
Next, in step S3 (first step), the
図3A中のエコー36〜40に対応して、図5Aには、雷雲の頂高度を示すエコー頂高度がエコー62〜78として表示され、一方で、図5Bには、雷雲中の雨水を鉛直方向に積算した鉛直積算雨水量がエコー82〜98として表示されている。なお、図3Aと図5A及び図5Bとでは、互いに対応していないエコーが存在しているが、これは、前述したように、図3Aは、同一高度のエコー36〜40であるため、他の高度で発生している雷雲については画面30に表示されることはなく、一方で、図5A及び図5Bの各エコー62〜78、82〜98は、観測可能範囲(円32)内に発生している全ての雷雲に対応する各エコーを示していることに起因している。
Corresponding to the echoes 36-40 in FIG. 3A, FIG. 5A shows echo top heights indicating the top altitude of thunderclouds as echoes 62-78, while FIG. 5B shows the rainwater in the thunderclouds vertically. The vertically accumulated rainwater amount accumulated in the direction is displayed as echoes 82-98. In FIG. 3A and FIGS. 5A and 5B, there are echoes that do not correspond to each other. As described above, FIG. 3A is the
次に、ステップS4において、解析用PC22は、図5Aのエコー頂高度の画像データ及び図5Bの鉛直積算雨水量の画像データ(プロダクト)を構成するメッシュに対して同時にラスタ走査を行い、対応するメッシュの情報から評価関数法による閾値判定を行う。
Next, in step S4, the
具体的には、図6の模式的説明図に示すように、エコー頂高度又は鉛直積算雨水量の画像データを示すプロダクト100を構成するメッシュ102(エコーを構成するメッシュ102a及びエコーを構成しないメッシュ102b)に対してラスタ走査を行い、エコー頂高度の画像データにおける走査箇所104でのメッシュの値と、鉛直積算雨水量の画像データにおける走査箇所104でのメッシュの値とについて、下記の(1)式に示す重み付けを行い、重み付けした各値を加算する。なお、第1の重み及び第2の重みの値は、雷雲等の気象現象の種類に応じて適宜設定する。
(1つのメッシュ102での鉛直積算雨水量)×(第1の重み)
+(対応するメッシュ102のエコー頂高度)×(第2の重み)>(閾値)
(1)
Specifically, as shown in the schematic explanatory diagram of FIG. 6, the mesh 102 (the
(Vertical accumulated rainwater volume with one mesh 102) x (first weight)
+ (Echo peak height of corresponding mesh 102) × (second weight)> (threshold)
(1)
加算結果と所定の閾値との間で(1)式の関係を満たすならば走査箇所104のメッシュ102の値を「1」とし、一方で、(1)式の関係を満たさなければメッシュ102の値を「0」とする2値化処理を行うことにより、エコー頂高度の画像データ、鉛直積算雨水量の画像データ及び(1)式に基づく新たな画像データ(2値化配列)を生成する。
If the relationship of the expression (1) is satisfied between the addition result and the predetermined threshold value, the value of the
次に、ステップS5において、解析用PC22は、ステップS4において生成した2値化配列の新たな画像データのメッシュに対してラスタ走査を行うことにより、下記のラベリングを行う。
Next, in step S5, the
具体的には、図7A〜図7Dの模式的説明図に示すように、解析用PC22は、2値化配列の画像データ105を構成するメッシュ106(「1」を示すメッシュ106a及び「0」を示すメッシュ106b)に対してラスタ走査を行い、走査箇所108がメッシュ106aに到達した際に(図7A参照)、そのメッシュ106aに対してラベリング番号を付与する。次に、解析用PC22は、ラベリング番号が付与されたメッシュ106a(メッシュ106c)の上下左右のメッシュが、「1」を示すメッシュ106aであるか否かを調べ(図7B参照)、メッシュ106aであれば、そのメッシュ106aに対してラベリング番号を付与する(図7C参照)。このようなラベリング処理を、メッシュ106aの存在を確認することができなくなるまで繰り返し行うことにより、画像データ105において、ラベリング番号が付与された複数のメッシュ106cからなるセル109が構成される(図7D参照)。なお、このセル109の面積が所定のセル化最低面積に到達しない場合(セル109を構成するメッシュ106cの数が所定数に到達しない場合)、解析用PC22は、該セル109を破棄する。
Specifically, as shown in the schematic explanatory diagrams of FIGS. 7A to 7D, the
次に、ステップS6において、解析用PC22は、今回のステップS5の処理により得られたセルと、前回のステップS5の処理により得られたセルとの相関を取り、相関性が高ければ、前回のセルと今回のセルとは同じセルであると共に、前回のセルの位置から今回のセルの位置に移動したものと判断して、解析用PC22に保存されているラベリングの情報を前回のセルの内容から今回のセルの内容に更新する。
Next, in step S6, the
具体的には、図8A及び図8Bの模式的説明図に示すように、解析用PC22は、図8Aの画像データ110中の前回のセル116(前回のステップS5の処理により得られたセル)を上下左右に移動させながら、図8Bの画像データ112中の今回のセル118(今回のステップS5の処理により得られたセル)との相関を取り、重複するメッシュ114の数が互いに最大となる(重複部分の面積が互いに最大となる)セル116、118を同一のセルとみなして、保存されているセル116のラベリング情報(ラベリング番号、鉛直積算雨水量、エコー頂強度、観測時刻等の各種情報)を、セル118のラベリング情報に更新する。すなわち、解析用PC22は、前回の処理(前回のPPI観測)から今回の処理(今回のPPI観測)までの間にセル116がセル118の位置にまで移動したものと判断し、該セル116の移動先に位置するセル118を今回の処理におけるセルとみなして更新処理を行う。なお、画像データ110、112内には、複数のセル120が存在するので、解析用PC22は、各セル120に対しても同様の相関判定を行う。
Specifically, as shown in the schematic explanatory diagrams of FIGS. 8A and 8B, the
ステップS7において、解析用PC22は、気象ドップラレーダ23がエコー強度データ中のエコー(例えば、図3Aのエコー36〜40)を追尾しているか否かを判定し、エコーの追尾をしていなければ、前記エコー強度データ中、注目すべきエコーが存在するか否かを判定し(ステップS8)、注目するエコーが存在すれば、該エコーを注目すべきエコー(注目エコー)として選定する(ステップS9、第2工程)。ステップS8、S9において、注目エコーが存在するか否かの判断処理及び注目エコーを選定するための判断処理としては、例えば、ステップS5、S6のラベリング及び相関判定によって生成されたセル109、116、118に対応するエコーが図3Aのエコー強度データの画像データ中に存在するか否かを判断し、そのようなエコー(図3Aのエコー36)が存在すれば、該エコー36を注目エコーとして選定する。
In step S7, the
次に、ステップS10において、解析用PC22は、注目エコー(エコー36)が、図3Aの観測時刻(例えば、PPI観測を完了した時刻)から図3Bの観測時刻までの間に移動する方向及び距離(移動先)を予測する。
Next, in step S10, the
具体的には、図9Aに示すように、先ず、エコー36に対応するセル130(ステップS5、S6の処理で用いられたセル)について、下記の(2)式及び(3)式に基づいて、前回のステップS6の処理におけるセル130の重心と、今回のステップS6の処理におけるセル130の重心とをそれぞれ算出し、算出した前記各重心より該セル130の移動方向及び移動距離を示す移動ベクトルを求める。
(重心のX座標)=Σ(Xn×メッシュ132aの個数)
/(X方向の座標数) (2)
(重心のY座標)=Σ(Yn×メッシュ132aの個数)
/(Y方向の座標数) (3)
Specifically, as shown in FIG. 9A, first, for the
(X coordinate of the center of gravity) = Σ (Xn × number of
/ (Number of coordinates in the X direction) (2)
(Y coordinate of the center of gravity) = Σ (Yn × number of
/ (Number of coordinates in Y direction) (3)
図9Aは、メッシュ132(セル130を構成するメッシュ132a及びセル130を構成しないメッシュ132b)について、各メッシュ132aの値が略同一である場合(全て「1」である場合)におけるセル130の重心の算出に関する模式説明図であり、X座標(X1〜X5)及びY座標(Y1〜Y5)でのメッシュ132aの個数に基づいて、上記(2)式及び(3)式により重心のX座標及びY座標を算出する。なお、(2)式及び(3)式において、Σは、n=1〜5(Xn=X1〜X5、Yn=Y1〜Y5)での総和を示す数学記号である。
FIG. 9A shows the center of gravity of the
次に、今回のセル130の重心に対応するエコー36の位置を、該エコー36の重心の位置とみなした後に、前記移動ベクトルの始点がエコー36の重心の位置と一致するように、前記移動ベクトルをエコー36に重ね合わせることにより該エコー36の移動先を予測する。すなわち、エコー36の重心の位置を始点とする前記移動ベクトルの終点が該エコー36の移動先として予測される。
Next, after considering the position of the
一方、図9Bは、メッシュ134(セル130を構成するメッシュ134a、134b及びセル130を構成しないメッシュ134c)について、メッシュ134a、134bの値を互いに異なる値に設定した場合(ステップS4〜S6において、2値化処理ではなく、複数の値に数値化した処理を行った場合)のセル130の重心の算出に関する模式説明図である。図9Bの場合では、図3A〜図4Dに示すように、エコー強度データが座標毎に異なるエコー強度を示すことを考慮して、X座標(X1〜X7)及びY座標(Y1〜Y7)について、重みWxn、Wynを用いた下記の(4)式及び(5)式により重心のX座標及びY座標を算出する。
(重心のX座標)=Σ(Xn×Wxn)/(X方向の座標数) (4)
(重心のY座標)=Σ(Yn×Wyn)/(Y方向の座標数) (5)
On the other hand, FIG. 9B shows a case where the values of the
(X coordinate of the center of gravity) = Σ (Xn × Wxn) / (number of coordinates in the X direction) (4)
(Y coordinate of the center of gravity) = Σ (Yn × Wyn) / (number of coordinates in the Y direction) (5)
なお、(4)式及び(5)式において、Σは、n=1〜7(Xn=X1〜X7、Yn=Y1〜Y7、Wxn=Wx1〜Wx7、Wyn=Wy1〜Wy7)での総和を示す数学記号である。また、エコー36の移動先の予測方法については、図9Aの場合と同様であるので、その詳細な説明については省略する。
In equations (4) and (5), Σ is the sum of n = 1 to 7 (Xn = X1 to X7, Yn = Y1 to Y7, Wxn = Wx1 to Wx7, Wyn = Wy1 to Wy7). Mathematical symbol shown. Further, the method of predicting the movement destination of the
次に、ステップS11において、解析用PC22は、エコー36が円32の外部(観測可能範囲外)にあるか否かを判定し、円32内にあれば、中心34と前記移動ベクトルの終点(移動先)とを結ぶ方向(図3Bに示す直線44の方向)をRHI観測の際の方位角とするRHI方位角書き替えコマンドを発行してデータベース24に格納する(ステップS12、第3工程)。
Next, in step S11, the
これにより、次回のサイクルのRHI観測において、気象ドップラレーダ23は、RHI方位角書き替えコマンドの方位角でRHI観測を行い、この結果、図3Bに示す画面30には、エコー36を囲む円46が表示されると共に、エコー36(に応じた気象現象)を前記RHI観測により追尾していることを示す直線44が表示される。
Thereby, in the RHI observation of the next cycle, the
従って、各サイクルで得られたエコー強度データがホストPC20から解析用PC22に入力される毎に、解析用PC22において上記のステップS1〜S12の処理を繰り返し行うことで、図3B〜図4Cに示すように、時間経過に伴い移動する注目エコー(エコー36)を確実且つ自動的に追尾することができる。なお、前述したように、直線44は、中心34と、解析用PC22にて予測した移動先とを結ぶ直線であり、一方で、円46は、前記移動先でのエコー36の重心を示している。そのため、時間経過に伴ってエコー36(に応じた気象現象)が成長し又は衰退しても、該エコー36を確実に捕捉することが可能である。また、図3B〜図4Cには、エコー36以外にも、エコー38〜60が表示されているが、これらのエコー38〜60は、衰退する雷雲、次回の観測時には気象ドップラレーダ23の観測範囲外となる雷雲、あるいは、成長する雷雲であってもエコー36に応じた雷雲程成長しないもの等である。
Therefore, whenever the echo intensity data obtained in each cycle is input from the
なお、図2のフローチャートにおいて、ステップS8で注目エコーが存在しないと判断した場合、解析用PC22は、ステップS9以降の処理、すなわち、RHI方位角書き替えコマンドの発行処理を行わない。
In the flowchart of FIG. 2, if it is determined in step S8 that there is no echo of interest, the
また、ステップS7で注目エコー(エコー36)を追尾中と判定した際に、解析用PC22は、追尾中のエコー36が消失したか否かを判定し、エコー36が消失していなければ、前回のステップS10の処理において予測したエコー36の移動先と、今回のエコー強度データに含まれるエコー36の位置とを比較して、移動予測誤差を算出し(ステップS14)、算出した移動予測誤差を次のステップS10の処理に反映することにより、エコー36の移動先を精度良く算出する。
Further, when it is determined in step S7 that the target echo (echo 36) is being tracked, the
さらに、ステップS13でエコー36が消失したと判定した際に、解析用PC22は、ステップS14以降の処理を行わず、従って、RHI方位角書き替えコマンドの発行を行わない。
Further, when it is determined in step S13 that the
さらにまた、図2及び図4Dに示すように、ステップS10の処理後、エコー36の重心が円32の外部(観測可能範囲外)に位置することにより、円46が円32の外部に位置することになる場合に、解析用PC22は、エコー36の追尾を停止する(ステップS15)。
Furthermore, as shown in FIG. 2 and FIG. 4D, the
以上説明したように、本実施形態に係る自動追尾システム10及び自動追尾方法によれば、気象ドップラレーダ23が観測した気象現象のエコー強度データに基づくエコー頂高度及び鉛直積算雨水量のプロダクト100を用いて、注目したいエコー36を追尾対象として自動的に決定し、さらに、決定したエコー36に応じた気象現象に対してRHI観測を行う際の方位角を自動的に算出するので、注目したいエコー36に応じた積乱雲を含む雷雲等の気象現象を確実且つ自動的に追尾することができる。
As described above, according to the
すなわち、図10Aに示すように、RHI観測の際の方位角が観測シーケンスにより固定とされている場合には、突発的に発生する雷雲等の気象現象を示すエコー140の方向に方位角を向けることができず、該エコー140を追尾することができない。
That is, as shown in FIG. 10A, when the azimuth angle at the time of RHI observation is fixed by the observation sequence, the azimuth angle is directed in the direction of an
これに対して、本実施形態では、図10Bに示すように、解析用PC22にてRHI方位角書き替えコマンドを発行した際に、気象ドップラレーダ23がRHI方位角書き替えコマンド中の方位角にてRHI観測を行うので、注目したいエコー140を確実に狙い打ちできるような追尾を行うことが可能となる。この結果、短時間で急激に発生する積乱雲のような鉛直構造の気象現象を取り逃がすことなく観測することができるので、前記積乱雲に起因した自然災害に対する有意な研究データを容易に収集することができる。また、気象ドップラレーダ23に上記の追尾機能を持たせることにより、前記気象現象の発生から衰退までをRHI観測することが可能となる。
In contrast, in the present embodiment, as shown in FIG. 10B, when the RHI azimuth rewriting command is issued by the
本実施形態に係る自動追尾システム10及び自動追尾方法は、上述した説明に限定されるものではなく、図11に示すように、例えば、3台の気象ドップラレーダ150a〜150cの観測可能範囲152a〜152cの重複部分154内に気象現象156がある場合に、解析用PC22(図1参照)は、各気象ドップラレーダ150a〜150cに対して、方位角の互いに異なるRHI走査方位角書き替えコマンドを発行し、各気象ドップラレーダ150a〜150cは、前記RHI走査方位角書き替えコマンドの方位角を示す直線44a〜44cに基づきRHI観測を行ってもよい。
The
この場合、各気象ドップラレーダ150a〜150cが同時にRHI観測を行っても、あるいは、PPI観測を行う気象ドップラレーダと、RHI観測を行う気象ドップラレーダとに分けて気象現象156に対する観測を行っても、複数の気象ドップラレーダ150a〜150cによる観測が可能となるので、気象現象156を効率よく且つ確実に追尾することができると共に、気象現象156の解析精度を向上することができる。また、PPI観測を行う気象ドップラレーダと、RHI観測を行う気象ドップラレーダとに分けて気象現象156に対する観測を行うことで、RHI観測によるデータと、PPI観測によるデータとを共に得ることが可能となる。
In this case, even if each of the
なお、この発明は、上述の実施形態に限らず、種々の構成を採り得ることは勿論である。 Of course, the present invention is not limited to the above-described embodiment, and various configurations can be adopted.
10…自動追尾システム 20…ホストPC
22…解析用PC 23…気象ドップラレーダ
24…データベース
36、38、40、60、62、78、82、98、140…エコー
44…直線 46…円
109、116、118、120、130…セル
10 ...
22 ... PC for
Claims (2)
求めた前記エコー頂高度及び前記鉛直積算雨水量に基づいて、前記エコー強度データ中の所定のエコーを、追尾すべき注目エコーとして決定する第2工程と、
前記注目エコーに応じた気象現象に対するRHI観測を前記気象レーダが行う際の方位角を算出する第3工程と、
を有することを特徴とする気象現象の自動追尾方法。 A first step of determining an echo top height and a vertical accumulated rainwater amount using echo intensity data indicating a meteorological phenomenon observed by a weather radar;
A second step of determining a predetermined echo in the echo intensity data as a target echo to be tracked based on the obtained echo peak altitude and the vertically integrated rainwater amount;
A third step of calculating an azimuth angle when the weather radar performs RHI observation on a meteorological phenomenon corresponding to the echo of interest;
A method for automatically tracking meteorological phenomena.
前記エコー強度データを用いてエコー頂高度及び鉛直積算雨水量を求め、求めた前記エコー頂高度及び前記鉛直積算雨水量に基づいて、前記エコー強度データ中の所定のエコーを、追尾すべき注目エコーとして決定し、決定した該注目エコーに応じた気象現象に対するRHI観測を前記気象レーダが行う際の方位角を算出し、算出した前記方位角を示す方位書き替えコマンドを発行するコマンド発行手段と、
を有し、
前記気象レーダは、前記コマンド発行手段が前記方位書き替えコマンドを発行した際に、前記注目エコーに応じた気象現象に対するRHI観測の方位角を、前記観測シーケンスに予め設定された方位角から、前記方位書き替えコマンドの方位角に変更して、前記気象現象に対するRHI観測を行う
ことを特徴とする気象現象の自動追尾システム。 A weather radar for acquiring echo intensity data indicating a weather phenomenon by performing PPI observation according to an observation sequence;
Using the echo intensity data, the echo top height and the vertical accumulated rainwater amount are obtained, and based on the obtained echo top altitude and the vertical accumulated rainwater amount, a predetermined echo in the echo intensity data to be tracked to be tracked Command issuing means for calculating an azimuth angle when the weather radar performs RHI observation for a meteorological phenomenon corresponding to the determined echo of interest, and issuing an azimuth rewriting command indicating the calculated azimuth angle;
Have
When the command issuing means issues the direction rewriting command, the weather radar determines an azimuth angle of RHI observation for a meteorological phenomenon corresponding to the echo of interest from an azimuth angle set in advance in the observation sequence. An automatic weather phenomenon tracking system, wherein the azimuth angle of the direction rewriting command is changed to perform RHI observation for the weather phenomenon.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009181539A JP2011033528A (en) | 2009-08-04 | 2009-08-04 | Method and system of automatically tracking weather phenomenon |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009181539A JP2011033528A (en) | 2009-08-04 | 2009-08-04 | Method and system of automatically tracking weather phenomenon |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2011033528A true JP2011033528A (en) | 2011-02-17 |
Family
ID=43762732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009181539A Pending JP2011033528A (en) | 2009-08-04 | 2009-08-04 | Method and system of automatically tracking weather phenomenon |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2011033528A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014048273A (en) * | 2012-09-04 | 2014-03-17 | Toshiba Corp | Weather radar device, phased array radar device and observation sequence preparation method |
JP2014173865A (en) * | 2013-03-06 | 2014-09-22 | Toshiba Corp | Weather radar device, observation sequence creation method and observation sequence creation program |
CN104613928A (en) * | 2015-02-09 | 2015-05-13 | 中国人民解放军63863部队 | Automatic tracking and air measurement method for optical pilot balloon theodolite |
JPWO2015005020A1 (en) * | 2013-07-11 | 2017-03-02 | 古野電気株式会社 | Weather information processing apparatus, weather radar system, and weather information processing method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07110379A (en) * | 1993-10-12 | 1995-04-25 | Tokyo Electric Power Co Inc:The | Thundercloud observation system |
JPH11271443A (en) * | 1998-03-25 | 1999-10-08 | Mitsubishi Electric Corp | Control method for weather radar |
JPH11281740A (en) * | 1998-03-30 | 1999-10-15 | Mitsubishi Electric Corp | Meteorological observation radar system |
-
2009
- 2009-08-04 JP JP2009181539A patent/JP2011033528A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07110379A (en) * | 1993-10-12 | 1995-04-25 | Tokyo Electric Power Co Inc:The | Thundercloud observation system |
JPH11271443A (en) * | 1998-03-25 | 1999-10-08 | Mitsubishi Electric Corp | Control method for weather radar |
JPH11281740A (en) * | 1998-03-30 | 1999-10-15 | Mitsubishi Electric Corp | Meteorological observation radar system |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014048273A (en) * | 2012-09-04 | 2014-03-17 | Toshiba Corp | Weather radar device, phased array radar device and observation sequence preparation method |
JP2014173865A (en) * | 2013-03-06 | 2014-09-22 | Toshiba Corp | Weather radar device, observation sequence creation method and observation sequence creation program |
US9329266B2 (en) | 2013-03-06 | 2016-05-03 | Kabushiki Kaisha Toshiba | Weather radar apparatus, observation sequence generation method, and observation sequence generation program |
JPWO2015005020A1 (en) * | 2013-07-11 | 2017-03-02 | 古野電気株式会社 | Weather information processing apparatus, weather radar system, and weather information processing method |
US10139487B2 (en) | 2013-07-11 | 2018-11-27 | Furuno Electric Co., Ltd. | Meteorological radar system and method and device for processing meteorological information |
CN104613928A (en) * | 2015-02-09 | 2015-05-13 | 中国人民解放军63863部队 | Automatic tracking and air measurement method for optical pilot balloon theodolite |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Heinselman et al. | High-temporal-resolution capabilities of the national weather radar testbed phased-array radar | |
US10915673B2 (en) | Device, method, apparatus, and computer-readable medium for solar site assessment | |
AU2011311805B2 (en) | System and method for generating derived products in a radar network | |
US8068050B2 (en) | Methods and systems for detection of hazard to aviation due to convective weather | |
CN109425855A (en) | It is recorded using simulated sensor data Augmented Reality sensor | |
US7486220B1 (en) | Storm top detection and prediction | |
CN109427214A (en) | It is recorded using simulated sensor data Augmented Reality sensor | |
JP7182869B2 (en) | Target detection device | |
JP5654223B2 (en) | Image processing apparatus, radar apparatus equipped with the same, image processing method, and image processing program | |
KR101590253B1 (en) | Method and device for simulation of sonar images of multi-beam imaging sonar | |
JP5394690B2 (en) | Wind prediction apparatus and program | |
CN110068836A (en) | A kind of laser radar curb sensory perceptual system of intelligent driving electric cleaning car | |
JP2010181406A (en) | System and method for displaying radar-measured turbulence intensity on vertical display | |
CN110297246B (en) | Networking X-waveband weather radar collaborative self-adaptive control method and system | |
JP6689396B2 (en) | Weather forecasting device, weather forecasting method, and program | |
CN108828691A (en) | It is short to face precipitation forecast method and device | |
JP2011033528A (en) | Method and system of automatically tracking weather phenomenon | |
CN115461262A (en) | Autonomous driving using surface element maps | |
KR101036823B1 (en) | Radar observational density analysis system and analysis method of the same | |
CN117269940B (en) | Point cloud data generation method and perception capability verification method of laser radar | |
JP5362493B2 (en) | Radar signal processing device | |
CN109444893A (en) | Phase identification product pattern splicing method and device based on dual polarization radar net | |
JP6316970B2 (en) | Weather radar control apparatus, weather observation system, weather radar control method, and weather radar control program | |
KR20130049521A (en) | Apparatus and method for predicting precipitation echo route | |
JP6139193B2 (en) | Radar display device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120803 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130725 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130730 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130919 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20131126 |