冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型_艾金泉

第35卷第4期

2015年2月

生态学报ACTA ECOLOGICA SINICA Vol．35，No．4Feb．，2015

http ：//www．ecologica．cn 基金项目：福建省科技计划项目（2010Ｒ1037）；福建师范大学研究生科技创新基金资助

收稿日期：2013-04-27；网络出版日期：2014-04-11

*通讯作者Corresponding author．E-mail ：whchenfz@sohu．com

DOI ：10．5846/stxb201304270834艾金泉，陈文惠，陈丽娟，张永贺，周毅军，郭啸川，褚武道．冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型．生态学报，

2015，35（4）：1175-1186．

Ai J Q ，Chen W H ，Chen L J ，Zhang Y H ，Zhou Y J ，Guo X H ，Chu W D．Hyperspectral remote sensing estimation models for foliar photosynthetic pigment contents at canopy level in an invasive species ，Spartina alterniflora ．Acta Ecologica Sinica ，2015，35（4）：1175-1186．

冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型

艾金泉，陈文惠*，陈丽娟，张永贺，周毅军，郭啸川，褚武道

福建师范大学地理科学学院，福州350007

摘要：以闽江河口鳝鱼滩湿地互花米草（Spartina alterniflora ）的实测冠层高光谱反射率和叶片光合色素含量（LPPC ）为数据源，

在分析LPPC 与原始光谱反射率、一阶导数光谱反射率、22种已报道光谱指数和14种新构建的植被指数相关性的基础上，利用

直线回归、指数回归、对数回归以及乘幂回归方法，系统地比较了36种植被指数在估算互花米草LPPC 中的表现。研究表明：

（1）一阶导数光谱反射率组合的植被指数用于估算互花米草的LPPC 优于原始光谱反射率；（2）红边区域一阶导数光谱是估测

互花米草LPPC 的最佳波段；（3）对于单一色素含量的估算，叶绿素a （Chla ）的最佳估算指数为FDNDVI ［

723，703］；叶绿素b （Chlb ）的最佳估算指数为FDＲVI ［723，525］；类胡萝卜素（Cars ）的最佳估算指数为FDNDVI ［723，703］；（4）对于使用统一参量

同时估算Chla 、Chlb 、Cars ，由FDＲVI ［723，703］建立的对数估算模型效果最佳。研究成果可为湿地植物生化参量反演提供参

考，也可为闽江河口湿地入侵种互花米草的动态监测和生态评估管理提供有力的科学依据。

关键词：光合色素；敏感波段；植被指数；红边；湿地植物

Hyperspectral remote sensing estimation models for foliar photosynthetic pigment contents at canopy level in an invasive species ，Spartina alterniflora

AI Jinquan ，CHEN Wenhui *CHEN Lijuan ，ZHANG Yonghe ，ZHOU Yijun ，GUO Xiaochuan ，CHU Wudao College of Geographical Science ，Fujian Normal University ，Fuzhou 350007，China

Abstract ：Foliar photosynthetic pigments are the most important biochemical parameters relative to the physiological function of wetland plants．Quantitative estimation of photosynthetic pigments can provide important information about the dynamics of the vegetation productivity ，vegetation stress ，or nutrient cycles within wetland ecosystems．However ，the estimation of foliar photosynthetic pigments is complicated because canopy reflectance in the visible and near infrared wavelengths is affected by confounding effects that come not only from foliar photosynthetic pigment contents variation but also from the changing environmental conditions of wetland ecosystem．Our objective was to address the question by establishing hyperspectral remote sensing estimation models for foliar photosynthetic pigments at canopy level in an invasive species ，Spartina alterniflora ．In this study ，the hyperspectral reflectance of canopy leaf and leaf photosynthetic pigment contents （LPPC ）from S．alterniflora in Min Ｒiver Shanyu beach were recorded as data source．The correlation between LPPC and raw spectral reflectance ，the first derivative reflectance ，22reported vegetable indices and 14new formed indices were determined．Based on the results of correlation analysis ，a total of 36indices were tested by linear regression ，exponential regression ，logarithm regression and the power of regression to explore their potentials in LPPC estimation in

6711生态学报35卷

S．alterniflora．The results showed that：（1）This study selected5wavebands in the region of400—900nm，which appeared to be the optimal wavebands for the S．alterniflora foliar photosynthetic pigments estimation．Of the selected wavebands，the most frequently occurring wavebands were723nm，703nm，525nm，752nm，668nm．（2）Vegetation index portfolio by the first derivative reflectance was evidently better than raw reflectance for estimating LPCC in S．alterniflora．Combining the optimal wavebands，the results indicated that the first derivative of reflectance in the red edge region（680—760nm）was the optimal band for estimating LPCC．（3）For a single pigment content，the best indexes for estimating chlorophyll a（Chla），chlorophyll b（Chlb）carotenoids（Cars）were FDNDVI［723，703］，FDＲVI［723，525］，and FDNDVI［723，703］，respectively．The three new formed indices were proved to have better linearity with corresponding photosynthetic pigment．（4）Using the same index to estimate different pigments，the best model was the logarithmic model using FDＲVI［723，703］，with high predicted correlation coefficientsＲ2of0.6997，0.7187，and 0.7132，respectively．This study would not only provide a good reference for hyperspectral remote sensing retrieval of biochemical variables in wetland vegetation，but also provide a strong scientific basis for the dynamic monitoring of S．alterniflora and management of ecological assessment in MinＲiver estuary．

Key Words：photosynthetic pigment；sensitive bands；vegetation index；red edge；wetland vegetation

