CN113108918B - 一种极轨气象卫星热红外遥感数据反演气温方法 - Google Patents

一种极轨气象卫星热红外遥感数据反演气温方法 Download PDF

Info

Publication number
CN113108918B
CN113108918B CN202110373988.6A CN202110373988A CN113108918B CN 113108918 B CN113108918 B CN 113108918B CN 202110373988 A CN202110373988 A CN 202110373988A CN 113108918 B CN113108918 B CN 113108918B
Authority
CN
China
Prior art keywords
data
model
air temperature
cloud
inversion
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
Application number
CN202110373988.6A
Other languages
English (en)
Other versions
CN113108918A (zh
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.)
Chengdu University of Information Technology
Original Assignee
Chengdu University of Information Technology
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 Chengdu University of Information Technology filed Critical Chengdu University of Information Technology
Priority to CN202110373988.6A priority Critical patent/CN113108918B/zh
Publication of CN113108918A publication Critical patent/CN113108918A/zh
Application granted granted Critical
Publication of CN113108918B publication Critical patent/CN113108918B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology

Landscapes

  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Radiation Pyrometers (AREA)

Abstract

本发明提供一种极轨气象卫星热红外遥感数据反演气温方法,包括以下步骤:步骤一:获取极轨卫星观测数据,匹配形成数据集;步骤二、将所述数据集分成有云的数据和晴空的数据;步骤三:利用逐步线性回归方法对所述有云的数据和晴空的数据中的变量进行变量分析,以筛选出影响气温模型构建的重要变量来作为气温反演模型的输入变量;步骤四:采用极端梯度提升算法,将筛选出的有云的数据作为云天的数据集构建云天反演模型;将筛选出的晴空的数据作为晴天的数据集构建晴空反演模型。本发明构建的模型可以快速估算高精度的有云和晴空条件下的近地面大气温度信息,且具有较高的空间分辨率。

Description

