基于ct图像的铝电解阴极炭块组分识别研究-刘庆生.pdf

《基于ct图像的铝电解阴极炭块组分识别研究-刘庆生.pdf》由会员分享，可在线阅读，更多相关《基于ct图像的铝电解阴极炭块组分识别研究-刘庆生.pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第卷 第4期 1 7年月 材 料 科 学 与 工 艺MAEIA IENE EHN Vl. N. 4A. 17I：1.111.n.1.1基于图像的铝电解阴极炭块组分识别研究刘庆生，曾少军，卢小能，宋翰林（江西理工大学冶金与化学工程学院，江西赣州41）摘 要：本文为铝电解用阴极炭块检测提供一种科学的方法，通过对铝电解用阴极炭块中炭骨料成分的颗粒级配计算，由骨料颗粒个数比例、面积比例分析铝电解用阴极炭块的质量.通过分析铝电解阴极炭块显微数字图像，基于数学形态学研发了一种实现其组分定量检测的技术方法.利用显微设备获取阴极炭块样品的数字图像及灰度值矩阵，基于显微成像原理区别出阴极炭块的组成成分为空隙、沥

2、青、炭骨料及杂质，通过对阴极炭块数字图像校正光密度、校正颜色、选择阀值、填充空隙、检测边缘、分离区域及进行腐蚀、膨胀、开运算与闭运算等其他数学形态学特征计算.实验得到样品实际组成与处理结果数值处于合理的偏差范围内，得出显微数字图像处理是研究铝电解阴极碳块组分识别及特征提取的简便可靠的方法.结果表明，基于断层扫描图片可以较准确地实现阴极炭块组成成分的虚拟识别.关键词：铝电解；显微；图像处理；阴极炭块；粒径级配中图分类号： 1文献标志码： A文章编号： 1（17）4 I nn，EN n， nn，N Hnln（l Mll n l Ennn，nx nv n n nl，n 41，n ）： In ， ln

3、 ll n lk n v n l n nl ll l n n n n n lk ln ll. l n lk ln ll n n l n l n n. l nl n lk， nv nn n ll. l n x n lk n n n . v， l， n n n n， nln ， ln n， l n， l ln， v lln， n， ln n xnn nl n， nn n n ln n l l n lln. l n nl l l n n nl vn n， nn l n l n ll nn nn n xn ln ll n lk. vn l ， n n l l vl nn nn n lk.： ln ll

4、； ； n； n lk； l n收稿日期： 11. 网络出版时间： 171.基金项目：国家自然科学基金资助项目（1411，141）.作者简介：刘庆生（17 ），男，博士，副教授.通信作者：刘庆生，El： l1 1.相比铝工业发达国家，我国的铝电解槽平均使用寿命要少很多，这使得我国铝工业生产过程中的成本增加，制约了铝电解工业的可持续发展.阴极炭块作为铝电解槽中最脆弱、也是最关键的部分，因其破损而引发的工业事故屡见不鲜.万方数据铝电解阴极炭块组分的差异会导致炭块性能有所不同，从而导致电解槽的使用寿命有所差异，因此需要通过科学的检测手段来检测炭块性能，并对配给方案进行优化延长电解槽使用寿命.现阶段对

5、碳素材料进行无损伤检测方法有超声波检测、射线检测和声发射检测.超声波检测是利用材料及其缺陷的声学性能差异对超声波传播波形反射情况和穿透时间的能量变化来检验材料内部缺陷的无损检测方法.然而由于碳素材料内部存在很多缺陷，造成超声波衰减十分严重，超声波检测无法准确地对碳素材料制品内部缺陷进行准确判断.采用射线能够检测到大规格阴极碳块的内部缺陷，并且能够准确反映缺陷的尺寸、形状和位置.但其局限在于射线机灵敏度，对尺寸过小的缺陷检测效果并不好.利用仪器接受声发射源释放的材料结构或缺陷信息，并对其进行分析称为声发射检测技术.传统的铝电解用炭块检测手段不仅实验流程繁琐费时费力，而且精确性不够，无法精确到颗粒

