多光谱图像的玉米叶片含水量检测

彭要奇，肖颖欣，傅泽田，董玉红，李鑫星，严海军，郑永军

1. 中国农业大学工学院食品质量与安全北京实验室，北京 100083 2. 中国农业大学土地科学与技术学院，北京 100083 3. 中国农业大学信息与电气工程学院，北京 100083 4. 中国农业大学水利与土木工程学院，北京 100083 5. 中国农业大学工学院，北京 100083

引 言

玉米是一种需水量较大的粮食作物，及时获取玉米生长过程中植株体的水分信息，对田间灌溉管理和旱情预警监测具有显著的科学意义和现实意义。 绿色植物进行光合反应、呼吸效应和表层蒸腾作用主要是借助叶片来完成，叶片运作过程受外界环境影响明显，作物在缺水状态下，叶片通常会出现蜷缩、枯黄和色斑等形状、颜色和纹理特征的改变，因此检测作物叶片含水量可以作为检测作物整体含水量的一个参考标准[1]。 目前常见的作物植株含水量测定方法主要有烘箱法[2]、涉及高危化学试剂的卡尔·费休法[3]、基于物理层面的电容法[4]和根据波长特征的微波法[5]等。 这些测定方法往往操作流程繁琐、费时，有些甚至需要危险化学试剂，不能满足目前测量工作中制定的安全、低成本和快速化的技术要求。

多光谱图像是指包含多个波段信息的图像，每个波段提取出对应的一幅灰度图像，它包含了该波段传感器的敏感度生成的场景亮度信息[6]。 多光谱图像比传统普通彩色图像包含更多的光谱段，又弥补了高光谱图像拍摄仪器价格昂贵且不易携带的缺点，更加适用于科学研究和实地应用。 多光谱图像在气象探测[7]、地貌绘制[8]、物料分类[9]、农林[10]、海洋[11]、医疗[12]、环境监测[13]等各领域获得了非常广泛的应用。

目前，多光谱图像在农业中的应用有了一定的进展。 陈思明[14]等利用多光谱线性分解方法对原始影像的像元进行分解，并构建了土壤有机质含量的反演模型，实现了湿地土壤有机质含量的准确实时监测。 杨敏[15]等开展基于多光谱图像数据的植被叶面积指数反演方法研究，提高了多光谱图像在植被影响土壤侵蚀方面评估的应用性。 Potgieter[16]等利用多光谱数据提取NDVI指数和EVI指数，分析作物叶片的叶绿素含量，对高粱育种计划的遗传改良进行评估。 Backoulou[17]等使用MS3100-CIR多光谱相机拍摄冬小麦的多光谱图像，以确定绿蝇虫害的分布领域，指导建立有害生物管理和监测系统。 多光谱图像在农业中的应用与常规手段相比，具有部署简单、快速高效，成本低廉和易于推广的优点。 本工作基于多光谱图像，应用改进的神经网络方法，建立玉米叶片含水量反演模型，实现玉米叶片含水量的快速检测。

1 实验部分

1.1 材料

实验材料为按照随机取样的方法采集于拔节期、大喇叭口期和吐丝期等三个关键需水生育期的玉米叶片。 采集地点为中国农业大学通州实验站5号试验田，采集时间为2018年7月19号、8月4号和9月13号三次进行，共55组叶片。

1.2 仪器设备

实验用Micasense公司生产的RedEdge-M型号多光谱相机，见图1。 该型号多光谱相机体积小，重量轻，可以手持操作也可搭载无人机使用，可以同时拍摄包含蓝(475 nm)、绿(560 nm)、红(668 nm)、近红外(840 nm)和红边(717 nm)等5个离散波段的图像。 每一个波段通道都有完整的成像仪窄带光学过滤器，可以获取高精度的多光谱图像。 拍摄前需要采集标准白板信息，对获取的多光谱图像进行校正，主要校正拍摄图像的反射率和辐射亮度信息，拍摄过程中通过直连下行光传感器能够有效消除太阳高度角对光谱反射的影响。

叶片含水量的常用指标是叶片等效水厚度，这一指标计算的原理如式(1)所示[18]。

(1)

式中，EWT为样本等效水厚度；m1为测量区域叶片的鲜重；m2为测量区域叶片的干重；s为测量区域的表面积。

采用YLS-D型号植物营养测定仪进行玉米叶片水分的测量，正常工作的量程在0～99.9RH%之间，工作误差为±0.5个标准单位，见图2。 能够满足测量叶片水分的要求。 使用前需要进行初始化校准。

图1 RedEdge-M多光谱相机

图2 YLS-D型号植物营养测定仪

1.3 样本叶片采集及数据分析

首先在实验站通过随机取样法采集55组玉米叶片，在光线良好的环境下，手持RedEdge-M多光谱相机垂直正对样本叶片进行拍摄，拍摄前需要先拍摄标准白板图像作为校正图像。 拍摄过程中使用WiFi模式直连手机端进行样本图像位置的实时预览，相机镜头距离叶片高度大约25 cm，保证样本位置都处于五个通道传感器的正中间时，按下拍摄按钮进行图像保存，图像格式为16位的TIFF图像。 继续使用初始化校准后的YLS-D型号植物营养测定仪通过五点取样法夹持玉米叶片，记录五个位置的叶片水厚度信息。

使用ENVI5.3软件绘制玉米叶片感兴趣区域(region of interest, ROI)，绘制过程中考虑到传感器镜头的边缘减光现象，所选区域应尽量远离图像边缘且尽可能大的包络玉米叶片，见图3。 以此借助软件可以得到ROI区域内的平均反射光谱，并根据标准白板的标定反射率、白板ROI范围内的平均反射光谱和玉米叶片样本白板ROI范围内的平均反射光谱，比值换算得到玉米叶片5个波段处的光谱反射率。

