基于机器视觉的大曲质量检测系统研究

2018-06-19张芯豪黄丹平2田建平2黄

食品与机械 2018年4期

张芯豪黄丹平,2田建平,2黄丹

(1. 四川理工学院机械工程学院, 四川自贡 643000； 2. 过程装备与控制工程四川省高校重点实验室，四川自贡 643000； 3. 四川理工学院生物工程学院，四川自贡 643000)

“曲乃酒之骨”，大曲质量直接影响白酒质量。现场要求曲胚几何参数标准长为(300±3) mm，宽为(200±3) mm，厚为(100±3) mm，标准曲坯应四角整齐、表面光滑、松紧一致、无缺边掉角、富有弹性。曲坯主要颜色为米黄色，其断面存在乳白色菌丝、红曲霉及其它污染霉菌[1]。目前大部分酒企仍依靠人工经验判断大曲质量，无量化指标，易受主观影响，判断结果不标准。同时大曲发酵时环境恶劣，高温度、高湿度、曲虫等因素导致酒企难以使用人工方法来完成该工作。因此酿酒行业迫切需要一种智能化曲房系统，对大曲生产过程进行智能化改造，从大曲的制作到大曲的发酵、出库等环节实现全面智能化生产，得到品质一致的高质量大曲。为实现智能化生产，检测大曲质量是最重要的一环，管理控制数据来源的选择也十分重要[2-4]。

影响大曲质量的因素有光滑度、水分、酸度、糖化力、液化力、穿衣、气味、菌丝、颜色等，可归纳为感官指标、理化指标和微生物指标。敖宗华等[5]研究了大曲感官指标、理化指标以及微生物指标，通过这些指标量化了大曲质量的判断标准。炊伟强等[6]研究了泸州老窖普级大曲和优级大曲各种指标的差别，进一步发现各指标与大曲质量关系。何宏魁等[7]将古井贡酒优级大曲和普级大曲进行比较分析，研究大曲感官特征与其内在微生物、理化指标间的关系。目前针对机器视觉技术领域的大曲质量自动在线检测技术研究，中国未见相关报道。

为解决白酒固态发酵自动化程度低和人工成本高的问题，同时为实现智能化曲房生产提供大曲质量数据来源，本研究通过机器视觉技术获取大曲视觉信息，动态阈值方法分割图像，测量大曲几何参数，通过Lab颜色空间提取大曲霉菌含量，通过CNN卷积神经网络检测大曲裂缝信息，建立视觉信息与大曲质量的对应关系表，综合判断大曲质量。

1 系统结构与工作原理

针对影响大曲质量的感官特征，研究其感官特征与视觉信息的对应关系，基于机器视觉、智能制造、运动控制等技术，研发一套大曲质量检测系统。如图1所示，大曲质量检测系统主要由面阵工业相机、环形LED光源、光电传感器、上位机、剔除机构、传输机构组成[8]。

1. 上位机 2. 传输机构 3. 大曲 4. 光电传感器 5. 环形LED光源 6. 面阵工业相机 7. 剔除机构

图1 大曲质量检测系统结构图

Figure 1 The structural diagram of Daqu quality inspection system

在机器视觉系统中工业相机、光源及打光方式的选择会直接影响视觉信息的提取效果，本系统工业相机选用德国Balser彩色面阵相机，型号为acA 1300-30gc，其分辨率为1 294×964 pixel，帧率为30 FPS。系统光源选用白色LED环形光源，其工作原理见图2。本系统中环形光源α为30°，可以有效消除大曲周围阴影，提高几何参数测量精度，同时可提高视觉信息中孔洞、裂缝、穿衣等细节信息[9]。该光源与面阵工业相机同轴安装，保证光源能均匀照射工业相机正下方的大曲。

该大曲质量检测系统工作原理：

(1) 传输机构2带动大曲3移动至光电传感器4正下方位置。

(2) 光电传感器4检测大曲3位置信息并传输给上位机1，该信息反馈至工业相机6作为工业相机6的触发信号。

α. 光源照射角度

(3) 工业相机6接收到触发信号，采集视觉信息并上传到上位机1中等待处理。

(4) 上位机1内软件对视觉信息进行处理，根据检测结果发送剔除信号给剔除机构7。

(5) 大曲3移动到剔除机构7区域，剔除机构7根据剔除信号剔除对应工位的不合格大曲3。

2 大曲质量检测算法

图3为系统软件流程图，主要针对大曲几何参数、颜色、裂缝3种特征进行检测。尺寸检测包含均值滤波、动态阈值分割、形态学处理等步骤；颜色检测首先将原图从RGB颜色空间转换到XYZ颜色空间，再由XYZ颜色空间转换到Lab颜色空间，最后通过特定阈值从不同通道中提取到特征目标；裂缝检测将原图灰度化后，裁剪ROI区域并分割为多张小图片，通过提前训练好的CNN神经网络判断该大曲是否存在裂缝。

图3 大曲质量检测流程图

2.1 大曲几何参数检测