叶片光合色素含量广泛用于诊断植被生长状况和健康状况，是植被环境胁迫、光合作用能力和植物发育阶段的指示器［1］。叶片光合色素含量与植被独特的光谱特征有密切的联系，在可见光波段内，叶片光谱特征主要受光合色素含量的影响，在近红外波段主要受叶片内部结构、叶片含水量及叶面积指数等影响；而高光谱遥感技术能通过微弱的光谱响应光合色素含量变化，因而，可用高光谱遥感估算叶片光合色素含量。王福民等［2］对水稻色素含量的估算研究表明，红边是估测水稻最敏感的波段。Wu等［3］基于高光谱植被指数建立了冬小麦叶绿素含量的高光谱估算模型，发现改进型的两种综合植被指数MCAＲI/OSAVI［705，750］和TCAＲI/ OSAVI［705，750］最适合叶绿素含量反演。Main等［4］研究了73种已见报道估算叶片叶绿素含量的高光谱指数的鲁棒性，结果表明基于红边位置导数光谱和红边拐点的导数光谱建立的指数用于估算叶片叶绿素含量适用性和稳定性最好。Sims等［5］比较常用的高光谱指数在不同植被类型、不同的叶片冠层结构反演叶片叶绿素含量中的表现，发现简单比值植被指数mSＲ705和归一化差值植被指数mND705与叶绿素含量相关性最好。这些研究表明，利用高光谱遥感技术估算陆地植被叶片光合色素含量是可行可靠的，只是对于不同的植物所用的最佳波段和最佳估算参量不同。

湿地植被是湿地生态系统最重要的组成部分，是衡量湿地生态系统健康的关键指标，对区域环境的稳定性具有很强的指示器作用［6］。湿地植被光合色素含量状况是湿地植被光谱变异最重要的影响因子，使用高光谱遥感无损估算其光合色素含量有利于了解其光谱特征变异机理、生理状况、物种组成和生物化学循环的研究［7］。但是，研究者也注意到，由于湿地生物多样性高、水文条件复杂等特殊的生态环境，给高光谱遥感技术应用于湿地植被生化参数的估测带来了诸多挑战［7］。很幸运的是，通过文献查找发现，仍有少数研究者已经开展了高光谱遥感在天然湿地植被光合色素含量的估算研究。卢霞等［8］研究了大米草（Spartina anglica）室内叶片光谱特征参数与叶绿素浓度关系，发现利用水分指数（WI）基于三次函数构建的叶绿素浓度的高光谱检测模型精度最高。李凤秀等［9］采用可见光-近红外波段对冠层水平湿地小叶章（Calamagrostis angustifolia）叶片Chla含量建立了高光谱遥感估算模型，发现比值微分光谱指数FＲVI［548，556］建立的Chla 的线性预测模型表现最佳。Ge等［10］对湿地植物入侵种芦竹（Arundo donax）冠层叶片的光合色素含量估算研究表明，693nm波段是色素含量最敏感的波段。这表明高光谱遥感技术用于估算天然湿地植被叶片光合色素含量成为可能。

自1979年以来，互花米草作为外来种被引入福建省沿海滩涂，它迅速扩张，已对当地生态安全产生了威胁。如何有效的对它的长势进行动态监测，具有重要的意义。本文以闽江河口鳝鱼滩湿地入侵物种互花米草

http：//www．ecologica．cn

http ：//www．ecologica．cn 为研究对象，选定闽江河口湿地生态定位观测站200m 附近的小区域作为实验区，以野外实测数据作为数据源，利用导数光谱技术和相关性分析法，确定互花米草光合色素含量的敏感波段，构建新的植被指数，并同已

报道的植被指数进行比较，

旨在确定具有适用性和可靠性强的光合色素含量估算模型。本研究可为入侵种互花米草的光谱特征变异机理、

长势动态监测、入侵机制和生态管理评估提供参考依据。1

材料与方法图1研究区地理位置与采样点Fig．1The location of the study area and sampling site

1．1研究区概况

实验区位于闽江河口鳝鱼滩湿地生态定位观测站

附近（图1），

地理坐标为26?00'36?—26?03'42?N ，119?34'12?—119?41'40?E ，地处中亚热带和南亚热带过渡

区。区内气候温暖湿润，

年均气温19．7?，年均降水量1380mm 左右，年均降水日数为153d ，受正规半日潮影

响［11］。该湿地是闽江河口面积最大和保护较好的一处

天然湿地，

湿地植物以芦苇（Phragmites australis ）、短叶茳芏（Cyperus malaccensis ）和互花米草为主，其中，互花

米草是外来入侵种。近年来互花米草蔓延严重，

2002年还是零星分布，到2008年已经超过200hm 2，给该湿

地生态系统造成严重的生态威胁

［12］。1．2

数据获取1．2．1冠层光谱采集冠层光谱反射率的测定采用美国Analytical Spectral Device （ASD ）公司生产的便携式野外地物光谱仪FieldSpec3，该光谱仪波段范围在350—2500nm ，其中，350—1000nm 光谱采样间隔为1．4nm ，光谱分辨率为3nm ；1000—2500nm 光谱采样间隔为2nm ，光谱分辨率为10nm 。在光谱输出时，仪器自动进行重采样（重采样间隔为1nm ），共输出2151个波段数。

本研究采样点位于鳝鱼滩观鸟台以西200m 范围左右互花米草入侵严重区域（图1），采样日期是2012年10月20日和2012年10月22日（8：00—15：00属于低潮期，植被不受潮汐影响，地面无积水），风力小于3级，无云，天气晴朗，时间是地方时10：00—13：30。总共选取36个样方，各个样方相隔10m 左右，所有样方互

花米草长势较好（没有其他植物），植株高约180—220cm ，植株密度约200—300株/m 2，覆盖度大于90%，叶

片颜色呈淡绿色，冠层结构接近呈平直型。传感器探头使用25?视场角，高度保持在冠层上方1．0m 附近。每个样方的大小为1m ?1m ，每次采集10条光谱曲线，且每组数据采集前都用白板对仪器进行一次校正，以

消除环境变化带来的影响，

总共采集36个有效样本360条冠层光谱曲线。数据预处理采用ViewSpePro 软件，对采集的冠层光谱曲线进行平滑、去噪、一阶导数处理，并以每个采样点10条光谱曲线的平均值作为该采样点的冠层反射率。考虑到湿地土壤含水量高，易对中红外波段产生较大的噪音和借鉴相关研究经验

