JPH10260267A

JPH10260267A - 局所雨域予測方法および装置

Info

Publication number: JPH10260267A
Application number: JP6466597A
Authority: JP
Inventors: Hidetomo Sakaino; 英朋境野; Tsutomu Horikoshi; 力堀越; Satoshi Suzuki; 智鈴木
Original assignee: Nippon Telegraph and Telephone Corp
Current assignee: Nippon Telegraph and Telephone Corp
Priority date: 1997-03-18
Filing date: 1997-03-18
Publication date: 1998-09-29
Anticipated expiration: 2017-03-18
Also published as: JP3377076B2

Abstract

(57)【要約】【課題】流体方程式を介して全体の移動速度を推定
し、同時に予測を行う、演算コストが大幅に軽減され
た、局所雨域予測方法および装置を提供する。【解決手段】画像入力部１１を用いて雨域のレーダー
画像を入力し、輪郭線・法線検出部１３で、画像蓄積部
１２から取り出した画像に前処理を行い、いくつかのブ
ロックを検出し、輪郭線を算出し、輪郭線の接線と法線
とを算出し、変位ベクトル検出部１４で、法線を基点か
ら輪郭線の外側ならびに内側方向に、後時点の雨域の輪
郭線と交わるまで延長し、法線ベクトルと法線の速度ベ
クトルとを算出し、流体方程式補間部１５で、流体方程
式に速度ベクトルを初期値として代入して、速度ベクト
ル値を補間し、繰り返して補間を行って予測速度ベクト
ルを算出し、雨域の形状濃淡値情報とともに濃度拡散・
移流方程式に代入することで、雨域の変化予測を行い、
出力部１６で予測画像を外部に出力する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、気象予測方法およ
び装置に関し、特に局所における雨域の短時間予測方法
および装置に関する。

【０００２】

【従来の技術】従来、気象予測を行う際、気象庁ではア
メダスや気象衛星から得られる２、３次元的な気温、気
圧、露点、風ベクトル等の様々な物理パラメーターを用
い、数理・物理方程式を介して、数十時間先もしくは数
日先の日本全国の気象現象の変化を予測している。一
方、近年になってより狭い地域での、分オーダーでの短
い時間における気象予測、特に降雨や降雪の情報が望ま
れてきているため、局地的に配置した気象レーダー装置
からの雨域の変化情報を用いた予測が重要となってきて
いる。

【０００３】雨域の変化は、降雨や降雪状態にある領域
の時間的な変化を意味し、その領域は見かけ上、雲の形
状をなしている。画像としての雨域は、様々な形状と濃
淡値の組み合わせとして表現されており、明瞭な特徴量
の記述が困難であることから画像処理手法に基づいた雨
域の動きの検出は容易ではない。そのため現在でも、主
として人の経験則に基づいた雨域の変化からの主観的な
読み取りが行われ、その情報が活用されている。

【０００４】しかし一方では、自動的かつ客観的な雨域
変化を計算機上で行うために、時空間的相関法による雨
域変化検出が試みられている。この方法は、２つの画像
フレーム間〔ｆ（ｔ）とｆ（ｔ＋１）：ｔは時間ステッ
プ〕の、濃淡値の類似性を算出した後、最も相関値が高
い点が移動量に対応していると仮定した上で、雨域の移
動量を推定する。推定された移動量に基づいて雨域を平
行移動させていく。このようにして推定された次のステ
ップの雨域画像を繰り返し、特定数のステップ先までの
予測が可能となる。

【０００５】従来より、相互相関法、重心法などの画像
処理を行うことによって、時系列的な雨域パターンから
移動ベクトルを求めて、その移動ベクトルを線形に外挿
する方法が中心的な予測方法であった。しかしながら、
これらの外挿法では大きな問題がある。それは、必ずし
も各移動ベクトルが精度よく推定されたわけではないの
と、物理的に複数の移動ベクトル間を補正する手段が考
慮されていないため、外挿を数ステップ反復するだけ
で、各ベクトルが不自然に交錯してしまうからである。
このとき、ある大きさの雨域パターンを平行移動させて
いるので、交錯時には、雨量の精度は重なった分だけ多
く予測されてしまう。また、方向も大きくずれて行って
しまうといった問題がある。

