中国生态农业学报  2017, Vol. 25 Issue (8): 1224-1235  DOI: 10.13930/j.cnki.cjea.170112
0

引用本文 

郑雯, 明金, 杨孟克, 周四维, 汪善勤. 基于波段深度分析和BP神经网络的水稻色素含量高光谱估算[J]. 中国生态农业学报, 2017, 25(8): 1224-1235. DOI: 10.13930/j.cnki.cjea.170112
ZHENG W, MING J, YANG M K, ZHOU S W, WANG S Q. Hyperspectral estimation of rice pigment content based on band depth analysis and BP neural network[J]. Chinese Journal of Eco-Agriculture, 2017, 25(8): 1224-1235. DOI: 10.13930/j.cnki.cjea.170112

基金项目

国家高技术研究发展计划(863计划)项目(2013AA102401-3)资助

通讯作者

汪善勤, 主要研究方向为农业遥感和数字土壤。E-mail:sqwang@mail.hzau.edu.cn

作者简介

郑雯, 研究方向为农业遥感。E-mail:trista_only@126.com

文章历史

收稿日期:2017-02-10
接受日期:2017-03-30
基于波段深度分析和BP神经网络的水稻色素含量高光谱估算*
郑雯1, 明金1, 杨孟克1, 周四维1, 汪善勤1,2     
1. 华中农业大学资源与环境学院 武汉 430070;
2. 农业部长江中下游耕地保育重点实验室 武汉 430070
摘要:该文以水稻田间氮肥水平试验为基础,采用单变量的线性和非线性回归方法,建立基于植被指数的水稻色素含量高光谱估算模型。各植被指数对色素含量的估计能力分析结果显示,植被指数在色素含量较大时存在饱和问题,为此尝试将波段深度分析(BDA)与BP神经网络结合,以提高利用高光谱技术对水稻叶片色素含量的估算精度。基于连续统去除处理的水稻冠层高光谱数据(400~750 nm),选取波段深度(BD)、波段深度比(BDR)、归一化波段深度(NBDI)和归一化面积波段指数(BNA)4种波段指数,在此基础上进行主成分分析(PCA)实现降维,然后采用反向传播(BP)神经网络方法对水稻叶片色素含量进行高光谱反演,探讨BDA与BP神经网络结合解决植被指数饱和问题的可能性和有效性。结果表明,波段深度分析突出了光谱吸收特征差异,挖掘了更多的潜在信息,使得光谱曲线的差异性得到增强。BD与BP结合的估算模型对水稻叶片中的类胡萝卜素含量估算精度最高(R2=0.61,RMSEP=0.128 mg·g-1),BNA与BP结合的估算模型对水稻叶片中的叶绿素含量估算精度最高(R2=0.73,RMSEP=0.343 mg·g-1)。对比分析BDA与BP结合的模型和植被指数最佳回归模型的精度,发现波段深度分析建立的BP神经网络模型能较好地解决饱和问题,提高水稻叶片色素含量的估算精度。
关键词:高光谱    水稻    色素    植被指数    波段深度分析    主成分分析    反向传播神经网络    
Hyperspectral estimation of rice pigment content based on band depth analysis and BP neural network*
ZHENG Wen1, MING Jin1, YANG Mengke1, ZHOU Siwei1, WANG Shanqin1,2     
1. College of Resources and Environment, Huazhong Agricultural University, Wuhan 430070, China;
2. Key Laboratory of Arable Land Conservation in Middle and Lower Reaches of Yangtze River, Ministry of Agriculture, Wuhan 430070, China
*This study was supported by the National High-tech R & D Program of China (863 Program) (2013AA102401-3)
** Corresponding author, WANG Shanqin. E-mail:sqwang@mail.hzau.edu.cn
Received Feb. 10, 2017; accepted Mar. 30, 2017
Abstract: The estimation accuracy of plant pigment content is low under higher pigment content since conventional vegetation indices tend to be less sensitive to the variance of pigment content. In order to improve estimation accuracy of rice carotenoid and chlorophyll contents with canopy reflectance during all growth stage, we explore the feasibility and effectiveness of combining the band depth analysis (BDA) and back propagation (BP) neural network to solve the problem of vegetation index saturation. With canopy hyperspectral data (400-750 nm), four band indices — band depth (BD), band depth ratio (BDR), normalized band depth index (NBDI) and band depth normalized to band area (BNA) — were calculated via continuum removal processing. Principal component analysis (PCA) was used to reduce the dimensions of hyperspectral data, and determined 10 principle components, which were introduced into BP neutral network as input variables. In the study, canopy hyperspectral reflectance and pigment content measurements were conducted in Meichuan Town of Hubei Province, China. Eight treatments of nitrogen fertilization (0, 45, 82.5, 127.5, 165, 210, 247.5 and 292.5 kg·hm-2) were applied to generate various indices of vegetation and pigment content. Linear and nonlinear regression models were used to quantitatively analyze the vegetation indices and measured pigment content. In addition, coefficient of determination (R2) and root mean square error (RMSE) were used to evaluate the models. All the hyperspectral indices were comparatively analyzed. As a result, BDA showed the differences in spectral absorption characteristics and revealed more potential information to enhance spectral difference. The estimation model combined band index BD and BP had the highest estimation accuracy for carotenoid content in rice leaves, with R2 = 0.61 and RMSE = 0.128 mg·g-1; while the estimation model combined band index BNA and BP had the highest estimation accuracy for chlorophyll content in rice leaves, with R2 = 0.73 and RMSE = 0.343 mg·g-1. Further comparison between BDA & BP models with the best regression model for vegetation index indicated that BP neutral network model based on BDA provided a better solution to saturation problem and a higher estimation precision of rice leaf pigment content.
Key words: Hyperspectral     Rice     Pigment     Vegetation index     Band depth analysis     Principal component analysis     Back propagation neural network    

