1、第40 卷第3 期2023年9 月doi:10.12084/j.issn.2096-3289.2023.03.009基于多尺度LBP算子的大理石纹理分割算法苏州科技大学学报(自然科学版)Journal of Suzhou University of Science and Technology(Natural Science Edition)Vol.40 No.3Sep.2023邹智伟，张战成*，姚浩男，徐少康（苏州科技大学电子与信息工程学院，江苏苏州2 150 0 9）摘要：纹理分割算法是仿大理石石材生产的重要步骤。大理石图像中各种纹理特征交杂在一起，当前主流纹理分割算法无法将其清楚、连续地
2、分割。受空洞卷积(Dilated Convolutions）的启发，提出一种多尺度LBP算子(MLBP),该算子将边缘和内部像素点组合排列，生成多个尺度算子，增大了特征提取的感受野，通过多个子尺度算子的融合能捕获到纹理的变化,可以适应大理石图像纹理风格多变的特点,增强特征提取的泛化能力。基于GMM聚类算法，在真实大理石数据集上验证了纹理分割效果，像素精度（PA）为9 3.2%，平均像素精度（MPA)为8 9.3%，平均交并比（M I o U)为7 1.2%，与FCN、K-m e a n s 和FCM等聚类算对比实验显示，MLBP特征算子优于传统的LBP特征算子。关键词：纹理分割；LBP算子；混
3、合高斯；大理石纹理中图分类号：TP391文献标志码：A文章编号：2 0 9 6-32 8 9(2 0 2 3)0 3-0 0 7 1-0 6“源于自然，超越自然”是建筑类人造石材装饰的价值追求。人造大理石的图案通常采集于天然石材，先经设计师修饰出图，然后经布料系统处理，分成多个颜色通道在模具中喷料制造。传统布料系统依赖人工对大理石图片进行纹理分割，成本高昂，同时分割结果也受个人主观因素影响，造成了新花色设计周期长、生产效果和原始设计效果偏差大等问题。利用图像处理技术对图片进行智能分割，代替传统的人工抠图，实现大理石的智能布料，是智能制造的发展方向。传统的图像分割算法只考虑了像素值信息而忽视图像
4、中形状信息和纹理信息，不符合纹理分割的要求。由于天然大理石纹理具有随机和不规则的特点，并且大理石纹理很难用高层语义特征描述，基于深度学习的FCN（Fu l l y C o n v o l u t i o n a l Ne t w o r k s）、A PC-Ne t（A d a p t i v e Py r a mi dContextNetwork）等场景分割算法很难直接用于大理石的纹理分割。目前纹理分割方法主要分为以下几类：基于灰度共生矩阵(Grey-Level Co-occurrenceMatrix,GLCM)的纹理分割算法2 、基于Gabor滤波器的纹理分割算法3和基于局部二值模式(Lo
5、cal BinaryPattern,LBP)的纹理分割算法困等。大理石图像纹理的底层视觉特征表示和聚类分割是实现分离的关键步骤。常用的特征描述算子有灰度共生矩阵、Gabor滤波器、LBP算子等。灰度共生矩阵是像素距离和角度的矩阵函数，它通过计算图像中一定距离和一定方向的两点灰度之间的相关性，来反映图像在方向、距离、变化幅度及快慢上的综合信息，由此代表纹理信息。Haralick等人0 提出了用灰度共生矩阵来描述纹理特征，用灰度和梯度的综合信息提取纹理特征，将图像的梯度信息加入到灰度共生矩阵中，使共生矩阵包含图像的纹理信息及排列信息，但计算过程太过穴余，计算效率低下。PratapT等人使用不同方向
6、的Gabor滤波器提取纹理特征，图像通过Gabor滤波器处理之后，在图像空间域上,将高斯核函数和正弦平面波相乘得到纹理特征信息。在频域上，纹理特征信息反映在不同方向和不同尺度上的多向性特征。Gabor能够提供良好的方向选择和尺度选择特性，而且对于光照变化不敏感，对光照变化有良好的适应性。如上所述，需要构造多尺度、多方向的Gabor滤波器，这会使图像尺寸、特征维度成倍增加,产生允余数据,使特征匹配更加困难。OjalaT等人8 则使用LBP算子提取纹理特征，该算子具有光照不敏感性、灰度不变性和旋转不变性，且鲁棒性强。不同于GLCM和Gabor滤波器计算收稿日期 2 0 2 1-12-17基金项目
7、国家自然科学基金面上项目（6 17 7 2 2 37)作者简介】邹智伟（1997 一）,男，江苏苏州人，硕士研究生,研究方向：计算机视觉与图像分割，E-mail;。*通信作者：张战成（197 7 一），男，博士，副教授，硕士生导师,E-mail:。72余、特征维数高，LBP算子计算时间快、消耗内存小,LBP算子的缺点是易受到噪声影响，为了克服该缺点出现多种LBP变体，如局部三值模式(Local TernaryPattern,LTP)、ELBPl 0 0 和CLBPl。但 LBP算子还是受算子尺度的影响，难以捕捉宏观纹理信息和纹理变化。常用的聚类算法主要有K均值聚类(K-means)12、模糊C
8、 均值(Fuzzy C-means,FCM)13和高斯混合模型(Gaussian Mixture Model,CMM)4I等。K-means算法的优势在于它的速度非常快，因为只需要计算点和簇中心之间的距离，它的缺点是必须人工选择聚类数量且聚类效果很大程度取决于聚类中心的初始位置；对于满足正态分布的数据,FCM算法的聚类效果更好，但对孤立点十分敏感且算法性能同样依赖初始点；GMM算法与常见的K-means聚类方法类似，但K-means属于硬聚类,结果中只包含每个像素点属于某个类的硬划分。在此基础上高斯混合模型则可以给出每个像素点属于多个类的概率，其包含的信息量远高于K-means和FCM算法，属
9、于软聚类。该模型优点在于，当某元素分配给两个簇的概率相同或非常接近时，可以通过主观选择或者设置可变值来获得更准确的结果。原始LBP算子在固定尺度上提取局部信息，存在两个问题：第一个问题是单一尺度不能有效的获取大尺寸且连续的纹理,解决该问题的一种方式是增大算子尺度15;第二个问题是在构造灰度特征向量时，仅将边缘像素和中心像素进行比较，而忽视边缘和中心之间的内部像素，丢失大量信息。受到深度学习中空洞卷积16 和特征金字塔17 的启发，将多分辨率特征联结，增大感受野从而捕获大尺度纹理特征。笔者提出MLBP（M u lt is c a le Lo c a l Bin a r y Pa t t e r
10、n)算子解决这两个问题。MLBP是多个子尺度LBP算子的叠加，不同于传统方法使用边缘的连续像素作为采样点，而使用子尺度采样点的组合，且选择大尺度算子的同时保留内部子尺度算子，减少特征丢失。在MLBP中，不同尺度的纹理特征会被不同的子尺度算子捕捉，不同子尺度算子能提取纹理的不同变化。1算法模型将MLBP算子提取的大理石图像纹理信息作为特征向量，通过高斯混合模型对特征向量分布进行分析，使用EM算法计算GMM参数，通过参数计算某个样本属于某个聚类的概率，挑选最大概率作为该样本的Label,则获得分割结果。算法整体流程如图1所示。苏州科技大学学报（自然科学版)Texturefeatureextract
11、ion2023年FeatureconnectionGMMClustering图1算法流程图MLBP算子基于LBP算子，如图2 所示，55MLBP算子由3个子尺度算子组成，分为外部子尺度算子和内部子尺度算子，不同于LBP算子使用16 个连续邻域像素作为采样点，而是将2 4个边缘采样点分解成3个邻域像素的组合，每个组合由8 个采样点组成，对应3个子尺度算子。计算过程如下：首先，将图像转为灰度矩阵，每个像素点是值为0,2 55 的一维向量，子尺度算子在灰度矩阵上对每个像素点进行计算，移动步长为1。计算公式为PLBPp.r(k)=Z gt(I(p)-(k)2,k e0,KP=1*）(*0,=0其中k代
12、表像素点，总像素点个数为K,I（*）是指该点的灰度值,p代表的是采样点（不同子尺度算子的采样点不同）,gt（*)代表门函数。将中心点和邻域采样点进行比较,大于中心点的为1，小于中心点的为0,类似(1)*0,=1(2)第3期门控逻辑。将8 位二进制数组合并转换为10 进制数，即该像素点的新灰度值，计算过程类似卷积操作，最后将多个子尺度算子提取的纹理向量级联得到最终纹理向量。邹智伟，等：基于多尺度LBP算子的大理石纹理分割算法73图2 MLBP55算子55MLBP是3个子尺寸算子的组合，7 7 尺度下是6 个（如图3所示），这不仅考虑边界像素，还考虑内部块的像素来构造特征向量，采用8 个采样点的组
13、合保证灰度值区间为0,2 55,通过不同大小的子尺度算子能增加感受野，使获取的特征更加完整连续，即单个子尺度算子提取不到的特征其他子尺度算子能够提取，领域像素的不同组合可以有效捕捉纹理特征的变化,脱离了原始LBP算子的局部性。图3MLBP77算子2实验结果及分析文中实验选取厂家提供的彩色大理石样板图作为数据集，将原图RGB颜色空间转为YUV颜色空间，选取Y通道进行纹理分割，包含两种像素尺寸规格：354517 6 18 和16 0 12 0 15，含有鱼肚白和北冰洋等不同风格的大理石图像。为了体现文中算法在纹理分割上的优越性,将该算法分割结果与基于LBP算子的图像分割算法FCM-LBP、K-m
14、e a n s-LBP和主流语义分割算法FCN对比,再详细对比不同纹理描述算子的灰度图结果来展现MLBP算子相比普通LBP算子的优势。742.1实验结果如图4与图5所示，大理石图像相比普通图像具有丰富的纹理语义信息和复杂的位置信息。如图4(c）和图5（c）,FC N语义分割算法很难将纹理特征从背景中分离，纹理特征发生错位，图像整体信息比较混乱，该语义分割算法不适合大理石纹理分割。LBP算子受限于固定尺寸，在捕获小纹理特征的同时会损失大纹理特征，导致提取的纹理特征不完整，如图4(d)-图4(h)中FCM+LBP、K-m e a n s+LBP和GMM+LBP分割算法结果中线条出现断续，纹理边缘特
15、征也有一定的损失；当分割目标差异较小时难以区分噪声和斑点类型纹理的区别,损失了大量的纹理信息；分割复杂纹理的边界时，在纹理和背景的界限上存在一定的过分割区域，分割时难以区分不同尺寸特征,噪声点和异常点偏多,鲁棒性较差。图4(f)-图4(g)受到聚类算法的影响,在数据较大时容易造成局部最优解，对孤立点非常敏感，分割结果存在许多噪声点，且分割效果非常依赖于聚类中心的初始化。文中基于MLBP提取纹理特征并利用GMM进行聚类，与FCM和K-means不同,GMM受噪声影响小,抗干扰能力强,如图4(h)-图4(i)。而MLBP加人了子尺度算子，能够提取不同尺度的纹理特征且捕获到纹理的变化，从而获得更加准
16、确的纹理特征，如图5(i)，清楚地将纹理线条和背景区分，纹理线条完整且连续，噪点少，近似标签图。(a)原图(b)Label(c)FCN(d)FCM+LBP(e)K-means+LBP(f)FCM+MLBP(g)K-means+MLBP(h)GMM+LBP(i)GMM+MLBP苏州科技大学学报（自然科学版）图4不同分割算法的分割结果2023年(a)原图如图6 所示，对比不同子尺度算子提取的纹理特征，可以看出每个子尺度算子捕获了不同尺度的纹理，大尺度算子获得较大的纹理信息，小尺度则捕捉小纹理信息。图6(a)中,LBP算子提取的特征局限于固定尺度，捕获的纹理特征并不完整，且纹理变化单一，而MLBP中
17、多个子尺度算子的结合完全提取图像的重要特征。在同一时刻,某一个子尺度捕获不到的纹理信息另一个子尺度可以捕获到,这是多分辨率提取特征的一个重要特性。观察图6(b)中各个结果,可以看出背景和前景纹理的变化,不同子尺度算子捕捉到纹理不同方向的变化，而LBP算子无法捕获纹理的变化。2.2实验分析为了更好地定量评价每个方法的分割效果，文中使用像素精度（Pixel Accuracy,PA）、均像素精度（MeanPixel Accuracy，M PA）、平均交并比（Mean Intersection over Union，M I o U)作为分割结果的评估标准。PA是计算标记正确的像素占总像素的比例；MPA
18、是一种逐像素标记的精度标准，它是PA像素精度的简单提升，计(b)Label(C)FCN()FCM+LBP(e)K-means+LBP()FCM+MLBP(g)K-means+MLBP(b)GMM+LBP()GMM+MLBP图5不同分割算法的细节对比第3期算每个类中被正确分类像素数的比例，之后求所有类的平均；MIoU指在每个类别上计算IoU值（真实值与预测值的交集/真实值与预测值的并集），最后进行平均操作。邹智伟，等：基于多尺度LBP算子的大理石纹理分割算法75(a)LBP从表1中可以看出，在大理石样板数据集的实验中,FCN语义分割算法精度较低,且空间位置上特征错位严重，不适合大理石纹理分割；将
19、LBP算子和FCM、K-me a n s、G M M 结合的算法优于FCN算法，但在细小纹理上分割上不够精确，难以将复杂的纹理从背景中分离出，导致纹理分割结果出现杂乱、断续等问题；文中提出的MLBP算子将多尺度纹理特征融合起来，能清晰、完整地将连续纹理特征提取，在不同纹理风格的大理石数据集中平均分割精度均是最优，基于MLBP的算法相比FCN分割算法均像素精度（MPA）提高了7.6%、8.1%、14%，对比基于LBP的算法提高了6.6%、6.4%、3.1%。最终文中算法比FCM-MLBP提高6.4%，比K-means-MLBP提高5.9%,均值正交比(MIoU)提高了8.6%、8%。表1三种算法
20、的客观指标评价的结果MethodPA/%FCN79.2FCM+LBP80.2K-means+LBP80.9FCM+MLBP86.8K-means+MLBP87.3GMM+LBP90.1GMM+MLBP93.23结语(b)MLBP(对应图（3)中各子尺度算子）图6 LBP、M LBP7 7 算子纹理灰度图对比MPA/%76.777.277.979.582.385.389.3MIoU/%45.656.558.762.663.268.971.2面向大理石图像的纹理特征提取,设计了改进的LBP算子MLBP,该算子通过融合多个子尺度算子提取的纹理信息，提高了LBP算子对纹理的描述性，同时保留LBP算子不
21、受图片的光照、角度等因素影响的优点。基于GMM对纹理图片的MLBP特征进行无监督聚类，实现了大理石纹理的自动分割。在真实的大理石图片上进行纹理分割实验，实验结果表明，基于MLBP方法结合GMM相比FCN、FC M 和K-means分割算法更加接近于标签分割结果，满足了工业生产的需求。文中算法仅在鱼肚花纹和麻石花纹上进行了实验，更多风格的纹理图像分割的效果有待进一步验证。参考文献：1王振华，胡伏原，吕凡，等.一种自适应增强图像纹理的边缘检测方法.苏州科技大学学报（自然科学版），2 0 18,35(2：6 8-7 2.2杨军，王恒亮.结合巴氏系数和灰度共生矩阵的遥感影像分割J.遥感信息，2 0 1
22、9,34(3：19-2 5.763王旭初，刘辉煌，牛彦敏。基于双流加权Gabor卷积网络融合的室内RGB-D图像语义分割.光学学报，2 0 2 0,40(19)：50-6 2.4 BASAR S,ALI M,OCHOA-RUIZ G,et al.A novel defocused image segmentation method based on PCNN and LBPJ.IEEE Access,2021,9:87219-87240.5彭丽丽，吴萍萍.基于改进水平集方法的腐蚀材料图像特征提取.苏州科技学院学报（自然科学版）,2 0 16,33(2)：45-50.6 AOUAT S,AIT-
23、HAMMI I,HAMOUCHENE I.A new approach for texture segmentation based on the gray level co-occurrence matrixJ.Multime-dia Tools and Applications,2021,80(2):24027-24052.7 PRATAP T,KOKIL P.Approximate optimization of gabor filter parameters in application to blood vessel segmentation in retinal imagesCJ/
24、IEEEInternational WIE Conference on Electrical and Computer Engineering.IEEE,2019:1-5.8 GARG M,DHIMAN G.A novel content-based image retrieval approach for classification using GLCM features and texture fused LBP variantsJ.Neural Computing and Applications,2021,33:1311-1328.9 SUKHIA K N,RIAZ M M,GHAF
25、OOR A,et al.Content-based remote sensing image retrieval using multiscale local termary pattemJ.DigitalSignal Processing,2020,104(5):102765.10 LIU L,ZHAO L,LONG Y,et al.Extended local binary paterns for texture classificationJ.Image and Vision Computing,2012,30(2):86-99.11 GUO Z,ZHANG L,ZHANG D.A co
26、mpleted modeling of local binary pattern operator for texture classificationJ.IEEE Transactions on ImageProcessing,2010,19(6):1657-1663.12 NITHYA A,APPATHURAI A,VENKATADRI N,et al.Kidney disease detection and segmentation using artificial neural network and multi-ker-nel k-means clustering for ultra
27、sound imagesJ.Measurement,2020,149:106952.13王荣淼,张峰峰,詹蔚，等。基于空间约束的模糊C均值聚类肝脏CT图像分割J.计算机应用,2 0 19,39(11):336 6-336 9.14 LI X,LI B,LIU F,et al.Segmentation of pulmonary nodules using a GMM fuzzy C-means algorithmJ.IEEE Access,2020,8:37541-37556.15 CAO Y,PRANATA S,YASUGI M,et al.Stagged multi-scale LBP fo
28、r pedestrian detectionCJ/IEEE International Conference on Image Process-ing.IEEE,2012:449-452.16 IM D,HAN D,CHOI S,et al.DT-CNN:Dilated and transposed convolution neural network accelerator for real-time image segmentation on mo-bile devicesCJ/IEEE International Symposium on Circuits and Systems.IEE
29、E,2019:1-5.17黄泽涛，刘洋，于成龙，等。时序特征融合的视频实例分割.中国图象图形学报，2 0 2 1,2 6(7)：16 92-17 0 3.Marble texture segmentation algorithm based on multi-scale LBP operator(School of Electronic&Information Engineering,SUST,Suzhou 215009,China)Abstract:Texture segmentation algorithm is an important step in the production of
30、 the ceramic tile which imi-tates the style of marble.Various texture features in marble images are intertwined,and existing mainstream tex-ture segmentation algorithms cannot segment them clearly and continuously.Inspired by Dilated Convolutions,amulti-scale LBP operator(MLBP)is proposed,which comb
31、ines the edge and internal pixels to generate multiplescale operators.Thus the receptive field of feature extraction is increased.Moreover,the texture variation can becaptured through the fusion of multiple sub-scale operators,accommodating to the polytropic texture style of mar-ble images and enhan
32、cing the generalization ability of feature extraction.Based on the GMM clustering algo-rithm,the texture segmentation effect is verified on the real marble dataset.The pixel accuracy(PA)is 93.2%,the average pixel accuracy(MPA)is 89.3%,and the average intersection ratio(MIoU)is 71.2%.Compared withFCN,K-means and FCM clustering on raw image and LBP,the proposed MLBP feature operator is better than thetraditional LBP feature operator.Key words:texture segmentation;LBP operator;GMM;marble texture苏州科技大学学报（自然科学版)ZOU Zhiwei,ZHANG Zhancheng,YAO Haonan,XU Shaokang2023年责任编辑：谢金春