一种消除植被冠层结构和地表背景影响的叶片反射率卫星遥感提取方法

摘要

本发明提供了一种消除植被冠层结构和地表背景影响的叶片反射率卫星遥感提取方法，属于植被遥感反演参数方法的研究领域。其步骤为：利用4‑Scale模型确定遥感影像像元光谱对应的阳叶的可视概率(PT)、光照背景的可视概率(PG)；计算光谱的角度指数(AI)，与PT、PG进行相关分析，建立基于AI的PT和PG的估算模型；以叶面积指数(LAI)、PT和PG作为3个检索项，利用4‑Scale模型模拟建立针对多次散射因子M的查找表(与波长相关)；最后，利用建立的估算模型和查找表，求得卫星遥感影像植被像元的平均叶片反射率。本发明方法可由遥感影像提取叶片的反射光谱，与优化迭代计算方法相比，计算效率提高，与需多步查找的查找表方法相比，计算过程更简单，效率更高。

说明书

一、技术领域

本发明是利用卫星遥感影像数据，定量反演像元的阳叶和光照背景可视概率，以多次散射反射因子描述阴影分量和多次散射反射的影响，利用基于4尺度几何光学模型理论发展的光谱降尺度模型提取像元平均叶片反射率的方法，属于植被遥感反演参数方法的研究领域。

二、背景技术

卫星观测的像元尺度反射信号受到叶片反射率、冠层结构和背景反射率等因子的影响，利用卫星遥感数据定量反演植被叶片和冠层生化参数，首先需要消除植被冠层结构和地表背景的影响。由冠层反射率光谱反演叶片反射率光谱，可以消除植被冠层结构和地表背景的影响，是利用卫星遥感数据高精度定量反演植被叶片和冠层生化参数的关键步骤。

为了消除植被冠层结构和地表背景对植被冠层生化参数遥感反演的影响，一些研究以几何光学模型为工具，根据冠层辐射传输原理，通过对冠层顶的反射光谱进行分解，获取叶片的反射率光谱。这些研究一般采用优化迭代计算或建立查找表的方法对冠层顶的反射光谱进行分解，例如，Zarco-Tejada等人在“Water content estimation in vegetationwith MODIS reflectance data and model inversion methods”文中利用SAILH冠层辐射传输模型研究叶面积指数(LAI)和观测几何角度的影响，并且利用优化迭代技术由冠层反射率光谱反演叶片反射率光谱；Zhang等人在“Leaf chlorophyll content retrievalfrom airborne hyperspectral remote sensing imagery”文中利用4尺度几何光学模型(4-Scale模型)研究冠层结构参数、观测几何角度以及地表背景反射率等因素的影响，引入了一个多次散射因子(M)，对4-Scale模型进行简化，提出了光谱降尺度模型，以太阳天顶角(SZA)、观测天顶角(VZA)和观测方位角(PHI)3个观测几何参数和叶面积指数(LAI)作为4个检索项来建立阳叶可视概率(PT)、光照背景可视概率(PG)与多次散射反射因子(M)的查找表，利用查找表来实现叶片反射率光谱的反演。

这两种方法都是基于遥感物理机理的物理模型，其中优化迭代计算方法的计算效率低，在大区域的卫星遥感影像应用中可操作性不强；查找表方法通过预先进行大量的模拟计算，建立查找表来简化计算过程，提高了反演计算的效率，但要求输入每个像元的冠层结构参数和观测几何角度作为检索项，并需要进行多步查找，过程繁琐，而且查找表是有局限性的，即使在其他查找项固定的情况下，PT和PG与LAI也不是一一映射的关系，查找表给出的是最大可能概率下的值。

三、发明内容

本发明的目的是提供一种适用于遥感影像的简单高效的模型方法，由遥感影像提取叶片的反射光谱，为有效提高叶片和冠层的生化参数的遥感反演精度提供基础。本发明的原理如下：

本发明分析利用叶片反射率和背景反射率在光谱形态上的差异构建角度指数(AI)，据其估算阳叶可视概率(PT)、光照背景可视概率(PG)，在此基础上结合叶面积指数(LAI)确定多次散射因子(M)，根据冠层反射率光谱和光谱降尺度模型求单叶片的反射率光谱。

本发明的技术方案主要包括以下步骤：

