融合卫星与浮标观测的海洋三维温度遥感超分重建方法

1.本发明涉及遥感信息处理与应用领域，具体涉及一种融合卫星与浮标观测的海洋三维温度遥感超分重建方法。


背景技术：

2.海洋内部动态过程复杂，argo浮标观测数据的时间序列太短，无法支持长期的海洋现象分析。因此，对全球海洋三维温度数据集进行长时序重建势在必行。20世纪70年代以来，卫星遥感技术提供了高时空分辨率、长时间序列的海洋观测，促进了大气和海洋科学的发展。但是，由于海洋的物理性质，卫星遥感的电磁波不能直接穿透海洋，因此无法直接获得海洋内部观测。


技术实现要素：

3.本发明的目的在于提供一种融合卫星与浮标观测的海洋三维温度遥感超分重建方法，该方法可以实现卫星观测从海表到海洋内部的扩展，提高海洋三维温度观测的空间分辨率。
4.为了实现上述目的，本发明采用的技术方案是：一种融合卫星与浮标观测的海洋三维温度遥感超分重建方法，包括以下步骤：
5.步骤s1：获取全球海洋多源卫星遥感观测数据与海洋内部温度数据并进行预处理；
6.步骤s2：根据argo实测格网数据的经纬度坐标，将每个格网点上的海洋内部温度量数据st与海表遥感及时空参量一一对应，分别标记为遥感参量特征矩阵x、实测温度st标签矩阵y，合并为矩阵a后划分为训练数据集、验证数据集和测试数据集；
7.步骤s3：基于cnn-lstm网络模型构建基于多源卫星遥感观测和海洋内部argo浮标实测数据的智能遥感时序超分重建模型，并输入1
°
分辨率的训练数据集进行模型训练，基于验证数据集调节模型的网络结构和超参数，从而得到最优的智能遥感时序超分重建模型，然后将测试数据集作为输入数据，用均方根误差衡量最终模型的泛化能力；
8.步骤s4：将设定年限内高空间分辨率0.25
°×
0.25
°
的海表遥感数据组为长时序特征矩阵x'，输入到训练好的智能遥感时序超分重建模型中，重建考虑时空特征的长时序高空间分辨率全球海洋三维温度数据；
9.步骤s5：选择r2和nrmse为评价指标，利用评价指标和空间分布格局对重建精度进行定量和定性分析。
10.进一步地，所述步骤s1具体包括以下步骤：
11.步骤s11：获取全球海洋多源卫星遥感观测数据与海洋内部温度数据，全球海洋多源卫星遥感观测数据包括海表高度ssh、海表温度sst和海表风场ssw，其中海表风场需分解成水平分量ussw和垂直分量vssw，海洋内部温度数据包括argo实测格网数据和en4廓线数据；
12.步骤s12：对获取的数据进行包括质量控制、去除异常值和陆地空值的预处理；基于矢量数据栅格化的思想将en4廓线数据处理为0.25
°
空间分辨率的稀疏格网化数据产品；将获取的所有数据的时间分辨率统一处理为月，将所有海表遥感观测数据的空间分辨率0.25
°×
0.25
°
利用最优插值法统一为argo格网观测的空间分辨率1
°×1°
，并保留0.25
°
的海表参量用于时序超分重建；基于argo观测数据的空间覆盖范围对1
°
海表卫星观测数据进行掩膜，使1
°
数据的空间覆盖统一；基于woa提供的0.25
°
掩膜数据对0.25
°
海表参量和海洋内部温度数据进行空间覆盖统一。
13.进一步地，所述步骤s2具体包括以下步骤：
14.步骤s21：根据argo实测格网数据的经纬度坐标，将海洋内部温度量数据st与输入模型的海表遥感及时空参量在1
°×1°
的格网上一一对应，分别标记为特征矩阵x、标签矩阵y；
15.步骤s22：按照cnn-lstm网络模型的数据接受形式将x和y进行纵向合并得到二维矩阵a，基于时间序列将矩阵a划分为训练数据集、验证数据集和测试数据集。
16.进一步地，所述步骤s3具体包括以下步骤：
17.步骤s31：基于cnn-lstm网络模型构建基于多源卫星遥感观测和海洋内部argo浮标实测数据的智能遥感时序超分重建模型；模型训练时输入1
°
分辨率的海表遥感参量和海洋内部温度数据，模型利用cov1d模块进行空间特征的提取，在训练和测试中，端到端映射以低分辨率图像为训练输入，重建时输出高分辨率图像，即训练时模型接受的单个特征通道大小为x*1，训练后的模型可输入同样降维处理的0.25
°
分辨率的海表遥感参量，但单个特征通道大小变为16x*1；
18.步骤s32：在智能遥感时序超分重建模型中，利用cnn模块从卫星海表数据中提取与格网中心点温度相关的空间特征，利用lstm模块学习每个网格点温度的时间相关性，通过验证数据集和格网搜索方法调节隐藏层层数、池化层类型、激活函数类型、丢弃率和学习率以获得最优的网络结构和超参数，得到最优的智能遥感时序超分重建模型；
19.步骤s33：在模型的网络结构和超参数调整后，将测试数据集输入最终得到的模型中，利用精度评价因子均方根误差rmse衡量最终模型的泛化能力。
20.进一步地，所述步骤s4具体包括以下步骤：
21.步骤s41：获取1993至2020年的高空间分辨率0.25
°×
0.25
°
的海表遥感观测和时空位置共7个参量的数据，标记为长时序特征矩阵x'；
22.步骤s42：将长时序特征矩阵x'输入步骤s3训练好的模型，利用模型在训练数据集中学习到的时空依赖关系，由x'直接预测得到长时序高分辨率全球海洋次表层温度矩阵y'，实现高时空分辨率卫星遥感观测向海洋内部映射，重建考虑时空特征的长时序高空间分辨率全球海洋三维温度数据，并命名为dors-0.25
°
。
23.进一步地，所述步骤s5具体包括以下步骤：
24.步骤s51：选择r2和nrmse为评价指标，其表达式如下：
[0025][0026]
[0027]
其中，yi为实测温度值，为模型估算的温度值，为实测温度的平均值，y
max
和y
min
为实测数据集中最大最小值，n为样本总数；
[0028]
步骤s52：利用评价指标和空间分布格局对dors-0.25
°
的重建精度进行定量和定性分析，并基于dors-0.25
°
数据集追踪近三十年来全球海洋上层2000m温度变化趋势。
[0029]
与现有技术相比，本发明具有以下有益效果：本发明利用现有的丰富的海表卫星遥感观测和海洋内部实测数据，结合深海遥感技术，基于顾及时序依赖性和空间特征的卷积长短期记忆神经网络cnn-lstm方法，重构长时序高空间分辨率全球海洋三维温度数据集，实现了卫星观测从海表到海洋内部的扩展，填补pre-argo时期海洋内部观测数据的空白，将全球海洋三维温度观测的空间分辨率提高至0.25
°
。
附图说明
[0030]
图1为本发明实施例的方法实现流程示意图；
[0031]
图2为本发明实施例的智能遥感时序超分重建模型的网络结构图；
[0032]
图3为本发明实施例的dors-0.25
°
与en4稀疏格网化数据的精度评价图；
[0033]
图4为本发明实施例的结果展示图，依次为2016年1月100m、500m深度层，argo st(左图)和dors-0.25
°
st(右图)；
[0034]
图5为本发明实施例的全球海洋上层2000m温度时序演变趋势展示图。
具体实施方式
[0035]
下面结合附图及实施例对本发明做进一步说明。
[0036]
应该指出，以下详细说明都是示例性的，旨在对本技术提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本技术所属技术领域的普通技术人员通常理解的相同含义。
[0037]
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本技术的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0038]
如图1所示，本实施例提供了一种融合卫星与浮标观测的海洋三维温度遥感超分重建方法，包括以下步骤：
[0039]
步骤s1：获取全球海洋多源卫星遥感观测数据与海洋内部温度数据并进行预处理。
[0040]
在本实施例中，所述步骤s1具体包括以下步骤：
[0041]
步骤s11：获取全球海洋多源卫星遥感观测数据与海洋内部温度数据，全球海洋多源卫星遥感观测数据包括海表高度ssh、海表温度sst和海表风场ssw，其中海表风场需分解成水平分量ussw和垂直分量vssw，上述数据的空间分辨率均为0.25
°×
0.25
°
；时空位置参量包括经度lon(longitude)、纬度lat(latitude)、时间doy(day ofyear)(当前时相处于一年之中的天数)；海洋内部温度数据包括由实测浮标插值得到的argo实测格网数据和来源于世界海洋数据库wod的en4廓线数据。
[0042]
步骤s12：对获取的数据进行预处理：首先对多源海表遥感和海洋内部格网观测数
据集进行质量控制、去除异常值和陆地空值处理；其次是基于矢量数据栅格化的思想将en4廓线数据处理为0.25
°
空间分辨率的稀疏格网化数据产品；然后是将上述数据的时间分辨率统一处理为月，将所有海表遥感观测数据的空间分辨率(0.25
°×
0.25
°
)利用最优插值法统一为argo格网观测的空间分辨率(1
°×1°
)，并保留0.25
°
的海表参量用于时序超分重建；最后基于argo观测数据的空间覆盖范围对1
°
海表卫星观测数据和海洋内部温度数据集进行掩膜，使1
°
数据的空间覆盖统一，基于woa提供的0.25
°
掩膜数据对0.25
°
海表参量和海洋内部温度数据进行空间覆盖统一。
[0043]
步骤s2：根据argo实测格网数据的经纬度坐标，将每个格网点上的海洋内部温度量数据st与海表遥感及时空参量一一对应，分别标记为特征矩阵x、标签矩阵y，合并为矩阵a后划分为训练数据集、验证数据集和测试数据集。
[0044]
在本实施例中，所述步骤s2具体包括以下步骤：
[0045]
步骤s21：将上述海表遥感参量和海洋内部温度数据降维至一维，并根据argo实测格网数据的经纬度坐标，将1
°×1°
的海洋内部温度量数据st与输入模型的海表遥感及时空参量在1
°×1°
的每个格网上一一对应，分别标记为一维的遥感参量特征矩阵x和实测温度st标签矩阵y。
[0046]
步骤s22：按照cnn-lstm网络模型的数据接受形式将x和y进行纵向合并得到二维矩阵a，特征通道数量为7，基于时间序列将矩阵a划分为训练数据集、验证数据集和测试数据集。
[0047]
步骤s3：基于cnn-lstm网络模型构建基于多源卫星遥感观测和海洋内部argo浮标实测数据的智能遥感时序超分重建模型，并输入1
°
分辨率的训练数据集进行模型训练，基于验证数据集调节模型的网络结构和超参数，从而得到最优的智能遥感时序超分重建模型，然后将测试数据集作为输入数据，用均方根误差衡量最终模型的泛化能力。本实施例中构建的智能遥感时序超分重建模型如图2所示。
[0048]
在本实施例中，所述步骤s具体包括以下步骤：
[0049]
步骤s31：在模型训练时输入一一对应的1
°
分辨率的海表遥感参量和海洋内部温度数据，cnn-lstm模型利用cov1d模块进行空间特征的提取，与cov2d相比虽降低了局部特征提取的能力但规避了空间二维数据分辨率不一致导致的输入与输出无法对应的问题。训练和测试中，端到端映射以低分辨率图像为训练输入，重建时输出高分辨率图像，即训练时模型接受的单个特征通道大小为x*1，训练后的模型可输入同样降维处理的0.25
°
分辨率的海表遥感参量，区别在于单个特征通道大小变为16x*1。
[0050]
步骤s32：在智能遥感时序超分重建模型中，利用cnn模块从卫星海表数据中提取与格网中心点温度相关的空间特征，利用lstm模块学习每个网格点温度的时间相关性，通过验证数据集和格网搜索方法调节隐藏层层数、池化层类型、激活函数类型、丢弃率和学习率以获得最优的网络结构和超参数，得到最优的智能遥感时序超分重建模型。
[0051]
步骤s33：在模型网络结构和超参数调整后，将测试数据集输入最终得到的模型中，利用精度评价因子均方根误差rmse衡量最终模型的泛化能力。
[0052]
步骤s4：将设定年限内高空间分辨率(0.25
°×
0.25
°
)的海表遥感数据组为长时序特征矩阵x'，输入到训练好的智能遥感时序超分重建模型中，重建考虑时空特征的长时序高空间分辨率全球海洋三维温度数据。
[0053]
在本实施例中，所述步骤s4具体包括以下步骤：
[0054]
步骤s41：获取1993至2020年的高空间分辨率(0.25
°×
0.25
°
)的海表遥感观测和时空位置共7个参量的数据，标记为长时序特征矩阵x'。
[0055]
步骤s42：将长时序特征矩阵x'输入步骤s3训练好的模型，利用模型在训练数据集中学习到的时空依赖关系，由x'直接预测得到长时序高分辨率全球海洋次表层温度矩阵y'，实现高时空分辨率卫星遥感观测向海洋内部映射，重建考虑时空特征的长时序高空间分辨率全球海洋三维温度数据，并命名为dors-0.25
°
。
[0056]
步骤s5：选择r2和nrmse为评价指标，利用评价指标和空间分布格局对重建精度进行定量和定性分析。
[0057]
在本实施例中，所述步骤s5具体包括以下步骤：
[0058]
步骤s51：选择r2和nrmse为评价指标，其表达式如下：
[0059][0060][0061]
其中，yi为实测温度值，为模型估算的温度值，为实测温度的平均值，y
max
和y
min
为实测数据集中最大最小值，n为样本总数。
[0062]
步骤s52：利用评价指标和空间分布格局对dors-0.25
°
的重建精度进行定量和定性分析。定量分析是将dors-0.25
°
与实测稀疏格网化en4数据集进行不同深度层位的精度对比(如图3所示)，定性分析是分析对比dors-0.25
°
与实测格网化argo数据集在100和500m深度层位的空间分布(如图4所示)，结果表明dors-0.25
°
各个层位重建结果均保持较高水平。经过上述精度对比后，基于遥感卫星视角超分重建的dors-0.25
°
数据集追踪近三十年来全球海洋上层2000m温度变化趋势，开展海洋内部暖化过程分析(如图5所示)。
[0063]
本发明提供了一种融合卫星与浮标观测的海洋三维温度遥感超分重建方法，该方法通过重建长时序高空间分辨率全球海洋三维温度数据集(dors-0.25
°
)，实现了卫星观测从海表到海洋内部的扩展，填补了pre-argo时期海洋内部观测数据的空白，将全球海洋三维温度观测的空间分辨率提高至0.25
°
，dors-0.25
°
可支持全球海洋温度异常和海洋热含量变化分析以及其他海洋气候变化研究，为海洋内部动力过程、全球变暖和可持续发展目标(sdgs)提供有价值和有前景的数据支持。
[0064]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0065]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产
生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0066]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0067]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0068]
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

