1.一种雷达回波零度层亮带智能识别预警方法,其特征在于:包括如下具体步骤:步骤1:读取一个雷达基数据文件,从该文件中提取基本反射率信息,记为其中, 表示雷达探测的仰角,γ[j]∈{γ_min,…,γ_max}表示雷达探测点到雷达中心点的距离,ω[k]∈{ω_min,…,ω_max}表示雷达探测点的方位角;
步骤2:定义一个集合H,将雷达探测区间按垂直高度划分为若干等份,每份垂直高度定义为Alt_N,集合H表示为:H∈{Alt_min,Alt_min+Alt_N,Alt_min+Alt_N×2,…,Alt_max}其中,Alt_min和Alt_max分别表示雷达探测的最小垂直高度和最大垂直高度;Alt_N为经验阈值,Alt_N为整数,取值范围为[5,20];
步骤3)定义数组VA[i],依次记录与探测点 具有相同方位角且位于同一垂线上的探测点的基本反射率数据,其具体步骤为:
步骤3.1)定义数组VA[i],令 中i=0,j=0,k=0,将
存储到VA[i],记为VA[0];
步骤3.2)令i不断自增1至仰角 找到所有与探测点 具有
相同方位角且位于同一垂线上的探测点的基本反射率数据,记为 并
将探测点的基本反射率数值存储到VA[i],其中, 的计算方法为:
步骤4)将步骤3.2)中数组VA[i]中数值最大的元素记为VA[i_max],记录其在数组中的位置为i_max,然后判断VA[i_max]是否满足以下公式1,如果满足,则进入步骤5);反之,则进入步骤8):式中,VA[i_max-1]和VA[i_max+1]分别表示数组VA中第i_max-1个和第i_max+1个元素; 表示集合 中元素的数量,即不同仰角的数量;Thd_A是一个经验阈值,取值范围为[5,20],单位为dBZ;
步骤5)存储VA[i]中数值最大的元素在数组中的位置i_max,清空数组VA[i],将步骤
3.1)中i和j重置为0,k自增1,即k=1,并将 存储到VA[i],记为VA[0];再按照步骤3.2)的方法,依次记录与探测点 具有相同方位角且位于同一垂线上的探测点的基本反射率数据,重新存储至数组VA[i];
步骤6)定义变量A(γ[j])和C(γ[j]),其中,A(γ[j])和C(γ[j])的初始值均为0;
取上步中的i_max,并以该i_max的值判断当前数组VA[i]是否满足公式1,如果满足,则将变量A(γ[j])和C(γ[j])分别自增1;反之,则仅将变量C(γ[j])自增1,A(γ[j])的值保持不变;
步骤7)定义二维数组E(δ,h),该二维数组中的各个元素初始值均为-1;
其中,δ∈{γ_min,…,γ_max};h∈H,h的计算方法包括如下步骤:首先,计算探测点 的垂直高度Vh,计算方法为:
然后,确定上述垂直高度Vh在步骤2中阈值范围H中所属区间,并将该区间的下界值赋给h;
最后采用步骤5和步骤6的方法,令k不断自增1至方位角ω[k]=ω_max,分析变量A(γ[j])和C(γ[j])是否满足以下公式2,如果满足,则将变量E(δ,h)标记为1;反之,则将变量E(δ,h)标记为0,得到距离库长γ[0]对应的判定值E(γ[0],h):其中,Thd_B是一个经验阈值,取值范围为[0.4,1];
步骤8)将步骤3.1)中i和k重置为0,j不断自增1,采用步骤3~步骤7的方法,直到γ[j]=γ_max,得到所有距离库长对应的判定值E(γ[j],h),即E(δ,h);
步骤9)以列为单位,找出步骤8中所有距离库长对应的判定值E(δ,h)中每一列元素中取值连续为1的且连续为1长度最长的一段序列,并将该序列长度记为Len_max(h);
再将所有满足Len_max(h)>Thd_C的序列长度Len_max(h)按大小进行排序,从中找出最大的前Trd_D个列,并记录下这些列的h,记为RH,RH是一个数组,该数组中的元素项记为RH(s),其中,s为自然数;Trd_C是一个经验阈值,取值范围为[5,MaxDis],MaxDis表示雷达探测的最大距离;Trd_D为一个经验阈值,取值范围为[1,4];
步骤10)如果RH(s)中没有任何高度信息,则可判定步骤1所取的雷达基数据文件中不存在零度层亮带;如果RH(s)中存在1个或多个高度信息,那么在这些高度层上,可判定为疑似存在零度层亮带。
2.根据权利要求1所述的雷达回波零度层亮带智能识别预警方法,其特征在于:对于所述疑似存在零度层亮带的高度层是否确定存在零度层亮带的判定方法如下:利用同一雷达上一时刻的天气雷达基数据文件,采用步骤1-步骤10的方法判定得到同一雷达上一时刻对应的RH’,如果RH’在RH疑似存在零度层亮带的高度层也疑似存在零度层亮带,则该高度层确定存在零度层亮带;反之,如果RH’在RH疑似存在零度层亮带的高度层不存在零度层亮带,则该高度层疑似存在零度层亮带。
3.根据权利要求2所述的雷达回波零度层亮带智能识别预警方法,其特征在于:对于所述疑似存在零度层亮带的高度层是否确定存在零度层亮带的判定方法如下:读取与当前雷达基数据文件探测时间相近的常规探空观测资料或数值模式资料,从这些资料中提取与步骤10中RH(s)中疑似存在零度层亮带的高度层的温度信息,记为DT(RH(s)),分析RH(s)是否满足公式3,如果满足,则判定当前雷达探测到在RH(s)中疑似存在零度层亮带的高度层确定存在零度层亮带;反之,在RH(s)中疑似存在零度层亮带的高度层不存在零度层亮带:Trd_E≤DT(RH(s))≤Thd_F 公式3其中,Trd_E和Thd_F均为经验阈值,Trd_E∈[-20℃,0℃],Trd_F∈[0℃,10℃]。
4.根据权利要求3所述的雷达回波零度层亮带智能识别预警方法,其特征在于:将疑似存在的零度层亮带或确定存在零度层亮带呈现给业务人员,用于雷达回波资料的零度层亮带智能识别和预警。