6、计算.本文提供了一种铝电解用阴极炭块质量检测技术，通过对铝电解使用的阴极炭块中炭骨料成分的颗粒级配计算，由骨料颗粒个数比例、面积比例分析铝电解用阴极炭块的质量.1 阴极炭块检测的基本思路图像处理流程见图1，其中，图像的组分识别和特征提取指对获得的炭块图像校正光密度和颜色、选择阀值、填充空隙、检测边缘、分离区域，进行腐蚀、膨胀、开运算与闭运算等，计算阴极炭块中炭骨料、杂质、空隙、沥青等具体含量；计算炭骨料颗粒级配.图像的获取 图像的预处理图像的组分识别和特征提取结果与讨论图 图像的处理流程.1 l n 图像的获取 阴极炭块的显微成像原理铝电解阴极炭块是由炭骨料、杂质、空隙、沥青等密度不同的成分组

7、成，利用样本中不同密度的部位对射线的吸收率不同，穿由透阴极炭块的不同密度部位的所需光强不同，结果反映为照射物质所对应的电压不同，最终在射线检测器上将阴极炭块成像1.同时，照射在阴极炭块表面的短波射线，仅有少量的射线被吸收，其他的则会穿过炭块原子间隙.阴极炭块的显微成像后各个像数点的吸收系数可通过测定阴极炭块经过显微成像器照射之前和之后的射线强度和进行计算得到.吸收定律的计算表达式为（，） In II （1）利用显微成像识别出铝电解阴极炭块组分为炭骨料、杂质、空隙、沥青.对图像进行形态学特征计算，得到铝电解阴极炭块中炭骨料的不同粒径比例.图为显微数字图像及对应的灰度值矩阵.图 显微数字图像（放大

8、）及对应炭块图像矩阵显示. M l （nl） n nn n lk x l 阴极炭块图像的采集阴极炭块样品的数字图像可以由nnVxl型显微设备获取，并可实现对铝电解阴极炭块的二维表征.选取铝电解阴极炭块试样时，要求被选用的铝电解用阴极炭块的质量较好，不存在较多缺陷及二次使用的情况.4材 料 科 学 与 工 艺 第卷 万方数据此外，为了更好地实现对铝电解阴极炭块的二维表征应选取圆柱试样，同时为有利于数字图像后续分析，应该获取侧面所投射的矩形图像.本文使用的试验样品是铝电解用圆柱形石墨质阴极炭块，尺寸为 . 阴极炭块图像的预处理对获取的铝电解用阴极炭块的数字图像进行矩形分割，运用高速高精度视觉处理技

9、术对AI（A In）计算范围进行自动检测.AI可以是任意的多边形形状，但为了便于后续图像处理，在操作中选择了隔离出静止数字图像内子集的矩形图像.利用Pnl等绘图软件环境下对获取的铝电解用炭块数字图像进行裁剪和消除黑色背景等操作，使阴极炭块图像的倾斜程度得以修正，并防止过度失真.图像修正结果如图所示.图 阴极炭块修正图像. n lk4 图像的特征提取和组分识别形态学运算部分的技术设计流程如图4所示.光密度校正颜色校正选择阀值孔洞填充粒径计算 腐蚀等处理 区域分离 边缘检测图 形态学运算技术的设计流程.4 nl n l ll n 光密度校正将彩色图像转化成为灰度图像，但实际上此操作效果不明显，因为

10、阴极炭块呈现黑色，彩色部分几乎没有.光密度是图像处理过程中所使用的单位，灰度值是数字图像形态学特征计算的中间变量，但是因为光密度难以被校正，不如校正上述步骤处理后的图像灰度值4.显微图像中亮度最高的位置光密度为时，没有光被吸收，透射率是1.炭块图像中亮度最低点的位置光密度趋近于.然而，图片中最亮的部位灰度值A不一定为，而是光密度为.因此为了防止测量数据出错，需要对图像的灰度值进行校正操作，将数字图像的光密度算出.灰度值与透过该物质的光强呈线性相关，而光密度和透过其光强呈指数关系，符合光吸收基本定律.其中图为灰度值直线，图为光密度曲线.250200150100500 50 100 150 200

11、 250灰度光强/cd图 灰度值直线. vl n ln2.52.01.51.00.5050 100 150 200 250光强/cd光密度图 光密度曲线. l n v 颜色校正颜色校正是调节修正阴极炭块图像在（亮度、对比度及伽玛值）、阴影、补偿视觉等方面的不足之处.显微设备显像管夹杂的散光会在一定程度上干扰数字图像的获取，所以使得阴极炭块数字图像的对比度和饱和度有所降低，通常将值定义为.4，但实际应用中，需根据具体情况通过P等数字图像软件调节值，解决阴极炭块图像局部太亮或者太淡等问题，获取清晰的数字图像；由于图像存在噪声及模糊的问题，需调节修正图像在阴影、补偿视觉处缺陷. 选择阀值阀值的选取方

