基于反射率-植被覆盖度二维空间的蒸散发遥感反演方法和应用

摘要

本发明属于地表蒸散发遥感估算技术领域，涉及基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法和应用。该方法包括：获取目标区域在预设时间尺度的相关数据；确定基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义，以及反射率－植被覆盖度二维空间中蒸发比的变化规律；根据长时间序列遥感短波红外反射率数据和植被覆盖度数据，确定反射率－植被覆盖度二维空间的干湿边，构建反射率－植被覆盖度的二维空间；根据扩展的P－T方程，以及构建的反射率－植被覆盖度二维空间中蒸发比的变化，确定与蒸发比相关的参数φ；基于所确定的参数φ，计算预设时间尺度的蒸散发。本发明对于准确地利用遥感手段获取地表蒸散发具有重要意义。

权利要求

1.一种基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，其特征在于，所述方法包括以下步骤：获取目标区域在预设时间尺度的相关数据；

确定基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义，以及反射率－植被覆盖度二维空间中蒸发比的变化规律；

根据长时间序列遥感短波红外反射率数据和植被覆盖度数据，确定反射率－植被覆盖度二维空间的干湿边，构建反射率－植被覆盖度的二维空间；

根据扩展的Priestley－Taylor方程，以及所构建的反射率－植被覆盖度二维空间中蒸发比的变化，确定与蒸发比相关的参数φ；

基于所确定的参数φ，计算所述时间尺度的蒸散发；

所述相关数据包括目标区域在预设时间尺度的气象数据、预设时间尺度的遥感反演地表净辐射和土壤热通量数据，以及长时间序列遥感短波红外反射率数据和植被覆盖度数据；

基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义包括：通过解译长时间序列的遥感短波红外反射率数据和植被覆盖度数据构建的梯形散点图，进行蒸散发的估算；所述方法由干边和湿边构成反射率－植被覆盖度特征空间的两个物理边界；

所述反射率－植被覆盖度二维空间中蒸发比的变化规律包括：所述方法中湿边像元处于潜在蒸发状态，蒸发比为最大值；所述方法中干边上的蒸散发主要来源于植被的蒸腾，干边上的蒸发比从0到最大值呈线性变化；

所述扩展的Priestley－Taylor方程的表达式如下：

其中，LE为潜热通量，φ为考虑空气动力学阻抗作用的综合参数，Rn为地表净辐射，G为土壤热通量，Δ为饱和水汽压曲线斜率，γ为干湿球常数；

根据蒸发比的定义，蒸发比与参数φ的关系为：

其中，EF为蒸发比，φ为考虑空气动力学阻抗作用的综合参数，Δ为饱和水汽压曲线斜率，γ为干湿球常数。

2.根据权利要求1所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，其特征在于，所述气象数据包括空气温度数据、大气压强数据和相对湿度数据。

3.根据权利要求1所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，其特征在于，根据所述长时间序列遥感短波红外反射率数据和植被覆盖度数据构成的二维空间散点图，确定所述二维空间的干湿边，构建反射率－植被覆盖度的二维空间，所述干湿边的公式为：STRd＝md+rd*Fr；

STRw＝mw+rw*Fr；

其中，STRd和STRw分别为干边和湿边上短波红外转换反射率，表达式为STR＝(1-R)2/(2R)，R为短波红外反射率；md和mw分别为极干裸土和极湿裸土处短波红外转换反射率；rd和rw分别为极干全植被覆盖区域和极湿全植被覆盖区域短波红外转换反射率与对应裸土处短波红外转换反射率的差值；Fr为植被覆盖度。

4.根据权利要求1所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，其特征在于，参数φ的确定是通过在梯形空间中双线性插值得到，具体步骤包括：假定湿边上φ值随着植被覆盖度的变化而保持不变，一直为最大值φmax，i；

假定干边最干燥裸土像元的φ值为全局最小值φmin，干边上φ值随着像元的植被覆盖度的增加而线性增加；

假定在某一给定像元的植被覆盖度情况下，φ值随着短波红外转换反射率STR的增加而从φmin，i到φmax，i线性增加。

5.根据权利要求4所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，其特征在于，在特征空间内任一像元(Fr，STR)的φ值计算公式为：φmin,i＝1.26*Fr；

φmax,i＝φmax＝1.26；

其中，STRd和STRw分别为干边和湿边上短波红外转换反射率，STR为待计算像元的短波红外转换反射率，φmin，i为植被覆盖度等于Fr时对应的φ最小值，φmax，i为植被覆盖度等于Fr时对应的φ最大值，Fr为植被覆盖度。