色素含量是描述作物生理状态的重要指标之一, 其变化可用于评估作物的光合能力和初级生产力。实时快速地获取色素含量信息成为监测水稻(Oryza sativa)生长、发育和产量的一种有效方法。传统的作物叶片色素含量测定常采用实验室理化分析方法, 费时、繁琐并具有破坏性。近年来, 高光谱技术以其波段连续性强、波谱分辨率高、光谱信息量大等优势为快速、高效、无损检测作物叶片色素含量提供了一种有效途径[1-2]。因此, 研究高光谱技术定量估算作物色素含量具有重要的实践意义。

在可见光范围内, 植被冠层反射光谱主要受叶片色素含量的影响; 而在近红外范围叶片组织结构、生物量、蛋白质、水分等起主要作用。大量研究表明, 可以用植被的冠层反射光谱来估算色素含量[3]。作物叶片中的色素主要包括叶绿素(叶绿素a、叶绿素b)和类胡萝卜素(胡萝卜素和叶黄素), 其中叶绿素是吸收光能的主要物质, 直接影响光合作用的光能利用。研究发现叶绿素a和叶绿素b的吸收峰波段分别为665 nm和643 nm, 类胡萝卜素的吸收峰所处的波段与叶绿素相重叠, 单独估计不同色素含量存在一定的难度。文献报道通过建立遥感原始波段或植被指数与生化组分含量之间的多元统计回归模型, 能较准确地估测植物叶片生化组分含量。王福民等[2]通过所有光谱波段的两两组合, 建立归一化比值色素与叶绿素、类胡萝卜素的统计模型, 得到最佳归一化色素指数。Sims等[4]通过寻找新的光谱指数来预测叶片色素, 但选取的特征光谱及参数因不同作物或不同试验、不同叶片结构而有明显差异。最明显的问题在于, 当叶片色素含量较高时, 常规植被指数将趋于饱和水平, 估算色素含量的精度大大降低[5]。Lichtenthaler等[6]指出NDVI对偏高叶绿素浓度不敏感, 原因在于NDVI对叶绿素吸收带的高度敏感性, 较低的叶绿素含量就能使得叶片对光吸收达到饱和状态。

针对上述问题, Gitelson等[7-8]提出红边位置700 nm处反射率比植被指数对叶绿素的吸收特征更加敏感。但是, 当叶片叶绿素含量改变时, 对应的吸收波段、红边位置与红边面积也随之改变[9]。连续统去除法可以矫正由于波段依赖而引起的波段反射率极值点的漂移, 将波段极值点调整到其真正的波段位置[10]。在此基础上, 对使用连续统去除之后的光谱进行波段深度分析, 增强吸收特征及减少无关的信息冗余[11]。Mutanga等[12]对实验室可控条件下的草地冠层光谱进行波段深度分析, 并与逐步多元线性回归结合估算草地地上生物量。Kokaly等[13]运用波段深度分析方法能够很好地估算干枯落叶的生理生化含量。Chen等[14]将波段深度分析与PLS结合来估算高覆盖草地的地上生物量。

目前在植被冠层光谱定量分析中常采用线性方法, 如逐步多元回归和PLS等。但是, 由于受天气、作物品种、土壤背景、试验条件等众多因素的影响, 冠层光谱与植物叶片组分含量之间的关系是非线性的, 现有方法很难拟合这种多因素影响下的非线性关系。人工神经网络因具有极强的非线性处理、自组织调整、自适应学习的能力而得到广泛应用。本文以提高水稻全生育期的冠层光谱估计类胡萝卜素和叶绿素的精度为主要目的, 采用波段深度分析(BDA)计算波段深度指数, 然后通过主成分分析(PCA)对波段深度指数进行降维, 建立BP神经网络模型预测水稻叶片的色素含量, 并探讨BDA与BP神经网络结合解决植被指数饱和问题的可能性和有效性。

1 材料与方法 1.1 试验材料

试验田位于湖北省武穴市梅川镇(30°06′N, 115°35′E), 土壤类型为水稻土, 0~20 cm耕层土壤pH 5.70, 有机质37.35 g·kg-1, 全氮1.60 g·kg-1, 速效磷5.83 mg·kg-1, 速效钾91.67 mg·kg-1。供试水稻品种为‘深两优5814’, 属杂交籼稻, 全生育期约138 d, 株型适中, 叶片挺直。

1.2 试验设计与方法

