GF-6 WFV光谱植被指数模型的叶面积指数反演

王枭轩，卢小平*，孟庆岩，李国清，王 俊，张琳琳，杨泽楠

1. 河南理工大学自然资源部矿山时空信息与生态修复重点实验室，河南 焦作 454000 2. 中国科学院空天信息创新研究院，北京 100101 3. 三亚中科遥感研究所，海南 三亚 572029 4. 河南省遥感测绘院，河南 郑州 450000

引 言

叶面积指数(LAI)是植被冠层结构最重要的生物物理参数之一[1]，它通常定义为单位地表面积上所有叶片面积之和的一半，它是反映植被个体特征和群体特征长势的关键指标，控制着地表植被的许多生物及物理过程[2]，因此，LAI的准确估算对作物长势、 产量预估和环境监测等具有重要的意义。

目前，国内外多数学者对LAI反演主要采用物理模型和传统光学植被指数模型两种方法。 其中，物理模型包含辐射传输的物理模型和计算机模型，辐射传输模型中具有代表性的是SUITA模型、 SAIL模型、 PROSPECT模型和PROSAIL模型等。 Jacquemoud等[3]通过分析冠层生物物理和辐射传输模型，对PROSAIL模型的历史进程进行了分析。 Duan等[4]基于高光谱无人机遥感影像数据，采用PROSAIL辐射传输模型结合LUT查表的方法对玉米、 马铃薯和向日葵三种作物叶面积指数进行反演； Li等[5]利用基于PROSAIL的检索表，GF-1遥感影像数据，为冬小麦品种特征和分布提供了先验知识，并利用光谱反射率结合15种植被指数反演得到不同物候期的叶面积指数，研究表明： 在不同生育期选择合适的反演模型和策略可以大大提高反演效果。 计算机模型具有代表性的如神经网络模型、 极限学习机模型和随机森林模型等，如刘俊等[6]采用BP神经网络模型反演了玉米LAI，证明了BP神经网络反演LAI的准确性； 贾洁琼等[7]利用随机森林结合PROSAIL模型反演LAI，证明了该模型能有效对农作物LAI进行监测； 雷宇斌等[8]采用极限学习机的方法反演道路植被LAI。 上述两种物理模型机理性强，使用的植被类型和空间范围更广，但模型参数较多且难获取，存在模型反演的病态问题。 传统光谱植被指数模型能够有效解决构造模型参数较多的问题，便于构建。 如Soudani等[9]基于IKONOS、 SPOT和ETM+等3种遥感影像数据结合5种植被指数模型研究温带针叶林和落叶林叶面积指数反演； Liu等[10]利用植被指数半经验模型综合分析了农作物在不同生育期的拟合精度； 贺佳等[11]通过不同植被指数回归模型反演冬小麦LAI； 赵虎等[12]利用5种植被指数回归模型对小麦LAI反演； 李军玲等[13]采用红边面积与蓝边面积的归一化指数反演冬小麦LAI。 上述传统光学植被指数的经验模型虽然能够及时有效、 无损害的监测作物长势，但机理性较弱，稳定性较差。

综上所述，虽然大部分国内外学者采用传统光谱植被指数模型或物理模型反演叶面积指数，但物理模型有反演病态的问题，传统光学植被指数模型会产生农作物机理的问题。 针对物理模型反演易产生病态问题，本工作基于我国国产自主卫星GF-6 WFV遥感影像，采用光谱植被指数模型反演乳塾期玉米LAI，同时引入红边波段参与建模，以增强模型稳定性。 红边波段范围在680～760 nm，是介于红波段到近红外波段，能有效指示农作物生长状况的敏感性光谱波段。

1 实验部分

1.1 实验区概述

实验区位于冀中平原西部的河北省石家庄市东南部栾城县(37°47′34″—37°59′20″N，140°28′36″—114°47′35″E)，地势平坦，整体由西北向东南倾斜，平原为主要地形，如图1所示。 这里气候温和，年平均气温12.8 ℃，年降水量423.7 mm，属温带大陆性季风气候，光热充足，雨热同季，有利于农作物生长。

图1 实验区采样点图Fig.1 Sampling sits of the experimental area

1.2 影像数据及预处理

高分6号多光谱遥感影像(GF-6 WFV)来源于中国资源卫星应用中心网，过境时间2019年10月08日与LAI采集时间基本一致，正好满足乳熟期玉米的时节，能更加有效的和准确的反应作物长势。 其包含8个波段，空间分辨率16 m，宽幅800 km，重访日期2 d，如表1所示。 同时对影像进行辐射定标，大气校正和正射校正。 辐射定标采用中国资源卫星应用中心下载的绝对定标系数，大气校正采用ENVI 5.3中的FLAASH大气校正模块，为获得地表真实反射率，正射校正采用DEM高程数据进行，使影像更加精确、 清晰。

1.3 实测数据

采用美国LI-COR LAI-2200C植物冠层分析仪，采集日期为2019年9月22日—29日上午07：30—09：30(基本与高分6号卫星过境时间一致) ，选取玉米样方的数量为50个，尺寸为16 m×16 m，并在其中挑选长势相同的3个点，所有样方都应当远离树木及建筑，还有道路，最少要保持30 m的距离，采样的时候需遮挡阳光，采用LAI-2000测量，具体来讲就是利用180°遮盖帽，之所以如此是为了确保测量不出现误差。 不仅如此，还要调整好实测地面和仪器之间的距离，按照标准应当保持5 cm，目的是确保玉米冠层和仪器的

表1 GF-6卫星WFV遥感影像Table 1 Remote sensing images of GF-6 WFV satellite

表2 光谱植被指数计算公式Table 2 Spectral vegetation indices

视角范围之间不产生误差。 为了防止叶面积指数测量误差，对各个样点的叶面积指数测量4次，取平均值。 同时，利用美国GPS定位仪实时记录每个玉米样点的经纬度。

1.4 方法

1.4.1 光学植被指数

光学植被指数就是根据地面实测的叶面积指数与光学遥感中的光谱数据或者变换行式(植被指数)建立回归分析模型。 采用这些植被指数可以用来估算叶面积指数，并监测植被长势。 通过选取6种典型抗饱和、 抗土壤和大气光学植被指数与实测叶面积指数建立回归模型。 同时考虑到GF-6 WFV遥感影像包含两个红边波段能更好的反演LAI，因此将其融入上述6种光学植被指数模型，构造出新的光学植被指数，植被指数计算公式如表2所示。

1.4.2 精度评定指数

通过精度评定可以有效的评价出植被指数与实测叶面积指数的拟合状况，得到最优的反演模型，即该模型反演的影像值和与其对应的实测值相一致。 选取决定系数(coefficient of determination，R2)和均方根误差(root mean square error，RMSE)作为评价指标，如式(1)和式(2)所示

(1)

式(1)中，y为实测生态参数值，yi为估算生态参数值，yc为实测生态参数平均值。

(2)

式(2)中，yi为预测值，yc为实测值，n为样本个数。

2 结果与讨论

2.1 相关性纷析

在拟合回归分析前，首先对光谱指数与LAI进行相关性分析，如表3所示。 根据分析结果可知，总体上20种光谱植被指数与LAI具有显著相关性，相关系数在0.4以上。 其中，具有红边参与构造的光谱指数相关性要高于无红边参与构造的光谱指数，尤其是NDSI3的相关性最高，达到0.881 0。

表3 光谱植被指数与乳熟期玉米LAI的相关系数Table 3 Correlation coefficient types between spectralvegetation index and LAI of maize at milk stage

2.2 LAI反演拟合模型分析

基于GF-6 WFV卫星遥感影像采用2.1节中提到的20种光谱指数与实测LAI建立光谱植被指数模型，如图2所示。

图2 光谱植被指数模型验证Fig.2 Verification of spectral vegetation index model

