CN103808263A - 谷粒粒形参数的高通量检测方法 - Google Patents

Abstract

本发明公开了一种谷粒粒形参数的高通量检测方法，包括以下步骤，搭建粒形检测平台；检测环境标定；影像采集与图像分割；轮廓提取与特征参数采集；利用粒形参数预设阈值进行粘连籽粒与杂质筛除；分析轮廓内图像，实现高通量颜色及纹理特征检测；输出谷物粒形参数检测结果。该粒形参数检测方法，其检测结果准确可靠，相比人工测量方法而言，极大地提高了检测效率和稳定性。能够满足基因研究中对谷粒粒形参数的检测需求。实现简单，检测高效准确，达到了应用的要求。

Description

谷粒粒形参数的高通量检测方法

技术领域

本发明涉及农作物领域，具体地，涉及一种谷粒粒形参数的高通量检测方法。

背景技术

目前，谷物籽粒的形状和大小是最重要的农业性状之一，它们不仅会影响产量和食用品质，还是市场价格的重要决定因素。因此，诸如遗传学、基因功能分析、分子育种和作物改良等植物学研究领域，均需要定量化评价种子的形态特征。高效，可靠，高通量的谷物粒形表现型测量方法是必需的。

一般来说，谷物籽粒的形态特征主要通过两种方法测量。最简易的是国标GB/T 17891-1999中介绍的方法，“随机数取完整无损的精米10粒，平放于测量板上，按照头对头、尾对尾，不重叠、不留隙的方式，紧靠直尺摆成一行，数出长度，求平均即为籽粒长度。同样的，按照同一方向肩靠肩排列，用直尺测量，读出宽度，求其平均即为籽粒宽度”，然后推算其他粒形参数。然而，这种手工方法能够获取的参数非常有限，同时测量的质量受籽粒形态变化以及测量人的主观因素等多种外界条件影响，并不稳定可靠。相比之下，利用数字成像技术获取籽粒的高分辨率图像，继而以此计算获取籽粒形态参数，不仅能够使谷物粒形参数实现自动且高效检测，还能够获取更多更复杂的形态、纹理及颜色参数等。

专利号“CN101275824”中公开了一种大米粒形检测方法。利用扫描仪获取大米籽粒图像，以灰度迭代法分割背景与籽粒，米粒轮廓进行逐点检测获取粒长、粒宽及长宽比。该技术方案中，灰度迭代法为自动分割技术，籽粒分割精度及粒形轮廓精度无法保持稳定，且所获取的粒形参数较少。专利号为“CN101905215A”公开的数字化水稻考种机，着重以风选法分离实粒和空瘪粒，并利用机器视觉获取粒数与长宽等初级形态参数。该发明以甄别空瘪粒为主要目的，机型庞大复杂，所能获取的粒形参数也较少。同时还有一系列玉米果穗装置及方法，均利用机器视觉技术获取玉米果穗形态及穗行数、行粒数等考种参数。上述技术方案均需由光源、专用CCD相机和图像采集卡构成的图像采集装置实现籽粒图像采集和形态分析功能。且分析的粒形参数不够全面。

发明内容

本发明的目的在于，针对上述问题，提出一种谷粒粒形参数的高通量检测方法，以实现检测结果准确可靠且谷粒粒形参数的检测齐全的优点。

为实现上述目的，本发明采用的技术方案是：

一种谷粒粒形参数的高通量检测方法，包括以下步骤，

步骤1、搭建粒形检测平台，

在水平工作台上，固定数码单反相机；

上述数码单反相机的成像平面与工作台平行，并设定数码单反相机的高度；

在水平工作台平面上放置单色背景的托盘，以承载待测籽粒；

将数码单反相机与计算机连接；

开启数码单反相机的实时取景模式，显示相机实时成像信息；

步骤2、检测环境标定；

以直尺作为参照物，放置于托盘内；

控制数码单反相机拍摄参照物影像，并将拍摄的影像回传至计算机；

计算机对拍摄的影像进行处理，即在影像中选择参照物标尺线段并输入实际尺寸，同时设定影像的像素数与实际尺寸的转换系数；

