芦笋基部最大直径检测

摘    要：为了解决当前芦笋尺寸测量效率低、精度差以及测量过程中易破损等问题，针对芦笋基部最大直径检测，提出基于视觉检测的外包裹矩形法，进而实现芦笋的自动分级。首先，运用传统图像处理方法（图像增强、分割和边界轮廓提取）得到芦笋轮廓的二值图像；随后，将芦笋基部外弯曲部分图像裁剪；最后，对裁剪后图像矩阵中的像素点进行搜索，确定四方位置上的首个像素点，并基于此绘制出芦笋基部的外包裹矩形，其宽度即为芦笋基部最大直径。试验结果表明，检测尺寸误差范围在0.01～0.50 mm，符合国家芦笋等级规格中的允许误差范围，以此方法实现芦笋分级的准确率达100%（依照国家芦笋等级规格）。此方法不仅高效准确且无损伤地完成芦笋的尺寸检测，也为其他农副产品以及不规则物体的尺寸检测提供了新的思路。

关键词：芦笋；分级；尺寸检测；图像处理；Matlab

中图分类号：S644.6 文献标志码：A 文章编号：1673-2871（2023）08-105-07

Abstract： In order to solve the current problems of low efficiency and poor accuracy of asparagus size measurement as well as the ease of breakage during the measurement process， and thus to achieve automatic grading of asparagus. In this paper， an outer wrapping rectangle method based on visual inspection is proposed for the detection of the maximum diameter at the base of asparagus. First， a binary image of the asparagus outline is obtained using traditional image processing methods（Image enhancement， segmentation and boundary contour extraction）; subsequently， the image of the outer curved part of the asparagus base is cropped; finally， the pixel points in the cropped image matrix are searched to determine the first pixel point in the quadrilateral position， and based on this， the outer wrapping rectangle of the asparagus base is drawn， the width of which is the maximum diameter of the asparagus base. The results of the experiment show that the detection error range is between 0.01 and 0.5 mm， which is in accordance with the national asparagus grading specifications， and that this method achieves 100% accuracy in grading asparagus （In accordance with the national asparagus grading specifications）. This method not only allows for efficient， accurate and non-destructive dimensional inspection of asparagus， but also brings new ideas to the dimensional inspection of other agricultural products and irregular objects.

Key words： Asparagus; Grading;  Size detection; Image processing; Matlab

芦笋又名石刁柏，质地柔软、口感细腻、风味鲜美，富含多种氨基酸、蛋白质和维生素，具有一定的药用价值，目前我国芦笋的种植面积已达到10万hm2，高居世界第一，除了西藏和青海外，全国各个地方均有种植[1]。

芦笋在大批量收割后，其尺寸参差不齐，如果不按照某一指标进行检测分级，种植者无法获得最大的经济效益。当前消费者和食品加工单位都对芦笋外部品质的要求不断提高，芦笋的外部尺寸检测分级显得举足轻重。农业农村部印发的《“十四五”全国农业机械化发展规划》指出，围绕果蔬等鲜活农产品保质增值，研发分级等初加工机械[2]。经过芦笋尺寸的检测分级，使得不同尺寸的芦笋物尽其用，像鲜芦笋、芦笋罐头、芦笋茶和芦笋美容品等[3]。综上所述，为实现芦笋的自动分级，对芦笋尺寸的检测至关重要。芦笋等级规格有基于长度和基部最大直径2种，芦笋的机械化收割主要是针对长度在18～23 cm的嫩茎[4]。笔者的研究主要通过视觉手段对芦笋（本试验针对绿芦笋）基部最大直径进行检测。

现阶段对于芦笋等农副产品尺寸的检测分级主要是通过人工或者传统机械（辊轮、筛网等装置）来实现的[5-11]。人工检测分级的主观性太强，导致精准程度较低；整个检测过程需要大量的劳动力，人工成本过高，造成经济效益低；同时人工检测还会造成过多的接触，会产生一定的食品安全隐患。传统机械检测主要是利用间距不同的辊轮、直径不同的圆孔以及滚筒的离心力进行分级，优点是可以通过上述参数的调整来实现对大多数的果蔬产品进行分级，适用度较高，并且分级效率高，但缺点也是显而易见的，一是分级过程中容易造成芦笋的破损和折断，破损的芦笋也会对其他芦笋造成影响；二是分级的精度不高；三是仍需要人工的辅助作业，自动化程度较低。

视觉检测主要是利用计算机来模拟人类的视觉功能，对于从外界采集到的图像进行分析和理解，来获取想要的信息以实现目标物尺寸的测量和分级等目的，具有精度高、检速快、自动化和实时性等优点[12]。国内外的众多学者很早就将该技术应用到各种农副产品的尺寸检测上，包括马铃薯、苹果、柑橘、火龙果和黄瓜等[13-21]，但上述的检测目标物大多是球形或者类球形的农产品，针对芦笋这样细长且具有一定柔软性目标物的研究还很少。独具一格的外形造成芦笋尺寸检测具有较高的难度，主要有以下几个方面：芦笋的基部形状并非圆形，难以通过常规方法对其进行测量；芦笋会出现部分的弯曲且具有一定的柔软性，增大测量的难度；由于收割等造成芦笋基部的不完整，其最大直径难以确定。

综上所述，笔者提出通过视觉检测的方式来对芦笋基部最大直径进行检测，首先对芦笋图像进行采集，经过图像预处理、图像分割和边界轮廓提取，提出一种新型方法从芦笋边缘二值图像上获取芦笋基部最大直径，即从图像四个方向上寻找其最外侧边缘点，以此绘制芦笋基部的外包裹矩形，其宽度即为芦笋基部的最大直径。

1 材料与方法