6.根据权利要求1～5任一项所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，其特征在于，基于所确定的参数φ，以及所述时间尺度的气象数据、遥感反演地表净辐射和土壤热通量数据，利用扩展的Priestley－Taylor方程计算所述时间尺度的蒸散发。

7.权利要求1～6任一项所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法在地表蒸散发研究领域中的应用。

说明书

基于反射率-植被覆盖度二维空间的蒸散发遥感反演方法和

应用

技术领域

[0001] 本发明属于地表蒸散发遥感估算技术领域，具体而言，涉及一种基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法和应用。

背景技术

[0002] 蒸散发(Evapotranspiration，ET)是地表水循环与能量平衡的重要组成部分，决定了土壤－植被－大气系统中的水、热传输，其定量估算在水文学、气象学、农学、地学等学科研究中，都发挥着重要的作用。而遥感技术从技术性、经济性和实效性等方面来看，被认为是获得地球表面区域尺度蒸散发分布最有效的方法。

[0003] 目前，已发展的遥感蒸散发估算方法有一源模型、二源模型和地表温度－植被指数三角形/梯形模型。地表温度－植被指数三角形/梯形模型进行地表蒸散发的遥感估算，通常只需要遥感数据和少量气象数据作为输入，方法简单，已得到了广泛的应用。但是传统的地表温度－植被指数三角形/梯形模型具有以下局限性：(1)需要可见光与热红外相配套的数据，这样直接排除了没有热红外波段而只有可见光波段卫星遥感传感器的应用；(2)由于地表温度不仅受土壤含水量的影响，还受太阳辐射、空气温度、相对湿度和风速等气象要素的影响，故无法仅利用卫星单一过境数据或长时间序列数据开展地表温度－植被指数三角形/梯形模型的干边和湿边确定，而需要耗时进行每一天参数的计算和校准，限制了该方法的应用。

[0004] 有别于热红外地表温度数据，遥感短波红外波段的反射率数据，只与下垫面水分含量有关，不随周围大气参数的变化而变化，因此，利用遥感短波红外反射率数据构建方法可以克服上述传统方法的两个局限性。

[0005] 鉴于此，特提出本发明。

发明内容

[0006] 本发明的第一目的在于提供一种基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，以克服上述问题或者至少部分地解决上述技术问题。

[0007] 本发明的第二目的在于提供一种所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法在地表蒸散发研究领域中的应用。

[0008] 为实现上述目的，本发明采用的技术方案为：

[0009] 根据本发明的一个方面，本发明提供一种基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，所述方法包括以下步骤：

[0010] 获取目标区域在预设时间尺度的相关数据；

[0011] 确定基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义，以及反射率－植被覆盖度二维空间中蒸发比的变化规律；

[0012] 根据长时间序列遥感短波红外反射率数据和植被覆盖度数据，确定反射率－植被覆盖度二维空间的干湿边，构建反射率－植被覆盖度的二维空间；

[0013] 根据扩展的Priestley－Taylor方程，以及所构建的反射率－植被覆盖度二维空间中蒸发比的变化，确定与蒸发比相关的参数φ；

[0014] 基于所确定的参数φ，计算所述时间尺度的蒸散发。

[0015] 作为进一步优选技术方案，所述相关数据包括目标区域在预设时间尺度的气象数据、预设时间尺度的遥感反演地表净辐射和土壤热通量数据，以及长时间序列遥感短波红外反射率数据和植被覆盖度数据；

[0016] 优选地，所述气象数据包括空气温度数据、大气压强数据和相对湿度数据。

[0017] 作为进一步优选技术方案，基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义包括：

[0018] 通过解译长时间序列的遥感短波红外反射率数据和植被覆盖度数据构建的梯形散点图，进行蒸散发的估算；

[0019] 所述方法由干边和湿边构成反射率－植被覆盖度特征空间的两个物理边界。

[0020] 作为进一步优选技术方案，所述反射率－植被覆盖度二维空间中蒸发比的变化规律包括：

[0021] 所述方法中湿边像元处于潜在蒸发状态，蒸发比为最大值；

[0022] 所述方法中干边上的蒸散发主要来源于植被的蒸腾，干边上的蒸发比从0到最大值呈线性变化。

[0023] 作为进一步优选技术方案，根据所述长时间序列遥感短波红外反射率数据和植被覆盖度数据构成的二维空间散点图，确定所述二维空间的干湿边，构建反射率－植被覆盖度的二维空间，所述干湿边的公式为：

[0024] STRd＝md+rd*Fr；

[0025] STRw＝mw+rw*Fr；

