基于图像处理的烟丝形貌特征的测量方法与流程

1.本发明涉及物理实验与烟草机械设备技术领域，尤其涉及基于图像处理算法的烟丝形貌特征测量方法。


背景技术：

2.烟丝的形貌特征，包括圆形度、矩形度、球形度，是卷烟工艺中一项重要的评价指标，对卷烟的感官质量及物理指标有明显的影响。在烟丝流化的研究中，不同的形貌特征影响烟丝的曳力。烟丝形状越偏离球形颗粒，表面受到的流化阻力越大。另一方面，烟叶经过片烟处理设备与切丝设备后，包括分切设备、叶片松散回溯机与切丝机，应该对烟丝的形貌进行分析。叶梗分离设备、供料成条机、卷接成型机中也应该采用形貌特征分析方法来判断国产烟草机械各机组的功能可靠性。传统测量方式多为人工评价测量，劳动强度大，采样周期长，评价客观性差，与烟草企业的现代化设计理念不相符。多种形态参数未形成统一的实验规范，因此，应该提出一种合适的烟丝图像测量方法，对混合烟丝形貌特征进行分析与评价。
3.在设计烟丝颗粒流化反应器之前，一般采用理论与数值模拟相结合的形式，初步分析烟丝的流化特征。受计算方法与数值模拟软件的限制，烟丝颗粒一般采用球形颗粒替代，加以形状因子来表征烟丝与规则形状的区别。因此，应该提出一种等效直径测量方法，计算混合烟丝的等效直径均值与偏差。
4.现有的烟丝物理性能测量效率低且精度不足，未能对烟丝矩形度、圆形度、球形度、等效直径等进行可视化表示。以往的测量方法均采用人工检测，依靠目测判断已不能满足烟丝标准规格测量要求。目前，对烟丝宽度的测量较为先进的方法为ccd图像测量，但仅对烟丝宽度进行了分析，尚未提出针对烟丝的矩形度、圆形度等特征属性的测量方法。


技术实现要素：