一种极轨气象卫星热红外遥感数据反演气温方法
技术领域
本发明涉及卫星气象技术技术领域,尤其涉及一种极轨气象卫星热红外遥感数据反演气温方法。
背景技术
目前,基于极轨卫星热红外数据遥感反演气温主要是在晴空条件下进行。根据相关文献表明卷云全球覆盖率大约为20%-30%,热带地区卷云覆盖率甚至高于70%。在有云的条件下,云的吸收、散射以及粒子微物理特性等会对卫星观测的辐射产生影响。因此在大气辐射传输物理过程中须考虑云的影响,以减小给卫星遥感物理反演气温的结果带来的误差。
极轨卫星提供的高空间分辨率的卫星数据能够用于近地面的气温反演,同时也能提供较高空间分辨率的气温产品。但是,由于极轨卫星通常对于同一地点的观测一天只有2次且观测区域不固定,同时极轨卫星资料的高分辨率数据量大。这些问题会给气温反演带来了一些困难。
发明内容
本发明的目的在于解决上述现有技术存在的缺陷,提供一种能够利用极轨卫星提供的高空间分辨率辐射数据,进行有云条件下的气温反演,为极轨卫星云天的气温反演提供新的方法。
一种极轨气象卫星热红外遥感数据反演气温方法,包括以下步骤:
步骤一:获取极轨卫星观测数据,匹配形成数据集;
步骤二、将所述数据集分成有云的数据和晴空的数据;
步骤三:利用逐步线性回归方法对所述有云的数据和晴空的数据中的变量进行变量分析,以筛选出影响气温模型构建的重要变量来作为气温反演模型的输入变量;
步骤四:采用极端梯度提升算法,将筛选出的有云的数据作为云天的数据集构建云天反演模型;将筛选出的晴空的数据作为晴天的数据集构建晴空反演模型;
步骤五:利用所述云天反演模型或晴空反演模型获得极轨气象卫星反演气温的估算结果;
步骤六:对所述估算结果的精度进行评估;
步骤七:根据评估结果确定极轨气象卫星反演的气温。
进一步地,如上所述的极轨气象卫星热红外遥感数据反演气温方法,所述极轨卫星观测数据包括:
极轨卫星观测热红外亮温、数值预报模式资料、数字高程模型、台站气温数据、和台站高度数据。
进一步地,如上所述的极轨气象卫星热红外遥感数据反演气温方法,所述数值预报模式资料包括:如GFS预报气温和高度数据。
进一步地,如上所述的极轨气象卫星热红外遥感数据反演气温方法,步骤二中所述将所述数据集分成有云的数据和晴空的数据包括:将所述数据集通过L2云检测,建立近地面气温反演数据集,该数据集包括云的数据和晴空的数据。
进一步地,如上所述的极轨气象卫星热红外遥感数据反演气温方法,所述步骤三包括:
将所述数据集作为机器学习模型的输入变量,将实际台站气温作为输出变量,分别建立训练集、验证集和测试数据集;
所述训练集用于建立云天反演模型或晴空反演模型;所述验证集评估模型的性能和调整模型参数;所述测试数据集用于评估模型的精度。
进一步地,如上所述的极轨气象卫星热红外遥感数据反演气温方法,所述评估模型的精度包括:
利用模型估算的气温与相同时间实际台站测得的气温进行对比,计算误差对模型进行精度评估:
Figure BDA0003010447260000031
Ta是模型反演的气温值,Tsta是气象台站实测得的气温,N是所选用的样本个数,RMSE为均方根误差,其能够反映出反演的台站气温的精确度。
有益效果:
本发明构建的模型可以快速估算高精度的有云和晴空条件下的近地面大气温度信息,且具有较高的空间分辨率。我国自主研发的极轨卫星FY-3D具有250m的空间分辨率的探测能力,能够获得极高空间分辨率的卫星热红外亮温数据,而数值预报本身能够预测气温但是空间分辨率通常高达20km以上,而机器学习能够学习不同数据的特征。因此,我们通过机器学习算法将热红外亮温卫星数据以及数值预报数据作为输入数据,采用机器学习算法,建立有云和晴空条件下气温反演模型,从而快速获得高精度的有云和晴空条件下的近地面大气温度信息。
附图说明
图1是XGB模型训练路线图;
图2是基于本发明算法FY-3D MERSI热红外数据气温估算与气象台站实测数据的对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明首要解决的问题是如何对极轨卫星观测资料进行时间和区域的固定,以及对大量的卫星数据进行去噪,去异常值等预处理。本发明将极轨卫星资料通过数据选择,数据清洗以及观测区域订正后进行下一步的气温反演工作。
在有云条件下,热红外通道的辐射会受到云层的吸收、散射等的影响,从而会对云天气温反演增加复杂度和困难度。本发明通过如下三个步骤进行有云条件下的气温反演研究。
(1)利用极轨卫星观测亮温(如FY-3D)、气象台站观测资料、数值预报模式资料(如GFS)以及数字高程模型(DEM),匹配形成数据集,再通过L2云检测,建立近地面气温反演数据集;
(2)利用逐步线性回归方法对反演数据集的变量进行变量分析,例如辐射亮温变量及其他要素等不同变量对近地面气温的影响程度进行分析,了解各个变量对于模型本身的重要性,筛选出影响气温模型构建的重要变量,从而作为气温反演模型的输入变量,这能够提高构建的气温反演模型的效率和有效性;
(3)利用极端梯度提升算法(eXtreme Gradient Boosting,缩写为XGB)建立气温反演模型。基于(2)步骤确定的极轨卫星观测热红外亮温(如FY-3D)、数值预报模式资料(如GFS预报气温和高度数据),数字高程模型(DEM)和台站高度数据等关键因子和气象台站的实测近地面大气温度信息。将极轨卫星观测热红外亮温(如FY-3D)、数值预报模式资料(如GFS预报气温和高度数据),数字高程模型(DEM),台站海拔高度作为机器学习模型的输入变量,将实际台站气温作为输出变量,分别建立3种数据集,即训练集、验证集(通常将原始数据80%划分为训练集,20%为验证集)和测试数据集(未被训练和验证过程使用的数据)。机器学习算法会在训练集上训练XGB模型,建立基于极轨气象卫星热红外亮温数据反演云天条件和晴空下的近地面气温模型。在验证集上评估模型性能和调整参数。调整到最优的参数后就在测试集上最后测试,评估XGB模型精度。极端梯度提升(XGB)算法是为每一个样本赋上一个权重值,通过每一步训练中得到的模型会对该节点数据进行估计,使得数据点的估计有对有错,估计后会增加分错的点的权重,减少分对点的权重,这样使得某些经常被分错的点会被赋上一个非常高的权重。然后进行n次迭代得到n个简单的分类器,加入正则项来控制模型的复杂度,防止过拟合,最后将它们组合起来得到一个最终的XGB模型。
通过云检测,将极轨气象卫星热红外亮温数据分成有云的数据和晴空的数据,采用极端梯度提升(XGB)算法,将有云的数据作为云天的数据集构建反演模型;如果判别晴空,则将识别为晴空的数据用来构建晴空反演模型,从而分别得到云天条件和晴空下的近地面气温2个模型。
(4)利用(3)步骤建立的气温反演模型和验证集,通过极轨卫星观测热红外亮温(如FY-3D)、数值预报模式资料(如GFS预报气温和高度数据)、数字高程模型(DEM)和台站海拔高度作为机器学习模型的输入变量,获得有云条件下极轨气象卫星反演气温的估算结果。利用模型估算的气温与相同时间实际台站测得的气温进行对比,计算误差对模型进行精度评估。精度评估主要通过计算统计参数进行评估,如均方根误差(RMSE)。
Figure BDA0003010447260000051
Ta是模型反演的气温值,Tsta是气象台站实测得的气温,N是所选用的样本个数。RMSE能够反映出反演的台站气温的精确度,较低的RMSE表示反演的台站气温与实际台站气温的差异较小,即表示反演气温的精度较高。
实施例说明
收集了2019年湖南区域97个气象台站的近地面大气温度实测数据,并同步收集了FY-3DMERSI卫星同步观测的热红外亮温数据、美国数值预报(GFS)的气温数据。
利用线性回归和随机森林分析方法分析MERSI热红外亮温数据、像元经纬度、观测角、GFS气温数据等因素对近地面大气温度估算的影响,选择气温估算精度的关键因子并评估了各个因子对整个模型性能的影响程度。表1给出了通过逐步回归分析确定的主要影响因子。
基于上一步骤确定的关键因子和气象台站的实测近地面大气温度信息,对输入、输出数据构建形成训练数据集,利用XGB算法构建有云条件下的近地面大气温度估算模型,整个气温估算流程如图1所示,极轨卫星观测热红外亮温(如FY-3D)、数值预报模式资料(如GFS预报气温和高度数据),数字高程模型(DEM)和台站高度数据作为输入数据,输出数据是反演的气象台站的近地面温度信息,通过云检测进行了数据晴空和云天数据集的分类,分别构建了晴空和云天的基于XGB算法的气温反演模型,两个模型,在模型训练和测试的时候,输入的数据和输出的数据都是一样的。最后选择2020年7月3日的FY-3DMERSI卫星同步观测的热红外亮温数据、美国数值预报(GFS)的气温等数据对算法进行了测试。测试结果表明本发明构建的模型可以快速估算高精度的有云和晴空条件下的近地面大气温度信息,且空间分辨率能够达到250m,如图2所示。
表1 XGB算法训练相关主要参数的信息
Figure BDA0003010447260000061
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (6)

