风暴单体识别与追踪(完成版)

1、汇报人：南京(nn jn)信息工程大学*2015年12月16日降水系统(xtng)中心识别与追踪1共二十八页Content01020304概述(i sh) SCIT算法(sun f)原理 SCIT算法性能评估 总结 2共二十八页一、概述(i sh) 基于(jy)雷达数据的风暴体识别、追踪及预警方法是最早出现的临近预报技术,也是天气雷达系统和强天气预警业务的基本组成部分。风暴体识别、追踪及预警方法可以分为三大类 : 持续性预报法( Persistence)、交叉相关法( Cross Correlation)和单体质心法( Centroid Tracking) 。三者都属于外推预报法( Extra

2、polation)。其中 ,持续性预报法目前已经被后两者取代 。3共二十八页简单地说, 交叉相关法就是把整个数据(shj)区域划分成若干小区域, 相邻时刻雷达回波图像的小区域之间计算相关系数 ,通过最大相关系数确定相邻时刻图像中的区域对应关系,进而确定回波区域的平均运动。 1. 交叉(jioch)相关法交叉相关法示意图 对两幅相邻时刻( t 1 , t 2 ) 的雷达反射率因子回波图像 ,分别划分成若干大小相同的子区域。对 t 1 时刻的一个子区域 A , 在 t 2时刻回波图像的候选区域( 大小由平均风速决定) 中, 将所有可能的子区域分别与 A 做相关计算 。t 2 时刻相关系数最大的子区

3、域 B即为与 A 匹配的子区域。将 A 和 B 的中心位置连接起来 , 即为回波移动矢量 。得到回波移动矢量之后,就可以进行外推预警。4共二十八页优点(yudin)： 算法简单； 对较大的回波,交叉相关法有较好的识别和追踪效果 。缺点： 只考虑了回波的水平移动,计算量大等； 另外, 对较小的单体 ,尤其是相互距离较近的多个单体 , 其识别和追踪效果较差 。1.1 交叉(jioch)相关法的优缺点：5共二十八页单体质心法首先要识别出单体,并计算其质心位置、体积和投影面积等特征,然后在前后两个时刻的扫描数据中进行单体匹配追踪,最后通过连续(linx)多时刻的匹配追踪结果外推预警。对风暴体的追踪 ,

4、 各种单体质心法的具体实现步骤不一样, 但一般都会利用后面的预警结果来辅助追踪。2. 单体质(tzh)心法6共二十八页单体质心法能较好地识别和追踪较小的孤立单体 , 并且能够提供单体的更详细的特征数据 ；与交叉相关法相比 ,单体质心法计算量小,能够更有效地识别和跟踪风暴单体( 如飑线中的各个单体) ；由于质心法将风暴单体都分离出来 ,因此(ync)可以计算每个单体的特征 ,并分析其演变过程 ,从而为人们提供了一个分析风暴的工具。2.1 单体质(tzh)心法特点7共二十八页1999 年 ,Wolfson 提出了一个“风暴增长衰减追踪器”，结合 ITWS 中的交叉相关追踪器 ,可以给出风暴系统的增

5、长和衰减趋势； 2002 年, Handwerker 开发了一个三维风暴跟踪器( TRACE3D)； 2003 年 , Lakshmanan 提出了一种新的风暴识别方法。他使用了 K 均值聚类方法,可以根据需要进行不同尺度的风暴识别 、追踪和预警；2005 年 Xu和2006 年Lakshmanan 在统计的框架(kun ji)内, 使用了 Bayes 学习、Gibbs抽样、 MCMC 等现代统计方法 , 提出了基于统计方法的时空动态模型( Spatio-temporalDynamical Model) 来预报风暴的移动、 演变。3. 近年来出现(chxin)的新方法8共二十八页二、SCIT算

6、法(sun f)原理将某个连续三维区域的体积V 和雷达(lid)反射率Z 大于给定的体积阈值和反射率阈值( TZ) ( 其中也可包含一定数量的比反射率阈值小不到5 dBZ 的点)定义为风暴单体。反射率阈值的大小，根据不同的风暴类型来给定，可分成以下几类: 超级单体风暴： TZ = 40 50 dBz 对流风暴： TZ = 30 40 dBz 中尺度对流复合体： TZ = 25 30 dBz 雪带： TZ = 15 20 dBz9共二十八页1、SCIT算法子(f zi)功能风暴单体段020304识别反射率因子的径向排列，并输出这些段上的信息到风暴单体质心子功能中。将段组合成二维分量，并使这些分量