［8，10］，本研究只选取400—900nm 之间的波段进行分析。

1．2．2光合色素含量测定野外光谱测量完后，采叶人员同时采集相应样本点冠层上部、中部、下部成熟叶片6—8叶，装入自封袋，标号，直接保存在放有冰袋的保鲜箱中，采样同一天拿回实验室进行光合色素含量测定。具体测定方法如下［13］，先将采样叶片全部剪碎均匀混合，称取剪碎混合后的新鲜叶片3份，每份0.200g ，分别放入研钵中，加入少量石英砂和碳酸钙粉及2—3mL 96%的乙醇研磨成匀浆，再加入乙醇10mL 继续研磨至组织变白，静置3—5min ，过滤到25mL 棕色容量瓶中，再用UV2450-紫外分光光度计测量665nm 、

649nm 下的光密度OD 值，然后根据下面公式计算叶绿素含量浓度：

77114期艾金泉等：冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型

Chla（mg/L）=13．95?OD

665－6．88?OD

649

（1）

Chlb（mg/L）=24．96?OD

649－7．32?OD

665

（2）

Cars（mg/L）=（1000?OD

470

－2．05?Chla－114．8?Chlb）/245（3）X（mg/g）=浓度（mg/L）?提取液体积（mL）/质量（g）/1000（4）式中，OD665为665nm处的吸光度；OD649为649nm处的吸光度；OD665为665nm处的吸光度。

1．3本文用到的高光谱植被指数

总结和筛选已报道用于冠层水平绿色植物叶片光合色素含量估算的高光谱指数，选择其中常用的22种分析（表1前22个指数）。所选指数主要有以下几类：比值植被指数，如ＲVI750/550；归一化植被指数，如

NDVI；综合植被指数，如MCAＲI/OSAVI；微分比值植被指数DＲ

715/DＲ

705

；微分归一化植被指数BmSＲ

705

。这

些已报道的高光谱指数中，所用到的估测光合色素含量的敏感波段集中在445、550、670、700、705、730、760、800nm附近，这为提出新的高光谱植被指数提供了参考。

本文充分借鉴前人研究经验，基于光谱特征分析和原始光谱、一阶导数光谱与光合色素含量间的相关性分析，选取得到的所有峰值点和拐点（图1，图2），通过波段组合，构建新的光谱指数并选取其中14个表现最好的指数纳入本文分析（表1中后14个指数）。

表1本文用于估算互花米草叶片光合色素含量的高光谱指数

Table1Hyperspectral indices used for estimation foliar pigment content of S．alterniflora in this study

光谱指数Spectral indices 计算公式

Algorithm formula

参考文献

Ｒeferences

编号

Eq．ID

ＲVI750/550Ｒ750/Ｒ550［14］（5）ＲVI750/700Ｒ750/Ｒ700［15］（6）Ｒ762/Ｒ734Ｒ762/Ｒ734［16］（7）

MSＲ［（Ｒ

750/Ｒ705）－1］/（Ｒ750/Ｒ705）+

槡1［17］（8）

NDVI（Ｒ800-Ｒ670）/（Ｒ800+Ｒ670）［15，18］（9）PＲI（Ｒ531-Ｒ570）/（Ｒ531+Ｒ570）［19］（10）NDVI705（Ｒ750-Ｒ705）/（Ｒ750+Ｒ705）［17，20］（11）TVI NDVI+

槡0.5［21］（12）mND705（Ｒ750-Ｒ705）/（Ｒ750+Ｒ705－2Ｒ445）［5］（13）Green-NDVI（Ｒ800-Ｒ550）/（Ｒ800+Ｒ550）［22］（14）SIPI（Ｒ800-Ｒ445）/（Ｒ800-Ｒ680）［23］（15）

CAＲI CAＲI=

|670a+Ｒ670+b|Ｒ700

Ｒ670a2+

槡1

a=（Ｒ700-Ｒ500）/150；b=Ｒ550－550a

［24］（16）

EVI2．5?（（Ｒ800-Ｒ670）/（Ｒ800－（6?Ｒ670）－（7．5?Ｒ475）+1））［25］（17）MCAＲI（（Ｒ700-Ｒ670）－0.2?（Ｒ700-Ｒ550））?（Ｒ700/Ｒ670）［26］（18）

MCAＲI/OSAVI （（Ｒ700-Ｒ670）－0.2?（Ｒ700-Ｒ550））?（Ｒ700/Ｒ670）

（1+0.16）?（Ｒ800-Ｒ670）/（Ｒ800+Ｒ670+0.16）

［26］（19）

Boochs A D703［27］（20）DＲ715/DＲ705DＲ715/DＲ705［16，28］（21）BmSＲ705（DＲ750－DＲ445）/（DＲ750+DＲ445）［29］（22）ＲEIP红边内最大一阶导数波长［15］（23）FDNDVI［750，550］（DＲ750－DＲ550）/（DＲ750+DＲ550）［30］（24）BmND705（DＲ750－DＲ445）/（DＲ750+DＲ705－2DＲ445）［29］（25）SＲE/SBE红边内最大一阶微分与蓝边内最大一阶微分的比值［31］（26）MCCV全波段中原始光谱反射率与光合色素含量相关性最大波长处的反射率值本文（27）

MCCFDV 全波段中一阶导数光谱反射率与光合色素含量相关性最大波长处的一阶导

数值

本文（28）

8711生态学报35卷

http：//www．ecologica．cn

续表

光谱指数Spectral indices 计算公式

Algorithm formula

参考文献

Ｒeferences

编号

Eq．ID

ＲVI668/723Ｒ668/Ｒ723本文（29）

ＲVI723/703Ｒ723/Ｒ703本文（30）

ＲVI752/703Ｒ752/Ｒ703本文（31）

ＲVI723/525Ｒ723/Ｒ525本文（32）

FDＲVI［723，668］DＲ723/DＲ668本文（33）

FDＲVI［723，703］DＲ723/DＲ703本文（34）

FDＲVI［752，703］DＲ752/DＲ703本文（35）

