一种小麦叶片衰老程度快速测量方法与流程

本发明涉及植物叶片衰老程度测量技术领域，特别涉及一种小麦叶片衰老程度快速测量方法。

背景技术：

衰老是植物生长发育的最后一个阶段，是普遍存在于生物界的生理现象。叶片是植物光合作用的重要场所，是衰老比较敏感的部位之一。小麦叶片衰老一般是指叶片代谢活动减弱、生理机能衰退的过程，是细胞死亡的一种类型。小麦叶片衰老最显著的变化就是叶片变黄或失绿，衰老的过程就是将可循环利用的营养物质转移到幼嫩叶片或种子中。盐胁迫、干旱胁迫、营养物质的缺乏和病虫害胁迫会引起小麦叶片早衰或加剧叶片衰老，导致严重减产，延缓衰老可以一定程度上提高产量。因此，小麦叶片衰老程度是一项重要指标。

长期以来，小麦叶片衰老程度的测量都是采用间接测量法，即通过测量与叶片衰老表现显著的生理指标来反映，如叶绿素含量、蛋白质含量、光合磷酸化能力、光合速率、呼吸速率以及各种活性酶和内源激素等。部分指标测定方法如下：

净光合速率：于每次采样当天9:00—12:00,用li-6400便携式光合仪测定待采叶片的净光合速率。

叶绿素含量：将叶样去除叶脉后剪成碎片,称取0.2g采用丙酮浸提比色法测定叶绿素含量。

总氮含量：将叶样于烘箱中105℃烘30min,80°c烘至恒重。采用凯氏定氮法测定叶片总氮含量。

气体交换参数：用li-6400便携式光合系统测定。

最大光化学效率（fv/fm）：用image-pam调制荧光仪测定光系统ⅱ（psⅱ）最大光化学效率。

可溶性糖和淀粉含量：用硫酸蒽酮比色法测定可溶性糖和淀粉含量。

由此可见，要研究小麦叶片的衰老问题，只能通过测定与衰老相关的指标，其指标数量多、测定过程繁琐，耗时费力，且无统一衰老程度评判标准。

技术实现要素：

本发明针对现有技术小麦叶征衰老过程评定方法的不便性，提供一种小麦叶片衰老程度快速测量方法，实现小麦叶片衰老程度的快速测量。

本发明的目的是这样实现的，一种小麦叶片衰老程度快速测量方法，包括如下步骤：

a.预备实验阶段：确定小麦衰老程度与叶片颜色变化的关系，包括以下分步：

a.1在小麦叶片衰老全过程中，使用常规方法测定与叶片衰老过程相关的反映光合能力的各项生理指标；

a.2结合a.1步测定的麦叶片衰老全过程的各生理指标及小麦叶片衰老过程，将小麦叶片衰老过程从开始衰老到完全衰老分为若干个等级；

a.3取各等级的叶片平铺在白色背景纸上获取叶片的rgb图像；

a.4利用matlab提取rgb图像中叶片图像的颜色特征参数的r、g、b值；

a.5根据小麦叶片颜色从开始衰老到最终死亡经历由绿转黄到枯黄最后变褐色，按下式确定反应小麦叶片颜色衰老变化的特征参数mt的变化规律；

a.6根据颜色特征参数mt与衰老等级的关系，构建颜色特征参数mt与叶片衰老程度等级的标准回归曲线方程；

b.待测叶片衰老程度测定流程包括以下分步：

b.1将待测小麦叶片平铺在一张白色背景底面上，拍摄小麦叶片的rgb图像；

b.2利用opencv库函数的imread读取叶片图像；

b.3通过灰度化cvtcolor和二值化threshold处理提取图片中的叶片的原rgb图像的r、g、b值；

b.4根据公式（1）计算叶片部位颜色特征值mt均值；

b.5将mt均值代入步骤a.6的标准回归曲线方程中，得到叶片的衰老等级。

本发明的上述小麦叶片衰老程度快速测量方法，在预备实验阶段，通过常规方法，监测小麦叶片衰老过程中各项生理指标，根据各生理指标将小麦叶片衰老过程划分为若干个等级，再结合不同衰老等级叶片颜色的图像数据，确定与衰老程度相关的颜色特征参数mt，再建立小麦衰老等级与颜色特征参数mt的标准回归曲线方程，通过前期预备实验阶段，总结出叶片衰老与颜色特征参数的曲线关系，以便于实际测量中快速应用；在实际叶片衰老程度检测和判定时，通过拍摄获取叶片图像，识别提取图像中叶片部位的r、g、b值，计算颜色特征参数mt后，代入标准回归曲线方程以确定小麦叶片的具体衰老等级。上述小麦叶片衰老程度快速测量方法，通过前期的预备总结得出回归曲线，为实际叶片快速测量做了充分总结，在实际测量中只需拍摄叶片图像，将颜色特征参数代入回归曲线方程中，即可很到叶片衰老等级，相比传统的需要测定多个与衰老程度相关的生理指标，具有省时省力，方便快速的优点，同时确定的测定结果准确率高，不受人为因素的影响。

为便于全面准确测量叶片的生理指标参数，a.1步中，反映光合能力的各项生理指标包括：光合磷酸化能力、净光合速率、蒸腾速率、呼吸速率、气体交换、最大光化学效率fv/fm。

为便于准确分级叶片的衰老程度，a.2步进行小麦叶片衰老分级方法为：根据a.1步测定的小麦衰老过程中各生理指标，利用神经网络adaboost分类算法将小麦叶片衰老程度训练为1～8个等级。

为准确确定mt与叶片颜色的关系，a.5中特征参数mt的变化规律按公式计算得出：

