一种基于机器视觉的胡萝卜表面缺陷识别量化方法与流程

本发明涉及农产品加工技术领域，具体涉及一种基于机器视觉的胡萝卜表面缺陷识别量化方法。

背景技术：

胡萝卜的分级销售有助于提升胡萝卜的市场竞争力和销售价格，提高胡萝卜生产企业利润。目前胡萝卜分选主要依靠人工完成，这种分选方式效率低，主观性强、标准不严格，具有人工分选的固有缺陷。并且随着用人成本的上升，进一步压缩了胡萝卜加工企业的利润。人工分选已经不能满足胡萝卜加工企业的生产要求。

机器视觉提供了一种高效、低成本和实时的果蔬在线检测、分级方法，采用机器视觉对胡萝卜进行在线分级，可以节省劳动力，提高识别准确性，保证胡萝卜分级的一致性，为胡萝卜生产企业降低生产成本，提高行业竞争力。目前有大量的运用机器视觉对果蔬分级的研究，但其主要适用于球形或类球形果蔬，如苹果、桃、梨、柑橘和马铃薯等。这些技术在长形果蔬上的应用具有局限性，因此需要研究开发适合胡萝卜分选的机器视觉方法。

目前国内还没有相对成熟的基于机器视觉的胡萝卜分选方法，根据中华人民共和国国内贸易行业标准(sb/t10450-2007)、胡萝卜购销等级要求(中华人民共和国商务部2007.12.28发布)，影响胡萝卜购销等级的关键指标有歪扭、弯曲、开裂、青头和病虫伤等，其衡量标准比较模糊，为胡萝卜缺陷检测带来了困难，需要对胡萝卜的缺陷指标进行量化。

技术实现要素：

现有的胡萝卜分级设备多根据胡萝卜的粗细不同，采用不同间隙的滚轴等机械方式使胡萝卜按照一定的粗度分级，虽然这种分级方式极大提高了生产效率，节省人工成本，但不能把有缺陷(青头、断裂、弯曲、开裂、须根等)的胡萝卜识别出来，影响分级效果，因此不适合大规模推广应用。本发明公开的胡萝卜表面缺陷检测方法可提供一种基于机器视觉的胡萝卜断裂、弯曲、开裂、须根、虫眼五种缺陷识别量化方法。

为达到以上目的，本发明采取的技术方案是：

一种基于机器视觉的胡萝卜表面缺陷识别量化方法，包括：胡萝卜图像的预处理；胡萝卜断裂缺陷识别量化方法、胡萝卜弯曲缺陷识别量化方法、胡萝卜开裂缺陷识别量化方法、胡萝卜须根缺陷识别量化方法和胡萝卜虫眼缺陷识别量化方法；

所述胡萝卜图像的预处理包括以下步骤：

(1)将胡萝卜的rgb图像转换为hsv图像，提取出h、s、v三个分量图像，对h、s、v三个分量图像分别选取阈值使背景变为黑色，得到消除背景的二值图像；

(2)从胡萝卜源图像中获取r、g、b三个分量图像，并根据公式①、②对r、g、b三个分量图像进行处理，得到灰度图像gray和灰度图像gray，再利用全局阈值对灰度图像gray和灰度图像gray分别进行二值化，消除复杂背景影响，得到灰度图像gray、灰度图像gray对应的二值图像；

gray＝g-b①

gray＝2*r-g-b②；

所述胡萝卜断裂缺陷识别量化方法，包括如下步骤：

(1)根据胡萝卜的h分量图像得到二值图像；

(2)获得二值图像的边缘；

(3)根据二值图像的边缘，获取二值图像的最小外接矩形，得到最小外接矩形的长l和宽w；

(4)计算胡萝卜的纵横比zr＝w/l；

(5)提取出二值图像两端的边缘点坐标，保存在矩阵q1和矩阵q2中；

(6)分别对矩阵q1、矩阵q2求斜率，得到斜率k1和斜率k2；

(7)根据zr、k1和k2，确定胡萝卜是否是真的断裂，防止误判并判断出断裂曲面的形状；

当zr大于阈值0.25时且k1或k2中出现多于5个无穷大时，胡萝卜为断裂，且zr越大断裂越严重，当k1或k2中出现无穷大次数越多时，表示断面越陡峭；

所述胡萝卜弯曲缺陷识别量化方法，包括如下步骤：

(1)根据胡萝卜的h分量图像得到二值图像；

(2)对得到的二值图像进行面积去噪，去掉小的斑点；

(3)用圆形腐蚀单元对二值图像进行腐蚀，消除边缘毛刺；