5.本发明的目的在于提出一种基于图像处理的烟丝形貌特征的测量方法，可以快速获得烟丝圆形度、矩形度、球形度和等效直径等形貌特征参数，可以快速测量大量随机混合的烟丝，具有较高的测量精度，并能根据烟草机械检测需求，扩展程序功能。
6.为了实现本发明的目的，本发明采取的技术方案为：
7.基于图像处理的烟丝形貌特征的测量方法，其特征在于：包括以下步骤：
8.步骤一.获取烟丝图像；所述烟丝图像中的烟丝相互不重合；
9.步骤二.对烟丝图像进行图像预处理分割得到各烟丝区域；
10.步骤三.统计各烟丝区域的实际周长和实际面积作为烟丝的实际周长和实际面积s，测量得到烟丝平均厚度h，以烟丝平均厚度h*实际面积s得到烟丝区域的体积v
p
；
11.步骤四.得到各烟丝区域的最小外接矩形，以最小外接矩形的对角线长度作为烟丝长度；通过各烟丝区域的周长和面积得到圆形度ψ；通过烟丝的平均厚度 h和烟丝的实际面积得到烟丝的等效直径d
p
；通过烟丝的平均厚度h、烟丝的实际面积和实际周长得到球形
度。
12.进一步的改进，步骤一中，获取烟丝图像的方法如下：
13.设置数码相机、led光源和试验台，所述数码相机电连接图像采集计算机，图像采集计算机电连接显示器；所述试验台位于数码相机正下方，试验台上方一侧固定有水平尺，另一侧固定有白板；所述led光源周向环绕设置在数码相机两侧；
14.数码相机将拍摄得到的烟丝图像上传至图像采集计算机。
15.进一步的改进，所述步骤三中，烟丝的实际周长和实际面积获取方法如下：
16.步骤a.白板上放置有一个已知长度的标定物l0，作为单位标定依据的标定物，拍摄图像p0；
17.步骤b.按照标准gb/t16447
‑
1996改变烟丝的含水率，随机打乱烟丝测量样品，取出部分烟丝作为检测对象，将所述烟丝平铺于白板，随机分散烟丝，防止烟丝之间的重叠；
18.步骤c.将玻璃板缓慢平压在所述烟丝与白板上，使所述烟丝平铺于水平面，防止烟丝的卷曲；
19.步骤d.在玻璃板平压所述烟丝3min后，拍摄静态图像p，重复步骤b到步骤c，拍摄n张烟丝图像；
20.步骤e.读取拍摄的图像p0，得到标定物的垂直投影像素长度l
p
；
21.步骤f.根据得到的像素长度l
p
，得到长度尺寸比lca：
[0022][0023]
步骤g.人工测量得到烟丝的平均厚度h；烟丝区域的像素周长*lca即烟丝区域的实际周长，烟丝区域的像素面积*lca2即烟丝区域的实际面积。
[0024]
进一步的改进，所述步骤四中，通过最小外接矩形的对角线长度*lca得到各区域的烟丝长度，圆形度ψ、等效直径d
p
和球形度的计算公式如下：
[0025]
圆形度ψ：
[0026][0027]
等效直径d
p
：
[0028][0029]
球形度
[0030]
s
f
＝2s+hl
[0031][0032]
矩形度r：
[0033][0034]
其中，h为烟丝的平均厚度、d为烟丝投影面积的等效直径、v
p
为烟丝区域的体积,s
f
为烟丝实际表面积，l为边缘轮廓提取后的烟丝实际周长，s为烟丝的实际面积，sr为烟丝最小外接矩形的实际面积。
[0035]
所述烟丝形貌特征综合性能评价指标pjzb计算公式如下：
[0036][0037]
其中，r
s
为烟丝理想的矩形度，d
ps
为烟丝理想的等效直径，ψ
s
为烟丝理想的圆形度，φ
s
为烟丝理想的球形度，ω1、ω2、ω3、ω4均为权重。
[0038]
进一步的改进，ω1为0.7，ω2为0.15,ω3为0.1,ω4为0.05；综合性能评价指标pjzb越接近于0，烟丝质量越好。
[0039]
与现有技术相比，本发明的优点在于：
[0040]
(1)提出了一种基于图像处理的烟丝形貌特征与等效直径测量方法，设计了一款未有的烟丝形貌特征测量系统。
[0041]
(2)基于图像处理算法，能够快速无损测量大量分散烟丝的几何形态属性。形态参数主要包括：，圆形度、矩形度、球形度等。
[0042]
(3)基于图像处理算法，能够统计实验区域烟丝等效直径。
附图说明
[0043]
图1为基于图像处理的烟丝形貌特征测量系统；
[0044]
图2为烟丝形貌特征测量技术方案；
[0045]
图3为烟丝图像处理测量系统；
[0046]
图4
‑
a为线段原始图像；
[0047]
图4
‑
b为图像处理后线段最小外接矩形；
[0048]
图5为以不同线段为标准计算线段长度的相对误差；
[0049]
图6为以不同线段为标准计算线段长度相对误差的均值及其均方差；
[0050]
图7
‑
a为烟丝原图；
[0051]
图7
‑
b为烟丝灰度化图；
[0052]
图7
‑
c为烟丝滤波去噪图；
[0053]
图7
‑
d为烟丝二值化图；
[0054]
图7
‑
e为烟丝取
‑
反、忽略小面积图；
[0055]
图7
‑
f为烟丝区域分割标记图；
[0056]
图8为球形度分布直方图；
[0057]
图9为烟丝长度分布直方图；
[0058]
图10为当量直径分布直方图；
[0059]
图11
‑
a为圆形度频率分布图；(横坐标表示圆形度，纵坐标表示频率；圆形度可描述为某一图形接近圆形的程度，频率是每个对象出现的次数与总次数的比值)；
[0060]
图11
‑
b为矩形度频率分布图；(横坐标表示矩形度，纵坐标表示频率；矩形度体现物体对其外接矩形的充满程度，频率是每个对象出现的次数与总次数的比值)；
[0061]
图11
‑
c为球形度频率分布图；(横坐标表示球形度，纵坐标表示频率；球形度是指物体的形状接近球体的程度，频率是每个对象出现的次数与总次数的比值)；
[0062]
图11
‑
d为圆形度概率密度分布拟合图；(横坐标表示圆形度，纵坐标表示概率密度；圆形度可描述为某一图形接近圆形的程度，概率密度表示瞬时幅值落在某指定范围内的概率)；
[0063]
图11
‑
e为矩形度概率密度分布拟合图；(横坐标表示矩形度，纵坐标表示概率密度；矩形度体现物体对其外接矩形的充满程度，概率密度表示瞬时幅值落在某指定范围内的概率)；
[0064]
图11
‑
f为球形度概率密度分布拟合图；(横坐标表示球形度，纵坐标表示概率密度；球形度是指物体的形状接近球体的程度，概率密度表示瞬时幅值落在某指定范围内的概率)；
[0065]
其中，显示器1、图像采集计算机2、led光源3、白板4、玻璃板5、数码相机6、试验台7、水平尺8、烟丝样本9。
具体实施方式
[0066]
为更进一步阐述本发明为达成预定目的所采用的技术手段及功效，以下结合附图及较佳实例，对依据本发明提出的基于图像处理的烟丝形貌特征的测量方法，详细说明如下：
[0067]
图1为基于图像处理的烟丝形貌特征与等效直径测量系统。提供了一种针对于烟草机械生产的成品烟丝检测方法，包括对烟丝矩形度、圆形度、球形度、等效直径等特征的可视化表达。整个实验装置包括显示器1、图像采集计算机2、 led光源3、白板4、玻璃板5、数码相机6、试验台7、水平尺8、烟丝样本9。
[0068]
本发明解决其技术问题所采用的技术方案依次为搭建物理系统、图片获取、单位标定、输入烟丝的厚度与像素尺寸比、图像预处理、烟丝多连通区域分割、烟丝长度与体积输出、烟丝圆形度、球形度输出、烟丝等效直径、球形度分布图计算等。具体实施步骤如下：
[0069]
步骤一：根据图1方案搭建烟丝形貌特征与等效直径测量实验系统，连接数码相机与计算机测量系统。两盏led光源左右对称放置，使烟丝样本照明均匀性好，提高了成像质量。光源下方放置白板，烟丝与背景完全分离，增加了烟丝之间的色差，提高了烟丝在ccd中的成像质量。图像采集计算机能够对采集的多组烟丝样本图像进行分析，基于自主创新设计的烟丝形貌特征与等效直径测量系统完成了烟丝等效直径、球形度等形状属性测量。
[0070]
步骤二：固定ccd相机与白板，并且在白板上放置一个已知面积或者长度的的标定物l0，作为单位标定依据，拍摄图像p0。
[0071]
步骤三：选取芙蓉王的成品卷烟采用的烟丝作为测试样品，按照标准 gb/t16447
‑
1996改变烟丝的含水率，然后随机打乱测量样品，取出部分烟丝作为检测对象。将一定数量的烟丝平铺于白板，随机分散烟丝，防止烟丝之间的重叠。
[0072]
步骤四：将玻璃板缓慢平压在烟丝样本与白板上，使烟丝样本平铺于水平面，防止烟丝的卷曲。
[0073]
步骤五：在玻璃板平压烟丝样本3min后，拍摄静态图像p。重复步骤三到步骤四，拍摄n张烟丝图像。
[0074]
步骤六：通过图像采集计算机读取数码相机拍摄的图像并保存；
[0075]
步骤七：在计算机单位标定系统中，读取拍摄的标定图片p0，得到物体的垂直投影像素长度l
p
。
[0076]
根据得到的像素长度l
p
，计算长度尺寸比lca：
[0077][0078]
步骤八：通过烟丝形貌特征测量分析系统读取图像p，记录尺寸比lca与平均厚度h。
[0079]
步骤九：基于烟丝形貌特征测量分析系统，统计出各个烟丝的等效直径、球形度、面积与对角线尺寸等数据。图像处理技术包括图像预处理、多连通域分割算法等，其中图像预处理包括图像灰度化、滤波去噪、图像二值化、背景去除、形态学闭运算等。
[0080]
基于图像处理的烟丝形貌特征测量系统，如图3所示。对烟丝形貌特征进行无损测量，避免了人为因素的干扰，可快速获取相关的数据，可靠性高。具体实施流程如下：
[0081]
①
运行烟丝图像处理测量系统，操控区内读取图片，对图像进行预处理。
[0082]
②
输入区内设置像素比为0.33，烟丝平均厚度为0.2mm。
[0083]
③
图片灰度化，将rgb图像转化为灰度化图像，减少图像内存；
[0084]
④
采用中值滤波技术去除烟丝灰度图中局部噪点；
[0085]
⑤
将灰度图转换为二值图，并进行背景去除。此时拍摄的烟丝为黑色，背景为白色；
[0086]
⑥
对图像进行二值图取反处理，使烟丝显示为白色，背景为黑色，此时的图像方便于后续形态学操作；
[0087]
⑦
通过形态学闭运算修正烟丝轮廓与中间孔缝；
[0088]
⑧
通过多连通域分割划分烟丝区域，并通过区域伪彩方式显示；
[0089]
⑨
提取各个区域的信息，由尺寸比lca完成像素转换，统计烟丝实际周长l，面积s；
[0090]
⑩
0计算各连通域中烟丝的对角线长度、圆形度ψ、矩形度r、等效直径d
p
、球形度φ等，计算公式如下：
[0091][0092]
s
f
＝2s+hl
[0093][0094][0095]
[0096]
式中，h为烟丝的平均厚度、d为烟丝投影面积的等效直径、vp为烟丝的体积,l为边缘轮廓提取后的烟丝周长、sr为烟丝最小外接矩形的实际面积。
[0097]
由于圆形度、矩形度、等效直径与球形度4个形貌参数的功能不一，除了等效直径为尺寸参数，其他三个都是无因次参数，不能直观的反应烟丝形貌的综合质量。为量化烟丝质量参数，提出综合性能评价指标pjzb计算式：
[0098][0099]
式中，ω
i
为不同指标的对应权重，r
s
为烟丝理想的矩形度，d
ps
为烟丝理想的等效直径，ψ
s
为烟丝理想的圆形度，φ
s
为烟丝理想的球形度。烟丝理想的矩形度、等效直径、圆形度和球形度应该根据烟草公司想要的烟丝尺寸来计算。由于矩形度对烟丝形貌质量参数影响相对较大，权重为ω1＝0.7。圆形度、相对等效直径与球形度影响较小，则权重依次为ω2＝0.15,ω3＝0.1,ω4＝0.05。
[0100]
综合性能评价指标pjzb越接近于0，表示烟厂加工出的实际烟丝形貌质量与理想烟丝形貌质量接近。
[0101]
其实施方式如下：
[0102]
1.单位标定
[0103]
在对烟丝进行实际测量之前，先对系统进行了标定试验，在同一测量条件下，对多个标准线条的图像像素长度进行测量，进行了像素尺寸比的验证。读取拍摄的标定图片p0，采用图像处理技术，得到图像中物体长度的像素点为l
p
，已知实际长度l0，根据公式
[0104][0105]
得到像素尺寸比lca。
[0106]
表1
[0107][0108]
如图4
‑
a和图4
‑
b所示，图中有11条线段，由于呈对称性，对称的线段长度像素点一样，则所有线段可用序号1、2、3、4、5、6表示。经过图像处理后获得各线段的最小外接矩形，分别以不同线段计算的尺寸比来验证其他线段的长度，并计算了它们的相对误差，其结果如图5所示，相对误差绝对值在0.5％以内，因此，该算法用来计算烟丝的尺寸等属性有效，都能够满足实际的要求。
[0109]
为了选取最优的像素尺寸比，计算了各个尺寸比相对误差的均值与均方差，如图6所示，线段5的相对误差均值与均方差皆为最小，因此将线段5的像素尺寸比lca：0.0918作
为本实验的计算尺寸比。
[0110]
2.烟丝形貌特征测量
[0111]
在烟丝形貌特征测量中，尺寸比为0.0918。基于烟丝形貌特征测量系统，读取图像p，输入尺寸比，计算出连通域中烟丝的对角线长度、面积、等效直径、球形度等参数，如图7
‑
11与表2所示。图7为烟丝图像处理流程图，分别对烟丝进行灰度化、滤波去噪、二值化、取反、忽略小面积烟丝、烟丝区域分割标记等，获取每个烟丝的基础信息。图8
‑
10分别绘制了烟丝群球形度、长度、当量直径的直方图，且图11绘制了烟丝的圆形度、矩形度、球形度的概率密度函数曲线。
[0112]
由图7
‑
a到图7
‑
f可以看出，烟丝原图中烟丝呈现暗黄色，由于光线并非完全均匀，中心成像清晰，两侧烟丝较暗。采用本发明专利提出的图像处理程序对烟丝图像进行灰度化与滤波去噪，图7
‑
c的图像识别度明显提高。进一步通过二值化操作，烟丝与背景完全分离，可以方便计算烟丝颗粒的面积。最后，对每幅待测图像进行取反、忽略小面积，与区域分割，每个烟丝颗粒被单独标记。基于颗粒像素，提取实际周长与面积，有利于后续烟丝形貌特征的分析。
[0113]
由图6球形度分布直方图与图9烟丝长度分布直方图可知，烟丝球形度集中在0.11附近，且左侧数据明显偏高。烟丝长度集中分布在23mm附近，由于图像测量考虑的样本较多，烟丝长度分布区间较为分散。
[0114]
基于公式：
[0115][0116]
计算了烟丝的当量直径(即等效直径)，当量直径分布直方图如图10所示。可以看出，经过多组烟丝样本的图像分析后，烟丝颗粒的当量直径呈近似正态分布，等效直径均值为1.7mm。
[0117]
烟丝的形貌特征分布包括圆形度、矩形度与球形度。图11
‑
a到图11
‑
c为形貌特征频率图。图11
‑
d到图11
‑
f为形貌特征概率密度分布拟合图，可以看出，圆形度、矩形度与球形度的均值依次为0.35、0.45、0.15。
[0118]
表2
[0119]
[0120]
通过本发明专利提出的烟丝形貌特征与等效直径测量系统，测量了多组烟丝图像。表2为烟丝形貌特征表，提供了一幅图像中52个烟丝颗粒的形貌特征与等效直径测量结果，从表中可以看出，各烟丝属性均在均值附近波动。