【０００６】

【発明が解決しようとする課題】これまで適用されてい
る時空間的相関法では、連続する２つの画像から１つの
動きベクトルを推定する。しかしながら、相関係数を算
出するためには画素ごとに演算を行う必要があるため
に、画像全体で数千回〜数万回オーダーの演算回数が要
求されるという問題がある。また、画像の濃淡値は、雨
域が生成・消滅を絶えまなく繰り返しているために、前
フレームの濃淡値と次フレームの濃淡値とが常に１対１
に対応していることはないので、濃淡値に基づいた相関
値への信頼性は必ずしも高くなくて、むしろ不安定であ
るという問題がある。さらに、相関係数が求められない
ような散漫な雨域画像の場合には、ほとんど動きベクト
ルを推定することができず、その間の補間方法について
も具体的な方法が見い出されていないのが現状である。

【０００７】また、相関法では、単位時間当たりの移動
量が大きいときは、大き目の窓関数を、移動量が小さい
ときは、小さ目の窓関数を適宜用意する必要があり、最
適な窓関数の大きさを決定することは容易ではなく、同
時に、決定しようとするならば、それに要する演算時間
は全体として数〜十倍増加するという問題が生ずる。

【０００８】さらに、画像を分割した数だけ、通常数十
個程度の、雨域の移動速度を推定することができない。

【０００９】特に、雨域が停滞性の場合、雨域の輸郭形
状の変化に比べると、雨域の濃淡値の変化が激しい。こ
のような雨域において、相関法をそのまま適用すると、
全体の雨域の移動はなくても、雨域の濃淡変化があるた
めに、何らかの動きベクトルが複数個得られてしまい、
その動きベクトルで雨域を移動させていた。そのため、
予測が進むほど、雨域の領域が実際からはずれていって
しまうという問題があった。

【００１０】雨域が停滞性であっても、フレーム間の雨
域パターンの精緻な対応づけは困難である。従って、人
の経験に基づいて雨域パターンの移動速度を推定するこ
とが実用上中心的となっている。

【００１１】また降水域の変化が一定という強い仮定が
おかれた単純なモデルであったために、対流に伴った降
水域の変化の程度を勘案した予測方法が試みられていな
かった。

【００１２】本発明は上述の問題点の解決を図り、相関
法を用いずに、レーダー画像から局所的な雨域の移動速
度を推定し、流体方程式を介して全体の移動速度を推定
し、同時に予測を行う、演算コストが大幅に軽減され
た、局所雨域予測方法および装置を提供することを目的
とする。

【００１３】

【課題を解決するための手段】本発明の局所雨域予測方
法および装置は、計算手段を有する局所雨域予測手段を
用いて局所雨域を予測する方法であって、局所雨域予測
手段の画像入力部を用いて予測の対象とする雨域のレー
ダー画像を入力し、画像蓄積部に入力された時系列画像
を蓄積し、輪郭線・法線検出部で、計算のために画像蓄
積部から取り出した時系列の雨域のレーダー画像に前処
理を行った後、それぞれから雲のかたまりであるいくつ
かの適当なブロックを検出し、検出されたブロックの輪
郭線を算出し、算出された輪郭線を所定の手順によって
複数の輪郭点で離散化し、輪郭点を基点として輪郭線の
接線と法線とを算出し、変位ベクトル検出部で、前時点
の法線を基点から輪郭線の外側ならびに内側方向に、後
時点の雨域の輪郭線と交わるまで延長し、該基点と該交
点の座標から法線ベクトルと法線の速度ベクトルとを算
出し、流体方程式補間部で、流体方程式の該当項に、部
分的に得られた速度ベクトルを初期値として代入して、
該画像のフレーム全体にわたって速度ベクトル値を補間
し、必要とする未来の時間まで繰り返して補間を行って
所望の予測速度ベクトルを算出し、算出された該速度ベ
クトルと雨域の形状濃淡値情報とを濃度拡散・移流方程
式の所定の項に代入することで、雨域の変化予測を行
い、出力部で、得られた変化予測に基づく予測画像を外
部に出力する。