FDＲVI［723，525］DＲ723/DＲ525本文（36）FDNDVI［723，668］（DＲ723－DＲ668）/（DＲ723+DＲ668）本文（37）FDNDVI［723，703］（DＲ723－DＲ703）/（DＲ723+DＲ703）本文（38）FDNDVI［752，703］（DＲ752－DＲ703）/（DＲ752+DＲ703）本文（39）FDNDVI［525，723］（DＲ525－DＲ723）/（DＲ525+DＲ723）本文（40）Ｒx指波长x nm处原始光谱反射率；DＲx指波长x nm处一阶导数光谱反射率

1．4模型构建、评价与验证

采用交叉验证法对数据建模，从所有的有效样本中，随机选取24个样本作为建模数据集，其他12个作为验证数据集。利用建模数据集，以各高光谱植被指数为自变量x，互花米草光合色素含量为因变量y建立估测模型。对同一组数据，建模过程中，选取直线回归、指数回归、对数回归以及乘幂回归中决定系数Ｒ2最大的模型进行建模，并相应给出模型的决定系数（C-Ｒ2）和均方根误差（ＲMSEC）作为精度评价标准。

模型建立后，利用验证数据集对其预测能力进行检验，并相应给出验证模型的决定系数（P-Ｒ2）、均方根误差（ＲMSEP）和斜率（Slope）作为精度评价标准。最终模型的确定取决于相对较高的决定系数（C-Ｒ2、P-Ｒ2）和相对较低的均方根误差（ＲMSEC、ＲMSEP）；此外结合验证模型的Slope的大小进行精度检验，Slope越接近1，模型的精度就越高［32-33］。其中，

ＲMSE=∑（y m－y c）2

槡/n（41）式中，y m是光合色素含量的实测值，y c是模型的估测值，n为样本的数量。

2结果与分析

2．1互花米草光合色素统计分析

为了减少误差和奇异值对模型拟合的影响，先对冠层水平互花米草叶片光合色素测定数据进行统计描述分析和奇异值去除，奇异值的判定方法见文献［34］。奇异值删除后，叶片光合色素与光谱配套的数据样本总数为36。随机选取24个样本作为建模数据集，其他12个作为验证数据集，各样本光合色素统计描述见表2。

表2互花米草叶片光合色素含量统计描述

Table2Statistic description of LPPC in S．alterniflora

色素Pigment 数据集

Datasets

样本数

Samples

平均值

Mean/（mg/g）

最大值

Max/（mg/g）

最小值

Min/（mg/g）

标准差

Standard

deviation

/（mg/g）

变异系数

Coefficient

of variation/%

Chla建模Calibration241．3502．4120.6820.4460.345验证validation121．3432．5770.7920.4720.351 Chlb建模240.4730.8950.2660.1550.328验证120.4510.8660.2960.1550.344 Cars建模240.2760.5000.1010.0970.351验证120.2750.4710.1410.0960.3499711

4期艾金泉等：冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型

http：//www．ecologica．cn

http ：//www．ecologica．cn 由表2可知，Chla 含量的最大值为2．577mg /g ，最小值为0.682mg /g ；Chlb 含量的最大值为0.895mg /g ，最小值为0.266mg /g ；Cars 含量的最大值为0.500mg /g ，最小值为0.101mg /g ；其光合色素含量有较大的

分布区间，

可满足建模要求。2．2互花米草光谱特征分析

图2是互花米草叶片原始光谱平均反射率与一阶导数光谱平均反射率曲线。从图2可知，互花米草原始光谱反射率曲线具有典型植物光谱特征。在绿光波段（520—600nm ），波长550nm 附近产生了一个明显反射

峰（绿峰），

主要是绿光对光合色素的强烈反射造成的；在红光波段，波长670nm 附近产生了一个较明显的吸收谷（红谷），主要是红光对叶绿素的吸收作用产生的。在近红外波段（700—900nm ），其光谱特征主要受叶

片结构和水分的影响，

最显著的特征是在700—800nm 之间有一反射的陡坡（红边）。一阶导数光谱反射率在525nm 和727nm 处出现两个极大值，代表光谱反射率增速最快的极点，

525nm 处反射率上升是由于叶绿素的反射作用增强引起的，

而727nm 处主要由于红边效应的作用（图2）

。图2

互花米草反射率光谱与一阶导数光谱曲线Fig．2Ｒeflected spectra and first derivative curves of S．alterniflora

2．3叶片光合色素含量与原始光谱及一阶导数光谱的相关性分析

图3显示互花米草叶片光合色素含量与原始光谱及一阶导数的相关性。由图3可知，Chla 、

Chlb 、Cars 与原始光谱的相关性变化趋势一致，都呈较强的负相关关系。在红至近红外波段（700—900nm ），原始光谱与

光合色素最大相关性出现在703nm 附近，Chla 、

Chlb 、Cars 与原始光谱达到了最大相关性值分别为－0.747、－0.682、－0.746；拐点出现在752nm 附近，从700到752nm 之间，相关系数值逐渐增大，而752nm 后相关系数值趋于稳定。在可见光波段（400—700nm ），互花米草叶片光合色素含量与一阶导数相关性波动性比较大，出现多个极值和拐点

。

图3

叶片光合色素与原始反射率及一阶导数光谱反射率的相关性（n =36）Fig．3Correlograms between LPPC and raw reflected spectral and its first derivative spectral

0811生态学报35卷

http ：//www．ecologica．cn

在绿光波段（520—600nm ），波长525nm 处是互花米草光合色素与一阶导数光谱的一个极值点，系蓝边的拐点产生；在红光波段（630—690nm ），Chla 、Chlb 、Cars 在668nm 处达到了最大的正相关关系，值分别为0.808，0.711，0.787，系红光吸收谷的位置；在691nm 处出达到了最大的负相关关系，值分别为－0.764，－0.665，－0.753；在红边区域（680—760nm ），723nm 附近出现拐点，波长小于723nm 时光谱一阶导数与叶片光合色素相关性达到0.01极显著相关水平的检验（r =0.4238，n =36），波长大于723nm 时没有通过检验（图3）。2．4