从图2可知，基于GF-6 WFV遥感影像构造的光谱植被指数与LAI呈现显著正相关，R2在0.51以上，说明整体拟合效果较好。 其中NDS13的拟合精度是最高的，R2为0.803，RMSE为0.301 2。 主要是由于该模型有两种红边波段构成，对LAI的解释能力较优。 综上所述，由红边参与构造的光谱植被指数与LAI拟合精度要高于未用红边参与构造的光谱植被指数，有利于叶面积指数的反演。

2.3 LAI反演空间分析

LAI对农作物长势分析、 产量预估和环境监测有着巨大的作用。 利用NDSI3光谱植被指数模型反演乳熟期玉米LAI，如图3所示。

从图3可知，从数值上分析，通过模型估算出来的LAI集中分布在2～3之间，接近于实测LAI值。 从空间分布上分析，玉米LAI的高值区主要分布在中部和南部地区，低值区分布在西北部，符合当地的实际情况。 但是东南部存在异常点，主要是受到城区道路和建筑物的影响，受到混合像元的影响。

图3 NDSI3光谱植被指数模型反演乳熟期LAIFig.3 NDSI3 spectral vegetation index model inversionfor LAI of corn at milk stage

2.4 模型验证

为了验证模型的准确性，将构建模型以外的20个实测LAI点与基于GF-6 WFV影像NDSI3模型反演的LAI进行回归分析，同时采用R2和RMSE进行精度评价，如图4所示。

图4 NDSI3模型验证Fig.4 NDSI3 model validation

从图4可知，实测LAI与NDSI3模型反演的LAI整体拟合较好，R2为0.804，RMSE为0.312 5，说明该模型能有效反演乳熟期玉米的生长状况。

2.5 MODIS LAI 验证

为了更加准确的验证模型的实用性，采用NDSI3模型反演得到乳熟期LAI的影像重采样到1 km，并随机抽取100个点。 将GF-6 WFV影像反演得到的叶面积指数值(LAIGF-6)与对应点的MODIS影像产品MOD15A2(LAIMODIS)进行回归拟合，并分析两者的相关性，如图5所示。

图5 交叉验证Fig.5 Cross validation

从图5可知，乳熟期玉米LAIGF-6与LAIMODIS有较好的一致性，R2为0.676 5，RMSE为0.248 6。 同时，对这100个点进行均值分析可知，LAIGF-6均值为2.1，LAIMODIS均值为2.9，说明定量化差异对本次实验影响较小，LAIMODIS均值要高于LAIGF-6，主要是由于MODIS影像像元混合严重和空间分辨率低有关。

3 结 论

利用实测LAI与GF-6遥感影像构建了5种光学植被指数和15种光谱植被指数对乳熟期玉米LAI进行了反演分析，得出以下结论：

(1)15种光谱植被指数和5种光学植被指数都与LAI达到了显著相关，红边参与构造的15种光谱植被指数相关性高于5种光学植被指数，其中归一化差值光谱指数3与乳熟期玉米LAI相关性最高。

(2) 通过拟合回归可以得出，15种光谱植被指数与乳熟期玉米LAI较5种光学植被指数拟合精度高。 通过反演结果图可知，并未出现过拟合的现象。 因此，可以有效证明归一化差值光谱指数3能有效反演乳熟期玉米LAI。

(3) 通过MODIS LAI产品验证可知，在反演精度要求不高时，MODIS产品可以有效反演玉米LAI，监测玉米长势。 当反演精度要求较高时，GF-6遥感影像可有效反演玉米LAI，成为监测玉米长势的重要数据源。

综上所述，由GF-6 WFV构建的归一化差值光谱指数3能有效反演乳熟期玉米LAI，为玉米长势监测提供依据。 但由于受影像获取时间和野外观测条件的限制，本研究只是对乳熟期玉米的叶面积指数进行了研究，接下来应尽可能从玉米全生育期的角度，分析该方法的适用性和精度验证。