JP3377078B2 - 対流場変化予測推定装置 - Google Patents

対流場変化予測推定装置

Info

Publication number
JP3377078B2
JP3377078B2 JP17718397A JP17718397A JP3377078B2 JP 3377078 B2 JP3377078 B2 JP 3377078B2 JP 17718397 A JP17718397 A JP 17718397A JP 17718397 A JP17718397 A JP 17718397A JP 3377078 B2 JP3377078 B2 JP 3377078B2
Authority
JP
Japan
Prior art keywords
image
moving
moving speed
convection
vorticity
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
Application number
JP17718397A
Other languages
English (en)
Other versions
JPH1123726A (ja
Inventor
英朋 境野
智 鈴木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP17718397A priority Critical patent/JP3377078B2/ja
Publication of JPH1123726A publication Critical patent/JPH1123726A/ja
Application granted granted Critical
Publication of JP3377078B2 publication Critical patent/JP3377078B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、時系列画像中の非
剛体の物体の移動速度を検出し、同時に、その予測を行
うことが必要とされる気象レーダーエコー画像からの降
水量の変化予測、流体工学における流体の挙動の解析
等、非剛体の動きを検出してその挙動を予測する分野に
属する。
【0002】
【従来の技術】従来、画像中の剛体系の物体の移動ベク
トルを検出する場合は、剛体系の物体上の照明変化がほ
とんどないモデルを採用していることが多い。その検出
方法には、オプティカルフローや相互相関法に基づいた
方法が中心的である(文献[1])。([1]Dona H.B
allard,“コンピュータビジョン”、日本コンピュータ
ビジョン協会)。
【0003】一方、非剛体の物体に対する適切な移動ベ
クトルはないと言える。これは、連続する画像間であっ
ても、非剛体の物体の輪郭線、濃淡値等の属性が同時に
変化するために、明確な対応づけを行えないことに起因
する。即ち、複雑な非剛体の物体に対しては、統計的な
類似性を追従していく相互相関法が適用されることが多
い。例えば、気象レーダーエコー画像を用いた場合、降
水パターンは非剛体的に変化するものの、刻々変化する
画像間に相互相関法が適用されている。
【0004】
【発明が解決しようとする課題】しかし、最近では、降
水パターンの変化予測が必要とされつつあるが、移動ベ
クトルを線形に外挿する方法をとっている。これによ
り、降水パターンの移動が表現できるが、非線形な軌跡
を示すことが多い場合には、この外挿方法によると未来
の移動ベクトルほど、精度は大きく下がってしまう問題
がある。
【0005】また、検出された流れ場から、対流変化を
予測する方法が取り入れられていないために、降水量の
変化、特に、湧きだしと吸い込み領域を予測することが
できなく、短時間気象予測で必要とされる擾乱の変化に
対応することができない。
【0006】本発明の目的は、非剛体で非線形な移動ベ
クトル変化から湧きだしと吸い込み領域を予測する対流
場変化予測推定装置を提供することである。
【0007】
【課題を解決するための手段】本発明の対流場変化予測
推定装置は、湧き出しと吸い込み領域を場として有する
対流場変化予測推定装置であって、画像入力する画像入
力手段と、入力された画像を時系列画像として蓄積する
画像蓄積手段と、画像蓄積手段に蓄積されている連続す
る2つ以上の2次元画像間で画像処理を行って、画像特
徴量を抽出しながら画像中の物体の移動速度を検出する
移動速度検出手段と、画像特徴量から検出され複数の
移動速度を初期移動速度として、移動平均フィルターを
反復的に適用して初期移動速度をゼロ速度領域へ伝播さ
せる反復的移動速度伝播手段と、移動速度を初期値とし
て、渦度と流れ関数を連立させて、渦度と流れ関数の変
化を算出する対流場算出手段と、渦度と流れ関数の密度
分布から流れ場の変動を予測する対流度予測手段と、対
流度予測手段の予測結果を出力する出力手段を有する。
【0008】また、画像入力手段は、気象レーダーエコ
ー画像とドップラー画像から観測された情報を入力する
画像入力手段を有する。
【0009】また、移動速度検出手段は、画像を複数の
小領域に分割して、連続する画像の各小領域間での相互
相関係数を解析して、最も相関係数の高い点を求めて小
領域の移動ベクトルとして検出する移動速度検出手段を
有する。
【0010】更に、移動速度検出手段は、渦度と流れ関
数の密度分布に基づいて、画像の濃淡値を変数にもつ流
体方程式の湧きだし項と吸い込み項に特定な比例定数を
含み、湧きだしと吸い込み領域の移動速度を検出する移
動速度検出手段を有する。
【0011】本発明では、限られた画像情報から局所的
に移動ベクトルを検出し、ゼロ領域へ反復法により局所
移動ベクトルを伝播させた後、流れ場の移動ベクトルを
渦度と流れ関数に置き換えて、渦度と流れ関数を連立し
て解き、その時間的な推移を求めながら、渦度と流れ関
数の密度分布の関係から、対流に伴った擾乱する領域の
予測を反復していく。
【0012】
【発明の実施の形態】次に、本発明の実施の形態につい
て、図面を参照して、詳細に説明する。
【0013】図1は、本発明の一実施形態の対流場変化
予測推定装置の構成図である。本装置は、画像入力する
画像入力部100と、入力された画像を時系列画像とし
て蓄積する画像蓄積部110と、画像蓄積部に蓄積され
ている連続する2つ以上の2次元画像間で画像処理を行
って、画像特徴量を抽出しながら画像中の物体の移動速
度を検出する移動速度検出部120と、画像特徴量から
検出され複数の移動速度を初期移動速度として、移動
平均フィルターを反復的に適用して初期移動速度をゼロ
速度領域へ伝播させる反復的移動速度伝播部130と、
移動速度を初期値として、渦度と流れ関数を連立させ
て、渦度と流れ関数の変化を算出する対流場算出部14
0と、渦度と流れ関数の密度分布から流れ場の変動を予
測する対流度予測部150と、その予測結果を出力する
出力部160により構成されている。
【0014】移動速度検出部120は、2つの異なる時
刻の画像だけを用いても移動量の検出等の画像処理を行
うことができるが、実際の降水パターンが複雑に、非定
常的に変化しているため平均処理を含める必要があり、
具体的には、例えば、過去6回の画像を前半、後半の3
画像の平均を求め、更に前半、後半の値の平均を求め、
2つ以上の2次元画像間で画像処理を行って、画像特徴
量を抽出しながら画像中の物体の移動速度を検出する。
【0015】図2は局所的に検出された移動ベクトルか
ら渦度・流れ関数による対流場を予測する流れについて
画像の特徴を示した推移図である。図の200は、従来
から用いられている小ブロック相互相関法を気象レーダ
ーエコー画像(降水量・降水域)に適用して局所的な移
動ベクトルを検出する。図の210と220は、移動平
均フィルターにより、水平方向と垂直方向とに独立に局
所移動ベクトルを伝播させた結果である。図の230と
240は、適当な反復回数の後、流れ関数、渦度関数
(ψ、ζ)を変数とする式(1)と(2)の各画素
(u、v)の移動ベクトル値を当てはめながら、連立方
程式を解く。
【0016】
【数1】 上式(1)(2)は、差分法に従って、時間項は前進差
分、1次微分項は中心差分、2次微分は位置に関して2
次オーダーの差分公式で離散化される。即ち、図の23
0と240が示すように、時間積分を施すことで、流れ
関数と渦度の時間的な変化が予測される。
【0017】図3は、渦度と流れ関数の密度分布に従っ
て、湧きだし領域と吸い込み領域が予測された結果を示
す図である。渦度の密度が高い領域300は湧きだしの
候補領域となり、逆に密度が低い領域310は吸い込み
の候補領域となる。
【0018】図4は、実際の降水パターンでの3時間の
湧きだし領域と吸い込み領域の予測的中率を示す。な
お、的中率は、面積と位置についての関数である。その
結果、湧きだし領域よりも、吸い込み領域を予測するこ
とが難しいことがわかる。しかしながら、従来このよう
な予測は試みられていないため、本発明の対流場変化予
測推定装置は、対流の急激な変化による異常降水現象の
予測への適用が非常に有効であると言える。
【0019】
【発明の効果】以上説明したように、本発明によれば、
限られた画像情報から局所的に移動ベクトルを検出し、
ゼロ領域へ反復法により局所移動ベクトルを伝播させた
後、流れ場の移動ベクトルを渦度と流れ関数に置き換え
て、渦度と流れ関数を連立して解き、その時間的な推移
を求めながら、渦度と流れ関数の密度分布の関係から、
対流に伴った擾乱する領域を予測していくことができる
と言う効果がある。
【図面の簡単な説明】
【図1】本発明の実施例の構成を示すブロック図であ
る。
【図2】局所移動ベクトルから渦度・流れ関数の算出と
予測による推移図である。
【図3】渦度分布と湧き出し・吸い込みの候補領域の予
測出力図である。
【図4】湧き出し領域と吸い込み領域の予測の評価結果
を示す図である。
【符号の説明】
100 画像入力部 110 画像蓄積部 120 移動速度検出部 130 反復的移動速度伝播部 140 流れ場算出部 150 対流度予測部 160 出力部
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.7 識別記号 FI G06T 7/20 G06F 15/70 400 (56)参考文献 特開 平8−304561(JP,A) 特開 平7−110378(JP,A) 篠沢一彦、藤井雅晴、曽根原登,“局 所並列計算による降雪レーダ予測方式の 検討”,電子情報通信学会論文誌,日 本,社団法人電子情報通信学会,1995年 7月25日,第J78−D−2巻、第7 号,p.1144−1149 境野英朋、堀越力、鈴木智,“レーダ ーエコー画像における移流・拡散方程式 を用いた降水パターン変化予測方法”, 電子情報通信学会総合大会講演論文集 情報・システム2,日本,社団法人電子 情報通信学会,1997年 3月 6日,D −11−197,P.197 境野英朋、末永康仁、石井健一郎, “力学的相互作用における弾性移動体の 計算法に関する検討”,電子情報通信学 会秋期大会−ソサイエティ先行大会−講 演論文集 情報・システム,日本,社団 法人電子情報通信学会,1994年 9月 5日,D−367,p.375 (58)調査した分野(Int.Cl.7,DB名) G01W 1/00 - 1/18 G01P 3/36 G01S 13/95 G06T 1/00 G06T 7/20 JICSTファイル(JOIS)

Claims (4)

    (57)【特許請求の範囲】
  1. 【請求項1】 湧き出しと吸い込み領域を場として有す
    る対流場変化予測推定装置であって、 画像入力する画像入力手段と、 入力された画像を時系列画像として蓄積する画像蓄積手
    段と、 前記画像蓄積手段に蓄積されている連続する2つ以上の
    2次元画像間で画像処理を行って、画像特徴量を抽出し
    ながら画像中の物体の移動速度を検出する移動速度検出
    手段と、 前記画像特徴量から検出され複数の移動速度を初期移
    動速度として、移動平均フィルターを反復的に適用して
    初期移動速度をゼロ速度領域へ伝播させる反復的移動速
    度伝播手段と、 移動速度を初期値として、渦度と流れ関数を連立させ
    て、渦度と流れ関数の変化を算出する対流場算出手段
    と、 渦度と流れ関数の密度分布から流れ場の変動を予測する
    対流度予測手段と、 前記対流度予測手段の予測結果を出力する出力手段を有
    する対流場変化予測推定装置。
  2. 【請求項2】 前記画像入力手段が、 気象レーダーエコー画像とドップラー画像から観測され
    た情報を入力する画像入力手段を有する請求項1記載の
    対流場変化予測推定装置。
  3. 【請求項3】 前記移動速度検出手段が、 画像を複数の小領域に分割して、連続する画像の各小領
    域間での相互相関係数を解析して、最も相関係数の高い
    点を求めて前記小領域の移動ベクトルとして検出する移
    動速度検出手段を有する請求項1記載の対流場変化予測
    推定装置。
  4. 【請求項4】 前記移動速度検出手段が、 渦度と流れ関数の密度分布に基づいて、画像の濃淡値を
    変数にもつ流体方程式の湧きだし項と吸い込み項に特定
    な比例定数を含み、湧きだしと吸い込み領域の移動速度
    を検出する移動速度検出手段を有する請求項1記載の対
    流場変化予測推定装置。
JP17718397A 1997-07-02 1997-07-02 対流場変化予測推定装置 Expired - Fee Related JP3377078B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP17718397A JP3377078B2 (ja) 1997-07-02 1997-07-02 対流場変化予測推定装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP17718397A JP3377078B2 (ja) 1997-07-02 1997-07-02 対流場変化予測推定装置

Publications (2)

Publication Number Publication Date
JPH1123726A JPH1123726A (ja) 1999-01-29
JP3377078B2 true JP3377078B2 (ja) 2003-02-17

Family

ID=16026640

Family Applications (1)

Application Number Title Priority Date Filing Date
JP17718397A Expired - Fee Related JP3377078B2 (ja) 1997-07-02 1997-07-02 対流場変化予測推定装置

Country Status (1)

Country Link
JP (1) JP3377078B2 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4538426B2 (ja) * 2006-04-12 2010-09-08 日本電信電話株式会社 移動ベクトル検出装置、移動ベクトル検出方法および移動ベクトル検出プログラム
JP4499775B2 (ja) * 2007-10-31 2010-07-07 日本電信電話株式会社 パターン予測装置、パターン予測方法及びパターン予測プログラム
CN102435411B (zh) * 2011-09-05 2013-11-27 中国人民解放军国防科学技术大学 Npls可压缩湍流雷诺应力的全场测量系统及方法
US11237299B2 (en) * 2017-05-01 2022-02-01 I.M. Systems Group, Inc. Self-learning nowcast system for modeling, recording, and predicting convective weather

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
境野英朋、堀越力、鈴木智,"レーダーエコー画像における移流・拡散方程式を用いた降水パターン変化予測方法",電子情報通信学会総合大会講演論文集 情報・システム2,日本,社団法人電子情報通信学会,1997年 3月 6日,D−11−197,P.197
境野英朋、末永康仁、石井健一郎,"力学的相互作用における弾性移動体の計算法に関する検討",電子情報通信学会秋期大会−ソサイエティ先行大会−講演論文集 情報・システム,日本,社団法人電子情報通信学会,1994年 9月 5日,D−367,p.375
篠沢一彦、藤井雅晴、曽根原登,"局所並列計算による降雪レーダ予測方式の検討",電子情報通信学会論文誌,日本,社団法人電子情報通信学会,1995年 7月25日,第J78−D−2巻、第7号,p.1144−1149

Also Published As

Publication number Publication date
JPH1123726A (ja) 1999-01-29

Similar Documents

Publication Publication Date Title
Salhi et al. Object tracking system using Camshift, Meanshift and Kalman filter
US8885876B2 (en) Visual tracking system and method thereof
US7536035B2 (en) Object velocity measuring apparatus and object velocity measuring method
US6628715B1 (en) Method and apparatus for estimating optical flow
Badenas et al. Motion-based segmentation and region tracking in image sequences
GB2320830A (en) Rain tracking radar
Han et al. Visual tracking by continuous density propagation in sequential Bayesian filtering framework
CN103604947B (zh) 一种时间分辨率自适应调整的流场状态测量方法
CN104219488A (zh) 目标图像的生成方法和装置、以及视频监控系统
Bae et al. Image tracking algorithm using template matching and PSNF-m
Moreno-Garcia et al. Video sequence motion tracking by fuzzification techniques
JP3377078B2 (ja) 対流場変化予測推定装置
JP4210292B2 (ja) 画像変化予測方法および画像変化予測装置
JP3953350B2 (ja) 映像シーン予測方法、装置、プログラム、および記録媒体
KR100635883B1 (ko) 실시간 물체 추적 장치
Richefeu et al. A new hybrid differential filter for motion detection
JPH1114760A (ja) 気象画像予測方法、装置および気象画像予測プログラムを記録した記録媒体
Srilekha et al. A novel approach for detection and tracking of vehicles using Kalman filter
JP2000348180A (ja) 移動領域追跡方法及び装置及び記録媒体
JP3377079B2 (ja) 流速推定装置
JP3351460B2 (ja) 移動体の挙動解析予測装置およびそのうねり対応づけ処理方法
JP3351461B2 (ja) 降水量最尤推定装置
Abbas et al. Real time fuzzy based traffic flow estimation and analysis
JP3390975B2 (ja) 画像流速推定方法、装置および画像流速推定プログラムを記録した記録媒体
JP3458940B2 (ja) 気象画像予測方法、装置および気象画像予測プログラムを記録した記録媒体

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20071206

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20081206

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091206

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101206

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101206

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111206

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111206

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121206

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121206

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131206

Year of fee payment: 11

LAPS Cancellation because of no payment of annual fees