试验田内共设置24个小区, 氮肥处理设0 kg·hm-2、45 kg·hm-2、82.5 kg·hm-2、127.5 kg·hm-2、165 kg·hm-2、210 kg·hm-2、247.5 kg·hm-2、292.5 kg·hm-28个水平, 分别用N0、N3、N5.5、N8.5、N11、N14、N16.5和N19.5表示。各处理3次重复, 随机区组排列。每个小区面积20.0 m2。磷、钾肥按P2O575 kg·hm-2、K2O 75 kg·hm-2施入。氮肥为尿素(含N 46%), 磷肥用过磷酸钙(含P2O5 12%), 钾肥为氯化钾(含K2O 60%), 所有肥料做基肥一次性施入。种植密度为255 000株·hm-2, 2015年5月25日播种, 6月28日移栽, 10月9日成熟。

1.3 测量项目与方法 1.3.1 水稻冠层高光谱反射率的测定

分别于2015年水稻分蘖期(7月10日)、拔节期(7月27日)、孕穗期(8月12日)、抽穗期(8月27日)和灌浆期(9月12日), 在相应时期内选择晴朗无云、无风天气, 每个小区选取5个有代表性观测点。在每个观测日的10:00—14:00, 使用FieldSpec Pro FRTM光谱仪(美国Analytical Spectral Devices公司生产)测量冠层光谱反射率。该光谱仪视场角为25°, 波谱范围350~2 500 nm, 其中350~1 000 nm区间光谱分辨率为3 nm, 采样间隔为1.4 nm; 在1 000~2 500 nm区间光谱分辨率为10 nm, 采样间隔为2 nm。各小区在测量前均进行标准白板校正, 测量时仪器探头垂直朝下, 距离水稻冠层高度约1.0 m。各观测点重复观测5次, 取25条光谱的算术平均值作为相应小区冠层光谱结果, 5个物候期24个小区共采集120条有效冠层光谱数据。图 1表示孕穗期不同氮素水平下的水稻冠层光谱反射率变化趋势, 其他关键物候期呈现相同规律。

图 1 不同氮素水平下水稻孕穗期冠层光谱反射率 Figure 1 Canopy spectral reflectance of rice at booting stage under different nitrogen applications
1.3.2 水稻叶片叶绿素及类胡萝卜素的测定

采用混合液(丙酮:无水乙醇=1:1) 提取法测定水稻叶片叶绿素(a+b)和类胡萝卜素含量, 每个小区对应测光谱的冠层, 与冠层光谱测量同步。选取有代表性植株4株, 从采样的水稻植株上等数量选取上中下完全展开叶, 剪碎, 混均匀后获取0.200 g, 加混合液50 mL, 于室温下遮光静置24 h左右至样品完全发白, 用分光光度计测定叶绿素、类胡萝卜素含量[15]

1.4 数据处理与分析 1.4.1 数据预处理

进行色素含量估算时, 仅选用信噪比较高的400~1 350 nm作为有效分析数据。为消除噪声, 采用Savitzky-Golay卷积平滑法进行光谱去噪, 移动窗口宽度为15, 多项式次数为2。将测得的色素含量与对应的冠层高光谱组成原始数据集, 随机分成两组, 第1组包含80个用于建模的样本数据, 第2组包含40个样本数据用于模型验证。

1.4.2 波段深度分析(BDA)

