一种基于Sentinel数据与深度学习模型的农田土壤水分模拟方法与流程

一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法
技术领域
1.本发明涉及微波遥感技术领域，尤其涉及一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法。


背景技术：

2.经检索，中国专利号cn105510173b公开了一种现场测定土壤含水量的方法，该发明方法通过称重法测定土壤含水量，其虽然在测定过程中可以有效避免土壤水分损失，但这类方法只能获取有限点位的土壤含水量信息，费时费力、成本较高，并且监测精度受到样点密度和空间分布的影响；土壤水分，是地表生态系统的重要组成部分，是农业、水文、气象等科学领域中的一个重要研究参数，其决定着农作物的生长状况；sentinel系列卫星，是欧盟委员会和欧洲航天局共同倡议的全球环境与安全监测系统的重要组成部分，其主要用于全球环境监测，最早发射的两颗卫星sentinel
‑
1a和sentinel
‑
1b分别于2014年4月3日和2016年4月25日升空，搭载了c波段合成孔径雷达，主要完成雷达成像任务；其中sentinel
‑
1a于2014年10月开始面向全球用户提供免费数据，可提供重复观测的c波段合成孔径雷达影像，双星星座的重访周期为6d，空间分辨率最优为5m，如此高空间分辨率和高时间分辨率的数据为区域尺度土壤水分的实时监测提供了良好的数据来源；而目前，传统的土壤水分监测手段主要是运用烘干称重法、中子水分仪法、时域反射仪等进行野外实地测量，但这类方法只能获取有限点位的土壤含水量信息，费时费力、成本较高，并且监测精度受到样点密度和空间分布的影响，而sentinel系列卫星具有大面积实时观测、效率高、成本低等优势；因此，发明出一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法变得尤为重要；
3.现有的农田土壤水分模拟方法大多通过传统的土壤水分监测手段实现的，该类方法只能获取有限点位的土壤含水量信息，费时费力、成本较高，并且监测精度受到样点密度和空间分布的影响，无法全天时、全天候的对农田土壤水分进行模拟，从而不利于实时针对农业生产进行保护；为此，我们提出一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法。


技术实现要素：