[0026] 其中，STRd和STRw分别为干边和湿边上短波红外转换反射率，表达式为STR＝(1-R)2/(2R)，R为短波红外反射率；md和mw分别为极干裸土和极湿裸土处短波红外转换反射率；rd和rw分别为极干全植被覆盖区域(表层土壤水分为0，但根区土壤水分无亏缺)和极湿全植被覆盖区域短波红外转换反射率与对应裸土处短波红外转换反射率的差值；Fr(Fraction cover)为植被覆盖度。

[0027] 作为进一步优选技术方案，所述扩展的Priestley－Taylor方程的表达式如下：

[0028]

[0029] 其中，LE(Latent heat flux)为潜热通量，φ为考虑空气动力学阻抗作用的综合参数，Rn为地表净辐射，G为土壤热通量，Δ为饱和水汽压曲线斜率，γ为干湿球常数。

[0030] 作为进一步优选技术方案，根据蒸发比的定义，蒸发比与参数φ的关系为：

[0031]

[0032] 其中，EF(Evapotranspiration Fraction)为蒸发比，φ为考虑空气动力学阻抗作用的综合参数，Δ为饱和水汽压曲线斜率，γ为干湿球常数。

[0033] 作为进一步优选技术方案，参数φ的确定是通过在梯形空间中双线性插值得到，具体步骤包括：

[0034] 假定湿边上φ值随着植被覆盖度的变化而保持不变，一直为最大值φmax，i；

[0035] 假定干边最干燥裸土像元的φ值为全局最小值φmin，干边上φ值随着像元的植被覆盖度的增加而线性增加；

[0036] 假定在某一给定像元的植被覆盖度情况下，φ值随着短波红外转换反射率STR(Shortwave Transformed Reflectivity)的增加而从φmin，i到φmax，i线性增加；

[0037] 优选地，在特征空间内任一像元(Fr，STR)的φ值计算公式为：

[0038]

[0039] φmin,i＝1.26*Fr；

[0040] φmax,i＝φmax＝1.26；

[0041] 其中，STRd和STRw分别为干边和湿边上短波红外转换反射率，STR 为待计算像元的短波红外转换反射率，φmin，i为植被覆盖度等于Fr时对应的φ最小值，φmax，i为植被覆盖度等于Fr时对应的φ最大值，Fr 为植被覆盖度。

[0042] 作为进一步优选技术方案，基于所确定的参数φ，以及所述时间尺度的气象数据、遥感反演地表净辐射和土壤热通量数据，利用扩展的Priestley－Taylor方程计算所述时间尺度的蒸散发。

[0043] 根据本发明的另一个方面，本发明提供一种以上所述的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法在地表蒸散发研究领域中的应用。

[0044] 与现有技术相比，本发明的有益效果在于：

[0045] 本发明的基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，可以利用遥感手段准确获取目标时间尺度的地表蒸散发，该方法利用地表短波红外反射率只与土壤状况有关，而不受大气状况影响的特征，可利用长时间序列的短波红外反射率和植被覆盖度数据构建唯一的一个二维特征空间，准确确定二维空间的干湿边，从而避免了仅利用所述时间尺度数据确定干湿边时，每个时间尺度数据均各自构成一个二维特征空间，由于所述时间尺度参数的计算和校准而造成的干湿边确定的不确定性。同时本发明也避免了传统的地表温度－植被指数三角形/梯形模型中的以下局限性：(1)需要可见光与热红外相配套的数据，直接排除了没有热红外波段而只有可见光波段卫星传感器的应用；(2)由于地表温度不仅受土壤含水量的影响，还受太阳辐射、空气温度、相对湿度和风速等气象要素的影响，故无法仅利用卫星单一过境数据或长时间序列数据开展地表温度－植被指数三角形/梯形模型的干边和湿边确定，而需要耗时进行每一天参数的计算和校准，限制了该方法的应用。本发明对于准确地利用遥感手段获取地表蒸散发具有重要意义。

[0046] 本发明能够高效准确地利用遥感手段准确获取目标时间尺度的地表蒸散发，具有方法简单、灵活、易于操作，可操作性强，结果可靠性好等特点。

[0047] 本发明的其他特征和优点将在随后的说明书中阐述，或者，部分特征和优点可以从说明书推知或毫无疑义地确定，或者通过实施本公开的上述技术即可得知。

[0048] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，作详细说明如下。

具体实施方式

[0049] 下面将结合实施方式和实施例对本发明的实施方案进行详细描述，但是本领域技术人员将会理解，下列实施方式和实施例仅用于说明本发明，而不应视为限制本发明的范围。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。未注明具体条件者，按照常规条件或制造商建议的条件进行。

