CN109829547B - 一种基于深度学习的一维综合孔径微波辐射计sst反演方法 - Google Patents
一种基于深度学习的一维综合孔径微波辐射计sst反演方法 Download PDFInfo
- Publication number
- CN109829547B CN109829547B CN201811547210.7A CN201811547210A CN109829547B CN 109829547 B CN109829547 B CN 109829547B CN 201811547210 A CN201811547210 A CN 201811547210A CN 109829547 B CN109829547 B CN 109829547B
- Authority
- CN
- China
- Prior art keywords
- sea surface
- surface temperature
- deep learning
- encoder
- self
- 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
Landscapes
- Radiation Pyrometers (AREA)
Abstract
本发明公开了一种基于深度学习的一维综合孔径微波辐射计SST反演方法,包括有以下步骤:首先:构建用于训练和验证深度学习模型的数据集,为海面温度反演提供数据支撑;其次:构建由自编码器和全连接层耦合的深度学习模型,以训练数据集训练自编码器,使自编码器达到降低数据误差的效果;之后:在自编码器的解码层后连接全连接层,构建完整的深度学习模型,将自编码器的输出作为全连接层的输入,通过最小化损失函数微调整个全连接层,从而反演出海面温度;最后:用验证数据集对整个深度学习模型检验,本发明解决了多入射角情况下,反演海面温度困难的问题,提高了反演的效率。
Description
技术领域
本发明涉及遥感技术领域,具体涉及一种基于深度学习的一维综合孔径微波辐射计SST反演方法。
背景技术:
海面温度(Sea Surface Temperature,SST)在全球气候变化和长期天气过程中起着重要作用,其中被动微波遥感能够进行全天候、全天时不间断地观测。被动微波遥感的代表仪器之一是实孔径微波辐射计,其能够提供包括海面温度在内的多种海洋环境要素产品。但是,由于实孔径微波辐射计的空间分辨率受制于天线的尺寸,导致其空间分辨率较低。针对这一缺点设计了一维综合孔径微波辐射计,与实孔径微波辐射计机械扫描成像方式不同,一维综合孔径微波辐射计采用小孔径天线阵列,解决了传统实孔径微波辐射计空间分辨率与天线物理口径之间的固有矛盾,有效的提高了观测的空间分辨率。但是,由于成像方式的差别,导致已有的实孔径微波辐射计海面温度反演算法无法应用于一维综合孔径微波辐射计。传统实孔径微波辐射计一般是以固定入射角扫描成像,而综合孔径微波辐射计对场景属于凝视成像,其入射角一般是在某一范围内变化,相比于实孔径微波辐射计更加复杂。
由于,目前尚没有用于海面温度(SST)观测的星载一维综合孔径微波辐射计,因此,需要设计一种用于星载一维综合孔径微波辐射计反演海面温度的方法。一维综合孔径微波辐射计在大气顶接收到的亮温是海水盐度、海面温度、海面风速、海面相对风向、大气水汽含量、云液态水含量和入射角的函数,本可以利用微波辐射传输正演模型反演出海面温度。但是,由于海洋大气环境的复杂性,这就增加了利用理论模型反演海面温度的难度。深度学习方法只需要利用实测数据集就可反演出海面温度,而不需要深究微波发射和传输过程中的深层物理机制,具有强大的非线性拟合能力,能够高精度地反演出海面温度,本发明就是提供这样一种基于深度学习的一维综合孔径微波辐射计SST反演方法,为后续星载一维综合孔径微波辐射计的载荷研制及应用提供技术支撑。
发明内容
本发明的目的在于提供一种基于深度学习的一维综合孔径微波辐射计SST反演方法,打破传统实孔径微波辐射计的思路,利用深度学习强大的非线性拟合性,解决了一维综合孔径微波辐射计多入射角反演海面温度困难的问题,能够高效、高精度地反演出海面温度。
一种基于深度学习的一维综合孔径微波辐射计SST反演方法,包括以下步骤:
步骤1:构建用于训练深度学习模型的数据集A,A包括用作模型输出的海面温度集Y,和用作模型输入的垂直极化亮温、水平极化亮温、入射角、海面风速、海面相对风向、大气水汽含量和云液态水含量集X,在X中加入随机噪声得到新的X′作为数据集B,将上述数据集A和B随机分为训练集和验证集两部分,训练集和验证集分别占数据总数的75%和25%;
步骤2:构建由自编码器和全连接层耦合的深度学习模型,将训练集中的X′输入自编码器,将训练集中的X作为自编码器输出的理想值,用最小化负对数似然的损失函数训练自编码器,直至解码层还原的数据与X的误差最小,从而达到降噪的效果;
步骤3:在自编码器的解码层后连接全连接层,构建完整的深度学习模型,将步骤2中自编码器的输出作为全连接层的输入,将训练集中的Y作为全连接层输出的理想值,通过最小化损失函数微调整个全连接层,从而反演出海面温度;
步骤4,将验证集中的X′作为深度学习模型的输入,计算对应的输出,并与验证集中的Y进行误差比较,检验方法的准确度。
优选的,所述步骤2中将训练集中X′输入到自编码器,采用如下算法:
X″=s(w×X′+b)
X″′=s(w′×X″+b′)
其中,X″是X′编码后的特征表达,w代表输入层到隐层的权重,b是隐层偏置,s是激活函数,w′是解码权重,b′是解码偏置,X″′是X″解码后的结果。
优选的,所述步骤2中通过最小化负对数似然的损失函数来训练这个模型,使得X″′和不含随机误差的X差距最小,认为这个模型达到了去除噪音的效果,具体公式为:L=-logP(X|X″′)。
优选的,所述步骤3中将步骤2中自编码器的输出作为全连接层的输入,将训练集中的Y作为全连接层输出的理想值,进行海面温度反演。
优选的,所述步骤3中的最小化损失函数为Levenberg-Marquardt最小化算法,通过大量的训练,得到精确的海面温度。
本发明的优点在于:该种基于深度学习的一维综合孔径微波辐射计SST反演方法,综合孔径微波辐射计与传统实孔径微波辐射计完全不同,其采用小孔径天线阵列,通过场景辐射亮温的空间频率采样得到可见度函数,再进行逆傅里叶变换等数学运算重构出场景的亮温图像,根据这一特点,构建用于模型训练的数据集A和B,B由A中X加入随机误差得到,加入误差之前的数据做真实数据。其次,采用最小化负对数似然的损失函数的方法,用构建好的数据训练自编码器,使其具有降噪的效果;之后,在训练好的自编码器之后连接全连接层,采用最小化损失函数和Levenberg-Marquardt最小化算法的方法,用大量数据训练全连接层。最终使深度学习模型可以解决综合孔径微波辐射计多入射角反演海面温度困难的问题,能够高效、高精度的反演出海面温度。
附图说明
图1为本发明中深度学习模型结构图。
图2为本发明中检验结果图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
如图1至图2所示,一种基于深度学习的一维综合孔径微波辐射计SST反演方法,包括以下步骤:
步骤1:为了得到较为准确的数据分布,从欧洲中期天气预报中心(ECMWF)获取2015年1月1日至12月31日每天的1°×1°海平面模式数据,包括海面温度、海面风速、海面风向、大气水汽含量和云液态水含量等要素。筛选出13836组数据,将数据输入微波辐射传输正演模型中,计算一维综合孔径微波辐射计接收到的垂直极化和水平极化亮温。然后用这些数据构建用于训练深度学习模型的数据集A,A包括用作模型输出的海面温度集Y和用作模型输入的垂直极化亮温、水平极化亮温、入射角、海面风速、海面相对风向、大气水汽含量和云液态水含量集X,在X中加入随机噪声得到新的X′作为数据集B,将上述数据集A和B随机分为训练集和验证集两部分,其中训练集占75%,验证集占25%。我们称X′为观测数据,X为真实数据;
步骤2:构建由自编码器和全连接层耦合的深度学习模型,使用自编码器的作用是为了降低观测数据集X′的误差,使X′与X的误差更小。将训练集中的X′输入自编码器,将训练集中的X作为自编码器输出的理想值,用最小化负对数似然的损失函数训练自编码器,直至解码层还原的数据与X的误差最小,从而达到降噪的效果;
步骤3:在自编码器的解码层后连接全连接层,构建完整的深度学习模型,将步骤2中自编码器的输出作为全连接层的输入,将训练集中的Y作为全连接层输出的理想值,通过最小化损失函数微调整个全连接层,从而反演出海面温度;
步骤4,在前面步骤的前提下,将验证集中的X′作为深度学习模型的输入,计算对应的输出,并与验证集中的Y进行误差比较,检验方法的准确度。
值得注意的是,所述步骤2中将训练集中X′输入到自编码器,采用如下算法:
X″=s(w×X′+b)
X″′=s(w′×X″+b′)
其中,X″是X′编码后的特征表达,w代表输入层到隐层的权重,b是隐层偏置,s是激活函数,w′是解码权重,b′是解码偏置,X″′是X″解码后的结果。
在本实施例中,所述步骤2中通过最小化负对数似然的损失函数来训练这个模型,使得X″′和不含随机误差的X差距最小,认为这个模型达到了去除噪音的效果,具体公式为:L=-logP(X|X″′)。
在本实施例中,所述步骤3中将步骤2中自编码器的输出作为全连接层的输入,将训练集中的Y作为全连接层输出的理想值,进行海面温度反演。
在本实施例中,所述海面温度反演主要包括:使用训练完成的自编码器消除观测数据的噪音后,将训练完成的数据输入到全连接层中,由于反演海面温度用到的数据为:一维综合孔径微波辐射计垂直极化和水平极化观测亮温、海面风速、海面相对风向、大气水含量、云液态水含量和入射角。一共7个气象海洋要素,所以输入层设为7个神经元加1个偏置,隐层层数设为2层,每层神经元个数分别为4个和2个,每层各包括1个偏置(本专利意在讲思路,隐层层数和神经元数需通过试验确定)输出层神经元数为1个(海面温度),各层内神经元之间互不连接,各层间神经元通过不同权值相互连接。然后按照设定好的激活函数,逐层向后传递,通过比较输出数据和输入数据的期望产生误差信号,将误差信号逐层反向传递,并不断修正各层的权值,得到海面温度。其中激活函数为sigmiod函数:
在本实施例中,所述步骤3中的最小化损失函数为Levenberg-Marquardt最小化算法,通过大量的训练,得到精确的海面温度。
由技术常识可知,本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此,上述公开的实施方案,就各方面而言,都只是举例说明,并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。
Claims (6)
1.一种基于深度学习的一维综合孔径微波辐射计海面温度反演方法,其特征在于,包括以下步骤:
步骤1:构建用于训练深度学习模型的数据集A,A包括用作模型输出的海面温度集Y,和用作模型输入的垂直极化亮温、水平极化亮温、入射角、海面风速、海面相对风向、大气水汽含量和云液态水含量集X,在X中加入随机噪声得到新的X′作为数据集B,将上述数据集A和B随机分为训练集和验证集两部分,训练集和验证集分别占数据总数的75%和25%;
步骤2:构建由自编码器和全连接层耦合的深度学习模型,将训练集中的X′输入自编码器,将训练集中的X作为自编码器输出的理想值,用最小化负对数似然的损失函数训练自编码器,直至解码层还原的数据与X的误差最小,从而达到降噪的效果;
步骤3:在自编码器的解码层后连接全连接层,构建完整的深度学习模型,将步骤2中自编码器的输出作为全连接层的输入,将训练集中的Y作为全连接层输出的理想值,通过最小化损失函数微调整个全连接层,从而反演出海面温度;
步骤4,将验证集中的X′作为深度学习模型的输入,计算对应的输出,并与验证集中的Y进行误差比较,检验方法的准确度。
2.根据权利要求1所述的一种基于深度学习的一维综合孔径微波辐射计海面温度反演方法,其特征在于:所述步骤2中将训练集中X′输入到自编码器,采用如下算法:
X″=s(w×X′+b)
X″′=s(w′×X″+b′)
其中,X″是X′编码后的特征表达,w代表输入层到隐层的权重,b是隐层偏置,s是激活函数,w′是解码权重,b′是解码偏置,X″′是X″解码后的结果。
3.根据权利要求2所述的一种基于深度学习的一维综合孔径微波辐射计海面温度反演方法,其特征在于:所述步骤2中通过最小化负对数似然的损失函数来训练这个模型,使得X″′和不含随机误差的X差距最小,认为这个模型达到了去除噪音的效果,具体公式为:L=-log P(X|X″′)。
4.根据权利要求1所述的一种基于深度学习的一维综合孔径微波辐射计海面温度反演方法,其特征在于:所述步骤3中将步骤2中自编码器的输出作为全连接层的输入,将训练集中的Y作为全连接层输出的理想值,进行海面温度反演。
6.根据权利要求5所述的一种基于深度学习的一维综合孔径微波辐射计海面温度反演方法,其特征在于:所述步骤3中的最小化损失函数为Levenberg-Marquardt最小化算法,通过大量的训练,得到精确的海面温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811547210.7A CN109829547B (zh) | 2018-12-18 | 2018-12-18 | 一种基于深度学习的一维综合孔径微波辐射计sst反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811547210.7A CN109829547B (zh) | 2018-12-18 | 2018-12-18 | 一种基于深度学习的一维综合孔径微波辐射计sst反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109829547A CN109829547A (zh) | 2019-05-31 |
CN109829547B true CN109829547B (zh) | 2020-10-09 |
Family
ID=66859812
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811547210.7A Active CN109829547B (zh) | 2018-12-18 | 2018-12-18 | 一种基于深度学习的一维综合孔径微波辐射计sst反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109829547B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114184999B (zh) * | 2020-09-15 | 2024-05-14 | 复旦大学 | 一种互耦小孔径阵列的生成式模型处理方法 |
CN112180369B (zh) * | 2020-09-23 | 2021-09-10 | 中国人民解放军国防科技大学 | 基于深度学习的一维综合孔径辐射计海面风速反演方法 |
CN113408742B (zh) * | 2021-06-24 | 2023-06-02 | 桂林理工大学 | 一种基于机器学习的高精度海表温度反演方法 |
CN114265087A (zh) * | 2021-12-24 | 2022-04-01 | 青岛镭测创芯科技有限公司 | 一种海浪参数反演方法及相关组件 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1301969A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 顺轨方向多基线干涉式综合孔径微波辐射计及其设计方法 |
CN201600451U (zh) * | 2009-09-28 | 2010-10-06 | 中国科学院空间科学与应用研究中心 | 采用干涉式孔径综合技术的大气温度探测仪的外定标装置 |
CN105760699A (zh) * | 2016-03-18 | 2016-07-13 | 中国科学院国家空间科学中心 | 一种海表盐度反演方法和装置 |
CN105988146A (zh) * | 2015-01-29 | 2016-10-05 | 中国科学院空间科学与应用研究中心 | 一种星载微波辐射计的应用数据处理方法 |
CN107300561A (zh) * | 2016-04-15 | 2017-10-27 | 北京空间飞行器总体设计部 | 基于多遥感器联合探测的海洋盐度卫星 |
CN108875905A (zh) * | 2018-04-09 | 2018-11-23 | 华中科技大学 | 一种大气温湿廓线的可见度函数直接反演方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060064279A1 (en) * | 2004-09-17 | 2006-03-23 | Eric Baum | Algorithm for retrieval of ocean surface temperature, wind speed and wind direction from remote microwave radiometric measurements |
-
2018
- 2018-12-18 CN CN201811547210.7A patent/CN109829547B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1301969A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 顺轨方向多基线干涉式综合孔径微波辐射计及其设计方法 |
CN201600451U (zh) * | 2009-09-28 | 2010-10-06 | 中国科学院空间科学与应用研究中心 | 采用干涉式孔径综合技术的大气温度探测仪的外定标装置 |
CN105988146A (zh) * | 2015-01-29 | 2016-10-05 | 中国科学院空间科学与应用研究中心 | 一种星载微波辐射计的应用数据处理方法 |
CN105760699A (zh) * | 2016-03-18 | 2016-07-13 | 中国科学院国家空间科学中心 | 一种海表盐度反演方法和装置 |
CN107300561A (zh) * | 2016-04-15 | 2017-10-27 | 北京空间飞行器总体设计部 | 基于多遥感器联合探测的海洋盐度卫星 |
CN108875905A (zh) * | 2018-04-09 | 2018-11-23 | 华中科技大学 | 一种大气温湿廓线的可见度函数直接反演方法 |
Non-Patent Citations (1)
Title |
---|
星载一维综合孔径微波辐射计海面温度反演研究;陈冠宇等;《第35届中国气象学会年会 S21 卫星气象与生态遥感》;20181024;摘要 * |
Also Published As
Publication number | Publication date |
---|---|
CN109829547A (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109829547B (zh) | 一种基于深度学习的一维综合孔径微波辐射计sst反演方法 | |
CN109725316B (zh) | 一种基于一维综合孔径微波辐射计海面温度物理反演方法 | |
CN109725317B (zh) | 一种基于一维综合孔径微波辐射计的海面亮温成像仿真方法 | |
CN113534159B (zh) | 一种基于sar卫星遥感数据的海面风速反演方法及系统 | |
CN104820250B (zh) | 一种极轨气象卫星virr海洋上云检测的处理方法 | |
CN103869296B (zh) | 一种基于成像面表征的极化sar地形辐射校正和几何纠正方法 | |
CN112180369B (zh) | 基于深度学习的一维综合孔径辐射计海面风速反演方法 | |
CN109655828B (zh) | 一种多频一维综合孔径微波辐射计sst深度学习反演方法 | |
CN110793649B (zh) | 一种使用后向散射截面校正粗糙海面亮温的方法 | |
CN104391295A (zh) | 一种图像熵最优的压缩传感sar稀疏自聚焦成像方法 | |
CN108646247B (zh) | 基于伽马过程线性回归的逆合成孔径雷达成像方法 | |
CN113050091B (zh) | 一种星载合成孔径雷达风速风向联合反演方法及系统 | |
CN105891832A (zh) | 一种基于交叉极化模型与cmod5n的海面风速反演方法 | |
CN114117886A (zh) | 一种用于多光谱遥感水深反演方法 | |
CN105353375A (zh) | 一种镜像综合孔径辐射计反演方法 | |
CN110991087B (zh) | 一种基于多入射角组网sar卫星数据的风场反演方法及系统 | |
CN110376572B (zh) | 一种机载激光雷达水下探测深度模拟装置 | |
Shangguan et al. | GF-3 polarimetric data quality assessment based on automatic extraction of distributed targets | |
CN101697009A (zh) | 海面波面还原方法 | |
Zhang et al. | A multiband model with successive projections algorithm for bathymetry estimation based on remotely sensed hyperspectral data in Qinghai Lake | |
CN115097451A (zh) | 一种基于sar卫星数据的海浪多参数反演方法及系统 | |
US11971309B2 (en) | Measurement data processing device | |
CN102073038A (zh) | 基于微小地形的遥感影像的地形校正的方法 | |
CN117237430B (zh) | 一种高精度多时序水深反演方法、计算设备及存储介质 | |
CN111090102B (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 |