[发明专利]一种基于反射光谱小波变换的植被参数遥感反演方法

说明书

技术领域

本发明涉及高光谱遥感监测和定量反演领域，特别涉及一种基于反射光谱小波变换的植被参数遥感反演方法。

背景技术

植被参数分类两类：生物物理参数和生物化学参数（生化参数）。生物物理参数包括生物量、叶面积指数、覆盖度、植株高度、叶倾角等；生化参数包括叶绿素、氮素、蛋白质、木质素等。不同的生化参数具有不同的吸收反射特性，表现为不同的吸收反射峰和谷，这些吸收反射峰谷具有大小、深浅等特征，即具有不同的尺度特性。目前的遥感反演是用原始光谱或是导数、伪吸收系数等变化光谱建立统计模型直接反演生化参数，或是利用机理模型反演生化参数。这种直接利用原始光谱及其变化形式的光谱作为反演生化参数的方法没有考虑生化参数本身吸收反射特征的尺度信息，因此不能较好的满足日益提高的精准农业对生化参数反演的精度要求。

发明内容

为了解决现有生化参数反演方法未考虑光谱尺度效应进行的生化参数反演精度不高的问题，本发明提供一种考虑光谱尺度特征、精度较高的基于反射光谱小波变换的植被生化参数遥感反演方法。

本发明解决其技术问题所采用的技术方案是： 

一种基于反射光谱小波变换的生化参数遥感反演方法，所述反演方法包括以下步骤：

1）光谱数据变换：获取不同条件的植被参数及其原始光谱，对原始光谱进行光谱变换；

2）小波变换：将原始光谱利用不同小波函数进行连续小波变换，生成具有不同频率的小波系数；原始光谱的连续小波变换公式如下：

               （1）

式中，a,b为实数，小波系数是连续的，连续小波变换在任何尺度上进行。

3）构建逐步回归模型：以不同尺度小波系数为自变量，以植被参数为因变量，进行逐步回归，选择植被参数反演所需的光谱波段，构建植被参数定量反演的模型，求算建模R²；

4）确定最佳小波分解尺度：比较不同小波分解尺度所构建模型的建模R²，将建模R²最大的模型确定为最佳模型。

作为优选的一种方案：所述反演方法还包括以下步骤：

5）模型验证：利用独立样本的植被参数及光谱数据，对光谱数据进行小波变换，生成不同尺度的小波系数，将变换系数代入步骤3）所建立的回归模型，验证回归模型的监测效果。

    进一步，所述步骤2）中的小波含量包括53种小波：Haar 小波； Daubechies 小波系列：‘db2’、‘db3’、……‘db10’，共9种；Symlets小波系列： ‘sym2’ ‘sym3’、……‘sym8’，共7种；Coiflets小波系列：‘coif1’、 ‘coif2’、……‘coif5’，共5种；Biorthogonal小波系列：‘bior1.1’、 ‘bior1.3’、……‘bior6.8’，共15种； Reverse biorthogonal 小波系列：‘rbio1.1’、‘rbio1.3’、……‘rbio6.8’，共15种； 以及 Discrete Meyerpseudo 小波。

再进一步，所述步骤2）中的植被参数包括生物物理参数，所述生物物理参数包括生物量和叶面积指数。

更进一步，所述植被参数还包括生物化学参数，所述生物化学参数包括叶绿素和氮素。

本发明的技术构思为：植被参数具有不同的吸收反射特点，在光谱上表现为不同的吸收反射峰，这些反射特征具有一定的尺度特性。小波分析可以将植被叶片或是冠层的反射光谱在不同尺度上分解，因此可以利用适合与目标参数的尺度光谱信息反演植被参数，以提高植被参数的反演精度。

本发明的有益效果主要表现在：（1）、能够显著的提高植被参数高光谱遥感反演精度，主要是因为小波分析可以将光谱在不同尺度上分析，利用合适尺度的光谱信息更能表现植被参数的吸收特征；（2）、针对不同植被参数具有广泛的适用性，即可以根据植被参数自身的吸收、反射特征自动确定合适的分解尺度；（3）、不仅适用于叶片或是冠层反射光谱，也适用于卫星遥感高光谱数据。

附图说明

图1是基于反射光谱小波变换的植被参数遥感反演方法的流程图；

图2是图2 小波变换的尺度对建模效果的影响。

具体实施方式

下面结合附图对本发明做进一步的说明。

参照图1和图2，一种基于反射光谱小波变换的植被参数遥感反演方法，所述反演方法包括以下步骤：

1）光谱数据变换：获取不同条件的植被参数及其光谱，对光谱进行各种光谱变换，包括导数变换、伪吸收变化等；