CN109917394B - 一种基于天气雷达的短临智能外推方法 - Google Patents
一种基于天气雷达的短临智能外推方法 Download PDFInfo
- Publication number
- CN109917394B CN109917394B CN201910187175.0A CN201910187175A CN109917394B CN 109917394 B CN109917394 B CN 109917394B CN 201910187175 A CN201910187175 A CN 201910187175A CN 109917394 B CN109917394 B CN 109917394B
- Authority
- CN
- China
- Prior art keywords
- storm
- monomer
- weather
- short
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于天气雷达的短临智能外推方法,所述方法包括对风暴单体历史个例的收集、识别和追踪,将获取的n+1项风暴单体轨迹序列作为范例存储于数据库中,当进行天气预测时,将获取的n项实时风暴单体序列与数据库中的所有范例进行相似度计算,找出相似度最高的范例,将该范例的第n+1项作为下个时次的风暴单体,以此计算出所有单体轨迹下个时次的状态,构成所预测的下个时次的天气状况,同时将新获取的实时风暴单体序列作为范例,存储于数据库中,以进行自学***。
Description
技术领域
本发明涉及气象学天气预测技术领域,特别涉及一种基于天气雷达的短临智能外推方法。
背景技术
对流天气的临近预报是指对未来0~2小时的对流天气***及其所伴随的灾害性天气的发生、发展、演变和消亡的预报,强对流天气是影响我国的主要灾害性天气之一,具有空间尺度小、生命史短、突发性强、发展演变迅速、破坏力大等特点,因此,针对强对流天气的临近预报一直是气象预报工作的重点和难点。
目前传统的天气预测方法有两种,即数值天气预报(NWP)法和基于光流矢量的天气预测方法。数值天气预报法是一种先对大气进行物理建模,再通过模拟物理模型进行预测的方法,该方法的优势体现在对于长时间范围内的天气预测比较准确,但是对0~2小时内的天气状况无法预测,因此不宜使用NWP方法进行短临预报。另一种在实际中得以应用的天气预测方法是基于光流矢量的方法,该方法是先通过两个雷达回波图像来估算光流矢量,该光流矢量可以理解为云的运动方向,然后根据最后一张雷达图进行合理的外推,得到最终的预测结果,该方法在0~2小时内的短临预测结果准确度更高。
以上基于雷达资料的短临外推预报方法的研究思路主要是围绕着单个天气变化过程中的对流单体演变过程,并未考虑不同天气变化过程中的不同对流单体间的类比,且缺乏自学习能力。
发明内容
发明目的:针对上述现有技术存在的问题和不足,本发明提供了一种基于天气雷达的短临智能外推方法,综合机器学***。
技术方案:本发明提供了一种基于天气雷达的短临智能外推方法,包括以下步骤:
步骤a:收集某个地区天气个例雷达基数据;
步骤b:提取每一时刻雷达基数据中的所有风暴单体,并提取、追踪其特征,获得若干条长度不等风暴单体轨迹序列;
步骤c:将获得的若干条长度不等风暴单体轨迹序列,切割成长度为n+1的风暴单体轨迹范例存储于数据库中;
步骤d:当进行天气预测时,读取当前n个时次的雷达基数据,利用SCIT识别、追踪,获得若干条长度为n的风暴单体轨迹序列,将该序列分别与数据库中所有历史风暴单体的前n项的轨迹范例进行相似度计算,找出相似度最高的范例,将该范例的第n+1项作为下个时次的风暴单体。按此方法计算出所有单体轨迹下个时次的状态,构成所预测的下个时次的天气状况,写为雷达基数据;
步骤e:按照步骤d进行迭代预测,直至推算出未来2小时全部的天气状况;
步骤f:将新获取的风暴单体序列切割成长度为n+1的范例,存储于数据库中,进行自学习,丰富范例库。
所述步骤b中提取风暴单体的方法为风暴单体识别追踪算法,提取的风暴单体特征包括当前季节、风暴体序列号、单体质心坐标、所在的仰角、质量权重体积、最大反射率、移动方向、移动速度。所述步骤d中,在进行相似度计算时,所筛选的历史风暴单体的形成位置位于地理相似判断范围内,所述地理相似判断范围为以待预测风暴单体的形成位置为圆心,半径不大于预设纬距的范围。
针对风暴单体的移动方向,设定与该移动方向偏差预设角度的范围,作为移向相似判断范围,针对风暴单体的移动速度,设定预设大小的上下波动范围,作为移速相似判断范围。
有益效果:与现有技术相比,本发明的优点在于:1、本发明提供了一种基于天气雷达的短临智能外推方法,不仅考虑天气***本身的演变信息,还将同一地区历史相似天气***的演变信息考虑在内,提高了强对流短临天气预报的预警能力和准确性。2、本发明所采用的预测方法只需要一台配置较高的小型服务器即可,相对于传统数值预报需要大型计算机进行辅助运行而言,大大缩短了预测所需要的时间,其次在预测过程中不需要进行大量的参数输入,大大提高了气象预测的工作效率。
附图说明
图1为基于范例推理的雷达基数据外推预报方法流程图;
图2为风暴单体识别流程图;
图3为质心法追踪流程图;
图4为二维风暴分量的识别示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。
如图1所示,基于机器学***稳季节ARIMA模型来进行拟合,更科学地实现不同序列片段间的相似性度量。各模块的方法如下:
(1)天气个例雷达基数据的收集、质量控制;
(2)基于雷达基数据的对流单体识别与跟踪;
(3)基于对流单体演变的时间序列进行特征选择;
(4)构建基于时序片段的对流单体演变范例库;
(5)设计适合对流单体演变序列的相似性度量算法;
(6)基于类比方法进行对流单体演变的短临外推;
(7)业务检验与对流单体演变范例库的自学习机制。
本发明设计的一种天气变化预测方法在实际应用过程当中,基于历史气象雷达基数据,针对风暴单体的移动路径进行预测,获得未来0~2小时预测的天气变化过程,具体包括如下步骤:
步骤a:收集某个地区天气个例雷达基数据。
步骤b:如图2所示,从一块地区的历史气象雷达基数据中识别风暴单体,具体步骤如下:
步骤b1:首先在体积扫描的每个径向上搜索强度大于一定阀值(默认为30dBz)的连续点,合并成有一定权重长度的一维风暴段;
步骤b2:如图4所示,在PPI(像素密度)层中按方位距离、长度重叠的相关性将风暴段合并成有一定权重面积(默认为10平方千米)的二维风暴分量;
步骤b3:最后按空间垂直相关将风暴分量合并成具有一定权重体积的三维风暴体,
识别完成后,获得风暴单体;
步骤b4:提取所识别风暴单体的特征,包括季节、质心位置、质量权重体积、最大反射率等,具体详见下表1所示:
表1风暴单体特征值
步骤c:如图3所示,采取质心追踪法追踪风暴单体的轨迹,所谓质心跟踪法,主要是根据前一时刻风暴的运动矢量和当前风暴与前一时刻风暴可能存在的位置关系对风暴单体的轨迹进行追踪,同时根据位置关系结果估计两个风暴单体的合并、***、新生和消亡。该算法的基本思想是:存在相似特征(包括体积、强度等)的风暴单体优先配对,近距离的风暴单体优先配对,但两者之间的最大速度不得超过限定值。风暴跟踪根据连续时间内的多个体扫识别出的风暴单体及其特征,通过对最后两个体扫匹配风暴体实现风暴跟踪。风暴单体的轨迹追踪的具体步骤如下:
步骤c1:用前一时刻t0每一个风暴体的运动矢量估计其在当前t1时刻可能的质量权重中心位置,如果该风暴体在t0时刻扫描体中是第一次测得,则用缺省的运动矢量,即:
其中Xexp、Yexp为缺省的运动矢量;vpast为上个时次的速度标量;θ为速度的方位角;Dt=t1-t0,表示两个风暴单体扫描间的时间间隔。
关键问题是如何给出t0时刻新生风暴体的初始速度,特别是当t0时刻是程序处理的第一个体扫描资料,无历史路径时,就必须较为准确的给出每一个风暴体的初始速度。通过分别计算两时刻风暴体的群体质量权重中心,来获取新生风暴体初始速度缺省值,方法如下:
t0时刻所有风暴体的群体质量权重中心为:
t1时刻所有风暴体的群体质量权重中心为:
则t0时刻任一风暴的速度缺省值为:
其中,Vz0i,Vz1i分别为t0和t1时刻任一风暴的反射率因子权重体积,Xz0i,Yz0i,Xz1i,Yz1i,分别为t0和t1时刻任一风暴的质量权重中心。
步骤c2:对于t0时刻探测的风暴,计算每个风暴在t1时刻的估计位置与t1时刻任一风暴体的质量权重中心之间的距离:
然后,按如下规则匹配两风暴:
(1)若t1时刻只有一个风暴与该风暴的距离小于等于相关风暴的最小距离阀值Td,且满足下列三个条件则认为这两个风暴是同一个风暴,实现匹配:
①两个风暴的反射率因子权重体积比大于或等于相关风暴最小体积比阈值;
②两个风暴的移动方向之间的夹角小于相关风暴最大移向偏差阈值;
③两个风暴的移动速度小于、等于相关风暴最大移速阈值。
其中(Xpast,Ypast)、(Xcur,Ycur)分别表示t0、t1时刻风暴的反射率权重中心。
(2)若t1时刻有多个风暴与该风暴的距离小于等于相关风暴的最小距离阀值Td,按以下情况处理:
如果t1时刻这几个风暴的质量权重体积之和与t0时刻风暴的质量权重体积的比值小于等于风暴***最小体积比阀值Tsplit。
其中M、Vcuri分别为t1时刻满足距离阀值的风暴数及任一风暴的质量权重体积。
说明可能存在风暴***,只要t1时刻这几个风暴满足(1)中的②和③两个条件,就认为它们都与t0时刻的风暴相关。
反之,如果t1时刻这几个风暴的质量权重体积之和与t0时刻风暴的质量权重体积的比值大于风暴***最小体积比阀值Tsplit,即:
则认为不可能存在风暴***,取t0时刻该风暴距离最小的一个与该风暴匹配,并判断是否满足(1)中的②和③两个条件,如果满足就认为其与t0时刻该风暴相关。
(3)若t1时刻没有风暴与该风暴的距离小于等于相关风暴的最小距离阀值Td,则认为t0时刻的该风暴在t1时刻消亡或发生合并。
步骤c3:若t0时刻的某一风暴在t1时刻没有风暴与之相关,就认为发生风暴的合并,若t1时刻的某一风暴在t0时刻没有风暴与之相关,则标记该风暴为新生风暴单体。
步骤c4:根据匹配结果,计算t1时刻每一风暴的运动矢量:
利用t1时刻风暴的运动矢量预报未来t2时刻的质量权重中心位置,如果该风暴是新生风暴,则不做预报,计算完成后,获得若干风暴单体的运动轨迹。
步骤d:将获得的若干条长度不等风暴单体轨迹序列,切割成长度为n+1的风暴单体轨迹范例存储于数据库中,自此范例库搭建完毕。
步骤e:此部分为天气预测部分,获取当前时段天气雷达基数据,执行步骤a至步c,获得若干风暴单体的运动轨迹,截取轨迹的前n项,得到若干个长度为n的风暴单体序列。依次将这些序列与数据库中范例的前n项进行对比,根据表2关联表的算法计算相似度(Hss),选择相似度最高的范例,将该范例的第n+1项风暴单体作为该风暴序列下个时次(未来第6分钟)的风暴单体状态,计算完成后,获得下个时次(未来第6分钟)的天气状况。
表2Contingency Table
其中,hits表示命中,即事件预测将发生,且确实发生;misses表示未命中,即事件预测不会发生,但确实发生;false alarms表示误报,即事件预测将发生,但没有发生;correct negatives表示正确的负面,即事件预测不会发生,也不会发生;Yes表示像素点有反射率,No表示无反射率。将两个风暴单体同一位置的像素点进行比较,当都显示为Yes时,hits+1,当都显示为No时,correct negatives+1。例如用Hss算法分别计算出前5对风暴单体的相似度,分别为S1、S2、S3、S4、S5,每个位置风暴单体的权重为W1~W5,距离当前时次越近,权重越高,默认权重为0.2、0.4、0.6、0.8、1.0,最后两条轨迹的总相似度为通过比较整个数据库中的范例找出与目标轨迹相似度最高的范例,以该范例的最后位的风暴单体作为预报结果。
步骤f:将预测的6分钟的天气状况作为已知数据,转步骤e,迭代预测出下下个时次(未来第12分钟)所有风暴单体的状态,重复此过程,直至预测未来2小时的天气状况,至此0~2小时外推结束。
Claims (4)
1.一种基于天气雷达的短临智能外推方法,包括以下步骤:
步骤a:收集天气个例雷达基数据;
步骤b:提取每一时刻雷达基数据中的所有风暴单体,并提取、追踪其特征,获得若干条长度不等风暴单体轨迹序列;具体包括:
步骤b1:首先在体积扫描的每个径向上搜索强度大于一定阈值的连续点,合并成有一定权重长度的一维风暴段;
步骤b2:在像素密度层中按方位距离、长度重叠的相关性将风暴段合并成有一定权重面积的二维风暴分量;
步骤b3:最后按空间垂直相关将风暴分量合并成具有一定权重体积的三维风暴体,识别完成后,获得风暴单体;
步骤b4:提取所识别风暴单体的特征,包括当前季节、风暴体序列号、单体质心坐标、所在的仰角、质量权重体积、最大反射率、移动方向、移动速度;
步骤c:将获得的若干条长度不等风暴单体轨迹序列,切割成长度为n+1的风暴单体轨迹范例存储于数据库中;
步骤d:在进行天气预测时,读取当前n个时次的雷达基数据,利用SCIT识别、追踪,获得若干条长度为n的风暴单体轨迹序列,将该序列分别与数据库中所有历史风暴单体的前n项的轨迹范例进行相似度计算,找出相似度最高的范例,将该范例的第n+1项作为下个时次的风暴单体,按此方法计算出所有单体轨迹下个时次的状态,构成所预测的下个时次的天气状况,写为雷达基数据;
步骤e:按照步骤d进行迭代预测,直至推算出未来2小时全部的天气状况;
步骤f:将新获取的风暴单体序列切割成长度为n+1的范例,存储于数据库中,进行自学习,丰富范例库。
2.根据权利要求1所述的一种基于天气雷达的短临智能外推方法,其特征在于,所述步骤b中提取风暴单体的方法为风暴单体识别追踪算法。
3.根据权利要求1所述的一种基于天气雷达的短临智能外推方法,其特征在于,所述步骤d中,在进行相似度计算时,所筛选的历史风暴单体的形成位置位于地理相似判断范围内。
4.根据权利要求3所述的一种基于天气雷达的短临智能外推方法,其特征在于,所述地理相似判断范围为以待预测风暴单体的形成位置为圆心,半径不大于预设纬距的范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910187175.0A CN109917394B (zh) | 2019-03-13 | 2019-03-13 | 一种基于天气雷达的短临智能外推方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910187175.0A CN109917394B (zh) | 2019-03-13 | 2019-03-13 | 一种基于天气雷达的短临智能外推方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109917394A CN109917394A (zh) | 2019-06-21 |
CN109917394B true CN109917394B (zh) | 2022-12-23 |
Family
ID=66964455
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910187175.0A Active CN109917394B (zh) | 2019-03-13 | 2019-03-13 | 一种基于天气雷达的短临智能外推方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109917394B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111366989A (zh) * | 2020-03-23 | 2020-07-03 | 上海眼控科技股份有限公司 | 天气预报方法、装置、计算机设备及存储介质 |
CN111815542B (zh) * | 2020-07-31 | 2024-05-07 | 南京信息工程大学 | 一种树木年轮图像髓心定位、年轮测量方法 |
CN112232674B (zh) * | 2020-10-16 | 2021-12-07 | 中国气象局气象探测中心 | 一种气象灾害评估方法、装置和*** |
CN113642813B (zh) * | 2021-10-18 | 2022-02-11 | 江苏铨铨信息科技有限公司 | 一种基于物理方程的降水外推预报方法 |
CN114740549B (zh) * | 2022-04-12 | 2024-02-06 | 北京墨迹风云科技股份有限公司 | 一种短临降水天气预测方法、装置、电子设备及可读介质 |
CN116933665A (zh) * | 2023-09-18 | 2023-10-24 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种基于深度学习的超级计算机预报模型的应用方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105488594A (zh) * | 2015-12-08 | 2016-04-13 | 南京信息工程大学 | 一种台风路径预测方法 |
CN107229084A (zh) * | 2017-06-08 | 2017-10-03 | 天津大学 | 一种自动识别,跟踪并预测对流***目标的方法 |
CN108535731A (zh) * | 2018-04-18 | 2018-09-14 | 青岛心中有数科技有限公司 | 短临降水预报方法与装置 |
CN108802733A (zh) * | 2018-06-15 | 2018-11-13 | 安徽诺安信息科技有限公司 | 一种基于气象雷达数据对冰雹探测的算法 |
CN108828691A (zh) * | 2018-04-18 | 2018-11-16 | 青岛心中有数科技有限公司 | 短临降水预报方法与装置 |
CN109100722A (zh) * | 2018-07-25 | 2018-12-28 | 南京信息工程大学 | 基于雷达回波图像扇区分量分析的风暴趋势预测方法 |
CN109164450A (zh) * | 2018-09-12 | 2019-01-08 | 天津大学 | 一种基于多普勒雷达数据的下击暴流预测方法 |
-
2019
- 2019-03-13 CN CN201910187175.0A patent/CN109917394B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105488594A (zh) * | 2015-12-08 | 2016-04-13 | 南京信息工程大学 | 一种台风路径预测方法 |
CN107229084A (zh) * | 2017-06-08 | 2017-10-03 | 天津大学 | 一种自动识别,跟踪并预测对流***目标的方法 |
CN108535731A (zh) * | 2018-04-18 | 2018-09-14 | 青岛心中有数科技有限公司 | 短临降水预报方法与装置 |
CN108828691A (zh) * | 2018-04-18 | 2018-11-16 | 青岛心中有数科技有限公司 | 短临降水预报方法与装置 |
CN108802733A (zh) * | 2018-06-15 | 2018-11-13 | 安徽诺安信息科技有限公司 | 一种基于气象雷达数据对冰雹探测的算法 |
CN109100722A (zh) * | 2018-07-25 | 2018-12-28 | 南京信息工程大学 | 基于雷达回波图像扇区分量分析的风暴趋势预测方法 |
CN109164450A (zh) * | 2018-09-12 | 2019-01-08 | 天津大学 | 一种基于多普勒雷达数据的下击暴流预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109917394A (zh) | 2019-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109917394B (zh) | 一种基于天气雷达的短临智能外推方法 | |
KR102129893B1 (ko) | 딥러닝 네트워크 및 평균 이동을 기반으로 하는 선박 자동추적 방법 및 시스템 | |
CN109738970B (zh) | 基于雷电数据挖掘实现雷电预警的方法、装置和存储介质 | |
CN110033051B (zh) | 一种基于多步聚类的拖网渔船行为判别方法 | |
CN110889324A (zh) | 一种基于yolo v3面向末端制导的热红外图像目标识别方法 | |
Jing et al. | AENN: A generative adversarial neural network for weather radar echo extrapolation | |
CN104156984A (zh) | 一种不均匀杂波环境下多目标跟踪的概率假设密度方法 | |
Li et al. | Deep spatio-temporal wind power forecasting | |
CN113283653B (zh) | 一种基于机器学习和ais数据的船舶轨迹预测方法 | |
Thongniran et al. | Spatio-temporal deep learning for ocean current prediction based on HF radar data | |
Zhuang et al. | Long-lead prediction of extreme precipitation cluster via a spatiotemporal convolutional neural network | |
CN115546710A (zh) | 水电厂的人员定位方法、装置、设备及可读存储介质 | |
Sun et al. | Image target detection algorithm compression and pruning based on neural network | |
Wu et al. | Inductive spatiotemporal graph convolutional networks for short-term quantitative precipitation forecasting | |
Viswambari et al. | Data mining techniques to predict weather: a survey | |
CN115048986B (zh) | 一种基于多分类器动态剪枝选择的地表冻融状态分类方法 | |
US10801841B1 (en) | Trajectory prediction via a feature vector approach | |
CN110515081A (zh) | 一种雷达回波零度层亮带智能识别预警方法 | |
Pillai et al. | Fine-Tuned EfficientNetB4 Transfer Learning Model for Weather Classification | |
Krishna | A review of weather forecasting models-based on data mining and artificial neural networks | |
Socaci et al. | XNow: A deep learning technique for nowcasting based on radar products’ values prediction | |
CN114137637A (zh) | 基于闪电和雷达数据的雷暴中心踪迹的集合概率预报方法 | |
Maceda et al. | Rainfall Classification Model for the Philippines using Optimized K-nearest Neighbor Algorithm with GridSearchCV Hyperparameter Tuning | |
Chawla et al. | Data efficient Random Forest model for avalanche forecasting | |
Wang et al. | Tiny-RainNet: A Deep CNN-BiLSTM Model for Short-Term Rainfall Prediction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |