结合形态学重建和超像素的多特征FCM 分割算法①
马喃喃,  刘　丛
(上海理工大学 光电信息与计算机工程学院, 上海 200093)通讯作者: 马喃喃摘　要: 针对现有模糊聚类分割算法对噪声的鲁棒性差且提取的图像特征不充分等问题, 本文提出了一种结合形态学重建和超像素的多特征模糊 C-均值(FCM)分割算法. 首先, 利用形态学闭合重建处理原图像, 提高了算法的鲁棒性和细节保护能力. 其次, 采用Mean-Shift 方法预分割重建图像, 获得一组超像素区域. 再次, 提取重建图像各像素的颜特征、纹理特征和梯度特征, 利用平均策略定义各超像素的颜特征、纹理特征和梯度特征, 组成多维特征向
量. 最后, 运用最大熵正则化的加权模糊 C-均值算法(EWFCM)的框架, 以超像素为单位, 以核诱导距离作为距离度量来聚类多维特征向量. 选取BSDS300数据集中的6幅图像完成实验对比. 结果表明, 本文算法具有更高的分割精度.关键词: 形态学重建; 超像素; Mean-Shift; EWFCM; 核诱导距离
引用格式:  马喃喃，刘丛.结合形态学重建和超像素的多特征FCM 分割算法.计算机系统应用,2021,30(2):194–200. /1003-3254/7783.html
Multi-Feature FCM Segmentation Algorithm Combining Morphological Reconstruction and Superpixels
MA Nan-Nan, LIU Cong
(School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093,China)
Abstract : Aiming at the problems in the existing fuzzy clustering segmentation algorithms, such as poor noise robustness and insufficient image feature extraction, we propose a multi-feature FCM segmentation algorithm combining morphological reconstruction and superpixels. First, the original image is subject to morphological closing reconstruction,which improves the robustness and detail-preserving ability of the algorithm. Secondly, the mean-shift method is employed to pre-segment the re
constructed image and obtain a set of superpixels. Thirdly, the color, texture and gradient features of each superpixel in the reconstructed image are extracted and defined by an averaging strategy to form the multi-dimensional feature vectors. Finally, these vectors are clustered by using the framework of the EWFCM algorithm,taking superpixels as the unit and the nuclear induced distance as the distance measure. Furthermore, six images in the BSDS300 data set are selected for the experimental comparison. The results show that the algorithm in this study has higher segmentation accuracy.
Key words : morphological reconstruction; superpixel; Mean-Shift; maximum-Entropy-regularized Weighted Fuzzy C-Means (EWFCM); nuclear induced distance
计算机系统应用 ISSN 1003-3254, CODEN CSAOBN
E-mail: Computer Systems & Applications,2021,30(2):194−200 [doi: 10.15888/jki.csa.007783] ©中国科学院软件研究所版权所有.
Tel: +86-10-62661041
① 基金项目: 国家自然科学基金(61703278)
Foundation item: National Natural Science Foundation of China (61703278)
收稿时间: 2020-06-22; 修改时间: 2020-07-21; 采用时间: 2020-07-27; csa 在线出版时间: 2021-01-27
194
1  引言
图像分割是图像处理到图像分析的一个关键步骤，其目的为将一幅图像划分为若干个具有不同特性且有意义的区域, 每个区域具有相似的特征. 图像分割算法大致可分为基于阈值的分割[1]、基于聚类的分割[2]、基于区域的分割[3]和基于图论的分割[4]等. 其中, 基于聚类的分割方法以其快速、高效的特点, 越来越广泛的用于遥感图像、医学图像和自然图像的分割.
模糊C均值聚类(FCM)[5]及其改进算法是各种聚类技术中使用最广泛的方法. FCM算法为每一个样本的归属引入了模糊性, 保留了原始图像更多的信息. 但是该算法存在3个缺点: (1)只考虑像素灰度忽视了像素间的空间信息, 算法缺乏抗噪性和鲁棒性; (2)没有充分利用不同特征的优势; (3)使用欧氏距离作为距离度量.
为了提高传统F C M算法的抗噪性和鲁棒性, Ahmed等[6]在FCM目标函数中增加了空间信息约束，提出了FCM_S算法, 由于每次迭代时都要计算邻域，因此该算法时间复杂度较高. 为了降低FCM_S算法的复杂度, Chen等[7]提出FCM_S1和FCM_S2算法, 由于均值/中值滤波图像可以提前计算, 因此在算法迭代前就获得了邻域信息, 降低了计算成本. Zhong等[8]利用熵的概念提出了自适应结合空间信息的AFCM_S1
算法, 同时解决了手动调参问题. Lei等[9]提出了结合形态学重建与隶属度滤波的FRFCM分割算法, 该算法利用形态学重建平滑原图像, 提高了算法的抗噪性和细节保护能力.
为了充分提取图像特征, 一些学者提出了融合图像多特征的算法. Belongie等[10]提出了一种基于颜和纹理特征的图像分割算法, 该算法将三种颜特征和三种纹理特征嵌入到单个特征向量中. Yu等[11]提出了融合多特征的图像分割算法, 该算法通过亮度相似度、纹理相似度和边缘相似度来定义区域差异性. 尽管这些算法得到了可靠的分割结果, 但不同特征之间的权重调节仍存在问题. Rajaby等[12]提出的WHIFCM 算法使用图像的相和强度分量, 并通过自适应调整的权重将其组合到模糊C均值目标函数中.
很多分割算法大都是使用欧氏距离计算图像中目标的相似性, 该距离是一种线性度量. 近年, 研究者发现使用非线性度量可以更多地发现图像中的分布信息. Gong等[13]提出了KWFLICM算法. 该算法使用像素的空间距离和灰度值来重新定义权重因子, 并引入核诱导距离进行距离度量, 实现了更好的分割效果. Shang等[14]提出CKS_FCM算法, 该算法将欧氏距离度量替换为核诱导距离度量, 提高了图像分割的精确度.
基于上述3个问题, 本文提出了结合形态学重建和超像素的多特征FCM分割算法(SRMFCM). 首先使用形态学闭合重建和Mean-Shift超像素分割方法[15]预处理原始图像; 然后提取重建图像各像素的颜、纹理和梯度特征, 利用平均策略定义各超像素的颜、纹理和梯度特征. 最后, 运用区域代表像素点, 用核诱导距离进行距离度量, 实现多特征加权聚类, 合并具有相同标签的区域来获得最终分割结果.
2  相关工作
2.1  形态学重建
RεRδ
形态学重建能在不知道噪声类型的情况下较好的去除噪声并保留物体轮廓[9]. 最基本的形态学腐蚀重建和形态学膨胀重建如式(1)所示.
εδf
gε(1)
f
(g)=ε(g)∨fε(i)g(f)=ε
(
ε(i−1)(g)
)
∨fδ(1)
f
(g)=δ(g)∧f δ(i)g(f)=δ
(
δ(i−1)(g)
)
∧f∨∧
g=ε(f)
g≤f g=δ(f)g≥f
g=ε(f)g=δ(f)
其中, 表示腐蚀操作, 表示膨胀操作. 表示原始图像, 表示标记图像. 式(1)满足的条件包括
, , ,
, 其中和分别代表逐点的最大值和最小值. 在实际运用中, 满足膨胀重建中的条件, 满足腐蚀重建中的条件.因此, 为了简洁有效, 总是利用和来形成标记图像.
RεRδ
R O R C R C
R C
R C
和可以结合成有较强滤波能力的形态学开启重建算子和形态学闭合重建算子. 因为更适于平滑图像纹理细节,
本文选择形态学闭合重建来处理原始图像[9]. 如式(2)所示, 其中B表示腐蚀和膨胀操作所需的结构元.
2.2  超像素分割
超像素分割利用图像的局部相似性, 将图像分割成若干个具有相似颜、亮度和纹理等特征的图像块.相比于以像素为单位的图像分割方法, 以超像素为单位的图像分割方法更有利于提取图像局部特征来获取
2021 年 第 30 卷 第 2 期计算机系统应用
195
M ×N (R 1,R 2,···,R n )R 1∪R 2∪···∪
R n =P i  j R i ∩R j =∅更有效的图像信息, 同时降低计算复杂度. 超像素分割的基本原理是将一幅大小为的图像P 划分为n 个非空区域, 这些区域满足且当时, . 现有的超像素分割方法
有基于图论的方法和基于梯度下降的方法. 本文采用的基于梯度下降的Mean-Shift 算法, 该算法生成的超像素块不规则, 边缘贴合度好, 具有很好的抗噪性.2.3  EWFCM
X ={x 1,
x 2,···,x n }⊂R m EWFCM [16]是最大熵正则化的加权模糊C-均值算法, 该算法通过最小化类内离散度同时最大化属性权重熵来获得最佳的聚类结果. 给定样本集合, m 是样本空间的维数, n 是样本个数.
EWFCM 算法通过最小化式(3)将样本X 分为c 个类.
u ki ∑c
k =1
u ki =10≤u ki ≤1αc k j w k j ∑m
j =1
w k j =10≤w k j ≤1(
x i j −c k j )2x i j c k j λ(λ>0)λ其中, 表示第i 个样本属于第k 类的程度, 满足条件
,
. 表示模糊因子, 一般取2.
表示第k 类中第j 维特征的聚类中心. 表示第
k 类中第j 维特征所占的权重, 满足条件
,. 表示样本点与中心点的欧
氏距离. 式(3)第一部分加权的欧氏距离度量用来控制集的形状和大小, 促进了集的聚集. 第二部分是属性权值的负熵, 它根据现有数据对所有属性权值的最优分布进行了正则化. 是正则化可调参数. 只要选择合适的, 就能平衡式(3)的两部分, 得到最优稳定解.
3  本文算法
3.1  图像预处理
P ′P ′R ={R 1,R 2,···,R n }R i (1≤i ≤n )|R i |首先, 利用形态学闭合重建处理原始图像P , 得到重建图像; 然后, 按照Mean-Shift 算法将重建图像预分割成n 块区域, 写成集合形式为，每个区域的总像素数表示为.3.2  特征提取
P ′P ′HSV 颜模型是一种与人类视觉感知很相似的颜模型, 本文分别从RGB 空间和HSV 空间中提取图像的颜特征向量. 在RGB 空间中, 图像的每个
P x ,y [P r x ,y ,P g x ,y ,P b x ,y ]P ′P x ,y [P h x ,y ,P s x ,y ,P v x ,y ]
像素
都会有一个3维颜特征向量.
同理, 在HSV 空间中, 图像每个像素也有一个
3维颜特征向量. 每个区域中所有像素在各颜通道上的均值作为本区域的颜特征. 例如，区域R i 在R 颜通道上的颜特征如式(4)所示.同理可得该区域在其他颜通道上的颜特征向量.P ′P ′
gray P ′
gray P x ,y [P t 1x ,y ,P t 2x ,y ,P t 3x ,y ,P t 4x ,y ,P t 5x ,y ,P t 6x ,y ,P t 7x ,y ,P t 8x ,y ]Gabor 是用于边缘提取的线性滤波器, 它具有良好的方向选择和尺度选择特性, 对光照变化不敏感, 因此十分适合纹理
分析. 本文选择1个尺度和8个方向的Gabor 滤波器来提取图像的纹理特征. 首先将图像转
为灰度图像, 然后利用Gabor 滤波器对进行卷积, 则每个像素都得到8维的纹理特征向量, 每个区域中所有
像素在各方向上的纹理特征均值作为本区域的纹理特征, 例如, 区域R i 在方向t 1上的纹理特征如式(5)所示.
同理可得该区域在其它方向上的纹理特征向量.G x G y P ′
gray 经典的图像梯度算法是考虑图像的每个像素的某个邻域内的灰度变化, 利用边缘临近的一阶或二阶导数变化规律, 对原始图像中像素某个邻域设置梯度算子, 通常运用小区域模板进行卷积来计算. 通常有Sobel 算子、Robinson 算子、Laplace 算子等. 本文利用Sobel 算子生成梯度图像, Sobel 算子有两个, 一个用于检测
水平边缘, 另一个用于检测垂直边缘. 如式(6)和式(7)所示. 和分别代表图像经横向及纵向边缘检测的图像
灰度值.
P ′
gray P G 图像的梯度图像如式(8)所示.
区域R i 中所有像素的梯度均值作为本区域的梯度
计算机系统应用
2021 年 第 30 卷 第 2 期
196
特征向量, 如式(9)所示.
3.3  目标函数
R ={R 1,R 2,···,R n }n ×15原图像经过预处理后得到n 块区域, 写成集合形式为. 每个区域R i 都有15维的特征向量, 分别是6维颜特征向量、8维纹理特征向量和1维梯度特征向量, 所以整体的数据集规模为.本文最初的目标函数如式(10)所示.
Φx i j −Φc k j  Φ()Φ(x i j )T Φ(c k j )
=K (x i j ,c k j )其中, 是利用核诱导距离表示的归一
化量度. 表示非线性映射, 在特征空间的内积表示
为. 核替换的过程如式(11)所示.
最终的目标函数如式(12)所示.其中,
u ki k j k j 式(13)表示高斯核函数
, σ表示函数的宽带参数, 且满足k (x ,x )=1. 使用拉格朗日乘子法最小化式(12), 整理得的迭代公式如式(14)~式(16)所示.
3.4  算法流程
本文提出的SRMFCM 算法流程如算法1.
算法1. SRMFCM 算法
正则化是最小化策略的实现λ输入: 原始图像P , 聚类数目C , 模糊因子m , 参数, 最大迭代次数
max_iter 和高斯核kernel .输出: 图像分割结果Segm_P .
Step 1. 根据式(1)和式(2)对图像P 进行形态学闭合重建, 生成重建图像P'.
Step 2. 利用Mean-Shift 方法将图像P'分割成N 块区域, 区域代表像素点完成后续聚类.
Step 3. 根据式(4)提取各区域的颜特征. 根据式(5)提取各区域的纹理特征. 根据式(6)~式(9) 提取各区域的梯度特征. 将所有特征向量整合为一个特征向量矩阵.
U 0Step 4. 初始化隶属度矩阵, 设置循环计算器t =0.
Step 5. 用式(15)更新聚类中心; 用式(14)更新隶属度矩阵; 用式(16)更新权重矩阵; 用式(12)计算目标函数F ; t ++.Step 6. 如果t >max_iter , 转向Step 7, 否则返回Step 5.Step 7. 根据隶属度矩阵得到最终分割结果Segm_P .
4  实验结果与分析
为了验证本文所提SRMFCM 算法的有效性, 本节测试了BSDS300数据集的6幅自然图像. 本文选择的
对比方法有4种, 第1种结合形态学重建和隶属度滤波的FRFCM 算法; 第2种是基于像素点的用熵来调节类内紧凑度和像素空间信息的AFCM_S1算法; 第3种是基于核诱导距离度量的KFCM 算法; 第4种是基于区域的融合图像颜和纹理特征的FCM 算法, 简写为FCM(RC+RT).
图1中Image1-Image6表示6幅测试图像, 每幅图像的尺寸大小都为321×481像素. Image1-Image6的主要
分割目标分别为海星, 马, 老虎, 鹰, 蜻蜓和草坪.本实验过程中, 只需要分割出上述主要目标物体即可，其他部分都作为背景区域对待, 因此聚类数目都设为2.
图2–图7展示了使用上述5种算法对这6幅自然图像的分割结果, 其中GT 表示真值图像.
首先分析图2, 由于FRFCM 算法、AFCM_S1算法和KFCM 算法都是以像素为单位, 所以海星缺失了很多像素, 虽然FCM(RC+RT)算法和SRMFCM 算法
2021 年 第 30 卷 第 2 期
计算机系统应用
197
的分割结果优于前三者, 但是FCM(RC+RT)算法分割的海星中仍缺少一些像素, 使用SRMFCM 算法分割的结果中仍存在部分噪声.
(a) Image1(b) Image2(c) Image3
(d) Image4(f) Image5(g) Image6
图1    测试图像
(a) GT (b) FRFCM
(c) AFCM_S1
(d) KFCM (e) FCM (RC+RT)(f) SRMFCM
图2    不同算法对Image1的分割结果
对于图3, 草地影响了AFCM_S1算法、KFCM 算法和FCM(RC+RT)算法对马的分割效果. 虽然FRFCM 算法中草地的影响较小, 但是分割的目标缺少部分像素. 相比于前4种算法, SRMFCM 算法分割结果比较精确.
(a) GT (b) FRFCM (c) AFCM_S1
(d) KFCM
(e) FCM (RC+RT)
(f) SRMFCM
图3    不同算法对Image2的分割结果
对于图4, FRFCM 算法的分割结果较差. 其它4种算法都较准确的分割出了老虎, 但使用SRMFCM 算法分割的结果噪声最少.
(a) GT (b) FRFCM
(c) AFCM_S1
(d) KFCM
(e) FCM (RC+RT)
(f) SRMFCM
图4    不同算法对Image3的分割结果
对于图5, 由于原图中的天空颜由内向外逐渐加
深, 所以影响了AFCM_S1算法和KFCM 算法的分割效果. FRFCM 算法、FCM(RC+RT)算法和SRMFCM 算法的分割结果较好, 但是相比于
FRFCM 算法和FCM(RC+RT)
算法
, SRMFCM 算法对于鹰的尾巴有更
好的细节保留.
(a) GT (b) FRFCM (c) AFCM_S1
(d) KFCM (e) FCM (RC+RT)(f) SRMFCM
图5    不同算法对Image4的分割结果
对于图6, 算法都可以将蜻蜓的主体分割出来,
对
蜻蜓的足部等细节也有很好的保留, 但使用
SRMFCM 算法得到的分割结果最好
.
(a) GT
(b) FRFCM
(c) AFCM_S1
(d) KFCM (e) FCM (RC+RT)(f) SRMFCM
图6    不同算法对Image5的分割结果
对于图7, AFCM_S1算法和KFCM 算法分割的草坪中有很多噪声, 而使用FRFCM 算法、FCM(RC+RT)
计算机系统应用
2021 年 第 30 卷 第 2 期
198

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。