互花米草叶片光合色素含量与高光谱指数的相关性分析

通过互花米草叶片光合色素与高光谱植被指数之间的相关性分析（表3）可知，

所有光谱指数中，基于原始光谱组合的植被指数相对于单波段的MCCV 没有明显优势，而对于导数光谱所组合的植被指数相对于导数单波段的MCCFDV 优势明显，

并且基于原始光谱反射率组合的植被指数与光合色素含量间的相关系数r 整体低于基于一阶导数光谱反射率组合的植被指数，这说明一阶导数光谱反射率组合的植被指数更适合用于估算互花米草光合色素含量。其中，FDNDVI ［752，703］与Chla 相关系数r 达到最大，值为0.842；FDＲVI ［723，525］与Chlb 的相关性达到最大，值为0.894；FDNDVI ［723，525］与Cars 的相关性达到最大，值为0.835。

表3叶片光合色素含量与高光谱植被指数之间的相关性分析（n =36）

Table 3

Correlation coefficients between hyperspectral vegetable indices and LPPC

指数Indices 相关系数r Correlation coefficient r

Chla Chlb Cars 指数Indices 相关系数r Correlation coefficient r Chla Chlb Cars MCCV －0.747－0.682－0.746ＲVI 752/7030.7650.8240.723ＲVI 750/5500.7030.7680.668ＲVI 723/5250.4860.5340.487ＲVI 750/7000.7150.7920.674MCCFDV 0.8080.7110.787Ｒ762/Ｒ7340.7290.7730.683Boochs A －0.651－0.517－0.657MSＲ0.6980.7690.658DＲ715/DＲ7050.7000.7480.624NDVI 0.5670.4390.580BmSＲ7050.7440.7690.705PＲI －0.1240.099－0.102ＲEIP

0.5830.5700.566NDVI7050.6580.7690.698FDNDVI ［750，550］0.7570.8110.715TVI 0.5480.6240.534BmND 7050.8330.8580.809mND 7050.7760.8440.801SＲE /SBE 0.8310.8830.812Green-NDVI 0.6640.7190.635FDＲVI ［668，723］0.8270.8640.831SIPI 0.4410.5060.444FDＲVI ［723，703］0.8390.8890.826CAＲI 0.5510.5540.522FDＲVI ［

752，703］0.8310.8660.791EVI 0.5550.4470.581FDＲVI ［723，525］0.8410.8940.820MCAＲI －0.609－0.489－0.604FDNDVI ［723，668］－0.821－0.856－0.829MCAＲI /OSAVI －0.749－0.659－0.726FDNDVI ［723，703］0.8360.8800.832ＲVI 668/723－0.326－0.398－0.332FDNDVI ［

752，703］0.8420.86500.811ＲVI 723/703

0.631

0.702

0.602

FDNDVI ［525，723］

－0.840

－0.881

－0.835

2．5互花米草叶片光合色素含量的高光谱遥感估算模型构建与验证

挑选表3中至少与一种光合色素相关性r 值大于0.8的12个植被指数利用建模数据集构建模型，以C-

Ｒ2、ＲMSEC 为评价指标，结果如表4所示。表4可知，除mND 705、MCCFDV 外其他指数对互花米草光合色素的估算都取得了较好的效果。其中，对Chla 的建模精度最高的4个指数依次为FDＲVI ［752，703］、FDNDVI ［752，703］、BmND 705、FDNDVI ［723，703］；对Chlb 的建模精度最高的4个指数依次为FDＲVI ［723，525］、SＲE /SBE 、FDNDVI ［723，525］、FDＲVI ［723，703］；对Cars 的建模精度最高的4个指数依次为FDNDVI ［723，525］、SＲE /SBE 、FDNDVI ［723，703］、FDＲVI ［723，703］。

为了验证模型的适用性和稳健性，利用验证数据集对模型预测能力进行验证，以P-Ｒ2、ＲMSEP 、Slope 为1

8114期艾金泉等：冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型

评价指标，其结果如表4所示。对Chla的验证精度最高的4个指数依次为FDＲVI［723，525］、mND705、FDＲVI ［723，703］、FDNDVI［723，703］；对Chlb的验证精度最高的4个指数依次为FDＲVI［723，525］、MCCFDV、FDNDVI［723，703］、FDＲVI［723，703］；对Cars的验证精度最高的4个指数依次为FDＲVI［723，703］、FDNDVI ［723，703］、mND

705

、FDＲVI［723，525］。

表4叶片光合色素含量估算模型及评价指标

Table4LPPC estimation models and evaluating indicator

指数Indices 光合色素

Pigment

模型函数

Model function

建模Calibration

决定系数

C-Ｒ2

均方根误差

ＲMSEC

验证Validation

决定系数

P-Ｒ2

均方根误差

ＲMSEP

斜率

Slope

mND705Chla y=5．871x－2．08540.58500.28150.71040.28180.4839 Chlb y=2．3211x－0.88540.75530.07520.70810.09610.5818

Cars y=1．2957x－0.48260.61130.05880.70210.05130.5508 MCCFDV Chla y=2．6972e3035．8x0.73700.23640.70090.28520.4663 Chlb y=0.7971e2032x0.50710.10530.74360.09660.3913

Cars y=0.5096e2414．3x0.50150.07120.57180.06690.3844 BmND705Chla y=2．6476x+2．11730.75730.21530.63430.30060.4575 Chlb y=0.9574x+0.75020.81810.06480.65450.10210.5118

Cars y=0.5487x+0.43460.69800.05190.61490.06170.4589

SＲE/SBE Chla y=1．8671ln（x）－1．81240.73880.22330.66460.28650.5125 Chlb y=0.6992ln（x）－0.71150.85590.05770.67120.09390.5874

Cars y=0.3972ln（x）－0.39710.71720.05020.65020.05850.5301

FDＲVI［723，703］Chla y=2．5189ln（x）+0.36220.73760.22380.70790.27990.4867 Chlb y=0.9384ln（x）+0.10470.84580.05970.73010.09380.5601

