董 领，卿粼波，何小海，黄帅坤，何海波

(1 四川大学 电子信息学院, 成都 610065；2 成都西图科技有限公司，成都 610065)

0 引 言

岩心颗粒图像的分割结果会影响后续磨圆度、面积等参数分析。一幅岩心图像通常包含多种类型的矿物颗粒，传统的手工分割虽然相对精确，但是操作繁琐耗时，不能满足实际需求，故需要利用计算机来实现岩心颗粒图像的自动分割[1]。

图像分割算法主要有阈值分割、边缘检测分割、聚类分割、图论分割等。但这些方法用来提取岩心颗粒图像都不能取得比较好的分割结果。首先，岩心图像里面的颗粒颜色不一，且有些颗粒表面存在复杂的纹理；其次，颗粒之间的填隙物成分复杂，以致于有些颗粒与背景之间、颗粒与颗粒之间的界限模糊，很难进行区分[2]。

基于超像素块的分割方法近年来已被广泛应用于图像分割领域。宋熙煜等人对比了多种超像素算法，其中SLIC算法展示了目前最优的边界依从性[3]；任大勇等人提出SLIC与分水岭算法相结合进行分割，但依旧存在比较明显的过分割现象[4]；郭艳婕等人在分割过程中加入纹理特征，并应用DBSCAN算法合并超像素块，但用于岩心图像处理时效率低且合并效果并不理想[5]。

本文针对岩心颗粒在颜色和纹理方面存在的复杂性，提出一种岩心颗粒图像分割算法，使用结合纹理特征的改进SLIC算法生成超像素块，根据最近邻接图算法合并超像素块，能够对岩心颗粒获得较为准确的提取结果。

1 算法概述

算法流程图如图1所示。本文算法包括对岩心颗粒图像进行增强、SLIC分割算法的改进以及对过分割之后得到的超像素块进行合并3个部分。首先，采用结合中值滤波的CLAHE算法对岩心颗粒图像进行预处理，增强图像中的颗粒部分，同时模糊背景部分；其次，利用LBP纹理特征改进SLIC算法，对岩心颗粒图像进行过分割，使得到的超像素块更加贴合目标边缘；最后，利用颜色均值对岩心颗粒图像进行两步合并，最终得到岩心颗粒图像的边缘分割结果。

图1 算法流程图

2 图像增强

由于岩心颗粒图像里面的颗粒表面和岩心背景可能存在一些其它颜色的斑点信息干扰，以及图像本身的噪声干扰，这些情况可能增加图像的分割难度，故在对图像进行分割处理之前，需要对岩心颗粒图像进行图像增强操作。普通的直方图均衡化可能会使图像中的部分区域由于对比度增强过大而成为噪点。故本文使用一种结合中值滤波的CLAHE算法来对岩心颗粒图像的局部进行对比度增强，算法步骤如下：

(1)将图像I分成N×N个子区域，计算每个区域的直方图h(x)，再根据阈值T0对直方图进行调整。对新直方图每个子区域的四邻域进行插值运算，得到对比度增强后的图像I'；

(2)根据公式(1)获得图像I'的背景二值图Ibg，其中Im由公式(2)给出并根据公式(3)对Ibg进行形态学处理：

Ibg(x,y)=220,Im(x,y)>T∂

(1)

Im=|I'-M5(I')|+|I'-M7(I')|

(2)

(3)

(4)

图像增强结果如图2所示。可以看出本文算法能较好地对岩心颗粒图像进行局部增强，使颗粒部分更突出，岩心背景部分更模糊，更有利于后续的分割操作。

(a)增强前 (b)增强后

3 基于纹理特征改进的SLIC超像素分割

3.1 SLIC算法

SLIC算法是先将彩色图像转化为CIELAB颜色空间和XY坐标下的5维特征向量，构造距离度量标准进行局部聚类。算法具体实现如下：

(2)距离度量。以2S×2S为搜索范围，计算种子点与搜索范围内每个像素点之间的距离D，式(5)～式(7)：

(5)

(6)

(7)

其中，dlab是颜色距离；dxy是空间距离；m是调节空间距离的权重系数；S是种子点间的距离。可以看出，D值越小，则像素之间越相似。

(3)重复上述步骤，迭代更新聚类中心，直到误差收敛满足停止条件。

(4)增强连通性。将出现多连通、尺寸过小、不连续等情况的超像素块分配给邻近的超像素。

3.2 LBP纹理特征

岩心颗粒图像由于存在十分丰富的纹理信息，所以纹理特征对图像分割结果有重要影响。LBP(Local Binary Pattern)算子是一种图像局部纹理特征提取算法，计算的是每个像素和周围像素间“相对”的关系[6]。LBP值计算公式(8)如下：

(8)

其中，ic是中心像素点灰度值；ip是相邻像素点的灰度值；s是一个符号函数，式(9)：

(9)

本文使用圆形LBP特征，设置采样点为8个，当半径radius取1、3时LBP特征图像如图3所示。

(a)原始图像，左为lenna图，右为岩心颗粒图像

由图3可以看出，半径越小，LBP图越精细，所以本文选取半径为1，采样点为8的圆形LBP特征。

3.3 改进的SLIC算法

SLIC算法虽然在空间和颜色上具有良好的复杂度，但是由于不包含纹理特征，当目标与背景包含不同纹理信息但颜色相近时，不能很好地分割出理想的超像素块[7]。将LBP纹理特征引入SLIC算法，可以有效解决分割精度不准确问题。

