遥感科学与应用技术

南方红壤侵蚀区长汀县不同生态恢复年限下芒萁叶绿素含量的高光谱估算模型

  • 邓超 , 1, 2 ,
  • 陈志彪 , 1, 2, * ,
  • 陈海滨 1, 2 ,
  • 陈志强 1, 2
展开
  • 1. 福建师范大学 湿润亚热带生态地理过程教育部重点实验室,福州 350007
  • 2. 福建师范大学地理科学学院,福州 350007
*通讯作者:陈志彪(1962-),男,福建平潭人,教授,博士,主要从事水土保持与资源环境等方面研究。E-mail:

作者简介:邓超(1992-),女,四川资阳人,硕士生,研究方向为侵蚀过程与生态调控。E-mail:

收稿日期: 2018-10-12

  要求修回日期: 2019-02-25

  网络出版日期: 2019-06-15

基金资助

国家重点研发计划项目(2016YFC0502905)

福建省社会发展引导性(重点)项目(2016Y0024)

福建省教育厅中青年教师教育科研项目(JT180292)

Hyperspectral Estimation Models of Chlorophyll Content for Dicranopteris Dichotoma Leaves at Different Ecological Restoration Stages in the Eroded Red Soil Areas of Southern China

  • DENG Chao , 1, 2 ,
  • CHEN Zhibiao , 1, 2, * ,
  • CHEN Haibin 1, 2 ,
  • CHEN Zhiqiang 1, 2
Expand
  • 1. Key Laboratory for Humid Subtropical Eco-geographical Processes of the Ministry of Education, Fujian Normal University, Fuzhou 350007, China
  • 2. College of Geographical Science, Fujian Normal University, Fuzhou 350007, China
*Corresponding author: CHEN Zhibiao, E-mail:

Received date: 2018-10-12

  Request revised date: 2019-02-25

  Online published: 2019-06-15

Supported by

National Key Research and Development Program of China, No.2016YFC0502905

Social Development Pilot Project of Fujian, No.2016Y0024

Education Research Project for Young and middle-aged Teachers of Fujian Education Department, No.JT180292

Copyright

《地球信息科学学报》编辑部 所有

摘要

芒萁是南方红壤侵蚀区生态恢复重要的地带性草本植物,对生态系统修复具有重要作用,监测其叶绿素含量能有效诊断生长健康状况。本文以福建省长汀县朱溪流域6个不同生态恢复年限下的芒萁叶片高光谱反射数据以及实测叶绿素含量为数据源,借助高光谱遥感技术分析不同恢复年限芒萁叶片原始光谱特征,筛选出光谱敏感波段并构建光谱指数,基于相关性分析,建立芒萁叶绿素单变量以及多元逐步回归模型,并确定最佳估算模型。结果表明:高光谱指数建立的单变量估算模型中,改进红边归一化植被指数(mNDVI705)、叶面叶绿素指数(LCI)、红边指数(Vog)、比值光谱指数(RVI603/407)、NDVI[603,407]高光谱指数建立的二次模型精度高,建模决定系数R2均超过了0.8,其中以高光谱指数为自变量建立的多元回归模型拟合R2值(0.886)最高。综合建模精度和模型验证精度,LCI指数构建的单变量模型以及基于高光谱指数的多元回归模型是估算芒萁叶片叶绿素含量最佳模型。本研究建立的叶绿素高光谱估算模型对快速、无损地监测水保植物芒萁生长具有重要意义。

本文引用格式

邓超 , 陈志彪 , 陈海滨 , 陈志强 . 南方红壤侵蚀区长汀县不同生态恢复年限下芒萁叶绿素含量的高光谱估算模型[J]. 地球信息科学学报, 2019 , 21(6) : 948 -957 . DOI: 10.12082/dqxxkx.2019.180513

Abstract

Dicranopteris dichotoma is an important zonal herbaceous plant for ecological restoration in the eroded red soil areas of southern China, and is an effective control of soil erosion. Using remote sensing techniques to monitor the chlorophyll content can help diagnose the vegetation growth and healthy condition of dicranopteris dichotoma. Based on hyperspectral reflectance data and the corresponding chlorophyll contents of dicranopteris dichotoma leaf from six different ecological restoration stages in Zhuxi small watershed of Changting County, Fujian Province, this study analyzed the hyperspectral curve properties of leaf, transformed the original spectral into the first derivative, and selected the sensitive wavebands to create ratio (RVI) and normalized (NDVI, FDNDVI) hyperspectral indices. Then correlation analysis was conducted for the chlorophyll contents and hyperspectral indices which were selected from reported indices and newly constructed indcies with sensitive wavelengths. Based on the correlation coefficients, we can chose the best indices to create estimation models. The linear, exponential, multiplicative, quadratic polynomial, logarithmic, and multivariate regression models were constructed for comparison. Furthermore, the optimal estimation model was determined by the accuracy of each estimation model. Results showed that the sensitive wavelengths of the original spectral for dicranopteris dichotoma leaf at different ecological restoration stages were 407 nm, 603 nm, and that the optimal wavebands of the first derivative were 463 nm, 554 nm, 674 nm, and 739 nm. The relationship between the chlorophyll content of dicranopteris dichotoma leaf and the hyperspectral indices of red edge position (λr), NDVI[603, 407], Modified Red Edge Normalized Difference Vegetation Index (mNDVI705), Vogelmann Index (Vog) were very significant, and the correlation coefficients were over 0.85. The estimation models of chlorophyll content established by hyperspectral indices of mNDVI705, Leaf Chlorophyll Index (LCI), Vog, RVI603/407, NDVI[603, 407] showed better test results, and the R2 were over 0.8. The model established by FDNDVI[739, 463] index had the highest verification accuracy, and the R2 reached 0.741. The multivariate regression model based on hyperspectral indices got highest test results with the highest R2. Therefore, the LCI index and the multivariate regression model based on hyperspectral indices have the strongest ability for predicting chlorophyll concentration, which provides scientific basis for dynamic monitoring of dicranopteris dichotoma in the eroded red soil regions of southern China. It is significant for monitoring soil and water conservation plants. Meanwhile, the objective of this research was to provide effective technical support for ecological restoration by building hyperspectral estimation models of chlorophyll content, with a rapid and non-destructive method for monitoring vegetation growth.

1 引言

