一种表面温度高温和低温数据集重构方法与流程

1.本发明涉及近地面气温数据重构技术领域，尤其涉及一种表面温度高温和低温数据集重构方法。


背景技术：

2.近地面气温是反映全球气候变化的重要变量，显著影响地球各圈层能量和物质的循环转换，特别是每日近地面最高气温，最低气温和平均气温对研究大气循环转换，气候变化以及极端天气事件的强度和频率具有重要意义；
3.目前地面观测数据在单点上具有时间连续性，但受到站点分布不均和对空间异质性表征差的影响缺乏空间连续性，卫星遥感数据可以进行大范围空间观测，但受到传感器过境时间和云雨天气的影响缺乏时间连续性，同化数据可以获得大区域长时间序列的气温数据，但其空间分辨率较低且目前缺乏近地面最高气温和最低气温数据，因此，本发明提出一种表面温度高温和低温数据集重构方法以解决现有技术中存在的问题。


技术实现要素：

4.针对上述问题，本发明的目的在于提出一种表面温度高温和低温数据集重构方法，该表面温度高温和低温数据集重构方法区分了晴天和非晴天天气状态，分别建立气温变化模型计算得到每日最高和最低气温，通过校正方程进一步提高数据精度，采用多种评价指标对近地面气温的时空变化趋势进行分析，以cmfd数据和era5数据集作为再分析数据源，具有较高的空间分辨率和可靠的精度，可准确捕捉区域气温变化情况，便于气候相关分析和其他地表驱动因子的研究。
5.为实现本发明的目的，本发明通过以下技术方案实现：一种表面温度高温和低温数据集重构方法，包括以下步骤：
6.步骤一、将每日天气划分为晴天天气和非晴天天气，并分析近地面气温的两种估算策略；
7.步骤二、先根据晴天状态下气温变化的周期性和非对称性建立晴天状态下的气温分段正弦函数，模拟近地面气温的变化情况估算近地面的最高气温和最低气温并校正，然后估算非晴天天气下的近地面的最高气温和最低气温并校正；
8.步骤三、将校正后的近地面最该气温和最低气温作为每日气温变化的限定阈值，并根据站点实测数据对近地面平均气温输出值进行多元线性回归校正以提高精度得到每日的近地面平均气温数据集；
9.步骤四、通过平均大气温度的垂直衰减率对各种气象站的近地面气温的数据进行校正；
10.步骤五、将决定系数、平均绝对误差和均方根误差三个指标作为评估指标，并利用评估指标对原始输出的气温数据集和校正后数据集与站点实测数据分别进行精度验证；
11.步骤六、通过对输出近地面气温数据与站点实测数据进行比较，评估各区域的近
地面最高气温、最低气温和平均气温的变化趋势，通过各气候变量进一步检验数据集的效果和区域适用性；
12.步骤七、对校正前和修正后的近地面最高气温、最低气温和平均气温及各站点实测数据进行精度验证，并对近地面的最高气温和最低气温的波动趋势进行分析。
13.进一步改进在于：所述步骤一中对于天气划分是借助modis自带的质量控制字段和利用站点实测数据进行判别。
14.进一步改进在于：所述步骤一中近地面气温的两种估算策略为当像元位置有对应气象站点及相邻站点间欧式距离小于0.3
°
时，利用多时相和空间相关性的方法对像元区域进行数据填补；针对站点分布较为稀疏且相邻两站点之间欧氏距离大于0.3
°
的地区，利用era5数据集逐小时数据对天气状态进行划分。
15.进一步改进在于：所述利用era5逐小时数据对天气状态进行划分时先基于era5数据集与cmfd的分析数据源逐三小时气温数据的空间相关关系构建了一个有效的降尺度过程，再将era5数据集在cmfd数据的辅助下进行空间降尺度。
16.进一步改进在于：所述步骤二中建立的晴天状态下的气温分段正弦函数由下式1和式2表示
[0017][0018][0019]
其中h
min
是近地面最低气温t
min
出现的时间，h
max
是近地面最高气温t
max
出现的时间，h
o
为输入时间，a
t
和b
t
是未知参数。
[0020]
进一步改进在于：所述步骤四中校正方法为先将观测气温按照统一的气温垂直递减率统一到海平面高度上，然后通过与海平面相对应的温度完成数据校正和插值过程，最后再校正到温度数据所对应的海拔高度。
[0021]
进一步改进在于：所述步骤六中采用线性回归、相关系数和t检验的方法对近地面平均气温进行分析，通过线性回归和相关系数得到近地面平均气温的年际变化速率和相关性。
[0022]
进一步改进在于：所述步骤七中对修正后的近地面最高气温、最低气温和平均气温及各站点实测数据进行精度验证时基于步骤四中的校正方法对对温度进行进一步的校正并将校正后的数据与实测数据对比后得到精度验证结果。
[0023]
本发明的有益效果为：本发明区分了晴天和非晴天天气状态，分别建立气温变化模型计算得到每日最高和最低气温，通过校正方程进一步提高数据精度，采用多种评价指标对近地面气温的时空变化趋势进行分析，以cmfd数据和era5数据集作为再分析数据源，具有较高的空间分辨率和可靠的精度，可准确捕捉区域气温变化情况，便于气候相关分析和其他地表驱动因子的研究。
附图说明
[0024]
图1为本发明操作流程图。
[0025]
图2为本发明天气划分流程图。
[0026]
图3为本发明空间降尺度过程图。
具体实施方式
[0027]
为了加深对本发明的理解，下面将结合实施例对本发明做进一步详述，本实施例仅用于解释本发明，并不构成对本发明保护范围的限定。
[0028]
实施例1
[0029]
根据图1、2、3所示，本实施例提供了一种表面温度高温和低温数据集重构方法，包括以下步骤：
[0030]
步骤一、借助modis自带的质量控制字段和利用站点实测数据将每日天气划分为晴天天气和非晴天天气，其中modis自带的质量控制字段用以判别2002
‑
2021年时间段的天气状态，第一种为当像元位置有对应气象站点及相邻站点间欧式距离小于0.3
°
时，利用多时相的方法对像元区域进行数据填补；
[0031]
站点实测数据用以判别1979
‑
2001年时间段的天气状态，并分析近地面气温的两种估算策略，针对站点分布较为稀疏且相邻两站点之间欧氏距离大于0.3
°
的地区，先基于era5数据集与cmfd的分析数据源逐三小时气温数据的空间相关关系构建了一个有效的降尺度过程，再将era5数据集在cmfd数据的辅助下进行空间降尺度，再利用era5数据集逐小时数据对天气状态进行划分；
[0032]
步骤二、先根据晴天状态下气温变化的周期性和非对称性建立晴天状态下的气温分段正弦函数，模拟近地面气温的变化情况估算近地面的最高气温和最低气温并校正，然后估算非晴天天气下的近地面的最高气温和最低气温并校正，其中建立的晴天状态下的气温分段正弦函数由下式1和式2表示
[0033][0034][0035]
其中h
min
是近地面最低气温t
min
出现的时间，h
max
是近地面最高气温t
max
出现的时间，h
o
为输入时间，a
t
和b
t
是未知参数；
[0036]
步骤三、将校正后的近地面最该气温和最低气温作为每日气温变化的限定阈值，并根据站点实测数据对近地面平均气温输出值进行多元线性回归校正以提高精度得到每日的近地面平均气温数据集；
[0037]
步骤四、通过平均大气温度的垂直衰减率对各种气象站的近地面气温的数据进行校正，先将观测气温按照统一的气温垂直递减率统一到海平面高度上，然后通过与海平面相对应的温度完成数据校正过程，最后再校正到温度数据所对应的海拔高度；
[0038]
步骤五、将决定系数、平均绝对误差和均方根误差三个指标作为评估指标，并利用
评估指标对原始输出的气温数据集和校正后数据集与站点实测数据分别进行精度验证；
[0039]
步骤六、通过对输出近地面气温数据与站点实测数据进行比较，评估各区域的近地面最高气温、最低气温和平均气温的变化趋势，通过各气候变量进一步检验数据集的效果和区域适用性，采用线性回归、相关系数和t检验的方法对近地面平均气温进行分析，通过线性回归和相关系数得到近地面平均气温的年际变化速率和相关性，由下式3和式4计算给出
[0040][0041][0042]
其中n表示时间序列长度的总年数，i表示年份，t
i
表示第i年的近地面平均气温，r表示气温与时间序列的相关关系，r值介于
‑
1～1之间；
[0043]
步骤七、对校正前和修正后的近地面最高气温、最低气温和平均气温及各站点实测数据进行精度验证，对修正后的近地面最高气温、最低气温和平均气温及各站点实测数据进行精度验证时基于步骤四中的校正方法对温度进行进一步的校正并将校正后的数据与实测数据对比后得到精度验证结果，并对近地面的最高气温和最低气温的波动趋势进行分析。
[0044]
实施例2
[0045]
根据图1、2、3所示，本实施例提供了一种表面温度高温和低温数据集重构方法，包括以下步骤：
[0046]
步骤一、借助modis自带的质量控制字段和利用站点实测数据将每日天气划分为晴天天气和非晴天天气，其中modis自带的质量控制字段用以判别2002
‑
2021年时间段的天气状态，第一种为当像元位置有对应气象站点及相邻站点间欧式距离小于0.3
°
时，利用空间相关性的方法对像元区域进行数据填补；
[0047]
站点实测数据用以判别1979
‑
2001年时间段的天气状态，并分析近地面气温的两种估算策略，针对站点分布较为稀疏且相邻两站点之间欧氏距离大于0.3
°
的地区，先基于era5数据集与cmfd的分析数据源逐三小时气温数据的空间相关关系构建了一个有效的降尺度过程，再将era5数据集在cmfd数据的辅助下进行空间降尺度，再利用era5数据集逐小时数据对天气状态进行划分；
[0048]
步骤二、先根据晴天状态下气温变化的周期性和非对称性建立晴天状态下的气温分段正弦函数，模拟近地面气温的变化情况估算近地面的最高气温和最低气温并校正，然后估算非晴天天气下的近地面的最高气温和最低气温并校正，其中建立的晴天状态下的气温分段正弦函数由下式1和式2表示
[0049]
[0050][0051]
其中h
min
是近地面最低气温t
min
出现的时间，h
max
是近地面最高气温t
max
出现的时间，h
o
为输入时间，a
t
和b
t
是未知参数；；
[0052]
步骤三、将校正后的近地面最该气温和最低气温作为每日气温变化的限定阈值，并根据站点实测数据对近地面平均气温输出值进行多元线性回归校正以提高精度得到每日的近地面平均气温数据集；
[0053]
步骤四、通过平均大气温度的垂直衰减率对各种气象站的近地面气温的数据进行校正，先将观测气温按照统一的气温垂直递减率统一到海平面高度上，然后通过与海平面相对应的温度完成数据插值过程，最后再校正到温度数据所对应的海拔高度；
[0054]
步骤五、将决定系数、平均绝对误差和均方根误差三个指标作为评估指标，并利用评估指标对原始输出的气温数据集和校正后数据集与站点实测数据分别进行精度验证；
[0055]
步骤六、通过对输出近地面气温数据与站点实测数据进行比较，评估各区域的近地面最高气温、最低气温和平均气温的变化趋势，通过各气候变量进一步检验数据集的效果和区域适用性，采用线性回归、相关系数和t检验的方法对近地面平均气温进行分析，通过线性回归和相关系数得到近地面平均气温的年际变化速率和相关性，由下式3和式4计算给出
[0056][0057][0058]
其中n表示时间序列长度的总年数，i表示年份，t
i
表示第i年的近地面平均气温，r表示气温与时间序列的相关关系，r值介于
‑
1～1之间；
[0059]
步骤七、对校正前和修正后的近地面最高气温、最低气温和平均气温及各站点实测数据进行精度验证，对修正后的近地面最高气温、最低气温和平均气温及各站点实测数据进行精度验证时基于步骤四中的校正方法对温度进行进一步的校正并将校正后的数据与实测数据对比后得到精度验证结果，并对近地面的最高气温和最低气温的波动趋势进行分析。
[0060]
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。