在实际采集大曲视觉信息过程中，由于光照、电流、机械运动等原因，视觉信息中会生成脉冲噪声、散粒噪声和高斯噪声。同时现场阳光、工人走动、发酵房暗室等因素可能造成照度不均匀，导致所采集视觉信息中有阴影，各处对比度不同，从而影响大曲几何参数的检测。因此本研究采用动态阈值分割法对大曲图像进行处理，首先对图像进行均值滤波，均值滤波的基本原理是把数字图像或数字序列中一点的值用该点部分邻域中各点值的均值代替，从而削弱噪声影响[10]。图4和图5所示为一张大曲图滤波前和滤波后的直方图。

式(1)中η为动态阈值法的阈值范围，表示原图g与均值滤波后的图f在同一阈值下像素点个数之差，本研究中设为20。

(1)

图4 原图直方图

图5 滤波后直方图

式中:

λ——灰度值，范围为0～255；

gR——均值滤波前R通道图；

gG——均值滤波前G通道图；

gB——均值滤波前B通道图；

fR——均值滤波后R通道图；

fG——均值滤波后G通道图；

fB——均值滤波后B通道图；

η——阈值差。

(2)

再将0～255所有灰度值条件下的二值图相加，如式(3)所示。

(3)

实际采集时大曲周围可能存在曲渣，同时大曲外围可能存在突出的谷壳或麦秆。因此对分割后的二值图像I去除小面积连通域，再通过开运算(先腐蚀后膨胀)去掉毛刺，最后求出连通域最小外接矩形，寻找该矩形边缘，将边缘像素长度乘以物像比例，得到实际长度。

R=As/Ar，

(4)

式中：

R——连通域矩形度，表示一个连通域接近矩形的程度，越接近1越像矩形；

As——连通域面积，像素点个数；

Ar——连通域最小外接矩形面积，像素点个数。

曲胚俯视图应十分接近矩形，因此若R<0.95，即判断大曲不完整[11]。

2.2 大曲颜色检测

本系统中颜色识别主要针对大曲断面的乳白色菌丝、红曲霉及可能存在的其它颜色污染菌霉。工业相机采集到的大曲图像是RGB彩色图，由于天气、光源衰减、光照角度等原因可能导致采集到的大曲图像质量不一，同时总结现场检测人员的判断经验，采用一种基于人对颜色感觉的Lab色彩空间进行颜色检测。

Lab色彩空间有颜色均匀性的优点，但与XYZ色彩空间不同，未包含人类能够感觉到的所有颜色，所以要先将RGB色彩空间转换为XYZ色彩空间，见式(5)[12]。

(5)

再由XYZ色彩空间转换为Lab色彩空间，见式(6)。

(6)

经过试验，本系统最终识别乳白色菌丝和红曲霉的方法[13]：

(7)

(8)

如式(7)所示在L通道以阈值206分割图像为二值图，可以有效提出乳白色菌丝部分。同理如式(8)在b通道以阈值167分割图像为二值图，提取红曲霉。最后计算乳白色菌丝与断面面积百分比N1，红曲霉与断面面积百分比N2。

2.3 裂缝检测

大曲表面裂缝形态多样，裂开部分内部与表面结构、颜色一致，因此难以通过颜色识别、阈值分割、纹理分析等方法实现对大曲表面裂缝的检测。根据裂缝检测特点本研究选用一种基于人工监督的卷积神经网络(CNN)进行裂缝检测，经过试验可以达到要求。

CNN提取特征主要依靠交替的卷积层和降采样层来完成，然后通过BP算法来调整网络结构，最后多次迭代训练直到网络达到收敛状态来提高网络的泛化性能。从本质上来看，CNN是一种输入到输出的映射，通过大量学习所有训练样本输入与其对应输出的映射关系来判断测试样本类型，因此输入与输出没有确定的数学表达式。卷积神经网络在图像处理领域应用广泛，具有两大特征[14]：

(1) 相邻两层神经元之间的连接方式采用的是局部连接而不是全连接。

(2) 在同一特征图中，神经元(激活单元)的权值是共享的。通过局部连接和权值共享2种方式，CNN在很大程度上减少权值数量，降低网络复杂度，起到一定的正则化作用。

CNN识别准确率需要调节学习率、学习批大小、卷积核数和训练次数等，通过多次调整来确定参数[15]。在开始训练前，程序将输入图像进行规整，否则图像过大会导致训练系数过多而无法收敛。同时将所有的权值设置为较小的不同随机数值来保证训练网络可以正常地学习。学习率是运行SGD(随机梯度下降法)算法重要参数，值设置过大算法可能无法优化，过小会使得算法优化速度太慢且陷入局部最优。学习批大小需要能整除训练样本数量，否则会使部分样本不参与训练。训练次数需要通过对同一个样本进行训练，通过观察得到测试误差最小时对应的训练次数。

本研究算法中1张大曲图片分割出10张小图片作为裂缝样本进行检测，每一张大小为300像素×300像素，剩余部分有背景影响，根据样本图片含有裂缝的张数判断大曲整体图片是否含有裂缝。

2.4 综合判断

经过研究，本系统中视觉信息与大曲质量对应关系见表1。

表1 大曲质量等级指标†