1.1 检测要求及试验环境

1.1.1 检测要求 芦笋外部形态和基部轮廓示意图如图1所示，本试验的检测对象为芦笋基部最大直径。

芦笋的尺寸等级标准按照中华人民共和国农业农村部2008年5月18日发布的《芦笋等级规格》，以芦笋基部最大直径作为划分等级规格的指标[22]，如表1所示。

1.1.2 试验环境 笔者试验所用的所有芦笋样本均来自于山东省菏泽芦笋种植基地，品种为绿芦笋。试验时间为2022年6－8月，试验地点为宁夏大学工程训练中心。所用试验环境：光源采用自然光源，电脑处理为Intel（R） Core（TM） i7-8550U CPU，操作系统为Windows 11家庭中文版，图像处理软件采用Matalb2020b。

1.2 图像采集

笔者试验的图像采集平台主要由相机、镜头以及夹持相机的移动导轨组成，固定高度和位置，控制该变量对成像的影响，考虑到未来实际应用中的光照情况，光源为自然光源。采集到的芦笋图像分辨率为600×800。

在对芦笋的图像采集之前，先使用标准件（笔者的研究使用第五套人民币一元硬币，直径为25 mm）对该试验采集装置进行标定，获取实际尺寸与像素尺寸的比值K。使用上述的采集平台对硬币的图像进行采集，通过一系列图像处理操作，测量硬币的面积，得到硬币的像素直径为64.701 3 pt，最后计算K值如式（1）所示。

1.3 芦笋图像处理

芦笋尺寸检测程序主要包括图像增强、灰度化、图像分割、芦笋边界轮廓提取、尺寸测量，程序基于Matlab软件进行编写。

1.3.1 图像增强 由于芦笋图像是在自然光照条件下进行采集的，因此存在光照不均匀的情况，图像会产生部分阴影或者模糊，如图2所示。使用J=imflatfield（I，Sigma）进行去模糊，其中I为原始图像；Sigma为去模糊所使用高斯平滑的标准差，Sigma值经过多次试验，最优值30；J为去模糊后返回的图像，如图3所示。

1.3.2 灰度化 由于检测的对象为芦笋，对采集到的芦笋图像（RGB图像）进行灰度化处理，可减少后续图像处理过程的计算量。笔者通过计算芦笋图像中每个像素点的R、G、B三个分量的加权和，将RGB彩色图像转化为灰度图像，采取权重平均法，如式（2）所示。

灰度化后的图像如图4所示。

1.3.3 图像分割 图像分割是将图像分离成前景（目标对象）和背景的过程，是图像处理过程中的基础性问题，分割的优劣直接关乎到后续的尺寸测量。在众多算法中，笔者选择了主动轮廓法[23-24]，又称snake法，是一种迭代式区域增长图像分割算法，snake表示平面中一组有序的像素点[xi=xiyiT]，其构成了感兴趣区域周围的一个边界曲线。该边界曲线沿着能量减小的方向进行移动，如式（3）：

其中，EI—内部能量，表示假定的边界曲线；EE—外部能量，表示当前snake所在点的图像。所以说，主动轮廓法的每一次迭代都要沿着能量降低的方向，即[Ek+1<Ek]，然后用新点代替旧点，继续迭代，直到E最小。考虑芦笋的大小形状以及反复试验，笔者选择80×480大小的矩形作为初始snake，迭代600次时目标和背景分割效果最佳。初始snake和分割后的芦笋图像如图5和图6所示。

使用其他图像分割算法（Otsu法阈值分割和自适应阈值分割）对芦笋图像进行分割，分割后的结果如图7和图8所示。从芦笋图像的分割结果可以看出，尽管Otsu法阈值分割相对于自适应阈值分割保留了更多的信息，但是两种阈值方式都没有完整的将芦笋分割，同时从分割结果来看，分割后的图像上还有少许的噪声点（图像上白色区域内的黑色竖线和黑点），不利于后续图像的边界提取以及尺寸测量。综合分析，使用主动轮廓法来进行芦笋图像的分割可以更好地完成芦笋尺寸的测量。

1.3.4 芦笋边界轮廓提取 分割后的芦笋图像为二值图像，已无噪声点，对其进行边界轮廓提取时选取低噪声情况下提取效果较好的Roberts算子[25]，Roberts算子采用对角线方向相邻两像素之差近似梯度幅值，检测图像中的边缘（灰度函数发生急剧突变的像素点之间的连线），其模板如式（4）所示。提取后芦笋轮廓如图9所示。

1.4 芦笋尺寸测量

尺寸测量即从图10的芦笋边界轮廓图像中获得预定的尺寸信息（芦笋基部最大直径），从二维平面上观察，即测量芦笋基部边缘点所构成的最大水平线段的长度。图像处理中测量线段的常用方法之一就是寻找角点，然后再计算两个角点之间的距离，Harris角点检测算法被公认为是较好的寻找角点的方法[26-27]。除此以外，还可以利用曲线拟合思想[28]，将芦笋边缘像素点近似拟合成为一条或者数条相连的曲线，通过计算曲线的长度或者曲线之间的距离来近似获得芦笋的基部最大直径。但是，上述两种方法不适用于芦笋基部最大直径的获取，原因主要是芦笋形状的特殊性，从芦笋的边缘轮廓图可以看出来，其基部位置多为曲线，且与两个边缘以曲线过渡，以及收割等原因导致芦笋基部不完整，需要将基部上完整部分的直径作为该芦笋的基部最大直径，从而造成曲线的难拟合或者角点的错误检测，两种方法所得结果如图10所示，从图上可以看到，由于芦笋基部的复杂性造成角点的错误检测，以及基部因不完整导致曲线拟合，难以计算正确距离。