————————————
作者简介:陈晓静(1983-),女,硕士,主研方向:图像处理,气象雷达模拟仿真;敬忠良,教授、博士生导师;张 军,博士后
收稿日期:2012-08-27 修回日期:2013-01-03 E-mail:fanzhongyan1000@163.com
基于卫星图像数据的大范围风场重建
陈晓静,敬忠良,张
军
(上海交通大学航空航天学院,上海 200240)
摘 要:在机载气象雷达早期开发中,数值天气模型是一个非常重要的环节,而天气模型中风场的准确性直接影响整个模型的准
确性。为解决数值天气模型中的连续图像运动分析问题,提出一种基于运动场 Helmholtz 分解的低维流体运动估计方法。通过少量
涡流粒子和源粒子的演化,得到光流场的低维参数化表达,即非旋转量和螺旋量基函数的线性组合,其中基函数由格林核梯度构
成,系数值和基函数参数通过最小化代价函数获得。实验结果表明,与传统光流法相比,该方法的计算速度快了近 4 倍,风场更
为准确地反映了实际的天气状况,在风场反演中更为可靠。
关键词:光流;图像运动分析;运动估计;机载气象雷达;数值天气;风场重建
Large-scale Wind Field Reconstruction Based on Satellite Image Data
CHEN Xiao-jing, JING Zhong-liang, ZHANG Jun
( School of Aeronautics and Astronautics, Shanghai Jiaotong University, Shanghai 200240, China)
【Abstract】In the early development of airborne meteorological radar, numerical weather is a very important part, while the accuracy of
wind field directly affects the accuracy of the entire numerical weather model. Aiming at sloving the motion analysis of successive images
in numerical weather model, low-dimensional fluid motion estimation method is used, and it is based on the Helmholtz decomposition of
motion field. Through the deformation of a small number of vortex and source particals, it can get the low-dimensional parametric
expression of optical flow field. The optical flow field consists of linear combinations of irrotational and solenoidal basis functions, which
are based on Green kernel gradient. The coefficient values and the basis function parameters are obtained by minimization of a function.
The experimental results show that compared with the traditional optical flow method, this method is nearly 4 times faster, and wind field
more accurately reflects the actual weather conditions. The method is more reliable in wind field reconstruction.
【Key words】
optical flow; image motion analysis; motion estimation; airborne meteorological radar; numerical weather; wind field
reconstruction
DOI: 10.3969/j.issn.1000-3428.2013.12.049
计 算 机 工 程
Computer Engineering
第 39 卷 第 12 期
Vol.39 No.12
2013 年 12 月
December 2013
文章编号:1000—3428(2013)12—0233—04
1
概述
数值天气在机载雷达研发模拟阶段至关重要,它作为
数据激励为机载气象雷达提供外部激励环境。数值天气反
演是根据飞行前方的连续卫星图片,通过数值计算,求解
描写天气演变过程的流体力学和热力学方程组,并通过合
适的云场模型来动态反演所在区域的天气信息
(
云场、风
场
)
,并把这些信息实时传送给雷达模拟器作为输入,雷达
对获取的信息进行处理,显示出前方天气信息。通过判断
雷达识别的天气信息与天气模型模拟的实际天气做对比,
来衡量雷达信号处理的准确性。在数值天气模型中,主要
包括风场重建、云微物理模型和云场动态模拟。本文的研
究重点在于利用卫星图片和一种基于光流的低微流体运动
估计器来重建风场。
光流法在运动图像分析、向量场分析中具有非常重要
的作用。目前常用的光流法是在基于
Horn
工作的基础上,
一种是对光流场加约束,改进光流的平滑;另一种是对灰
度场加以约束
[1]
。然而,针对卫星云图的运动场
(
风场
)
估计
的特殊问题,光流法存在一些问题:抗噪性不够好,计算
量很大,多数算法计算复杂耗时;图像边界处,图像序列
间没有相互匹配的部分
[2]
,不再满足光流约束方程,因此,
边界处计算的运动矢量与实际的相差比较大
[3]
。
鉴于传统光流的种种缺点,
Anne Cuzol
、
Pierre Hellier
等人提出了低维流体运动估计方法。这种方法是基于运动
场
Helmholtz
分解的,把向量场分为自由散度分量和自由涡
度分量,目的在于提供一个光流的低维参数化表达,通过
少量的涡流粒子和源粒子的变形来描述流体的运动。
2
种分
量都通过正则狄拉克测度对涡度和散度图离散化来近似。
结果得到非旋转量和螺旋量基函数的线性组合,这些基函
数通过格林核梯度和涡度图或者散度图做卷积分别得到。
系数值和基函数参数通过最小化一个函数得到,而这个函
数依赖于流体力学质量守恒原理的积分形式
[4]
。