作物叶片的色素含量主要与可见光波段紧密相关, 因此深度分析主要处理400~750 nm的光谱数据, 此波段包含了色素强吸收的蓝、红光区及“红边”波段范围。首先, 对光谱数据进行连续统去除变换, 校正由于波段依赖而引起的波段反射率极值点的偏移, 有效增强吸收特征, 增大各光谱曲线之间的差异[10]。“连续统线”定义为连接局部原始光谱反射率峰值点之间的线段(图 2)。将反射率曲线上每个波长的光谱反射率(R)除以相应波长处连续统线上的值(Rc)可以得到连续统去除光谱(R')。

图 2 400~750 nm连续统去除变换光谱曲线 Figure 2 Continuum removed reflectance spectra from 400 to 750 nm

然后, 对连续统去除变换光谱按照表 1中的公式计算得到波段深度(band depth, BD)、波段深度比(band depth ratio, BDR)、归一化波段深度指数(normalized band depth index, NBDI)、归一化面积波段深度(band depth normalized to band area, BNA), 结果能够进一步反映色素的吸收特征。

表1 吸收特征光谱变换的波段深度分析公式 Table 1 Band depth analysis formulas of absorption feature transform spectrum
1.4.3 建模方法

主成分分析(PCA)是一种通过降维技术把多个变量化为少数几个主成分的多元统计方法。本文采用PCA分别对原始波段(R)、波段深度(BD)、波段深度比(BDR)、归一化波段深度指数(NBDI)、归一化面积波段深度(BNA)提取主成分, 为了后续构建BP神经网络的一致性, 各取前10个主成分, 对应的累积贡献率分别为99.49%、99.85%、99.55%、90.46%和99.97%。经过PCA处理后, 将光谱维数从351减少到10, 既实现了高光谱数据的降维处理, 又保留了原光谱数据绝大部分的特征信息。然后将主成分代入BP神经网络模型进行训练。

BP神经网络模型是一种由输入层、隐含层和输出层组成的前向网络, 学习过程由信号的正向传播与误差的反向传播两个过程组成。正向传播时, 输入的样本以S型传递函数从输入层经各神经元按层处理, 通过隐含层再传向输出层。如果输出层的实际输出与期望输出不符, 则转入误差反向传播, 即将输出误差以某种形式通过隐含层向输入层反传, 并将误差分摊给各层所有单元, 进行权值修正, 从而使误差信号趋于最小, 满足预先设定的限差要求[16]

以10个主成分作为BP神经网络的输入矢量, 隐含层节点数直接影响网络的非线性预测性能, 采用“试错法”确定隐含层节点数, 即尝试不同值的对比, 获得最佳节点数, 类胡萝卜素和叶绿素对应模型的隐含层节点数分别为10、11。1个输出层为水稻叶片色素含量。网络训练的迭代次数为10 000, 训练目标参数为0.001, 学习误差为0.01。训练函数用L-M优化算法函数Trainlm(图 3)。

图 3 估计水稻色素含量的BP神经网络模型 Figure 3 BP neural network model of estimating pigment content in rice
1.4.4 色素光谱参数计算

参照已有的色素光谱参数计算方法, 选择常用的、认可度较高且与本研究相关的17个色素光谱指数(表 2), 分别进行了线性和非线性模型拟合, 挑选色素最佳估计模型。

表2 常用估算水稻叶片色素含量的高光谱参数 Table 2 Hyperspectral indices for estimating leaf pigment content of rice

采用拟合决定系数(R2)、拟合均方根误差(RMSEC)评价模型的建模精度; 估计决定系数(R2)、估计均方根误差(RMSEP)和相对分析误差(RPD)评价模型的验证精度。R2反映模型建立和验证的稳定性, R2越接近于1, 说明模型的稳定性越好、拟合程度越高。RMSE和RPD用来评价模型的误差和估计精度, RMSE值越小说明模型估算能力越好, RPD值大于1.4时说明模型可用, 大于2.0时模型非常优异[17]。数据处理在Matlab 2010b中编程实现。

均方根误差(root mean square error):

${\text{RMSE}} = \sqrt {\frac{{\sum\nolimits_{i = 1}^n {{{({y_i} - {{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{y} }_t})}^2}} }}{n}} $ (1)

相对分析误差(relative percent deviation):

${\text{RPD}} = \frac{{\sqrt {\frac{{\sum {{{({y_i} - \bar y)}^2}} }}{n}} }}{{{\text{RMSE}}}}$ (2)

式中: ${{y_i}}$代表实测值, ${{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{y} }_t}}$代表预测值, ${\bar y}$代表实测值的均值, n为样本数。

2 结果与分析 2.1 类胡萝卜素、叶绿素的含量统计

建模集类胡萝卜素含量均值和标准差分别为0.72 mg·g-1和0.24 mg·g-1; 验证集类胡萝卜素含量均值和标准差分别为0.71 mg·g-1和0.19 mg·g-1 (表 3)。建模集叶绿素含量均值和标准差分别为2.16 mg·g-1和0.79 mg·g-1; 验证集叶绿素含量均值和标准差分别为2.21 mg·g-1和0.66 mg·g-1(表 3)。试验所划分的建模集和验证集样本均值差异较小, 标准差较低, 两者具有很好的相似性, 说明样本的划分满足模型的构建和检验需要。

表3 水稻样本叶片色素含量的统计特征 Table 3 Statistical characteristics of pigment content of leaves of rice samples

对各生育期不同施肥水平下的水稻叶片类胡萝卜素和叶绿素含量进行单因素方差分析, 结果表明水稻分蘖期到拔节期, 不同施肥水平处理下, 叶片类胡萝卜素含量差异不显著; 孕穗后叶片类胡萝卜素开始有显著性差异。水稻各生育期, 不同施氮处理下的叶片叶绿素含量差异显著。同时类胡萝卜素和叶绿素含量的总体变化趋势一致, 即分蘖期到孕穗期含量逐渐增加, 到抽穗期有明显减少的过程, 灌浆期的色素含量明显低于其他各期(图 4)。水稻分蘖期到孕穗期, 叶片中的类胡萝卜素和叶绿素含量逐渐增加, 孕穗前水稻植株处于营养生长阶段, 叶片还在继续发育, 类胡萝卜素和叶绿素含量不断增加。水稻孕穗期到灌浆期, 叶片中的类胡萝卜素和叶绿素含量逐渐减少, 反映出水稻抽穗初期, 叶片内的营养物质向穗部转移, 类胡萝卜素和叶绿素的合成减弱。而在灌浆时期, 叶片慢慢衰老变黄, 叶片作为光合器官的作用进一步削弱, 合成的类胡萝卜素和叶绿素越来减少。

图 4 不同施氮水平下水稻各生育期叶片类胡萝卜素(a)和叶绿素(b)含量的变化(n=120) Figure 4 Carotenoid (a) and chlorophyll (b) contents of rice leaves at different growth stages at different nitrogen applications (n=120)
2.2 各光谱指数反演模型的比较分析 2.2.1 光谱指数与色素含量相关性分析

