一种结合微波和红外遥感影像的土壤水分反演方法技术

本发明专利技术是一种结合微波和红外遥感影像融合的土壤水分反演方法。本发明专利技术结合PDI指数和VSWI指数的适用性特点，根据地表植被覆盖度的不同，建立了基于综合干旱指数（CDI）的遥感监测模型。先使用红外遥感数据得到综合干旱指数结果，然后结合微波遥感数据，求得综合干旱指数反演结果与微波遥感数据的相关系数，建立回归方程，影像融合同时对微波数据进行降尺度，结果可将红外数据反演的综合干旱指数结果转化为真实的土壤水分值，并且用微波遥感数据对红外数据的空值区域进行填补，即可得到完整的高分辨率土壤水分产品。该方法简单实用，提高了大范围土壤水分反演的效率。

Soil moisture retrieval method combining microwave and infrared remote sensing image

The invention relates to a method for soil moisture inversion combined with microwave and infrared remote sensing image fusion. According to the applicability of PDI index and VSWI index, according to the difference of vegetation coverage, a remote sensing monitoring model based on Integrated Drought Index (CDI) was established. The first use of infrared remote sensing data to obtain the comprehensive drought index results, and then combined with the microwave remote sensing data, the correlation coefficient inversion results comprehensive drought index and microwave remote sensing data, the regression equation, image fusion of the microwave and data downscaling, results the comprehensive drought index of red data inversion results into the real value of soil moisture and, using microwave remote sensing data of the infrared data null region to fill, you can get high resolution soil moisture products complete. The method is simple and practical, and it can improve the efficiency of soil moisture retrieval.

【技术实现步骤摘要】