(4)在二值图像中提取出连通像素面积在某一范围内的区域，即为胡萝卜主体区域；

(5)抽取胡萝卜主体区域的骨架，获得骨架图像；

(6)对骨架图像进行hough变换，得到骨架图像的两条拟合直线，两条拟合直线的夹角分别为θ1，θ2；

(7)根据θ1，θ2计算得到弯曲度br，br＝180/π*|θ1-θ2|；

当弯曲度br大于阈值5时，表示胡萝卜弯曲，且弯曲度br越大，表示胡萝卜弯曲程度越大；

所述胡萝卜开裂缺陷识别量化方法，包括如下步骤：

(1)根据胡萝卜的r、g、b三个分量图像，通过式①得到灰度图像gray，对灰度图像gray进行二值化，得到二值图像；

(2)在二值图像中提取出连通像素面积在某一范围(如[200050000])内的区域；

(3)获得该区域的等价椭圆，得到等价椭圆的长半轴a和短半轴b以及区域面积a2；

(4)根据椭圆面积公式s＝π*a*b计算出等价椭圆的面积；

(5)计算区域面积和等价椭圆面积之比i＝a2/s，长半轴和短半轴之比a/b；

(6)开裂度用长半轴与短半轴构成三角形的夹角表示，kr＝2*180/π*arctan(b/a)；

当a/b大于5时且i＝a2/s大于阈值0.8时，判断其为开裂区域，开裂度kr越大说明开裂越严重；

所述胡萝卜须根缺陷识别量化方法，包括如下步骤：

(1)根据胡萝卜的s分量图像，得到二值图像；

(2)对二值图像进行内部孔洞的填充；

(3)获取二值图像的边缘点坐标，保存在矩阵p中；

(4)从矩阵p中提取出上边缘的坐标存在矩阵p1中，从矩阵p提取出下边缘坐标存在矩阵p2中；

(5)分别对矩阵p1和矩阵p2每隔50个像素进行取样，取样后得到矩阵p3和矩阵p4；

(6)分别对矩阵p3、矩阵p4求斜率，得到斜率k3、斜率k4，求取k3和k4的峰值，当有一个峰值大于阈值1.5时，则胡萝卜存在一处须根，当有若干个峰值大于阈值1.5时，则胡萝卜存在若干处须根，找到峰值对应的点，并记录在二值图像中，求取坐标，该坐标对应胡萝卜须根的位置；

所述胡萝卜虫眼缺陷识别量化方法，包括如下步骤：

(1)根据胡萝卜的r、g、b三个分量图像，通过式②得到灰度图像gray，对灰度图像gray进行二值化，得到二值图像，二值图像的面积为a；

(2)在二值图像中提取出面积在某一范围(如[100020000])的连通区域；

(3)获得连通区域的等价椭圆，得到等价椭圆的长半轴a、短半轴b和周长c；

(4)根据椭圆面积公式s＝π*a*b计算出等价椭圆的面积；

(5)获得连通区域的圆形度e＝4*π*s/c²；

(6)计算虫眼比cr＝s/a；

当圆形度e大于阈值0.6时，判断其为虫眼，并计算出虫眼比。

附图说明

本发明有如下附图：

图1本发明的缺陷识别量化流程图。

图2本发明中须根识别量化流程图。

图3本发明中断裂识别量化示意图。

图4本发明中弯曲识别量化示意图。

图5本发明中开裂识别量化示意图。

图6本发明中须根识别示意图。

具体实施方式

以下结合附图对本发明作进一步详细说明。

实施例，请参照图1-图6，图1给出了本发明所述的胡萝卜缺陷识别量化主流程图；图2给出了本发明中胡萝卜须根识别量化检测的流程图；图3示出了本发明中断裂识别量化检测的示意图；图4示意性的给出了本发明中弯曲识别量化检测的原理图；图5展示了本发明中开裂识别量化检测的示意图；图6示出了本发明中胡萝卜须根识别量化检测的原理图。

首先采用ccd摄像头拍摄胡萝卜图像，图像分辨率为2592*1944且胡萝卜为水平放置。接着依据图1所示流程依次对所采集的图像进行如下处理：

对图像进行预处理，分别获取胡萝卜图像的h、s、v图像和r、g、b图像，然后对h分量图像采用阈值为0.2进行二值化，获得消除背景的二值图像；对s分量图像进行二值化，阈值为0.5，背景设为黑色；对v分量图像采用阈值为0.4进行二值化，获得二值图像；利用公式①得到灰度图像并二值化，阈值为100；利用公式②得到灰度图像，用全局阈值(100)进行二值化，得到突出红色部分的二值图像。接下来根据前面所提到的方法进行胡萝卜断裂、弯曲、开裂、须根和虫眼的量化识别。

