专利名称：多普勒雷达遥感强风暴的预警方法
技术领域：
本发明属于大气环境的遥感监测技术，具体涉及一种多普勒雷达遥感强风暴的自动预警方法。
背景技术：
强对流天气，也称强风暴，是主要的灾害性天气之一，包括雷电、冰雹、雷雨、大风等，其发生时会产生巨大的社会�：ΑＲ虼耍郧糠绫┑募嗖夂驮ぞ哂邢灾纳缁岷途靡庖�。由于强风暴空间尺度较小、持续时间较短，运用常规观测手段难以对其进行有效监测。而多普勒天气雷达目前是强风暴的主要遥感监测手段，具有很高的时空分辨率，并且可以观测风暴的三维空间结构。强风暴在天气雷达上以三维回波图像的形式出现。当回波分布稀疏，且变化缓慢时，较容易进行预警。但对回波分布密集的风暴簇，或回波变化剧烈，频繁出现分裂与合并等复杂情况，会给预警带来很大困难，现有的方法无法处理。但随着气候变暖，极端强对流灾害天气频发，这些回波分布复杂情况的出现也愈发频繁。因此，针对这些回波分布密集、分裂合并频繁等复杂情况，提出能够处理这些问题的风暴自动预警方法是十分重要的。
发明内容
本发明的目的，是提供一种强风暴的自动预警方法，以弥补现有技术的不足。本发明首先在风暴的识别上，采用基于数学形态学的方法，可以较好地解决虚假合并问题，并能从复杂的风暴簇中分离出风暴单体。然后，在风暴的追踪上，引入序贯蒙特卡洛方法，使得追踪过程大为简化。不仅可以在追踪的过程中同时处理风暴的分裂与合并，还可以处理风暴的漏测情况。最后，对风暴的预警，结合光流法得到的运动矢量场，进行线性拟合外推预警。本发明的强风暴预警方法步骤如下I.数据预处理将质量控制后的雷达基数据，插值到三维地理坐标下。2.强风暴的识别本发明提出了基于数学形态学的风暴识别方法，首先使用第一级阈值进行单阈值识别；其次，对识别得到的风暴执行腐蚀操作，以消除虚假合并；然后，使用高一级阈值进行识别，并对识别得到的风暴.体进行膨胀操作，当风暴的边界在膨胀的过程中相互之间接触，或接触到了原来较低阈值识别的风暴的边界时，则停止膨胀过程；最后，逐次使用更高级别的阈值进行识别，并在每一个级阈值的识别过程中执行腐蚀和膨胀操作。3.强风暴的追踪本发明将序贯蒙特卡洛方法应用于风暴的追踪，其基本过程是一个随着时间序列进行迭代的过程，每次迭代包括三步抽样，预测和测量，从而实时地获得所需要的动态息。本发明对相邻时刻风暴的追踪采用雷达组合反射率因子图像。抽样所得的每一个粒子都要通过一阶线性系统模型，进行一步预测。
(I)相邻时刻风暴的追踪。定义^为事件t时刻的第m个风暴由t-ι时刻的第k个风暴发展而来。当观测Yt到达后，计算事件发生的概率Pierk IY,) == = I-JVm
mk其中，Am是t时刻第m个风暴的面积，Ak是t_l时刻第k个风暴的面积，NFm, k是指从t-Ι时刻第k个风暴抽样出来的所有粒子中，经一步预测之后落入t时刻第m个风暴的粒子数目，如果I X)大于阈值 ；，本发明设为O. 5，则t时刻的第m个风暴被认为是由t-i时刻的第k个风暴发展而来的。(2)风暴漏测的处理。使用延迟逻辑来解决这个问题，即在当前时刻不做出决策，而是延迟一个或多个时刻，直到有充分多的信息时才做出决策。首先，在t-2时刻，对风暴 做均匀抽样并进行一步预测；然后，在t-Ι时刻，对一步预测得到的新的粒子更新权重w，方法如下j = I.. NP，其中，I」,η是在粒子j所处位置实际观测到的反射率因
rmm
子的值，NP是在t-2时刻从风暴中抽样的粒子数目，Tzmin是风暴算法中的第I级反射率因子阈值。正规化这些粒子的权重，使得1二 wV,-. =1。按照= Σ: MVlXy来估计的风暴虚拟质心位置。然后，对粒子集
进行重抽样并进行一步预测，根据丨Y,)即可判断t-2，t_l，和t时刻风暴的关系，即是否出现了漏测。(3)分裂、合并的处理。根据Ρ(Τ· Y,)的值，并从粒子与风暴边界的相对位置，容易判断如果t时刻的两个风暴都是由t-i时刻的同一个风暴发展而来的，则认为t-i时刻的这个风暴发生了分裂；类似地，可以判别合并的情况4.强风暴的预警首先，对最近时刻及前一时刻的雷达数据，计算其组合反射率因子图像，使用光流法计算得到运动矢量�。蝗缓蠖悦恳桓鍪侗鸬姆绫扑闫渌哺乔虻钠骄硕噶�(VaJ ;最后，使用V·和风暴的质心位置，预警未来时刻的风暴位置；对风暴的顶高、体积等属性参数，使用一阶线性模型预警。5.结果输出不难将以上方法编制成软件系统在计算机上，以图像形式显示出来。本发明的优点可以有效的消除风暴的虚假合并，而且能够在回波分布较密集的情况下正确地识别风暴；在风暴追踪过程中充分考虑了可能的风暴漏测情况，有效解决了风暴的分裂与合并问题；结合使用回波运动矢量场和一阶线性模型进行预警，得到更稳定的结果。
图I是本发明的强风暴的预警方法的流程示意图。
具体实施例方式如图I是本发明的强风暴预警方法的实现步骤如下
I.数据预处理。首先进行质量控制，去除质量有严重问题的雷达基数据。然后采用径向和方位上的最近邻居法和垂直线性内插法，将极坐标系下的基数据插值到三维地理坐标下。2.强风暴的识别。主要步骤如下I)使用第一级阈值Tz,min，使用三维聚类法进行风暴的识别。2)对第I)步识别得到的三维风暴进行一次腐蚀操作，消除相邻风暴之间的微弱连接，从而解决虚假合并的问题。3)依次使用阈值Tza = Tz,min+(i_l) X5dBZ, i = 2…Nttoesh，进行风暴的识别，并对识别的结果进行腐蚀和膨胀操作，从而逐级地分离出风暴簇中强弱不均的各个风暴单体。3.强风暴的追踪将序贯蒙特卡洛方法的原理应用于风暴的追踪，其基本过程是一个随着时间序列进行迭代的过程，每次迭代包括三步抽样，预测和测量。(I)对相邻时刻风暴的追足示a)抽样。对t-Ι时刻的每一个被识别出来的风暴，在其覆盖的区域内进行重要性抽样，抽样时考虑每一个样本所处位置处的反射率因子的大�。檠牧Ｗ邮狽Pk = floor (Ak), k = I. . . N^1其中，Ak为第k个风暴的投影面积，大小在30 1000km2。Nw为t_l时刻检测到的风暴总数。Floor(Ak)是指取不超过Ak的最大整数。对每一个粒子状态以它所在的位置坐标X = (X，y)表示。同时，为每个粒子分配其所属风暴的标识号和一个相同的权重W。暴体的标识号用来在后继的预测和测量阶段区分属于不同风暴的粒子。对t-i时刻的第k个风暴，从中抽样得到的粒子集合可以表示为
权利要求
1.一种强风暴的自动预警方法，其步骤如下 (1)、数据预处理首先进行质量控制，去除质量有严重问题的雷达基数据，然后采用径向和方位上的最近邻居法和垂直线性内插法，将极坐标系下的基数据插值到三维地理坐标下； (2)、强风暴的识别本发明提出了基于数学形态学的风暴识别方法，首先使用第一级阈值进行单阈值识别；其次，对识别得到的风暴执行腐蚀操作，以消除虚假合并；然后，使用高一级阈值进行识别，并对识别得到的风暴.体进行膨胀操作，当风暴的边界在膨胀的过程中相互之间接触，或接触到了原来较低阈值识别的风暴的边界时，则停止膨胀过程；最后，逐次使用更高级别的阈值进行识别，并在每一个级阈值的识别过程中执行腐蚀和膨胀操作； (3)、强风暴的追踪本发明将序贯蒙特卡洛方法应用于风暴的追踪，其基本过程是一个随着时间序列进行迭代的过程，每次迭代包括三步抽样，预测和测量，从而实时地获得所需要的动态信息；本发明对相邻时刻风暴的追踪采用雷达组合反射率因子图像，抽样所得的每一个粒子都要通过一阶线性系统模型，进行一步预测； 相邻时刻风暴的追踪定义^^为事件t时刻的第m个风暴由t-Ι时刻的第k个风暴发展而来，当观测Yt到达后，计算事件发生的概率
2.如权利要求I所述的强风暴预警方法，其特征是步骤(2)中风暴的三维识别采用了基于数学形态学的方法，即通过在多阈值识别过程中，综合使用膨胀和腐蚀操作，可以成功地识别出风暴的虚假合并,并从风暴簇中分离出相距较近的风暴。
3.如权利要求I所述的强风暴预警方法，其特征是步骤(3)中使用了序贯蒙特卡洛方法，使用ρ(6Γ Y,)作为关键指标，来判断不同时刻风暴的对应关系，简化了追踪过程，可以同时处理风暴的分裂与合并，还可以处理风暴的漏测情况。
4.如权利要求I所述的强风暴预警方法，其特征是步骤(4)中结合光流法得到的运动矢量�。邢咝阅夂贤馔圃ぞ�。
全文摘要
本发明涉及一种强风暴的自动预警方法。本发明利用多普勒天气雷达数据。首先在风暴的识别上，采用基于数学形态学的方法，即通过在多阈值识别过程中，综合使用膨胀和腐蚀操作，可以成功地识别出风暴的虚假合并，并从风暴簇中分离出相距较近的风暴。然后，在风暴的追踪上，使用序贯蒙特卡洛方法，使得追踪过程大为简化，还可以处理风暴的分裂合并，和风暴漏测的情况。最后，对风暴的预警，结合光流法得到的运动矢量�。邢咝阅夂贤馔圃ぞ�。本发明可以处理雷达回波分布密集、分裂合并频繁等复杂情况。本发明的用途，是在各省市气象台安装预警软件系统，提供强对流天气的预警服务，满足国家防灾减灾的迫切需要。
文档编号G01W1/10GK102721987SQ20121019092
公开日2012年10月10日 申请日期2012年6月12日 优先权日2012年6月12日
发明者卞继方, 宋珍妮, 韩雷 申请人:中国海洋大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