改进的SLIC算法在初始化聚类中心时，需要加入LBP值，即Ck=[lk,ak,bk,xk,yk,LBPk]T，在距离度量步骤加入纹理距离，式(10)：

(10)

最终距离公式改为式(11)：

(11)

其中，n为调节纹理距离的权重系数。

加入纹理特征能够生成更好的贴合提取目标边缘轮廓的超像素块。如图4所示，(a)为使用原始SLIC算法分别在eagle图和岩心颗粒图像上超像素分割的结果，(b)为本文算法分别在两张图的分割结果。可见用原始SLIC算法进行分割，在目标和背景存在明显差异的情况下，某些部分依旧没有被分割出来，而本文算法所产生的超像素块比较有效地解决了分割中出现的欠分割问题，分割的结果更完整。

(a)原始SLIC算法

4 基于超像素合并的岩心颗粒边缘分割

4.1 最近邻接图构建

图像在经过超像素分割之后，被划分为K个区域，每个区域的像素具有同一个标签值。K个区域用区域邻接图(Region Adjacency Graph，RAG)表示为G=(V,E)，G是一个无向图，V是顶点集，每个顶点代表一个区域，E是边集，表示超像素块之间的相邻关系，边权重代表两个区域之间的相似度[8]。

区域邻接图由于计算效率低，本文用RAG的改进算法即最近邻接图(NNG)来搜索合并相似区域。NNG为一个有向图Gm=(Vm,Em)，其中Vm=V，Em为有向边的集合。边指向的是权重最小的相邻节点。RAG和NNG如图5所示。

图5 RAG和NNG示意图

NNG去掉了冗余的边信息，降低了算法的复杂度[9]。合并算法步骤如下：

(1)输入预分割图像；

(2)构建RAG，根据边权重生成对应的NNG；

(3)将每个环两端的区域进行合并，所有环合并完成后，更新NNG。

4.2 基于超像素块颜色均值的合并

颜色特征是基于像素点的特征，也是衡量图像相似度常用的特征之一。相邻超像素块之间的相似度可以用面积大小和颜色差来度量，计算公式(12)～(14)：

(12)

γ(A,B)=(vA-vB)2

(13)

γ(A,B)<ε

(14)

其中，vi表示超像素块中第i个颜色值；sizek表示区域面积；vk为区域的颜色均值；γ(A,B)表示区域相似度；ε表示阈值。

超像素分割之后，先计算每个超像素块的vk，再对超像素分割图像进行一处或多处背景选点，根据选点对应的标签值去找选点区域所属的超像素块，将其它超像素块的vk与背景点进行是否相似比较，如果相似就将其标签直接划分为背景进行背景区域融合，本文取背景阈值为0.1。对剩下的非背景区域，根据相似度γ(A,B)，利用最近邻接图进行合并，得到最终的合并结果。

5 实验结果与分析

5.1 实验结果

本文选用BSD数据集中的100张图像以及30张岩心颗粒图像进行实验，对比了使用SLIC+DBSCAN的算法[5]以及结合SLIC的分水岭算法[4]，分别选取2张BSD数据集图像和2张岩心颗粒图像作为展示，其分割结果如图6所示。

图6 分割结果对比

由图6可以看出，对于目标单一且背景差异明显的图，如图6中的第一张图，3种算法均能分割出目标；对于目标丰富、背景复杂的图像，如图6中的第二张图以及两张岩心颗粒图像，SLIC+DBSCAN算法存在明显的欠分割情况，有些地方边界缺失严重；结合SLIC的分水岭算法最终的融合效果不理想，某些地方过分割明显，且两种算法在岩心颗粒图像上效果均不好；而本文算法能在保留好边界的情况下，去除杂质，不仅能在公共图像上比较好地提取目标，也适用于岩心图像，提取的结果更完整。

5.2 实验分析

本文采用经典性能指标评价，查全率R(Recall)和查准率P(Precision)进行评估，计算公式(15)和公式(16)如下：

(15)

(16)

一般R和P值越大，表示分割性能越好，但是由于R和P值是反比关系，所以可以用F指数来衡量，F指数定义(17)如下：

(17)

由于本文算法提取的边界只占有一个像素，用来对比计算的人工分割结果是纯手工标注，具有一定的误差，故在计算客观评价标准时，允许4个像素的误差存在。本文从3个指标上和SLIC+DBSCAN的算法以及结合SLIC的分水岭算法进行对比，计算结果见表1，可以看出在公共数据图上，本文算法分割效果与对比算法相似，F指数只是略高于对比算法；而对于岩心颗粒图像，本文算法的分割结果明显优于其他算法，F指数也明显高于对比算法。因此，本文所提出的算法适用于岩心颗粒图像的提取，分割结果更准确。

表1 不同算法分割结果评估

6 结束语

本文提出了一种基于改进的SLIC岩心颗粒图像提取算法。在对岩心颗粒图像进行分割前，先结合中值滤波的CLAHE预处理，对图像进行局部增强，使岩心颗粒目标更明显，同时模糊岩心背景区域，让颗粒目标与背景更具有区分度；利用纹理特征对图像分割的影响，在对岩心颗粒图像进行超像素分割时加入LBP纹理特征，产生更加贴合目标边界的超像素块；最后利用最近邻接图对颜色特征相似的区域进行合并，最终分割出岩心颗粒。实验结果表明，本文算法对岩心颗粒图像的提取结果比较准确和完整，能有效减少后续对岩心颗粒进行分析的复杂度。