叶绿素含量是植物的重要生长参数,是植物进行光合作用重要的物质[1]。叶绿素含量的变化可以反映植物叶片光合作用功能的强弱以及表征逆境胁迫下植物组织、器官的衰老状况,对植物的光合能力、发育阶段及营养状况有重要的指示作用[2],通过研究植物的叶绿素含量能有效的监测植物生长发育状况。传统的化学分析法得到植物的叶绿素含量虽然结果可靠,但是对植物具有破坏性且不能对同一植物进行长期的动态监测[3]。高光谱遥感技术具有波段多且窄的特点[4],并具有省时、减少破坏的优点,能定量分析植被光谱差异,识别植物类型,反演植物叶绿素、水分等生理生化参数[5],感知环境因子对植物的影响[6]。高光谱技术在生态环境监测,尤其是在植被叶绿素含量变化研究中发挥着越来越重要的作用[7,8]。在湿地生态系统[9]以及草原生态系统[10]植被监测的相关研究中取得了较好的效果。同时,通过建立苹果树[11]、白毛杨[12]、马尾松[13]等的叶绿素高光谱估算模型能有效监测单一植被物种的生长状况。
严重的水土流失造成土壤养分流失、土壤质量下降,加剧生态系统的退化[14]。芒萁作为水土保持的先锋植物,是严重退化生态系统最后退出,也是治理中最早侵入的地带性草本之一。在区域生态恢复过程中,芒萁能很快替代其他人工种植的草本,成为优势草种,其生长的好坏可反映对水土流失的控制及生态恢复的效果。叶绿素是诊断植被健康状况的重要指标,通过建立叶绿素含量的高光谱估算模型可快速监测芒萁的生长状况。选择土壤侵蚀区不同生态恢复年限下芒萁进行研究,对南方红壤侵蚀区生态恢复与重建具有重要意义[15,16]。近年来,对水稻、玉米等农作物和对杨树、梧桐等树种的生理生化参数进行深入的研究,而利用高光谱遥感技术对严重退化区生态恢复具有重要作用的水土保持植物芒萁叶片叶绿素含量研究还尚未见报道。
本文以福建长汀县土壤侵蚀区朱溪流域的地带性草本芒萁为研究对象,分析不同生态恢复年限下芒萁叶片光谱特征以及通过高光谱指数构建叶绿素含量的估算模型,为南方红壤侵蚀区芒萁的长势监测提供技术支撑。

2 研究区概况与数据源

2.1 研究区概况

福建省长汀县是中国南方红壤侵蚀区的典型代表地区,是水土流失最严重的县份之一,由于不合理的土地利用方式,该地区土壤和植被遭到严重毁坏,生物多样性减少,土地生产力下降,导致严重的水土流失和生态系统退化[17]。其地理位置为116°23'30E"-116°30'30"E,25°38'15N"-25°42'55"N。流域地势自东北向西南倾斜,海拔270~680 m,地貌以低山、丘陵为主,地形较为破碎,在河流沿岸及支流有河谷盆地分布。气候属于中亚热带季风性湿润气候,流域常年平均气温约为18.4 ℃,1月平均气温大约为8 ℃,7月平均气温约为27.2 ℃,年降雨量介于1400~2450 mm之间,全年降水量呈典型的双峰型。土壤以湿热的气候下风化而发育成的红壤为主,风化壳深厚,一般为10~20 cm,节理发育,结构松软,抗侵蚀性差,土层保水保肥能力差。

2.2 样品测定

在研究区内分别选择生态恢复年限为6a、12a、18a、36a的来油坑(6a)、龙颈(12a)、游坊(18a)、八十里河(36a)的4个生态恢复样地(图1),以及未治理对照地和风水林对照样地各一个,其中每一个样地选择上、中、下坡为采样点,并以东南西北4个方向采样,将采集好的叶片迅速装入自封袋并编号,放入保鲜箱中,野外采样完成后尽快带回实验室。
Fig. 1 Overview of the study area and laboratory experiment

图1 研究区样地概况和室内实验

2.3 光谱数据测定

测定芒萁叶片光谱采用光谱仪Field Spec3,该光谱仪波段范围为350~2500 nm,光谱测量在室内进行,能有效避免测量多个样品时因角度不同造成光谱信息误差。在测定前进行标准白板校正,每个样品重复测定10次光谱,剔除光谱反射率大于1的异常数据,对处理后的光谱数据进行分析,求均值并作为该样品最终光谱反射率。测定72个芒萁叶片样本共720条光谱数据,利用ViewSpePro软件分析光谱数据。考虑到光谱曲线末端噪音较大以及植物色素的影响集中于可见光波段[18],本文选取400~1350 nm的波段范围进行分析。

2.4 叶绿素含量测定

将测定完光谱数据的芒萁叶片从保鲜袋中取出,继续测定叶片的叶绿素含量。首先将叶片剪碎、混匀,每份称取0.2 g,倒入研钵中,分别加入少量的无水C2H5OH及适量的石英砂和CaCO3粉,研磨成匀浆,再用无水C2H5OH洗入离心管中,离心之后将上清液倒入25 ml容量瓶,再用无水C2H5OH定容至25 ml,摇匀之后,用紫外可见分光光度计测量665 nm和649 nm处的吸光度[19],并采用叶绿素浓度计算公式计算芒萁叶片浓度含量。
叶绿素a浓度(Chla)=13.95×D665-6.88×D649 (1)
叶绿素b浓度(Chlb)=24.96×D665-7.32×D649 (2)
叶绿素总浓度(Chl)=(Chla+Chlb)×(V/M)(3)
叶绿素含量=叶绿素总浓度×浸提液体积×稀释倍数/样品质量(4)
式中:D665D649分别表示在665 nm和649 nm处的吸光度值;V为定容体积/ml;M为样品重量/g。叶绿素浓度单位为mg/L,含量单位为mg/g。

3 研究方法

目前用于估算叶绿素含量的高光谱指数,主要包括比值植被指数和归一化植被指数。其中,红边参数能有效反映叶绿素含量的变化,尤其是定量描述红边特征的红边位置(λr)、红边幅值(Dr)和红边峰值面积(SDr)。另外,对光谱曲线进行导数处理可有效消除干扰信息。本文在原始光谱曲线特征、一阶微分光谱与叶绿素含量的相关性分析基础上,挑选出相关系数高的敏感波段,并通过波段组合构建比值和归一化光谱指数,波段组合在反映植物信息方面比单波段的数据更可靠。从各指数(表1)中挑选与叶绿素含量相关性较高的14种高光谱指数建立单变量回归模型。为充分利用波段信息来估算芒萁叶绿素含量,进一步通过多元回归分析[20]进行波段选择与重组。多元回归分析的一般形式为:
Y = b 0 + b 1 x 1 + b 2 x 2 + + b i x i (5)
式中:Y代表芒萁的叶绿素含量;bi代表第i个特征变量的回归系数;b0为回归常数;xi表示第i个高光谱指数。
Tab. 1 Hyperspectral indices for estimating the chlorophyll content of dicranopteris dichotoma leaves

表1 估算芒萁叶片叶绿素含量的高光谱指数