12、法有很多7，选取不同阙值能第4期刘庆生，等：基于图像的铝电解阴极炭块组分识别研究万方数据区分出炭块各个组成成分.将预处理后的铝电解用炭块的图像的灰度值的直方图绘制出，发现具其有明显的双峰特性，截取谷底的灰度值作为阈值，背景灰度值可看作固定的数值，被拍摄的炭块图像中各个成分有类似的对比度，根据峰谷的范围，选定一个阈值范围分割阴极炭块数字图像.以图阴极炭块修正图像为例，其灰度值直方图如图7所示.1614121086420 100 200灰度像素个数/103图 灰度值直方图.7 vl 在IPP软件中选择阀值后对图像分割4次.以图阴极炭块修正图像为例，第1次选定灰度值为，分割出铝电解用炭块中的空洞区域

13、；第次选定灰度值为1，分割出铝电解用炭块中的沥青区域；第次选定灰度值为11，分割出铝电解用炭块中的炭骨料区域；第4次选定灰度值为1，分割出铝电解用炭块中的杂质区域. 空隙填充通过填充空隙获取没有空隙的炭骨料颗粒图像.选取粒径较小的黑色空隙进行填充，获取满足炭骨料颗粒级配的计算素材.为填充较大的空隙选取1 ，并且发现图像填充的效果良好.其中，图为二值化处理后未进行空隙填充的炭骨料颗粒图像，图为填充空隙后的炭骨料图像.图 未进行空隙填充的炭骨料颗粒图像. Pl n vlln图 空隙填充处理后的炭骨料颗粒图像. I n l v lln 边缘检测图像的分界处有少量颜色或亮度不相同的点，使得数字图像边缘

14、变得模糊.为了使这些边缘相邻点之间的有更加清晰明确的分界线，基于梯度值的变化进行边缘检测的计算1.基于l法的基本原理11，先定义一个向量（，），向量的方向就是对应函数增大的最大值方向，梯度的幅值为（，j） （）其中，和计算公式为 （ 4） （ 7 ）（） （ 1 ） （ 4）（4）可利用卷积模版来演示和，图1为l算子的领域像素点标记示意图.图 算子的领域像素点标记示意图.1 Pxl k k l l卷积模版中心的像素点是l算法的关键.通过搜索整个图像来获取数字图像的边缘曲线，通过识别、绘制炭骨料图像的具体形状来实现图像的分割.根据n法和l法进行边缘检测，炭块图像的分割基于n法中的l法.图11为l

15、法分割中的图像边界.图 法处理中的图像边界.11 n n l n材 料 科 学 与 工 艺 第卷 万方数据 区域分离通过灰度值可以分离炭块中的的炭骨料、杂质、空隙、沥青区域，即AI（需要计算的区域）的选择由灰度值决定1，以图阴极炭块修正图像为例的局部图像1如图1所示.图 阴极炭块修正图像的局部截图.1 A ll n n lk通过选定灰度值在1的图像分离出炭骨料，并用纯白色标定选定区域，如图1所示.为了利于观察及分析，将炭骨料标定为白色，背景标定为黑色.图14为炭骨料颗粒图像与区域分离后的炭骨料颗粒局部图像.对阴极炭块图像进行形态学计算，将不同粒径炭骨料进行区域分离，然后统计处理后的颗粒图像颗粒

16、面积和粒径.炭骨料图像可划分为4个级配：粉料、小颗粒、中颗粒及大颗粒，图1为其显示效果.图 选定阀值后的炭块数字图像（局部）.1 n lk l l l （ll）(a) (b)图 炭骨料颗粒图像.14 l ： （） vll ； （） ll n(a) (b) (c)图 区域分离处理后的炭骨料颗粒图像.1 n l nl n：（） l l； （） l； （） ll l 腐蚀、膨胀、开闭运算等形态学运算在数字图像处理过程中腐蚀发挥较大作用14，不仅能消除图像目标区域的边界点，使边界点向AI区域收缩，而且还能消去炭骨料图像中颗粒之间相互粘连.无数个点组成图像，当放大足够倍数时，便可观察到这些点，假定炭骨料