7、垂直相关构成三维单体，再计算这些单体的属性并输出至后两个子功能。依据风暴移动的历史来预报风暴将来的质心位置。风暴单体质心风暴单体跟踪风暴位置预报01通过将当前体积扫描发现的单体与前次体积扫描的单体作匹配来监视单体的移动。10共二十八页1.1 SCIT单体(dn t)识别方法-风暴单体段定义：一段连续的、沿着一个径向的一系列距离库，其反射率因子超过(chogu)一个规定的阈值。阈值（REFLECTIVITY） ：30,35,40,45,50,55, 60 dBZ (系列算法只使用一个阈值（30dBZ)-差别可调参数：DROPOUT REF DIFF:样本体积内，有效反射率与阈值的差值。默认5dB

8、Z （0-10dBZ）DROPOUT COUNT: 小于R的距离库数量。 默认2 （0-5）SEGMENT LENGTH 1-7 ：默认 1.9km （1.9-9.9km）大于该值的部分被保存。11共二十八页图2-1 沿径向(jn xin)的50,45,40,35及30dBZ阈值风暴段的识别图2-2 风暴(fngbo)段例子REFLECTIVITY：50dBZ、45dBZ、30dBZDROPOUT REF DIFF ：5DROPOUT COUNT ：212共二十八页1.2 SCIT单体识别方法-风暴(fngbo)单体分量当最后一个仰角扫描的段被处理完后，在空间上相邻的风暴段被结合成一个二维的风

9、暴单体分量(fn ling)。可调参数：AZIMUTHAL SEPARATION: 两个风暴段之间的偏离量 1.5 （1.5-3.5）SEGMENT OVERLAP : 两个风暴段的重叠部分 2km （0-5km）NUMBER OF SEGMENTS : 2D风暴段必须包含的1D的数目。 2 （1-4）COMPONENT AREA 1-7： 2D风暴段的面积。 10km2 (10-30km2)13共二十八页1.3 单体(dn t)识别方法-三维风暴单体的确定当当前体扫的所有仰角扫描都处理完之后，这些分量将按照质量的大小从大到小排列，然后做垂直相关分析。每一个确定的三维风暴单体由两个或更多个在相

10、继仰角上的二维分量构成。由于只用最强的（最高反射率因子阈值）二维分量来确定三维风暴单体，最终产品实际是一个三维风暴单体质心。这些风暴单体按它们(t men)的基于单体的垂直累积液态含水量（VIL）值排列。14共二十八页垂直(chuzh)相关分析 从第2个高度层开始每个高度层的风暴分量都和它下面相邻的高度层的风暴分量进行比较，可分成以下几个步骤:(1) 寻相邻仰角层的风暴分量的反射率质心之间的水平距离小于5 km的风暴分量为垂直相关的风暴分量。(2) 步骤(1)后，如果没有垂直相关的风暴分量，则把搜寻半径改为7.5km，继续搜寻。(3) 步骤(2)后，如果还没有垂直相关的风暴分量，则把搜寻半径改

11、为10km，继续搜寻。*将垂直相关的分量保存下来，不相关的当做2D风暴单体保存。1.3单体识别方法-三维风暴(fngbo)单体的确定15共二十八页两个很近3D单体结合成一个单体：（1）这两个单体在同一仰角上没有(mi yu)2D单体；（2）两个单体的质心的水平距离小于 HORIZONTAL MERGE 10km；（3）两单体间垂直空隙小于HEIGHT MERGE 4.0km EVEVATION MERGE 3.0两个很近3D单体结合成一个单体：（惠惠你看要不要加一个192上的东西？）1.3单体识别方法-三维风暴(fngbo)单体的确定16共二十八页图2-3 所有30dBZ的风暴段被结合(jih

12、)成一个二维分量，假定他们满足方位角距离和径向距离的判据如果对应(duyng)低阈值的分量的质心落于高阈值分量的范围内，则低阈值的分量被丢弃。图2-4 风暴单体段和它们构成的圆形二维分量以及由多个二维分量构成的三维风暴单体17共二十八页2、单体(dn t)追踪方法对在两个(lin )相继体扫中识别的风暴单体进行时间相关处理以确定每一个被识别的风暴单体的路径。检查体扫之间的时间间隔A利用前一个体扫单体质心位置确定目前体扫中该单体的初猜位置B计算当前体扫中识别的每个单体与其对应的初猜位置之间的距离C计算所有已经被识别的目前单体的新的运动向量D18共二十八页图2-5 时间相关联过程中的一个例子。代表