以上述转换系数标定数码单反相机与托盘的位置关系，准确托盘内物体的实际距离；

步骤3、影像采集与图像分割；

将待测籽粒均匀地散布于托盘内，避免大片籽粒相粘连；

开启数码单反相机实时取景模式，显示相机实时成像的籽粒图像；

控制数码单反相机拍摄谷物籽粒图像，实时回传至计算机；

使用计算机对采集的图像进行分割，将采集的谷物籽粒图像转换到HSL色彩空间，调整色调、亮度及饱和度阈值，实时观察分割效果，并设定分割阈值进行背景分割； 

步骤4、轮廓提取与特征参数采集；

利用八邻域寻迹算法获取各籽粒轮廓，以籽粒轮廓信息提取籽粒形态特征参数；

步骤5、利用粒形参数预设阈值进行粘连籽粒与杂质筛除；

步骤6、分析轮廓内图像，实现高通量颜色及纹理特征检测；

步骤7、输出谷物粒形参数检测结果。

根据本发明的优选实施例，上述步骤3中的图像分割，可直接载入预设阈值进行分割。

根据本发明的优选实施例，所述步骤4中所述特征参数采集具体包括：

求取籽粒轮廓周长、面积及重心坐标；

获取籽粒周长参数p及面积参数a，

求取复杂度com，                                                ，

形状参数sha ，

，

离散度dis ，；

以籽粒轮廓点集求取凸包点集；

以轮廓凸包点集求取籽粒轮廓最小外接圆，获取最小外接圆半径r，

求取圆形度cir， 

，

紧密度tig，

；

以轮廓凸包点集求取籽粒轮廓最小外接矩形，获取最小外接矩形面积；

求取矩形度

；

以轮廓凸包点集求取籽粒长宽；

搜索凸包点集，以最长轴为籽粒长度l并将凸包点集分为左右两部分，左、右点集分别求取至最长轴的最大距离

和

，则籽粒宽度W，

，

继而求取长宽比lwr 

。

本发明的技术方案具有以下有益效果：

本发明的技术方案，提出的粒形参数检测方法，其检测结果准确可靠，相比人工测量方法而言，极大地提高了检测效率和稳定性。能够满足基因研究中对谷粒粒形参数的检测需求。实现简单，检测高效准确，达到了应用的要求。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

图1为本发明实施例所述的谷粒粒形参数的高通量检测结构示意图；

图2为本发明实施例所述的谷粒粒形参数的高通量检测方法流程图；

图3为图像处理软件的像素距离到实际距离标定对话框示意图；

图4为图像分割阈值设定UI界面示意图；

图5为籽粒轮廓提取与标号结果示意图；

图6为杂点清除设置界面示意图；

图7a至图7d为谷物籽粒形态特征提取效果示意图；

图8a为籽粒样本人工测量值与检测值籽粒长度比较示意图；

图8b为籽粒样本人工测量值与检测值籽粒宽度比较示意图。

结合附图，本发明实施例中附图标记如下：

1-翻拍架；2-数码单反相机；3-相机连接线；4-计算机；5-参照物；6-待测籽粒；7-托盘；。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

实施例一、

一种谷粒粒形参数的高通量检测方法，包括以下步骤，

步骤1、搭建粒形检测平台，如图1所示，

在水平工作台上，借助翻拍架1或三脚架等固定数码单反相机2；

上述数码单反相机2的成像平面与工作台平行，并设定相机的高度；设定的高度以获得较理想的成像为标准。

在水平工作台平面上放置单色背景的托盘7，以承载待测籽粒；托盘7颜色以与籽粒颜色有较大反差为宜；

将数码单反相机2与计算机4通过相机连接线3连接；

开启数码单反相机的实时取景模式，显示相机实时成像信息；

步骤2、检测环境标定；

以直尺作为参照物5，放置于托盘内；

控制数码单反相机拍摄参照物影像，并将拍摄的影像回传至计算机；