【００１４】輪郭線・法線検出部の前処理には、雨域画
像の２値化、膨張収縮、孤立点除去、ラベリングが含ま
れることが好ましく、変位ベクトル検出部の輪郭線の接
線と法線との算出には、輪郭点を中心とする中心差分法
が用いられることが好ましい。

【００１５】また、流体方程式補間部で用いられる流体
方程式がナビエ・ストークス方程式であってもよく、バ
ーガーズ方程式であってもよい。

【００１６】さらに、流体方程式補間部で用いられる濃
度拡散・移流方程式には、湧き出し項、吸い込み項、消
散項が含まれることが好ましい。

【００１７】本発明の局所雨域予測装置は、局所雨域を
予測するための計算手段を有する局所雨域予測装置であ
って、予測の対象とする雨域のレーダー画像を入力する
ための画像入力部と、入力された時系列画像を蓄積する
画像蓄積部と、計算のために画像蓄積部から取り出され
た時系列の雨域のレーダー画像に前処理を行った後、そ
れぞれから雲のかたまりであるいくつかの適当なブロッ
クを検出し、検出されたブロックの輪郭線を算出し、算
出された輪郭線を所定の手順によって複数の輪郭点で離
散化し、輪郭点を基点として輪郭線の接線と法線とを算
出する輪郭線・法線検出部と、前時点の法線を基点から
輪郭線の外側ならびに内側方向に、後時点の雨域の輪郭
線と交わるまで延長し、該基点と該交点の座標から法線
ベクトルと法線の速度ベクトルとを算出する変位ベクト
ル検出部と、流体方程式の該当項に部分的に得られた速
度ベクトルを初期値として代入して、該画像のフレーム
全体にわたって速度ベクトル値を補間し、必要とする未
来の時間まで繰り返して補間を行って所望の予測速度ベ
クトルを算出し、算出された該速度ベクトルと雨域の形
状濃淡値情報とを濃度拡散・移流方程式の所定の項に代
入することで、雨域の変化予測を行う流体方程式補間部
と、得られた変化予測に基づく予測画像を外部に出力す
る出力部とを備える。

【００１８】本発明では、時間的に異なるフレームにお
ける雨域を入力し、入力データより過去の雨域の輪郭線
を求めた後、その輸郭線にまばらに設定された法線を輪
郭線の内外に対して延長し、現在のフレームにおける雨
域との交点を求め、輸郭線の変位から速度ベクトルを推
定することにより局所的な雨域の移動速度が推定され
る。続いて、速度ベクトルの全体的な整合をとるため
に、求まっている局所的な速度ベクトルを初期値として
流体方程式に代入する。同時に、流体方程式は非定常な
表現があるために、時間について方程式を積分し続けて
反復的に補間することによりより密に分布した速度ベク
トルに関する予測ができる。さらに、算出された速度ベ
クトルを濃度拡散・移流方程式の速度（移流）項に代入
し雨域の形状濃淡値情報を濃度拡散・移流方程式の拡散
項、移流項、湧きだし項、吸い込み項、消散項に代入
し、時間積分を行うことで、次の時間の濃淡値分布の変
化を求める。この時間積分を繰り返すことで濃淡値変化
の予測結果が得られ、画素単位の移流速度に基づいて、
雨域を変化させていくことができる。

【００１９】本発明は、気象レーダーエコー画像に基づ
いた狭い地域の雨域の時系列変化を短時間予測し、気象
の変化に敏感に影響されるような狭い地域における様々
な経済活動が円滑に営まれるように、高速にかつ簡易に
局所的な雨域の移動速度を推定し、同時にその変化の予
測を雨域の変化に反映させて、画像の形で２次元的に予
測する方法と装置である。