4.本发明的目的是为了解决现有技术中存在的缺陷，而提出的一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法。
5.为了实现上述目的，本发明采用了如下技术方案：
6.一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法，该模拟方法具体步骤如下：
7.步骤一、sentinel数据采集；利用sentinel
‑
1搭载的合成孔径雷达对待研究目标区域进行遥感图像拍摄，得到sentinel数据；
8.步骤二、数据预处理：通过专用处理软件对步骤一所述sentinel数据进行预处理；
9.步骤三、土壤水分分析：提取经过预处理后所述sentinel数据中的雷达后向散射系数，并基于alpha模型进行分析，同时根据介电模型求解土壤体积含水量；
10.步骤四、土壤水分模拟：获取步骤三所述雷达后向散射系数和土壤体积含水量，并基于人工神经网络进行深度学习，构建预测模型，最后基于预测模型并通过3d技术对所述待研究目标区域进行土壤水分变化模拟。
11.进一步地，步骤一所述sentinel数据包括2个时相t1和t2的sar影像。
12.进一步地，步骤二所述预处理包括但不限于辐射校正、地形辐射校正、滤波和地形校正；所述专用处理软件具体为snap或sentinel
‑
1toolbox中的一种。
13.进一步地，步骤三所述土壤水分分析过程如下：
14.s1：将2个时相的后向散射系数做比值处理，得到alpha模型，其具体公式如下：
[0015][0016]
式中：σ0表示雷达后向散射系数；θ为雷达入射角度；ε
s
为土壤相对介电常数；t1和t2表示雷达数据获取的时间；α
pp
为极化幅度，是雷达入射角度和土壤介电常数的函数；pp表示极化方式，为hh或者vv；所述极化幅度α
pp
可表达为如下公式：
[0017][0018][0019]
s2：基于alpha模型构建观测方程，其具体公式如下：
[0020][0021]
s3：获取n个时间序列的sar影像，并基于式(1)构建矩阵方程组，其具体公式如下：
[0022][0023]
s4：采用边界约束最小二乘对步骤s3所述方程组进行求解，得到极化幅度α
pp
的值，并基于式(2)和式(3)求得土壤的介电常数ε
s
；
[0024]
s5：最后，采用介电模型并根据步骤s4所述土壤的介电常数求解土壤体积含水量。
[0025]
进一步地，所述介电模型的具体公式如下：
[0026]
当频率在1.4ghz
‑
1.8ghz时，所述介电常数的模型表达式为：
[0027][0028]
当频率在1.4ghz
‑
1.8ghz时，所述介电常数的模型表达式为：
[0029]
ε＇
m
＝1.15ε
m
‑
0.68
ꢀꢀꢀ
(7)
[0030]
式中：p
b
为土壤容重；p
s
为土壤比重；m
v
为土壤水分；β为指数；ε
fw
为自由水介电常数。
[0031]
进一步地，步骤四所述预测模型是基于雷达后向散射系数和土壤体积含水量形成的回归模型，其通过深度学习技术进行构建，其具体构建过程如下：
[0032]
ss1：首先，提取大量sentinel数据中的雷达后向散射系数以及经过土壤水分分析后得到的土壤体积含水量；
[0033]
ss2：然后，对所述雷达后向散射系数和土壤体积含水量进行向量化处理，形成样本集，并将其划分为70％的训练集和30％的测试集；
[0034]
ss3：接着，构建基于人工神经网络的回归模型，并将所述70％的训练集作为输入数据输入其中进行训练，即得到预测模型；
[0035]
ss4：最后，将所述30％的测试集输入进行步骤ss3所述预测模型，若该模型符合预期，则输出该模型，否则，进行重采样，直至达到期望值。
[0036]
进一步地，所述雷达后向散射系数包括但不限于土壤介电常数、地表粗糙度和植被覆盖因素。
[0037]
相比于现有技术，本发明的有益效果在于：
[0038]
本技术公开了一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法，其利用sentinel
‑
1卫星获取的多时相c波段合成孔径雷达数据，并基于深度学习模型进行回归学习，通过研究雷达后向散射系数与土壤体积含水量之间的关系，从而实现了对农田土壤水分的模拟，并且本发明考虑到了雷达后向散射系数中土壤介电常数、地表粗糙度和植被覆盖等因素的影响，通过alpha模型进行消除，从而保证了对农田土壤水分计算精度，其相较于现有的农田土壤水分模拟方法具有大面积实时观测、效率高、成本低等优势，能够全天时、全天候的对农田土壤水分进行模拟，从而有利于实时针对农业生产进行保护。
附图说明
[0039]
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。
[0040]
图1为本发明提出的一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法的整体流程图。
具体实施方式
[0041]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。
[0042]
在本发明的描述中，需要理解的是，术语“上”、“下”、“前”、“后”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以
特定的方位构造和操作，因此不能理解为对本发明的限制。
[0043]
参照图1，本实施例公开了一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法，该模拟方法具体步骤如下：
[0044]
步骤一、sentinel数据采集；利用sentinel
‑
1搭载的合成孔径雷达对待研究目标区域进行遥感图像拍摄，得到sentinel数据；
[0045]
具体的，该sentinel数据包括2个时相t1和t2的sar影像。
[0046]
步骤二、数据预处理：通过专用处理软件对步骤一sentinel数据进行预处理；
[0047]
具体的，该预处理包括但不限于辐射校正、地形辐射校正、滤波和地形校正；专用处理软件具体为snap或sentinel
‑
1toolbox中的一种。
[0048]
步骤三、土壤水分分析：提取经过预处理后sentinel数据中的雷达后向散射系数，并基于alpha模型进行分析，同时根据介电模型求解土壤体积含水量；
[0049]
具体的，该土壤水分分析过程如下：
[0050]
s1：将2个时相的后向散射系数做比值处理，得到alpha模型，其具体公式如下：
[0051][0052]
式中：σ0表示雷达后向散射系数；θ为雷达入射角度；ε
s
为土壤相对介电常数；t1和t2表示雷达数据获取的时间；α
pp
为极化幅度，是雷达入射角度和土壤介电常数的函数；pp表示极化方式，为hh或者vv；极化幅度α
pp
可表达为如下公式：
[0053][0054][0055]
s2：基于alpha模型构建观测方程，其具体公式如下：
[0056][0057]
s3：获取n个时间序列的sar影像，并基于式(1)构建矩阵方程组，其具体公式如下：
[0058][0059]
s4：采用边界约束最小二乘对步骤s3方程组进行求解，得到极化幅度α
pp
的值，并基于式(2)和式(3)求得土壤的介电常数ε
s
；
[0060]
s5：最后，采用介电模型并根据步骤s4土壤的介电常数求解土壤体积含水量；
[0061]
具体的，该介电模型的具体公式如下：
[0062]
当频率在1.4ghz
‑
1.8ghz时，介电常数的模型表达式为：
[0063][0064]
当频率在1.4ghz
‑
1.8ghz时，介电常数的模型表达式为：
[0065]
ε＇
m
＝1.15ε
m
‑
0.68
ꢀꢀꢀ
(7)
[0066]
式中：p
b
为土壤容重；p
s
为土壤比重；m
v
为土壤水分；β为指数；ε
fw
为自由水介电常数。
[0067]
步骤四、土壤水分模拟：获取步骤三雷达后向散射系数和土壤体积含水量，并基于人工神经网络进行深度学习，构建预测模型，最后基于预测模型并通过3d技术对待研究目标区域进行土壤水分变化模拟；
[0068]
具体的，该雷达后向散射系数包括但不限于土壤介电常数、地表粗糙度和植被覆盖因素。
[0069]
参照图1，本实施例公开了一种基于sentinel数据与深度学习模型的农田土壤水分模拟方法，除与上述实施例相同的结构外，本实施例将具体介绍预测模型；
[0070]
具体的，该预测模型是基于雷达后向散射系数和土壤体积含水量形成的回归模型，其通过深度学习技术进行构建，其具体构建过程如下：首先，提取大量sentinel数据中的雷达后向散射系数以及经过土壤水分分析后得到的土壤体积含水量；然后，对雷达后向散射系数和土壤体积含水量进行向量化处理，形成样本集，并将其划分为70％的训练集和30％的测试集；接着，构建基于人工神经网络的回归模型，并将70％的训练集作为输入数据输入其中进行训练，即得到预测模型；最后，将30％的测试集输入进行步骤ss3预测模型，若该模型符合预期，则输出该模型，否则，进行重采样，直至达到期望值。
[0071]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。