Cars y=0.5337ln（x）+0.06630.71050.05080.71770.05560.5098

FDＲVI［668，723］Chla y=16．59x+2．74810.69920.26870.67530.23970.5619 Chlb y=6．179x－0.99350.80130.09280.66820.06780.6341

Cars y=3．5859x－0.57780.70090.05310.68690.05160.6021

FDＲVI［752，703］Chla y=1．337ln（x）+2．2380.77360.20790.63280.29560.4744 Chlb y=0.4785ln（x）+0.79060.81860.06470.64940.10050.5238

Cars y=0.2763ln（x）+0.45910.70910.05090.60590.06120.4715

FDＲVI［723，525］Chla y=0.36x+0.59370.70770.22460.72950.26190.5417 Chlb y=0.1363x+0.26370.85790.05220.80580.08090.5988

Cars y=0.0706x+0.13470.67660.05070.69920.05420.5533 FDNDVI［723，668］Chla y=－8．2384ln（x）+2．74250.69890.23980.67330.26930.5617 Chlb y=－3．0685ln（x）+0.99140.80100.06780.66590.09310.6336

Cars y=－1．7812ln（x）+0.57660.70090.05160.68560.05320.6024 FDNDVI［723，703］Chla y=5．3318x+0.32370.73790.22370.69970.28260.4869 Chlb y=1．9849x+0.09070.84480.05980.71870.09530.5596

Cars y=1．1311x+0.05780.71260.05060.71320.05590.5127 FDNDVI［752，703］Chla y=2．8975x+2．26340.76550.21160.65800.28940.4726 Chlb y=1．0408x+0.80080.81590.06520.68020.09680.5258

Cars y=0.5966x+0.46360.69630.05190.61580.06080.4628 FDNDVI［723，525］Chla y=－8．7581x－4．61240.73790.22390.60370.34250.4996 Chlb y=－3．2737x－1．7560.85410.05810.62460.12250.5829

Cars y=－1．8677x－0.9960.72190.04970.60580.06920.5173

综合比较，虽然已见报道的指数中BmND705、SＲE/SBE等模型（图4）对互花米草的光合色素的建模精度很高，但其验证模型的精度都相对较低；而mND705预测精度很高，但其建模模型精度却较低，说明已见报道的2811生态学报35卷

http：//www．ecologica．cn

http ：//www．ecologica．cn 植被指数建立的互花米草光合色素模型的稳定性不强，不是最佳估算模型。为了找到更好的估算模型，本文选取建模精度和预测精度排名都能在所有估算模型前4名的模型作为最佳模型，这样既保证了模型的精度又

使模型的可靠性更强。发现对于单一色素含量估算而言，

Chla 和Cars 的最佳估算指数为FDNDVI ［723，703］，Chlb 的最佳估算指数为FDＲVI ［723，525］，其直观形状如图4所示。其中FDNDVI ［

723，703］建立的Chla 估算模型的C-Ｒ2、ＲMSEC 、P-Ｒ2、ＲMSEP 、Slope 分别为0.7379、0.2237、0.6997、0.2826、0.4869；FDＲVI ［723，525］建立的Chlb 估算模型的C-Ｒ2、ＲMSEC 、P-Ｒ2、ＲMSEP 、Slope 分别为0.8579、0.0522、0.8058、0.0809、0.5988；FDNDVI ［723，703］建立的Cars 估算模型的C-Ｒ2、ＲMSEC 、P-Ｒ2、ＲMSEP 、Slope 分别为0.7126、0.0506、0.7132、0.0559、0.5127。

整体上看，如果只用一个植被指数来估算3种不同的光合色素含量，那么以FDＲVI ［723，703］建立的对

数估算模型效果最佳，它是所有指数中唯一的对Chla 、

Chlb 、Cars 估算精度C-Ｒ2、P-Ｒ2均能达到0.7以上的指数，

且其ＲMSEC 、ＲMSEP 、Slope 精度排名位于所有模型前4以上（图4）。结合单一色素估算的最佳波段，不难发现红边区域一阶导数光谱是互花米草色素含量估算的最佳波段

。

图4

叶片光合色素含量实测值与预测值拟合散点图Fig．4Scatter plots describing the measured and predicted LPPC

3

讨论3．1冠层水平互花米草叶片光合色素含量估算的最佳波段

冠层水平，

植物的光谱反射率由绿色植物和非绿色植物的光谱组成，它受植物冠层结构（叶片的方向性和阴影等）和土壤背景等外部因素的影响［35］。为了消除外部因素对植物光谱产生的噪音对估测精度的影响，

38114期艾金泉等：冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型

4811生态学报35卷

研究者通常使用导数光谱技术、相关性分析、主成分分析、植被指数构建、红边位移效应分析等方法确定光合

色素估算的最佳波段。杨杰等［36］通过任意两波段组合的比值光谱指数和归一化光谱指数与水稻叶绿素含量

相关性分析表明，709、724、745nm和749nm波段组合的指数能较好的预测Chla、Chlb的含量。Zarco-Tejadaa

等［37］通过叶片辐射传输模型PＲOSPECT-5模拟发现由515nm和570nm两波段组合的比值植被指数Ｒ515/Ｒ

用于人工种植葡萄的Cars的估算精度最高。Kochubey等［38］对冬小麦、玉米、甜菜和野生葡萄四种农作物570

叶片叶绿素含量变化引起一阶导数光谱反射率值变化的研究发现，当农作物叶片的叶绿素含量从小到大变化

时，红边区域（690—730nm）一阶导数光谱反射率值的变化也出现从小到大的相应变化，725、715nm及702

nm是表现最敏感的波段。

本文充分借鉴以上研究成果，利用导数光谱技术和相关系数法（图2中的峰值点或拐点）得到了互花米

草光合色素敏感的波段，并对其构建光谱指数进行验证（表4），表明723、703、525、752、668nm的一阶导数

光谱是估测互花米草叶片光合色素含量的敏感波段。由此可知，冠层水平，红边导数光谱（680—760nm）是

