本发明属于水文遥感技术领域，具体涉及一种基于多源遥感数据融合的土壤水分测量方法及装置。
背景技术：
土壤水分作为陆面资源形成、转化、消耗过程研究中的基本参数，是土壤-植被-大气系统能量与物质交换的重要参数之一，也是地表能量平衡的重要决定因素。
常规的土壤水分地面测量方法，如重量法测量和时域反射仪法都是基于点的测量，仅能代表空间范围很小(0.0025m2)的土壤水分信息，采样点有限且成本高。
随着遥感技术的发展，光学和微波遥感技术成为估算土壤水分的重要手段。常用的光学遥感手段，如可见光和近红外波段可通过识别土壤颜色和地表温度等因素反演土壤水分，但易受云雾的影响。微波遥感，如被动微波遥感反演土壤水分的机理在于地表反射的微波后向散射系数与土壤介电常数密切相关，而土壤的介电常数主要由土壤水分决定，但其空间分辨率较低。
全球卫星导航系统(globalnavigationsatellitesystem,gnss)除向全球空间用户提供导航定位功能外，其反射信号(gnssreflectometry,gnss-r)也可用于地表参数监测。相比于其他遥感手段，利用gnss-r遥感数据估算土壤水分具有低成本、高时空分辨率、覆盖范围广、精度高等诸多优点。随着星载旋风全球导航卫星系统(cygnss)卫星的成功发射，可为大范围、高空间分辨率的土壤水分研究提供更多数据源。但星载gnss卫星数据为不连续分布的离散点，且密度较低。
因此，无论是常规的土壤水分地面测量方法还是目前的遥感技术的估算方法，得到的土壤水分的结果与实际均存在较差的差距。
因此，需要提供一种针对上述现有技术中不足的改进技术方案。
技术实现要素：
本发明的目的是提供一种基于多源遥感数据融合的土壤水分测量方法及装置，用以克服上述现有技术中存在土壤水分测量结果差的问题。
为了实现上述目的，本发明提供如下技术方案：
一种基于多源遥感数据融合的土壤水分测量方法，用于对目标区域土壤所含的水分进行测量，其特征在于，
通过旋风全球导航卫星系统获取所述目标区域土壤的原始数据，通过所述原始数据和所述目标区域土壤的实测数据确定所述目标区域的土壤水分估算结果；
通过中分辨率成像光谱分析仪获取所述目标区域土壤的第一地表温度数据和通过中国气象局陆面数据同化系统获取所述目标区域土壤的第二地表温度数据，将所述第一地表温度数据和所述第二地表温度数据进行时空融合，得到地表温度估算结果；
获取影响土壤水分的辅助参数；所述辅助参数包括植被覆盖指数、地面高程、植被覆盖度、土壤粗糙度、地面坡度、降水；
将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行数据融合，得到多源遥感数据融合模型，通过所述多源遥感数据融合模型得到土壤水分测量结果。
进一步的，所述土壤水分估算结果的确定过程为：
通过cygnss获取所述目标区域土壤的原始数据，并对所述原始数据进行预处理，得到标准数据；
利用双基雷达方程对所述标准数据进行处理，得到地表有效反射率；
获取所述目标区域土壤的实测数据，通过所述实测数据和所述地表有效反射率确定所述土壤水分估算结果。
进一步的，所述预处理为采用数据质量控制方法对所述原始数据进行处理。
进一步的，所述地表有效反射率为：
其中，sr为地表有效反射率；snr是峰值功率减去噪声后的信噪比值；prt是发射功率；gt是gps发射天线的增益；gr是镜面反射点的天线增益；λ表示gpsl1频段信号的波长；rts是从gnss卫星到镜面反射点的范围；rsr是从镜面反射点到接收器的距离；θ为入射角。
进一步的，所述地表温度估算结果的确定过程为：
通过光谱分析仪获取所述目标区域土壤的第一地表温度数据和通过中国气象局陆面数据同化系统获取所述目标区域土壤的第二地表温度数据；
将所述第一地表温度数据和所述第二地表温度数据进行时空融合处理，得到所述地表温度估算结果。
进一步的，所述时空模型是通过像元反射率的时间变化特征的改进型自适应遥感图像时空融合方法构建得到的，所述时空模型为estarfm时空模型。
进一步的，通过所述时空模型进行的融合处理包括邻近相似像元搜索、权重值计算、计算转换系数。
进一步的，所述辅助参数通过权重对比方法得到的。
进一步的，执行将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行数据融合，得到多源遥感数据融合模型，通过所述多源遥感数据融合模型得到土壤水分测量结果的过程如下：
将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行空间匹配；
将进行空间匹配后的土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行处理，得到多源遥感数据融合模型；
通过所述多源遥感数据融合模型得到土壤水分测量结果。
为了实现上述目的，本发明提供一种基于多源遥感数据融合的土壤水分测量装置，包括处理器、存储器以及存储在所述存储器中并可在所述处理器上的运行的计算机程序，所述处理器在执行所述计算机程序时实现如上述的基于多源遥感数据融合的土壤水分测量方法的步骤。
与最接近的现有技术相比，本发明提供的技术方案具有如下优异效果：
本发明基于多源遥感数据融合的土壤水分测量方法在土壤水分测量过程中，基于bp神经网络，融合星载cygnss数据和光学modis数据，即利用空间上广泛分布、时间上连续运行的星载cygnss数据作为主要数据源进行地表土壤水分反演，探测结果同样可以为利用l波段的土壤水分遥感卫星提供标定；最终获得高质量、高时空分辨率、且空间连续分布的土壤水分数据。通过本申请的多源遥感数据模型能够获得更加准确的土壤水分测量结果。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。其中：
图1为本发明多源遥感数据融合的土壤水分测量方法的总体流程图；
图2为本发明多源遥感数据融合的土壤水分测量方法的详细流程图；
图3为本发明cygnss卫星信号传播和散射示意图；
图4为本发明实施例中利用bp神经网络融合cygnss、modis及辅助数据的示意图；
图5为本发明实施例中2018年4月1～7日估算土壤水分空间分布图；
图6为本发明实施例中估算土壤水分与实测土壤水分线性拟合的示例。
具体实施方式
下面将参考附图并结合实施例来详细说明本发明。各个示例通过本发明的解释的方式提供而非限制本发明。实际上，本领域的技术人员将清楚，在不脱离本发明的范围或精神的情况下，可在本发明中进行修改和变型。例如，示为或描述为一个实施例的一部分的特征可用于另一个实施例，以产生又一个实施例。因此，所期望的是，本发明包含归入所附权利要求及其等同物的范围内的此类修改和变型。
本发明基于多源遥感数据融合的土壤水分测量方法，主要用于对一定区域土壤所含的水分进行测量，主要技术构思如图1所示，将cygnss土壤水分估算结果、时空融合方式得到的地表温度估算结果和辅助参数这三种类型的数据通过神经网络实现多源数据融合，通过多源数据融合的处理得到土壤水分估算结果，本发明通过如图2所示基于多源遥感数据融合的土壤水分测量方法得到的土壤水分估算结果更加准确。
方法实施例：
步骤s1，通过星载旋风全球导航卫星系统cygnss获取所述目标区域土壤的原始数据，通过所述原始数据和所述目标区域土壤的实测数据确定所述目标区域的土壤水分估算结果。
在执行通过星载旋风全球导航卫星系统cygnss获取原始数据，通过所述原始数据和所述土壤的实测水分确定土壤水分估算结果时，包括如下过程：
步骤s11，通过cygnss获取所述目标区域土壤的原始数据，并对所述原始数据进行预处理，得到标准数据；
通过星载旋风全球导航卫星系统cygnss获取目标区域土壤的原始数据，所述原始数据包括信噪比值(snr)、发射功率(prt)、gps发射天线的增益(gt)、镜面反射点的天线增益(gr)、gpsl1频段信号的波长(λ)、从gnss卫星到镜面反射点的范围(rts)、从镜面反射点到接收器的距离(rsr)、入射角(θ)。
如图3所示为cygnss卫星信号传播和散射示意图，通过星载旋风全球导航卫星系统(cygnss)获取的原始数据受系统噪声、天线增益异常、等效全向辐射功率异常等问题影响较大，若直接使用原始数据进行估算，会导致土壤水分的估算结果准确性差；另外，当原始数据的地面高程大于600米时，使用数据压缩法也会原始数据的质量产生影响，因此，在本申请中，为了提高通过cygnss获取的原始数据的质量，采用数据质量控制方法对原始数据进行预处理，将预处理后的原始数据称为标准数据；
步骤s12，利用双基雷达方程对所述标准数据进行处理，得到地表有效反射率；
在土壤水分计算中，假设地表反射信号是一个连续的观测量，因此，多基于星载反射信号gnss-r估算土壤水分研究，将研究目标区域集中在植被覆盖度较小或地表粗糙度较低的区域。利用双基雷达方程可知：
式中，pr为特定极化方式下经地面反射后到达接收机的反射信号功率，可视为snrp；pt为gps卫星发射的右旋圆极化(rhcp)功率；gt为gps发射天线的增益；rts为gps卫星与镜面反射点之间的距离；rsr为镜面反射点与镜面反射点之间的距离；gr为接收天线的增益；λ是gps波长；γp是特定极化方式下的地表反射率。
假如忽略gnss卫星发射功率的波动及不同gnss卫星的反射功率差异，将pt与gt视为常量，那么主要影响pr的分量可以概括为rsr、rts、以及gr。此外，延迟多普勒图(ddm)同时受到系统基底噪声的影响，利用非相干累加的方法在一定程度消弱了噪声的影响，但不能完全消除。因此，需综合考虑反射信号功率、发射机及接收机到镜面反射点的距离，以及接收天线增益对ddm峰值功率的影响。
在本申请实施例中，利用双基雷达方程对信噪比值(snr)、发射功率(prt)、gps发射天线的增益(gt)、镜面反射点的天线增益(gr)、gpsl1频段信号的波长(λ)、从gnss卫星到镜面反射点的范围(rts)、从镜面反射点到接收器的距离(rsr)、入射角(θ)的标准数据进行处理，得到地表有效反射效率。
在本申请实施例中，得到的地表有效反射效率的模型为：
其中，sr为地表有效反射率；snr是峰值功率减去噪声后的信噪比值；prt是发射功率；gt是gps发射天线的增益；gr是镜面反射点的天线增益；λ表示gpsl1频段信号的波长；rts是从gnss卫星到镜面反射点的范围；rsr是从镜面反射点到接收器的距离；θ为入射角。
步骤s13，获取所述目标区域土壤的实测数据，通过所述实测数据和所述地表有效反射率确定土壤水分估算结果；
步骤s131，构建土壤水分估算模型；
在本申请实施例中，土壤水分估算模型为：
smcygnss＝(sr-a)/b
式中，smcygnss为估算的土壤水分，sr为地表有效反射率，a和b分别为通过实测数据与sr线性拟合，得到的第一相关参数和第二相关参数。
步骤s132，通过所述实测数据、所述地表有效反射率和所述土壤水分估算模型，确定土壤水分估算结果；
所述土壤水分估算结果为cygnss土壤水分估算结果，通过土壤水分估算模型进行估算，得到的cygnss土壤水分估算结果精度高，时间分辨率为每日，但结果呈现离散分布的情况。
步骤s2，通过中分辨率成像光谱分析仪获取所述目标区域土壤的第一地表温度数据和通过中国气象局陆面数据同化系统获取所述目标区域土壤的第二地表温度数据，将所述第一地表温度数据和所述第二地表温度数据进行时空融合，得到融合的地表温度数据估算结果。
步骤s21，通过光谱分析仪获取所述目标区域土壤的第一地表温度数据和通过中国气象局陆面数据同化系统获取所述目标区域土壤的第二地表温度数据；
在本申请实施例中，通过光谱分析仪获取所述目标区域土壤的第一地表温度数据时，具体是通过中分辨率成像光谱仪(modis)实现的，那么通过中分辨率成像光谱仪获取的地表温度数据为modislst数据；在本申请实施例中，modislst数据为所述第一地表温度数据。
在本申请实施例中，采用中分辨率成像光谱仪即可满足要求，作为其他实施方式，光谱分析仪也可以采用高分辨率成像光谱分析仪，成本相对较高。可以理解的是，以上描述仅为示例性的，本申请实施例对此并不进行限定。
在本申请实施例中，中国气象局陆面数据同化系统存储有所有陆面的数据信息，通过中国气象局陆面数据同化系统能够直接获取所述目标区域土壤的地表温度数据，即cldaslst数据；在本申请实施例中，cldaslst数据为第二地表温度数据。
步骤s22，将所述第一地表温度数据和所述第二地表温度数据进行时空融合处理，得到所述融合的地表温度数据估算结果；
在本申请实施例中，在构建时空融合模型时考虑像元反射率的时间变化特征的改进型自适应遥感图像时空融合方法所构建的时空模型，即构建的是estarfm(enhancedspatialandtemporaladaptivereflectancefusionmodel)时空融合模型。时空融合模型的构建过程为现有技术，本申请实施例中不再详细介绍。
在本申请实施例中，利用estarfm时空融合模型对cldaslst数据和modislst数据进行时空融合处理包括：选取两对时间相同且无云的modislst数据和cldaslst数据，及预测目标日期的cldaslst数据影像，通过estarfm模型进行邻近相似像元搜索、权重值计算、计算转换系数、预测期中心像元值计算等步骤，得到逐日的融合后的lst影像。
对于处理的过程多为现有技术，以下本申请主要介绍一下modislst数据和cldaslst数据之间的加权函数；
其中，l表示modis影像的反射率；m表示cldas影像的反射率；(xw/2，yw/2)为预测日期中心像元坐标，tp为获取日期，b为影像波段，w表示移动窗口大小，窗口大小由目标区域地表的均匀性决定；n表示相似像元的个数；(xi，yi)为第i个相似像元的坐标；wi代表第i个像元的权重值；vi为计算得到的第i个像元的转换系数。
具体过程如下：
s231，分别从第一地表温度数据中和从第二地表温度数据中选取相同时间对应的modislst数据和cldaslst数据；
在本申请实施例中，以天为单位，要测量的是目标区域一段时间内的土壤水分测量结果，在执行该步骤时，分别先输入的是第一天的第一地表温度数据(即modislst数据)和第一天的第二地表温度数据(即cldaslst数据)。
s232，输入目标区域预测日期的cldaslst数据影像；
s233，基于所述estarfm时空融合模型对所相同时间对应的modislst数据和cldaslst数据、所述cldaslst数据影像进行处理，得到一段时间融合后的lst影像，所述一段时间融合后的lst影像为所述融合的地表温度数据估算结果；
所述处理包括：进行邻近相似像元搜索、权重值计算、计算转换系数、预测期中心像元值计算等。通过estarfm模型对modislst数据和cldaslst数据进行邻近相似像元搜索、权重值计算、计算转换系数、预测期中心像元值计算等步骤，分别得到每一天的lst影像，汇集形成融合后的所述融合的地表温度数据估算结果。
在本申请实施例中，使用两对在同一天获取的modis影像(云雾覆盖度<15％)与cldas影像，以及预测日期的cldas影像来获取高空间分辨率、受云雾影响小的影像。
本申请采用estarfm时空数据融合模型对modis和cldas进行融合，获取融合后的高分辨率、不受云雾影响的，且时空连续的地表温度lst数据，并通过实际同期的modislst评价融合结果的可靠性，在此基础上求得完整的时间序列lst数据。
也就是说，在本申请实施例中，通过采用estarfm时空数据融合模型对modislst数据和cldaslst数据进行融合，能够获得空间完整且分布率较高的lst分布。为了保证融合的地表温度数据结果的lst分布的准确性，再采用实际同期的modislst数据对融合后的地表温度估算结果进行评价，从而保证融合结果的可靠性，即地表温度估算结果具有完整的时间序列lst数据。
步骤s3，获取影响土壤水分的辅助参数；所述辅助参数包括植被覆盖指数(ndvi)、地面高程(dem)、植被覆盖度、土壤粗糙度、地面坡度、和降水；
对于土壤水分的含量受很多因素的影响，例如：地面坡度陡的相对于地面平缓的来讲，地面坡度陡的不易存储水源，水会随着坡度向下流走；降水量大的区域与降水量小的区域对比，降水量小的区域，其土壤水分含量较少；相类似的，植被覆盖度密集的区域，其土壤水分含量较多。
在本申请实施例中，从许多影响土壤水分的影响因素中，采用权重对比方法，对不同辅助参数的结果进行精度评估和筛选，获取对土壤水具有较强相关性的参数组合。综合考虑，选择了包括植被覆盖指数(ndvi)，地面高程(dem)，植被覆盖度、土壤粗糙度、地面坡度、降水等辅助参数，通过选取的这些辅助参数，能够在融合过程中提高土壤水分的估算结果。
在本申请实施例中，辅助参数中，植被覆盖指数(ndvi)、地面高程(dem)、土壤粗糙度、地面坡度和降水是通过数据官方网站免费获取；植被覆盖度和地面坡度通过已有成熟的经验模型获得。
在本申请实施例中，步骤s1、步骤s2与步骤s3之间并不存在先后时序关系，可以先执行步骤s1，再执行步骤s2，后执行步骤s3；也可以先执行步骤2，再执行步骤s1，后执行步骤s3；或者同时执行步骤s1、步骤s2和步骤s3。可以理解的是，以上描述仅为示例性的，本申请实施例对此并不进行限定。
步骤s4，将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行数据融合，得到多源遥感数据融合模型，通过所述多源遥感数据融合模型能够得到土壤水分测量结果。
如图4所示为对应的利用bp神经网络融合cygnss、modis及辅助数据的示意图，具体过程如下：
步骤s41，将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行空间匹配；
通过空间经纬度统一的方式，先将所述土壤水分估算结果、所述地表温度数据估算结果和所述辅助参数的空间进行匹配，便于所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数在统一空间下进行处理。
步骤s42，将进行空间匹配后的土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行处理，得到多源遥感数据融合模型；
在本申请实施例中，通过引入bp神经网络算法，将土壤水分估算结果、所述地表温度估算结果和所述辅助参数作为bp神经网络的输入因子，通过反复实验确定bp神经网络模型的隐含层节点与层数，确定激励函数、学习率等参数；进而得到bp神经网络模型。该bp神经网络模型为多源遥感数据融合模型。
在申请实施例中，bp神经网络可用类似的机器学习方法替代，例如随机森林等。可以理解的是，以上描述仅为示例性的，本申请实施例对此并不进行限定。
在本申请实施例中，bp神经网络采用非线性拟合，隐含层节点为8，层数为2，激励函数为s型函数，最大训练次数设置为6000次，误差度量最小化为rmse，误差阈值设置为0.001，学习率设置为0.05。
步骤s43，通过所述多源遥感数据融合模型能够得到土壤水分测量结果。
在本申请实施例中，将多源遥感数据代入至所述多源遥感数据融合模型，就能得到对应土壤水分的测量结果，高空间分辨率为1km*1km。如图5所示为本发明实施例中2018年4月1～7日估算土壤水分空间分布图；
在本申请实施例中，已应用该模型得到中国东南部地区土壤水分的逐日连续分布图，该分布图与实测结果吻合度较高，如图6所示为估算土壤水分与实测土壤水分线性拟合的结果。因此，通过本申请的多源遥感数据模型能够获得更加准确的土壤水分测量结果。
装置实施例：
本申请还提出一种基于多源遥感数据融合的土壤水分测量装置，包括处理器、存储器以及存储在所述存储器中并可在所述处理器上的运行的计算机程序，其特征在于，所述处理器在执行所述计算机程序时实现上述的基于多源遥感数据融合的土壤水分测量方法的步骤。
由于基于多源遥感数据融合的土壤水分测量方法的步骤已在方法实施例中详细介绍，此处不再进行赘述。
与现有技术相比，本发明基于多源遥感数据融合的土壤水分测量方法在土壤水分测量过程中，基于bp神经网络，融合星载cygnss数据和光学modis数据，即利用空间上广泛分布、时间上连续运行的星载cygnss数据作为主要数据源进行地表土壤水分反演，探测结果同样可以为利用l波段的土壤水分遥感卫星提供标定；最终获得高质量、高时空分辨率、且空间连续分布的土壤水分数据。通过本申请的多源遥感数据模型能够获得更加准确的土壤水分测量结果。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均在本发明待批权利要求保护范围之内。
1.一种基于多源遥感数据融合的土壤水分测量方法，用于对目标区域土壤所含的水分进行测量，其特征在于，
通过旋风全球导航卫星系统获取所述目标区域土壤的原始数据，通过所述原始数据和所述目标区域土壤的实测数据确定所述目标区域的土壤水分估算结果；
通过中分辨率成像光谱分析仪获取所述目标区域土壤的第一地表温度数据和通过中国气象局陆面数据同化系统获取所述目标区域土壤的第二地表温度数据，将所述第一地表温度数据和所述第二地表温度数据进行时空融合，得到地表温度估算结果；
获取影响土壤水分的辅助参数；所述辅助参数包括植被覆盖指数、地面高程、植被覆盖度、土壤粗糙度、地面坡度、降水；
将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行数据融合，得到多源遥感数据融合模型，通过所述多源遥感数据融合模型得到土壤水分测量结果。
2.根据权利要求1所述的基于多源遥感数据融合的土壤水分测量方法，所述土壤水分估算结果的确定过程为：
通过cygnss获取所述目标区域土壤的原始数据，并对所述原始数据进行预处理，得到标准数据；
利用双基雷达方程对所述标准数据进行处理，得到地表有效反射率；
获取所述目标区域土壤的实测数据，通过所述实测数据和所述地表有效反射率确定所述土壤水分估算结果。
3.根据权利要求2所述的基于多源遥感数据融合的土壤水分测量方法，所述预处理为采用数据质量控制方法对所述原始数据进行处理。
4.根据权利要求3所述的基于多源遥感数据融合的土壤水分测量方法，所述地表有效反射率为：
其中，sr为地表有效反射率；snr是峰值功率减去噪声后的信噪比值；prt是发射功率；gt是gps发射天线的增益；gr是镜面反射点的天线增益；λ表示gpsl1频段信号的波长；rts是从gnss卫星到镜面反射点的范围；rsr是从镜面反射点到接收器的距离；θ为入射角。
5.根据权利要求1所述的基于多源遥感数据融合的土壤水分测量方法，所述地表温度估算结果的确定过程为：
通过光谱分析仪获取所述目标区域土壤的第一地表温度数据和通过中国气象局陆面数据同化系统获取所述目标区域土壤的第二地表温度数据；
将所述第一地表温度数据和所述第二地表温度数据进行时空融合处理，得到所述地表温度估算结果。
6.根据权利要求5所述的基于多源遥感数据融合的土壤水分测量方法，所述时空模型是通过像元反射率的时间变化特征的改进型自适应遥感图像时空融合方法构建得到的，所述时空模型为estarfm时空模型。
7.根据权利要求5或6所述的基于多源遥感数据融合的土壤水分测量方法，通过所述时空模型进行的融合处理包括邻近相似像元搜索、权重值计算、计算转换系数。
8.根据权利要求1所述的基于多源遥感数据融合的土壤水分测量方法，所述辅助参数通过权重对比方法得到的。
9.根据权利要求1所述的基于多源遥感数据融合的土壤水分测量方法，执行将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行数据融合，得到多源遥感数据融合模型，通过所述多源遥感数据融合模型得到土壤水分测量结果的过程如下：
将所述土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行空间匹配；
将进行空间匹配后的土壤水分估算结果、所述地表温度估算结果和所述辅助参数进行处理，得到多源遥感数据融合模型；
通过所述多源遥感数据融合模型得到土壤水分测量结果。
10.一种基于多源遥感数据融合的土壤水分测量装置，包括处理器、存储器以及存储在所述存储器中并可在所述处理器上的运行的计算机程序，其特征在于，所述处理器在执行所述计算机程序时实现如权利要求1-9任一项所述的基于多源遥感数据融合的土壤水分测量方法的步骤。
技术总结