基于图像处理的农作物植株水分胁迫表型检测方法
【专利摘要】本发明提供了一种基于图像处理的农作物植株水分胁迫表型检测方法，该方法包括：S1：将农作物的图像分割为植株图案和背景图案；S2：对所述植株图案进行分割，得到茎秆图案和叶片图案；S3：将茎秆图案细化为茎秆直线，将叶片图案细化为叶片曲线，计算直线OS与直线OP1所形成的夹角SOP1的角度值以及直线OS与直线OP2所形成的夹角SOP2的角度值；S4：计算夹角SOP1的角度值与夹角SOP2的角度值的比值，并根据所述比值对所述农作物的叶片受水平胁迫时的变化进行评估。本发明提供的基于图像处理的农作物植株水分胁迫表型检测方法，能够对作物细微变化进行精确感知和监测，并准确反映作物与环境之间的互作关系。
【专利说明】
基于图像处理的农作物植株水分胁迫表型检测方法
技术领域
[0001]本发明涉及农业信息技术领域，尤其涉及一种基于图像处理的农作物植株水分胁 迫表型检测方法。
【背景技术】
[0002]植物表型（phenotype)被定义为植物基因型（genotype)和所处环境决定的形状、 结构、大小、颜色等全部可测的生物体外在表现。即表型是一个基因型与环境互作产生的全 部或部分可辨识特征和性状。从功能基因组学到作物栽培生理的各个层面的研究都涉及到 对植株的各种特征和性状即表型的鉴别与分析，以及对复杂的植物生长环境的监测与控 制。基于图像的植物表型检测能够直接反应物体的形貌特征，具有非侵入式和高通量等特 点，已经逐渐应用于植物表型分析中，然而，目前的植物表型检测方法需要一定的人工交 互，自动化程度不够，对作物细微变化无法精确感知和监测，不能准确反映作物与环境之间 的互作关系。
【发明内容】
[0003] (一)要解决的技术问题
[0004] 本发明要解决的技术问题是:提供一种基于图像处理的农作物植株水分胁迫表型 检测方法，能够对作物细微变化进行精确感知和监测，并准确反映作物与环境之间的互作 关系。
[0005] (二)技术方案
[0006] 为解决上述技术问题，本发明的技术方案提供了一种基于图像处理的农作物植株 水分胁迫表型检测方法，包括：
[0007] S1:将农作物的图像分割为植株图案和背景图案；
[0008] S2:对所述植株图案进行分割，得到茎杆图案和叶片图案；
[0009] S3:将茎杆图案细化为茎杆直线，将叶片图案细化为叶片曲线，记所述茎杆直线与 所述叶片曲线的交叉点为位置点〇,在茎杆直线上选取位置点S，在叶片曲线上选取位置点 P4P位置点p2,计算直线0S与直线(^所形成的夹角SOPi的角度值以及直线0S与直线0P2所形 成的夹角S0P 2的角度值；
[0010] S4:计算夹角SOP^角度值与夹角S0P2的角度值的比值，并根据所述比值对所述农 作物的叶片受水平胁迫时的变化进行评估。
[0011] 优选地，采集不同时间所述农作物的图像，对每一副所述农作物的图像分别进行 步骤S1~S4，以衡量所述农作物应对环境变化的自身调整能力。
[0012] 优选地，步骤S2包括：
[0013] S21:扫描所述植株图案的第一行像素，将超过预设个数植株像素的联通区域作为 候选区域，并对每一个候选区域进行标记，赋予不同候选区域不同的标记值；
[0014] S22:对所述植株图案的第二行像素至最末行像素依次执行以下处理：
[0015] 扫描该行像素，将超过所述预设个数植株像素的联通区域作为新候选区域，计算 每一个新候选区域与每一个已标记的候选区域的联通关系，对于每一个新候选区域，若其 与1个已标记的候选区域存在联通关系，则该新候选区域继承所述1个已标记的候选区域的 标记值，若其与0个已标记的候选区域存在联通关系，则对该新候选区域赋予新的标记值， 若其与2个已标记的候选区域存在联通关系，则对该新候选区域进行拆分，并将拆分后得到 的一个区域赋予所述2个已标记的候选区域中一个已标记的候选区域的标记值，将拆分后 得到的另一个区域赋予所述2个已标记的候选区域中另一个已标记的候选区域的标记值；
[0016] S23:判断每一个标记值对应的所有像素所组成的形状，从中选取对应的所有像素 所组成形状符合矩形分布的标记值，将所述选取的标记值对应的所有像素所组成的区域作 为茎杆图案，将余下的标记值所对应的像素进行合并处理，得到叶片图案。
[0017]优选地，所述对该新候选区域进行拆分包括：
[0019] 其中，CL为该新候选区域的像素数量，CLp1、CLp2分别为所述2个已标记的候选区域 的像素数量，CLi、CL 2分别为拆分后两个区域的像素数量。
[0020] 优选地，在步骤S3中，位置APi与所述茎杆直线的距离为，位置点p2与所述茎杆
直线的距离为·
，其中，L为所述叶片曲线上对应叶尖的位置点与所述茎杆直线的距离。 [0021](三)有益效果
[0022]本发明提供的基于图像处理的农作物植株水分胁迫表型检测方法，能够对作物细 微变化进行精确感知和监测，并准确反映作物与环境之间的互作关系。
【附图说明】
[0023]图1是本发明实施方式提供的一种基于图像处理的农作物植株水分胁迫表型检测 方法的流程图；
[0024] 图2是本发明实施方式提供的对农作物图像进行植株和背景分割后的示意图；
[0025] 图3是本发明实施方式提供的对植株图案进行器官分割后的示意图；
[0026] 图4是本发明实施方式提供的叶片角度计算示意图；
[0027] 图5是本发明实施方式提供的植株叶片角度随时间变化曲线示意图；
[0028] 图6是本发明实施方式提供的一天相对湿度变化曲线示意图。
【具体实施方式】
[0029]下面结合附图和实施例，对本发明的【具体实施方式】作进一步详细描述。以下实施 例用于说明本发明，但不用来限制本发明的范围。
[0030]参见图1，图1是本发明实施方式提供的一种基于图像处理的农作物植株水分胁迫 表型检测方法的流程图，包括：
[0031] S1:将农作物的图像分割为植株图案和背景图案；
[0032] S2:对所述植株图案进行分割，得到茎杆图案和叶片图案；
[0033] S3:将茎杆图案细化为茎杆直线，将叶片图案细化为叶片曲线，记所述茎杆直线与 所述叶片曲线的交叉点为位置点〇,在茎杆直线上选取位置点S，在叶片曲线上选取位置点 Pi和位置点P 2,计算直线0S与直线(^所形成的夹角SOPi的角度值以及直线0S与直线0P2所形 成的夹角S0P 2的角度值；
[0034] S4:计算夹角SOPj^角度值与夹角S0P2的角度值的比值，并根据所述比值对所述农 作物的叶片受水平胁迫时的变化进行评估。
[0035]优选地，可以采集不同时间所述农作物的图像，对每一副所述农作物的图像分别 进行上述步骤S1~S4,以衡量所述农作物应对环境变化的自身调整能力，例如，可以采集一 天时间内不同时间所述农作物的图像。
[0036] 其中，步骤S2可以包括：
[0037] S21:扫描所述植株图案的第一行像素，将超过预设个数植株像素的联通区域作为 候选区域，并对每一个候选区域进行标记，赋予不同候选区域不同的标记值；
[0038] S22:对所述植株图案的第二行像素至最末行像素依次执行以下处理：
[0039]扫描该行像素，将超过所述预设个数植株像素的联通区域作为新候选区域，计算 每一个新候选区域与每一个已标记的候选区域的联通关系，对于每一个新候选区域，若其 与1个已标记的候选区域存在联通关系，则该新候选区域继承所述1个已标记的候选区域的 标记值，若其与0个已标记的候选区域存在联通关系，则对该新候选区域赋予新的标记值， 若其与2个已标记的候选区域存在联通关系，则对该新候选区域进行拆分，并将拆分后得到 的一个区域赋予所述2个已标记的候选区域中一个已标记的候选区域的标记值，将拆分后 得到的另一个区域赋予所述2个已标记的候选区域中另一个已标记的候选区域的标记值； [0040] S23:判断每一个标记值对应的所有像素所组成的形状，从中选取对应的所有像素 所组成形状符合矩形分布的标记值，将所述选取的标记值对应的所有像素所组成的区域作 为茎杆图案，将余下的标记值所对应的像素进行合并处理，得到叶片图案。
[0041 ]其中，在步骤S22中，所述对该新候选区域进行拆分包括：
[0043]其中，CL为该新候选区域的像素数量，CLp1、CLp2分别为所述2个已标记的候选区域 的像素数量，CLi、CL2分别为拆分后两个区域的像素数量。
[0044] 其中，在步骤S3中，位置APi与所述茎杆直线的距离为\
，位置点P2与所述茎杆直 线的距离为
，其中，L为所述叶片曲线上对应叶尖的位置点与所述茎杆直线的距离。
[0045] 本发明实施方式提供的基于图像处理的农作物植株水分胁迫表型检测方法，能够 对作物细微变化进行精确感知和监测，并准确反映作物与环境之间的互作关系。
[0046] 例如，采用上述方法可以对玉米进行植株水分胁迫表型检测，在水分胁迫条件下， 玉米叶片含水量下降，细胞膨压发生变化，植株产生应激反应，叶片气孔的保卫细胞伸长， 气孔开度降低，叶片发生卷曲减小蒸散面积降低水分散失。叶片形态变化反应了对水分胁 迫应激反应的灵敏程度，从一定程度上可以表征植株抗水分胁迫的能力，检测方法包括：
[0047] 1、连续图像序列获取与预处理；
[0048]检测对象可以为转抗旱基因的玉米植株，获取图像前14天无灌溉，土壤体积含水 量18.4%，可将工业相机安装在一个支架上，被测植株处于相机与幕布之间，尽量使相机图 像平面、玉米植株平面、幕布三者平行。从上午8:00开始到午后17:00为止，间隔10分钟程序 控制自动获取一幅图像，共获取54张，图像分辨率4000 X 3000像素，然后对图像序列采用基 于机器学习的背景分割方法，将图像中的内容分割为植株和背景，分割结果如图2所示。
[0049] 2、植株器官分割
[0050] 对于预处理后的每一副图像，检测叶片形态变化，首先需要分割植株上叶片和茎 杆等器官，玉米叶片与莖杆的颜色纹理十分相近，图割(graph cut)和水平集（level set) 等广义的图像分割方法很难获得理想的分割结果。玉米植株的形态有一定规律，叶片围绕 茎杆在植株两侧交替分布，利用这种形态特点作为先验知识，提出基于扫描线的叶片-茎杆 分割算法，算法具体步骤如下：
[0051] 1)扫描图像第一行像素，将超过3个植株像素的联通区域标记为叶片/茎杆候选区 域，并对每一个候选区域进行标记，赋予不同候选区域不同的标记值；
[0052] 2)扫描第二行像素，获得新候选区域，计算每一个新候选区域与步骤1)中已标记 的候选区域的联通关系(相邻关系），对于每一个新候选区域，若其仅与1个已标记的候选区 域存在联通关系，则说明该新候选区域与该1个已标记的候选区域属同一个器官，继承该1 个已标记的候选区域的标记值，若其与0个已标记的候选区域存在联通关系，即没有已标记 的候选区域与其存在联通关系，则说明该新候选区域是一个独立的器官，赋予该候选区域 一个新的标记值，若其与2个已标记的候选区域存在联通关系，说明该区域为两个器官的连 接处，采用下式对该新候选区域进行像素拆分后分别标记：
[0054]其中，CL表示该新候选区域的像素数量，0^、0^2表示与该新候选区域存在联通 关系的2个候选区域的像素数量，CL^CLs为拆分后两个区域的像素数量；
[0055] 3)重复步骤2)，直到扫描完整个图像，判断每个标记值所对应的像素集合中所有 像素组成的形状，将其中符合矩形分布的一组像素集合作为为茎杆区域，该集合中像素的 X-方向均值与植株上所有像素X-方向均值的距离应在特定阈值内（如30pixel);
[0056] 4)对余下的标记值所对应的像素集合进行合并，由于玉米叶片边缘常伴有一些凸 起和褶皱，导致一些本应属于一片叶子的像素集合，被标记为几个不同的叶片像素集合，例 如，可以按照联通的像素比例对不同像素集合进行合并，若一个标记值所对应的像素集合 中超过设定比例的像素与另一个标记值所对应的像素集合相联通，则合并这2个标记值所 对应的像素集合，例如，一个标记值所对应的像素集合由10行像素组成，如果有7行(70 % ) 像素与一个标记值所对应的像素集合相联通，则合并这2个标记值所对应的像素集合，重复 上述操作直到所有能够合并的像素集合被合并为止，分割结果如下图3所示。
[0057] 3、植株水分胁迫表型指标计算
[0058]采用两个角度值做为植株水分胁迫表型指标，如图4所示，将茎杆图案和叶片图案 分别细化为茎杆直线和叶片曲线，茎杆直线0S，叶片尖端点距离茎杆直线的距离记为UPi、 p2点分别分别为L的1/4和2/3，S是茎杆直线上一点，0是茎杆直线与叶片曲线的交点，计算 S0P^S0P 2的角度值，并以S0PjPS0P2的比值作为叶片形态变化的衡量指标，SOPi/SOP〗表征 了叶片受到水分胁迫时的弯曲变化程度，该值计算较为简便，对于不同茎叶夹角的叶片该 值具有可比性。
[0059] 4、植株水分胁迫抗表型指标与环境因子的协同变化趋势
[0000] SOP1/SOP2表征了叶片一天内弯曲程度的变化，该值可做为植株水分胁迫表型指 标，用来衡量植株应对环境变化的自身调整能力。如图5所示，一天内SOPi/SOP〗变化趋势，随 着时间推移，中午达到最大数值，午后数值逐渐下降。
[0061] 玉米植株的水分胁迫表型是基因型与环境互作的结果，在土壤含水量一定的条件 下，空气相对湿度是最主要的影响因素，当天的空气相对湿度变化曲线如图6所示，相对湿 度同水分胁迫表型参数指标(S0PVS0P2)呈现较为一致的协同变化趋势，证明S0PVS0P2是 能够表征植株水分胁迫的表型指标。
[0062] 通过上述1天内的连续图像序列，发现水分胁迫玉米植株叶片形态细微变化，采用 图像处理算法计算了表征这种变化的玉米植株水分胁迫表型参数，与同步测量的环境数据 对照(空气相对湿度），二者具有较为一致的协同变化规律。
[0063] 以上实施方式仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通 技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有 等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定。
【主权项】
1. 一种基于图像处理的农作物植株水分胁迫表型检测方法，其特征在于，包括： S1:将农作物的图像分割为植株图案和背景图案； S2:对所述植株图案进行分割，得到茎杆图案和叶片图案； S3:将茎杆图案细化为茎杆直线，将叶片图案细化为叶片曲线，记所述茎杆直线与所述 叶片曲线的交叉点为位置点0,在茎杆直线上选取位置点S，在叶片曲线上选取位置点Pi和 位置点P2,计算直线0S与直线OPi所形成的夹角SOPi的角度值W及直线0S与直线0P2所形成 的夹角S0P2的角度值； S4:计算夹角SOPi的角度值与夹角S0P2的角度值的比值，并根据所述比值对所述农作物 的叶片受水平胁迫时的变化进行评估。2. 根据权利要求1所述的基于图像处理的农作物植株水分胁迫表型检测方法，其特征 在于，采集不同时间所述农作物的图像，对每一副所述农作物的图像分别进行步骤S1~S4， W衡量所述农作物应对环境变化的自身调整能力。3. 根据权利要求1所述的基于图像处理的农作物植株水分胁迫表型检测方法，其特征 在于，步骤S2包括： S21:扫描所述植株图案的第一行像素，将超过预设个数植株像素的联通区域作为候选 区域，并对每一个候选区域进行标记，赋予不同候选区域不同的标记值； S22:对所述植株图案的第二行像素至最末行像素依次执行W下处理： 扫描该行像素，将超过所述预设个数植株像素的联通区域作为新候选区域，计算每一 个新候选区域与每一个已标记的候选区域的联通关系，对于每一个新候选区域，若其与1个 已标记的候选区域存在联通关系，则该新候选区域继承所述1个已标记的候选区域的标记 值，若其与0个已标记的候选区域存在联通关系，则对该新候选区域赋予新的标记值，若其 与2个已标记的候选区域存在联通关系，则对该新候选区域进行拆分，并将拆分后得到的一 个区域赋予所述2个已标记的候选区域中一个已标记的候选区域的标记值，将拆分后得到 的另一个区域赋予所述2个已标记的候选区域中另一个已标记的候选区域的标记值； S23:判断每一个标记值对应的所有像素所组成的形状，从中选取对应的所有像素所组 成形状符合矩形分布的标记值，将所述选取的标记值对应的所有像素所组成的区域作为茎 杆图案，将余下的标记值所对应的像素进行合并处理，得到叶片图案。4. 根据权利要求3所述的基于图像处理的农作物植株水分胁迫表型检测方法，其特征 在于，所述对该新候选区域进行拆分包括：其中，化为该新候选区域的像素数量，CLpi、CLp2分别为所述2个已标记的候选区域的像 素数量，化1、CL2分别为拆分后两个区域的像素数量。5. 根据权利要求1所述的基于图像处理的农作物植株水分胁迫表型检测方法，其特征 在于，在步骤S3中，位置点Pi与所述茎杆直线的距离为位置点P2与所述茎杆直线的距离 4 为^，其中，L为所述叶片曲线上对应叶尖的位置点与所述茎杆直线的距离。
【文档编号】G06T7/00GK106097372SQ201610465992
【公开日】2016年11月9日
【申请日】2016年6月23日
【发明人】王传宇, 郭新宇, 杜建军, 温维亮, 吴升
【申请人】北京农业信息技术研究中心