(1)首先，将地面样区的地理坐标与遥感影像匹配，提取对应像元的光谱，利用4-Scale模型确定遥感光谱对应的阳叶的可视概率(PT)、光照背景的可视概率(PG)。冠层反射率光谱的准确模拟是进行冠层光谱分解，获取叶片的反射率光谱的基础和前提条件。4-Scale模型是冠层几何光学物理模型，根据像元内树冠和背景，以及光照面和阴影面定义四分量：阳叶、阴叶、光照背景和阴影背景。4-Scale模型的输入参数包括：植被冠层结构参数，观测几何参数、叶片反射率、叶片透射率和背景反射率等。其中一组植被冠层结构参数和观测几何参数输入值代表一个观测场景。根据地面观测站点的实验数据设置4-Scale模型的输入值：参照地面观测实验得到的冠层结构参数统计特征，设置模型的植被冠层结构参数输入值范围；按照传感器的地面像元大小和成像条件来设置模型的观测几何参数输入值范围；利用地面观测实验得到的叶片反射率、叶片透射率和背景反射率光谱作为4-Scale模型的光谱数据输入。模拟过程分为两个部分：先根据输入的植被冠层结构参数，观测几何参数等建立模拟观测场景，计算模拟场景中四分量的可视概率：阳叶的可视概率(PT)、阴叶的可视概率(ZT)、光照背景的可视概率(PG)和阴影背景的可视概率(ZG)。然后根据输入的叶片反射率，叶片透射率和背景反射率等光谱数据计算多次散射反射率等，模拟计算四分量的反射率(阳叶的反射率(ρ_PT)、阴叶的反射率(ρ_ZT)、光照背景的反射率(ρ_PG)和阴影背景的反射率(ρ_ZG))，最后利用公式(1)模拟计算冠层反射率(ρ)。在模拟过程中，根据遥感与模拟的冠层光谱差别最小的原则(两者之间比较计算得到的调整后的R平方最大，归一化均方根误差nRMSE最小)对冠层结构参数进行寻优，选取与遥感光谱一致性最好的模拟光谱所对应的阳叶的可视概率(PT)、光照背景的可视概率(PG)、叶面积指数(LAI)和其他模型模拟参数值作为遥感光谱对应的植被冠层结构参数和四分量参数的值。

ρ＝ρ_PTPT+ρ_ZTZT+ρ_PGPG+ρ_ZGZG>

(2)角度指数(AI)估算阳叶可视概率(PT)、光照背景可视概率(PG)。叶片反射率和背景反射率在光谱形态上差异明显，而冠层反射率的光谱形态主要由阳叶和光照背景这两个分量决定，所以可根据冠层反射率的光谱形态特征来估算这两个分量的可视概率PT和PG。利用红光(645nm)、近红外(858nm)和短波红外(1640nm)波段和对应的反射率建立的角度指数(AI)来刻画冠层反射率的光谱形态特征。以3个波段波长(λ₁、λ₂和λ₃)和对应的反射率(f₁、f₂和f₃)构成的二维空间里的3个点(λ₁，f₁)、(λ₂，f₂)和(λ₃，f₃)来定义以(λ₂，f₂)为顶点的角度指数(AI)：

设置λ₁＝645nm，λ₂＝858nm，λ₃＝1640nm，根据公式(2)，计算遥感影像像元光谱的角度指数(AI)，并与步骤(1)中得到的阳叶的可视概率(PT)、光照背景的可视概率(PG)进行相关分析，建立基于AI的PT和PG的估算模型。

(3)多次散射因子(M)查找表的建立。多次散射因子(M)是为了简化4-Scale模型，实现将冠层反射率数据降尺度得到叶片反射率光谱而提出的变量，包括阴影分量和多次散射反射的影响。理论上，M可用下式计算：

利用4-Scale模型模拟建立针对M的查找表，首先，根据传感器的地面像元大小和成像条件来设置若干组符合客观现实的观测几何参数、植被冠层结构参数，再选择任意一组背景反射率ρ_PG和叶片反射率ρ_L，利用4-Scale模型模拟得到冠层反射率ρ、PT和PG，代入上式计算M。然后，根据模拟结果以叶面积指数(LAI)、阳叶可视概率(PT)和光照背景可视概率(PG)作为3个检索项，建立多个M的查找表(与波长相关)。

(4)像元平均叶片反射率的提取。以4尺度几何光学模型为基础构建光谱降尺度模型：

在背景反射率ρ_PG已知的情况下，根据冠层反射率光谱ρ计算AI，利用基于AI的PT和PG的估算模型求得阳叶的可视概率PT、光照背景的可视概率PG，结合LAI数据查找给定波长的多次散射因子M，所有结果代入公式中，求得单叶片的反射率光谱ρ_L。