计算机对拍摄的影像进行处理，即在影像中选择参照物标尺线段并输入实际尺寸，同时设定影像的像素数与实际尺寸的转换系数；

以上述转换系数标定数码单反相机与托盘的位置关系，准确托盘内物体的实际距离；

步骤3、影像采集与图像分割；

将待测籽粒均匀地散布于托盘内，避免大片籽粒相粘连，以提高籽粒轮廓提取的准确性；

开启数码单反相机实时取景模式，显示相机实时成像的籽粒图像；

控制数码单反相机拍摄谷物籽粒图像，实时回传至计算机；

使用计算机对采集的图像进行分割，将采集的谷物籽粒图像转换到HSL色彩空间，调整色调、亮度及饱和度阈值，实时观察分割效果，并设定分割阈值进行背景分割； 

步骤4、轮廓提取与特征参数采集；

利用八邻域寻迹算法获取各籽粒轮廓，以籽粒轮廓信息提取籽粒形态特征参数；

步骤5、利用粒形参数预设阈值进行粘连籽粒与杂质筛除；

步骤6、分析轮廓内图像，实现高通量颜色及纹理特征检测；

步骤7、输出谷物粒形参数检测结果。

根据本发明的优选实施例，上述步骤3中的图像分割，可直接载入预设阈值进行分割。

步骤4所提到特征参数求包括：

（1）求取籽粒轮廓周长、面积及重心坐标；

获取籽粒周长参数p及面积参数a，

求取复杂度com， ，

形状参数sha ，

，

离散度dis ，

；

（2）以籽粒轮廓点集求取凸包点集；将轮廓点集转化为凸包点集可以极大地减少后续分析过程的计算量。

（3）以轮廓凸包点集求取籽粒轮廓最小外接圆，获取最小外接圆半径r，

求取圆形度cir， 

，

紧密度tig，；

（4）以轮廓凸包点集求取籽粒轮廓最小外接矩形，获取最小外接矩形面积

；

求取矩形度

；

（5）以轮廓凸包点集求取籽粒长宽；

搜索凸包点集，以最长轴为籽粒长度l并将凸包点集分为左右两部分，左、右点集分别求取至最长轴的最大距离

和，则籽粒宽度W，

，

继而求取长宽比lwr 。

为了更好地实现上述方案，各环节计算结果均实时显示于粒形参数检测软件的主界面中，便于控制检测精度。粒形参数检测完成后，将同时保存原始图像、分割结果、轮廓提取及标号结果以及粒形参数检测结果等信息。结果保存后，计算机切换至实时取景模式，等待下次采集的籽粒图像进行分析处理。

实施例二、

1、按照本发明技术方案提供的规范构建拍摄环境，开启计算机内图像处理软件（下文简称软件）并与相机建立连接，软件切换至“实时取景”模式。

2、将标定参照物放置于谷物托盘中，控制相机拍摄参照物图像；参照物图像回传至计算机，软件切换至“分析”模式。

3、选择标定功能，在参照物图像中左键点击标尺起始点；移动鼠标至标尺结束点，再次点击鼠标左键；此时参照物图像中绘制了起始点与结束点像素连线并显示了两点间的像素数。如果不满意本次选取，可以点击右键取消，重新点击距离标定按钮进行新标定；如果满意选取则点击鼠标左键，软件开启距离标定设定UI。

4、在距离标定UI如图3所示，中填写参照物起始点与结束点所代表的实际距离（mm），软件将计算像素到实际尺寸的转换率，点击确定以设定该转换率。最后将参照物取出。

5、开启“实时取景”模式，在谷物托盘放入待测籽粒样本；观察软件主显示区域的实时图像采集结果，如果成像区域内籽粒清晰、粘连少、且无籽粒处于图像采集的边缘，则控制相机进行籽粒图像的采集；图像回传至计算机，软件将切换至“分析”模式，主显示界面呈现所采集的籽粒图像。

6、选择“HSL图像分割”功能，开启HSL色彩空间滤色器UI界面如图4所示。分别设置色相、饱和度和亮度阈值，通过预览窗口观察籽粒目标分割效果，满意后保存阈值并进行图像分割。