本专利技术属于遥感图像处理领域，特别涉及一种结合微波和红外遥感影像融合的土壤水分反演方法。
技术介绍
陆表土壤含水量不仅是旱情监测的重要指标，还是气候、水文、生态、农业等领域的重要参数，也是全球气候变化的重要组成部分。遥感技术的发展为获取大范围地表土壤水分信息提供了有效手段。光学/热红外遥感主要基于土壤及土壤上植被的光谱反射特性及表面发射特性来估算土壤水分。当研究区域较大时，尽管光学/热红外遥感反演地表土壤水分有着较高的空间分辨率，但受云层的影响严重，造成了结果在空间上的不连续，不能有效对土壤水分进行业务化的监测。被动微波遥感有着较好的物理基础，受云层的影响不大，但星载被动微波传感器的空间分辨率相对较低，应用受到一定限制。本专利技术结合红外和微波遥感数据，旨在得到完整的高分辨率土壤水分产品，能提高大范围区域的陆表土壤水分监测效率，便于土壤水分产品的业务化生产。
技术实现思路
本专利技术根据土壤植被覆盖度的不同，综合使用垂直干旱指数(PDI)和植被供水指数(VSWI)建立综合干旱指数(CDI)监测模型。先使用红外遥感数据根据综合干旱指数模型得到反演结果，然后结合微波遥感数据，求得综合干旱指数模型反演结果与微波遥感土壤水分数据的相关系数，建立回归方程，进行影像融合同时对微波数据进行降尺度，将红外遥感数据反演的结果转化为真实的土壤水分值，并且用微波遥感数据对红外数据的空值区域进行填补，从而得到完整的高分辨率土壤水分产品。采用的具体技术方案如下：一种结合微波和红外遥感影像的土壤水分反演方法，包括以下步骤：步骤1：干旱指数PDI适合低植被覆盖区，VSWI适合高植被覆盖区，故基于红外遥感影像，以NDVI指数作为区别植被覆盖程度的阈值，建立综合PDI指数与VSWI指数的综合干旱指数；步骤2：进行综合干旱指数与微波遥感土壤水分数据的相关性分析，为得到最佳拟合方程，以完成两种数据的融合；步骤3：根据步骤2得到的最佳拟合方程，融合微波和红外影像的结果。本融合算法融合同时可以进行对微波数据的降尺度，并且用微波数据补全综合干旱指数的空值区，最终得到完整土壤水分产品。所述步骤1在基于红外遥感影像建立综合干旱指数时，具体包括子步骤：步骤1.1，计算影像的垂直干旱指数，其计算公式为： P D I = 1 M 2 + 1 ( R r e d + MR n i r ) , ]]>式中，Rred为经过大气校正的红光波段反射率，Rnir为经过大气校正的近红外波段反射率，M为土壤线的斜率；步骤1.2，计算影像的植被供水指数，其计算公式为：VSWI＝NDVI/Ts，式中，Ts为植被的冠层温度，NDVI是归一化植被指数。VSWI是根据植被覆盖状况的变化来进行反演的，所以不适用于低植被覆盖地区。步骤1.3，PDI适用于低植被覆盖区，VSWI适用于高植被覆盖区，以NDVI为区分植被覆盖区域高低的阈值，综合上两种指数建立综合干旱指数(CDI，compositive drought index)。综合干旱指数的定义公式如下：式中，CDIi表示在任意像元点i的综合干旱指数，PDIi表示在任意像元点i的垂直干旱指数，VSWIi表示在任意像元点i的植被供水指数，PDImax和PDImin分别表示垂直干旱指数最大值和最小值，VSWImax和VSWImin分别表示植被供水指数最大值和最小值。所述步骤2在进行综合干旱指数与微波遥感数据相关性分析时，进一步包括子步骤：步骤2.1：任取50组CDI像元值，求各组像元的平均值。红外遥感数据分辨率较高，此处取的每组包含若干个CDI像元(具体个数根据所采用红外和微波影像的分辨率比例来确定)，取均值是为了让均值代表此组内的CDI像元值；步骤2.2：取50个与步骤2.1对应位置的微波遥感土壤水分数据的像元值，与步骤2.1中的50个CDI均值建立拟合关系；步骤2.3：选取相关性最好的拟合模型函数。所述步骤3在融合微波/红外影像的结果时，进一步包括子步骤：步骤3.1，微波遥感数据质量控制。微波土壤水分数据会由于各种原因造成某个像元为空值，则该空值像元使用其周围8个像元值的平均值作为该像元值。按行列顺序逐一计算，计算公式如下： x θ = Σ i = 0 7 x i / ( 8 - k ) , ]]>式中，xθ表示该空值像元，xi表示该像元的邻接像元值，k(k∈[0,7]且k∈Z)为周围8个像元中空值像元的个数。步骤3.2，影像融合同时对微波数据进行降尺度，并把CDI的空值区补全，得到高分辨率的完整土壤水分产品。所述步骤3.2在融合CDI和微波土壤水分数据时，进一步包括子步骤：步骤3.2.1：由步骤2所得，微波遥感数据与综合干旱指数的最佳拟合关系为：θM＝f(θCDI)，此处假设为线性关系，即θM＝a+b×θCDI，式中，θM为红外影像网格尺度的一个像元代表的实际土壤含水量，θCDI为综合干旱指数上的一个像元值，a,b为系数；步骤3.2.2：取红外遥感影像网格上相邻i个像元的土壤含水量的平均值等于对应的微波遥感影像网格上一个像元的土壤含水量的值： θ S = ( θ M 1 + θ M 2 + ....... + θ M i ) / i = f ( θ M i ) , ]]>式中，θS为微波遥感网格尺度上的一个像元的土壤含水量，为红外遥感影像网格尺度上一个像元的土壤含水量；步骤3.2.3：把步骤3.2.2中的公式代入步骤3.2.1中，有 θ本文档来自技高网...

【技术保护点】
一种结合微波和红外遥感影像的土壤水分反演方法，其特征在于:包括以下步骤：步骤1：基于红外遥感影像建立以NDVI为阈值的综合PDI指数与VSWI指数的综合干旱指数；步骤2：进行综合干旱指数与微波遥感土壤水分数据的相关性分析；步骤3：融合微波和红外影像的结果，得到完整土壤水分产品。