对于卫星观测的遥感影像，先进行遥感影像的预处理，包括大气校正、几何校正、像元分类和非植被像元掩膜处理等，得到植被像元的冠层反射率光谱ρ。然后根据地表背景类型，选择土壤光谱数据库中对应的典型地物反射率作为背景反射率ρ_PG。再利用基于AI的PT和PG的估算模型估算植被像元的PT和PG。同时，根据已经利用地面数据校准过的LAI的遥感估算模型，进行LAI的遥感反演，或者采用已有的经过系统验证的LAI标准产品。进而基于PT、PG和LAI确定M。最后求得植被像元平均叶片反射率。

四、附图说明

图1为2005年3月云南勐遮星地同步实验研究区的卫星遥感影像灰度图

(a)为Hyperion高光谱影像数据的灰度图

(b)为MODIS多光谱影像数据的灰度图

图2为4尺度模型模拟的冠层反射率光谱与Hyperion影像上提取的像元光谱的比较示例图

图3为4尺度模型模拟的冠层反射率光谱与MODIS影像上提取的像元光谱的比较示例图

图4为云南研究区站点对应的Hyperion影像像元光谱计算得到的角度指数(AI)与4-Scale模型反演得到的PT和PG的相关关系

(a)为AI与PT的相关关系

(b)为AI与PG的相关关系

图5为云南研究区站点对应的MODIS影像像元光谱计算得到的角度指数(AI)与4-Scale模型反演得到的PT和PG的相关关系

(a)为AI与PT的相关关系

(b)为AI与PG的相关关系

图6为Hyperion影像数据反演的植被像元的平均叶片反射率光谱与地面实测的叶片反射率光谱的比较示例图

图7为MODIS影像数据反演的植被像元的平均叶片反射率光谱与地面实测的叶片反射率光谱的比较示例图

五、具体实施方式

以下通过具体实例对本发明作进一步描述：

以2005年3月云南勐遮星地同步实验研究区的Hyperion高光谱影像(图1(a))和MODIS多光谱影像(图1(b))为实例，根据技术方案步骤(1)所述，模拟地面观测站点的Hyperion和MODIS像元光谱(冠层反射率光谱)，图2对比了模拟的冠层反射率光谱与Hyperion影像上提取的光谱(将地面观测站点的地理坐标与影像匹配)，图3对比了模拟的冠层反射率光谱与MODIS影像上提取的光谱(将地面观测站点的地理坐标与影像匹配)，结果表明模拟的冠层光谱与卫星遥感影像像元光谱具有很好的一致性。在准确模拟冠层反射率光谱的基础上，确定Hyperion和MODIS遥感影像像元光谱对应的阳叶的可视概率(PT)、光照背景的可视概率(PG)。

根据技术方案步骤(2)所述，以Hyperion影像上提取的光谱和其对应的阳叶的可视概率(PT)、光照背景的可视概率(PG)作为训练数据，建立适用于Hyperion影像的基于AI的PT的估算模型(图4(a))和基于AI的PG的估算模型(图4(b))。同时，以MODIS影像上提取的光谱和其对应的阳叶的可视概率(PT)、光照背景的可视概率(PG)作为训练数据，建立适用于MODIS影像的基于AI的PT的估算模型(图5(a))和基于AI的PG的估算模型(图5(b))。

根据技术方案步骤(3)所述，分别建立适用于Hyperion影像的多次散射因子M的查找表(如表1)和适用于MODIS影像的多次散射因子M的查找表(如表2)。

根据技术方案步骤(4)所述，分别求得Hyperion和MODIS遥感影像植被像元的平均叶片反射率。图6是Hyperion影像数据反演的植被像元的平均叶片反射率光谱与地面实测的叶片反射率光谱的对比，图7是MODIS影像数据反演的植被像元的平均叶片反射率光谱与地面实测的叶片反射率光谱的对比。

根据本发明提出的方法，在Hyperion和MODIS不同空间分辨率传感器的像元上，角度指数(AI)与PT和PG都具有极好的相关性，由影像光谱反演的叶片光谱与实测数据高度相似，说明了构建的方法有应用到多种卫星传感器的潜力。

表1 适用于Hyperion数据的以PT、PG和LAI作为检索项的多次散射因子M查找表示例

表2 适用于MODIS数据的以PT、PG和LAI作为检索项的多次散射因子M查找表示例