高光谱指数类型 高光谱指数 定义或计算公式 参考文献
常用光谱指数 红边位置(λr) 红边位置 [15],[21]-[22]
红边振幅(Dr) 红边内一阶微分的最大值
红边面积(SDr) 红边内一阶微分的总和
蓝边位置(λb) 蓝边位置
蓝边振幅(Db) 蓝边内一阶微分的最大值
蓝边面积(SDb) 蓝边内一阶微分的总和
红边内一阶微分总和与蓝边内一阶微分总和比值(SDr/SDb) 红边内一阶微分总和与蓝边内一阶微分总和比值
红边内一阶微分归一化总和与蓝边内一阶微分总和的归一化值 (SDr-SDb)/(SDr+SDb)
归一化植被指数(NDVI) (R800-R670)/(R800+R670) [23]
红边归一化植被指数(NDVI705 (R750-R705)/(R750+R705) [24]
改进红边归一化植被指数(mNDVI705 (R750-R705)/(R750+R705-2R445) [25]
改进红边比值植被指数
(mSR705
(R750-R445)/(R750+R445) [26]
结构不敏感色素指数(SIPI) (R800-R445)/(R800+R680) [27]
植被衰减指数(PSRI) (R680-R500)/R750 [28]
叶面叶绿素指数(LCI) (R850-R710)/(R850+R680) [11]
归一化比值叶绿素指数(PSND) (R810-R674)/(R810+R674) [29]
红边指数(Vog) R740/R720 [30]
敏感波段构建的光谱指数 RVI750/700 R750/R700 [31]
RVI603/407 R603/R407 本文
NDVI[603, 407] (R603-R407)/(R603+R407)
FDRVI[554, 463] DR554/DR463
FDRVI[739, 674] DR739/DR674
FDRVI[674, 554] DR674/DR554
FDRVI[739, 554] DR739/DR554
FDNDVI[739, 463] (DR739-DR463)/(DR739+DR463)
FDNDVI[739, 554] (DR739-DR554)/(DR739+DR554)
FDNDVI[554, 463] (DR554-DR463)/(DR554+DR463)
FDNDVI[674, 554] (DR674-DR554)/(DR674+DR554)
FDNDVI[674, 463] (DR674-DR463)/(DR674+DR463)

注:Rx表示波长在x处的原始光谱反射率;DRx同表示波长在x处的一阶微分光谱反射率;FDRVI[x1, x2]表示在波长x1, x2处一阶微分的比值植被指数;FDNDVI[x1, x2]为在波长x1, x2处一阶微分的归一化植被指数;(DRα-DRβ)/(DRα+DRβ)为在波长α, β处的归一化植被指数。

将挑选出的70组芒萁叶片样本随机分成两组,其中54组样本数据用于建模,16组样本用于验证。将高光谱指数作为自变量,芒萁叶片叶绿素含量作为因变量,分别建立线性、指数、乘幂、二次多项式、对数模型,决定系数(R²)最大的模型选定为每一指数的最优模型。在模型建立的同时通过验证数据集对模型进行检验,结合验证模型的决定系数、均方根误差(RMSE)(式(6))和相对误差(RE)(式(7))评价实测值与估算值的拟合效果,并绘制二者之间的关系图,筛选出芒萁叶片叶绿素含量最佳的预测模型,估算模型具体流程如图2所示。
Fig. 2 Flowchart of determining the best dicranopteris dichotoma chlorophyll estimation model

图2 芒萁叶绿素估算模型流程

RMSE = i = 1 n ( y i - y i ' ) 2 n (6)
RE = 1 n i = 1 n y i ' - y i y i × 100 % (7)
式中:yi表示实测值;yi'代表预测值;n为样本数量。

4 结果及分析

4.1 不同生态恢复年限下芒萁叶片原始光谱特征

不同生态恢复年限下的芒萁光谱曲线均呈现出典型植物光谱特征,在波长500 nm处是蓝紫光对叶绿素的吸收谷,550 nm附近存在一个明显的反射峰;在波长680 nm处是红光吸收谷,由于红光对叶绿素的强烈吸收,不同生态恢复年限下的芒萁叶片反射率在两处吸收谷差异较小,在可见光区域反射率受各种色素影响较大;在近红外波段,光谱曲线受到水分和叶片结构的影响较大,尤其在700~750 nm之间反射率急剧增加,形成一个陡坡称为“红边”,风水林对照地的红边斜率最大,龙颈(12a)的红边斜率最小。在800~1350 nm之间出现一个高反射率平台,反射率值仅出现小范围波动。不同生态恢复年限的芒萁光谱反射率值存在差异。从图3可看出,不同生态恢复年限下的芒萁叶片绿峰值不同,风水林对照地的芒萁叶片反射率最大,此处芒萁叶片叶绿素含量实测值也最大。在700~1350 nm间,风水林对照地的芒萁叶片反射率远远高于其他样地,光谱反射率最大达到0.6,生态恢复12年的芒萁光谱反射率最小,差异性可能是受到细胞结构的影响[25]
Fig. 3 The characteristics of original spectral curves for dicranopteris dichotoma leaf at different ecological restoration stages in Zhuxi small watershed of Changting County

图3 长汀县朱溪流域不同生态恢复年限下的芒萁叶片原始光谱反射率曲线变化特征
注:WZL指未进行生态修复的未治理对照样地;FSL指生态恢复较好的风水林对照样地;6a、12a、18a、36a代表生态恢复年限分别为6年、12年、18年和36年。

4.2 芒萁叶片叶绿素含量与原始光谱反射率及一阶导数的相关分析

将不同生态恢复年限下的芒萁叶片叶绿素含量与原始光谱数据进行相关性分析。从图4可看出:原始光谱反射率与叶绿素含量在407 nm处相关系数最大(+0.520),随着波长增加相关系数值逐渐减小,在波长为603 nm处达到最大负相关(-0.430)。699 nm之后随着波长增加,相关系数值快速增大,到760 nm处相关系数值趋于稳定。可见,407 nm和603 nm波段可作为估算芒萁叶片叶绿素含量的敏感波长。芒萁光谱一阶微分与叶绿素含量的相关性分析(图5)发现,相关系数波动性较大,出现多个极点和拐点,波长674 nm处相关系数值最大可达0.803;在绿光波段范围内,波长554 nm处达到最大负相关(-0.801)。在350~700 nm可见光波段叶绿素含量与一阶微分存在较大的相关性,在463 nm、739 nm附近均出现较高的相关性,相关系数值分别为-0.737、0.725,当波长大于1000 nm之后,相关性较弱,较少能通过0.01(r=0.303,n=54)极显著检验水平。因此,可以选取一阶微分与芒萁叶片叶绿素含量相关性较大的4个波长作为估算叶绿素含量的敏感波长。
Fig. 4 Correlation coefficients between chlorophyll content and the original spectral data of dicranopteris dichotoma in Zhuxi small watershed of Changting County

图4 长汀县朱溪流域芒萁叶绿素含量与原始光谱反射率相关系数
注:图中虚线表示显著性达0.01检验水平时相关系数的临界值。

Fig. 5 Correlation coefficients between chlorophyll content and the first derivative spectral data of dicranopteris dichotoma in Zhuxi small watershed of Changting County

图5 长汀县朱溪流域芒萁叶绿素含量与一阶微分光谱的相关系数
注:图中虚线表示显著性达0.01检验水平时相关系数的临界值。

4.3 高光谱指数与叶绿素含量的相关性

将选取的30个高光谱指数与叶片叶绿素含量进行相关性分析(表2)。芒萁叶片叶绿素含量与红边位置(λr)呈极显著相关,相关系数达0.862,但与红边幅值、红边面积相关系数值都未能通过0.01显著性检验,说明芒萁叶片叶绿素含量与红边位置关系较为密切,可将红边位置作为参量来构建估算模型。蓝边(490~530 nm)内最大一阶微分值(Db),蓝边位置(λb)与叶绿素含量相关性不显著;红边面积与蓝边面积比值(SDr/SDb)与叶绿素含量呈极显著正相关,相关系数值达0.758。
Tab. 2 Correlation coefficients between the hyperspectral indices and chlorophyll content of dicranopteris dichotoma in Zhuxi small watershed of Changting County

表2 长汀县朱溪流域叶绿素含量与高光谱指数的相关性分析

高光谱指数 相关系数 高光谱指数 相关系数
λr 0.862 ** SRI 0.834 **
Dr -0.023 PSND 0.377 **
SDr 0.424 ** Vog 0.860 **
λb -0.083 RVI750/700 0.808 **
Db -0.168 RVI603/407 -0.754 **
SDb -0.168 NDVI[603, 407] -0.874 **
SDr/SDb 0.758 ** FDRVI[554, 463] -0.292
(SDr-SDb)/(SDr+SDb) 0.114 FDRVI[739, 674] 0.588 **
NDVI 0.443 ** FDRVI[674, 554] -0.350 **
NDVI705 0.824 ** FDRVI[739, 554] 0.669 **
mNDVI705 0.857 ** FDNDVI[739, 463] 0.755 **
mSR705 0.162 FDNDVI[739, 554] 0.724 **
SIPI -0.718 ** FDNDVI[554, 463] 0.393 **
PSRI -0.330 ** FDNDVI[674, 554] 0.278
LCI 0.848 ** FDNDVI[554, 688] -0.759 **

注:** 表示在0.01水平上显著。

4.4 基于高光谱指数的叶绿素含量估算模型建立与验证

4.4.1 单变量线性和非线性估算模型
在构建的植被指数中,NDVI[603, 407]、FDNDVI[739, 463]、FDNDVI[739, 554]与叶绿素含量的相关系数都较高,其中NDVI[603, 407]是相关系数最大的指数(表2)。挑选出14个相关性系数大于0.700的高光谱指数作为自变量,芒萁叶片叶绿素含量作为因变量分别建立线性、二次、对数、幂函数以及指数5种回归模型,再根据R²最大为优选原则确定出每一种指数的最佳拟合模型(表3),从表中可以看出以λr和SDr/SDb建立的最佳估算模型都为线性模型;SIPI的最佳估算模型为乘幂模型;SRI和RVI750/700最佳估算模型为指数模型;基于原始光谱数据与一阶微分光谱选取的敏感波长而构建的植被指数的最佳估算模型均为二次多项式模型。大部分模型的R2值都大于0.700,其中mNDVI705、LCI、Vog、RVI603/407、NDVI[603, 407]指数构建的最佳模型都为二次多项式,决定系数R2均超过了0.800。为了验证模型的适用性,以预测决定系数(R2)、均方根误差(MERS)、相对误差(RE)为评价指标,分析验证数据集检验模型的预测能力发现:预测R2值达0.870以上的高光谱指数共有6个,分别为FDNDVI[739, 463]、NDVI705、SRI、RVI750/700、Vog、LCI;MERS最小的4个指数分别是FDNDVI[739, 463]、NDVI705、Vog、LCI;RE最小的4个指数分别是FDNDVI[739, 463]、mNDVI705、NDVI705、LCI。综合3个指标确定模型验证精度,精度最高的指数为FDNDVI[739, 463],其次为LCI和NDVI705。同时,建模精度较高的指数分别为LCI和mNDVI,而FDNDVI[739,463]和NDVI705指数建模的R²值都处于中等水平,分别为0.741、0.750。红边位置的建模精度和验证精度都处于中等水平,通过敏感波段构建的植被指数中FDNDVI[739, 463]为芒萁叶绿素估算模型的最佳指数,已报道的高光谱指数中LCI、mNDVI都能较好的估算芒萁叶片叶绿素含量。LCI指数构建的二次模型的拟合R2、预测R2值都较高,是最佳单变量估算模型。
Tab. 3 Chlorophyll content estimation models of dicranopteris dichotoma leaf and evaluation indicators in Zhuxi small watershed of Changting County

表3 长汀县朱溪流域芒萁叶片叶绿素含量估算模型及评价指标

高光谱指数 模型建立 模型验证
回归方程 拟合R2 预测R2 MERS RE%
λr y=0.096x-66.543 0.779 0.859 0.363 22.87
SDr/SDb y=0.018x-0.962 0.592 0.776 0.332 19.92
NDVI705 y=1.794-9.103x+19.72x2 0.750 0.886 0.256 16.20
mNDVI705 y=2.831-13.123x+21.146x2 0.849 0.860 0.281 16.14
SIPI y=2.71x(-25.09) 0.577 0.683 0.420 21.59
LCI y=5.212-21.726x+26.671x2 0.865 0.872 0.273 16.38
SRI y=0.2e0.772x 0.743 0.882 0.287 17.12
Vog y=9.969-18.516x+9.133x2 0.820 0.878 0.261 17.20
RVI750/700 y=0.275e0.469x 0.693 0.877 0.306 17.40
RVI603/407 y=4.693-1.506x+0.139x2 0.833 0.782 0.326 18.38
NDVI[603, 407] y=4.583-9.133x+5.149x2 0.835 0.740 0.359 18.37
FDNDVI[739, 463] y=13.333-35.175x+24.404x2 0.741 0.908 0.235 15.47
FDNDVI[688, 554] y= -9.885+18.412x-6.812x2 0.728 0.689 0.410 29.43
FDNDVI[739, 554] y=1.370-3.818x+4.408x2 0.777 0.729 0.386 28.72
4.4.2 多元逐步回归模型
虽然以敏感波长为自变量、芒萁叶绿素含量为因变量建立多元回归方程拟合R2值达0.700以上,但模型的预测R2仅为0.634,预测精度较低(表4)。而基于高光谱指数建立的回归方程拟合R2值(0.886)高于敏感波段建立的多元回归模型和高光谱指数建立的单变量回归模型,预测R2值达0.800以上,精度较高(图6)。
Tab. 4 Multivariate regression models for the chlorophyll content estimation of dicranopteris dichotoma in Zhuxi small watershed of Changting County

表4 长汀县朱溪流域芒萁叶片叶绿素含量多元回归模型

光谱形式 回归方程 拟合R2 预测R2 RMSE
敏感波长 y=1.443-2359.116×R463'-1793.879×R554'-237.549×R674'+227.782×R739' 0.757 0.634 0.598
高光谱指数 y=-25.830+0.03x1+0.579x2-1.662x3-6.195x4+6.593x5 0.886 0.838 0.349

注:表中R'x代表在波长x处的一阶微分值;x1,x2,x3,x4,x5分别代表高光谱指数 λr,FDNDVI[688, 554],NDVI[603, 407],mNDVI705,Vog。

Fig. 6 Scatter plots for the relationship between predicted and measured values of the chlorophyll content in Zhuxi small watershed of Changting County

图6 长汀县朱溪流域芒萁叶片叶绿素含量实测值与预测值拟合散点图

最后综合建模精度和验证模型精度,就单变量指数来看,LCI指数建立的二次模型是不同生态恢复年限下芒萁叶片叶绿素含量的最佳估算模型;多元回归模型中,以λr,FDNDVI[688, 554],NDVI[603, 407],mNDVI705,Vog指数为自变量建立的模型拟合R2值最高,预测精度较高,是最佳的估算模型。

5 讨论

5.1 芒萁叶绿素含量的敏感波长

通过波段组合构建的光谱指数能更准确地反映植被信息。研究者利用相关系数法确定原始光谱曲线和一阶微分光谱中的敏感波段,构建了互花米草、水稻、杉木等[32]植被的光谱指数,发现红边是估测水稻的敏感波段,而互花米草的敏感波段主要集中在红边导数(680~760 nm),杉木敏感波段部分出现在红边位置附近,苹果树敏感波段不仅出现在红边附近还出现在500~600 nm范围内。因此,在借鉴前人研究方法的基础上本文选取敏感波长时分别基于原始光谱和一阶微分来确定,敏感波段主要为407、603、463、554、674和739 nm,芒萁叶片叶绿素含量的敏感波段与其他植物的敏感波段不同,这可能由于研究的植被种类不同而有差异,但红边都是重要的敏感波段选取处。

5.2 芒萁叶绿素含量的最佳光谱指数

一阶微分和原始光谱选取的波段构建的植被指数与叶绿素含量都表现出较好的相关性,这与前人的研究结果[9,11,25]类似。文中通过一阶微分构建的FDNDVI[739, 554]虽然预测效果较LCI指数略低,但比NDVI705、SIPI等常见光谱指数好。不同的植被指数具有不同的特点,本文系统分析了在估算叶绿素含量方面效果较好的高光谱指数,发现芒萁叶片最佳的高光谱指数为LCI、mNDVI,同时这些指数在小麦[7]、水稻[26]、苹果树[11]等其他植被叶绿素含量的估算中也有较好的表现。由于易受到营养胁迫的影响,在红壤侵蚀区,土壤养分流失会影响到芒萁的生长,而芒萁本身对养分需求较低,因此,为避免单一的植被指数不能准确的估算叶绿素含量的变化,针对红壤侵蚀区不同生态恢复条件下的芒萁叶片叶绿素,本文选择了单变量模型的同时还建立多元回归模型,发现基于高光谱植被的多元回归模型精度大于单变量模型,这与王家强[8]在胡杨叶绿素含量的估测模型研究中多元回归模型精度均大于一元线性模型的结果一致。

6 结论

本文基于芒萁高光谱数据和叶绿素实测数据,分析了南方红壤侵蚀区不同生态恢复阶段芒萁叶片叶绿素含量的原始光谱特征以及构建了单变量估算模型和多元逐步回归模型,并确定了芒萁叶绿素的最佳估算模型,结论如下:
(1) 不同生态恢复年限下的芒萁叶片原始光谱反射率在可见光范围内主要受到植被色素的影响,反射率存在一定的差异。生态系统恢复越好,光谱反射率的红边斜率越大。
(2) 叶片原始光谱与叶绿素含量在可见光区域主要表现出负相关,一阶微分光谱与叶绿素含量在可见光范围内相关性显著,并出现多个极值。原始光谱曲线敏感波长为407 nm和603 nm,一阶微分曲线的敏感波长是463 nm、554 nm、674 nm和739 nm。
(3) 在建立的单变量回归模型中,LCI构建的一元二次模型建模精度和验证模型精度均较高,都能较好的预测芒萁叶片叶绿素含量;多元回归模型中,以高光谱指数λr,FDNDVI[688, 554],NDVI[603, 407],mNDVI705,Vog为自变量的多元回归模型是最佳的预测模型。
本文研究发现LCI指数在芒萁叶绿素含量的估算中具有优势,与其他植物最佳估算指数结果不同。同时,多元回归模型精度优于单变量模型,这与前人研究结果类似。芒萁叶绿素估算模型是在某一红壤侵蚀区实验数据的基础上构建和验证,今后还可通过不同土壤侵蚀区,不同季节下的芒萁数据来对模型进行检验和完善,以便进一步提高模型的稳定性和区域拓展性。通过卫星图像的反演,可将实验室的实测数据与星载或机载高光谱遥感影像,实现对水保植物芒萁的大面积监测。

The authors have declared that no competing interests exist.

[1]
Reich P B.Effects of low concentrations of o(3) on net photosynthesis, dark respiration, and chlorophyll contents in aging hybrid poplar leaves[J]. Plant Physiology, 1983,73(2):291-296.

[2]
王会霞,石辉,李秧秧.城市大气环境下绿化植物叶片比叶重和光合色素含量[J].中国环境科学,2011,31(7):1134-1142研究了大叶女贞(Ligustrum lucidum Ait.)和小叶女贞(Ligustrum quihoui Carr.)在不同大气环境条件下叶片中Pb、Cd、叶绿素a (Chl a)、叶绿素b(Chl b)、类胡萝卜素(Car)的含量及比叶重(LMA),分析了因素之间的相互关系.不同环境条件下叶片中Pb、Cd、Chl a、Chl b、Chl (a+b)、Car的含量及LMA具有显著差异.随污染程度的加剧,叶片中Pb、Cd的含量、LMA和Chl a/Chl b升高,而Chl a、Chl b和Chl (a+b)含量则呈现相反的变化趋势.大叶女贞叶片中Car含量变化不明显,而小叶女贞叶片中Car含量在污染严重的采样点明显较其他采样点高.2种植物叶片中Pb、Cd的含量对叶片LMA和光合色素含量均具有明显影响,Chl a、Chl b及Chl (a+b)含量与叶片中Pb、Cd的含量呈显著负相关(P<0.05).叶片中Cd含量的增大导致Car含量增大,而Pb的作用与物种有关.叶片中Pb含量的增大导致2物种LMA增大,而Cd含量导致大叶女贞LMA增大而小叶女贞LMA变化不明显.污染胁迫下植物叶片会增大LMA、Chl a/Chl b,减小Chl a、Chl b和Chl (a+b)含量,据此可对大气污染进行生物监测.

[ Wang H X, Shi H, Li Y Y.Leaf mass per area and photosynthetic pigments of greening plant species under different urban atmospheric environment[J]. China Environmental Science, 2011,31(7):1134-1142. ]

[3]
Hansen P M, Schjoerring J K.Reflectance measurement of canopy biomass and nitrogen status in wheat crops using normalized difference vegetation indices and partial least squares regression[J]. Remote Sensing of Environment, 2003,86(4):542-553.

[4]
刘伟东,项月琴,郑兰芬,等.高光谱数据与水稻叶面积指数及叶绿素密度的相关分析[J].遥感学报,2000,4(4):279-283分析了北京大屯科技站水稻叶面积指数(LAI)、叶绿素密度(CH.D)与高光谱分辨率遥感数据在整个生育期内的变化过程。利用微分技术处理水稻群体反射光谱以减少土壤等低频背景光谱噪音的影响。通过单相关分析和逐步回归方法研究水稻LAI、CH.D分别与光谱反射率、反射率的一阶微分光谱的相关关系, 并建立预测回归方程。结果表明, 微分技术能够改善光谱数据与LAI、CH.D的相关性, CH.D与光谱数据的相关明显优于同LAI的。

DOI

[ Liu W D, Xang Y Q, Zheng L F, et al.Relationships between rice LAI, CH.D and hyperspectra data[J]. Journal of Remote Sensing, 2000,4(4):279-283. ]

[5]
王润生. 高光谱遥感的物质组分和物质成分反演的应用分析[J].地球信息科学学报,2009,11(3):261-267高光谱遥感的最大特点是可以提取和重建像元光谱,从而依据光谱特征直接识别地物类型、地物组成,乃至地物的成分,反演地物的物理、化学参量。本文阐述了高光谱的光谱识别方法所能探测的岩石、土壤、植被和人工建筑物的物质组分或成分;另外,与多光谱相比较的角度,讨论了高光谱图像地物特征识别和物质成分反演等技术方法的应用。

[ Wang R S.Spectral identification and inversion of composition and component of objects with hyperspectral remote sensing[J]. Journal of Geo-information Science, 2009,11(3):261-267. ]

[6]
彭羽,米凯,覃亚,等.半干旱草原不同生态恢复阶段优势植物的高光谱特征分析[J].光谱学与光谱分析,2014,34(11):3090-3094为了将高光谱遥感技术应用到生态恢复监测, 采用ASD Fields HH便携式野外光谱仪测量半干旱草原区内蒙古和林格尔县不同生态恢复阶段优势植物物种的高光谱数据, 通过ViewSpecPro软件对原始光谱数据进行降维和去噪等处理, 进行原始反射率、 一阶微分分析, 采用排序软件Canoco 4.5进行去趋势典范对应分析。 结果表明, 生态恢复早期的优势植物物种为狗尾草、 小叶锦鸡儿, 生态恢复5年的优势植物物种为猪毛菜、 小叶锦鸡儿, 生态恢复后期的优势植物物种为樟子松、 猪毛菜, 与实测结果一致。 DCCA分析结果表明, 生态恢复早期植被的冠层光谱值, 特征波段主要为短波光, 且各物种光谱值差异较大; 生态恢复5年主要受1 000~1 050 nm波段反射率影响, 生态恢复后期主要受1 040~1 075 nm波段的近红外反射率影响。 采用DCCA方法能够直观反映不同生态恢复阶段优势物种的光谱特征。

[ Peng Y, Mi K, Qin Y, et al.Spectral reflectance characteristics of dominant plant species at different eco-restoring stages in the semi-arid grassland[J]. Spectroscopy and Spectral Analysis, 2014,34(11):3090-3094. ]

[7]
梁亮,杨敏华,张连蓬,等.基于SVR算法的小麦冠层叶绿素含量高光谱反演[J].农业工程学报,2012,28(20):162-171,294为给小麦的长势监测与农艺决策提供科学依据,利用高光谱技术实现了小麦冠层叶绿素含量的估测。通过分析18种高光谱指数对叶绿素的估测能力,筛选出可敏感表征叶绿素含量的指数REP,利用地面光谱数据为样本集,以最小二乘支持向量回归(least squares support vector regression,LS-SVR)算法建立了小麦冠层叶绿素含量反演模型,其校正决定系数C-R2与预测决定系数P-R2分别为0.751与0.722,在各指数中反演精度最高。进一步分析表明,REP对叶绿素含量以及 LAI值较高与较低的样本均具备良好的预测能力,可有效避免样本取值范围以及冠层郁闭度等因素对叶绿素含量估测的影响。利用LS-SVR反演模型完成了OMIS影像叶绿素含量的遥感填图,并以地面实测值进行检验,其拟合模型R2与RMSE值分别为0.676与1.715。结果表明,高光谱指数REP所建立的LS-SVR模型实现了叶绿素含量的准确估测,可用于小麦叶绿素含量信息的快速、无损获取。

[ Liang L, Yang M H, Zhang L P, et al.Chlorophyll content inversion with hyperspectral technology for wheat canopy based on support vector regression algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012,28(20):162-171,294. ]

[8]
王家强,伍维模,李志军,等.基于高光谱指数的天然胡杨叶绿素遥感建模研究[J].干旱区资源与环境,2014,28(10):95-99

[ Wang J Q, Wu W M, Li Z J, et al.Remote sensing modeling research of natural Populus Euphratica chlorophyll based on hyperspectral indices[J]. Journal of Arid Land Resources and Environment, 2014,28(10):95-99. ]

[9]
艾金泉,陈文惠,陈丽娟,等.冠层水平互花米草叶片光合色素含量的高光谱遥感估算模型[J].生态学报,2015,35(4):1175-1186以闽江河口鳝鱼滩湿地互花米草(<i>Spartina alterniflora</i>)的实测冠层高光谱反射率和叶片光合色素含量(LPPC)为数据源,在分析LPPC与原始光谱反射率、一阶导数光谱反射率、22种已报道光谱指数和14种新构建的植被指数相关性的基础上,利用直线回归、指数回归、对数回归以及乘幂回归方法,系统地比较了36种植被指数在估算互花米草LPPC中的表现。研究表明:(1)一阶导数光谱反射率组合的植被指数用于估算互花米草的LPPC优于原始光谱反射率;(2)红边区域一阶导数光谱是估测互花米草LPPC的最佳波段;(3)对于单一色素含量的估算,叶绿素a (Chla)的最佳估算指数为FDNDVI[723,703];叶绿素b (Chlb)的最佳估算指数为FDRVI[723,525];类胡萝卜素(Cars)的最佳估算指数为FDNDVI[723,703];(4)对于使用统一参量同时估算Chla、Chlb、Cars,由FDRVI[723,703]建立的对数估算模型效果最佳。研究成果可为湿地植物生化参量反演提供参考,也可为闽江河口湿地入侵种互花米草的动态监测和生态评估管理提供有力的科学依据。

DOI

[ Ai J Q, Chen W H, Chen L J, et al.Hyperspectral remote sensing estimation models for foliar photosynthetic pigment contents at canopy level in an invasive species[J]. Acta Ecologica Sinica, 2015,35(4):1175-1186. ]

[10]
肖汉,陈秀万,杨振宇,等.基于光谱分析的草地叶绿素含量估测植被指数[J].光谱学与光谱分析,2014,34(11):3075-3078对现有叶绿素遥感估测研究方法进行比较, 确定植被指数法是其中最实用、 普适性最强的研究方法。 近年来, 草地退化问题日益严峻, 需要进一步从光谱分析、 植物生化参数估测的角度加以研究, 因而亟需建立一种用于反演草地植被叶绿素含量的植被指数。 首先对四川省松潘草原和内蒙古自治区贡格尔草原的草地实测反射率光谱曲线及其一阶微分曲线进行分析, 通过这两种光谱与叶绿素含量的相关性分析, 找到红边区域(red-edge position, REP)与草地叶绿素含量之间的规律, 即叶绿素含量越高, 反射率一阶微分曲线的红边拐点(red-edge inflection point, REIP)取值越高, 由此构建草地叶绿素含量估测植被指数(grassland chlorophyll index, GCI), 选取最适宜反演的波段, 最后采用卫星高光谱影像计算GCI, 将计算结果与野外试验观测的叶绿素含量数据进行精度分析验证。 结果证明, 对于草地叶绿素含量来说, GCI比其他叶绿素指数的敏感性更强, 具有较高的草地叶绿素含量估测精度。 GCI是第一个针对草地叶绿素含量估测而被提出的植被指数, 其对遥感反演草地叶绿素含量具有广泛应用潜力。 同时这种基于光谱分析的草原植被叶绿素含量估测方法为其他的草原植被生化参数估测、 草原植被生长状况评价以及草地生态环境变化大面积监测提供了新的研究思路。

[ Xi H, Chen X W, Yang Z N, et al.Vegetation index estimation by chlorophyll content of grassland based on spectral analysis[J]. Spectroscopy and Spectral Analysis, 2014,34(11):3075-3078. ]

[11]
梁爽,赵庚星,朱西存.苹果树叶片叶绿素含量高光谱估测模型研究[J].光谱学与光谱分析,2012,32(5):1367-1370叶片叶绿素含量是评估果树长势和产量的重要参数, 实现快速、 无损、 精确的叶绿素含量估测具有重要意义。 本研究以山东农业大学苹果园为试验区, 采用高光谱分析技术探索苹果树叶片叶绿素含量的估测方法。 通过分析叶片高光谱曲线特征, 对原始光谱分别进行一阶微分、 红边位置以及叶面叶绿素指数(LCI)变换, 分别将其与叶绿素含量进行相关分析及回归分析, 建立叶绿素含量估测模型并进行检验, 从中筛选出精度最高的模型。 结果显示, 以LCI为变量的估测模型以及以一阶微分521和523 nm组合为变量的估测模型拟合精度最高, 决定系数R2分别为0.845和0.839, 均方根误差RMSE分别为2.961和2.719, 相对误差RE%分别为4.71%和4.70%。 因此LCI及一阶微分是估测苹果树叶片叶绿素含量的重要指标。 该模型对指导苹果树栽培生产具有积极意义。

[ Li S, Zhao G X, Zhu X C.Hyperspectral estimation models of chlorophyll content in apple leaves[J]. Spectroscopy and Spectral Analysis, 2012,32(5):1367-1370. ]

[12]
张永贺,郭啸川,褚武道,等.基于红边位置的木荷叶片叶绿素含量估测模型研究[J].红外与激光工程,2013,42(3):798-804利用红边参数反演作物参数是定量遥感研究的一个热点, 红边参数中红边位置与作物生化组分强相关, 为监测作物胁迫提供了一个非常敏感的指标。准确估测植被叶绿素含量,对于研究森林健康和胁迫、森林生产力的估计, 碳循环的研究有着重要的意义。介绍几种红边位置算法, 并对这些算法及其应用进行了比较,通过选取红边位置的不同敏感波段来估测植被叶片叶绿素含量。经室内光谱获取叶片的光谱数据,采用一阶光谱导数法、平滑处理后一阶光谱导数法、线性四点内插法、五次多项式拟合法四种算法处理光谱数据,获得红边位置变量,并与叶绿素含量进行拟合,构建估测木荷叶片叶绿素含量的回归模型。结果表明:各种算法获取的红边位置变量所构建的回归模型估测叶绿素含量是可行的;五次多项式拟合法估算精度是最高的,其获取红边位置计算相对复杂;线性四点内插法估算精度次之,但计算较简便。

[ Zhang Y H, Guo X C, Chu W D, et al.Estimation model of schima superba leaf chlorophyll content based on red edge position[J]. Infrared and Laser Engineering, 2013,42(3):798-804. ]

[13]
杜华强,葛宏立,范文义,等.马尾松针叶光谱特征与其叶绿素含量间关系研究[J].光谱学与光谱分析,2009,29(11):3033-3037以马尾松针叶野外高光谱数据为基础, 分析了马尾松光谱变化, 构建或借助不同光谱特征参数, 在理论和实践分析的基础上, 建立了马尾松针叶叶绿素含量与光谱反射率及9个特征参数之间的关系。 研究结果表明: (1) 马尾松叶绿素含量在527, 703, 1 364及1 640 nm四个波长附近, 与其反射率具有较好的线性关系, 为马尾松遥感监测在波段选择上提供了依据; (2)红边位置、 红边平均反射率、 红边位置附近平均反射率、 红边斜率、 红边面积、 红谷吸收深度、 绿峰反射高度、 红边归一化植被指数、 红边植被胁迫指数等9个马尾松反射光谱特征参数均与叶绿素含量间存在指数函数关系, 相关系数绝对值在0.5~0.7之间; (3) 采用9个光谱特征参数建立了马尾松针叶叶绿素含量预测模型, 且所建立的基于高斯核函数变换的偏最小二乘回归模型对叶绿素含量的预测精度远远大于传统线性回归模型, 模型的均方误差为0.0088, 平均绝对百分误差为0.7617%。

[ Du H Q, Ge H L, Fan W Y, et al.Study on relationships between total chlorophyll with hyperspectral features for leaves of pinus massoniana forest[J]. Spectroscopy and Spectral Analysis, 2009,29(11):3033-3037. ]

[14]
孙小涛,周忠发,陈全,等.重点生态功能区水土流失敏感性评价与分布研究——以贵州省雷山县为例[J].水土保持学报,2016,30(6):73-78

[ Sun X T, Zhou Z F, Chen Q, et al.Sensitivity evaluation and the spatial distribution of soil erosion in keyecological function areas: A case of Leishan city of Guizhou[J]. Journal of Soil and Water Conservation, 2016,30(6):73-78. ]

[15]
Chen Z, Chen Z, Bai L, et al.A catastrophe model to assess soil restoration under ecological restoration in the red soil hilly region of china[J]. Pedosphere, 2017,27(4):778-787.

[16]
Chen Z, Chen Z, Yan X, et al.Stoichiometric mechanisms of dicranopteris dichotoma, growth and resistance to nutrient limitation in the Zhuxi watershed in the red soil hilly region of China[J]. Plant and Soil, 2016,398(1-2):367-379.

[17]
陈志强,陈志彪.南方红壤侵蚀区土壤肥力质量的突变——以福建省长汀县为例[J].生态学报,2013,33(10):3002-3010以南方红壤侵蚀区典型区域福建省长汀县为研究区,将土壤肥力质量10个因子作为内部因子,坡度、植被覆盖度、水土流失强度等作为外部因子,构建土壤肥力质量演变的尖点突变模型,并分析土壤肥力质量演变分别与土壤肥力质量等级、水土流失强度、坡度和植被覆盖度的关系。研究结果表明:1)90个样点中,突变的样点共27个,占30%,稳定的样点共63个,占70%;2)突变样点主要对应于土壤肥力质量等级1和等级2(分别占突变样点总数的48.15%和33.33%)、水土流失微度和轻度(分别占37.04%和44.44%)、坡度5-10°和10-15°(分别占37.04%和40.74%)、植被覆盖度>0.4和0.3-0.4(分别占48.15%和37.04%);分叉集Δ与土壤肥力质量等级、水土流失强度、坡度和植被覆盖度都呈极显著(<em>P</em><0.01)相关关系,皮尔逊相关系数绝对值的大小顺序为:水土流失强度>植被覆盖度>土壤肥力质量等级>坡度;突变主要发生于土壤肥力质量等级较高、水土流失强度较轻、坡度中等、植被覆盖较好的地点;3)土壤肥力质量演变时间相对较长,应根据中间过渡状态来判断是否产生突变;土壤肥力质量处于突变状态时可用较小投入产生较大效益,在关注严重水土流失区生态恢复与重建的同时,不应忽视突变区的治理。

DOI

[ Chen Z Q, Chen Z B.Mutation of soil fertility quality in the red eroded area of southern China: A case study in Changting County, Fujian Province[J]. Acta Ecologica Sinica, 2013,33(10):3002-3010. ]

[18]
褚武道,陈文惠,艾金泉,等.基于偏最小二乘法的樟树叶片叶绿素含量高光谱估算[J].福建师范大学学报(自然科学版),2014,30(1):65-70,90

[ Chu W D, Chen W H, Ai J Q, et al.Hyperspectral estimation of camphor tree leaf chlorophyll content based on partial least squares regression[J]. Journal of Fujian Normal University (Natural Science Edition), 2014,30(1):65-70,90. ]

[19]
宫兆宁,赵雅莉,赵文吉,等.基于光谱指数的植物叶片叶绿素含量的估算模型[J].生态学报,2014,34(20):5736-5745叶绿素是光合作用能力和植被发育阶段的指示器,是监测湿地植被生长健康状况的重要指标之一;高光谱遥感技术可以为植物叶绿素含量的定量化诊断提供简便有效、非破坏性的数据采集和处理方法。为保证被探测叶片面积相同,消除背景反射、叶片表面弯曲造成的光谱波动及叶片内部变异造成的影响,研究采用Field Spec 3光谱仪加载手持叶夹式叶片光谱探测器,测定野鸭湖湿地典型植物的叶片高光谱反射率数据,同时通过分光光度计室内测定相应叶片的叶绿素含量。采用相关性及单变量线性拟合分析技术,建立二者的关系模型,包括叶绿素含量与"三边"参数的相关模型以及比值光谱指数(SR)模型和归一化差值光谱指数(ND)模型,并采用交叉检验中的3K-CV方法对估算模型进行模型精度检验。结果表明:植物叶片叶绿素含量与"三边"参数大多都呈极显著相关,相关系数最大达到0.867;计算光谱反射率组成的比值(SR)和归一化(ND)光谱指数与叶绿素含量的决定系数,总体相关性比较高,较好的波段组合均为550-700nm与700-1400nm以及550-700nm与1600-1900nm,与叶绿素含量相关性最好的指数分别是SR(565nm,740nm)和ND(565nm,735nm)。并通过选取相关性最佳的光谱特征参数,分别基于"三边"参数和ND模型指数构建了植物叶片叶绿素含量的估算模型。其中,基于红边位置(WP_r)光谱特征参数和ND(565nm,735nm)光谱指数建立的叶绿素含量估算模型,取得了较好的测试效果,检验拟合方程的决定系数(<em>R</em><sup>2</sup>)都达到0.8以上,估算模型分别为<em>y</em>=0.113<em>x</em>-78.74,<em>y</em>=5.5762<em>x</em>+4.4828。通过3K-CV方法进行测试和检验,植物叶绿素含量估算模型均取得了较为理想的预测精度,预测精度的分别为93.9%及90.7%。高光谱遥感技术对植被进行微弱光谱差异的定量分析,在植被遥感研究与应用中表现出强大优势,为植物叶绿素含量诊断中的实际应用提供了重要的理论依据和技术支持。

DOI

[ Gong Z N, Zhao Y L, Zhao W J, et al.Estimation model for plant leaf chlorophyll content based on the spectral index content[J]. Acta Ecologica Sinica, 2014,34(20):5736-5745. ]

[20]
李春昉,郭际,赵绍丰.多元统计分析之因子分析浅析[J].价值工程,2010,29(15):128-129

[ Li C F, Guo J, Zhao S F [ The analysis of factor analysis in multivariate statistical analysis[J]. Value Engineering, 2010,29(15):128-129. ]

[21]
Horler D N H, Dockray M, Barber J. The red edge of plant leaf reflectance[J]. International Journal of Remote Sensing, 1983,4(2):273-288

[22]
姚付启,张振华,杨润亚,等.基于红边参数的植被叶绿素含量高光谱估算模型[J].农业工程学报,2009,25(S2):123-129

[ Yao F Q, Zhang Z H, Yang R Y, et al.Hyperspectral models for estimating vegetation chlorophyll content based on red edge parameter[J]. Transactions of the CSAE, 2009,25(s2):123-129. ]

[23]
Tucker C J.Red and photographic infrared linear combinations for monitoring vegetation[J]. Remote Sensing and Environment, 1979,8(2):127-150.

[24]
Gitelson A, Merzlyak M N.Spectral reflectance changes associated with autumn senescence of Aesculus hippocastanum L. and Acer platanoides L. leaves. spectral features and relation to chlorophyll estimation[J]. Journal of Plant Physiology, 1994,143(3):286-292.

[25]
Datt B.A new reflectance index for remote sensing of chlorophyll content in higher plants: Tests using Eucalyptus leaves[J]. Journal of Plant Physiology, 1999,154(1):30-36.

[26]
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[J]. Remote Sensing of Environment, 2003,84(4):526-537.

[27]
Peñuelas J, Gamon J A, Fredeen A L, et al.Reflectance indices associated with physiological changes in nitrogen and water-limited sunflower leaves[J]. Remote Sensing of Environment, 1994,48(2):135-146.

[28]
Merzlyak M N, Chivkunova O B, Gitelson A A, et al.Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening[J]. Physiologia Plantarum, 2010,106(1):135-141.

[29]
Blackburn G A.Quantifying chlorophylls and caroteniods at leaf and canopy scales: An evaluation of some hyperspectral approaches[J]. Remote Sensing of Environment, 1998,66(3):273-285.

[30]
Vogelmann J E, Rock B N, Moss D M.Red edge spectral measurements from sugar maple leaves[J]. International Journal of Remote Sensing, 1993,14(8):1563-1575.

[31]
Haboudane D, Miller J R, Tremblay N, et al.Integrated narrow-band vegetation indices for prediction of crop chlorophyll content for application to precision agriculture[J]. Remote Sensing of Environment, 2002,81(2-3):416-426.

[32]
焦全军,张霞,张兵,等.基于叶片光谱的森林叶绿素浓度反演研究[J].国土资源遥感,2006,18(2):26-30<p>准确估测森林叶片叶绿素浓度有助于进一步理解和模拟森林生态系统。基于江西省千烟州试验区主要树种光谱数据和相应叶绿素浓度等化学参量数据,使用高光谱指数和偏最小二乘(PLS)回归方法进行森林叶绿素浓度的反演,对具代表性的几种叶绿素相关的光谱指数进行反演能力评价,构造出表征叶绿素吸收特征的叶绿素吸收面积指数(CAAI),发现TCARI(改进型的叶绿素吸收比率指数)和CAAI能较好地估测试验区森林主要树种的叶绿素浓度。针对叶片生化参量之间强相关现象,首先使用了叶片生化参量吸收特征分析方法,选定特定波段。然后对叶绿素浓度进行PLS回归估测,并从PLS获得估测模型的结果来评价和解释几种高光谱指数的叶绿素反演能力。</p>

DOI

[ Jiao Q J, Zhang X, Zhang B, et al.The petrieval of forest chlorophyll concentration based on foliar spectra[J]. Remote Sensing for Land and Resources, 2006,18(2):26-30. ]

文章导航

/