【技术特征摘要】
1.一种结合微波和红外遥感影像的土壤水分反演方法，其特征在于:包括以下步骤：步骤1：基于红外遥感影像建立以NDVI为阈值的综合PDI指数与VSWI指数的综合干旱指数；步骤2：进行综合干旱指数与微波遥感土壤水分数据的相关性分析；步骤3：融合微波和红外影像的结果，得到完整土壤水分产品。2.如权利要求书1所述的一种结合微波和红外遥感影像的土壤水分反演方法，其特征在于：所述步骤1在基于红外遥感影像建立综合干旱指数，具体包括以下子步骤：步骤1.1：计算红外影像的垂直干旱指数，其计算公式为： P D I = 1 M 2 + 1 ( R r e d + MR n i r ) , ]]>式中，Rred为经过大气校正的红光波段反射率，Rnir为经过大气校正的近红外波段反射率，M为土壤线的斜率；步骤1.2:计算红外影像的植被供水指数，其计算公式为：VSWI＝NDVI/Ts，式中，Ts为植被的冠层温度，NDVI是归一化植被指数；步骤1.3:PDI适用于低植被覆盖区，VSWI适用于高植被覆盖区，以NDVI为阈值，综合此两种指数建立综合干旱指数CDI，综合干旱指数的定义公式如下： CDI i = 1 - PDI i - PDI min PDI max - PDI min ( NDVI i ≤ 0.3 ) VSWI i - VSWI min VSWI max - VSWI min ( NDVI i > 0.3 ) , ]]>式中，CDIi表示在任意像元点i的综合干旱指数，PDIi表示在任意像元点i的垂直干旱指数，VSWIi表示在任意像元点i的植被供水指数，PDImax和PDImin分别表示垂直干旱指数最大值和最小值，VSWImax和VSWImin分别表示植被供水指数最大值和最小值。3.如权利要求书2所述的一种结合微波和红外遥感影像的土壤水分反演方法，其特征在于：所述步骤2在进行综合干旱指数与微波遥感数据相关性分析，进一步包括以下子步骤：步骤2.1：任取50组综合干旱指数CDI像元值；步骤2.2：取对应的微波遥感土壤水分数据的像元值，建立拟合关系；步骤2.3：选取相关性最好的拟合模型函数。4.如权利要求书3所述的一种结合微波和红外遥感影像的土壤水分反演方法，其特征在于：所述步骤3在融合微波和红外影像的结果，进一步包括以下子步骤：步骤3.1：微波遥感数据质量控制：若微波土壤水分数据中的个像元为空值，则该空值像元使用其周围8个像元值的平均值作为该像元值；计算公式如下： x θ = Σ i = 0 7 x i / ( 8 - k ) , ]]>式中，xθ表示该空值像元，xi表示该像元的邻接像元值，k为周围8个像元中空值的个数,k∈[0,7]且k∈Z；步骤3.2：影像融合同时进行降尺度，以获得高分辨率的完整土壤水分产品。5.如权利要求书4所述的一种结合微波和红外遥感影像的土壤水分反演方法，其特征在于：所述步骤3.2在影像融合同时进行降尺度，进一步包括以下子步骤：步骤3.2.1：由步骤2所得，微波遥感数据与综合干旱指数的最佳拟合关系为：θM＝f(θCDI)，此处假设为线性关系，即θM＝a+b×θCDI，式中，θM为红外影像网格尺度上的一个像元代表的实际土壤含水量，θCDI为综合干旱指数上一个像元的值，a,b为系数；步骤3.2.2：取红外遥感影像网格上相邻i个像元的土壤含水量的平均值为微波遥感影像网格上一个像元的土壤含水量的值： θ S = ( θ M 1 ...

【专利技术属性】
技术研发人员：张文，任燕，谢文君，胡一杰，孟令奎，余长慧，马晓琳，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42