[发明专利]一种对流天气的临近预报方法及系统有效
申请号: | 201510368821.5 | 申请日: | 2015-06-29 |
公开(公告)号: | CN104977584B | 公开(公告)日: | 2017-08-15 |
发明(设计)人: | 曹春燕;陈元昭;张文海;陈训来;徐文文;李程;刘东华 | 申请(专利权)人: | 深圳市气象台 |
主分类号: | G01S13/95 | 分类号: | G01S13/95 |
代理公司: | 深圳中一专利商标事务所44237 | 代理人: | 张全文 |
地址: | 518040 广东省深*** | 国省代码: | 广东;44 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 对流 天气 临近 预报 方法 系统 | ||
1.一种对流天气的临近预报方法,其特征在于,包括:
获取探测区域内的雷达拼图;
基于所述雷达拼图,采用局部约束的光流法计算得到雷达回波的光流场,对于所述光流场的光流(u,v),所述u,v为所述雷达拼图其中一点的光流在x,y方向上的分量,该点在时刻t的灰度值为I(x,y,t),相应的光流约束方程为:
Ixu+Iyv+It=0,
其中:若该点在(t+Δt)时运动到(x+Δx,y+Δy),引入进一步的约束条件来求解所述光流约束方程:
依据所述雷达回波的光流场进行雷达回波外推,完成对流天气的临近预报。
2.如权利要求1所述的方法,其特征在于,在所述获取探测区域内的雷达拼图之后,所述基于所述雷达拼图,采用局部约束的光流法计算得到雷达回波的光流场之前,所述方法还包括:
采用中值滤波对所述雷达拼图的数据做平滑处理。
3.如权利要求2所述的方法,其特征在于,所述采用中值滤波对所述雷达拼图的数据做平滑处理包括:
采用预设结构的二维滑动模板,将板内像素按照像素值的大小进行排序,生成单调上升或下降的二维数据序列g(x,y)=med{f(x-k,y-1),(k,1∈ w)},其中,所述f(x,y),和g(x,y)分别为原始的所述雷达拼图的数据和处理后的所述雷达拼图的数据,所述w为所述预设结构的二维滑动模板。
4.如权利要求3所述的方法,其特征在于,所述预设结构的二维滑动模板的形状包括线状,圆形,十字形或圆环形。
5.如权利要求1所述的方法,其特征在于,所述基于所述雷达拼图,采用局部约束的光流法计算得到雷达回波的光流场还包括:
若所述雷达拼图某点的光流矢量与其相邻点的光流矢量平均值的偏差或夹角超过预设阈值,则在所述雷达回波的光流场中剔除该点的光流矢量。
6.一种对流天气的临近预报系统,其特征在于,所述系统为基于地理信息系统的可视化平台,包括服务器和客户端,所述客户端用于显示所述服务器的临近预报结果,所述服务器包括:
获取单元,用于获取探测区域内的雷达拼图;
运算单元,用于基于所述雷达拼图,采用局部约束的光流法计算得到雷达回波的光流场,对于所述光流场的光流(u,v),所述u,v为所述雷达拼图其中一点的光流在x,y方向上的分量,该点在时刻t的灰度值为I(x,y,t),相应的光流约束方程为:
Ixu+Iyv+It=0,
其中:若该点在(t+Δt)时运动到(x+Δx,y+Δy),引入进一步的约束条件来求解所述光流约束方程:
预报单元,用于依据所述雷达回波的光流场进行雷达回波外推,完成对流天气的临近预报。
7.如权利要求6所述的系统,其特征在于,所述服务器还包括:
中值滤波单元,用于采用中值滤波对所述雷达拼图的数据做平滑处理。
8.如权利要求7所述的系统,其特征在于,所述中值滤波单元具体用于:
采用预设结构的二维滑动模板,将板内像素按照像素值的大小进行排序,生成单调上升或下降的二维数据序列g(x,y)=med{f(x-k,y-1),(k,1∈w),其中,所述f(x,y),和g(x,y)分别为原始的所述雷达拼图的数据和处理后的所述雷达拼图的数据,所述w为所述预设结构的二维滑动模板。
9.如权利要求8所述的系统,其特征在于,所述预设结构的二维滑动模板的形状包括线状,圆形,十字形或圆环形。
10.如权利要求6所述的系统,其特征在于,所述运算单元还用于:
若所述雷达拼图某点的光流矢量与其相邻点的光流矢量平均值的偏差或夹角超过预设阈值,则在所述雷达回波的光流场中剔除该点的光流矢量。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于深圳市气象台,未经深圳市气象台许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201510368821.5/1.html,转载请声明来源钻瓜专利网。
- 上一篇:无线通信同步系统
- 下一篇:用于一调频连续波雷达感测系统的信号处理方法及装置