1. 引言
筒子纱作为纺织生产过程中重要的生产原材料,目前仍然采用人工巡检更换方式进行,这种方式用工成本高,不利于规模化生产。换筒机器人可以有效提高生产效率,开发换筒机器人的首要工作是设计纱筒定位算法,针对企业实际生产中纱筒架上的不同类型纱筒,本文开发了一种基于机器视觉的纱筒定位算法。
在纺织企业中,纱筒区域一般包含20到30排纱筒架,每排包括5组纱筒架,每个纱筒架按列分为A、B架,每个架子分为四层并且每层装有一个纱筒,即每排纱筒架装有40个纱筒如图1(a),每组纱筒架装有八个纱筒,且纱筒角度关于加强度对称,如图1(b)所示。这种纱筒安装方法导致在自动检测中出现了几个需要解决的问题:1) 同一个纱筒架上的每个纱筒相对于相机的空间位置不同,最终导致其成像结果存在较大差异;2) 受实际生产情况的限制,打光光源只能采用单一方向上的光源,但每个纱筒存在轴方向上的旋转,因此难以取得较好的打光效果;3) 获取目标纱筒架上的纱筒图像时,会同时获取到后排多排纱筒架以及生产设备,对纱筒图像分割产生较大干扰。
Figure 1. (a) Image of yarn bobbin; (b) Cone yarn bobbin; (c) Cylinder yarn bobbin; (d) Background yarn bobbin; (e) Background equipment
图1. (a) 纱筒图像;(b) 圆锥形纱筒;(c) 圆柱形纱筒;(d) 背景区域后排纱筒;(e) 背景区域中的生产设备
图像获取系统获得的纱筒图像如图1(a)所示,可以看到纱筒架上存在圆锥形纱筒(图1(b))和圆柱形纱筒(图1(c))两种类型纱筒,其中圆锥形纱筒内筒面积小,且其成像特征与后排纱筒内筒相近,在进行分割定位时背景区域中的纱筒(图1(d))会对算法带来较大的干扰,此外在部分纱筒图像中还会存在生产设备等背景,这种背景区域成像特点随机且灰度特征与纱筒内筒相似,因此在对纱筒进行分割定位前需要设计滤波器对图像进行降噪,减少背景区域对定位算法的干扰。
常见的降噪算法包括:中值滤波 [1] [2] [3] 、高斯滤波 [4] [5] 、非局域均值滤波 [6] 以及双边滤波 [7] 等。魏松涛和何文思 [8] 提出了一种结合梯度影响因子的均值滤波算法,将梯度信息加入到均值滤波的卷积核里,这样可以在滤波时尽量的保存区域边缘信息,对非边缘区域的像素进行较大程度的滤波。这样达到了减少图像整体噪声而且较好的保存了边缘信息的效果。刘保中和刘建宾 [9] 提出的改进型非局部均值图像降噪,通过采用噪声模型对图像整体灰度特性进行估计,再利用基于块匹配的滤波思想,能够在较大程度上保持各个像素块的信息,这对小型斑点噪声有着较好的效果,但是难以将各个像素块的特征信息强化并表达,且滤波后图像的整体复杂度并没有得到一定程度上的减小。Boughida和Boubekeur [10] 提出基于图像像素样本直方图与协方差矩阵的贝叶斯协同降噪算法,通过对图像的三通道灰度直方图与协方差矩阵的计算得到全局贝叶斯滤波器,该方法对图像像素样本体量有较大依赖,对不同体量的像素样本滤波效果相差较大;花思齐等 [11] 提出了一种基于改进型5 × 5平均中值滤波的背景抑制算法,通过构建5 × 5的滤波核并设置邻域滤波计算以达到在保护图像边缘与细节信息同时提高算法检测的实时性。
为了解决上述问题与挑战,本文提出了一种改进型各向异性扩散滤波结合卷积滤波的纱筒定位检测算法。主要工作包括:搭建图像采集系统并进行相机标定;然后设计滤波器对图像进行滤波,增强纱筒区域特征的同时抑制纱筒背景和纱筒纱线纹理等高频噪声;最后根据纱筒成像特征设计圆形卷积核对图形进行卷积,对卷积图像进行Blob分析,通过特征筛选提取纱筒内筒轮廓,实现对纱筒内筒的定位。
2. 材料与方法
2.1. 图像采集系统
纺织企业中,每排纱筒架中间存在1.2 M左右的间隔,这使得自动换筒系统可以自由移动至每个纱筒架的正面,负责纱筒图像获取的工业相机安装在换筒机器臂的上方,如图2(a)所示。系统相机选用海康威视MV-CA050-10GM搭配8 mm焦距镜头,得到分辨率为3072 × 2048的彩色纱筒图像如图2(b)所示。
(a) (b)
Figure 2. (a) Image acquisition system; (b) Image acquired
图2. (a) 图像获取系统;(b) 获取到的图像
2.2. 滤波器设计
2.2.1. 直方图均衡化
直方图均衡化属于间接对比度增强算法,利用图像直方图对纱筒图像的对比度进行调整以达到增强图像的局部对比度的同时不影响整体图像的对比度。设纱筒图像在(x,y)处的灰度值为f,均衡化后的图像为g,则对图像的增强过程可以表述为将(x,y)处的灰度f映射为g,对于一幅纱筒图像其灰度级fk出现的概率Pr为
(1)
其中,fk为第k级灰度;nk为图像中灰度级为fk的像素个数;n为像素总个数;L为图像的灰度级数。
直方图变换函数为:
(2)
在对纱筒图像进行直方图均衡化过程中,先对纱筒图像的灰度情况进行统计并计算出原始直方图分布,然后根据计算出的累计直方图分布求出fk与gk的灰度映射关系,按照映射关系对纱筒图像的各个像素进行灰度转换即可完成对纱筒图像的直方图均衡化。
2.2.2. 优化扩散方程的八方向各向异性扩散滤波
由于纱筒是由纱线缠绕而成的,纱筒的正面以及筒身存在明显的纹理信息。直方图均衡化使得纱筒图像整体对比度增强,但是也增强了纱线的纹理特征。为了在保留纱筒正面边缘信息的同时抑制纱线带来的纹理噪声,选择各向异性扩散滤波(Perona-Malik Anisotropic Diffusion Filter, P-M Diffusion Equation)可以在降低图像噪声的同时保留、增强边缘信息。
各向异性扩散滤波中首先要考虑扩散函数,扩散函数体现扩散方程对图像区域的扩散力度,对滤波平滑程度有着重大影响。PM扩散方程的理想滤波效果是平滑灰度梯度较小区域,保留图像边缘信息,因此扩散函数要对低灰度梯度区域扩散力度大,对高灰度梯度区域扩散力度小,常见扩散函数有式(3)、(4),本文在此基础上对式(3)进行改进得到新扩散函数如式(5),
(3)
(4)
(5)
其中k为梯度阈值常量,e为自然指数。三种扩散函数的曲线对比图如图3所示,可以看到新扩散函数相比于另外两个传统扩散函数收敛速度更快,对应扩散方程滤波时间更少、速度更快。
传统PM扩散方程是对像素点东南西北四个方向上的像素点进行散度计算,为了更好保留边缘细节信息,本文增加东北、东南、西北、西南四个方向上的散度计算,八方向散度计算示意图如图4所示。
Figure 4. Schematic diagram of the eight-way dispersion calculation
图4. 八方向散度计算示意图
图像I中像素点I(x,y)的八方向散度计算公式如下:
(6)
对于图像I在离散条件下的迭代模型如下:
(7)
2.2.3. 基于圆形卷积滤波的纱筒内筒定位算法
纱筒图像背景复杂,直接通过阈值分割、Blob分析等方法难以准确定位纱筒内筒位置,考虑纱筒内外筒在几何特征上呈现出圆形特征,设计圆形卷积滤波核,对图像进行卷积操作,设计圆形滤波核如图5:
在设计卷积核时要考虑实际环境中打光效果给图像带来的影响,本文处理对象中由于打光侧位于纱筒左侧,因此在设计时要对纱筒左上侧光照较强的区域进行抑制。卷积后图像保留下的信息主要是前后排纱筒的内筒区域,通过快速阈值分割得到含有目标纱筒内筒信息的ROI区域,在对ROI区域进行Blob分析,先进行连通域分析再根据区域面积进行筛选,得到纱筒内筒空心区域。
为了分割提取纱筒内筒区域,对卷积滤波前后图像执行图像相减操作,得到的图像中背景及纹理信息被有效抑制,接着通过分割以及曲线拟合得到纱筒内筒外轮廓,具体方法如下:
1) 根据纱筒内筒空心区域位置信息建立圆形掩膜,分割出包含纱筒内筒的区域。
2) 为减少内筒空心区域给分割算法带来的干扰,内筒区域需要去除内筒空心区域,得到包含纱筒内筒和纱线信息的图像I。
3) 采用梯度、灰度阈值相结合的分割方法对图像I进行分割,得到纱筒内筒大部分区域,计算区域外轮廓,并建立二元二次模型,选择边缘曲率变换平稳点作为拟合点,通过最小二乘法拟合得到纱筒内筒外轮廓。
3. 实验
3.1. 相机标定
根据仿真标定的分析结果认为网格尺寸在6~14 mm时标定误差较小 [12] ,以MATLAB标定工具箱为基础进行试验论证,对标定格尺寸在2~30 mm范围内的标定板进行实验,证明棋盘格尺寸在10~15 mm时标定效果最好 [13] 。同时,文献也对标定图片数量与标定效果关系进行大量实验,实验证明当图片数量达到8张时标定结果接近理想值,当标定图片继续增多后,结果趋于平稳。综上,本文选择尺寸15 mm、10 × 10的标定板,选取12个不同的角度以张正友标定法进行相机标定如图6,图7。
相机标定时的相对位置以及每幅标定图像的标定误差如图8所示,标定得到相机的内参矩阵和畸变矩阵D1如表1所示。
Figure 6. Camera calibration. (a) Multi-angle calibration images; (b) calibration image corners
图6. 相机标定。(a) 多角度标定图;(b) 标定角点图
Figure 7. Camera calibration result. (a) Camera calibration error; (b) Camera calibration angle
图7. 相机标定结果。(a) 相机标定误差;(b) 相机标定角度
3.2. 算法实验结果
3.2.1. 直方图均衡化效果
根据上述直方图均衡化模型对纱筒三通道图像分别进行滤波,并将滤波后的三个通道再次叠加得到直方图均衡化后的彩色纱筒图像,滤波效果及其灰度图像灰度直方图如图8所示。可以看到经过均衡化滤波后的图像灰度特征得到加强,纱筒内筒与纱线区域灰度特征分化明显,特别是高光纱筒滤波效果更突出。
(a) Gray Image (b) Histogram equalisation effect(c) Grayscale histogram of the original image (d) Grayscale histogram after equalisation
Figure 8. (a) Gray image; (b) Histogram equalisation effect; (c) Grayscale histogram of the yarn bobbin original image; (d) Grayscale histogram after equalisation
图8. (a) 灰度图像;(b) 直方图均衡化;(c) 纱筒原图灰度直方图;(d) 均衡化后灰度直方图
3.2.2. 扩散滤波效果
对图像采用优化后的扩散方程模型对图像进行滤波并与传统扩散模型进行对比如图9所示,可以看到,相比于传统的扩散方程滤波,优化后的滤波模型对纹理、纱线印痕等高频噪声抑制效果更好,对图像整体噪声抑制效果更好,同时保留图像的边缘信息。
3.2.3. 高斯卷积滤波效果
对滤波后的图像应用本文设计的圆形卷积核进行卷积滤波效果如图10所示,可以看到图像中的大部分纱线区域以及背景纱筒区域被有效抑制,特别是后排纱筒架等灰度特征与纱筒内筒相近的噪声也得到了有效的抑制,大幅降低了图像的噪声强度,减少背景区域对全局分割算法的干扰。
(a) Original image (b) Convolution filtering effect
Figure 10. Filter result of yarn bobbin. (a) Original image; (b) Convolution filtering effect
图10. 纱筒滤波效果。(a) 原图像;(b) 卷积滤波效果
卷积后即可对纱筒图像进行分割得到纱筒内筒区域,首先进行快速阈值分割,在进行连通域分析,通过区域面积、圆度等区域特征进行筛选即可得到纱筒内筒区域,考虑到纱筒轮廓呈现椭圆特征,建立二元二次拟合模型对轮廓进行拟合,并计算拟合曲线的最小外接矩形,完成对纱筒的定位,实验效果如图11所示。对卷积滤波后的纱筒图像按照上述图像运算得到目标纱筒内筒区域提取效果如图11(a)所示,其中可以看到部分纱筒区域结果中包含部分干扰区域,对其计算最小外接矩形并标注如图11(b)所示。定位结果与误差数据如表2。
(a) Region Extraction (b) Region Positioning Frame
Figure 11. Contour extraction and positioning result of yarn bobbin. (a) Yarn bobbin area extraction. (b) Yarn bobbin positioning frame.
图11. 纱筒轮廓提取及定位效果。(a) 纱筒区域提取;(b) 纱筒定位框
Table 2. Yarn bobbins positioning test results
表2. 纱筒定位检测结果
4. 结束语
本文针对企业生产场景下获取到的纱筒图像,设计了一种基于卷积滤波的纱筒定位检测算法。首先对图像进行滤波,通过对图像直方图均衡化增强图像区域特征,然后对传统各向异性扩散滤波进行优化,提出一种更符合企业使用要求的扩散函数并增加四个方向上的散度计算,增强滤波模型对纹理等高频信息的抑制效果,同时保留图像区域边缘信息,接着根据纱筒内筒几何特征与打光成像效果设计圆形卷积核对图像进行卷积,抑制纱筒图像大部分背景,最后对卷积图像分割,并进行Blob分析,根据圆度、区域面积等特征筛选得到纱筒内筒区域,再建立二元二次拟合模型拟合纱筒内同外轮廓,计算外轮廓最小外接矩形,实现对纱筒的定位检测。实验证明本文方法在保证检测效率的同时具有较高的准确率,能够对一个纱筒架上的八个不同类型纱筒同时定位检测。该论文对在实际生产情况下的纱筒纱线余量检测具有一定的启示意义,为传统算法在复杂条件下的目标检测提供支持。