BP神经网络是一种根据研究样本误差，使用逆传播算法反复训练完善的多层前馈网络, 在数据预测、样本检测和图像融合等多领域具有广泛的实施操作性[19]。 而早期经典的BP神经网络算法在训练过程存在收敛速度较慢、受到个别样本异常值扰动的情况下易出现局部过拟合和网络结构最优化难以确定等缺点； 引入Levenberg-Marquardt方法进行改进，能够进一步减小振荡、加快网络的收敛速度。

图3 玉米叶片ROI的绘制

2 结果与讨论

2.1 玉米叶片含水量反演模型的建立

基于多光谱图像的玉米叶片含水量反演模型的输入神经元数目为5个，即蓝、绿、红、近红外和红边等5个波段对应的反射率； 输出神经元为1个，即玉米叶片等效水厚度。 随机选取43组玉米叶片样本的光谱反射率作为建模数据样本，部分用于建模的训练数据列于表1。

表1 部分建模数据

由于训练的各样本数据均有最大值和最小值，并且单位不统一，采用MATLAB程序对输入数据和输出数据进行归一化处理至同一量纲，经处理后，建模训练样本值归一化范围在(0.1，0.9)之间。

目前对于隐含层单元数的确定尚没有明确最佳的方法，主要依据Kolmogorov定理、经验公式和试凑法来确定[20]。 根据蓝、绿、红、近红外和红边等5个波段图像对应的反射率作为5个输入神经元，玉米叶片含水量作为1个输出神经元，初求解出隐含层神经元数量在2～11个之间，将每个隐含层神经元数目对应的网络结构反复训练20次，记录拟合精度最高的一次，以此确定出最理想的隐含层神经元数目。 隐含层神经元数目训练求解过程的主要参数记录如表2所示。

根据表2训练结果可知，当隐含层神经元数目在9个时，反演模型训练过程中迭代次数最少和均方差最小。 因此确定中BP神经网络中隐含层单元数的数目为9个，并确定基于多光谱图像的玉米叶片含水量反演模型的BP神经网络结构为5-9-1的结构，如图4所示。

表2 隐含层神经元求解主要参数

图4 含水量反演模型神经网络结构

2.2 模型准确度分析

确定完BP神经网络结构后，将43组玉米叶片样本返回至5-9-1的网络结构中继续大量训练，参考网络输出值与目标参考值的相关性达到最佳效果后停止。 基于莱文贝格-马夸特方法改进的玉米叶片含水量反演模型的训练性能和训练拟合相关度信息如图5所示。

从图5可以看出，基于莱文贝格-马夸特方法改进的玉米叶片含水量反演模型训练迭代次数少，收敛速度快，小误差样本较为集中，均方差MSE值是0.009 689 9，整体网络结构训练效果好。 另外对神经网络模型进行相关回归性分析可以看出，基于多光谱图像的玉米叶片含水量反演模型的网络输出值与目标参考值的相关系数达到了0.896 37，玉米叶片含水量反演模型整体具有较高的精度。

2.3 检测结果验证

最后对剩余的12组玉米叶片样本进行反演值和参考值的仿真验证。 将12组玉米叶片样本的5个波段图像的反射率代入反演模型计算出玉米叶片含水量反演值，与参考值的对比结果如表3所示。

根据表3中的反演值和参考值信息绘制仿真结果相关关系图，如图6所示，反演值和参考值具有较高的相关性，二者相关系数r能够达到0.894 8，可见基于多光谱图像的玉米叶片含水量反演模型较为准确可靠，能够满足农业环境参数的检测需求。

图5 反演模型训练信息

表3 玉米叶片含水量反演值和参考值结果

图6 反演值和参考值验证

3 结 论

以55组不同生育期的玉米叶片为研究对象，获取玉米叶片的光谱信息和水分状况，结合莱文贝格-马夸特方法改进神经网络模型构建基于多光谱图像的玉米叶片含水量反演模型，实现了对玉米叶片含水量的快速检测。

首先使用ENVI5.3软件作玉米叶片感兴趣区域的构建，得到ROI区域内的平均反射光谱，参照标准白板出厂时提供的专属标定反射率、白板ROI范围内的平均反射光谱和玉米叶片样本白板ROI范围内的平均反射光谱，比值换算得到玉米叶片5个波段的光谱反射率。 随机选择43组样本作为训练集进行反演模型构建，其余12个样本作为测试集用于仿真验证。 结果表明，经莱文贝格-马夸特方法改进后的玉米叶片含水量反演模型，网络输出值与目标参考值的相关系数达到了0.896 37，反演模型整体具有较高的精度。 对剩余的12组玉米叶片样本导入反演模型计算并进行仿真验证，反演值和参考值具有较高的一致性，二者相关系数r能够达到0.894 8。 因此将多光谱图像用于玉米叶片含水量的检测，可以实现对玉米叶片水分信息的快速获取，为高效的农场管理模式提供了一种全新的解决方案和参考依据。

多光谱对于玉米叶片含水量的检测具有较好效果，对于反演模型的构建，本研究仅仅选取了改进的神经网络进行学习训练，今后可以尝试多种方法包含有多元线性回归、偏最小二乘法、逐步线性回归、主成分分析和支持向量机的方法进行建模，比较各种方法的优劣，选出最适用于作物叶片含水量反演模型构建的方法。