CN113051824B - 基于剖面观测数据的海洋内部温度场重建方法 - Google Patents
基于剖面观测数据的海洋内部温度场重建方法 Download PDFInfo
- Publication number
- CN113051824B CN113051824B CN202110324151.2A CN202110324151A CN113051824B CN 113051824 B CN113051824 B CN 113051824B CN 202110324151 A CN202110324151 A CN 202110324151A CN 113051824 B CN113051824 B CN 113051824B
- Authority
- CN
- China
- Prior art keywords
- data
- grid
- profile
- ocean
- training
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K3/00—Thermometers giving results other than momentary value of temperature
- G01K3/08—Thermometers giving results other than momentary value of temperature giving differences of values; giving differentiated values
- G01K3/14—Thermometers giving results other than momentary value of temperature giving differences of values; giving differentiated values in respect of space
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明涉及一种基于剖面观测数据的海洋内部温度场重建方法,包括以下步骤:步骤S1:获取历史海洋观测的剖面数据,并对剖面数据进行格网化预处理,得到格网化剖面观测数据;步骤S2:通过预设比例随机选取格网化剖面观测数据作为交叉验证集,其余作为输入数据;步骤S3:根据输入数据集生成训练集和测试集;步骤S4:基于卷积自编码器神经网络对训练集进行迭代训练,构建卷积自编码器模型;步骤S5:将测试集输入到自编码器模型中,重建长时间序列的海洋温度数据。本发明提高了海洋内部温度时空观测分辨率与完整性,在海洋遥感及全球气候变化领域有较高的应用价值。
Description
技术领域
本发明涉及遥感信息处理与应用领域,具体涉及一种基于剖面观测数据的海洋内部温度场重建方法。
背景技术
全球变暖的大背景下,人类越来越重视海洋内部热力结构变化。海洋温度和海洋热含量是表征海洋热力结构的重要参量,其中海洋温度是对全球海洋热力结构变化最直接的反映,是海气交互作用强烈程度的指标。海洋内部温度是海洋物理性质的关键变量,反映着海洋密度和层化结构,是全球洋流循环的主要驱动因素之一。近年来卫星遥感技术在海洋领域应用成效显著,但目前遥感在海洋领域大多应用于海表参量反演,无法直接观测海洋内部信息。纵观人类对海观测历史,主要的剖面观测仪器主要是一下几种,1967年以前主要以南森瓶(Nansen bottle)和机械深海温盐测量器(Mechanical bathythermographs,MBTs)为主;1968年后至21世纪初主要以抛弃式测温仪器(Xpendable Bathy Thermograph,XBT)为主;21世纪初期至今,以Argo浮标为主。受不同历史时期技术条件的限制,海洋观测剖面空间分布不均,时空分辨率较低。海洋内部动力过程具有时空多维度多尺度特点,基于现有历史观测资料填补历史观测的空白,重建长时序高时空覆盖率的海洋温度场对了解全球海洋的变化机制和过程非常重要。
发明内容
有鉴于此,本发明的目的在于提供一种基于剖面观测数据的海洋内部温度场重建方法,为海洋环境与全球气候变化研究提供有力的数据支撑。
为实现上述目的,本发明采用如下技术方案:
一种基于剖面观测数据的海洋内部温度场重建方法,包括以下步骤:
步骤S1:获取历史海洋观测的剖面数据,并对剖面数据进行格网化预处理,得到格网化剖面观测数据;
步骤S2:通过预设比例随机选取格网化剖面观测数据作为交叉验证集,其余作为输入数据;
步骤S3:根据输入数据集生成训练集和测试集;
步骤S4:基于卷积自编码器神经网络对训练集进行迭代训练,构建卷积自编码器模型;
步骤S5:将测试集输入到自编码器模型中,重建长时间序列的海洋温度数据。
进一步的,所述步骤S1具体为:
获取海洋EN4 profile、Argo浮标数据、IAP时间序列格网温度和海表温度遥感数据,并进行空间范围统一处理;
统计IAP数据不同深度不同时相的极值,对EN4 profile进行极值检查;
将海表温度遥感数据、IAP数据、EN4 profile和Argo数据的时间分辨率统一为每月,利用最邻近插值将海表温度遥感数据空间分辨率统一为1°×1°;
剖面数据格网化处理中,0m、5m、和10m的格网深度范围拓展2m,20m和30m格网深度范围拓展5m,50~1000m的格网深度拓展10m。将检索到的格网中观测值取平均值代表格网温度值。
进一步的,所述步骤S2具体为:
每一时相二维稀疏分布的格网化剖面观测,构成三维时间序列的格网化剖面观测数据;
通过matlab中find算法提取出时间序列格网化剖面观测数据中非空值的格网索引,再以randperm算法预设比例从上述非空值格网索引中随机选取,保留为交叉验证集;
交叉验证集保存以后,将其对应的时间序列格网化剖面观测数据中的值设为空值,生成输入数据集。
进一步的,所述步骤S3具体为:将基于输入数据集生成的归一化温度随机掩膜设为空值,并随机加入5%高斯噪声,生成训练数据集,基于完整的输入数据集生成归一化的输入变量构成测试集。
进一步的,所述模型训练和测试过程中,根据输入的位置、时间、协方差,通过二维卷积学习主要的时空特征。
进一步的,所述步骤S4具体为:
将训练集分为多个mini-batch,再输入到卷积自编码器神经网络中进行训练;
参数调优防止过拟合或欠拟合,通过调节L2正则化、丢弃正则化和学习率进行模型参数优化。
本发明与现有技术相比具有以下有益效果:
本发明以时间序列数据为基础,充分利用的剖面数据含有多种海洋观测仪器的数据,如MBTs、XBT、CTD及Argo等,基于剖面观测数据的时空分布特征,采用卷积自编码器神经网络方法,建立高效的时序重建模型,时空重建了0-1000m海洋温度,重建结果有良好的精度,能准确、直观的呈现全球海洋0-1000m温度结构的空间分布特征和时间序列变化。
附图说明
图1为本发明的方法流程示意图;
图2 为本发明实施例的精度评价图;
图3~6为本发明实施例的结果展示图。其中图3~图6依次为2006年5月、2019年9月30m、200m和600m、1000m深度层,Argo格网温度(左图)和模型重建的温度(右图Reconstruction)对比图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
请参照图1,本发明提供一种基于剖面观测数据的海洋内部温度场重建方法,包括以下步骤:
步骤S1:获取历史海洋观测的剖面数据EN4 profiles,并对剖面数据进行格网化预处理,得到格网化剖面观测数据(EN4 Gridded measure data,EN4-GMD)。获取Argo浮标数据、中国科学院IAP格网化温度数据和遥感海表温度数据。
步骤S11:获取海洋EN4 profile、Argo浮标数据和IAP时间序列格网温度数据,并将这些数据的空间范围统一(经度:180°W~180°E,纬度:77.5°S~78.5°N);
步骤S12:统计IAP数据不同深度不同时相的极值,对EN4 profile进行极值检查,超出极值的认定为异常值,将其去除;
步骤S13:将海表温度遥感数据、IAP数据、EN4 profile和Argo数据的时间分辨率统一为每月,利用最邻近插值将海表温度遥感数据空间分辨率统一为1°×1°;
步骤S14:剖面数据格网化处理中,0m、5m、和10m的格网深度范围拓展2m(如5m检索3m-7m的剖面),20m和30m格网深度范围拓展5m,50~1000m的格网深度拓展10m。将检索到的格网中观测值取平均值代表格网温度值,其中0m深度由1955年1月至1981年12月的EN4剖面数据和1982年1月至2019年12月的遥感海表温度数据共同组成。
步骤S2:将每一时相非空的EN4-GMD随机选取1%作为交叉验证集(EN4-GMD-VS),其余作为输入数据(EN4-GMD-IN)。EN4-GMD-VS不参与整个重建运算,仅用于重建结果精度验证。
步骤S21:将每一时相空间稀疏分布的格网化剖面观测数据,构成三维时间序列的EN4-GMD;
步骤S22:通过matlab中find算法提取出时间序列EN4-GMD中有值格网的索引,再以randperm算法从上述每一时相有值格网索引中随机选取1%,保留为交叉验证集(EN4-GMD-VS);
步骤S23:EN4-GMD-VS保存以后,将其对应的时间序列EN4-GMD中的值设为空值,生成输入数据集(EN4-GMD-IN);
步骤S3:基于输入的时空数据EN4-GMD-IN生成训练集和测试集,将EN4-GMD-IN生成的归一化温度随机掩膜其中一部分设为空值,并随机加入5%高斯噪声,生成训练数据集,基于完整的EN4-GMD-IN生成归一化输入变量构成测试集;
步骤S31:EN4-GMD-IN数据以1955年1月为基准,计算每一时相对应的时间序列数,构成时间time输入变量。经度0-360°维度-77.5-78.5,通过meshgrid方法构成二维lon和lat变量;
步骤S32:将经纬度、温度及时间等变量进行归一化处理,温度距平减掉时空平均值,时间time变量的归一化以时序数比12的正余弦值表示;
步骤S33:上述归一化输入变量直接构成测试集,训练数据集则需将归一化温度随机掩膜(设为空值)其中一部分,并加入5%的高斯噪声。
步骤S4:使用卷积自编码器神经网络对训练集进行迭代训练,建立卷积自编码器模型;
步骤S41:因海洋时空数据量较大,需要设置batch-size将完整训练集分为多个mini-batch,再输入到卷积自编码器神经网络中进行训练,这样可以降低模型训练时所占内存;
步骤S42:参数调优防止过拟合或欠拟合,通过调节L2正则化、丢弃正则化和学习率等进行模型参数优化。
步骤S5:将测试集输入到自编码器模型中,重建长时间序列的海洋温度数据,得到0-1000m海洋内部温度分布。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。
Claims (2)
1.一种基于剖面观测数据的海洋内部温度场重建方法,其特征在于,包括以下步骤:
步骤S1:获取历史海洋观测的剖面数据,并对剖面数据进行格网化预处理,得到格网化剖面观测数据;
步骤S2:通过预设比例随机选取格网化剖面观测数据作为交叉验证集,其余作为输入数据集;
步骤S3:根据输入数据集生成训练集和测试集;
步骤S4:基于卷积自编码器神经网络对训练集进行迭代训练,构建卷积自编码器模型;
步骤S5:将测试集输入到卷积自编码器模型中,重建长时间序列的海洋内部温度数据;
所述步骤S1具体为:
获取海洋EN4 profile、Argo浮标数据、IAP时间序列格网温度和海表温度遥感数据,并进行空间范围统一处理;
统计IAP时间序列格网温度不同深度不同时相的极值,对EN4 profile进行极值检查;
将海表温度遥感数据、IAP时间序列格网温度、EN4 profile和Argo浮标数据的时间分辨率统一为每月,利用最邻近插值将海表温度遥感数据空间分辨率统一为1°×1°;
剖面数据格网化处理中,0m、5m、和10m的格网深度范围拓展2m,20m和30m格网深度范围拓展5m,50~1000m的格网深度拓展10m;
将检索到的格网中观测值取平均值代表格网温度值;
所述步骤S2具体为:
将每一时相空间稀疏分布的格网化剖面观测数据,构成三维时间序列的格网化剖面观测数据;
通过matlab中find算法提取出时间序列格网化剖面观测数据中非空值的格网索引,再以randperm算法根据 预设比例从上述非空值的格网索引中随机选取,保留为交叉验证集;
交叉验证集保存以后,将其对应的时间序列格网化剖面观测数据中的值设为空值,生成输入数据集;
所述步骤S3具体为:将基于输入数据集生成的归一化的温度随机掩膜设为空值,并随机加入5%高斯噪声,生成训练数据集,基于完整的输入数据集生成归一化的输入变量构成测试集;
所述步骤S4具体为:
将训练集分为多个mini-batch,再输入到卷积自编码器神经网络中进行训练;
参数调优防止过拟合或欠拟合,通过调节L2正则化、丢弃正则化和学习率进行模型参数优化。
2.据权利要求1所述的基于剖面观测数据的海洋内部温度场重建方法,其特征在于,所述模型训练和测试过程中,根据输入的位置、时间、协方差,通过二维卷积学习时空特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110324151.2A CN113051824B (zh) | 2021-03-26 | 2021-03-26 | 基于剖面观测数据的海洋内部温度场重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110324151.2A CN113051824B (zh) | 2021-03-26 | 2021-03-26 | 基于剖面观测数据的海洋内部温度场重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113051824A CN113051824A (zh) | 2021-06-29 |
CN113051824B true CN113051824B (zh) | 2022-06-14 |
Family
ID=76515370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110324151.2A Active CN113051824B (zh) | 2021-03-26 | 2021-03-26 | 基于剖面观测数据的海洋内部温度场重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113051824B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115096274B (zh) * | 2022-06-16 | 2023-08-18 | 中国人民解放军61540部队 | 海洋观测系统和mapping方案评价方法及系统 |
CN115238569B (zh) * | 2022-06-21 | 2023-04-28 | 中国人民解放军61540部队 | 基于信息嵌合深度ffnn模型的温盐重构方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013089840A1 (en) * | 2011-12-14 | 2013-06-20 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Multi-layer ocean vertical gradient constraints |
CN109543768A (zh) * | 2018-11-30 | 2019-03-29 | 福州大学 | 基于多源卫星遥感的海洋内部温盐信息智能提取方法 |
CN109543356A (zh) * | 2019-01-07 | 2019-03-29 | 福州大学 | 考虑空间非平稳性的海洋内部温盐结构遥感反演方法 |
CN110516856A (zh) * | 2019-08-12 | 2019-11-29 | 中国海洋大学 | 基于卷积神经网络的估算海洋次表层温度的方法 |
-
2021
- 2021-03-26 CN CN202110324151.2A patent/CN113051824B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013089840A1 (en) * | 2011-12-14 | 2013-06-20 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Multi-layer ocean vertical gradient constraints |
CN109543768A (zh) * | 2018-11-30 | 2019-03-29 | 福州大学 | 基于多源卫星遥感的海洋内部温盐信息智能提取方法 |
CN109543356A (zh) * | 2019-01-07 | 2019-03-29 | 福州大学 | 考虑空间非平稳性的海洋内部温盐结构遥感反演方法 |
CN110516856A (zh) * | 2019-08-12 | 2019-11-29 | 中国海洋大学 | 基于卷积神经网络的估算海洋次表层温度的方法 |
Non-Patent Citations (2)
Title |
---|
Estimating Subsurface Thermohaline Structure of the Global Ocean Using Surface Remote Sensing Observations;Su Hua et.al;《REMOTE SENSING》;20190826;全文 * |
基于LightGBM的全球海洋次表层温盐遥感预测;张天一;《遥感学报》;20201031;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113051824A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110361778B (zh) | 一种基于生成对抗网络的地震数据重建方法 | |
CN113051824B (zh) | 基于剖面观测数据的海洋内部温度场重建方法 | |
GB2547816B (en) | Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation | |
CN108981957A (zh) | 基于自组织神经网络及经验正交函数的水下温度场重构方法 | |
Konomi et al. | Adaptive Bayesian nonstationary modeling for large spatial datasets using covariance approximations | |
CN112287294B (zh) | 一种基于深度学习的时空双向土壤含水量插值方法 | |
Chen | Logistic models of fractal dimension growth of urban morphology | |
Heaton et al. | Constructing valid spatial processes on the sphere using kernel convolutions | |
CN116152465A (zh) | 一种基于多级插值法的海洋环境场重建方法 | |
Pauthenet et al. | Four-dimensional temperature, salinity and mixed-layer depth in the Gulf Stream, reconstructed from remote-sensing and in situ observations with neural networks | |
Lubkov et al. | Method for reconstructing the monthly mean water transparencies for the northwestern part of the Black Sea as an example | |
Chen et al. | TemproNet: A transformer-based deep learning model for seawater temperature prediction | |
CN117271979A (zh) | 一种基于深度学习的赤道印度洋表层海流流速预测方法 | |
CN116050460B (zh) | 基于注意力神经网络的气温数据空间插值方法 | |
CN116861640A (zh) | 基于时空插值的高分辨温盐数据生成方法及装置 | |
CN115564989A (zh) | 面向土地利用分类的随机森林算法 | |
CN114155354B (zh) | 基于图卷积网络的电容层析成像重建方法与装置 | |
CN110837924B (zh) | 一种水质浊度预测方法 | |
CN115186755A (zh) | 一种考虑相关性的温盐数据垂向梯度误差检测方法 | |
Berretta et al. | Adaptive Environmental Sampling: The Interplay Between Geostatistics and Geometry. | |
CN113011086A (zh) | 一种基于ga-svr算法森林生物量的估测方法 | |
Chen et al. | The improved kriging interpolation algorithm for local underwater terrain based on fractal compensation | |
Amendola et al. | A neural network ensemble downscaling system (SIBILLA) for seasonal forecasts over Italy: winter case studies | |
Zheng et al. | A multi-factor data-driven prediction model for cyanobacteria blooms in lakes and reservoirs | |
CN117494573B (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 |