胡萝卜断裂识别量化步骤如下：

1)对胡萝卜的h分量图像进行阈值分割，阈值为0.2，得到二值图像；

2)用面积去噪的方式去除噪声，获得二值图像的边缘；

3)根据二值图像的边缘，获取二值图像的最小外接矩形，得到最小外接矩形的长l和宽w；

4)计算胡萝卜的纵横比zr＝w/l；

5)提取出二值图像两端的边缘坐标，保存在q1和q2中(如图3(q)所示；

6)计算出q1和q2的斜率k1、k2(如图3(k1)所示)；

当纵横比zr大于0.25且k1或k2中有不少于5个无穷大时，即判断为断裂，本实施例中的纵横比zr为0.27且k1中有15个无穷大点，所以本实施例可以判定为断裂，且断裂面很陡峭。

胡萝卜弯曲识别量化实施步骤如下：

1)对胡萝卜的h分量图像进行阈值分割，阈值为0.2，得到二值图像；

2)对得到的二值图像进行面积去噪，去掉小的斑点；

3)用半径为3的圆形腐蚀单元对二值图像进行腐蚀处理，消除边缘毛刺；

4)在二值图像中提取出面积大于100000个像素面积的区域，即为胡萝卜主体区域；

5)抽取胡萝卜主体区域的骨架，获得骨架图像；

6)对骨架图像进行hough变换，获得如图4所示的两条拟合直线，两条拟合直线的夹角分别为θ1，θ2；

6)根据θ1，θ2计算得到弯曲度br＝180/π*|θ1-θ2|；

本实施例中，弯曲度br为11.24，弯曲度br越大，表示胡萝卜弯曲程度越严重。

胡萝卜开裂识别量化结果如图5所示，具体步骤如下：

1)获取胡萝卜rgb图像后，用公式①计算得到灰度图像，并二值化，得到二值图像，阈值为100；

2)在二值图像中提取出面积在[200050000]的连通区域；

3)获得连通区域的等价椭圆和连通区域的面积a2；

4)获得等价椭圆的长半轴a和短半轴b，椭圆面积为s＝π*a*b；

5)计算i＝a2/s、长半轴和短半轴之比a/b和开裂度kr＝2*180/π*arctan(b/a)；

本例中a/b＝6.53且i＝0.97，所以判断其为开裂，开裂度kr＝17.41。

胡萝卜须根识别量化步骤如下：

1)根据图2所示流程图，选取s分量图像进行二值化，得到二值图像，阈值为0.5；

2)对二值图像进行内部孔洞的填充；

3)提取出二值图像的边缘点坐标，保存到p中，如图6(p)所示；

4)按顺时针从p中提取出上边缘坐标(图6(p1))和下边缘坐标，分别存在p1和p2(图6(p2))；

5)对p1和p2每隔50个像素进行取样，取样后得到p3和p4，

6)分别对矩阵p3、矩阵p4求斜率，得到斜率k3、斜率k4，求取k3和k4的峰值，如图6(k)所示；

本例中有一个峰值为2.07，大于1.5，所以胡萝卜有一个须根，须根的位置在图上的坐标为(1417，377)；

胡萝卜虫眼识别量化步骤如下：

1)获取胡萝卜的rgb图像，通过公式②获得强调红色的灰度图像，并二值化，得到二值图像，阈值为100；二值图像的面积为a＝480936；

2)通过面积去噪的方式，在二值图像中获得面积在[100020000]的连通区域，即是虫眼区域；

3)获得虫眼区域的等价椭圆并得到等价椭圆的长半轴a＝34.16、短半轴b＝28.57、周长c＝201.87；

4)根据椭圆面积公式s＝π*a*b计算出等价椭圆的面积，即为虫眼面积；

5)圆形度计算公式e＝4*π*s/c²；

6)计算虫眼比cr＝s/a；

本例中椭圆面积s为3065.8个像素面积，虫眼比cr＝0.0064，虫眼圆形度为0.945。

本发明的优点：

本发明所公开的胡萝卜断裂、弯曲、开裂、侧根和虫眼等缺陷识别量化检测方法利用ccd相机获得胡萝卜实时图像，基于图像处理技术对获得的实时图像进行缺陷识别量化，克服了人工检测的主观性，提高了缺陷识别量化的客观性和准确性。应用于农业生产加工和分级，能够提高生产效率，降低生产成本。

本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。