[发明专利]一种强对流单体样本的标记方法

权利要求书

1.一种强对流单体样本的标记方法，其特征在于，包括以下步骤：

步骤一、输入对流天气实况和一组连续体扫的雷达基数据，所述对流天气实况包括对流天气类型、发生结束时间、发生结束地点；

以所述该组雷达基数据为原始数据，求取获得一组组合反射率图，记为第一组合反射率图；同时获得一组去除了0.5°仰角反射率后得到的组合反射率图，记为第二组合反射率图；并对上述得到的两组组合反射率图分别进行单体分割，从而得到单体；

步骤二、对存在于第一组合反射率图而不存在于第二组合反射率图的单体进行超折射滤除，并以第一组合反射率图中保留的单体作为最终的单体分割结果；

步骤三、通过光流法实现前后体扫时刻的单体跟踪，从而确定前后体扫时刻的单体是否关联；

步骤四、对强对流单体进行标记：

用户根据对流天气实况对一个体扫时刻的所有单体分别进行类型的标记，根据步骤三的跟踪结果自动给出与所标记单体相关联的所有单体的类型的标记，单体的类型包括强对流单体类型和非强对流单体类型，其中，所述强对流单体的类型包括：冰雹、短时强降水、雷雨大风、冰雹+短时强降水、冰雹+雷雨大风和冰雹+短时强降水+雷雨大风；

将自动给出的标记结果中不符合对流天气实况的标记结果进行手动更正。

2.根据权利要求1所述强对流单体样本的标记方法，其特征在于：步骤一中，单体分割的步骤如下：

步骤1)筛选n个种子区域：将以45dBZ为阈值，以80km²为面积阈值提取的连通域作为种子区域，分别标记为S₁、S₂、……、S_n-1、S_n；

步骤2)筛选最大关注区域：用阈值分割法得到所有取值大于30dBZ的区域作为最大关注区域，分别标记为N₁、N₂、……、N_m-1、N_m，N_i中包含种子区域的个数为j，j≤n，1≤i≤m，若j≤1，则最大关注区域N_i表示为一个单体；否则，对最大关注区域N_i进行如下操作；

步骤2-1)令该最大关注区域N_i中的种子区域长大形成k个单体核H，包括：

从最大关注区域N_i的外围历经30dBZ、35dBZ和40dBZ三个值域到达种子区域边界，将其中阈值在40dBZ、35dBZ和30dBZ以上的区域分别记为种子区域的内核、中核和外核；

单体核H的形成根据以下情况确定：

1)若最大关注区域N_i的外核包含的种子区域数目小于2，则该最大关注区域N_i的外核形成了单体核H；

2)若最大关注区域N_i的外核包含的种子区域数目不小于2，且最大关注区域N_i的中核包含的种子区域数目小于2，则该最大关注区域N_i的中核形成了单体核H；

3)在上述情况1)和2)均不成立的情况下，则该最大关注区域N_i的内核形成了单体核H；

步骤2-2)膨胀-避让式分割单体：

①最大关注区域N_i的最小外包矩形宽为w、高为h，将k个单体核H记为H₁、H₂、......、H_k-1、H_k，设一个数组A的大小为w×h，并设该数组A所有元素的初值为0；

②用数组A将最大关注区域N_i中的所有像素点做确定性标记和尚不确定性标记：即将数组A中与第1个单体核对应的元素赋值为1，以此表示此像素点属于第1个单体；将数组A中与最大关注区域N_i中非单体核区域对应的元素赋值为k+1；

③用数学形态学的膨胀法对最大关注区域N_i中的所有单体核H₁、H₂、......、H_k-1、H_k进行膨胀，得到膨胀模板H^*₁、H^*₂、......、H^*_k-1、H^*_k；遍历中的新增点P，p＝1，2，3，....，k-1,k；若A(P)＝m+1，则将P点拾入单体p，即A(P)＝p；

④逆转p的顺序，重复上述③的操作直至数组A中的所有值均不为m+1。

3.根据权利要求1所述强对流单体样本的标记方法，其特征在于：步骤三的具体内容如下：

步骤1)利用光流法求得光流数：利用步骤二得到前体扫时刻单体和后体扫时刻单体；然后输入该两帧组合反射率图，并通过光流法计算出光流场数据；遍历光流场数据点，对于X方向速度或者Y方向速度大于30的光流场数据点赋值为0；通过光流场数据计算得到前体扫时刻所有像素在后体扫时刻的流入位置；对比前后体扫时刻各单体的范围，计算出光流数；光流数记为SUM(A→B)，用以表示由前体扫时刻单体A流入后体扫时刻单体B的像素数；用光流数量化前体扫时刻单体A和后体扫时刻单体B的关联程度，

步骤2)结合前后体扫时刻单体的面积S，综合考虑前后体扫时刻单体出现分裂或合并的情况，制定如下三个条件，若满足其一则认为前体扫时刻单体A和后体扫时刻单体B存在关联，

条件一、SUM(A→B)>S(A)×50％；

条件二、SUM(A→B)>S(B)×50％；

条件三、设SUM(A→B)＝c，对前体扫时刻的l个单体流入后体扫时刻单体B的所有像素求和，记为SUM，c/SUM>50％，且c/S(B)>33％。