估测互花米草叶片光合色素含量的最佳波段，这一研究结论与冯伟等［31］，Main等［4］研究结果类似，只是所选取的具体波段不一样。这可能有以下几方面的原因［4，7］：（1）研究表明，红边是主要由植物的色素含量、水分含量等固有的生化参量决定的，并且相对于其他波段（蓝光或红外波段）不易受土壤等背景影响的波段，而运用红边导数光谱处理进一步的消除了土壤等背景光谱的影响，故用红边一阶导数波段组合的光谱指数反演叶绿素含量效果比较好；（2）估算的具体植物类别不同，他们的冠层结构可能有差异，且湿地植物光谱变异性强，所以导致具体使用的光谱波段不一样。

3．2互花米草光合色素含量高光谱遥感估测的可行性

利用高光谱遥感技术精确估测湿地植物互花米草的光合色素，对于监测其生长状况和生态评估管理具有

重要的意义。虽然本文对互花米草光合色素含量估算取得了较好的效果，但本研究取得的研究成果是对实验

条件进行了严格的控制下取得的。首先采样点位于闽江河口湿地生态定位站观鸟台以西200m附近的小区

域，湿地土壤单一，地形平坦；其次，互花米草的叶片冠层光谱采集是基于模拟陆地植物光谱条件下采集的，包

括采集的时间是互花米草不受潮汐影响的时间段，地面没有明显积水以及尽量选取接近平直型的冠层进行光

谱采集。但真实的湿地环境复杂，不但湿地植物生物多样性高，而且受到潮汐的影响，而本研究只研究了秋季

互花米草色素含量的高光谱估算模型，并没有对不同季节、不同湿地环境条件下互花米草光合色素估测模型

进行研究，这也是今后需进一步深入研究的方向。此外，在冠层水平下证明用高光谱能够准确测定互花米草

的光合色素含量，在卫星遥感也不一定行，还需研究相应的反演模型进一步能使本研究的结论推广至航空或

航天高光谱遥感。

4结论

本研究通过导数光谱技术、相关性分析法和构建新的植被指数等高光谱处理方法，利用地面实测数据，比

较22种已报道的植被指数和14种新构建的植被指数在估算冠层水平互花米草叶片光合色素含量中的表现，

得出的主要结论如下：

（1）不同色素含量与原始光谱反射率和一阶导数光谱反射率在各波长范围相关性变化趋势具有一致性，但其相关系数大小有所差异。

（2）已报道指数和新构建的指数均表明，一阶导数光谱反射率组合的植被指数用于估算互花米草光合色素的含量优于原始光谱反射率。

（3）通过导数光谱技术、相关性分析和最优估算模型精度比较得出，互花米草光合色素估算的敏感波段依次是723、703、525、752、668nm的一阶导数光谱。

（4）对于单一色素的估算，Chla、Chlb、Cars的最佳估算模型分别是由FDNDVI［723，703］、FDＲVI［723，

525］、FDNDVI［723，703］建立的线性估算模型。

http：//www．ecologica．cn

http ：//www．ecologica．cn （5）对于不同色素使用同一植被指数估算而言，指数FDＲVI ［723，703］建立的对数估算模型效果最好。参考文献（Ｒeferences ）：

［1］Xue L H ，Yang L Z．Deriving leaf chlorophyll content of green-leafy vegetables from hyperspectral reflectance．ISPＲS Journal of Photogrammetry

and Ｒemote Sensing ，2009，64（1）：97-106．

［2］王福民，黄敬峰，刘占宇，王秀珍．水稻色素含量估算的最优比值色素指数研究．浙江大学学报：农业与生命科学版，2009，35（3）：

321-328．

［3］Wu C Y ，Niu Z ，Tang Q ，Huang W J．Estimating chlorophyll content from hyperspectral vegetation indices ：Modeling and validation．Agricultural

and Forest Meteorology ，

2008，148（8/9）：1230-1241．［4］Main Ｒ，Cho M Z ，Mathieu Ｒ，O'Kennedy M M ，Ｒamoelo A ，Koch S．An investigation into robust spectral indices for leaf chlorophyll estimation．

ISPＲS Journal of Photogrammetry and Ｒemote Sensing ，2011，66（6）：751-761．

［5］Sims D A ，Gamon J A．Estimation of vegetation water content and photosynthetic tissue area from spectral reflectance ：a comparison of indices

based on liquid water and chlorophyll absorption features．Ｒemote Sensing of Environment ，2003，84（4）：526-537．

［6］Leonard L A ，Luther M E．Flow hydrodynamics in tidal marsh canopies．Limnology and Oceanography ，1995，40（8）：1474-1484．

［7］Adam E ，Mutanga O ，Ｒugege D．Multispectral and hyperspectral remote sensing for identification and mapping of wetland vegetation ：A review．

Wetlands Ecology and Management ，

2010，18（3）：281-296．［8］卢霞，刘付程，田慧娟．大米草室内叶片光谱特征参数与叶绿素浓度关系研究．测绘科学，2010，35（6）：99-102．

［9］李凤秀，张柏，刘殿伟，王宗明，宋开山，靳华安，刘焕军．湿地小叶章叶绿素含量的高光谱遥感估算模型．生态学杂志，2008，27（7）：

1077-1083．

［10］Ge S K ，Carrutherrs ＲI ，Spencer D F ，Yu Q．Canopy assessment of biochemical features by ground-based hyperspectral data for an invasive

species ，giant reed （Arundo donax ）．Environmental Monitoring and Assessment ，

2008，147（1/3）：271-278．［11］曾从盛，张林海，仝川．闽江河口湿地芦苇和互花米草氮、磷养分季节动态．湿地科学，

2009，7（1）：16-25．［12］高元竞．闽江河口湿地生态服务功能价值评价［D ］．福州：福建农林大学，2009．

［13］林辉，臧卓，刘秀英，孙华．森林树种高光谱遥感研究．北京：中国林业出版社，2011：23-23．