技术特征：
1.一种融合卫星与浮标观测的海洋三维温度遥感超分重建方法，其特征在于，包括以下步骤：步骤s1：获取全球海洋多源卫星遥感观测数据与海洋内部温度数据并进行预处理；步骤s2：根据argo实测格网数据的经纬度坐标，将每个格网点上的海洋内部温度量数据st与海表遥感及时空参量一一对应，分别标记为遥感参量特征矩阵x、实测温度st标签矩阵y，合并为矩阵a后划分为训练数据集、验证数据集和测试数据集；步骤s3：基于cnn-lstm网络模型构建基于多源卫星遥感观测和海洋内部argo浮标实测数据的智能遥感时序超分重建模型，并输入1
°
分辨率的训练数据集进行模型训练，基于验证数据集调节模型的网络结构和超参数，从而得到最优的智能遥感时序超分重建模型，然后将测试数据集作为输入数据，用均方根误差衡量最终模型的泛化能力；步骤s4：将设定年限内高空间分辨率0.25
°×
0.25
°
的海表遥感数据组为长时序特征矩阵x'，输入到训练好的智能遥感时序超分重建模型中，重建考虑时空特征的长时序高空间分辨率全球海洋三维温度数据；步骤s5：选择r2和nrmse为评价指标，利用评价指标和空间分布格局对重建精度进行定量和定性分析。2.根据权利要求1所述的融合卫星与浮标观测的海洋三维温度遥感超分重建方法，其特征在于，所述步骤s1具体包括以下步骤：步骤s11：获取全球海洋多源卫星遥感观测数据与海洋内部温度数据，全球海洋多源卫星遥感观测数据包括海表高度ssh、海表温度sst和海表风场ssw，其中海表风场需分解成水平分量ussw和垂直分量vssw，海洋内部温度数据包括argo实测格网数据和en4廓线数据；步骤s12：对获取的数据进行包括质量控制、去除异常值和陆地空值的预处理；基于矢量数据栅格化的思想将en4廓线数据处理为0.25
°
空间分辨率的稀疏格网化数据产品；将获取的所有数据的时间分辨率统一处理为月，将所有海表遥感观测数据的空间分辨率0.25
°×
0.25
°
利用最优插值法统一为argo格网观测的空间分辨率1
°×1°
，并保留0.25
°
的海表参量用于时序超分重建；基于argo观测数据的空间覆盖范围对1
°
海表卫星观测数据进行掩膜，使1
°
数据的空间覆盖统一；基于woa提供的0.25
°
掩膜数据对0.25
°
海表参量和海洋内部温度数据进行空间覆盖统一。3.根据权利要求1所述的融合卫星与浮标观测的海洋三维温度遥感超分重建方法，其特征在于，所述步骤s2具体包括以下步骤：步骤s21：根据argo实测格网数据的经纬度坐标，将海洋内部温度量数据st与输入模型的海表遥感及时空参量在1
°×1°
的格网上一一对应，分别标记为特征矩阵x、标签矩阵y；步骤s22：按照cnn-lstm网络模型的数据接受形式将x和y进行纵向合并得到二维矩阵a，基于时间序列将矩阵a划分为训练数据集、验证数据集和测试数据集。4.根据权利要求1所述的融合卫星与浮标观测的海洋三维温度遥感超分重建方法，其特征在于，所述步骤s3具体包括以下步骤：步骤s31：基于cnn-lstm网络模型构建基于多源卫星遥感观测和海洋内部argo浮标实测数据的智能遥感时序超分重建模型；模型训练时输入1
°
分辨率的海表遥感参量和海洋内部温度数据，模型利用cov1d模块进行空间特征的提取，在训练和测试中，端到端映射以低分辨率图像为训练输入，重建时输出高分辨率图像，即训练时模型接受的单个特征通道大
小为x*1，训练后的模型可输入同样降维处理的0.25
°
分辨率的海表遥感参量，但单个特征通道大小变为16x*1；步骤s32：在智能遥感时序超分重建模型中，利用cnn模块从卫星海表数据中提取与格网中心点温度相关的空间特征，利用lstm模块学习每个网格点温度的时间相关性，通过验证数据集和格网搜索方法调节隐藏层层数、池化层类型、激活函数类型、丢弃率和学习率以获得最优的网络结构和超参数，得到最优的智能遥感时序超分重建模型；步骤s33：在模型的网络结构和超参数调整后，将测试数据集输入最终得到的模型中，利用精度评价因子均方根误差rmse衡量最终模型的泛化能力。5.根据权利要求1所述的融合卫星与浮标观测的海洋三维温度遥感超分重建方法，其特征在于，所述步骤s4具体包括以下步骤：步骤s41：获取1993至2020年的高空间分辨率0.25
°×
0.25
°
的海表遥感观测和时空位置共7个参量的数据，标记为长时序特征矩阵x'；步骤s42：将长时序特征矩阵x'输入步骤s3训练好的模型，利用模型在训练数据集中学习到的时空依赖关系，由x'直接预测得到长时序高分辨率全球海洋次表层温度矩阵y'，实现高时空分辨率卫星遥感观测向海洋内部映射，重建考虑时空特征的长时序高空间分辨率全球海洋三维温度数据，并命名为dors-0.25
°
。6.根据权利要求5所述的融合卫星与浮标观测的海洋三维温度遥感超分重建方法，其特征在于，所述步骤s5具体包括以下步骤：步骤s51：选择r2和nrmse为评价指标，其表达式如下：和nrmse为评价指标，其表达式如下：其中，y
i
为实测温度值，为模型估算的温度值，为实测温度的平均值，y
max
和y
min
为实测数据集中最大最小值，n为样本总数；步骤s52：利用评价指标和空间分布格局对dors-0.25
°
的重建精度进行定量和定性分析，并基于dors-0.25
°
数据集追踪近三十年来全球海洋上层2000m温度变化趋势。

技术总结
本发明涉及一种融合卫星与浮标观测的海洋三维温度遥感超分重建方法，包括：获取海洋多源卫星遥感观测数据与海洋内部温度数据并预处理；根据Argo实测格网数据的坐标，将海洋内部温度量数据与海表遥感及时空参量对应，标记、合并为相应矩阵后划分训练、验证和测试数据集；基于CNN-LSTM网络模型构建基于多源卫星遥感观测和海洋内部Argo浮标实测数据的智能遥感时序超分重建模型，输入1


技术研发人员：苏华 姜进文 王安
受保护的技术使用者：福州大学
技术研发日：2023.07.14
技术公布日：2023/10/15