公式(1)；

所述特征参数mt为归一化参数，属于叶片衰老颜色时0<mt<1，不属于叶片衰老颜色时mt=0或1。

为便于快速确定小麦叶片衰老等级，a.6步中的标准回归曲线方程为：y=7.0679x+0.989，其中x值为小麦叶片图像mt均值，y为衰老等级。

为便于识别图像中的叶片并对叶片图像进行数据提取，b.3步中，利用opencv库函数读取叶片的rgb图像，通过灰度化cv_rgb2gra和二值化thresh_binary处理提取图片中的叶片的原rgb图像的r、g、b值，具体为，将三通道彩色图像转换成二通道灰度图像，标准化灰度图像值为0-255，选择阈值为200，将灰度图像转换成只有0或1的二值图像，灰度图像值为0-200的转换成0，200≤灰度图像值＜255的转换成1，转换后的值为0是背景，值为1是叶片，并记录转换后值为1的矩阵位置，再提取值为1叶片位置的原rgb彩色图像的r、g、b值。

为便于叶片图像的拍摄获取，a.3步和b.1步小麦叶片图像通过手机或高清数码相机拍摄。

附图说明

图1为衰老等级为1～8级小麦叶片的图像（已经转化为二通道灰度图像）。

图2为小麦叶片衰老程度标准回归曲线。

具体实施方式

本实施例的小麦叶片衰老程度快速测量方法，包括如下步骤：

a.预备实验阶段：确定小麦衰老程度与叶片颜色变化的关系，包括以下分步：

在小麦叶片衰老全过程中，使用常规方法测定与叶片衰老过程相关的反映光合能力的各项生理指标；这些生理指标包括：光合磷酸化能力、净光合速率、蒸腾速率、呼吸速率、气体交换、最大光化学效率fv/fm。

结合上步测定的麦叶片衰老全过程的各生理指标及小麦叶片衰老过程，将小麦衰老过程中各生理指标参数，利用神经网络adaboost分类算法将小麦叶片衰老程度训练为1～8个从开始衰老到完全衰老的8等级。

上述衰老全过程测量过程中，取各阶段等级的叶片平铺在白色背景纸上获取叶片的rgb图像，如图1所示（已经转化处理为二通道灰度），利用matlab软件提取rgb图像中叶片图像的颜色特征参数的r、g、b值；根据小麦叶片颜色从开始衰老到最终死亡经历由绿转黄到枯黄最后变褐色，总结确定反应小麦叶片颜色衰老变化的特征参数mt的公式：

公式（1）；

该方法计算的特征参数mt为归一化参数，属于叶片衰老颜色时0<mt<1，不属于叶片衰老颜色时mt=0或1；因此，通mt参数既可应用于区分小麦叶片衰老程度又可应用于辨别小麦叶片是否衰老。

根据上述测量过程中得到的颜色特征参数mt与衰老等级的关系，构建颜色特征参数mt与叶片衰老程度等级的标准回归曲线方程；y=7.0679x+0.989，其中x值为小麦叶片图像mt均值，y为衰老等级,所得的小麦叶片衰老程度标准回归曲线如图2所示。

通过上述预处理测量阶段，通过常规方法，监测小麦叶片衰老过程中各项生理指标，根据各生理指标将小麦叶片衰老过程划分为若干个等级，再结合不同衰老等级叶片颜色的图像数据，确定与衰老程度相关的颜色特征参数mt，再建立小麦衰老等级与颜色特征参数mt的标准回归曲线方程，通过前期预备实验阶段，总结出叶片衰老与颜色特征参数的曲线关系，以便于实际测量中通过叶片图像数据快速确定叶片的衰老等级。前述预处理所得的标准回归曲线，可用于同类生长特性的小麦叶片衰老等级的确定中。

通过上述预处理结果进行小麦生长各阶段待测叶片衰老程度测定流程包括以下分步：

首先，将待测小麦叶片平铺在一张白色背景底面上，通过数码相机或手机等高清拍摄装备拍摄小麦叶片的rgb图像；利用opencv库函数的imread读取叶片的rgb图像，将三通道彩色图像转换成二通道灰度图像，标准化灰度图像值为0-255，选择阈值为200，将灰度图像转换成只有0或1的二值图像，灰度图像值为0-200的转换成0，200≤灰度图像值＜255的转换成1，转换后的值为0是背景，值为1是叶片，并记录转换后值为1的矩阵位置，再提取值为1叶片位置的原rgb彩色图像的r、g、b值；再根据公式（1）计算叶片各部位颜色特征值mt，并求得mt均值，最后将mt均值代入预处理测量阶段总结的标准回归曲线方程中，得到叶片的衰老等级，例如若最后算得均值mt=0.63，则计算得出衰老等级为y=5.44，介于衰老等级5与等级6之间。

本发实施例中上述小麦叶片衰老程度快速测量方法，通过前期的预备总结得出回归曲线，为实际叶片快速测量做了充分总结，在实际测量中只需拍摄叶片图像，将颜色特征参数代入回归曲线方程中，即可很到叶片衰老等级，相比传统的需要测定多个与衰老程度相关的生理指标，具有省时省力，方便快速的优点，同时确定的测定结果准确率高，不受人为因素的影响。

本发明的小麦叶片衰老程度快速测量方法的并不局限于上述实施方式，同时也可应用于水稻、大麦等农植物叶片衰老过程的评价中。不同植物的衰老等级可能因植物生长特性而有所不同，不同植物的mt的计算公式和叶片衰老程度标准回归曲线会存在一定的差异，但以本发明为基础的修改和变化均属于本发明的保护范围。