13、当前体扫单体的初猜位置，短箭头(jintu)代表当前体扫中单体位置和相应的初猜位置间的距离一旦完成追踪过程，将计算风暴单体的属性并将其直到前10个体扫的时间序列都存储(cn ch)起来。19共二十八页3、风暴位置(wi zhi)预报根据该风暴单体过去(guq)移动的记录来预报该风暴单体将来的之心位置。通过外推法来预报下一时刻的风暴单体的位置。位置预报最少做15min，最长做60min。20共二十八页三、SCIT算法(sun f)性能评估1）最大反射率=30dBZ；2）在至少一个仰角扫描范围内不小于30dBZ的回波强度水平尺度至少是52；3）估计区域的3D质心与另一个区域的质心至少远离5km；4

14、）估计区域要被一个局部极小反射率分开，该反射率比区域内的最大反射率至少低10dBZ。1）体扫模式追踪所有被识别单体的生命周期；2）确定不正确的时间关联的数目；3）注意它们的特点；4）计算用算法正确追踪的单体的百分率。结果显示：预报误差会随着时间而增加。单体识别单体追踪单体定位01020321共二十八页表3-1 对于各种风暴类型SCIT算法(sun f)的正确识别率表3-2 对于各种( zhn)风暴类型系列算法的正确识别率SCIT算法：40dBZ以上的单体被探测到的可能性有68%，50dBZ以上的单体被探测到的可能性有96%.而系列算法：40dBZ以上的单体被探测到的可能性有24%，50dBZ以

15、上的单体被探测到的可能性有41%.同时，算法结果对于孤立单体优于其它类型的降水系统。关于层状降水类型需要注意的一点是：对于这种降水系统算法性能并不太好。算法遗漏的两种类型是：1) 没有识别出大的成熟单体附近的衰亡或新生的单体。2） 没有识别在只有一个仰角扫描范围内符合尺度和强度标准的单体。22共二十八页表3-1 对于(duy)每一个例SCIT算法正确时间关联的百分比表3-2 SCIT算法对不同的主导时间的平均(pngjn)预报误差大多数时间关联误差发生于两个不同的紧密相连的单体的生长和消散期。预报误差会随着时间而增加。23共二十八页表3-1 a.体扫197，第四个仰角上，一个新的单体在单体1的

16、附近生成；b.体扫198，第四个仰角上，该新生单体被不正确的与单体1相关联；c.在添加(tin ji)了一个方向性的阈值（90）之后，图像上显示一个正确的、新识别的单体29.表3-1 a.组合反射率因子图（第167个体(gt)扫模式）其中一个新单体在17号单体预报路径上生成; b.下一个体(gt)扫（168）SCIT算法探测到该单体，但是错误的将它与17号单体时间关联。24共二十八页4.、产品(chnpn)及总结WSR-88D的9.0版本，由SCIT算法生成表格、字母和图表产品。这些产品提供关于风暴(fngbo)结构和追踪信息的风暴(fngbo)单体属性来帮助用户进行预测和预警过程。风暴(fn

17、gbo)追踪属性可用于每一种格式，而风暴(fngbo)结果属性仅仅用于字符格式。风暴(fngbo)结构趋势可以以图形方式显示在PUP上面。25共二十八页4、产品(chnpn)及总结相对于使用一个反射率阈值的风暴系列算法， SCIT算法采用七个反射率阈值，大大提高了风暴单体的识别；SCIT算法不能识别小的（长度或面积上），浅的（没有通过两个连续仰角扫描的阈值），或者最大反射率小于30dBZ（例如发展中的降水系统，高耸的Cu）的单体。另外，当降水风暴单体形成时极其接近(jijn)彼此，上述阈值可能将这些单体合并在一起，作为一个单体识别。这种情况往往发生在大面积相似的反射率飑线上；对孤立的风暴单体和组织程度高的飑线、多单体风暴中的单体可以以最大的精度追踪；短期风暴的位置预测往往比长期预测更准确。26共二十八页请大家(dji)批评与指正！27共二十八页内容