第7期 任会峰等.基于自适应谷底检测的浮选泡沫形态特征提取 785 基于自适应谷底检测的浮选泡沫形态特征提取 任会峰阳春华周 璇桂卫华鄢锋 (中南大学信息科学与工程学院,长沙410083) 摘要针对矿物浮选泡沫大小形状不一、光照不均等特点,提出一种基于适应度反馈微粒群优化算法 自适应优选谷底检测阈值的泡沫形态特征提取方法。将微粒群算法的惯性权重和加速因子设置为全局 最优点适应度的函数,优选谷底检测阈值,采用局部灰度极小选择边界检测模板,根据多角度逻辑规则 比较,获取气泡边界,分割结果的统计分析表明泡沫形态特征服从gamma分布。 关键词浮选泡沫形态特征微粒群谷底检测模糊划分熵 中图分类号TP391 文献标识码 A 文章编号 1000—3932(2011)07-0785-05 矿物浮选泡沫的视觉表观特征是矿物的物理 1 矿物泡沫图像的获取及特点 性质、浮选药剂、机械和操作条件等工况的综合反 泡沫浮选是气、液、固三相交互的复杂物理化 映…,其中泡沫的形态特征反映泡沫的矿物携带 学反应过程。在浮选机中添加各种药剂,与送人 量,直接关系到回收率和品位等重要的经济技术 浮选机的微粒矿物、工业用水等混合形成矿浆。 指标 。Moolman D W等采用灰度阈值法提取气 皮带轮带动叶轮高速旋转,风机压人的空气流在 泡中心的亮点来预测泡沫大小分布 J,但无法获 旋涡作用以及矿浆与气流垂直交叉运动的剪切作 得准确的泡沫边缘,且浮选泡沫噪声较大,亮点分 用下被分割成众多微小气泡,气泡在上升过程中 布不均,存在严重的过分割;Wang W X等采用基 与疏水的有价矿物粘附,在浮选机液面表层聚集 于流域的边缘探索法分割泡沫图像 j,该方法的 成为矿化泡沫,在刮板的作用下回收富集为精矿。 准确性受气泡类型和光照条件的影响;Kaartinen J 等采用分水岭算法分割锌浮选泡沫图像 J,该方 而脉石等亲水性矿物由于不能与气泡粘附而滞留 在矿浆中,作为尾矿排出。气泡越小,与微细矿粒 法依然受到亮点分布不均的影响,易造成过分割, 并且对于气泡大小不一的泡沫图像,该方法效率 的碰撞机会越多,表面势能越小,泡沫属于相对稳 急剧下降;阳春华等在此基础上先后提出聚类预 定的热力学体系,泡沫含矿量越高,越有利于矿物 分割及其改进算法 ’ ,有效降低了过分割与欠 的富集回收。亮度高且色温稳定的荧光灯作为光 分割的风险,但是计算量较大,缺乏分割效果的客 源为摄像机照明,同时降低环境光照的影响,光源 观评价标准;其它经典的边缘检测方法 都受到 尽量靠近摄像机,其安装高度略高于摄像机镜头 气泡亮点分布不均以及气泡之间灰度变化不明显 的高度,在泡沫表层上方1 10era处安装彩色CCD 的影响,导致边缘检测模糊移位、边缘丢失等问 摄像机拍摄泡沫视频,视野为12×9cm ,拍摄速 题。 率为7.5帧/s。图像视频信号通过光纤传输到工 笔者以矿物浮选泡沫为研究对象,首先通过 业计算机的图像采集卡,其将连续的模拟信号转 形态开闭重构预处理泡沫图像,降低噪声的影响, 换成离散的数字信号实时显示。 然后选择合适的边界检测模板,利用改进的微粒 图1是现场采集的一幅泡沫图像,从中可见, 群算法自适应优选谷底检测阈值,根据逻辑规则 摄像机视野内充满挤压层叠的泡沫,浮选气泡大 获取气泡边界,再通过细化、滤波等后期处理得到 小不一、形状各异,泡沫表面亮点分布不均、边界 最终的分割结果,并对所得的气泡形状和大小形 模糊,泡沫图像噪声较大。 态特征进行统计分析,证明泡沫形态特征服从 gamma分布。所提方法可以有效提高泡沫分割的 精度和效率,对混合尺寸和形状的泡沫分割具有 收稿日期:2011 5之4(修改稿) 基金项目:国家自然科学基金资助项目(60874069, 较强的鲁棒性。 61071176);国家杰出青年科学基金资助项目(61025015) 786 化工自动化及仪表 第38卷 图1 泡沫图像 2适应度反馈微粒群优化谷底检测的矿物泡沫 图像分割 泡沫图像含有的噪声是进行分割处理的障 碍,通过面积重构开闭运算等预处理 降低噪声 的影响。图2是泡沫灰度图像纵切面示意图,灰 度较低的像素点构成气泡的边缘,每个气泡的边 缘具有不同的灰度值,构成谷底分布特点。 图2 泡沫灰度图像纵切面示意图 2.1谷底检测算法 设 i√)表示矿物泡沫图像在像素点(i, )的 灰度值,i=1,2,…,M;j=1,2,…,Ⅳ。如图3所 示,定义边缘检测模板 有3 X 3个子模板 ,m =0,1,…;8表示子模板的序号,且 位于模板 的几何中心。每个子模板包含后X 个像素n…矗 =1,2,3,…;n=0,1,…,(k X k一1)。 4 3 x 5 x o x1 x b X1 X8 图3边缘检测模板 设g(i,J)∈{0,1}表示像素(i, )的状态,g ( )=0表示边界像素,g(i,J)=1表示非边界像 素。定义gn。(i, )表示f(i, )在0。的状态,F 为 子模板Xm的平均灰度值,Fm 磊。 。在 点(i√)的0。、45。、90。、135。4个方向分别进行灰 度值比较: F5>Fo—t<F1,g0。(i√,k,t)=0 > —t< ,g4s (i,j,k, )=0 一F7>一Fo—f<一F3,g90 ( √, ,£):0 (1) Fs>Fo—t< ,gl35。(iJ, ,f)=0 根据式(1)的逻辑规则确定气泡边界。其中 t为阈值,t太大,导致边缘不连续或丢失;t太小, 导致边缘太粗或者伪边缘,因此t合适的取值大 小是算法的关键。 2.2适应度反馈微粒群优化谷底检测阈值 针对传统的边缘检测算法的阈值选取问题, Kapur J N等提出一维直方图的最大熵法 J,但是 泡沫图像的信噪比较低,一维最大熵阈值分割的 效果不佳。Abutaleb A S提出二维最大熵法搜索 阈值 ,能够有效地抑制噪声,但运算速度慢、复 杂性高,很难应用于实时处理。 阈值优选的本质可以看作目标函数极值问题 的求解,Pai P Y等采用遗传算法动态获取最优的 分割阈值…J。笔者针对标准微粒群算法(Stand— ard Particle Swarm Optimization,SPSO)容易早熟, 并且惯性权重09线性递减导致进化后期算法收 敛速度明显下降¨ 的问题,采用适应度反馈的形 式,将惯性权重 和加速因子C 、C:设置为全局 最优点适应度的函数: 1/{1 exp[一 F(G^)] } (2) c,=1+I/{1+exp[一 ×F(G )] } =1,2 其中F(G )是算法第k次迭代所得的全局最 优点适应度,Ot=1/F(G ),n∈(1,2),一般n取值 为1.5即可。对粒子进行更新操作: “=∞ +c r。( 一 )+C2/'2( 一 ) … 【3) : + 其中k=1,2,…,m为迭代次数,惯性权重oJ 根据迭代次数线性递减,加速因子c 和c 取值为 2,r 和r:是均匀分布于[0,1]之间的两个随机 数。 由式(2)、(3)可见,如果全局适应度很大,粒 子距离全局最优点会很远,因为cc,和c数值较大, 所以粒子以较大的飞行速度探索整个解空间;反 之,当粒子接近全局最优点时,因为∞和c数值较 小,粒子以较小的移动步长开发局部解空间。 定义泡沫图像的背景和目标分别为: 第7期 任会峰等.基于自适应谷底检测的浮选泡沫形态特征提取 787 其中t为分割阈值,P 表示图像中灰度级i出 最优阈值取值,否则转步骤e; 现的概率,0≤i≤L一1,P = p ,其模糊划分 e.按式(2)计算 和Cl,C2; 熵¨ 为: f.按式(3)更新粒子的速度和位置,返回b。 珥 H =1og2pT +筹+ (5) 其中步骤C的具体过程为: FOR i=1:M 其中Ht=一 。log2P ,HT=一 P log2P。。 IF f(Xi)<f(Pi)THEN Pi=Xi,f(Pi)=f(Xi) IF f(Xi)<f(Ps)THEN Pg=Xi,f(Ps)=f(Xi) 当珥取得最大值时,对应的阈值就是最佳阈 END FOR 值,记为t ,因此选择适应度函数为: fitness(t )=max{HT} (6) 2.3泡沫图像分割结果 0≤I≤£一l 对适应度反馈微粒群优化谷底检测的矿物泡 令 、P 分别表示当前第i个粒子经历的最 优点和全局最优点,则采用FFPSO算法优选t的 沫边界进行细化和滤波,删除孤立点。图4是对 步骤如下: 图1泡沫图像采用不同分割方法的处理结果。从 a。初始化种群规模 、终止条件,为每个粒 中可见,利用阈值分割算法 分割泡沫图像很难 子随机产生一个位置 和速度vi,i=1,2,…, ; 确定气泡的边缘信息,更无法合理地统计出泡沫 b.按照式(6)计算粒子的适应度 X ); 的个数和面积。采用Vincent提出的经典分水岭 C.更新个体极值和全局极值; 分割法 ,易受干扰信号影响,造成过分割和欠 d.若满足终止条件,则P 对应的位置就是 分割。 a.阀值分割 b.分水岭分割 c优化的谷底检测 图4不同方法的泡沫分割结果 对分割后的图像进行像素标定得到气泡的实 表1 泡沫形态特征gamma分布统计参数 际面积,并以与气泡区域具有相同标准二阶中心 矩的椭圆的离心率描述气泡形状特征。图5是泡 沫大小和形状的统计分布直方图,利用极大似然 法进行gamma拟合,直方图的包络曲线为gamma 拟合效果,可见泡沫大小服从gamma分布,将泡 2.4泡沫图像分割效果评估 沫离心率向量 对数化以提高其对称性,证明泡 笔者采用一种区域均匀化测度 方法评价 沫形状也服从gamma分布。泡沫形态特征的统 分割效果,该方法无需先验分割图像,且评判更加 计分布参数见表1。 客观,具体计算方法为: ( 一 ,) ,u=1—2c — (7) 互 N( …一-,…厂 式中C——阈值数目; 『__一第 个气泡区域; 泡沫图像的像素总数; O 0 0.2 0.4 0.6 0.8 l O 气泡离心率(一log ) . 像素i的灰度级; 『__一第 个气泡区域的灰度均值; 图5泡沫形态特征的统计分布 ——泡沫图像的最大灰度值; 788 化工自动化及仪表 第38卷 ——泡沫图像的最小灰度值。 Vision based Control of Zinc Flotation‘a Case Study 显然,U数值越大表明分割效果越好。表2 [J].Control Engineering Practice,2006,14(12):1455 是采用文献[9]的阈值分割法,与标准的谷底检 —1466. 测及笔者提出的适应度反馈微粒群优化的谷底检 [6] 阳春华,杨尽英,牟学民等.基于聚类预分割和高低 测方法比较,可知笔者所提方法更加优越。 精度距离重构的彩色浮选泡沫图像分割[J].电子 与信息学报,2008,30(6):1286—1290. 表2泡沫图像分割效果评估 [7] 阳春华,周开军,牟学民等.基于计算机视觉的浮选 泡沫颜色及尺寸测量方法[J].仪器仪表学报, 2009,30(4):717—721. [8] 段瑞玲,李庆祥,李玉和.图像边缘检测方法研究综 述[J].光学技术,2005,31(3):415-419. 3结束语 [9] Kapur J N,Sahoo P K,Wong A K-C.A New Method for 浮选泡沫谷底检测阈值的选取受到泡沫形态 Graylevel Picture Thresholding Using the Entropy of 时变,以及气泡尺寸和形状混合分布特点的制约。 the Histogram[J].Computer Vision Graphics Image 通过适应度反馈的形式改进微粒群算法优选谷底 Process,1985,29(3):273—285. 检测阈值,提高了气泡边界检测的精度和效率,适 [10] Abutlaeb A S.Automatic Thresholding of Graylevel 合复杂多变的泡沫图像分割。对谷底检测气泡边 Pictures Using Two-dimension Entorpy[J].Comput- 界的泡沫形态特征进行统计分析,证明了泡沫的 er Vision Graphics and Image Process,1989,47(1): 22—32. 形态特征服从gamma分布,可以选择形状参数和 Pal P Y,Chang C C,Chan Y K et a1.An Adaptable 尺度参数描述泡沫形态特征,进而用于判断浮选 Threshold Detector,Information Sciences,201 1,181 工况,以及指导回收率或品位等经济技术指标的 (8):1463—1483. 预测,这也是下一步要研究的问题。 [12] 刘波,王凌,金以慧等.微粒群优化算法研究进展 参考文献 [J].化工自动化及仪表,2005,32(3):1—6. [13] Maitra M.Chatterjee A.A Hybrid Cooperative—COIn— Al&ich C,Marals C,ghean B J et a1.Online Monito. prchensivc Learning Based PSO Algorithm for Image ring and Control of Froth Flotation Systems with Ma— Segmentation Using Multilevel Thresholding[J].Ex— chine Vision:A Review[J].International Journal of pert Systems with Applications,2008,34(2):1341— Mineral Processing,2010,96(4):1~l3. 1350. [2] Bouchard J,Andr6 D,del Villar R et a1.Column Flo- [14] 路迈西,王勇,王凡等.利用阈值分割技术提取煤 ration Simulation and Control:An Overview[J].Min— 泥浮选泡沫图像的物理特征[J].煤炭科学技术, ends Engineering,2009,22(6):519~529. 2002,30(8):34—37. [3] Moolman D W,Aldrich C,Schmitz G P J.The lnterre. [15] Vincent L and SoiUe P.Watersheds in Diigtal 1ationship between Surface Froth Characteristics and Spaces:An Efifcient Algorithm Based on Immersion Industiral notation Performance『J].Minerals Engin. Simulations『J].IEEE Trans on Pattenr Analysis and erring,1996,9(8):837~854. Machine Intelligence,1991,13(6):583~598. [4] Wang W X,Bergholm F,Yang B.Froth Delineation [16] Zhang H,Fritts J E,Goldman S A.Image Scgmenta・ Based on Image Classiifcation[J].Mineral Engineer— tion Evaluation:A Survey of Unsupervised Methods ing,2003,16(11):1183~1192. [J].Computer Vision and Image Understanding, [5] Kaartinen J,Hatonen J,Hyotyniemi H et a1.Machine 2008,110(2):260—280. Morphological Characte rjstjcs Extraction of Mineral Flotation Froth Based on Adaptive Valley Edge Detection REN Hui—feng,YANG Chun—hua,ZHOU Xuan,GUI Wei—hua,YAN Feng (School ofInformation Science and Engineering,Central South University,Changsha 410083,China) (Continued on Page 806) 8O6 化工自动化及仪表 第38卷 [4] 李凤苏.苯及苯系物对作业人员健康状况的影响 [J].职业与健康,2008,24(17):1760~1762. 空气中挥发性有机物的测定研究[J].化工自动化 及仪表,2010,37(6):55—57. n H,Kim S J et a1.Rapid Determina- [12] Oh S Y,Du Shition of Floral Aroma Compounds of Lilac Blossom by [5] Badjagbo K,Loranger S,Moore S et a1.BTEX Expo— sures among Automobile Mechanics and Painters and Their Associated Health Risks[J].Hum Ecol Risk As— sess,2010,16(2):301—316. Fast Gas Chromatography Combined with Surface A— coustic Wave Sensor[J].Chromatographia,2008,79 (1/2):170~178. r J D,Bailey W F.Fehy R A et a1.Method for [13] Fai[6] 郭建民,傅国林,王晓宁.如何保证在线气相色谱仪 的分析准确性[J].化工自动化及仪表,2000,27 (5):59~60. Rapid on—site Identiifcation of VOCs[J].Envion r[7] 郭冰.气相色谱仪及其应用[J].石油化工自动化, Sic—China,2009,21(7):1005—1008. 2007,(5):88—89. [8] 傅若农.气相色谱近年的发展[J].色谱,2009,27 [14] Watson G W,Staples E J,Viswanathan S.Perform- ance Evaluation of a Surface Acoustic Wave Analyzer (5):584~591. [9] 刑家续,李东伟.931C在线气相色谱仪的应用[J]. to Measure VOCs in Air and Water[J].Enviorn 化工自动化及仪表,2001,28(4):77~78. Prog,2003,22(3):215—226. [10] Zhang S,Zhao T B,Xu X et a1.Determination of [{5] Staples E J。Viswanathan S.Development of a N0ve】 BTEX Compounds in Solid—Liquid Mixing Paint Using Odor Measurement System Using Gas Chromatography the Combination of Solid Phase Extraction,Thermal 山Surface Acoustic Wave Sensor[J].Air Waste Desorption and GC-FID[J].Chromatographia,2010, Manage,2008,58(12):1522~1528. 71(11/12):i131~l135. [16] 毛友安,刘巍,钟科军.电子鼻技术应用研究进展 [11] 周烈兴,钱天才,王绍华.光离子化气相色谱法对 [J].化学传感器,2009,(2):l2—17. Benzenes Determination with Electronic Nose GUO Xia,TIAN Sen.1in。NING Ping.JIANG Lei (Faculty ofEnvironmental Science and Engineering,Kunming University ofScience and Technology,Kunming 650093,China) Abstract The electronic nose(zNose4200,GC/SAW)and gas chromatography(GC522,FID)were respec— tively employed to analyze feasibility of detecting Benzenes(BTEX)with electronic nose.The results show that the sensitivity of electronic nose outperforms the general gas chromatography more than 10 times,and even up to 140 times;and its reproducibility and accuracy are also better than the generla gas chromatography and it can benefit volatile and semi-volatile organic compounds analysis. Key words electronic nose,generla gas chromatography,SAW,Benzenes (Continued from Page 788) Abstract Considering the different size and shape and poor illumination of mineral flotation froth,the mot- phologieal characteristics extraction of mineral flotation froth was proposed based Oil the adaptive valley edge detection of fitness feedback particles swarm optimization(FFPSO).Having the inertia weight and accelera- tion factors of PSO taken as the global minima function to choose best detection threshold.and having the local gray scale minima based to select edge detection template,and as well as having the multi-angle logical rules compared to detect froth boundary.Statistical analysis of segmentation results shows that the froth morphologi- cal characteristics are gamma distirbution. Key words flotation froth,morphological characteristics,particles swarnl optimization,valley edge detec— tion,fuzzy partition entropy