【００２０】

【発明の実施の形態】次に、本発明の実施の形態につい
て図面を参照して説明する。図１は本発明の実施の形態
の局所雨域予測方法と装置の説明図であり、（ａ）は局
所雨域予測方法のフローチャート、（ｂ）は局所雨域予
測装置のブロック構成図である。図中符号１１は画像入
力部、１２は画像蓄積部、１３は輪郭線・法線検出部、
１４は変位ベクトル検出部、１５は流体方程式補間部、
１６は出力部、Ｓ１０１〜Ｓ１１３は各ステップであ
る。

【００２１】本発明の局所雨域予測装置は、気象レーダ
エコー画像などから所望の地域の所望の雨域画像を入力
する画像入力部１１と、画像入力部１１から入力した画
像データを蓄積する画像蓄積部１２と、画像蓄積部１２
から予測のために比較する現在と過去の画像データを取
り出し、比較のための前処理を行い、画像データの雨域
を雲の塊ごとにブロックとして分割し、ブロックごとに
輪郭線とその法線を算出する輪郭線・法線検出部１３
と、各ブロックの前時点の法線を延長し、後時点の輪郭
線との交点を算出し、速度ベクトルを算出する変位ベク
トル検出部１４と、所定の流体方程式の所定の項に算出
された局所的な速度ベクトルを初期値として代入し、必
要時間分だけ時間積分を繰り返して予測速度値を計算
し、濃度拡散・移流方程式の所定の項に、算出された予
測速度値と形状濃度情報を代入して雨域変化を予測する
流体方程式補間部１５と、流体方程式補間部１５で予測
された雨域変化を画像として出力する出力部１６とから
構成される。

【００２２】次に図１（ａ）を参照して本発明の局所雨
域予測方法の各過程を説明する。予測を開始すると（Ｓ
１０１）、画像入力部１１を用いて外部の気象レーダエ
コー画像などから所望の地域の所望の雨域のレーダー画
像を入力し（Ｓ１０２）、画像蓄積部１２に蓄積する
（Ｓ１０３）。画像蓄積部１２には過去の時系列画像が
蓄積される。

【００２３】次に画像蓄積部１２から現在の画像と、比
較の対象とする過去の画像を取り出し（Ｓ１０４）、輪
郭線・法線検出部ｌ３で、それぞれの画像に対して雨域
画像の２値化、膨張収縮、孤立点除去、ラベリングとい
った一連の前処理を行い（Ｓ１０５）、雨域をいくつか
の雲の塊のブロックに分割し（Ｓ１０６）、同一のラベ
リングをもつブロックごとに輪郭線とその法線を求める
（Ｓ１０７）。必要によりこれを現在を含め比較の対象
とするいくつか過去のフレームについて行う。各ブロッ
クの雨域の画像のフレーム間での移動速度を推定するた
めに、変位ベクトル検出部１４で、まず、各ブロックの
前時点の法線を内側および外側に延長し、後時点の雨域
の対応するブロックの輪郭線との交点を算出し、法線長
さと経過時間から速度ベクトルを算出する（Ｓ１０
８）。次に流体方程式補間部１５で、局所（部分）的に
推定された速度ベクトルから、その近傍の速度ベクトル
を推定するために、流体方程式の所定の項に、算定され
た局所的な速度ベクトルを初期値として代入して（Ｓ１
０９）、必要な時間だけ時間積分を繰り返して予測速度
値を計算する（Ｓ１１０）。さらに、湧きだし項、吸い
込み項、消散項をも有する濃度拡散・移流方程式の所定
の項に予測速度値と形状濃淡情報を代入して、雨域の変
化を予測する（Ｓ１１１）。予測の結果を出力部１６で
外部へ画像として出力して（Ｓ１１２）、工程を終了す
る（Ｓ１１３）。