分析色素含量与光谱参数的相关性发现, 在全生育期, 所有光谱参数均与叶片类胡萝卜素和叶绿素含量之间存在显著相关关系, 一些与红边相关的光谱参数, 如CTR2、PSSNRb等均与叶片类胡萝卜素和叶绿素含量达到较好的相关性(0.7以上)。CTR2与类胡萝卜素呈现显著的负相关, 相关系数达-0.77。PRI与叶绿素的含量相关性最大, 相关系数达0.8以上(表 4)。叶片光化学指数PRI的定义为531 nm和570 nm处反射率的归一化植被指数, 这两个波段位置的反射率受到叶黄素循环的影响并和叶片的光能利用率密切相关[28]。Filella等[29]指出叶片叶绿素含量和叶黄素含量之间有显著的相关性, 叶绿素含量可以间接反映叶黄素的变化。而且植被光能利用率与叶绿素进行光合作用的能力息息相关, 因而能够很好地估算叶绿素含量。

表4 高光谱特征参数与水稻叶片色素含量之间的相关系数分析(n=120) Table 4 Correlation analysis between pigment centent and hyperspectral indices of rice leaves (n=120)
2.2.2 光谱指数对色素含量估计能力

应用17种有代表性的光谱参数, 在以拟合R2最大的优选原则基础上, 构造对类胡萝卜素和叶绿素的最佳预测模型, 用来与BDA和BP神经网络结合建立的水稻叶片色素含量估算模型进行精度比较(表 5, 表 6)。结果表明类胡萝卜素含量的光谱指数模型拟合度较好, R2均在0.6附近, RMSE均相对较小, 而模型验证结果可行度差(R2 < 0.50, RPD < 1.4)。叶绿素光谱指数模型中, PRI、PSRI和SIPI拟合度和验证精度均较高(R2 > 0.60, RPD > 1.4), 可以用来估算叶绿素含量。其余光谱指数建模结果不理想, 模型的验证集精度不高, RMSEP较大, RPD较小。上述结果说明色素光谱指数模型不能很好地估算水稻叶片类胡萝卜素含量, 而植被指数如PRI、PSRI和SIPI使用了与叶绿素作用相关的蓝红光吸收和绿峰相对反射的波段, 可以较好地估算叶片叶绿素含量。

表5 水稻叶片类胡萝卜素含量(Cars)与高光谱参数的线性与非线性回归模型 Table 5 Comparison of linear vs non-linear spectral index regression models of rice leaf carotenoid content (Cars)
表6 水稻叶片总叶绿素含量(Chl)与高光谱参数的线性与非线性回归模型 Table 6 Comparison of linear vs non-linear spectral index regression models of rice leaves chlorophyll content (Chl)
2.2.3 各光谱指数对不同色素含量样本的饱和现象分析

图 5显示了类胡萝卜含量与估算精度较高的CTR2、PSSNRa、PSSNRb的散点关系以及叶绿素含量与估算精度高的PRI、PSRI、SIPI的散点关系。当类胡萝卜素含量大于0.8 mg·g-1时, CTR2、PSSNRa、PSSNRb与类胡萝卜素含量的散点分布呈现随机性; 当叶绿素含量大于2 mg·g-1时, PRI、PSRI、SIPI与叶绿素含量的散点分布也呈现不规律变化, 即色素含量增加到一定程度, 高光谱参数不能有规律地反映叶绿素含量的变化, 表现出上下浮动的平稳性。通过散点图确定变化的拐点, 将类胡萝卜素含量分为2个子集(子集Ⅰ:类胡萝卜素含量小于0.8 mg·g-1, 为中低含量样本集; 子集Ⅱ:类胡萝卜素含量大于0.8 mg·g-1, 为高含量样本集), 叶绿素含量分为2个子集(子集Ⅰ:叶绿素含量小于2 mg·g-1, 为中低含量样本集; 子集Ⅱ:叶绿素含量大于2 mg·g-1, 为高含量样本集)。然后在子集中随机选取40个样本, 分别进行类胡萝卜素含量、叶绿素含量的预测, 并以预测结果的RMSE作为评价指标。若光谱指数对色素含量不同的样本具有良好的预测能力, 则各子集的RMSE应趋于一致。

图 5 水稻叶片类胡萝卜素含量、叶绿素含量与高光谱参数的散点关系变化图 Figure 5 Relationships between spectral parameters and carotenoid content and chlorophyll content of rice leaves 各参数的意义见表 2 Meanings of parameters are shown in the Table 2.

根据CTR2、PSSNRa、PSSNRb估算类胡萝卜素含量并利用PRI、PSRI、SIPI估算叶绿素含量预测结果的RMSE值可知(表 7), 当色素含量高时, RMSE值越大, 对高含量样本预测精度相对较低, 存在饱和现象。其中与类胡萝卜素相关系数最好的CTR2和与叶绿素相关系数最好的PRI对不同样本含量预测精度均高于其他指数。

表7 光谱指数对水稻叶片不同色素含量样本子集预测结果的均方根误差 Table 7 RMSE of prediction results for samples subsets of different pigment contents of rive leaves
2.3 BDABP结合的色素含量估算