7、选择“特征提取”功能，将对分割得到的籽粒目标区域进行标号并获取其轮廓；通过轮廓点集，继而得到籽粒的形态特征参数。边界寻迹结果将叠加到原图显示到主工作区，籽粒的右下角显示区域标号如图5所示，将对应分析结果文件中的每条结果记录。当鼠标移动到籽粒区域时，籽粒区域的边界将高亮显示，点击选定籽粒则显示该籽粒的详细形态参数等。

8、选择“清除杂点”功能，可以根据籽粒的像素信息或形态信息进行阈值筛选如图6所示，排除灰尘、背景亮斑、粘连籽粒等杂点。

9、选择“保存”功能，保存原始图像、分割结果、轮廓提取及标号结果以及粒形参数检测结果等信息。

10、软件切换至“实时取景”模式，等待进行下一组样本检测。

实施例三、

1、利用以标定好的谷物粒形参数检测环境，将样本籽粒散布于谷物托盘中，通过“实时取景”功能观察待测样本，当样本分布符合检测要求时，控制相机采集籽粒图像并回传至计算机，检测软件切换至“分析”模式。

2、选择“参数检测”功能，软件将载入软件保存的转换率和分割阈值，进行图像分割；分割得到的籽粒目标区域进行标号并获取其轮廓，继而获取籽粒的形态特征参数；并根据设定阈值进行杂点筛选。

3、参数检测结果将与采集得到的籽粒图像叠加显示到主工作区，籽粒的右下角显示区域标号，将对应分析结果文件中的每条结果记录。当鼠标移动到籽粒区域时，籽粒区域的边界将高亮显示，点击选定籽粒则显示该籽粒的详细形态参数等，便于对检测准确性进行评估。如果不满意自动检测结果，还可以按照实例二的6-8步骤进行有人工干预的参数检测。

4、选择“保存”功能，保存原始图像、分割结果、轮廓提取及标号结果以及粒形参数检测结果等信息，软件切换至“实时取景”模式，等待进行下一组样本检测。

现代植物基因组学的飞速发展，对籽粒形态特征提取的便捷性和精度性提出了更高的要求。本发明技术方案以现代实验室常备仪器设备为基础，通过操作方法和检测算法的明创新，实现了准确谷物粒形参数的高通量检测方法。本发明技术方案的关键点包括：

利用参照物（如直尺等）对检测人员自己构建的检测环境进行标定，使软件能够准确、可靠的获取相机成像平面内的籽粒实际尺寸。按照本发明制定的规范进行设备组装和检测环境标定方法，能够满足谷物籽粒形态特征检测工作的稳定性要求。

本发明提出利用数码单反相机所提供的控制接口，通过计算机直接控制相机进行谷粒粒形图像的拍摄与粒形参数的检测分析工作。以“实时取景”模式进行采样准备，拍摄籽粒图像并回传至计算机后，切换至“分析”模式，进行图像分割、轮廓提取与标号和粒形参数检测。提高了检测效率。

将拍摄的籽粒图像转换到HSL色彩空间，再根据阈值进行图像分割运算，减弱了光线变化对图像分割精度的影响，使本发明无需额外配置光源环境，在室内条件下即能完成粒形参数的检测工作。同时，图像分割阈值设定提供了良好的交互环境，可实时观测分割效果，提高了籽粒目标的分割精度。分割结果将进行若干次膨胀和腐蚀运算，去除小花柄、细小杂质或籽粒轻微粘连的影响。

通过实现上述两项技术途径，本发明提供的检测方法无需设计制造专用的粒形检测装置，利用实验室环境中常见的仪器设备进行简单的组装，并按照检测标准进行粒形检测，既可实现粒形参数的高通量检测工作。