【００２４】図２は、本発明の雨域の個々のブロックか
ら法線・接線を求める方法（Ｓ１０７）を説明するため
の模式的平面図であり、図中符号２１は輪郭線、２２は
輪郭点、２３は接線、２４は法線である。ここでは、雨
域の輪郭線２１は仮に円パターンとなっているものとす
る。円パターンを適当な数からなる輪郭点２２によって
離散化する。図の輪郭点２２上の接線２３の方向は、こ
の輪郭点２２を中心としてその両隣の輪郭点２２を用い
た中心差分値の成分である。従って、法線２４の方向
は、この成分との内積がゼロとなるように求めた結果に
等しい。このような計算を離散化したすべての輪郭点２
２の点列について行う。なお、実際の複雑な輪郭形状を
もつ雨域パターンに適用するためには、前処理として、
２値化、膨張収縮、孤立点除去、ラベリングを施してお
く必要がある。

【００２５】図３は、本発明で雨域のブロックの２つの
画像のフレーム間で輪郭線がひろがった例を示す模式的
平面図であり、（ａ）は輪郭線が等方的に拡がった状態
を、（ｂ）は輪郭線が異方的に拡がった状態を示す。図
中符号３１、３６は過去の輪郭線、３４、３７は法線、
３５、３８は現在の輪郭線である。

【００２６】（ａ）と（ｂ）では小円の大円に対する配
置のずれ方が異なっている。ここでは、小円を過去のあ
る雨域、大円を現在の大きさの雨域であることを想定し
ている。（ａ）の場合は、大小２つの円の変形量は、明
らかに等方的に得られる。一方、（ｂ）の場合は、上方
に大きく、下方に小さい変形量として得られる。このよ
うに輪郭線の変化は通常一様ではない。

【００２７】図４は、本発明で実際の雨域について輪郭
線の変形量を算定した結果を示す模式的平面図であり、
図中符号４１は過去の輪郭線、４４は法線、４５は現在
の輪郭線である。小さい（網点の）雨域のブロックが、
画像の数フレーム先で大きく広がった場合の例である。
網点の輪郭線４１の法線４４を内外に伸ばして、大きい
方の輪郭線４５との交わったところまでの変形量によっ
て、相対的な移動速度が求まる。この場合、全体的に雨
域は上方に広がっているが、下方への移動もあるため
に、雨域の下部は、凸形状と凹形状とに変形しているこ
とがわかる。本発明では法線を外部と内部に延長して後
時点の輪郭線との交点を求めるので、このような凹形状
の変化にも対応できる。

【００２８】図５は、本発明で局所的な変形量（移動速
度）から、全体の移動速度を推定した結果を示す摸式的
平面図であり、図中符号５１は過去の輪郭線、５４は速
度ベクトル、５５は現在の輪郭線である。部分的に求ま
っている初期の移動速度をナビエ・ストークス方程式に
代入することで、全体の移動速度、即ち移流速度が予測
される。時間積分を続けると、次の画像のフレームに、
予測された移流速度を容易に得ることができる。

【００２９】図６は、実際に局所雨域予測装置の流体方
程式補間部で過去と現在の輪郭点を算出し変位方向をデ
ィスプレー上に表示した中間調画像であり、図中符号６
２は過去の輪郭点、６４は速度ベクトル、６７は現在の
輪郭点である。

【００３０】本発明では、速度ベクトルを輪郭線の変位
から粗く求めた後、移動平均を反復的に繰り返すため
に、各速度ベクトルは滑らかなつながりをもち、雨域パ
ターンは流体方程式に従って、物理的に変形していく。
また、移動平均による反復的な補間、補正に加えて、輪
郭線の変位から求まった速度ベクトルを初期ベクトルと
して、ナビエ・ストークス方程式に流速ベクトルとして
代入して、次のステップの速度ベクトルを物理的に補
間、補正することもできる。ナビエ・ストークス方程式
を用いていると、特に、渦パターンの生成・消滅予測が
可能である。

【００３１】実施の形態の説明では流体方程式補間部で
用いられる流体方程式としてナビエ・ストーク方程式を
用いる例で説明したが、バーガーズ方程式を用いても類
似の結果が得られる。