[0050] 在至少一个实施例中提供一种基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，所述方法包括以下步骤：

[0051] 获取目标区域在预设时间尺度的相关数据；

[0052] 确定基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义，以及反射率－植被覆盖度二维空间中蒸发比的变化规律；

[0053] 根据长时间序列遥感短波红外反射率数据和植被覆盖度数据，确定反射率－植被覆盖度二维空间的干湿边，构建反射率－植被覆盖度的二维空间；

[0054] 根据扩展的Priestley－Taylor方程，以及所构建的反射率－植被覆盖度二维空间中蒸发比的变化，确定与蒸发比相关的参数φ；

[0055] 基于所确定的参数φ，计算所述时间尺度的蒸散发。

[0056] 可以理解的是，本发明中的蒸发比为本领域的常用术语，通常是指蒸散发与可利用能量的比值。另外，本发明中其他未详细解释说明的术语，均按照本领域的常用术语理解即可，本发明对此不作特殊限制。

[0057] 本发明能够高效准确地利用遥感手段准确获取目标时间尺度的地表蒸散发，具有方法简单、灵活、易于操作，可操作性强，结果可靠性好等特点，对于准确地利用遥感手段获取地表蒸散发具有重要意义。

[0058] 在本发明的一种优选实施方式中，所述方法包括以下步骤：

[0059] 获取目标区域在预设时间尺度的气象数据、所述预设时间尺度的遥感反演地表净辐射和土壤热通量数据、以及长时间序列遥感短波红外反射率和植被覆盖度数据；所述气象数据包括：空气温度数据、大气压强数据、相对湿度数据；

[0060] 确定基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义、以及反射率－植被覆盖度二维空间中蒸发比(蒸散发与可利用能量的比值)的变化规律；

[0061] 根据长时间序列遥感短波红外反射率数据和植被覆盖度数据，确定反射率－植被覆盖度二维空间的干湿边，构建反射率－植被覆盖度的二维空间；

[0062] 根据扩展的Priestley－Taylor方程，根据所述构建的反射率－植被覆盖度二维空间中蒸发比的变化，确定与蒸发比相关的参数φ；

[0063] 根据所述确定的参数φ，所述时间尺度的遥感地表净辐射和土壤热通量数据，所述时间尺度的气象数据，利用所述扩展的 Priestley－Taylor方程计算所述时间尺度的蒸散发。

[0064] 需要说明的是，本发明的方法在计算过程中所需要的气象数据、遥感反演地表净辐射和土壤热通量数据、以及长时间序列遥感短波红外反射率和植被覆盖度数据，但本发明不对这些具体的数据做具体限定；以及空气温度数据、大气压强数据、相对湿度数据等气象数据，本发明在此不一一列举，对于这些具体的气象数据也不做具体限定。

[0065] 应当理解的是，本发明所述方法的说明中未详细描述的内容，均是本领域技术人员容易想到的常用参数或操作方式，或是本领域常规的条件，因此可以省略对其的详细说明。

[0066] 为便于对本实施例进行理解，下面对本发明实施例所公开的一种一种估算蒸散发的反射率－植被覆盖度梯形方法进行详细介绍。

[0067] 实施例

[0068] 一种基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，利用该方法估算蒸散发包括以下步骤：

[0069] 步骤S101：获取目标区域在预设时间尺度的气象数据，所述预设时间尺度的遥感反演地表净辐射和土壤热通量数据，以及长时间序列遥感短波红外反射率和植被覆盖度数据；

[0070] 其中，所述气象数据包括：空气温度数据、大气压强数据和相对湿度数据。

[0071] 步骤S102：确定基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法的原理与定义、以及反射率－植被覆盖度二维空间中蒸发比 (蒸散发与可利用能量的比值)的变化规律。

[0072] 所述基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，即通过解译长时间序列的遥感短波红外反射率数据和植被覆盖度数据构建的梯形散点图(以反射率为纵坐标)，进行蒸散发的估算。所述方法由干边和湿边构成反射率－植被覆盖度特征空间的两个物理边界。对一个给定的植被覆盖度，土壤水分从干边到湿边不断增加，地表短波红外反射率由最小值增加到最大值。方法中湿边像元处于潜在蒸发(腾发)状态，即蒸发比(蒸散发与可利用能量的比值)EF为最大值；而干边上的蒸散发主要来源于植被的蒸腾，干边上的蒸发比 EF从0(裸土)到最大值线性变化。