†X表示1张完整大曲图片分割出的10张裂缝样本图片中判断为含有裂缝的图片张数；R为矩形度，表示大曲面积与其最小外接矩形面积的比；N1表示乳白色菌丝面积占断面面积的百分比；N2表示红曲霉面积占断面面积的百分比。

3 试验结果及分析

3.1 几何参数检测结果

现场采用半自动化压曲机制作大曲，因此在测量时大曲尺寸只可能比标准值小。图6为几何参数测量示意图。

图6(a)中大曲已发酵完闭2个月，因此其内部水分散失，尺寸缩小，同时其背景光照不均，有明显亮斑。图6(d)中显示大曲几何参数，其中长为292.04 mm，宽为170.1 mm，矩形度R为0.953 259，符合实际数据，可实现预期功能，用于测量大曲几何参数。

3.2 颜色检测结果

图7为大曲断面图，含有红曲霉和乳白色菌丝。红曲霉为污染霉，该大曲红曲霉与断面面积百分比为4.54%。乳白色菌丝表示微生物发酵成熟，该大曲乳白色菌丝与断面面积百分比为50.50%。

图6 几何参数测量示意图

图7 颜色检测示意图

3.3 裂缝检测结果

图8为部分裂缝样本图，采用CNN卷积神经网络检测大曲表面是否存在裂缝。本研究使用4 800个有裂缝样本和5 600个无裂缝样本进行训练，使用1 000个样本进行测试，其中有裂缝样本300个，无裂缝图700个。

本方法中有3个重要参数：学习批大小，表示一次参数更新运算所需的样本数量；学习率，表示参数到达最优值过程的速度快慢；学习次数，表示训练过程中数据迭代次数。将学习批大小设置为4，为寻找最佳学习率和学习次数，进行试验。由表2可知，学习次数为10时，最佳学习率为0.10。

图8 裂缝样本组合图

学习率有裂缝样本正确率/%无裂缝样本正确率/%总体样本正确率/%0.0139.093.877.40.0588.793.391.90.1092.096.395.00.1588.789.489.20.2058.788.786.70.5039.079.067.0

† 学习次数为10。

由表3可知，学习率为0.10时，最佳学习次数最低为30。针对本算法中裂缝样本特征，得到本研究最佳参数：学习批大小4，学习率0.10，学习次数30，该条件下大曲表面裂缝的识别正确率可达99.0%。

表3 学习次数对样本正确率的影响†

† 学习率为0.10。

4 结论

本试验研究了大曲质量检测系统通过动态阈值分割法、Lab颜色空间、CNN卷积神经网络等机器视觉技术处理采集到的大曲视觉信息，综合判断大曲质量信息。通过实验验证，本研究具有实时性好、非接触、易实现制曲工艺自动化等优点，对大曲几何参数检测精度可达到±1 mm，对红曲霉和乳白色菌丝可精准提取其数量，试验内裂缝样本识别正确率可达到99.0%，均满足现场制曲工艺要求，可实现对大曲质量检测并提高大曲质量检测的客观性和准确度。

本试验尚未研究大曲视觉信息和理化特性的联系，不能得到更全面的大曲参数，因此，后续将展开对此类问题的相关研究，以便实现制曲工艺的全面自动化。

[1] 泸州老窖集团有限责任公司. 泸型酒技艺大全[M]. 北京: 中国轻工业出版社, 2011: 82-85.

[2] 沈怡方. 白酒生产技术全书[M]. 北京: 中国轻工业出版社, 1998: 54-57.

[3] 李安军, 何宏魁. “基于物联网技术的数字化曲房综合管控系统”在浓香型白酒生产过程中的应用[J]. 酿酒科技, 2014(2): 63-64, 67.

[4] 敖宗华, 陕小虎, 沈才洪, 等. 国内主要大曲相关标准及研究进展[J]. 酿酒科技, 2010(2): 104-108.

[5] 肖冬光, 赵树欣, 陈叶福, 等. 白酒生产技术[M]. 2版. 北京: 化学工业出版社, 2013: 70-74.

[6] 何宏魁, 李兰, 汤有宏, 等. 古井贡酒大曲感官评价与其微生物、理化指标的关系研究[J]. 酿酒, 2017, 44(1): 58-63.

[7] 炊伟强, 敖宗华, 张春林, 等. 泸州老窖大曲感官特征与微生物、理化指标和生化性能的关联研究[J]. 食品与生物技术学报, 2011, 30(5): 761-766.

[8] 黄丹平, 廖世鹏, 于少东, 等. 基于机器视觉瓦楞纸板自动计数系统研究[J]. 机械工程学报, 2017, 53(10): 79-86.

[9] 郭康, 黄丹平, 王磊, 等. 基于机器视觉的曲面字符缺陷检测系统研究[J]. 现代制造工程, 2016(12): 122-127.

[10] 黄河, 李庆武, 范习健. 采用局部动态阈值的图像分割算法[J]. 光电子技术, 2011, 31(1): 10-13.

[11] RAFAEL C G, RICHARD E W, STEVEN L E. 数字图像处理[M]. 阮秋琦, 阮宇智, 译. 3版. 北京: 电子工业出版社, 2011: 487-490.