基于自适应加权图像块的广义模糊C均值算法 下载: 855次
1 引言
在图像处理领域,图像分割是具有一定难度的基础性任务,在过去的几十年里,有众多的基于图像分割的研究成果。图像分割算法主要有阈值门限法[1-2]、聚类法[3-6]、分水岭方法[7]等。
模糊C均值(FCM)算法[8]作为一种软聚类方法,被广泛地应用于样本分类和图像分割中。然而,FCM算法并没有考虑图像像素的邻域信息,因此具有一定的噪声敏感性。为了提升FCM的鲁棒性,许多改进方法将邻域信息加入聚类的进程,如Ahmed等[3]将空间邻域项引入FCM的目标函数(FCM_S),随后,提出FCM_S的两种变体算法[4],思路在于用均值滤波和中值滤波替代空间邻域项。另外,基于灰度级的模糊C均值算法也得到了广泛的应用,其重点是利用邻域信息构造原图像的加权和图像[9-11],如加强型模糊C均值(EnFCM)算法[9]、快速广义模糊C均值(FGFCM)算法[10]等。上述算法在利用邻域信息时,共有的缺陷是需要利用经验来设置一些参数,如邻域项权重等。于是,为避免引入额外参数,Krinidis 等[12]提出基于局部信息模糊C均值(FLICM)算法,并在此基础上,提出基于FLICM算法的改进版本[13-14],但是这些算法的时间复杂度都较高。
此外,针对提升算法收敛速度的研究成果,如Zhu 等[15]构建广义的模糊C均值(GFCM)算法,与其加速原理相似的是,为了提升GFCM算法鲁棒性,Zhao等[16-17]提出了基于邻域广义模糊C均值(GFCM_S)算法以及相应的核化版本。需指出的是,收敛加速的原因是在GFCM、GFCM_S或核化版本的目标函数上引入新颖的关于隶属度的限制项,在迭代过程中奖励较高隶属度且惩罚较低隶属度。
图像块广泛地应用于图像降噪方面,一般来说,基于图像块的降噪方法比基于单个像素的降噪方法效果好,因为相比单个像素,图像块包含更多的信息和能够更好地描述图像[18]。已有少量文献从图像块角度构建模糊聚类目标函数[19],但是在构建过程中仍然存在其他问题待解决,比如图像块内各像素权重自适应确定问题以及构建的目标函数迭代加速求解问题。本文以图像块为基础,提出一种加速的自适应加权图像块模糊C均值算法。本文算法以图像块为基本单位构建新的目标函数,图像块内各像素权重由邻域像素与中心像素空间位置相对关系、图像块像素灰度关系自适应确定,同时在目标函数中引入GFCM算法中的隶属度限制项以加快算法迭代速度。利用模拟图像和实际图像对本文算法进行分割实验,取得了良好的分割效果。
2 广义的模糊聚类算法
最初GFCM算法由Zhu等[15]提出,其目标是加快聚类的迭代速度并保持较好的分类效果。假设I={xj|j=1,2,…,n}表示一幅灰度图像,n表示像素总数,xj为第j个像素。GFCM算法的目标函数为
式中:c为设定的聚类数;m为模糊因子;vi为第i个聚类中心;uij为隶属度,表征像素xj划分为第i个类的模糊度;U={uij}为c×N矩阵;V={vi}为c×1矩阵。
利用拉格朗日乘子法最小化(1)式,可得隶属度和聚类中心的迭代公式。
式中:α (0≤α<1)控制算法的收敛速度,当α=0时,GFCM算法退化为标准的FCM算法。
(3)、(4)式中修正系数aj的本质是相对地提升较高隶属度而减弱较低隶属度,在迭代过程中随聚类中心进行调整,GFCM算法的迭代流程可以参考文献[ 15]。
3 自适应加权图像块广义模糊C均值(GFCM_WP)算法
3.1 图像块各像素权重
对图像块做如下设定:Nj表示以像素xj为中心的图像块,其对应的灰度值为pjr(r∈Nj,表示图像块内第r个像素)。设图像块大小(长度)为q,显然,Nj是q×q的图像块,则pjr含q×q个像素灰度值(r=q×q)。
在每个图像块内,由于中心像素与邻域像素的位置相对关系不尽相同,同时各个像素灰度值大小都可能不同,那么各个像素在聚类中所起的作用也不相同,则必须为每个像素赋予不同的权重。图像块内各个像素的权重公式为
式中:wscn表示由图像块内中心像素和邻域像素空间位置关系(SCN)得到的权重;wgip表示由图像块内各像素灰度关系(GIP)得到的权重。二者综合作用得到图像块内各像素的权重。
首先对wscn进行说明。图像块内邻域像素与中心像素位置关系通常表示为
式中:djk表示空间欧氏距离,且djk只能表示像素xj和xk的相对关系。按照(6)式的定义,对任意图像块来说,wscn为固定矩阵,这显然是不合适的。事实上,在考虑djk的基础上再考虑像素xj和xk对应的图像块Nj和Nk之间的关系更能如实反映两像素的空间关系[14],因为如果图像块相似(同),则说明像素xj和xk在空间上连续。设计新的wscn表达式为
如(7)式所示,新的表达式添加了cos_sim(Nj,Nk)项,该项表示图像块Nj和Nk的余弦相似度,表达式为
由(8)式可知,根据图像块Nj和Nk的相关性强弱,cos_sim(Nj,Nk)的值在[0,1]之间波动。如果相关性越强,其值越大,则(7)式中分母的值越小,得到像素xj和xk的空间关系值wscn就越大,否则,反之。如此,由(7)式可自适应地对邻域像素与中心像素的空间关系赋予权重。
然后对wgip进行说明。对任意图像块Nj,对该图像块内各像素的变化情况进行衡量,公式为[11]
其中
式中:
式中:
由(12)式可知,当图像块内像素为边缘(噪声)点时,wgip值接近0,反之,wgip值较大,满足本文对权重赋值的期望。通过对wscn和wgip的说明,得到(5)式所示的综合权重,显然,需要对(5)式中总权重wjr进行归一化,
根据(13)式所示的权重的分配方法,给出几个具体的实例进行验证,选择有代表性的三个图像块A、B、C,如
图 1. 图像块权重分配实例。(a)选择的三个图像块;(b)图像块A的灰度值;(c)图像块B的灰度值(d) 图像块C的灰度值;(e)图像块A的wjr;(f)图像块B的wjr;(g)图像块C的wjr
Fig. 1. Examples of weight allocation of image patch. (a) Three selected image blocks; (b) gray value of image block A; (c) gray value of image block B; (d) gray value of image block C; (e) wjr of image block A; (f) wjr of image block B; (g) wjr of image block C
由
所以,本文的权重设置方法可以有效减小噪声和边缘等因素带来的冲击。另外,除了图像块大小q外,权重的计算并未引入其他参数,显示了算法的自适应性。需指出的是,
3.2 GFCM_WP算法目标函数及求解
如上节所示,将图像块内每个像素自适应分配权重后,构造的基于图像块的目标函数为
式中:vir为图像块内相应位置像素的聚类中心。需指出的是,(14)式的约束条件与(1)式相同,但是两式中的aj本质相同,而表达形式不同,(14)式中的aj稍后描述。
与其他模糊聚类算法相似,仍然利用拉格朗日乘子法对目标函数求解。首先构建(14)式的辅助函数:
分别对uij、vir求偏导,并令其为0,可得
由(16)式和约束条件
遵循GFCM算法的快速迭代思想,将(18)式中的表达成如下形式:
对比(19)式和(4)式可以看出,aj的表达形式由于不同的算法而有不同的表达形式。
3.3 GFCM_WP算法执行步骤
算法的执行流程可以总结如下。
输入:聚类数c,最大迭代次数T,终止条件ε,设置控制收敛速度系数α,图像块大小q,模糊指数m。
1) 根据(5)式求解每个图像块的权重wjr;
2) 随机初始化V(0)=[
3) For t=1 to T;
4) 将V(t-1)代入(19)式得到相应的
5) 将V(t-1)和
6) 将U(t)代入(17)式更新V(t);
7) If ‖V(t)-V(t-1)‖<ε或者q>T,迭代结束,Then执行步骤10;Else t=t+1,转步骤4;
8) Endif;
9) Endfor;
10) 输出U,由此得到图像块中心像素对应的类别。
4 实验结果与分析
4.1 实验说明和参数设置
为了展示GFCM_WP算法的分割表现,采用合成图像和实际图像进行分割实验。在本文中,应用GFCM[15]、KGFCM_S1[17]、KGFCM_S2[17]、EnFCM[9]、FGFCM[10]、NDFCM[6]、WIPFCM[19]与GFCM_WP进行对比。实验环境为:Matlab (R2014a)、3.40 GHz Intel© CoreTM i7-3770处理器,12 GB内存,Windows10企业版操作系统。为了对比公平,算法的共用参数应设置一致,如
表 1. 算法参数
Table 1. Parameters of algorithms
|
分割的评价指标如下。
1)当图像具有标准分割结果时,采用分割准确率(SA)[10]和调整兰德指数(ARI)[11]来评价,这两个指标均越大越好,表达式为
式中:Ai为算法得到的第i类的集合;Ci为标准分割结果中第i类的集合。
设R和T分别为算法得到的分割集合和标准分割结果,那么a、b、c、d分别为R∩T、
2)当图像没有标准分割结果时,采用一种基于熵信息的评价指标[14],其表达式为
式中:Hr(I)和Hl(I)分别为预期区域熵和布局熵。该有效性指标的思想是:分割应使每个分割区域内像素的均匀性最大化,并使区域间的均匀性最小,E越小,表明算法分割效果越好。
4.2 合成图像分割实验
人工合成了如
图 2. 合成图像上不同算法的分割结果。(a)原合成图像;(b)添加SPN(0.1)图像;(c) GFCM算法;(d) KGFCM_S1算法;(e) KGFCM_S2算法;(f) EnFCM算法;(g) FGFCM算法;(h) NDFCM算法;(i) WIPFCM算法;(j) GFCM_WP算法
Fig. 2. Segmentation results of different algorithms on the synthetic image. (a) Original synthetic image; (b) image with SPN(0.1); (c) GFCM algorithm; (d) KGFCM_S1 algorithm; (e) KGFCM_S2 algorithm; (f) EnFCM algorithm; (g) FGFCM algorithm; (h) NDFCM algorithm; (i) WIPFCM algorithm; (j) GFCM_WP algorithm
图 3. 合成图像上不同算法的分割结果。(a)添加WGN(0,0.01) & SPN(0.1)图像;(b) GFCM算法;(c) KGFCM_S1算法;(d) KGFCM_S2算法;(e) EnFCM算法;(f) FGFCM算法;(g) NDFCM算法;(h) WIPFCM算法;(i) GFCM_WP算法
Fig. 3. Segmentation results of different algorithms on the synthetic image. (a) Synthetic image with WGN (0,0.01) & SPN (0.1); (b) GFCM algorithm; (c) KGFCM_S1 algorithm; (d) KGFCM_S2 algorithm; (e) EnFCM algorithm; (f) FGFCM algorithm; (g) NDFCM algorithm; (h) WIPFCM algorithm; (i) GFCM_WP algorithm
表 2. 不同算法在合成图像上的分割结果
Table 2. Segmentation results of different algorithms on the synthetic image
|
4.3 实际图像分割实验
采用四幅图像进行分割评价,如
图 4. 待分割实际图像。(a) Bird图;(b) Bird标准分割图;(c) House图;(d) Coins图;(e) Rocks图
Fig. 4. Actual images to be segmented. (a) Bird image; (b) standard segmentation image of Bird; (c) House image; (d) Coins image; (e) Rocks image
由
图 5. 不同算法在Bird图像的分割结果。(a)添加SPN(0.2)图像;(b) GFCM算法;(c) KGFCM_S1算法;(d) KGFCM_S2算法;(e) EnFCM算法;(f) FGFCM算法;(g) NDFCM算法;(h) WIPFCM算法;(i) GFCM_WP算法
Fig. 5. Segmentation results of different algorithms on the Bird image. (a) Image with SPN(0.2); (b) GFCM algorithm; (c) KGFCM_S1 algorithm; (d) KGFCM_S2 algorithm; (e) EnFCM algorithm; (f) FGFCM algorithm; (g) NDFCM algorithm; (h) WIPFCM algorithm; (i) GFCM_WP algorithm
图 6. 不同算法在House图像的分割结果。(a)添加WGN(0,0.002) & SPN(0.05)图像;(b) GFCM算法;(c) KGFCM_S1算法;(d) KGFCM_S2算法;(e) EnFCM算法;(f) FGFCM算法;(g) NDFCM算法;(h) WIPFCM算法;(i) GFCM_WP算法
Fig. 6. Segmentation results of different algorithms on the House image. (a)Image with WGN(0,0.002) & SPN(0.05); (b) GFCM algorithm; (c) KGFCM_S1 algorithm; (d) KGFCM_S2 algorithm; (e) EnFCM algorithm; (f) FGFCM algorithm; (g) NDFCM algorithm; (h) WIPFCM algorithm; (i) GFCM_WP algorithm
图 7. 不同算法在Coins图像的分割结果。(a)添加SPN(0.1)图像;(b) GFCM算法;(c) KGFCM_S1算法;(d) KGFCM_S2算法;(e) EnFCM算法;(f) FGFCM算法;(g) NDFCM算法;(h) WIPFCM算法;(i) GFCM_WP算法
Fig. 7. Segmentation results of different algorithms on the Coins image. (a)Image with SPN(0.1); (b) GFCM algorithm; (c) KGFCM_S1 algorithm; (d) KGFCM_S2 algorithm; (e) EnFCM algorithm; (f) FGFCM algorithm; (g) NDFCM algorithm; (h) WIPFCM algorithm; (i) GFCM_WP algorithm
图 8. 不同算法在Rocks图像的分割结果。(a)添加SPN(0.1)图像;(b) GFCM算法;(c) KGFCM_S1算法;(d) KGFCM_S2算法;(e) EnFCM算法;(f) FGFCM算法;(g) NDFCM算法;(h) WIPFCM算法;(i) GFCM_WP算法
Fig. 8. Segmentation results of different algorithms on the Rocks image. (a)Image with SPN(0.1); (b) GFCM algorithm; (c) KGFCM_S1 algorithm; (d) KGFCM_S2 algorithm; (e) EnFCM algorithm; (f) FGFCM algorithm; (g) NDFCM algorithm; (h) WIPFCM algorithm; (i) GFCM_WP algorithm
表 3. 不同算法对真实图像的分割结果
Table 3. Segmentation results of different algorithms on the real images
|
5 结论
GFCM_WP算法直接以图像块为单位构建目标函数进而将邻域信息直接引入迭代进程。图像块内各像素权重由邻域像素和中心像素的空间关系以及图像块内各像素灰度关系综合确定,权重的确定过程体现了自适应性。利用新的目标函数得到的聚类中心和隶属度设计了算法的执行步骤。结合合成图像和实际图像对GFCM_WP算法和对比算法进行了分割性能测试,分割效果证明了所提算法的有效性。
[2] 聂方彦, 李建奇, 张平凤, 等. 一种基于Tsallis相对熵的图像分割阈值选取方法[J]. 激光与光电子学进展, 2017, 54(7): 071002.
[4] Chen S C, Zhang D Q. Robust image segmentation using FCM with spatial constraints based on new kernel-induced distance measure[J]. IEEE Transactions on Systems, Man, and Cybernetics. Part B, Cybernetics, 2004, 34(4): 1907-1916.
[5] 贾洪, 郑楚君, 李灿标, 等. 基于局部线结构约束的FCM聚类视网膜血管分割[J]. 光学学报, 2020, 40(9): 0910001.
[6] Guo F F, Wang X X, Shen J. Adaptive fuzzy c-means algorithm based on local noise detecting for image segmentation[J]. Iet Image Processing, 2016, 10(4): 272-279.
[7] 黄鸿, 金莹莹, 李政英, 等. 基于分水岭及半监督最小误差重构的荧光微球分割及分类方法[J]. 中国激光, 2018, 45(3): 0307013.
[8] Bezdek J C. A convergence theorem for the fuzzy ISODATA clustering algorithms[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1980, 2(1): 1-8.
[9] SzilagyiL, BenyoZ, Szilagyi SM, et al.MR brain image segmentation using an enhanced fuzzy C-means algorithm[C]//Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (IEEE Cat. No.03CH37439). September 17-21, 2003, Cancun, Mexico.New York: IEEE Press, 2003: 724- 726.
[10] Cai W L, Chen S C, Zhang D Q. Fast and robust fuzzy c-means clustering algorithms incorporating local information for image segmentation[J]. Pattern Recognition, 2007, 40(3): 825-838.
[11] 朱占龙, 王军芬. 基于自适应模糊C均值与后处理的图像分割算法[J]. 激光与光电子学进展, 2018, 55(1): 011004.
[12] Krinidis S, Chatzis V. A robust fuzzy local information C-means clustering algorithm[J]. IEEE Transactions on Image Processing, 2010, 19(5): 1328-1337.
[15] Zhu L, Chung F L, Wang S T. Generalized fuzzy C-means clustering algorithm with improved fuzzy partitions[J]. IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics), 2009, 39(3): 578-591.
[16] ZhaoF, Jiao LC. Spatial improved fuzzy c-means clustering for image segmentation[C]//Proceedings of 2011 International Conference on Electronic & Mechanical Engineering and Information Technology. August 12-14, 2011, Harbin, China.New York: IEEE Press, 2011: 4791- 4794.
[18] Liu Y L, Wang J, Chen X, et al. A robust and fast non-local means algorithm for image denoising[J]. Journal of Computer Science and Technology, 2008, 23(2): 270-279.
[20] 陆海青, 葛洪伟. 自适应灰度加权的鲁棒模糊C均值图像分割[J]. 智能系统学报, 2018, 13(4): 584-593.
Lu H Q, Ge H W. Adaptive gray-weighted robust fuzzy C-means algorithm for image segmentation[J]. CAAI Transactions on Intelligent Systems, 2018, 13(4): 584-593.
Article Outline
朱占龙, 董建彬, 李明亮, 郑一博, 王远. 基于自适应加权图像块的广义模糊C均值算法[J]. 激光与光电子学进展, 2020, 57(24): 241006. Zhanlong Zhu, Jianbin Dong, Mingliang Li, Yibo Zheng, Yuan Wang. Generalized Fuzzy C-Means for Image Segmentation Based on Adaptive Weighted Image Patch[J]. Laser & Optoelectronics Progress, 2020, 57(24): 241006.