【００３２】即ち、従来予測方法では、せいぜい３０分
程度先の予測には有効な場合もあるが、本発明のよう
に、数時間以上先まで安定にかつ精度よく予測すること
は不可能である。さらに、どの従来方法も移動ベクトル
の推定の仕方に差異があるものの、雨域の形状変化と濃
淡変化を同時に物理的に予測する方法は存在しなかっ
た。従って対流に伴った、雨域パターンの生成・消滅が
全く考慮されていなかった。本発明ではこれらのことが
考慮されており、細かい雨域パターンが速度ベクトルに
乗って近づいて融合する状態を含めて流体方程式系に基
づいて予測することが最も自然かつ容易な方法である。

【００３３】本発明では、曲線の変位から速度ベクトル
を推定している。しかしながら、従来の発明では、速度
ベクトルを求める際、本発明の方法のように、法線を輪
郭線の外側と内側に探索することが考慮されていなかっ
た。

【００３４】図７は、本発明と従来法による予測結果の
比較グラフである。雨域の面積の的中率として広く用い
られている、ＣＳＩ（ＣｒｉｔｉｃａｌＳｕｃｃｅｓ
ｓＩｎｄｅｘ）を用いると、従来法（相互相関法＋外挿
法やニューラルネットモデル）よりも、流体方程式に基
づいた方法の方が経過時間による的中率の低下が抑制さ
れているのがわかる。特に、本発明を適用して、局所的
な移流速度を流体方程式に与えた場合は、移流速度を雨
域場に大局的移流動速度として一様に与えた場合よりも
的中率が向上していることがわかる。

【００３５】以上の性能評価結果から、本発明の予測の
精度の高さが容易にわかる。即ち、従来方法よりも、高
速に、しかも、３時間以降でも実用的な的中率が本発明
により達成できた。

【００３６】

【発明の効果】以上説明したように、本発明によれば、
局所的な移流速度を相関法を用いずに、輪郭線とそれに
対する法線を利用することで、推定でき、さらに流体方
程式に代入することで、全体の速度ベクトルが推定で
き、その結果高速かつ簡易に局所的な雨域の移動速度が
推定でき、かつ雨域の予測精度を大幅に向上させること
ができるという効果がある。

【図面の簡単な説明】

【図１】本発明の実施の形態の局所雨域予測方法と装置
の説明図である。（ａ）は局所雨域予測方法のフローチ
ャートである。（ｂ）は局所雨域予測装置のブロック構
成図である。

【図２】本発明の雨域の個々のブロックから法線・接線
を求める方法を説明するための模式的平面図である。

【図３】本発明で雨域のブロックの２つの画像のフレー
ム間で輪郭線がひろがった例を示す模式的平面図であ
る。（ａ）は輪郭線が等方的に拡がった状態を示す。
（ｂ）は輪郭線が異方的に拡がった状態を示す。

【図４】本発明で実際の雨域について輪郭線の変形量を
算定した結果を示す模式的平面図である。

【図５】本発明で局所的な変形量（移動速度）から、全
体の移動速度を推定した結果を示す摸式的平面図であ
る。

【図６】実際に局所雨域予測装置の流体方程式補間部で
過去と現在の輪郭点を算出し変位方向をディスプレー上
に表示した中間調画像である。

【図７】本発明と従来法による予測結果の比較グラフで
ある。

【符号の説明】

１１画像入力部１２画像蓄積部１３輪郭線・法線検出部１４変位ベクトル検出部１５流体方程式補間部１６出力部２１輪郭線２２輪郭点２３接線２４、３４、３７、４４法線３１、３６、４１、５１過去の輪郭線３５、３８、４５、５５現在の輪郭線５４速度ベクトル６２過去の輪郭点６４速度ベクトル６７現在の輪郭点Ｓ１０１〜Ｓ１１３各ステップ

Claims

【特許請求の範囲】

