高光谱成像用于定量判断茶树叶片受冻程度_毛艺霖.pdf

第 卷，第期 光谱学与光谱分析 ，年月 ，高光谱成像用于定量判断茶树叶片受冻程度毛艺霖，李赫，王玉，范凯，孙立涛，王会，宋大鹏，申加枝，丁兆堂，青岛农业大学茶叶研究所，山东 青岛 山东省农业科学院茶叶研究所，山东 济南 日照市农业科学研究院茶叶研究所，山东 日照 摘要低温冻害是茶园中最常见的自然灾害之一。茶树叶片低温胁迫的定量监测对于评估茶园冻害程度和及时采取措施具有重要意义。茶树低温胁迫的传统检测方法，主要是通过人工观察和理化指标的测定，存在精度低、效率低和主观性强等问题，严重影响了灾害后期的茶树管理。该研究提出了一种基于高光谱成像的茶树冻害程度定量判断方法。首先，利用高光谱成像设备采集自然环境中茶树叶片在无冻害发生、冻害初期和冻害后期三个阶段的光谱数据，提取叶片的平均反射率；测定相应叶片中的相对电导率（）、叶绿素（）和丙二醛（）等生理生化指标。其次，利用多元散射校正（）、一阶导数（）和平滑滤波（）算法对采集的原始高光谱数据进行预处理，并利用无信息变量消除（）和竞争性自适应重加权（）算法筛选预处理后高光谱数据的特征波段。最后，利用卷积神经网络（）、支持向量机（）和偏最小二乘法（）建立 、和 含量的定量预测模型。结果表明：（）经 算法预处理的光谱曲线比原始光谱曲线的波峰和波谷更加突出，提高了光谱的分辨率和灵敏度，有利于提高后期回归模型的精度；（）算法筛选的特征波段数量最多，后期建模效果更优；算法筛选的特征波段数量最少，更适合用传统的机器学习方法建立回归模型；（）、和 的最佳预测模型分别为 （，）、（，）和 （，）。利用高光谱成像技术与多种算法相结合，可以实现对茶树叶片低温胁迫程度的无损、准确和定量监测，对快速预测茶园冻害发生和采取必要措施具有重要意义。关键词茶树；冻害；高光谱成像；深度学习中图分类号：文献标识码：（）收稿日期：，修订日期：基金项目：国家自然科学基金项目（），山东省农业科技创新项目（，），山东省农科院创新工程项目（，）资助作者简介：毛艺霖，女，年生，青岛农业大学茶叶研究所硕士研究生 ：通讯作者 ：；引言茶树（）是重要的经济作物，栽培历史悠久，分布广泛。茶树的起源中心位于气候温暖湿润的中国西南部，最适宜生长温度为 。低温已经成为影响茶树越冬和春茶生长的重要环境因素之一，严重影响茶叶生产的经济效益。因此，如何快速预测冻害的发生并采取必要措施，是茶叶生产中减轻冻害影响的重要问题。然而，传统的低温胁迫检测主要是通过茶树生理生化指标的测定和灾后实地表观的调查，存在耗时费力、准确性低和预防滞后等问题。因此，迫切需要提出一种无损、快速和准确的低温诱导成分监测方法。高光谱成像（）是一种新兴的光学技术，已成为取代传统破坏性取样和实验室测试的经济型方法。目前，高光谱技术已经应用于检测小麦、甜菜等作物的内含成分。此外，该技术也应用于茶叶的理化成分分析。例如，等利用高光谱成像技术监测红茶加工过程中茶多酚、游离氨基酸和咖啡碱的含量，从而定量判断萎凋和发酵的程度；等利用高光谱监测茶树叶片中干旱诱导成分，从而快速评价茶树的干旱程度。这些研究表明了高光谱成像技术在估算茶树内含成分方面的潜力。然而，基于高光谱成像监测茶树低温诱导成分和冻害状况的研究仍然缺乏。在高光谱数据的实际应用中，通常将其与波段选择算法和机器学习算法相结合，以提高模型的运算速度和鲁棒性。高光谱图像中存在着与研究目标无关的冗余信息，因此确定哪些信息对研究目标贡献更大尤为重要。因此，许多研究人员将无 信 息 变 量 消 除（）和 竞 争 性 自 适 应 重 加 权（）等波段选择方法应用于光谱数据，以降低维度并去除冗余光谱。目前，机器学习方法已经广泛应用于将高通量数据与作物性状相关联。然而，与传统的机器学习方法相比，深度学习方法在分类和回归问题上具有更好的效果。例如，等提出了一种基于 、小波变换和四通道残差网络（）的茶树叶病虫害识别框架，结果表明 模型的精度要优于支持向量机（）模型。等利用深度神经网络（）、偏最小二乘法（）和 网络建立了玉米叶面积预测模型，结果表明 模型的精度优于其他传统的机器学习方法。目前，深度学习模型在茶树低温胁迫表型方面的预测结果尚不清楚。基于上述分析，提出了一种基于高光谱成像定量判断茶树叶片受冻程度的方法。首先，利用高光谱成像设备采集不同 冻 害 程 度 茶 树 叶 片 的 光 谱 数 据，利 用 多 元 散 射 校 正（）、一阶导数（）和平滑滤波（）算法对高光谱数据进行预处理，通过 和 算法对预处理后光谱数据进行特征波段的筛选。然后，测定叶片中的相对电导率（）、叶绿素（）和丙二醛（）含量。最后，利用卷积神经网络（）的深度学习方法、和 两种机器学习方法建立茶树叶片理化参数的预测模型。实验部分 试验设计试验地点位于中国山东省日照市茶叶研究所（北纬 ，东经 ，海拔 ）。试验茶园分为 个小区（），每个小区取 片成熟叶片。根据园内气象站的测量，该地块 年 月至 年月的空气温度（地上处）和土壤温度（地下 处）变化如图所示。分别于 年 月 无冻害发生期（）、年月 冻害发生初期（）和 年月 冻害发生后期（）采集样品。共采集 个样本。图试验期间研究区域的气温和地温变化 数据采集 高光谱数据采集高光谱数据的采集过程如图所示。高光谱成像采集系统包括台成像光谱相机（，无锡双利合谱科技有限公司）、台对称分布的卤素灯线光源（，中国）、计算机和其他部件。采集系统的外部由黑色暗箱封闭。采集的高光谱图像的空间分辨率为 （空间光谱），光谱通道数为 ，光谱范围为 。为了减少暗电流噪声和外部环境的影响，采集的光谱数据需要进行黑白校正，如式（）。（）式（）中，是校正图像，是原始图像，是通过用反射率约为的镜头盖覆盖而获得的黑色参考图像，是使用反射率约为 的纯白色标准白板获得的白色参考图像。图高光谱数据的采集过程 理化指标的测定测 定 了 叶 片 的 三 个 生 理 生 化 指 标：、和。其中，叶 片 值是由便 携 式 叶 绿 素 仪（，日本）测定。叶片 是由电导率仪（，中国）测定，如式（）。叶片 含量是根据苏州格锐思生物科技有限公司的试剂盒测定，如式（）。含量试剂盒的方法为微板法，产品编号为 。（）（）（）式（）中，为叶片煮沸前外渗液的电导值，为叶片煮沸后外渗液的电导值。式（）中，和 分别是样液在 和 处读取的吸光光度值，是样本质量。光谱数据的标准化和特征提取在 （，中国）软件中打开黑白校正后的高光谱图像，并使用镜头、反射率校准等工具进一步校正图像，以获得标准化的高光谱图像。在 （，美国）软件中打开标准化后的高第期 毛艺霖等：高光谱成像用于定量判断茶树叶片受冻程度光谱图像，将整个茶树叶片区域定义为感兴趣区域（），并提取所有样本的 。然后计算每个样本 的平均反射率。数据分析 光谱数据的预处理由于高光谱设备和环境因素的影响，茶树叶片光谱受散射效应和噪声的影响，削弱了茶叶内部理化指标的光谱信号，影响回归模型的精度。因此，我们利用 、和对高光谱数据进行预处理。特征波段的筛选高光谱数据虽然有助于我们更加细致的对物品分类和识别，但是随着波段的增多，必然会导致信息的冗余和数据处理的复杂性。因此，我们利用 和 算法来选择全波段光谱数据中的代表性波段作为“特征波段”。模型的建立与评价为找到茶树冻害评估的最优预测模型，我们利用 、和 建立了茶树叶片光谱数据与其理化参数之间的回归模型，并比较了模型的效果。为进一步保证算法的准确性，采用十折交叉验证，即将数据集分成 份，轮流将其中份作为训练数据，份作为测试数据，重复三次，取结果的平均值。为了更精确的评价回归模型 的性能，我 们采用 决 定 系 数（）和 均 方 根 误 差（）来评估模型的性能。越大，越小，表明模型的性能越好。和 的公式如式（）和式（）（?）（?）（）（?）（）其中，为对应数据集中的样本个数，?和分别为茶叶样本的预测值和测量值，?为样本的平均实测值。结果与讨论 理化指标的分析 、和 的含量可以作为反映茶树冻害状况的良好指标。它们的测定结果如图所示。结果表明，随着冻害的发生，值逐渐降低，而 和 趋于增加。有趣的是，在 至 时期，每个指标变化相对缓慢，但在 时期变化显著（）。这可能是因为茶树在低温伤害初期具有抵御逆境的自我保护功能。然而，当低温胁迫超过茶树的抵御能力时，叶绿素会大幅度降解、电解液图三个时期的理化指标数据分布（）：值；（）：；（）：含量 （）：；（）：；（）：图原始光谱与预处理后的光谱 光谱学与光谱分析第 卷泄漏率和渗透调节物质会显著增加，导致茶叶受到低温伤害。因此，、和 含量的变化可以用来量化茶树冻害状况。即 值越低，和 越高，则表明茶树冻害越严重。光谱数据的预处理为了去除原始光谱数据的基线漂移、噪声等信息，建立稳定、可靠的定量分析模型，我们结合、和对光谱数据进行预处理。原始平均反射率光谱图和预处理后光谱曲线如图所示。预处理后，可以清晰地观察到光谱的吸收峰和反射谷更加突出，提高了光谱的灵敏度。特征波段的选择和分析为了消除无关波段对模型精度的影响，我们使用 和 算法选择特征波段，如表和图所示。结果表明，在 的筛选算法中，算法筛选出的特征波段最多，为 ；算法筛选的特征波段最少，为。在 的筛选算法中，算法筛选的特征波段最多，为；算法筛选的特 征 波 段 最 少，为 。在 的 筛 选 算 法 中，算法筛选的特征波段最多，为 ；算法筛选的特征波段最少，为。表特征波段筛选的结果 指标筛选算法特征波段的数量 图 和 算法筛选的特征波段（）：；（）：；（）：；（）：；（）：；（）：（）：；（）：；（）：；（）：；（）：；（）：不同模型的比较分析我们通过模型的和 衡量预测模型的性能（表和图），筛选最佳模型。结果表明，在 的预测模型中，模 型 的 精 度 最 高（，），而 的精度最低（，）；在 的预测模型中，模型的精度最高（，），模型的精度最低（，）；在 的 预 测 模 型 中，模型的精度最高（，），模型的精度最低（，）。类似的，基于高光谱的 模型在玉米幼苗的冷害诊断中也取得了较好的效果。总的来说，算法和 算法相结合的模型具有较高的精度，而与 算法相结合的模型精度最低。这可能是因为深度学习可以直接从数据中获取高级特征，其性能将随着数据规 模 的 增 加 而 不 断增长。因而 包 含 更 多 信 息 的 模型的预测精度更高。算法筛选到的信息量较少，并且可能缺少一些关键信息，导致模型精度降低。另一方面，深度学习非常适合大量的数据集，本工作只有 个样本，这在一定程度上限制了深度学习的实际应用。结论提出一种基于高光谱成像茶树受冻程度定量判断的方法。首先，采集冻害未发生、初期和后期的茶树叶片样本，利用高光谱成像获取叶片的光谱数据，并测定相应叶片的 值、和 含量。其次，利用 、和算法对采集的原始高光谱数据进行预处理，并通过 和第期 毛艺霖等：高光谱成像用于定量判断茶树叶片受冻程度表不同模型的评价结果 指标筛选算法建模算法训练集测试集 算法对预处理后高光谱数据进行特征波段的筛选。最后，利 用、和 算 法 建 立 茶 树 叶 片 、和 含量的预测模型。结果表明：（）算法预处理后的光谱曲线比原始光谱曲线波峰和波谷更加突出，提高了光谱的分辨率和灵敏度，有利于提高后期建立回归模型的精度；（）算法筛选的特征波段数量最多，后期建模效果更优，而 算法筛选的特征波段数量最少，更适合与传统机器学习的方法建立回归模型；（）模 型 在 预 测 叶 片 值（，）和 含量（，）方面具有最高的精度，模型在预测叶片 （，图实际值和预测值的散点图（）：；（）：；（）：（）：；（）：；（）：）方面具有最高的精度。高光谱成像技术和多种算法相结合建立的模型具有准确的预测结果，可用于测定 、和 含量。本研究为茶树叶片低温胁迫程度的快速、准确和定量地监测奠定了基础，对快速预测茶园冻害发生并采取必要措施具有重要意义。，（）：，：，（）：，（）：，（）：，（）：，：，：，（）：，（）：，：，（任东，沈俊，任顺，等）（光谱学与光谱分析），（）：，：光谱学与光谱分析第 卷 ，：，（李 彦 东，郝 宗 波，雷航）（计