在特征检测运算中，利用八邻域寻迹算法进行籽粒轮廓提取，提高了轮廓提取的精细度。在形态参数计算过程中，先求取轮廓点集的凸包，继而以凸包点集进行最小外接矩形、最小外接圆以及籽粒长宽等形态参数的计算，可以在不影响结果准确性的前提下，极大地降低运算量、提高运算速度如图7a至图7d所示。其中图7a为籽粒轮廓，图7b为凸包，图7c为最小外接矩形，图7d为最小外接圆。进行上述运算后，不仅能够提供更为准确的长、宽以及长宽比数据，还能提供籽粒面积、周长、圆度、矩形度、紧密度、复杂度、形状参数和离散度等更为复杂的形态指标。本发明可以对成像区域内的所有籽粒进行粒形参数的分析检测并输出，便于对谷物样本进行更细致深入地统计分析。

本发明利用现代实验室中常备的数码单反相机、三脚架、直尺和计算机等仪器设备，构建谷物粒形参数高通量检测环境。能够减少设备投资，提高现有设备的利用率和使用灵活性。

本发明的检测结果与用游标卡尺人工测量结果的比较分析如表1所示，对于长宽结果，两者间具有极显著（P<0.01）的相关关系，且两种方法所取得的结果差异不显著（P>0.05）。本发明的检测结果与人工测量结果间的关系如图8a和图8b所示。

表1、本发明技术方案检测值与人工测量值比较分析结果表：

。

最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种谷粒粒形参数的高通量检测方法，其特征在于，包括以下步骤，

步骤1、搭建粒形检测平台，

在水平工作台上，固定数码单反相机；

上述数码单反相机的成像平面与工作台平行，并设定数码单反相机的高度；

在水平工作台平面上放置单色背景的托盘，以承载待测籽粒；

将数码单反相机与计算机连接；

开启数码单反相机的实时取景模式，显示相机实时成像信息；

步骤2、检测环境标定；

以直尺作为参照物，放置于托盘内；

控制数码单反相机拍摄参照物影像，并将拍摄的影像回传至计算机；

计算机对拍摄的影像进行处理，即在影像中选择参照物标尺线段并输入实际尺寸，同时设定影像的像素数与实际尺寸的转换系数；

以上述转换系数标定数码单反相机与托盘的位置关系，准确托盘内物体的实际距离；

步骤3、影像采集与图像分割；

将待测籽粒均匀地散布于托盘内，避免大片籽粒相粘连；

开启数码单反相机实时取景模式，显示相机实时成像的籽粒图像；

控制数码单反相机拍摄谷物籽粒图像，实时回传至计算机；

使用计算机对采集的图像进行分割，将采集的谷物籽粒图像转换到HSL色彩空间，调整色调、亮度及饱和度阈值，实时观察分割效果，并设定分割阈值进行背景分割； 

步骤4、轮廓提取与特征参数采集；

利用八邻域寻迹算法获取各籽粒轮廓，以籽粒轮廓信息提取籽粒形态特征参数；

步骤5、利用粒形参数预设阈值进行粘连籽粒与杂质筛除；

步骤6、分析轮廓内图像，实现高通量颜色及纹理特征检测；

步骤7、输出谷物粒形参数检测结果。

2.根据权利要求1所述的谷粒粒形参数的高通量检测方法，其特征在于，上述步骤3中的图像分割，可直接载入预设阈值进行分割。

3.根据权利要求1或2所述的谷粒粒形参数的高通量检测方法，其特征在于，所述步骤4中所述特征参数采集具体包括：

求取籽粒轮廓周长、面积及重心坐标；

获取籽粒周长参数p及面积参数a，

求取复杂度com，                                                

，

形状参数sha ，

，

离散度dis ，；

以籽粒轮廓点集求取凸包点集；

以轮廓凸包点集求取籽粒轮廓最小外接圆，获取最小外接圆半径r，

求取圆形度cir， ，

紧密度tig，

；

以轮廓凸包点集求取籽粒轮廓最小外接矩形，获取最小外接矩形面积

；

求取矩形度

；

以轮廓凸包点集求取籽粒长宽；

搜索凸包点集，以最长轴为籽粒长度l并将凸包点集分为左右两部分，左、右点集分别求取至最长轴的最大距离

和

，则籽粒宽度W，

，

继而求取长宽比lwr 

。

Images