【請求項１】計算手段を有する局所雨域予測手段を用
いて局所雨域を予測する方法であって、前記局所雨域予測手段の画像入力部を用いて予測の対象
とする雨域のレーダー画像を入力し、画像蓄積部に入力された時系列画像を蓄積し、輪郭線・法線検出部で、計算のために前記画像蓄積部か
ら取り出した時系列の雨域のレーダー画像に前処理を行
った後、それぞれから雲のかたまりであるいくつかの適
当なブロックを検出し、検出された前記ブロックの輪郭
線を算出し、算出された前記輪郭線を所定の手順によっ
て複数の輪郭点で離散化し、前記輪郭点を基点として前
記輪郭線の接線と法線とを算出し、変位ベクトル検出部で、前時点の前記法線を前記基点か
ら前記輪郭線の外側ならびに内側方向に、後時点の雨域
の輪郭線と交わるまで延長し、該基点と該交点の座標か
ら法線ベクトルと法線の速度ベクトルとを算出し、流体方程式補間部で、流体方程式の該当項に、部分的に
得られた前記速度ベクトルを初期値として代入して、該
画像のフレーム全体にわたって速度ベクトル値を補間
し、必要とする未来の時間まで繰り返して補間を行って
所望の予測速度ベクトルを算出し、算出された該速度ベ
クトルと前記雨域の形状濃淡値情報とを濃度拡散・移流
方程式の所定の項に代入することで、雨域の変化予測を
行い、出力部で、得られた変化予測に基づく予測画像を外部に
出力する、ことを特徴とする局所雨域予測方法。
【請求項２】前記輪郭線・法線検出部の前記前処理に
は、雨域画像の２値化、膨張収縮、孤立点除去、ラベリ
ングを含む、請求項１に記載の局所雨域予測方法。
【請求項３】前記変位ベクトル検出部の輪郭線の接線
と法線との算出には、前記輪郭点を中心とする中心差分
法が用いられる、請求項１に記載の局所雨域予測方法。
【請求項４】前記流体方程式補間部で用いられる前記
流体方程式がナビエ・ストークス方程式である、請求項
１に記載の局所雨域予測方法。
【請求項５】前記流体方程式補間部で用いられる前記
流体方程式がバーガーズ方程式である、請求項１に記載
の局所雨域予測方法。
【請求項６】前記流体方程式補間部で用いられる濃度
拡散・移流方程式には、湧き出し項、吸い込み項、消散
項を含む、請求項１に記載の局所雨域予測方法。
【請求項７】局所雨域を予測するための計算手段を有
する局所雨域予測装置であって、予測の対象とする雨域のレーダー画像を入力するための
画像入力部と、入力された時系列画像を蓄積する画像蓄積部と、計算のために前記画像蓄積部から取り出された時系列の
雨域のレーダー画像に前処理を行った後、それぞれから
雲のかたまりであるいくつかの適当なブロックを検出
し、検出された前記ブロックの輪郭線を算出し、算出さ
れた前記輪郭線を所定の手順によって複数の輪郭点で離
散化し、前記輪郭点を基点として前記輪郭線の接線と法
線とを算出する輪郭線・法線検出部と、前時点の前記法線を前記基点から前記輪郭線の外側なら
びに内側方向に、後時点の雨域の輪郭線と交わるまで延
長し、該基点と該交点の座標から法線ベクトルと法線の
速度ベクトルとを算出する変位ベクトル検出部と、流体方程式の該当項に部分的に得られた前記速度ベクト
ルを初期値として代入して、該画像のフレーム全体にわ
たって速度ベクトル値を補間し、必要とする未来の時間
まで繰り返して補間を行って所望の予測速度ベクトルを
算出し、算出された該速度ベクトルと前記雨域の形状濃
淡値情報とを濃度拡散・移流方程式の所定の項に代入す
ることで、雨域の変化予測を行う流体方程式補間部と、得られた変化予測に基づく予測画像を外部に出力する出
力部と、を備えたことを特徴とする局所雨域予測装置。