17、图像中每个点表示一个二维变量，这个二维变量由（m，）表示，元素的横纵坐标分别由m和表示.其中，腐蚀的表达式为 m， m （）式中：代表经过腐蚀处理后炭骨料图像的集合；代表腐蚀处理前的炭骨料图像的集合；m表示腐蚀操作后炭骨料图像的结构元素.在数字图像处理过程中通过膨胀来扩充物体的边界点1，扩充边界点对AI区域的选择，可用于填充空隙.假定炭骨料图像中每个点表示一个二维变量，这个二维变量由（m，）表示，元素的横纵坐标分别由m和表示.其中，膨胀的表达式为 m， m （）7第4期刘庆生，等：基于图像的铝电解阴极炭块组分识别研究万方数据 先对阴极炭块图像进行腐蚀处理，再进行膨胀处理，即开运算实质1.闭运算

18、处理步骤正好与开运算相反，即先膨胀后腐蚀.开运算用于处理图像中较小的突出部分，平滑炭骨料图像，而闭运算用于填充图像的微小的空隙，使边界平滑.通过腐蚀、膨胀、开运算、闭运算等操作使炭块数字图像中相互粘连的炭骨料颗粒分隔.图像的中心相互粘连的炭骨料颗粒经过腐蚀后明显分离开.图1为处理前后的局部图像.(a) (b)图 炭骨料颗粒局部图像.1 l n l：（） n n； （） n n 结果与讨论基于形态学的图像处理原理17，计算图像的颗粒粒径.通过nnk模型应证颗粒级配的理论1.nnk模型为 （7）式中：表示阴极炭块成分中炭骨料颗粒最大粒径；表示炭骨料颗粒的最小粒径；为炭骨料颗粒的分布系数；表示累计百

19、分比.炭块组分（炭骨料、杂质、空隙、沥青）的所占面积百分比由阀值计算得出.以图为例将其灰度值划分为4个级段： 、 1、1，计算出各阶段图像面积所占比例.将图像像素与图像中标尺按比例转化，可以计算出炭块各成分面积.其中像素面积由各灰度值所占方格数统计计算得出，具体面积则由像素面积转化而来.将阴极炭块实际各组分所占的面积比例与实验测得的数据进行对比，其中，实验测得孔隙、沥青、炭骨料、杂质含量分别为.1、1.、7.、.；实际组成成分中孔隙、沥青、炭骨料、杂质含量分别为.、1.1、7.、1.由统计结果可知，实验结果与实际组成的偏差值处于可允许的范围，因此这种计算方法是可行有效的.阴极炭块成分的计算结果

20、与实际成分的对比如图17所示.利用I软件可对每个颗粒进行面积或粒径统计计算1，根据粒径不同划分区间，然后统计不同区间内的颗粒数量及其面积，并将计算结果与实际结果进行对比.结果发现，实验测定的粉料、小颗粒、中颗粒、大颗粒所占面积比例分别为7.4、.、4.、.7；实际组成粉料、小颗粒、中颗粒、大颗粒所占面积比例分别为.、.、4.4，、.然而由于在实际测定中由于炭块是深黑色的，会对判断并统计相互黏连小颗粒产生了误差，因此在实际组成成分中，大颗粒的比例要小于实验测量中大颗粒的比例，而小颗粒恰好相反.图1和图1分别为统计的结果及对比结果.80706050403020100孔隙 沥青 炭骨 杂质实验测定实

21、际成分体积分数/%图 阴极炭块计算结果与实际成分比例对比图.17 n lln l n lk l n 250200150100500粉料 小骨料 中骨料 大骨料面积/(10-6m2)0.469789150.692319560.532564233.297953图 种炭骨料颗粒的面积分布示意图.1 kn n l n 材 料 科 学 与 工 艺 第卷 万方数据6050403020100粉料 大颗粒 中颗粒 小颗粒体积分数/%实验测定实际成分图 种类型炭骨料颗粒的计算个数与实际个数比例对比图.1 n ll n n ln n l 结 论1）对采集的铝电解阴极炭块显微数字图像利用数学形态学进行检测分析，利用