1.一种极轨气象卫星热红外遥感数据反演气温方法,其特征在于,包括以下步骤:
步骤一:获取极轨卫星观测数据,匹配形成数据集;
步骤二、将所述数据集分成有云的数据和晴空的数据;
步骤三:利用逐步线性回归方法对所述有云的数据和晴空的数据中的变量进行变量分析,以筛选出影响气温模型构建的重要变量来作为气温反演模型的输入变量;
步骤四:采用极端梯度提升算法,将筛选出的有云的数据作为云天的数据集构建云天反演模型;将筛选出的晴空的数据作为晴天的数据集构建晴空反演模型;
步骤五:利用所述云天反演模型或晴空反演模型获得极轨气象卫星反演气温的估算结果;
步骤六:对所述估算结果的精度进行评估;
步骤七:根据评估结果确定极轨气象卫星反演的气温;
所述步骤四包括:利用极端梯度提升算法建立气温反演模型依赖于极轨卫星观测热红外亮温、数值预报模式资料,将极轨卫星观测热红外亮温、数值预报模式资料,数字高程模型,台站海拔高度作为机器学习模型的输入变量,将实际台站气温作为输出变量,分别建立3种数据集,即训练集、验证集和测试数据集;机器学习算法会在训练集上训练XGB模型,建立基于极轨气象卫星热红外亮温数据反演云天条件和晴空下的近地面气温模型。
2.根据权利要求1所述的极轨气象卫星热红外遥感数据反演气温方法,其特征在于,所述极轨卫星观测数据包括:
极轨卫星观测热红外亮温、数值预报模式资料、数字高程模型、台站气温数据、和台站高度数据。
3.根据权利要求2所述的极轨气象卫星热红外遥感数据反演气温方法,其特征在于,所述数值预报模式资料包括:如数值预报模式资料预报气温和高度数据。
4.根据权利要求1所述的极轨气象卫星热红外遥感数据反演气温方法,其特征在于,
步骤二中所述将所述数据集分成有云的数据和晴空的数据包括:将所述数据集通过L2云检测,建立近地面气温反演数据集,该数据集包括云的数据和晴空的数据。
5.根据权利要求1所述的极轨气象卫星热红外遥感数据反演气温方法,其特征在于,所述步骤四包括:
将所述数据集作为机器学习模型的输入变量,将实际台站气温作为输出变量,分别建立训练集、验证集和测试数据集;
所述训练集用于建立云天反演模型或晴空反演模型;所述验证集评估模型的性能和调整模型参数;所述测试数据集用于评估模型的精度。
6.根据权利要求5所述的极轨气象卫星热红外遥感数据反演气温方法,其特征在于,所述评估模型的精度包括:
利用模型估算的气温与相同时间实际台站测得的气温进行对比,计算误差对模型进行精度评估:
Figure FDA0003625317730000021
Ta是模型反演的气温值,Tsta是气象台站实测得的气温,N是所选用的样本个数;RMSE为均方根误差,其能够反映出反演的台站气温的精确度。
CN202110373988.6A 2021-04-07 2021-04-07 一种极轨气象卫星热红外遥感数据反演气温方法 Active CN113108918B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110373988.6A CN113108918B (zh) 2021-04-07 2021-04-07 一种极轨气象卫星热红外遥感数据反演气温方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110373988.6A CN113108918B (zh) 2021-04-07 2021-04-07 一种极轨气象卫星热红外遥感数据反演气温方法

