CN115691049A - 一种基于深度学习的对流初生预警方法 - Google Patents

一种基于深度学习的对流初生预警方法 Download PDF

Info

Publication number
CN115691049A
CN115691049A CN202211076007.2A CN202211076007A CN115691049A CN 115691049 A CN115691049 A CN 115691049A CN 202211076007 A CN202211076007 A CN 202211076007A CN 115691049 A CN115691049 A CN 115691049A
Authority
CN
China
Prior art keywords
convection
remote sensing
time
image
resolution
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.)
Pending
Application number
CN202211076007.2A
Other languages
English (en)
Inventor
杨春蕾
谢梦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Institute Of Technical Physics
Yunyao Power Technology Suzhou Co ltd
Original Assignee
Suzhou Institute Of Technical Physics
Yunyao Power Technology Suzhou Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Suzhou Institute Of Technical Physics, Yunyao Power Technology Suzhou Co ltd filed Critical Suzhou Institute Of Technical Physics
Priority to CN202211076007.2A priority Critical patent/CN115691049A/zh
Publication of CN115691049A publication Critical patent/CN115691049A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Image Processing (AREA)

Abstract

本发明涉及一种基于深度学习的对流初生预警方法,包括:建立深度超分辨率模型以提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息,并对遥感影像进行训练以重构所述遥感影像的图像分辨率和时间分辨率;获取重构的遥感影像按照预定的对流云判识策略采用云顶亮温阈值法进行对流云的初判识;基于面积重叠法利用计算机视觉算法提取对流云以进行对流云追踪;基于GOES‑R算法对最终的对流云遥感影像进行判识、依据预设的对流初生判定条件以判识对流初生;以及,结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估并结合雷达数据对对流初生反演结果进行检验。本发明能够实现强对流天气监测预警,提高对流初生的判识精度。

Description

