作者简介: 李 哲, 女, 1993年生, 新疆大学资源与环境科学学院硕士研究生 e-mail: 13999940867@163.com
叶片叶绿素能够有效监测植被的生长状况, 利用光谱指数反演植被叶绿素含量是目前的通用方法。 实测了盐生植物光谱反射率和叶片叶绿素含量。 对SPAD值进行变换, 对比Pearson与VIP方法探讨盐生植被叶片叶绿素含量与植被指数的相关性并进行精度验证, 从中选出最佳拟合模型。 研究表明, 通过对Pearson与VIP相关性分析, 最终选定VIP方法建立植被指数的叶片叶绿素估算模型, NDVI705, ARVI, CIred edge, PRI, VARI, PSRI和NPCI的VIP值均大于0.8, 因此选定这七个植被指数为最优植被指数; 预测结果显示, 所有模型的相关性都在0.7以上, 预测值与实测值相关性最好的是经过倒数变换的SPAD值, R=0.816, RMSE=0.007。 基于VIP方法的反演模型能较好地估算研究区植被叶绿素含量, 该方法为植物叶绿素含量诊断的实际应用提供了重要的理论依据和技术支持。
The chlorophyll can effectively monitoring vegetation growth status, currently the hyperspectral vegetation index (VI) is one of the common methods that have been widely applied to estimate the leaf chlorophyll content (LCC) inversion. Non-destructive rapid estimation of chlorophyll using hyperspectral remote sensing technology is a prerequisite to dynamically monitor chlorophyll content, and it is an important research issue of vegetation remote sensing. The author measured the leaf spectral reflectance and chlorophyll relative content of desert plants, analyzed the spectral curves of different desert plants under the same chlorophyll content, then transformed the SPAD value, compared Pearson and VIP methods to study the correlation between chlorophyll content and vegetation index of desert vegetation. Finally, the author selected the best fitting model from accuracy test. The results showed that: Based on the comparative analysis between Pearson and VIP, established the chlorophyll content estimation model by VIP method, selected 7 vegetation indices, which was NDVI705, ARVI, CIred edge, PRI, VARI, PSRI and NPCI respectively, the value of VIP all greater than 0.8, thus these 7 vegetation indices were the optimal vegetation indices. The prediction results indicated that the correlation of all models was more than 0.7, the best correlation between the predicted value and the measured value was the SPAD value of the reciprocal transformation, R=0.816, RMSE=0.007. The inversion model based on VIP method can estimate the chlorophyll content of vegetation in the study area, it provides an important theoretical basis and technical support for the practical application in the diagnosis of plant chlorophyll content.
叶绿素是植物生长过程的一个重要生化参数, 是植物光合作用能力的指示器, 已被广泛用于植物的养分循环、 发育情况、 生产力和疾病监测等研究[1], 叶片叶绿素含量(leaf chlorophyll content, LCC)能够反映植物的生长状况[2], 因此, 在生态系统较为脆弱的西北荒漠地带, 监测植被叶片的叶绿素含量对干旱区生态系统的恢复与重建和生态系统生产力的评估具有重要指导意义。
Le Maire等[3]总结了2002年之前所有发表的有关叶绿素含量预测的光谱指数。 在电磁波的蓝光和红光范围, 叶绿素表现为强烈吸收, 光谱指数基于那些对叶绿素敏感而对其他因素不敏感的窄波段光谱比值构建。 应用于叶片尺度的植被光谱指数大多数结构相对简单, 常用比值型、 差值型、 归一化差值型、 新二重差值型及改进版本等。 较为典型的有改进的叶绿素吸收指数MCARI和土壤调节指数, 其优点在于消除冠层参数, 如叶面积指数LAI和土壤背景变化等影响[4]。 Dash等[5]针对高浓度的植被叶绿素含量, 发展了新型植被指数MTCI; Daughtry等[6]根据前人经验, 并综合OSAVI指数的优点, 构建MCARI指数的同时, 进一步提出了MCARI/OSAVI和TCARI/OSAVI指数, 提高了反演叶绿素的准确性; 李鑫川等利用可见光-近红外与短波红外的特征建立了新型植被指数来估算农作物的叶面积指数; 廖钦洪等[7]有效利用遥感影像波段对各观测角的玉米冠层反射率进行分析, 对比新型植被指数HD-TCARI和TCARI反演叶绿素的特点。
由于受西北干旱区水盐胁迫的影响, 使得研究区的典型植被如柽柳(Tamarix)、 梭梭(Haloxylon ammodendron)、 白刺(Nitraria tangutorum)等与普通的植被相比形态和生理特征都具有差异性。 本工作选取多种已发表的叶绿素植被指数, 通过数学变换手段研究这些指数在盐生植被特定季节的适应性, 并对结果进行分析评价, 从中筛选出适应荒漠植被特定季节的几种光谱指数。
艾比湖湿地自然保护区位于新疆博尔塔拉蒙古自治州(44° 30'— 45° 09'N, 82° 36'— 83° 50'E), 与精河县、 阿拉山口市和博乐市相邻。 总面积2 670.85 km2, 其中水域面积为683.81 km2。 艾比湖湿地是我国典型的沙漠湖泊湿地, 湖边植物具有良好的防风固沙的特点, 保护区夏季干旱少雨, 冬季寒冷, 年日照时数大于2 500 h, 年均温度8.3 ℃, 年均降水90 mm, 蒸发量高于3 000 mm, 属于中温带大陆性干旱气候[8]。 保护区内有各类盐生植物占全国盐生植物种类的60%左右, 物种多样性低, 优势种主要为耐盐碱性的梭梭、 柽柳、 盐穗木(Halostachys caspica)、 盐节木(Halocnemum strobilaceum)、 花花柴(Karelinia caspica)、 芦苇(Phragmites australis)等。
于2016年10月对艾比湖湿地国家自然保护区进行实地考察, 保护区布设35个样点, 其中每个样点采集两种荒漠植被。 用新鲜的冰块对采集回来的样本进行冷却保存, 并放入密封的塑料袋, 确保植被样本不被损坏。 使用便携式地物光谱(Field Spec HiRes, ASD, 美国)对植物光谱曲线进行采集, 该仪器光谱范围为350~2 500 nm, 光谱分辨率3 nm。 测量前用白板定标, 仪器距样点上方5 cm, 每个点位进行5次反射光谱测量, 最后取其平均值作为该点的原始光谱反射率。 然后利用光谱仪的后处理软件分析采集的植物光谱曲线, 将每个点号的5条光谱数据置于Process下Statistics内的Mean命令进行计算平均值, 然后在Origin8.0软件中进行微分变换与多元统计分析。
采用日本KONICA MINOLTA SPAD-502叶绿素计对盐生植物叶片(梭梭、 胡杨(Populus euphratica)、 柽柳、 芦苇、 白刺、 花花柴、 盐节木)的相对叶绿素含量快速无损监测。 通过测量叶片的吸收率计算出SPAD值, 得到的数据应与叶片内部结构的叶绿素浓度成正比[9]。 通常情况下, 叶片叶绿素浓度的相对值就是SPAD的值。 盐生植物叶片的SPAD值与光谱测量同步进行, 每个样点选取2种植被, 每种四棵植株, 每株植被选择健康未枯萎, 并且是第2、 3片叶片进行测量; 每叶片测得10个SPAD值, 最后取其平均值作为该叶片的SPAD值。
研究[10]表明, 光谱指数和叶绿素含量之间的经验模型是田间和实验室测量叶绿素含量进行验证的有效估算方法。 因此, 选取了15个前人构建的高光谱植被指数(表1)应用于荒漠植被叶片叶绿素含量的反演研究。
![]() | 表1 光谱指数计算公式 Table 1 Formula for the spectral index |
Pearson相关系数, 通常用来度量两个变量X和Y之间的相互关系, 计算公式为[17]
式中, Cov(X, Y)代表X与Y的协方差, Var(X)和Var(Y)代表X和Y的方差。 当相关系数为1时, X与Y的关系为Y=aX+b, a> 0; 当相关性为-1时, X与Y的关系为Y=aX+b, a< 0。 如果X与Y相互独立, 相关系数为0。
通过计算变量投影重要性指标VIP来分析各高光谱植被指数对盐生植被叶片叶绿素浓度预的贡献。 VIP的计算公式为
式中,
将SPAD值(R)进行了LOG(R), SQRT(R), 1/R, LOG(1/R), SQRT(1/R), 1/LOG(R)六种数学变换, 对比两种相关分析方法以研究其相关性, 然后选用相关性最好的一种方法建立叶片叶绿素含量估算模型(n=70), 并随机选取23个独立样本进行预测验证。 首先采用Pearson相关分析分析研究多种变换下SPAD值与植被指数的相关性(见图1)。
由图1可知, 在未变换的Pearson相关矩阵图中, 与SPAD值呈现较好相关性的共有4个植被指数, 相关性由大到小分别是: PSRI(-0.597* * ), NPCI(-0.535* * ), VARI(0.494* )和ARVI(0.421* ), 其中PSRI和NPCI在0.01水平上显著相关, VARI和ARVI在0.05水平上显著相关。 在其余变换的Pearson相关矩阵图中, 与SPAD值呈现较好相关性的共有5个植被指数, 分别是PSRI, NPCI, VARI, ARVI和PRI。 在LOG(R)变换的Pearson相关矩阵图中, 相关性由大到小分别是: PSRI(-0.634* * ), NPCI(-0.577* * ), VARI(0.562* * ), ARVI(0.441* )和PRI(-0.433* ), 其中PSRI, NPCI, VARI在0.01水平上显著相关, ARVI和PRI在0.05水平上显著相关。 在SQRT(R)变换的Pearson相关矩阵图中, 相关性由大到小分别是: PSRI(-0.620* * ), NPCI(-0.560* * ), VARI(0.531* * ), ARVI(0.435* )和PRI(-0.426* ), 其中PSRI, NPCI, VARI在0.01水平上显著相关, ARVI和PRI在0.05水平上显著相关。 在1/R变换的Pearson相关矩阵中, 相关性由大到小分别是: PSRI(0.629* * ), VARI(-0.595* * ), NPCI(0.585* * ), ARVI(-0.425* )和PRI(0.418* ), 其中PSRI, VARI和NPCI在0.01水平上显著相关, ARVI和PRI在0.05水平上显著相关。 在LOG(1/R)变换Pearson相关矩阵图中, 相关性由大到小分别是: PSRI(0.634* * ), NPCI(0.577* * ), VARI(-0.562* * ), ARVI(-0.441* )和PRI(0.433* ), 其中PSRI, VARI和NPCI在0.01水墙上显著相关, ARVI和PRI在0.05水平上显著相关。 在SQRT(1/R)变换Pearson相关矩阵图中, 相关性由大到小分别是: PSRI(0.637* * ), VARI(0.585* * ), NPCI(-0.583* * ), ARVI(-0.437* )和PRI(0.430* ), 其中PSRI, VARI和NPCI在0.01水平上显著相关, ARVI和PRI在0.05水平上显著相关。 在1/LOG(R变换Pearson相关矩阵图中, 相关性由大到小分别是: PSRI(0.635* * ), VARI(-0.585* * ), NPCI(0.585* * ), ARVI(-0.434* )和PRI(0.427* ), 其中PSRI, VARI和NPCI在0.01水平上显著相关, ARVI和PRI在0.05水平上显著相关。
通过Pearson相关性分析, 筛选出了五个植被指数, 分别是PSRI, NPCI, VARI, ARVI和PRI。 显著性在0.01水平的植被指数分别是PSRI, VARI和NPCI, 表明该三种植被指数与SPAD值的相关性较好。 将变换前后的SPAD值与筛选的植被指数用SPSS软件进行多元回归线性分析, 得到如下结果(见表2)。
![]() | 表2 SPAD值与植被指数拟合评定指标(Pearson) Table 2 The fit evaluation indicators of the SPAD with vegetation index (Pearson) |
笔者在Pearson相关分析评定指标结果的基础上, 又采用了VIP相关分析研究多种变换下SPAD值与植被指数的相关性(见图2)。
基于光谱测量数据和PLSR模型的VIP方法计算分析结果, 显示不同SPAD变换应用中, 各植被光谱指数在估算盐生植物叶片叶绿素时的重要性(图2)。 利用VIP方法进行重要性评价时, 当VIP值大于1时, 说明该植被指数在估算植物叶片叶绿素时具有重要作用; 当植被指数的VIP值在0.5~1之间时, 说明该植被指数的预测重要性还需要进一步分析, 需要通过增加样本或根据其他条件进行判断分析, 而一般情况下, VIP值大于0.8的对应植被指数是预测叶片叶绿素值的重要植被指数[19]。
从图2中可以看出, NDVI705, ARVI, CIred edge, PRI, VARI, PSRI和NPCI这七个指数的VIP值始终大于0.8, 其中ARVI, CIred edge, PRI, VARI, PSRI和NPCI的值始终大于1, 说明这些植被指数在估算植物叶片叶绿素时具有重要作用。 因此, 根据VIP图选取了不同的植被指数, 将变换前后的SPAD值与筛选的植被指数用SPSS软件进行多元回归线性分析, 得到如表3所示结果。
![]() | 表3 SPAD值与植被指数拟合评定指标(VIP) Table 3 The fit evaluation indicators of the SPAD with vegetation index (VIP) |
由表3得知, 1/R, SQRT(1/R)和1/(LOGR)这三种SPAD变换与对应筛选的植被指数相关性最好, 相关性达到0.8, 根据R最大, RMSE最小的原则, 相关性最优的叶绿素变换是SPAD倒数变换, 相关性最小的是未经变换的SPAD值, R仅为0.73, RMSE最大达到11.28。
通过对比Pearson与VIP的SPAD值与植被指数的拟合评定指标表, 最终采用VIP筛选的植被指数建立叶片叶绿素的估算模型(见表4)。
![]() | 表4 SPAD值与植被指数最佳估算模型 Table 4 The optimal estimation model of SPAD and vegetation index |
为了进一步验证模型的准确性和稳定性, 采用预测集的23个样本对叶绿素含量估算模型的预测能力进行验证(图3)。
![]() | 图3 预测值与实测值的散点图Fig.3 Scatter plot between the observed LCC and the predicted values with different vegetation indices |
从图3中可知, 经过倒数变换的SPAD值其预测值与实测值相关性最好, R达到0.816, RMSE为0.007, 所有模型的相关性都在0.7以上, 模型相对较差的是未经变换的SPAD值, R为0.742, RMSE为6.36。 可以看出, 基于不同植被指数的反演模型都能较好地得到研究区植被叶绿素含量, 经过验证发现这七种模型对盐生植被叶绿素含量的预测是可行的。
采用野外实测光谱及其室内实测叶绿素数据, 对SPAD值进行数学变换, 探索高光谱遥感识别植被叶片叶绿素含量的敏感波段, 通过对比验证Pearson与VIP模型的优劣, 探索盐生植被叶绿素含量与植被指数的相关性, 建立最佳估算模型, 并对模型进行了检验和预测, 这是一个新的尝试。 本研究仍有些许不足之处, 受时间和自然条件限制, 采集的样本数量较少, 需要适当增加样本数量。 结果利用利用高光谱植被指数反演秋季的盐生植被叶绿素是可行的, 但对于全年的植被叶绿素含量反演仍需要进一步的研究探讨, 以提高光谱指数的普适性。
通过对比作者最终采用VIP筛选的植被指数建立叶片叶绿素的估算模型, 筛选了NDVI705, ARVI, CIred edge, PRI, VARI, PSRI和NPCI七个指数。 基于VIP方法的SPAD值与植被指数的预测模型相关性都在0.7以上, 根据R最大RMSE最小的原则, 预测值与实测值相关性最好的是经过倒数变换的SPAD值, R达到0.816, RMSE为0.007。 因此基于VIP方法的反演模型被认为是最佳估算荒漠植被叶片叶绿素含量模型。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|