表 8显示BDA与BP结合在建模和验证时均较大程度地提高了水稻叶片色素含量估算精度, 能较好地克服饱和问题。对于类胡萝卜素的预测, BD、BDR、BNA的验证效果所得到的R2、RPD均高于原始波段R的验证效果, RMSE则低于原始波段R验证的对应值(图 6, 图 7), 说明BDA的反演精度要优于原始波段R, 能够更好地挖掘与色素关系密切的信息。BD与BP结合的模型建模集R2最大, 验证集R2为0.61, 虽然比BNA与BP结合的模型验证R2略小, 但有最小的RMSEP和最大的RPD(达1.67), 所以BD与BP结合的模型能够较好地估算水稻叶片类胡萝卜素含量。对于叶绿素的预测, BDA的验证效果所得到的R2、RPD均高于原始波段R的验证效果, RMSE则低于原始波段R验证的对应值, 说明BDA比原始波段R能够更好地预测叶绿素含量。BNA与BP结合的模型建模集R2最大, 验证集R2为0.73, 虽然比BD与BP结合的模型验证R2略小, 差异不大, 但有最小的RMSEP和最大的RPD(达1.91), 所以BNA与BP结合的模型能够较好地估算水稻叶片总叶绿素含量。

表8 BDA与BP结合的模型估算水稻叶片类胡萝卜素和叶绿素含量的效果 Table 8 Performance of combination of BDA and BP for predicting rice leaves carotenoid and chlorophyll contents
图 6 BD与BP结合的模型估算类水稻叶片胡萝卜素含量的建模集(a)和验证集(b) Figure 6 Calibration (a) and validation (b) results for predicting carotenoid contents of rice leaves using the combination of BD and BP
图 7 BNA与BP结合的模型估算水稻叶片叶绿素含量的建模集(a)和验证集(b) Figure 7 Calibration (a) and validation (b) results for predicting chlorophyll contents of rice leaves using the combination of BNA and BP

与基于光谱参数建立的最佳回归模型相比, BDA与BP结合建立的模型均较大程度地提高了水稻叶片色素含量的估算精度。对比类胡萝卜素和叶绿素的估算模型, 叶绿素的估算模型精度均高于类胡萝卜素的模型精度, 由于类胡萝卜素和叶绿素的吸收峰所处波段重叠, 且叶片中叶绿素含量远高于类胡萝卜素含量, 因此类胡萝卜素含量的估算较为困难。

3 结论与讨论

本研究为了解决高光谱植被指数反演水稻叶片色素含量存在的饱和问题, 尝试将波段深度分析与BP神经网络结合建立水稻叶片色素含量估算模型。并与17种光谱参数分别建立的模型进行了精度比较。

研究表明, 进行连续统去除变换后, 更加突出了水稻叶片的吸收特征的差异, 挖掘了更多的潜在信息, 使得不同色素含量光谱曲线的差异性特征得到一定程度上的增强, 尤其是红吸收谷区域的波段深度特征蕴藏了较多与水稻叶片色素含量相关的信息, 能够更好地估算水稻叶片色素含量[6]。主要原因可能是波段深度分析处理可以有效减弱土壤背景、大气散射和吸收对目标光谱特征的影响。Kokaly等[13]指出, 波段深度分析方法可以降低土壤背景和大气吸收的影响, 有效估算出叶片的生化参数。Mutanga等[12]表明, 波段深度分析可以解决现有指数的饱和问题, 能够估算高密度植被区域的生物量。

大量研究结果表明, 由于受天气、日照、叶片几何结构、土壤背景、人为操作等因素的影响, 植被的生理生化参量与冠层光谱的关系是非线性的[30-31]。而神经网络对复杂且非线性问题的拟合有着极大的优势, 能够很好的预测非线性函数逼近等问题, 成为当前解决非线性问题的重要方法。

植被指数在色素含量较高时存在饱和问题, 构建植被指数与水稻叶片色素含量的最佳回归模型估算精度不高。马文勇等[32]通过植被指数反演草地叶绿素, 结果表明光谱指数构造形式多样, 且与草地叶绿素含量关系复杂, 在一定程度上影响叶绿素的估算精度。姜海玲等[33]通过地面高光谱数据重采样模拟不同传感器, 利用光谱指数反演植被叶绿素含量, 指出不同传感器的光谱指数估算精度和稳定性均不同。本文将波段深度分析与BP神经网络结合建立的水稻叶片色素含量模型能较好地解决饱和问题, 针对植被指数所囊括的波段具有一定的局限性, 能够实现不同波段之间的优势互补, 提高估算精度。

通过对可见光波段数据的波段深度分析, 利用主成分分析进行降维, 有效保留了信息, 增强了网络性能, 防止过度拟合, 提高了色素含量反演精度。然而本文只针对单一品种和小区试验条件下, 尚需要通过不同地区和品种类型等田间试验做进一步的验证和完善, 辐射传输模型更具有普适性, 仍需我们开展相关研究, 建立经验-物理耦合模型[34], 进一步提高模型估算的准确性和适用性。同时, 在更大尺度、更广范围内实现所构建模型与空间遥感信息的结合, 可为大区域估算水稻叶片色素含量提供理论依据和实践参考。

