CN113570642B - 基于背景场资料和机器学习的静止轨道卫星对流初生预警方法 - Google Patents
基于背景场资料和机器学习的静止轨道卫星对流初生预警方法 Download PDFInfo
- Publication number
- CN113570642B CN113570642B CN202110650036.4A CN202110650036A CN113570642B CN 113570642 B CN113570642 B CN 113570642B CN 202110650036 A CN202110650036 A CN 202110650036A CN 113570642 B CN113570642 B CN 113570642B
- Authority
- CN
- China
- Prior art keywords
- convection
- early warning
- satellite
- primary
- target
- 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
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000010801 machine learning Methods 0.000 title claims abstract description 23
- 230000003068 static effect Effects 0.000 title claims abstract description 8
- 230000003287 optical effect Effects 0.000 claims abstract description 28
- 238000007781 pre-processing Methods 0.000 claims abstract description 18
- 238000007637 random forest analysis Methods 0.000 claims abstract description 14
- 230000004044 response Effects 0.000 claims abstract description 6
- 238000012549 training Methods 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 18
- 238000001514 detection method Methods 0.000 claims description 16
- WWVKQTNONPWVEL-UHFFFAOYSA-N caffeic acid phenethyl ester Natural products C1=C(O)C(O)=CC=C1C=CC(=O)OCC1=CC=CC=C1 WWVKQTNONPWVEL-UHFFFAOYSA-N 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 11
- SWUARLUWKZWEBQ-UHFFFAOYSA-N phenylethyl ester of caffeic acid Natural products C1=C(O)C(O)=CC=C1C=CC(=O)OCCC1=CC=CC=C1 SWUARLUWKZWEBQ-UHFFFAOYSA-N 0.000 claims description 11
- 230000035945 sensitivity Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- 238000005381 potential energy Methods 0.000 claims description 7
- 238000007710 freezing Methods 0.000 claims description 4
- 230000008014 freezing Effects 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 4
- 238000007499 fusion processing Methods 0.000 claims description 3
- 230000035800 maturation Effects 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 239000005436 troposphere Substances 0.000 claims description 3
- 238000003908 quality control method Methods 0.000 claims description 2
- 239000013589 supplement Substances 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 description 23
- 238000009826 distribution Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 230000001186 cumulative effect Effects 0.000 description 6
- 230000006872 improvement Effects 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 5
- 238000012795 verification Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 230000000977 initiatory effect Effects 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000002329 infrared spectrum Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000208818 Helianthus Species 0.000 description 1
- 235000003222 Helianthus annuus Nutrition 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000007668 thin rolling process Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/269—Analysis of motion using gradient-based methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30192—Weather; Meteorology
-
- 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
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于背景场资料和机器学习的静止轨道卫星对流初生预警方法，所述方法包括：对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；根据判别跟踪确定的对流目标的多光谱响应，结合数值预报背景场，构造敏感因子；将敏感因子输入预先建立和训练好的基于随机森林方法的对流初生预警模型，得到分类结果，从而实现对流初生预警。本发明能够实现在保证提前于雷达回波35dBZ的前提下，对虚假目标进行过滤，从而降低对流初生预警的虚警率。
Description
技术领域
本发明涉及强对流天气预警技术领域，尤其涉及基于背景场资料和机器学习的静止轨道卫星对流初生预警方法。
背景技术
美国阿拉巴马大学(University of Alabama in Huntsville)开发的多光谱的对流初生识别预警算法，用来监测未来两个小时内可能发展成为回波强度≥35dBz强对流系统的发生，目前的版本是v2.0。该算法包括对流体判识、对流体跟踪和多光谱对流初生识别三大部分。主要判据为12条。利用MSG资料进行的算法验证结果表明，产品精度为80.75％。GOES-R对流初生算法的核心来自SATCAST(Satellite Convection Analysis AndTracking)，这是从2002年开始研发的算法，由NASA(National Aeronauticas and spaceAdministration)、NOAA(National Ocean and Atmospheric Administration)和NFS(National Science Foundation)共同资助，业务应用于美国飞控中心、一些国家天气预报服务部门、GOES-R测试平台、航空预报中心和海洋预报中心，目前版本v2.0。
美国威斯康辛大学的UWCI(University of Wisconsin Convective Initiation)算法(Sieglaff,et al,2011)的物理基础是：快速发展的对流因为云顶升高而导致云顶温度下降，因此具有较大的云顶降温率。在技术方法上，该算法提出了块平均(box-average)的概念，利用当前时刻的块平均云顶亮温与前一时刻的块平均亮温差来粗判对流，通过7条判据，进一步剔除因平流、薄卷云等影响的非对流像元，最后符合云顶降温率的像元即对流初生。该算法采用单通道数据和云类型产品作为辅助数据，因此简单、快速和高效，但缺点是在有云水平移动、薄卷云移动和复杂多层云的情况下，会产生误判。
国家卫星气象中心基于美国阿拉巴马大学的SATCAST算法，针对风云四号气象卫星也开发了对流初生预警算法。该算法主要采用了动态阈值的对流初判识、光流与面积重叠方法进行对流体跟踪，以及BP_Adaboost机器学习等技术方法，并且目前该算法产品已经进入业务使用。
值得指出的是，目前这些对流初生预警方法只采用了红外通道信息，并未引入可见光\近红外的信息。目前国家卫星气象中心业务对流初生产品解决了提前量与命中率问题，但在有云水平移动、薄卷云移动和复杂多层云的情况下，仍然受到干扰会产生虚警。
发明内容
本发明的目的在于克服现有技术缺陷，提出了基于背景场资料和机器学习的静止轨道卫星对流初生预警方法。
为了实现上述目的，本发明提出了一种基于背景场资料和机器学习的静止轨道卫星对流初生预警方法，所述方法包括：
对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；
对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；
计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；
根据判别跟踪确定的对流目标的多光谱响应，结合数值预报背景场，构造敏感因子；
将敏感因子输入预先建立和训练好的基于随机森林方法的对流初生预警模型，得到分类结果，从而实现对流初生预警。
作为上述方法的一种改进，所述对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；具体包括：
对每帧卫星云图进行通道分离处理，将卫星云图中不同通道的测值分离，每一个通道形成一幅圆盘图；再进行定位，定标以及质量控制处理得到该卫星云图预处理后的图像数据。
作为上述方法的一种改进，所述对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；具体包括：对预处理后的每帧图像均进行以下处理：
根据10.8μm亮温的云类型阈值选择目标，然后进行图像分割，对连通区域进行连续编号，形成标号图像；
对于标号图像，进行尺度检测，剔除小于预设的最小面积阈值的目标，再进行最小峰值检测，使得在目标区内的最小值与周围象元的温差在预设的温差阈值以上，得到检测后的标号图像；
对于检测后的标号图像，采用单时刻多光谱阈值法剔除其中的非对流单元，从而得到对流目标。
作为上述方法的一种改进，所述计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；具体包括：
计算两帧图像光流，根据光流对对流目标的位置进行订正，订正后计算对流目标前后时刻的重合面积，并结合前一时刻对流目标面积计算比值得到重合率，判断重合率是否大于预设比值，判断为是，则后一帧的对流目标是前一帧的后续目标。
作为上述方法的一种改进，所述敏感因子具体包括：
10.8μm亮温，用于表示云顶高；
10.8-6.2μm亮温差，用于表示云顶相对于上层对流层的高度；
10.8-12.0μm亮温差，用于表示分裂窗云顶相态；
8.5-10.8μm亮温差，用于表示云顶相态；
2*10.8-12.0-8.5μm，用于表示云顶相态；
10.8μm亮温时间变化率，用于表示云的抬升；
10.8-6.2μm亮温差时间变化率，用于表示云的抬升；
10.8-12.0μm亮温差时间变化率，用于表示云顶相态变化；
8.5-10.8μm亮温差时间变化率，用于表示云顶冻结变化；
0.64μm反照率，用于表示有效粒子半径或光学厚度；
0.64与1.6μm通道比值，用于表示云顶相态；
BLI指数，用于表示对流抬升趋势；和
CAPE对流有效位能，用于表示对流抬升的势能。
作为上述方法的一种改进，所述方法还包括对流初生预警模型的训练步骤；具体包括：
获取卫星资料、数值预报场背景资料和地基雷达资料建立训练集；
将训练集输入基于随机森林方法建立的对流初生预警模型，调整分类树数和预测变量总数，使用袋外数据样本计算估计误差，当估计误差符合预设的阈值要求时，确定预测变量数，得到训练好的对流初生预警模型。
作为上述方法的一种改进，所述获取卫星资料、数值预报场背景资料和地基雷达资料建立训练集；具体包括：
获取卫星资料、数值预报背景场资料和地基雷达资料建立从对流初生到成熟阶段整个生命史的时间序列数据集；
对时间序列数据集进行时空匹配融合处理，确认对流初生发生的位置与时间；当卫星提前于雷达进行对流初生预警时，以卫星识别的对流初生为实际的对流初生；当雷达识别的对流提前于对流初生预警时，对卫星对流观测进行反向跟踪，找到雷达对应的对流单元补充为实际的对流初生；
然后再进行反向追踪查询到对流初生触发前的背景条件，由对流触发前的背景条件和时间序列数据集构造案例数据集；
结合灾害类型数据与案例数据集进行时空匹配得到标签数据，从而建立训练集。
一种基于背景场资料和机器学习的静止轨道卫星对流初生预警系统，所述系统包括：对流初生预警模型、预处理模块、粗判别模块、判别跟踪模块、敏感因子构造模块和对流初生预警输出模块；其中，
所述预处理模块，用于对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；
所述粗判别模块，用于对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；
所述判别跟踪模块，用于计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；
所述敏感因子构造模块，用于根据判别跟踪确定的对流目标的多光谱响应，结合数值预报背景场，构造敏感因子；
所述对流初生预警输出模块，用于将敏感因子输入预先建立和训练好的基于随机森林方法的对流初生预警模型，得到分类结果，从而实现对流初生预警。
与现有技术相比，本发明的优势在于：
本发明在静止轨道气象卫星业务对流初生算法基于红外通道的基础上，引入了数值预报的最佳抬升指数与对流有效位能以及可见光通道资料，能够实现在保证提前于雷达回波35dBZ的前提下，对虚假目标进行过滤，从而降低对流初生预警的虚警率。
附图说明
图1是本发明的基于背景场资料和机器学习的静止轨道卫星对流初生预警方法总体框图；
图2是根据光流将对流进行位置订正示意图；
图3是光流订正之后面积重叠对流跟踪算法，考虑合并、分裂状况；
图4是多光谱对流判据敏感性实验与分类数据集制作模型流程图；
图5是目标识别、跟踪主体框架以及对流完整时间序列案例数据集构建示意图；
图6是基于背景场资料以及机器学习的风云四号对流初生预警技术对流识别的案例，其中，图6(a)是2020年7月6日世界时02:45案例；图6(b)是对流初生预警技术识别结果雷达回波；
图7(a)是数值预报背景场中BLI对流初生识别正负样本的概率分布图；
图7(b)是数值预报背景场中CAPE对流初生识别正负样本的概率分布图；
图8根据数值预报背景场中BLI与CAPE阈值化后对流初生识别结果对比，其中，图8(a)是基于CAPE阈值化对流初生识别结果，图8(b)是基于BLI阈值化对流初生识别结果；图8(c)是未经过BLI或CAPE阈值化后对流初生识别结果；图8(d)是雷达回波累积频数分布图；
图9展示了0.64μm与1.6μm反照率强度比分布与对流识别对应的关系，其中，图9(a)是北京时2020年8月10日11：00时0.64μm与1.6μm反照率强度比，图9(b)是北京时2020年8月10日11：00-11:30雷达回波累积分布图；图9(c)是北京时2020年8月10日12：00时0.64μm与1.6μm反照率强度比，图9(d)是北京时2020年8月10日12：00-12:30雷达回波累积分布图。框图中为0.64μm与1.6μm反照率强度比小值，而雷达回波未有反应。
具体实施方式
下面结合附图和实施例对本发明的技术方案进行详细的说明。
实施例1
如图1所示，本发明的实施例1提出了基于背景场资料和机器学习的静止轨道卫星对流初生预警方法。以风云四号气象卫星为例，进行说明。本方法主要由对流单元的初识别与跟踪模块和机器学习模型训练与验证模块两个模块构成。对流单元的初识别与跟踪模块由数据实时获取与预处理、对流初识别、对流跟踪与对流初生预警敏感性因子构造以及机器学习前向模型计算等子模块构成。机器学习模型训练与验证模块由模型训练与验证子模块构成。
步骤1)数据预处理
数据预处理进行通道分离，定位、定标等文件数据提取，质量检验和容错处理等工作。通道分离是将原始展宽图像数据中不同通道的测值分离开来，每一个通道形成一幅圆盘图，便于在后续计算中数据的反复调用。定位、定标工作的质量对卫星对流云判识的质量有非常大的影响。没有高质量的定位、定标预处理工作，对流云判识的质量是无法保证。
步骤2)对流快速粗判别
传统的地面系统对流判识算法首先根据10.8μm亮温的云类型阈值选择目标，然后进行图像分割。分割后的目标确保与其他目标分离，连续编号，形成标号图像，同时，与标号目标中象元的原始行列号也记录下来。对于所得到的标号图像，分别进行尺度检测和最小峰值检测。尺度检测检测是为了剔除过小的目标。参考RDT(Rapid DevelopmentThunderstorm)算法，取49象元为最小面积阈值，按FY-4A AGRI4km×4km空间分辨率计算，大约为784km2。最小峰值检测则要求在目标区内必须有最小值，同时，最小值与周围象元的温差在3K以上。当完成尺度检测和最小峰值检测后，再次对对流进行单时刻多光谱阈值法(PCA以后筛选出的判据条件)进行非对流单元的剔除，多光谱阈值由线下训练模块提供。
步骤3)对流跟踪
对流跟踪方法首先需要估计数据帧间的光流，对判识的目标利用光流进行位置的订正，如图2所示，再采用面积重叠的方法进行判别跟踪，如图3所示。这种方法假设在连续的两张卫星图像上，若前后两时刻的云团的位置有重合，且达到一定重合率，则可认为后一时刻的云团可能是前一时刻的云团发展而来的。
图3为面积重合法的示例，其中左侧云团为T时刻的云团，右侧为T+1时刻的云团，其重合部分如图。如果重合率(重合面积和前一个时刻的云团面积比值)大于40％，则T+1时刻的这个云团就认为是从T时刻的云团演变而来的。T时刻的云团就认为是T+1时刻云团的先行目标，T+1时刻的云团则为T时刻云团的后续目标。
步骤4)对流初生敏感因子构造
Mecikalski et al.(2010)详细描述了多通道红外判据在检测对流云发展的研究和应用情况。Siewert et al.(2009)进一步讨论了利用MSG SEVIRI(Meteosat SecondGeneration，第二代气象卫星)数据和不同的对流追踪技术检测对流初生。但他只是说明了在对流初生算法中使用多种光谱信息检测的重要性。Mecikalski et al.(2010)给出了所有可能的光谱检测并将它们从物理意义上分为三类：1)云厚度，2)云顶冻结以及3)上升气流强度。从这三个物理分类进行检测，确定哪些光谱检测是多余的，哪些包含了最重要的信息。因此，研究结果给出了所有对流事件的光谱平均值，以及哪个光谱检测最重要。对于风云四号卫星的对流初生算法，这就是多光谱检测的判据。Mecikalski et al.(2010)使用MSG其他光谱通道作为三类物理检测的扩展，其结果是从平均值中分离出标准差，然后构建临界值。由于之前的研究是使用15分钟分辨率的数据，因此在使用5分钟时间分辨率的数据时要相应做一些变化。
Mecikalski and Bedka(2006)的算法在业务中只用了当前GOES系列卫星的红外光谱通道。当前的算法可以探测出90％(当满足所有的判据阈值时)的对流初生情形，但同时有较高的虚警率(Mecikalski et al.2008)。高误报率源于基于单个像元的目标追踪和验证，如果采用目标验证方法，则可提高精度。本发明在红外光谱通道的基础上，创新地提出使用可见光与数值预报背景场的不稳定因子作为输入(具体参考表1)，进一步剔除虚假目标，从而减低对流初生识别的虚警率。
表1根据FY-4A/AGRI与背景场资料构造的对流初生的敏感因子。
| 对流初生敏感因子 | 物理意义 | 数据来源 |
| 10.8μm亮温 | 云顶高 | FY-4A/AGRI |
| 10.8-6.2μm亮温差 | 云顶相对于上层对流层的高度 | FY-4A/AGRI |
| 10.8-12.0μm亮温差 | 分裂窗云顶相态 | FY-4A/AGRI |
| 8.5-10.8μm亮温差 | 云顶相态 | FY-4A/AGRI |
| 2*10.8-12.0-8.5μm亮温差 | 云顶相态 | FY-4A/AGRI |
| 10.8μm亮温时间变化率 | 云的抬升 | FY-4A/AGRI |
| 10.8-6.2μm亮温差时间变化率 | 云的抬升 | FY-4A/AGRI |
| 10.8-12.0μm亮温差时间变化率 | 云顶相态变化 | FY-4A/AGRI |
| 8.5-10.8μm亮温差时间变化率 | 云顶冻结变化 | FY-4A/AGRI |
| 0.64μm反照率* | 有效粒子半径\光学厚度 | FY-4A/AGRI |
| 0.64与1.6μm通道比值 | 云顶相态 | FY-4A/AGRI |
| 最佳抬升指数(BLI) | 对流抬升趋势 | 数值预报预报场 |
| 对流有效位能(CAPE) | 对流抬升的势能 | 数值预报预报场 |
*表示太阳天顶角小于60°时使用。
步骤5)训练与验证数据集构造
多光谱对流判据敏感性实验与分类数据集制作。针对每一个强对流天气个例，其生命每个时刻都匹配上相应的参数来描述对流系统的特征来形成一个观测记录样本。这些参数来自卫星观测和数值预报分析场，都是中国地区的对流强降水预报中常用的参数，其中有一部分参数也被其他研究用于对流天气系统的经验预报。天气多源数据与候选对流判据时空匹配；敏感性试验：主成分量分析，挑选出候选判据(IFs)里面敏感的部分，作为实际的判据；按照下垫面、季节等因素进行分类，制作数据集。图4是多光谱对流判据敏感性实验与分类数据集制作模型流程图；图5是目标识别、跟踪主体框架以及对流完整时间序列案例数据集构建示意图。
相对于SATCAST的模糊检验方法，本发明首次提出基于目标全周期跟踪的精确检验方法。分别根据卫星资料与地基雷达资料建立对流目标检测识别与跟踪主体框架，构建包含从对流初生到成熟阶段整个生命史时间序列数据集，构造途径包括快速发展对流监测系统与人工标注方法。在分别建立的对流时间序列数据集的基础上，对卫星与地基雷达识别的对流系统进行时空匹配融合处理，确认对流初生发生的位置与时间。对于为卫星提前于雷达(35dBz)对流初生预警的情形，此时卫星识别的对流初生认定为实际的对流初生；如果出现雷达(35dBz)对流提前于初生预警的情形，需要对卫星对流观测进行反向跟踪，找到雷达35dBz对应的对流单元补充为实际的对流初生。确定对流初生发生的位置与时间后仍需要反向追踪进行延伸，以便查询到对流初生触发前时刻的背景场的条件信息。这样完成了对流触发前背景条件、对流初生、对流发展与成熟等对流完整时间序列案例数据集的构造。在此基础上，构建的灾害类型数据与完整时间序列对流案例数据集时空匹配，作为灾害类型预警的模型训练与检验标签数据
步骤6)基于随机森林的对流初生识别算法
随机森林方法的基本机制是首先根据训练数据建立多个分类树，然后通过多数投票得到最终的预测结果。随机森林算法通过用学习样本替换引导样本来开始随机选择。接下来，从每个引导示例创建树时，将在每个节点上创建一个拆分，以最恰当地分隔两个类。如果创建的引导样本数为500，则最终的分类树数也为500，形成一个林。例如，如果分类树数为500，400棵树预测是对流初生结果，剩下的100棵树预测不是对流初生，那么最终的预测结果就是对流初生。一旦建立了分类树，就可以随时将FY-4A/AGRI与背景场资料构造的对流初生的敏感因子输入随机森林模型，得到分类的结果。
在执行随机森林模型之前，需要调整两个参数。第一个参数是分类树数，分类树数的增加会减少估计误差，但计算量也会增加。有结果已经证明，500棵树足以产生稳定的对流初生的判识估计。第二个参数是预测变量总数(Q)和从Q预测变量中随机选择的预测变量数(P)。在随机森林法的训练过程中，只使用训练数据集。在确定解释变量的最佳组合或P的最佳值时，使用袋外数据(Out-Of-Bag,OOB)样本来计算估计误差。在确定了最佳解释变量组合后，才确定P。
实现效果：
图6显示了发生在北京西北部的2020年7月6日基于背景场资料以及机器学习的风云四号对流初生预警技术对流识别案例。其中，图6(a)是2020年7月6日世界时02:45案例；图6(b)是对流初生预警技术识别结果雷达回波。从图中可以看到，该事件预警的时间为时间时2:45，而雷达反应时间为4:00，与雷达相比，基于背景场资料以及机器学习的风云四号对流初生预警技术的提前量为1个小时15分钟。
图7显示了数值预报背景场中最佳抬升指数(BLI)与对流有效位能(CAPE)对流初生判识正负样本的概率分布图。其中，图7(a)是数值预报背景场中BLI对流初生识别正负样本的概率分布图；图7(b)是数值预报背景场中CAPE对流初生识别正负样本的概率分布图；根据该图可以发现，通过BLI与CAPE的阈值法处理可以剔除一部分虚假目标。
图8给出了根据数值预报背景场中最佳抬升指数或者对流有效位能阈值化后对流初生识别与未经过背景资料资料处理的对流初生识别对比结果。其中，图8(a)是基于CAPE阈值化对流初生识别结果，图8(b)是基于BLI阈值化对流初生识别结果；图8(c)是未经过BLI或CAPE阈值化后对流初生识别结果；图8(d)是雷达回波累积频数分布图；对比框结果可知，加入背景场资料能够有效剔除虚假的对流初生目标。
图9展示了0.64μm与1.6μm反照率强度比分布与对流识别对应的关系。其中，图9(a)是北京时2020年8月10日11：00时0.64μm与1.6μm反照率强度比，图9(b)是北京时2020年8月10日11：00-11:30雷达回波累积分布图；图9(c)是北京时2020年8月10日12：00时0.64μm与1.6μm反照率强度比，图9(d)是北京时2020年8月10日12：00-12:30雷达回波累积分布图。框图中为0.64μm与1.6μm反照率强度比小值，而雷达回波未有反应。对比从累积的雷达回波频数可以看出，0.64μm与1.6μm反照率强度较弱的没有雷达回波的分布，可以证明0.64μm与1.6μm反照率强度比在剔除虚假的对流初生目标方面能提供信息。
实施例2
本发明的实施例2提出了一种基于背景场资料和机器学习的静止轨道卫星对流初生预警系统，该系统包括：对流初生预警模型、预处理模块、粗判别模块、判别跟踪模块、敏感因子构造模块和对流初生预警输出模块；其中，
所述预处理模块，用于对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；
所述粗判别模块，用于对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；
所述判别跟踪模块，用于计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；
所述敏感因子构造模块，用于根据判别跟踪确定的对流目标的多光谱响应，结合数值预报背景场，构造敏感因子；
所述对流初生预警输出模块，用于将敏感因子输入预先建立和训练好的基于随机森林方法的对流初生预警模型，得到分类结果，从而实现对流初生预警。
需要说明的是，本发明以风云四号静止轨道气象卫星为例，但不限于风云四号，对于新一代静止轨道气象卫星，例如日本的葵花八号卫星、美国的GOES-R卫星及欧洲的MTG卫星都适用。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种基于背景场资料和机器学习的静止轨道卫星对流初生预警方法，所述方法包括：
对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；
对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；
计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；
根据判别跟踪确定的对流目标的多光谱响应，结合数值预报背景场，构造敏感因子；
将敏感因子输入预先建立和训练好的基于随机森林方法的对流初生预警模型，得到分类结果，从而实现对流初生预警；
所述方法还包括对流初生预警模型的训练步骤；具体包括：
获取卫星资料、数值预报场背景资料和地基雷达资料建立训练集；
将训练集输入基于随机森林方法建立的对流初生预警模型，调整分类树数和预测变量总数，使用袋外数据样本计算估计误差，当估计误差符合预设的阈值要求时，确定预测变量数，得到训练好的对流初生预警模型；
所述获取卫星资料、数值预报场背景资料和地基雷达资料建立训练集；具体包括：
获取卫星资料、数值预报背景场资料和地基雷达资料建立从对流初生到成熟阶段整个生命史的时间序列数据集；
对时间序列数据集进行时空匹配融合处理，确认对流初生发生的位置与时间；当卫星提前于雷达进行对流初生预警时，以卫星识别的对流初生为实际的对流初生；当雷达识别的对流提前于对流初生预警时，对卫星对流观测进行反向跟踪，找到雷达对应的对流单元补充为实际的对流初生；
然后再进行反向追踪查询到对流初生触发前的背景条件，由对流触发前的背景条件和时间序列数据集构造案例数据集；
结合灾害类型数据与案例数据集进行时空匹配得到标签数据，从而建立训练集。
2.根据权利要求1所述的基于背景场资料和机器学习的静止轨道卫星对流初生预警方法，其特征在于，所述对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；具体包括：
对每帧卫星云图进行通道分离处理，将卫星云图中不同通道的测值分离，每一个通道形成一幅圆盘图；再进行定位，定标以及质量控制处理得到该卫星云图预处理后的图像数据。
3.根据权利要求1所述的基于背景场资料和机器学习的静止轨道卫星对流初生预警方法，其特征在于，所述对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；具体包括：对预处理后的每帧图像均进行以下处理：
根据10.8μm亮温的云类型阈值选择目标，然后进行图像分割，对连通区域进行连续编号，形成标号图像；
对于标号图像，进行尺度检测，剔除小于预设的最小面积阈值的目标，再进行最小峰值检测，使得在目标区内的最小值与周围象元的温差在预设的温差阈值以上，得到检测后的标号图像；
对于检测后的标号图像，采用单时刻多光谱阈值法剔除其中的非对流单元，从而得到对流目标。
4.根据权利要求1所述的基于背景场资料和机器学习的静止轨道卫星对流初生预警方法，其特征在于，所述计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；具体包括：
计算两帧图像光流，根据光流对对流目标的位置进行订正，订正后计算对流目标前后时刻的重合面积，并结合前一时刻对流目标面积计算比值得到重合率，判断重合率是否大于预设比值，判断为是，则后一帧的对流目标是前一帧的后续目标。
5.根据权利要求1所述的基于背景场资料和机器学习的静止轨道卫星对流初生预警方法，其特征在于，所述敏感因子具体包括：
10.8μm亮温，用于表示云顶高；
10.8-6.2μm亮温差，用于表示云顶相对于上层对流层的高度；
10.8-12.0μm亮温差，用于表示分裂窗云顶相态；
8.5-10.8μm亮温差，用于表示云顶相态；
2*10.8-12.0-8.5μm，用于表示云顶相态；
10.8μm亮温时间变化率，用于表示云的抬升；
10.8-6.2μm亮温差时间变化率，用于表示云的抬升；
10.8-12.0μm亮温差时间变化率，用于表示云顶相态变化；
8.5-10.8μm亮温差时间变化率，用于表示云顶冻结变化；
0.64μm反照率，用于表示有效粒子半径或光学厚度；
0.64与1.6μm通道比值，用于表示云顶相态；
BLI指数，用于表示对流抬升趋势；和
CAPE对流有效位能，用于表示对流抬升的势能。
6.一种基于权利要求1的背景场资料和机器学习的静止轨道卫星对流初生预警方法的预警系统，其特征在于，所述系统包括：对流初生预警模型、预处理模块、粗判别模块、判别跟踪模块、敏感因子构造模块和对流初生预警输出模块；其中，
所述预处理模块，用于对静止轨道气象卫星连续采集的两帧卫星云图进行预处理；
所述粗判别模块，用于对预处理后的两帧图像分别进行对流快速粗判别得到对应的对流目标；
所述判别跟踪模块，用于计算预处理后的两帧图像间的光流，根据光流对对流目标进行位置的订正，再采用面积重叠的方法进行判别跟踪，确定对流目标；
所述敏感因子构造模块，用于根据判别跟踪确定的对流目标的多光谱响应，结合数值预报背景场，构造敏感因子；
所述对流初生预警输出模块，用于将敏感因子输入预先建立和训练好的基于随机森林方法的对流初生预警模型，得到分类结果，从而实现对流初生预警。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110650036.4A CN113570642B (zh) | 2021-06-10 | 2021-06-10 | 基于背景场资料和机器学习的静止轨道卫星对流初生预警方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110650036.4A CN113570642B (zh) | 2021-06-10 | 2021-06-10 | 基于背景场资料和机器学习的静止轨道卫星对流初生预警方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113570642A CN113570642A (zh) | 2021-10-29 |
| CN113570642B true CN113570642B (zh) | 2024-01-05 |
Family
ID=78161949
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110650036.4A Active CN113570642B (zh) | 2021-06-10 | 2021-06-10 | 基于背景场资料和机器学习的静止轨道卫星对流初生预警方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113570642B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104502999A (zh) * | 2014-12-10 | 2015-04-08 | 中国科学院遥感与数字地球研究所 | 一种基于静止卫星多通道数据的昼夜云检测方法与装置 |
| KR101855652B1 (ko) * | 2018-02-13 | 2018-06-25 | 대한민국 | 기계학습을 이용한 위성기반 전운량 산출 시스템 및 그 방법 |
| CN109740639A (zh) * | 2018-12-15 | 2019-05-10 | 中国科学院深圳先进技术研究院 | 一种风云卫星遥感影像云检测方法、系统及电子设备 |
| CN110567886A (zh) * | 2019-09-10 | 2019-12-13 | 西安电子科技大学 | 基于半监督空谱特征的多光谱云检测方法 |
| CN111274865A (zh) * | 2019-12-14 | 2020-06-12 | 深圳先进技术研究院 | 一种基于全卷积神经网络的遥感图像云检测方法及装置 |
| CN111638565A (zh) * | 2020-06-09 | 2020-09-08 | 云南电网有限责任公司电力科学研究院 | 一种山区暴雨监测方法 |
| CN111666903A (zh) * | 2020-06-10 | 2020-09-15 | 中国电子科技集团公司第二十八研究所 | 一种卫星云图中雷暴云团的识别方法 |
| KR20210060841A (ko) * | 2019-11-19 | 2021-05-27 | 연세대학교 산학협력단 | 위성 자료를 이용한 해상 초미세먼지 산출 방법 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10303942B2 (en) * | 2017-02-16 | 2019-05-28 | Siemens Aktiengesellschaft | Short term cloud forecast, improved cloud recognition and prediction and uncertainty index estimation |
| MX2021004475A (es) * | 2018-10-19 | 2021-06-04 | Climate Llc | Tecnica de aprendizaje automatico para identificar nubes y sombras de nubes en imagenes satelitales. |
-
2021
- 2021-06-10 CN CN202110650036.4A patent/CN113570642B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104502999A (zh) * | 2014-12-10 | 2015-04-08 | 中国科学院遥感与数字地球研究所 | 一种基于静止卫星多通道数据的昼夜云检测方法与装置 |
| KR101855652B1 (ko) * | 2018-02-13 | 2018-06-25 | 대한민국 | 기계학습을 이용한 위성기반 전운량 산출 시스템 및 그 방법 |
| CN109740639A (zh) * | 2018-12-15 | 2019-05-10 | 中国科学院深圳先进技术研究院 | 一种风云卫星遥感影像云检测方法、系统及电子设备 |
| CN110567886A (zh) * | 2019-09-10 | 2019-12-13 | 西安电子科技大学 | 基于半监督空谱特征的多光谱云检测方法 |
| KR20210060841A (ko) * | 2019-11-19 | 2021-05-27 | 연세대학교 산학협력단 | 위성 자료를 이용한 해상 초미세먼지 산출 방법 |
| CN111274865A (zh) * | 2019-12-14 | 2020-06-12 | 深圳先进技术研究院 | 一种基于全卷积神经网络的遥感图像云检测方法及装置 |
| CN111638565A (zh) * | 2020-06-09 | 2020-09-08 | 云南电网有限责任公司电力科学研究院 | 一种山区暴雨监测方法 |
| CN111666903A (zh) * | 2020-06-10 | 2020-09-15 | 中国电子科技集团公司第二十八研究所 | 一种卫星云图中雷暴云团的识别方法 |
Non-Patent Citations (3)
| Title |
|---|
| Wubin Huang ; Yingsha Jiang ; Xinwei Liu ; Yongjie Pan ; Xia Li ; Runxia Guo ; Yuxia Huang ; Bolong Duan.Classified Early-warning and Nowcasting of Hail Weather Based on Radar Products and Random Forest Algorithm.《2019 International Conference on Meteorology Observations (ICMO)》.2019,全文. * |
| 机器学习在卫星资料中的应用;郝伊一;唐伟;周勇;王喆;;卫星应用(01);全文 * |
| 覃丹宇 ; 方宗义.利用静止气象卫星监测初生对流的研究进展.《气象》.2014,全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113570642A (zh) | 2021-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102129893B1 (ko) | 딥러닝 네트워크 및 평균 이동을 기반으로 하는 선박 자동추적 방법 및 시스템 | |
| CN113705478B (zh) | 一种基于改进YOLOv5的红树林单木目标检测方法 | |
| CN111983732B (zh) | 一种基于深度学习的降水强度估计方法 | |
| CN110738127A (zh) | 基于无监督深度学习神经网络算法的安全帽识别方法 | |
| CN111507371B (zh) | 自动评估对训练图像的标签可靠性的方法和装置 | |
| CN106595551A (zh) | 基于深度学习的输电线路覆冰图像中覆冰厚度检测方法 | |
| CN110533100B (zh) | 一种基于机器学习进行cme检测和跟踪的方法 | |
| CN110728252B (zh) | 一种应用于区域人员运动轨迹监控的人脸检测方法 | |
| CN107578441B (zh) | 热带气旋中心定位的红外亮温偏差角梯度分布均匀性方法 | |
| CN108711148A (zh) | 一种基于深度学习的轮胎缺陷智能检测方法 | |
| CN109407177B (zh) | 基于机器学习与常规气象观测的大雾识别系统和应用方法 | |
| CN115437036A (zh) | 一种基于葵花卫星的对流初生预报方法 | |
| CN113076899B (zh) | 一种基于目标跟踪算法的高压输电线路异物检测方法 | |
| WO2020111934A1 (en) | A method and system for detection of natural disaster occurrence | |
| CN113838010B (zh) | 基于红外热波及卷积神经网络的复合材料损伤智能检测方法 | |
| CN111582069B (zh) | 一种空基监视平台的轨道障碍物零样本分类方法及装置 | |
| CN109960972B (zh) | 一种基于中高分辨率时序遥感数据的农林作物识别方法 | |
| CN115691049A (zh) | 一种基于深度学习的对流初生预警方法 | |
| CN112132157A (zh) | 一种基于树莓派的步态人脸融合识别方法 | |
| CN117523420B (zh) | 一种基于雷达产品数据的闪电落区识别方法及系统 | |
| CN115542430A (zh) | 一种卫星定量降水估计方法 | |
| Jones et al. | A semi-Lagrangian method for detecting and tracking deep convective clouds in geostationary satellite observations | |
| CN114627286A (zh) | 基于PSPNet和改进YOLOv4的地磅人员入侵的检测方法 | |
| CN113570642B (zh) | 基于背景场资料和机器学习的静止轨道卫星对流初生预警方法 | |
| CN116994161A (zh) | 一种基于改进YOLOv5的绝缘子缺陷检测方法 |
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 |