22、显微设备获取阴极炭块样品的数字图像及灰度值矩阵，由于射线穿透阴极炭块中空隙、沥青、炭骨料及杂质等不同密度组成成分时的光强不同，会导致深黑色的炭块图像中各组分的灰度值有所差异，为深层次炭骨料组分分析提供有效途径.）基于图像处理的形态学计算得出炭骨料颗粒不同面积的比例和级配个数的比例，将研究结果与实际组成数据对比，得到样品实际组成与处理结果数值处于合理的偏差范围内，因此数字图像处理技术可作为阴极炭块组分识别的有效途径.）计算机智能处理阴极炭块图像极大地提高了检测精度，同时大大降低了人为因素的误差，可快速对人工所无法完成的高精度的炭骨料颗粒的粒径和面积进行计算，并能以数据结果为基础进行后续的各类分析

23、和处理，且该识别及测试方法无需消耗大量的实验材料，节省了实验的成本.参考文献：1 AHIAV. n n n nn . I Vn n，1， 1（）： 1. 方黎勇，段建华，陈浩，等.基于显微图像的岩芯孔隙分形特征.强激光与粒子束，1，7（）：4.AN n， AN n， HEN H， l.l n . H P n Pl ，1，7（）： 4. 朱小洁，钱付平，张浩. 射线显微的应用现状及发展.化工新型材料，11（增刊1）：.H ， IAN n， HAN H.Aln n xnn M. N l Ml， 11（ 1）： .4 丁卫华，雷曼，郭瑞.混凝土图像的几何校正.理论与应用研究，1（）：1.IN W，

24、EI Mn， . n n n n . n Aln，1（）：1. HH P， HAAAI A. nn nlk ll EPM vln： n nn l n nl nnn n vn ll ln.En Pl nl， ， （）： 14 14. 姜维.浅谈数字图像处理技术及其应用.信息与电脑：理论版，1（）：111.IAN W. l n n ln. Inn n ： ，1（）：111.7 MAA ， HWAN W. nl n nn vl . IEEE nn nInn ， 1， （）： 174. AMIEN . Innnv n nnn nn. In， 7， 4（）：4. 李靖宇，穆伟斌，金成，等.图像分割在医学

25、图像处理中的应用研究.微型机与应用，1（）：1.I n， M Wn， IN n， l. nn n l n l n . M I Aln，1（）：1.1杨丽雯，曾朝阳，张永继.一种基于数学形态学的灰度图像边缘检测方法.国外电子测量技术，1（）：7.AN n， EN n， HAN n.M n l n l l . nEln Mn nl， 1 （ ）：7.11范彩霞，朱虹.基于物理模型的目标边缘检测方法西安理工大学学报，14（4）：4744.AN x， H Hn. n n l l. nl n nv nl， 14（4）：4744.1徐瑞.图像分割方法及性能评价综述.宁波工程学院学报，11（）：77. .A

26、nvv nn nn n vln . nl Nnnv nl，11（）：77.1王江涛，石红岩，李文.彩色图像分割算法综述.第4期刘庆生，等：基于图像的铝电解阴极炭块组分识别研究万方数据信息安全与技术，1（4）：7.WAN n， HI Hnn， I Wn. A v nl nn n. Inn n nl， 1（4）：7.14吴集林.论数字图像的膨胀与腐蚀.广东广播电视大学学报，7（1）： 117.W ln. l ln n n .nl n nv nn， 7（1）：117.1 IA ， IAN ， N ， l.n nl ln vl ln l nl . Innnlnl Mn l Mn： n， n Aln，1

27、，41（1）： 1111.1李忠民.基于图像处理的沥青混合料检测系统设计研究.西安：长安大学，1：4.17潘艳珠，吴文亮，王瑞宜.基于数字图像处理技术的沥青混合料级配离析评价方法.中外公路，11，1（4）：14.PAN n，W Wnln，WAN . nl n nl l xn n vln . nl n n H， 11，1（4）：14.1彭勇，孙立军，杨宇亮，等.一种基于数字图像处理技术的沥青混合料均匀性研究新方法.公路交通科技，4（11）：11PEN n， N n， AN ln， l. An n n l l nl. nl H nnn n vln， 4（11）：11.1张季如，祝杰，黄丽，等.土壤微观结构定量分析的IPP图像技术研究.武汉理工大学学报， （4）：.HAN ， H ， HAN ， l. IPP n nv nl l . nl Wn nv nl， （4）：.（编辑 程利冬）7材 料 科 学 与 工 艺 第卷 万方数据