［14］Gitelson A A ，Kaufman Y J ，Merzlyak M N．Use of a green channel in remote sensing of global vegetation from EOS-MODIS．Ｒemote Sensing of

Environment ，

1996，58（3）：289-298．［15］Haboudane D ，Millera J Ｒ，Tremblay N ，Zarcote-Jada P J ，Dextraze L．Integrated narrow-band vegetation indices for prediction of crop chlorophyll

content for application to precision agriculture．Ｒemote Sensing of Environment ，

2002，81（2/3）：416-426．［16］王强，易秋香，包安明，罗毅，赵金．基于高光谱反射率的棉花冠层叶绿素密度估算．农业工程学报，2012，28（15）：125-132．

［17］Sims D A ，Gamon J A．Ｒelationships between leaf pigment content and spectral reflectance across a wide range of species ，leaf structures and

developmental stages．Ｒemote Sensing of Environment ，2002，81（2/3）：337-354．

［18］Tucker C J．Ｒed and photographic infrared linear combinations for monitoring vegetation．Ｒemote Sensing of Environment ，1979，8（2）：127-150．［19］Gamon J A ，Serrano L ，Surfus J S．The photochemical reflectance index ：an optical indicator of photosynthetic radiation use efficiency across

species ，functional types ，and nutrient levels．Oecologia ，1997，112（4）：492-501．

［20］Gitelson A A ，Merzlyak M N．Spectral reflectance changes associated with autumn senescence of Aesculus hippocastanum L．and Acer platnoides

L．leaves．spectral features and relation to chlorophyll estimation．Journal of Plant Physiology ，

1994，143（3）：286-292．［21］Yang X H ，Huang J F ，Wang J W ，Wang X Z ，Liu Z Y．Estimation of vegetation biophysical parameters by remote sensing using radial basis

function neural network．Journal of Zhejiang University Science A ，

2007，8（6）：883-895．［22］Gitelson A A ，Merzlyak M N．Signature analysis of leaf reflectance spectra ：algorithm development for remote sensing of chlorophyll．Journal of

Plant Physiology ，1996，148（3-4）：494-500．

［23］Pe uelas J ，Gamon J A ，Fredeen A L ，Merino J ，Field C B．Ｒeflectance indices associated with physiological changes in nitrogen-and water-

limited sunflower leaves．Ｒemote Sensing of Environment ，

1994，48（2）：135-146．［24］Kim M S ，Daughtry C S T ，Chappelle E W ，McMurtrey III J E ，Walthall C L．The Use of High Spectral Ｒesolution Bands for Estimating Absorbed

Photosynthetically Active Ｒadiation （A par ）．France ：Proceedings of the 6th Symposium on Physical Measurements and Signatures in Ｒemote

Sensing Vai D'Isere ，

1994：299-306．［25］Huete A Ｒ，Liu H Q ，Batchily K ，Vanleeuwen W．A comparison of vegetation indices over a global set of TM images for EOS-MODIS．Ｒemote

Sensing of Environment ，1997，59（3）：440-451．

［26］Daughtry C S T ，Walthall C L ，Kim M S ，de Colstoun E B ，McMurtrey J E．Estimating corn leaf chlorophyll concentration from leaf and canopy

5

8114期艾金泉等：冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型

6811生态学报35卷

reflectance．Ｒemote Sensing of Environment，2000，74（2）：229-239．

［27］Boochs F，Kupfer G，Dockter K，Kühbauch W．Shape of the red edge as vitality indicator for plants．International Journal ofＲemote Sensing，1990，11（10）：1741-1753．

［28］Vogelmann J E，Ｒock B N，Moss D M．Ｒed edge spectral measurements from sugar maple leaves．International Journal ofＲemote Sensing，1993，14（8）：1563-1575．

［29］Le Maire G，Fran ois C，Dufrêne E．Towards universal broad leaf chlorophyll indices using PＲOSPECT simulated database and hyperspectral reflectance measurements．Ｒemote Sensing of Environment，2004，89（1）：1-28．

［30］蒋金豹，陈云浩，黄文江．用高光谱微分指数估测条锈病胁迫下小麦冠层叶绿素密度．光谱学与光谱分析，2010，30（8）：2243-2247．［31］冯伟，朱艳，田永超，马吉锋，庄森，曹卫星．基于高光谱遥感的小麦冠层叶片光合色素密度监测．生态学报，2008，28（10）：4902-4915．

［32］任红艳，史学正，庄大方，江东，徐新良，黄耀欢，刘磊，施润和．土壤全氮含量与碳氮比的高光谱反射估测影响因素研究．遥感技术与应用，2012，27（3）：372-379．

［33］梁亮，杨敏华，张连蓬，林卉，周兴东．基于SVＲ算法的小麦冠层叶绿素含量高光谱反演．农业工程学报，2012，28（20）：162-171．［34］黄敬峰，王福民，王秀珍．水稻高光谱遥感实验研究．杭州：浙江大学出版社，2010：136-161，284-300．

［35］Blackburn G A．Quantifying chlorophylls and carotenoids at leaf and canopy scales：an evaluation of some hyperspectral approaches．Ｒemote Sensing of Environment，1998，66（3）：273-285．

［36］杨杰，田永超，姚霞，曹卫星，张玉森，朱艳．水稻上部叶片叶绿素含量的高光谱估算模型．生态学报，2009，29（12）：6561-6571．［37］Zarco-Tejadaa P J，Guillén-Climenta M L，Hernández-ClementebＲ，Catalinac A，Gonzálezc MＲ，Martínc P．Estimating leaf carotenoid content in vineyards using high resolution hyperspectral imagery acquired from an unmanned aerial vehicle（UAV）．Agricultural and Forest Meteorology，2013，171-172：281-284．

［38］Kochubey S M，Kazantsev T A．Changes in the first derivatives of leaf reflectance spectra of various plants induced by variations of chlorophyll content．Journal of Plant Physiology，2007，164（12）：1648-1655．

http：//www．ecologica．cn

本文来源：https://www.bwwdw.com/article/ouqq.html