Publications (2)

Publication Number Publication Date
CN113108918A CN113108918A (zh) 2021-07-13
CN113108918B true CN113108918B (zh) 2022-06-21

Family

ID=76714539

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110373988.6A Active CN113108918B (zh) 2021-04-07 2021-04-07 一种极轨气象卫星热红外遥感数据反演气温方法

Country Status (1)

Country Link
CN (1) CN113108918B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114354510B (zh) * 2022-01-12 2024-06-04 国家卫星气象中心(国家空间天气监测预警中心) 一种基于极轨气象卫星mersi的白天云相态反演方法
CN115032720A (zh) * 2022-07-15 2022-09-09 国网上海市电力公司 基于随机森林的多模式集成预报在地面气温预报中的应用
CN116450700B (zh) * 2023-03-09 2024-02-13 中国农业科学院农业资源与农业区划研究所 极轨卫星地表温度时间归一化方法、装置及电子设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830448A (zh) * 2012-08-21 2012-12-19 北京琨奇电子系统有限公司 一种微波高光谱晴空定标方法、装置及系统
CN103293084A (zh) * 2013-05-08 2013-09-11 南京大学 基于多光谱气象卫星信息的海雾全天时全天候反演方法
CN104820250A (zh) * 2015-04-14 2015-08-05 国家卫星气象中心 一种极轨气象卫星virr海洋上云检测的处理方法
WO2018145229A1 (zh) * 2017-02-10 2018-08-16 广西壮族自治区气象减灾研究所 一种近地面气温的大面积精确反演方法
CN109580003A (zh) * 2018-12-18 2019-04-05 成都信息工程大学 一种静止气象卫星热红外数据估算近地面大气温度方法
CN109871637A (zh) * 2019-03-06 2019-06-11 成都信息工程大学 一种云天条件下近地面气温估算方法
CN110688606A (zh) * 2019-08-30 2020-01-14 中国科学院遥感与数字地球研究所 一种热红外遥感反演云微物理参数的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106019314B (zh) * 2016-08-05 2018-11-02 成都信息工程大学 一种卷云条件下短波红外卫星遥感二氧化碳的方法
CN109459392B (zh) * 2018-11-06 2019-06-14 南京农业大学 一种基于无人机多光谱影像的水稻地上部生物量估测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830448A (zh) * 2012-08-21 2012-12-19 北京琨奇电子系统有限公司 一种微波高光谱晴空定标方法、装置及系统
CN103293084A (zh) * 2013-05-08 2013-09-11 南京大学 基于多光谱气象卫星信息的海雾全天时全天候反演方法
CN104820250A (zh) * 2015-04-14 2015-08-05 国家卫星气象中心 一种极轨气象卫星virr海洋上云检测的处理方法
WO2018145229A1 (zh) * 2017-02-10 2018-08-16 广西壮族自治区气象减灾研究所 一种近地面气温的大面积精确反演方法
CN109580003A (zh) * 2018-12-18 2019-04-05 成都信息工程大学 一种静止气象卫星热红外数据估算近地面大气温度方法
CN109871637A (zh) * 2019-03-06 2019-06-11 成都信息工程大学 一种云天条件下近地面气温估算方法
CN110688606A (zh) * 2019-08-30 2020-01-14 中国科学院遥感与数字地球研究所 一种热红外遥感反演云微物理参数的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Retrieval of Total Precipitable Water from Himawari-8 AHI Data: A Comparison of Random Forest, Extreme Gradient Boosting, and Deep Neural Network;Yeonjin Lee 等;《remote sensing》;20190724;全文 *
基于FY-3C MWHTS的台风降水反演算法研究;李娜 等;《遥感技术与应用》;20191031;全文 *