参考文献
[1] 唐延林, 黄敬峰, 王人潮. 水稻不同发育时期高光谱与叶绿素和类胡萝卜素的变化规律[J]. 中国水稻科学, 2004, 18(1): 59–66.
Tang Y L, Huang J F, Wang R C. Change law of hyperspectral data with chlorophyll and carotenoid for rice at different developmental stages[J]. Chinese Journal of Rice Science, 2004, 18(1): 59–66.
[2] 王福民, 黄敬峰, 王秀珍. 水稻叶片叶绿素、类胡萝卜素含量估算的归一化色素指数研究[J]. 光谱学与光谱分析, 2009, 29(4): 1064–1068.
Wang F M, Huang J F, Wang X Z. Normalized difference ratio pigment index for estimating chlorophyll and carotenoid contents in leaves of rice[J]. Spectroscopy and Spectral Analysis, 2009, 29(4): 1064–1068.
[3] Blackburn G A. Relationships between spectral reflectance and pigment concentrations in stacks of deciduous broadleaves[J]. Remote Sensing of Environment, 1999, 70(2): 224–237. DOI:10.1016/S0034-4257(99)00048-6
[4] Sims D A, Gamon J A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages[J]. Remote Sensing of Environment, 2002, 81(2/3): 337–354.
[5] Datt B. Remote sensing of chlorophyll a, chlorophyll b, chlorophyll a+b, and total carotenoid content in eucalyptus leaves[J]. Remote Sensing of Environment, 1998, 66(2): 111–121. DOI:10.1016/S0034-4257(98)00046-7
[6] Lichtenthaler H K, Gitelson A, Lang M. Non-destructive determination of chlorophyll content of leaves of a green and an aurea mutant of tobacco by reflectance measurements[J]. Journal of Plant Physiology, 1996, 148(3/4): 483–493.
[7] Gitelson A A, Merzlyak M N, Lichtenthaler H K. Detection of red edge position and chlorophyll content by reflectance measurements near 700 nm[J]. Journal of Plant Physiology, 1996, 148(3/4): 501–508.
[8] Gitelson A A, Merzlyak M N. Remote estimation of chlorophyll content in higher plant leaves[J]. International Journal of Remote Sensing, 1997, 18(12): 2691–2697. DOI:10.1080/014311697217558
[9] 姚付启, 张振华, 杨润亚, 等. 基于红边参数的植被叶绿素含量高光谱估算模型[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.
[10] 黄敬峰, 王福民, 王秀珍. 水稻高光谱遥感实验研究[M]. 杭州: 浙江大学出版社, 2010.
Huang J F, Wang F M, Wang X Z. Hyperspectral Experiment for Paddy Rice Remote Sensing[M]. Hangzhou: Zhejiang University Press, 2010.
[11] Mutanga O, Skidmore A K. Integrating imaging spectroscopy and neural networks to map grass quality in the Kruger National Park, South Africa[J]. Remote Sensing of Environment, 2004, 90(1): 104–115. DOI:10.1016/j.rse.2003.12.004
[12] Mutanga O, Skidmore A K. Hyperspectral band depth analysis for a better estimation of grass biomass (Cenchrus ciliaris) measured under controlled laboratory conditions[J]. International Journal of Applied Earth Observation and Geoinformation, 2004, 5(2): 87–96. DOI:10.1016/j.jag.2004.01.001
[13] Kokaly R F, Clark R N. Spectroscopic determination of leaf biochemistry using band-depth analysis of absorption features and stepwise multiple linear regression[J]. Remote sensing of Environment, 1999, 67(3): 267–287. DOI:10.1016/S0034-4257(98)00084-4
[14] Chen J, Gu S, Shen M G, et al. Estimating aboveground biomass of grassland having a high canopy cover: An exploratory analysis of in situ hyperspectral data[J]. International Journal of Remote Sensing, 2009, 30(24): 6497–6517. DOI:10.1080/01431160902882496
[15] 黄秋婵, 韦友欢, 韦方立, 等. 常见夹竹桃科植物叶片颜色及类胡萝卜素含量的比较分析[J]. 安徽农业科学, 2011, 39(25): 15203–15204.
Huang Q C, Wei Y H, Wei F L, et al. The comparative analysis of leaf color and carotenoids content of the common apocynaceae plants[J]. Journal of Anhui Agricultural Science, 2011, 39(25): 15203–15204. DOI:10.3969/j.issn.0517-6611.2011.25.007
[16] 贾方方, 张黎明, 任天宝, 等. 基于BP神经网络的烟草叶片质体色素高光谱反演[J]. 烟草科技, 2016, 49(7): 8–13.
Jia F F, Zhang L M, Ren T B, et al. Hyperspectral inversion to estimate plastid pigment contents in tobacco leaves based on BP neural network[J]. Tobacco Science & Technology, 2016, 49(7): 8–13.
[17] Wang S Q, Li W D, Li J, et al. Prediction of soil texture using FT-NIR spectroscopy and PXRF spectrometry with data fusion[J]. Soil Science, 2013, 178(11): 626–638. DOI:10.1097/SS.0000000000000026
[18] Gitelson A, Merzlyak M N. Spectral reflectance changes associated with autumn senescence of Aesculus hippocastanum L. leaves. Spectral features and relation to chlorophyll estimation[J]. Journal of Plant Physiology, 1994, 143(3): 286–292. DOI:10.1016/S0176-1617(11)81633-0
[19] 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. DOI:10.1016/S0034-4257(98)00059-5
[20] Carter G A. Ratios of leaf reflectances in narrow wavebands as indicators of plant stress[J]. International Journal of Remote Sensing, 1994, 15(3): 697–703. DOI:10.1080/01431169408954109
[21] Rouse Jr J W, Haas R H, Deering D W. Monitoring the vernal advancement and retrogradation (green wave effect) of natural vegetation[R]. Greenbelt, MD, USA: NASA/GSFC, 1974
[22] 朱西存, 赵庚星, 王瑞燕, 等. 苹果叶片的高光谱特征及其色素含量监测[J]. 中国农业科学, 2010, 43(6): 1189–1197.
Zhu X C, Zhao G X, Wang R Y, et al. Hyperspectral characteristics of apple leaves and their pigment contents monitoring[J]. Chinese Journal of Rice Science, 2010, 43(6): 1189–1197. DOI:10.3864/j.issn.0578-1752.2010.06.011
[23] Daughtry C S T, Walthall C L, Kim M S, et al. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance[J]. Remote Sensing of Environment, 2000, 74(2): 229–239. DOI:10.1016/S0034-4257(00)00113-9
[24] Gamon J A, Peñuelas J, Field C B. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency[J]. Remote Sensing of Environment, 1992, 41(1): 35–44. DOI:10.1016/0034-4257(92)90059-S
[25] Merzlyak M N, Gitelson A A, Chivkunova O B, et al. Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening[J]. Physiologia Plantarum, 1999, 106(1): 135–141. DOI:10.1034/j.1399-3054.1999.106119.x
[26] Gitelson A A, Kaufman Y J, Merzlyak M N. Use of a green channel in remote sensing of global vegetation from EOS-MODIS[J]. Remote Sensing of Environment, 1996, 58(3): 289–298. DOI:10.1016/S0034-4257(96)00072-7
[27] Penuelas J, Baret F, Filella I. Semi-empirical indices to assess carotenoids/chlorophyll a ratio from leaf spectral reflectance[J]. Photosynthetica, 1995, 31(2): 221–230.
[28] 薛惠云, 张永江, 刘连涛, 等. 干旱胁迫与复水对棉花叶片光谱、光合和荧光参数的影响[J]. 中国农业科学, 2013, 46(11): 2386–2393.
Xue H Y, Zhang Y J, Liu L T, et al. Responses of spectral reflectance, photosynthesis and chlorophyll fluorescence in cotton during drought stress and rewatering[J]. Scientia Agricultura Sinica, 2013, 46(11): 2386–2393. DOI:10.3864/j.issn.0578-1752.2013.11.024
[29] Filella I, Porcar-Castell A, Munné-Bosch S, et al. PRI assessment of long-term changes in carotenoids/chlorophyll ratio and short-term changes in de-epoxidation state of the xanthophyll cycle[J]. International Journal of Remote Sensing, 2009, 30(17): 4443–4455. DOI:10.1080/01431160802575661
[30] Curran P J, Dungan J L, Peterson D L. Estimating the foliar biochemical concentration of leaves with reflectance spectrometry: Testing the Kokaly and Clark methodologies[J]. Remote Sensing of Environment, 2001, 76(3): 349–359. DOI:10.1016/S0034-4257(01)00182-1
[31] 杨可明, 孙阳阳, 刘飞, 等. 叶绿素含量BP反演模型的光谱信息输入因子构建研究[J]. 科学技术与工程, 2015, 15(15): 82–87.
Yang K M, Sun Y Y, Liu F, et al. Research on constructing the input factors of BP model for inversing leaf chlorophyll content based on the spectral information[J]. Science Technology and Engineering, 2015, 15(15): 82–87. DOI:10.3969/j.issn.1671-1815.2015.15.015
[32] 马文勇, 王训明. 基于高光谱分析的草地叶绿素含量估算研究进展[J]. 地理科学进展, 2016, 35(1): 25–34.
Ma W Y, Wang X M. Progress on grassland chlorophyll content estimation by hyperspectral analysis[J]. Progress in Geography, 2016, 35(1): 25–34.
[33] 姜海玲, 杨杭, 陈小平, 等. 利用光谱指数反演植被叶绿素含量的精度及稳定性研究[J]. 光谱学与光谱分析, 2015, 35(4): 975–981.
Jiang H L, Yang H, Chen X P, et al. Research on accuracy and stability of inversing vegetation chlorophyll content by spectral index method[J]. Spectroscopy and Spectral Analysis, 2015, 35(4): 975–981.
[34] 徐晋, 蒙继华. 农作物叶绿素含量遥感估算的研究进展与展望[J]. 遥感技术与应用, 2016, 31(1): 74–85.
Xu J, Meng J H. Overview on estimating crop chlorophyll content with remote sensing[J]. Remote Sensing Technology and Application, 2016, 31(1): 74–85.