一种基于深度学习的对流初生预警方法
技术领域
本发明涉及强对流天气预警技术领域,特别是涉及一种基于深度学习的对流初生预警方法。
背景技术
强对流天气预警是灾害防御的第一道防线,也是最重要的一环。但是,目前各国对于强对流天气信息的掌握都还不够全面。在实际业务中通常应用雷达反射率(基于雷达数据的线性外推方法)进行强对流的辨别,但是高原、沙漠、海上等地没有雷达覆盖,所以很难做出时间地点精准的预报。
由于天气***有着发展复杂且快速变化的特点,准确判断并完成追踪等过程变成了薄弱环节。随着对强对流天气生成条件的不断了解,对流初生是强对流发生的重要发展阶段,对流初生是指气象雷达综合反射率图像上首次出现大于35dBZ的回波。研究发现通过监测对流初生过程能够准确预警出强对流天气的出现,对能够提早预报预警起到重要作用。
目前,国际上主要的对流初生算法有4种,分别为RDT(Rapidly DevelopingThunderstorms),ForTraCC(Forecasting and Tracking the Evolution of CloudClusters),GOES-R (The Geostationary Operational Environment Satellite R-Series Program)和UWCI(The University of Wisconsin Convective Initiation)。在对流判识过程中,RDT算法更加关注垂直形态的峰值检测,并使用速度外推追踪对流,其业务算法结合多种资料如数值预报产品和卫星反演产品;ForTraCC算法对活跃对流更加关注,并综合考虑了对流合并与***的情况,能实现外推预报;GOES-R算法使用多光谱判识技术监测对流云,主要针对静止气象卫星;UWCI 算法提出了块平均(box average)概念,并使用云类型产品和单通道影像作为辅助数据,技术简单、快速。
目前,针对强对流预警的研究中,中国气象卫星中心开发的对流初生产品由于算法不公开,产品发布具有滞后性(通常为15-20分钟),无法在实际业务中服务,且还需要完成阈值优选、检验验证等大量工作。而中央气象台和清华大学联合开发了一种基于深度神经网络的雷达回波外推方法,预报准确率提高约40%,但其因雷达分布不均,无法满足海上等无雷达地区的强对流预警需求。同时,上述传统算法都只使用了近红外通道,由于近红外通道影像分辨率较低,导致难以实现精细化预警。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种基于深度学习的对流初生预警方法,用于解决现有技术中对流初生预警滞后且精度低的问题。
为实现上述目的及其他相关目的,本发明提供如下技术方案:
一种基于深度学习的对流初生预警方法,包括:
建立深度超分辨率模型以提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息,并对遥感影像进行训练以重构所述遥感影像的图像分辨率和时间分辨率;
获取重构的遥感影像按照预定的对流云判识策略采用云顶亮温阈值法进行对流云的初判识;
基于面积重叠法利用计算机视觉算法提取对流云以进行对流云追踪;
基于GOES-R算法对最终的对流云进行判识、依据预设的对流初生判定条件以判识对流初生;以及,
结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估,并结合雷达数据对对流初生反演结果进行检验。
作为发明的一种优选方案,所述建立深度超分辨率模型以提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息具体包括:
利用深度学习的非线性映射能力和对栅格数据的信息提取能力建立深度超分辨率模型;
利用深度超分辨率模型提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息;以及,
利用获取的有效信息将时空低分辨率遥感影像重构为时空高分辨率遥感影像。
作为发明的一种优选方案,所述对遥感影像进行训练以重构所述遥感影像的图像分辨率和时间分辨率具体包括:
采用ESRGAN模型对不同图像分辨率的遥感影像进行训练以将所述遥感影像重构为高图像分辨率的遥感影像,所述高图像分辨率的遥感影像保留原遥感影像的图像分辨率的通道信息;以及,
采用Super SloMo模型对不同时间分辨率的遥感影像进行训练以将所述遥感影像重构为高时间分辨率的遥感影像,所述高时间分辨率的遥感影像与雷达数据的时间分辨率保持一致。
作为发明的一种优选方案,所述预定的对流云判识策略具体为:对流云体的最冷温度阀值达到-55℃、对流云体的最暖温度阀值达到-10-5℃、对流云体在垂直方向上的温度等值线间隔达到1℃、对流云体极值达到3℃、且对流判识的最小面积阀值达到一个红外像元大小。
作为发明的一种优选方案,所述采用云顶亮温阈值法进行对流云的初判识时利用RDT算法,并采用自适应阈值选取长波红外通道阈值。
作为发明的一种优选方案,所述预设的对流初生判定条件具体为:10.7μm亮温≤0℃、 10.7μm亮温时间变化率≤[-4℃(15min)-1]、-35℃<6.5和10.7μm亮温差<-10℃、-25℃<13.3 和10.7μm亮温差<-5℃、6.5和10.7μm亮温时间变化率>[3℃(15min)-1]、且13.3和10.7μm 亮温时间变化率>[3℃(15min)-1]。
作为发明的一种优选方案,所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估具体包括:
采用峰值信噪比和结够相似性参数评价时空降尺度的准确性;其中,
所述峰值信噪比
Figure BDA0003829620310000031
其中,
Figure BDA0003829620310000032
表示图像点颜色的最大数值;
所述结构相似性
Figure BDA0003829620310000033
其中μx、μy为均值,σx、σy为标准差,σxy为协方差,c1=(0.01L)2,c2=(0.03L)2,L对于8bit灰度图取值为255。
作为发明的一种优选方案,所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估还包括:
采用均方根误差和空间相关系数判断预测结果的精度;其中,
所述均方根误差
Figure BDA0003829620310000034
所述空间相关系数
Figure BDA0003829620310000035
其中,n是像素数,Oi和Si为观测值和预测值,
Figure BDA0003829620310000036
Figure BDA0003829620310000037
分别为相应的平均值。
作为发明的一种优选方案,所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估之后,还包括:
根据检验评估结果优化对流初生判识的参数和阈值。
如上所述,本发明的基于深度学习的对流初生预警方法,具有以下有益效果:
本发明采用深度学习技术对风云卫星影像进行时空降尺度,空间降尺度是在空间范围内将影像的低分辨率降尺度为高分辨率,提高地物信息的分辨能力,反映出更细节的地面信息,时间降尺度是将粗时间分辨率的影像降尺度成更细时间分辨率的影像。由于风云卫星主要用于天气监测与预报、防灾减灾、气候变化、生态环境监测、服务“一带一路”建设等领域,因此,通过深度学习提高风云卫星的空间分辨率和时间分辨率,可为后续的降水预报、天气监测等提供更加精细化的数据;并且解决了多源数据时空分辨率不匹配问题,对气象的模式预报和提供时空分辨率统一的高精度降水服务具有深远的意义;而且时空降尺度得到更高时空分辨率的遥感影像,可以提升中小尺度云团连续快速监测能力,对强对流天气监测预警,尤其对突发性中小尺度天气***监测具有重要意义。
本发明采用时空降尺度方式进行资料融合,实现云图细节信息有效利用,解决小尺度漏报和大尺度延迟的问题,大大提高了预报准确性。其中,确定空间尺度从4km降至1km,时间尺度从15min降至6min。同时,在实际的对流初生反演过程中,涉及到很多气象参数和阈值的设定,如对流云的识别,通过云顶温度的经验值(~241K)判别,本发明通过深度学习的方法,结合RDT算法和GOES-R算法,不断优化算法在目标判识、目标追踪、对流初生判识的各个过程中判别条件权重和阈值的优化,对不同地区有针对性的改变阈值,提高对流初生的精度;同时,本发明结合如面积重叠法、交叉相关法、计算机视觉等多种云团追踪技术,其在处理线性运动目标跟踪方面有独到的优势。
附图说明
图1显示为本发明实施例一的方法流程图。
图2显示为本发明实施例一中ESRGAN模型结构示意图。
图3显示为本发明实施例一中Super SloMo模型结构示意图。
图4显示为本发明实施例一中应用RDT算法得到的云团识别图。
图5a显示为本发明实施例二中4km×4km的低分辨率影像。
图5b显示为本发明实施例二中双三次插值得到的高分辨率影像。
图5c显示为本发明实施例二中基于ESRGAN模型使用预训练权重到的高分辨率影像。
图5d显示为本发明实施例二中基于ESRGAN模型迁移学习到的高分辨率影像。
图5e显示为本发明实施例二中1km×1km的实际高分辨率影像。
图6a显示为本发明实施例三中基于Super SloMo模型预测的影像。
图6b显示为本发明实施例三中真实遥感影像。
图7显示为本发明实施例四中对流初生预警的影像。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。
在本发明的说明书和权利要求书及上述附图中的描述的一些流程中,包含了按照特定顺序出现的多个操作,但是应该清楚了解,这些操作可以不按照其在本文中出现的顺序来执行或并行执行,操作的序号如101、102等,仅仅是用于区分开各个不同的操作,序号本身不代表任何的执行顺序。另外,这些流程可以包括更多或更少的操作,并且这些操作可以按顺序执行或并行执行。需要说明的是,本文中的“第一”、“第二”等描述,是用于区分不同的消息、设备、模块等,不代表先后顺序,也不限定“第一”和“第二”是不同的类型。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
请参阅图1,本发明提供一种基于深度学习的对流初生预警方法,包括:
S100、建立深度超分辨率模型以提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息,并对遥感影像进行训练以重构所述遥感影像的图像分辨率和时间分辨率。
具体的,建立深度超分辨率模型以提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息具体包括:
S101、利用深度学习的非线性映射能力和对栅格数据的信息提取能力建立深度超分辨率模型;
S102、利用深度超分辨率模型提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息;以及,
S103、利用获取的有效信息将时空低分辨率遥感影像重构为时空高分辨率遥感影像。
对遥感影像进行训练以重构所述遥感影像的图像分辨率和时间分辨率具体包括:
S104、采用ESRGAN模型对不同图像分辨率的遥感影像进行训练以将所述遥感影像重构为高图像分辨率的遥感影像,所述高图像分辨率的遥感影像保留原遥感影像的图像分辨率的通道信息,本实施例中,通过采用ESRGAN模型将遥感影像的图像分辨率从重构到1km,并且能保留4km图像分辨率的通道信息。
如图2所示,图2示出了ESRGAN模型的具体网络模型结构示意图。本实施例中采用ESRGAN的生成器实现图像超分辨率重建,ESRGAN(Enhanced Super-ResolutionGenerative Adversarial Networks)模型是超分辨率生成对抗网络,能够在单幅图像超分辨率期间生成真实感纹理。其基于SRGAN改进,主要体现在以下三个方面:(1)使用对抗损失和感知损失; (2)引入Residual-in-Residu Dense Block(RRDB),结合多层残差网络和密集连接,该结构更易训练;(3)使用激活前的VGG特征来优化感知损失,可以提供对亮度一致性和纹理恢复更强的监督力。
S105、采用Super SloMo模型对不同时间分辨率的遥感影像进行训练以将所述遥感影像重构为高时间分辨率的遥感影像,所述高时间分辨率的遥感影像与雷达数据的时间分辨率保持一致,本实施例中,通过采用Super SloMo模型将遥感影像的时间分辨率从15min重构到 6min。
如图3所示,图3示出了Super SloMo模型结构示意图。Super SloMo模型是一种端到端的卷积神经网络,能在两个输入图像之间生成任意多的中间视频帧,其整个框架是由两个全卷积神经网络U-Net组成。首先,用一个U-Net计算相邻输入图像之间的双向光流;再在每个时间步长上来线性拟合光流,进而逼近中间帧的双向光流,为解决运动边界出现伪影的问题,使用另一个U-Net来优化近似的光流;最后,将输入的两张图像进行扭曲和线性融合,从而形成中间帧。此外,Super SloMo模型的光流计算网络和插值网络的参数,都不依赖于***值帧的特定时间步长(时间步长被为网络的输入),因此,可以并行地在两帧之间的任意时间步长上插帧,从而突破了很多单帧插值方法的局限性。
S200、获取重构的遥感影像按照预定的对流云判识策略采用云顶亮温阈值法进行对流云的初判识。
具体的,采用云顶亮温阈值法进行对流云的初判识时利用RDT算法,考虑对流云体在垂直方向上的高度变化,采用自适应阈值选取长波红外通道阈值,RDT算法由欧洲气象卫星开发组织EUMETSATSAF开发;预定的对流云判识策略具体为:对流云体的最冷温度阀值达到-55℃、对流云体的最暖温度阀值达到-10-5℃、对流云体在垂直方向上的温度等值线间隔达到1℃、对流云体极值达到3℃、且对流判识的最小面积阀值达到一个红外像元大小,同时满足上述条件即可判定位对流云。如图4所示,图4示出了应用RDT算法得到的云团识别图。
S300、基于面积重叠法利用计算机视觉算法提取对流云以进行对流云追踪,其中,计算机视觉算法具体采用Boosting算法,通过选取重叠区域值,考虑相关值,利用boosting算法可从一系列的图像中自动提取,进行分析和理解,实现对对流云的追踪。
S400、基于GOES-R算法对最终的对流云遥感影像进行判识、依据预设的对流初生判定条件以判识对流初生。
本实施例中具体参考当前GOES卫星业务CI算法采用的判据,并使用GOES-R ABI(Advanced Baseline Imager)仪器的3个通道数据,对流初生判定条件具体图下表所示:
Figure BDA0003829620310000071
S500、结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估,并结合雷达数据对对流初生反演结果进行检验。
具体的,所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估具体包括:
S501、采用峰值信噪比(Peak Signal to Noise Ratio,PSNR)和结构相似性(Structure Similarity Index Measure,SSIM)参数评价时空降尺度的准确性;其中,
所述峰值信噪比
Figure BDA0003829620310000072
其中,
Figure BDA0003829620310000073
表示图像点颜色的最大数值,若每个采样点用8位表示,则为255;
所述结构相似性
Figure BDA0003829620310000074
其中μx、μy为均值,σx、σy为标准差,σxy为协方差,c1=(0.01L)2,c2=(0.03L)2,L对于8bit灰度图取值为255。
通常情况下,PSNR的数值越高,表示重建图像与原始图像越接近。为更加准确评价图像质量,再使用结构相似性(SSIM)评价图像质量,结构相似性是超分辨图像与实际高分辨率图像的结构相似程度,取值不大于1,其值越接近1则表明两者结构越相似,影像重建效果越好。
S502、采用均方根误差(Root Mean Square Error,RMSE)和空间相关系数(Correlation Coefficient,CC)判断预测结果的精度;其中,
所述均方根误差
Figure BDA0003829620310000081
所述空间相关系数
Figure BDA0003829620310000082
其中,n是像素数,Oi和Si为观测值和预测值,
Figure BDA0003829620310000083
Figure BDA0003829620310000084
分别为相应的平均值。
RMSE是衡量预测值与实测值之间的偏差,RMSE值越小,偏差越小,代表模型精度越高;CC值越接近1则说明它们相关性越高。
所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估之后,还包括:
S503、根据检验评估结果优化对流初生判识的参数和阈值。
本发明使用四种评分标准对对流初生反演结果的预报能力进行评估,分别为公正预报率 (equitable threat score,ETS)、频率偏差Frequency Bias、命中率(probabilityof detection, POD)和假警率(false alarm ratio,FAR)公正预报率ETS是指评估去掉随机因素的发生的事件被正确预报出的比率,频率偏差Frequency Bias是指评估预报发生的事件和实际发生的事件的比率(大于1表示高估了事件发生的频率,小于1代表低估),命中率POD是指评估发生的事件被正确预报出的比率,假警率FAR是指评估实际没发生但预报有的事件的比率,它们的最优值分别是1、1、1、0,变化范围分别为-1/3~1、0~1、0~1、0~∞,其计算公式具体如下:
Figure BDA0003829620310000085
Figure BDA0003829620310000086
Figure BDA0003829620310000087
Figure BDA0003829620310000088
其中,H表示正确预报的格点数,M表示漏报格点数,F表示空报格点数,C表示正确预报出未存在对流初生事件格点数,具体如下表所示。
Figure BDA0003829620310000091
本发明在传统对流初生反演算法的基础上,基于国产风云卫星进行算法研发及改进,为满足灾害性预警的精细化需求,重点处理了风云卫星数据部分,采用深度学习技术对风云卫星影像进行时空降尺度,空间降尺度是在空间范围内将影像的低分辨率降尺度为高分辨率,提高地物信息的分辨能力,反映出更细节的地面信息,时间降尺度是将粗时间分辨率的影像降尺度成更细时间分辨率的影像。由于风云卫星主要用于天气监测与预报、防灾减灾、气候变化、生态环境监测等领域,因此,通过深度学习提高风云卫星的空间分辨率和时间分辨率,可为后续的降水预报、天气监测等提供更加精细化的数据;并且解决了多源数据时空分辨率不匹配问题,对气象的模式预报和提供时空分辨率统一的高精度降水服务具有深远的意义;而且时空降尺度得到更高时空分辨率的遥感影像,可以提升中小尺度云团连续快速监测能力,对强对流天气监测预警,尤其对突发性中小尺度天气***监测具有重要意义。
本发明采用时空降尺度方式进行资料融合,实现云图细节信息有效利用,解决小尺度漏报和大尺度延迟的问题,大大提高了预报准确性。其中,确定空间尺度从4km降至1km,时间尺度从15min降至6min。同时,在实际的对流初生反演过程中,涉及到很多气象参数和阈值的设定,如对流云的识别,通过云顶温度的经验值(~241K)判别,本发明通过深度学习的方法,结合RDT算法,不断优化算法在目标判识、目标追踪、对流初生判识的各个过程中判别条件权重和阈值的优化,对不同地区有针对性的改变阈值,提高对流初生的精度;同时,本发明结合如面积重叠法、交叉相关法、计算机视觉等多种云团追踪技术,其在处理线性运动目标跟踪方面有独到的优势。
本发明将实现高时空精度的对流初生预警***,同时本发明使用的数据为国产的风云卫星,将提高在航空、气象、农业等方面的预警能力。
1、航空方向:为航空公司提供对流初生算法,实现可以根据实时的卫星数据生成对流的范围、强度信息,为飞机的航行提供气象保障。具体对标美国气象公司(The WeatherCompany) 中的航空产品实时湍流报告(WSI Total Turbulence),在航路上集成对流初生预警产品。
据统计,我国每年因航班延误造成的经济损失就高达500亿元人民币,只要能正确掌握和运用气象条件能够避免2/3的损失。因此,准确监测和预报强对流天气的发生及利用强对流气象预报产品指导航空业务运行,可减小因天气问题造成的经济和生命损失。
通过提供强对流信息着重为中国航空公司飞行、尤其国际飞行提供气象保障,实现便捷、高效的气象服务。并通过提供有助于客户控制成本、增强飞行安全和提高效率的关键天气信息,结合中国的卫星、雷达、地面站点数据,为全国乃至全球航空飞行提供气象服务。
2、气象方向:各部门联动,加强服务平台、灾害防御体系建设,并妥善做好减灾工作是气象部门的主要任务之一。而基于深度学习的灾害性天气预警的目的是从源头上解决问题,比地面雷达更早地预警降雨***的发展,提早发出强对流天气的预警,从技术上提升预警能力。其在气象业务方面的应用效果突出,产生了良好的效益。
3、其他方面:在农业、物流、航运等方面提前预警,有助于物流、航运公司及时掌握灾害天气预警等级,提高灾害天气应对能力,并高效采取措施,将损失和风险降到最低。此外,成功的预警强对流,可以有效减少因此类天气造成的工农业生产的损失,保障人民的生命财产安全。
实施例二
本实施例公开了一空间降尺度实例,具体为:选用2018年7月16日4张影像数据作为测试数据(尺寸为1200×784),使用不同的方法对同一个测试集进行测试,并计算PSNR、SSIM、RMSE和CC指标。源资料风云4A影像空间分辨率为4km×4km,目标资料为1km ×1km分辨率的双三次插值、基于ESRGAN模型预训练权重方法和迁移方法模型的精度指标评价,具体精度指标评价对比如下表所示:
Figure BDA0003829620310000101
根据对比结果可知,基于ESRGAN模型分别使用预训练权重方法和迁移方法,它们在 PSNR、SSIM和EMSE指标上都比双三次插值方法取得更好的空间降尺度效果,表明基于深度学习的方法具备提取历史空间相关性的能力。基于ESRGAN模型的迁移方法较之使用预训练权重方法在PSNR、SSIM、RMSE和CC指标上均提高,并分别提高了0.827、0.009、0.002、0.001,说明迁移方法有效地学习了风云4A号遥感影像的空间特征。
图5a~图5e为选取的2018年7月16日8点19分影像空间降尺度的效果,其中图5a为4km×4km的低分辨率影像,图5b为双三次插值得到的高分辨率影像,图5c为基于ESRGAN 模型使用预训练权重到的高分辨率影像,图5d为基于ESRGAN模型迁移学习到的高分辨率影像,图5e为1km×1km的实际高分辨率影像。由图5a~图5e可明显看出相对比双三次插值方法,在4倍分辨率提升场景下基于深度学习的方法也能够更准确地重构出较高分辨率的遥感影像,并且基于ESRGAN模型使用迁移方法得到的高分辨率影像边缘清晰,总体效果上与真实高分辨率影像更接近。
实施例三
本实施例公开了一时间降尺度实例,具体为:根据风云4A影像2018年7月16日8点15分和8点30分的影像,设置时间步长,输入模型,测试出在该时间段任意的影像。图6a 和图6b是深度学习预测出的8点23分影像,其中,图6a为(a.基于Super SloMo模型预测的影像,图6b为真实遥感影像。Super SloMo模型预测的结果与实际影像的精度计算结果如下表所示:
Figure BDA0003829620310000111
由上表可知,基于Super SloMo模型进行预训练测试出的遥感影像与真实的遥感影像计算的RMSE为0.211、CC为0.9815、PSNR为33.508、SSIM为0.7815,并且在总体效果上与真实高分辨率影像更接近,且云的运动方向保持一致,时间分辨率实现从间隔15min变成6min。
实施例四
本实施例以一个具体实例为例进行说明,本实施例提供一次北京对流初生预警个例研究,如图7所示。北京2018.7.16暴雨强度仅次于2011年“6.23暴雨”,为北京市极端暴雨事件,持续约58小时,此次由强对流引起的暴雨发生于在7.15,在7.16达到最大。本实施例利用卫星数据进行实时检测,在11:30识别到对流云,在11:40对对流云进行追踪和条件判识时,发现其满足了6个判识条件中的5个条件,具体信息如下表。因此判定该处有对流初生,即将发生强对流。并通过雷达图像进行检验发现,雷达在12:30才监测到强对流,因此通过该对流初生的追踪和判识,可实现提前50分钟预警强对流。
Figure BDA0003829620310000121
应该理解的是,虽然附图的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,其可以以其他的顺序执行。而且,附图的流程图中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,其执行顺序也不必然是依次进行,而是可以与其他步骤或者其他步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

Claims (9)

1.一种基于深度学习的对流初生预警方法,其特征在于,包括:
建立深度超分辨率模型以提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息,并对遥感影像进行训练以重构所述遥感影像的图像分辨率和时间分辨率;
获取重构的遥感影像按照预定的对流云判识策略采用云顶亮温阈值法进行对流云的初判识;
基于面积重叠法利用计算机视觉算法提取对流云以进行对流云追踪;
基于GOES-R算法对最终的对流云进行判识、依据预设的对流初生判定条件以判识对流初生;以及,
结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估,并结合雷达数据对对流初生反演结果进行检验。
2.根据权利要求1所述的基于深度学习的对流初生预警方法,其特征在于,所述建立深度超分辨率模型以提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息具体包括:
利用深度学习的非线性映射能力和对栅格数据的信息提取能力建立深度超分辨率模型;
利用深度超分辨率模型提取具备不同图像分辨率和时间分辨率的遥感影像中的有效信息;以及,
利用获取的有效信息将时空低分辨率遥感影像重构为时空高分辨率遥感影像。
3.根据权利要求1或2所述的基于深度学习的对流初生预警方法,其特征在于,所述对遥感影像进行训练以重构所述遥感影像的图像分辨率和时间分辨率具体包括:
采用ESRGAN模型对不同图像分辨率的遥感影像进行训练以将所述遥感影像重构为高图像分辨率的遥感影像,所述高图像分辨率的遥感影像保留原遥感影像的图像分辨率的通道信息;以及,
采用Super SloMo模型对不同时间分辨率的遥感影像进行训练以将所述遥感影像重构为高时间分辨率的遥感影像,所述高时间分辨率的遥感影像与雷达数据的时间分辨率保持一致。
4.根据权利要求1所述的基于深度学习的对流初生预警方法,其特征在于,所述预定的对流云判识策略具体为:对流云体的最冷温度阀值达到-55℃、对流云体的最暖温度阀值达到-10-5℃、对流云体在垂直方向上的温度等值线间隔达到1℃、对流云体极值达到3℃、且对流判识的最小面积阀值达到一个红外像元大小。
5.根据权利要求1所述的基于深度学习的对流初生预警方法,其特征在于,所述采用云顶亮温阈值法进行对流云的初判识时利用RDT算法,并采用自适应阈值选取长波红外通道阈值。
6.根据权利要求1所述的基于深度学习的对流初生预警方法,其特征在于,所述预设的对流初生判定条件具体为:10.7μm亮温≤0℃、10.7μm亮温时间变化率≤[-4℃(15min)-1]、-35℃<6.5和10.7μm亮温差<-10℃、-25℃<13.3和10.7μm亮温差<-5℃、6.5和10.7μm亮温时间变化率>[3℃(15min)-1]、且13.3和10.7μm亮温时间变化率>[3℃(15min)-1]。
7.根据权利要求1所述的基于深度学习的对流初生预警方法,其特征在于,所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估具体包括:
采用峰值信噪比和结构相似性参数评价时空降尺度的准确性;其中,
所述峰值信噪比
Figure FDA0003829620300000021
其中,
Figure FDA0003829620300000022
表示图像点颜色的最大数值;
所述结构相似性
Figure FDA0003829620300000023
其中μx、μy为均值,σx、σy为标准差,σxy为协方差,c1=(0.01L),c2=(0.03L),L对于8bit灰度图取值为255。
8.根据权利要求7所述的基于深度学习的对流初生预警方法,其特征在于,所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估还包括:
采用均方根误差和空间相关系数判断预测结果的精度;其中,
所述均方根误差
Figure FDA0003829620300000024
所述空间相关系数
Figure FDA0003829620300000025
其中,n是像素数,Oi和Si为观测值和预测值,
Figure FDA0003829620300000026
Figure FDA0003829620300000027
分别为相应的平均值。
9.根据权利要求8所述的基于深度学习的对流初生预警方法,其特征在于,所述结合高时空精度的卫星通道图像对卫星时空降尺度结果进行检验评估之后,还包括:
根据检验评估结果优化对流初生判识的参数和阈值。
CN202211076007.2A 2022-09-02 2022-09-02 一种基于深度学习的对流初生预警方法 Pending CN115691049A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211076007.2A CN115691049A (zh) 2022-09-02 2022-09-02 一种基于深度学习的对流初生预警方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211076007.2A CN115691049A (zh) 2022-09-02 2022-09-02 一种基于深度学习的对流初生预警方法

Publications (1)

Publication Number Publication Date
CN115691049A true CN115691049A (zh) 2023-02-03

Family

ID=85060796

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211076007.2A Pending CN115691049A (zh) 2022-09-02 2022-09-02 一种基于深度学习的对流初生预警方法

Country Status (1)

Country Link
CN (1) CN115691049A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117518299A (zh) * 2024-01-05 2024-02-06 南京大学 一种分类强对流临近概率预报方法、***、设备及终端
CN117592002A (zh) * 2024-01-18 2024-02-23 国家卫星气象中心(国家空间天气监测预警中心) 一种初生对流的判识方法及装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117518299A (zh) * 2024-01-05 2024-02-06 南京大学 一种分类强对流临近概率预报方法、***、设备及终端
CN117518299B (zh) * 2024-01-05 2024-03-22 南京大学 一种分类强对流临近概率预报方法、***、设备及终端
CN117592002A (zh) * 2024-01-18 2024-02-23 国家卫星气象中心(国家空间天气监测预警中心) 一种初生对流的判识方法及装置
CN117592002B (zh) * 2024-01-18 2024-04-26 国家卫星气象中心(国家空间天气监测预警中心) 一种初生对流的判识方法及装置

Similar Documents

Publication Publication Date Title
US20220043182A1 (en) Spatial autocorrelation machine learning-based downscaling method and system of satellite precipitation data
CN115691049A (zh) 一种基于深度学习的对流初生预警方法
CN107728234B (zh) 一种基于大气电场数据的雷电强度值预测方法
CN109814175B (zh) 一种基于卫星的强对流监测方法及其应用
KR102006847B1 (ko) 다중규모 광학 흐름을 이용한 레이더 초단기 강수 예측을 위한 장치 및 방법
Liang et al. A composite approach of radar echo extrapolation based on TREC vectors in combination with model-predicted winds
CN108646319B (zh) 一种短时强降雨预报方法及***
CN111983732B (zh) 一种基于深度学习的降水强度估计方法
CN113496104A (zh) 基于深度学习的降水预报订正方法及***
CN110942111B (zh) 一种识别强对流云团的方法及装置
CN113469278B (zh) 一种基于深度卷积神经网络的强天气目标识别方法
CN104237890B (zh) 一种由“列车效应”引起的暴雨识别及预报方法
CN109407177B (zh) 基于机器学习与常规气象观测的大雾识别***和应用方法
CN108254750B (zh) 一种基于雷达资料的下击暴流智能识别预警方法
CN113298295B (zh) 一种面向电力生产的气象预报***
CN110849807A (zh) 一种基于深度学习的适用于道路能见度的监测方法与***
CN116384733A (zh) 一种强对流天气下基于气象雷达的输电线路风险预警方法
CN114742206B (zh) 综合多时空尺度多普勒雷达数据的降雨强度估计方法
CN114048944A (zh) 一种暴雨诱发地质灾害下应撤离人口及损毁房屋的预估方法
CN115907574A (zh) 一种流域性暴雨洪涝承灾体重置成本遥感模拟方法
CN117710508A (zh) 基于改进条件生成对抗网络的近地表温度反演方法和装置
CN109283505B (zh) 一种修正雷达回波外推图像发散现象的方法
CN114325879A (zh) 一种基于分级概率的定量降水订正方法
CN113821895A (zh) 输电线路覆冰厚度预测模型的构建方法、装置及存储介质
CN116413835A (zh) 融合站点降水和雷达回波的ResUNet定量降水估计方法及***

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