Also Published As

Publication number Publication date
CN113108918A (zh) 2021-07-13

Similar Documents

Publication Publication Date Title
CN113108918B (zh) 一种极轨气象卫星热红外遥感数据反演气温方法
CN111651941B (zh) 一种全球电离层电子总含量预测的算法
CN109580003B (zh) 一种静止气象卫星热红外数据估算近地面大气温度方法
CN110058328B (zh) 东北夏季降水多模式组合降尺度预测方法
CN109784552B (zh) 一种基于Re-ESF算法的空间变系数PM2.5浓度估算模型的构建方法
WO2020043027A1 (zh) 一种大气污染预测中对区域边界优化的方法
CN110174359B (zh) 一种基于高斯过程回归的航空高光谱影像土壤重金属浓度评估方法
CN108875254A (zh) 一种大气温湿廓线的一维变分反演方法
CN113348471B (zh) 一种大气污染预测中对区域边界优化的方法
CN111737913B (zh) 一种基于云水含量反演的mwhts晴空观测亮温选择方法
CN111797131A (zh) 一种基于遥感降水产品的极端降水区域频率分析方法
CN113743013A (zh) 一种基于XGBoost气温预测数据修正的方法
CN111091082A (zh) 一种基于高分辨率遥感数据的流域污染溯源方法
CN113408111B (zh) 大气可降水量反演方法及系统、电子设备和存储介质
CN105974495A (zh) 利用分类拟合法预判目标区域未来平均云量的方法
CN114880933A (zh) 一种基于再分析资料的无探空站点地基微波辐射计大气温湿廓线反演方法及系统
CN113468799A (zh) 一种静止气象卫星观测中获取近地面pm2.5浓度的方法及系统
Ruiz-Arias et al. Direct normal irradiance modeling: Evaluating the impact on accuracy of worldwide gridded aerosol databases
CN114595876A (zh) 一种区域风场预测模型生成方法和装置、电子设备
Chen et al. Proactive quality control: Observing system simulation experiments with the Lorenz’96 model
CN112560270A (zh) 一种化学危害同化系统
CN117409334B (zh) 一种基于静止卫星遥感数据的气溶胶类型识别方法
CN116432032A (zh) 基于多源数据和机器学习的气象数据异常事件识别方法
CN113311509B (zh) 一种基于神经网络的mwhts对海面气压的敏感性测试方法
CN110222301B (zh) 一种雾霾条件下地表太阳短波辐射计算方法

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