[0073] 步骤S103：根据所述长时间序列遥感短波红外反射率数据和植被覆盖度数据，确定反射率－植被覆盖度二维空间的干湿边，构建反射率－植被覆盖度的二维空间。

[0074] 由于地表短波红外反射率只与下垫面水分含量有关，而不受大气状况影响，因此可以由长时间的短波红外反射率和植被覆盖度数据形成一个二维空间。

[0075] 根据所述长时间序列遥感短波反射率数据和植被覆盖度数据构成的二维空间散点图，确定所述二维空间的干湿边，构建构建反射率－植被覆盖度的二维空间，所述干湿边的公式为：

[0076] STRd＝md+rd*Fr；

[0077] STRw＝mw+rw*Fr；

[0078] 其中，STRd和STRw分别为干边和湿边上短波红外转换反射率，表达式为STR＝(1-R)2/(2R)，R为短波红外反射率；md和mw分别为极干裸土和极湿裸土处短波红外转换反射率；rd和rw分别为极干全植被覆盖区域(表层土壤水分为0，但根区土壤水分无亏缺)和极湿全植被覆盖区域短波红外转换反射率与对应裸土处短波红外转换反射率的差值；Fr为植被覆盖度。

[0079] 步骤S104：根据扩展的Priestley－Taylor方程，根据所述构建的反射率－植被覆盖度二维空间中蒸发比的变化，确定与蒸发比相关的参数φ。

[0080] 所述扩展的Priestley－Taylor方程的表达式如下：

[0081]

[0082] 其中，LE为潜热通量，φ为考虑空气动力学阻抗作用的综合参数，Rn为地表净辐射，G为土壤热通量，Δ为饱和水汽压曲线斜率，γ为干湿球常数。

[0083] 根据蒸发比的定义(蒸散发与可利用能量的比值，可利用能量为地表净辐射与土壤热通量的差值)，蒸发比与参数φ的关系为：

[0084]

[0085] 其中，EF为蒸发比，φ为考虑空气动力学阻抗作用的综合参数，Δ为饱和水汽压曲线斜率，γ为干湿球常数。

[0086] 参数φ的确定是通过在梯形空间中双线性插值得到，具体步骤包括：

[0087] (1)：假定湿边上φ值随着植被覆盖度的变化而保持不变，一直为最大值1.26(φmax，i＝1.26)；

[0088] (2)：假定干边最干燥裸土像元的φ值为全局最小值(φmin＝0)，干边上φ值随着Fr的增加而线性增加；

[0089] (3)：假定在某一给定的Fr情况下φ值随着STR的增加而从φmin，i到φmax，i线性增加。

[0090] 即，在特征空间内任一像元(Fr，STR)的φ值计算公式为：

[0091]

[0092] φmin,i＝1.26*Fr；

[0093] φmax,i＝φmax＝1.26；

[0094] 其中，STRd和STRw分别为干边和湿边上短波红外转换反射率，STR 为待计算像元的短波红外转换反射率，φmin，i为植被覆盖度等于Fr时对应的φ最小值，φmax，i为植被覆盖度等于Fr时对应的φ最大值，Fr 为植被覆盖度。

[0095] 步骤S105：根据所述确定的参数φ，所述时间尺度的遥感地表净辐射和土壤热通量数据，所述时间尺度的气象数据，利用所述扩展的Priestley－Taylor方程计算所述时间尺度的蒸散发。

[0096] 本发明实施例提供的一种基于反射率－植被覆盖度二维空间的蒸散发遥感反演方法，有效的解决了传统的地表温度－植被指数三角形/ 梯形模型的以下局限性：(1)需要可见光与热红外相配套的数据，直接排除了没有热红外波段而只有可见光波段卫星遥感传感器的应用； (2)由于地表温度不仅受土壤含水量的影响，还受太阳辐射、空气温度、相对湿度和风速等气象要素的影响，故无法仅利用卫星单一过境数据或长时间序列数据开展地表温度－植被指数三角形/梯形模型的干边和湿边确定，而需要耗时进行每一天参数的计算和校准，限制了该方法的应用。本方法选择了地表短波红外反射率，只需要可见光波段数据即可进行蒸散发的估算；同时本方法利用了地表短波红外反射率只与下垫面水分含量有关，而不受大气状况影响的特征，利用长时间序列的短波红外反射率和植被覆盖度数据构成唯一的二维特征空间，准确确定二维空间的干湿边，避免了仅利用所述时间尺度数据确定干湿边时，由于所述时间尺度参数的计算和校准而造成的干湿边确定的不确定性。本发明对于准确地利用遥感手段获取地表蒸散发具有重要意义。

[0097] 最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。