激光与光电子学进展, 2020, 57 (24): 241006, 网络出版: 2020-12-02   

基于自适应加权图像块的广义模糊C均值算法 下载: 855次

Generalized Fuzzy C-Means for Image Segmentation Based on Adaptive Weighted Image Patch
朱占龙 1,2,3董建彬 1,2,3李明亮 1,2,3郑一博 2,3,*王远 2,3
作者单位
1 河北地质大学信息工程学院, 河北 石家庄 050031
2 河北省光电信息与地球探测技术重点实验室, 河北 石家庄 050031
3 河北省智能传感物联网工程研究中心, 河北 石家庄 050031
摘要
广义模糊C均值算法是一种比模糊C均值算法收敛速度更快的算法,然而它在分割灰度图像时对噪声敏感。为了改善其鲁棒性,提出基于图像块的像素灰度值加权的广义模糊C均值算法。该算法利用图像块代替单个像素构建目标函数,图像块内各像素的权重由邻域像素和中心像素空间关系及图像块内各像素灰度关系综合确定。以新目标函数为基础,利用拉格朗日乘子法推导出含图像块形式的隶属度和聚类中心表达式。通过这种方式,将邻域信息融入进聚类进程,提升算法的鲁棒性。利用合成图像和实际图像进行分割实验,结果表明:所提算法具有较强的鲁棒性和良好的分割性能。
Abstract
Generalized fuzzy C-means algorithm is a faster convergence algorithm than fuzzy C-means algorithm. However, it is sensitive to noise when segmenting gray images. In order to improve its robustness, a generalized fuzzy C-means algorithm based on the weighting of pixel gray value in image patch is proposed. In this algorithm, instead of a single pixel, the image patch is used to construct the objective function. The weight of each pixel in the image patch is determined by the spatial relationship between neighboring pixels and central pixel and the gray relationship of each pixel in the image patch. The expressions of membership and cluster center, in the form of image patch, are derived by using Lagrange multiplier method based on the new objective function. In this way, the neighborhood information is integrated into the clustering process, and then improves the robustness of the algorithm. Segmentation experiments are carried out with a synthetic image and several real images, and the segmentation results show that the proposed algorithm has strong robustness and good segmentation performance.

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算法的目标函数为

J(U,V)=i=1cj=1nuijmxj-vi2+j=1naji=1cuij(1-uijm-1),s.t.i=1cuij=1,uij(0,1),0j=1nuijn,(1)

式中:c为设定的聚类数;m为模糊因子;vi为第i个聚类中心;uij为隶属度,表征像素xj划分为第i个类的模糊度;U={uij}为c×N矩阵;V={vi}为c×1矩阵。

利用拉格朗日乘子法最小化(1)式,可得隶属度和聚类中心的迭代公式。

uij=xj-vi2-ajs=1c(xj-vs2-aj)-1/(m-1),(2)vi=j=1nuijmxj/j=1nuijm,(3)aj=α·min{xj-vl2|l{1,,c}},(4)

式中:α (0≤α<1)控制算法的收敛速度,当α=0时,GFCM算法退化为标准的FCM算法。

(3)、(4)式中修正系数aj的本质是相对地提升较高隶属度而减弱较低隶属度,在迭代过程中随聚类中心进行调整,GFCM算法的迭代流程可以参考文献[ 15]。

3 自适应加权图像块广义模糊C均值(GFCM_WP)算法

3.1 图像块各像素权重

对图像块做如下设定:Nj表示以像素xj为中心的图像块,其对应的灰度值为pjr(rNj,表示图像块内第r个像素)。设图像块大小(长度)为q,显然,Njq×q的图像块,则pjrq×q个像素灰度值(r=q×q)。

在每个图像块内,由于中心像素与邻域像素的位置相对关系不尽相同,同时各个像素灰度值大小都可能不同,那么各个像素在聚类中所起的作用也不相同,则必须为每个像素赋予不同的权重。图像块内各个像素的权重公式为

wjr=wscn·wgip,(5)

式中:wscn表示由图像块内中心像素和邻域像素空间位置关系(SCN)得到的权重;wgip表示由图像块内各像素灰度关系(GIP)得到的权重。二者综合作用得到图像块内各像素的权重。

首先对wscn进行说明。图像块内邻域像素与中心像素位置关系通常表示为

wscn=1/(1+djk),(6)

式中:djk表示空间欧氏距离,且djk只能表示像素xjxk的相对关系。按照(6)式的定义,对任意图像块来说,wscn为固定矩阵,这显然是不合适的。事实上,在考虑djk的基础上再考虑像素xjxk对应的图像块NjNk之间的关系更能如实反映两像素的空间关系[14],因为如果图像块相似(同),则说明像素xjxk在空间上连续。设计新的wscn表达式为

wscn=11+djk-cos_sim(Nj,Nk)(7)

如(7)式所示,新的表达式添加了cos_sim(Nj,Nk)项,该项表示图像块NjNk的余弦相似度,表达式为

cos_sim(Nj,Nk)=r=1q×q(pjr×pkr)r=1q×q(pjr2)·r=1q×q(pkr2)(8)

由(8)式可知,根据图像块NjNk的相关性强弱,cos_sim(Nj,Nk)的值在[0,1]之间波动。如果相关性越强,其值越大,则(7)式中分母的值越小,得到像素xjxk的空间关系值wscn就越大,否则,反之。如此,由(7)式可自适应地对邻域像素与中心像素的空间关系赋予权重。

然后对wgip进行说明。对任意图像块Nj,对该图像块内各像素的变化情况进行衡量,公式为[11]

cjr=pjr-rNjβrpjr/rNjβr,(9)

其中

βr=1,pjr-p-jrσj0,pjr-p-jr>σj,(10)

式中: p-jr表示图像块Nj的均值;σj表示该图像块的均方根。由(9)式可知,当图像块内各像素灰度值变化不大时,cjr值较小;当该图像块内像素位于边缘(噪声)区域时,cjr值较大。对图像块各像素对应的cjr取平均值[20],得

c-=rNjcjr/q2,(11)

式中: c-表示图像块内各像素变化的均值,这里用表达式(cjr- c-)衡量图像块各像素相对变化量。如果该表达式的值较小,则说明该值为非边缘点,该点对邻域的影响较小,希望增加其权重;否则,该值为边缘(噪声)点,希望减小其权重。采用下式赋予权重,

wgip=exp[-(cjr-c-)](12)

由(12)式可知,当图像块内像素为边缘(噪声)点时,wgip值接近0,反之,wgip值较大,满足本文对权重赋值的期望。通过对wscnwgip的说明,得到(5)式所示的综合权重,显然,需要对(5)式中总权重wjr进行归一化,

wjr=wscn·wgiprNjwscn·wgip(13)

根据(13)式所示的权重的分配方法,给出几个具体的实例进行验证,选择有代表性的三个图像块A、B、C,如图1(a)所示,图1(b)~(d)分别为对应的灰度值。

图 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

下载图片 查看所有图片

图1(b)~(d)可知,图像块A内的像素含椒盐噪声,但中心像素不是噪声点;图像块B内的像素为边缘点;图像块C内的像素含椒盐噪声,中心像素为盐噪声。同时,由图1(b)~(d)也可以看出,该示例所选择的三个图像块大小均为q=3,由(13)式计算得到三个图像块内各像素权重大小,如图1(e)~(g)所示。通过权重大小可以看出该图像块内对聚类结果影响较大的像素,以图像块A为例,对图像块有较大影响的像素是24、24和25。

所以,本文的权重设置方法可以有效减小噪声和边缘等因素带来的冲击。另外,除了图像块大小q外,权重的计算并未引入其他参数,显示了算法的自适应性。需指出的是,图1(e)~(g)中wjr为近似值,而且权重wjr可以事先计算,在聚类进程中采用权重值即可。

3.2 GFCM_WP算法目标函数及求解

如上节所示,将图像块内每个像素自适应分配权重后,构造的基于图像块的目标函数为

J(U,V)=i=1cj=1nuijmrNjwjrpjr-vir2+j=1naji=1cuij(1-uijm-1),(14)

式中:vir为图像块内相应位置像素的聚类中心。需指出的是,(14)式的约束条件与(1)式相同,但是两式中的aj本质相同,而表达形式不同,(14)式中的aj稍后描述。

与其他模糊聚类算法相似,仍然利用拉格朗日乘子法对目标函数求解。首先构建(14)式的辅助函数:

L(U,V)=J(U,V)+j=1nλji=1cuij-1(15)

分别对uijvir求偏导,并令其为0,可得

uij=-ajmrNjwjrpjr-vir2-aj1/(m-1),(16)vir=j=1nuijmwjrpjr/j=1nuijmwjr(17)

由(16)式和约束条件 i=1cuij=1可得:

uij=rNjwjrpjr-vir2-ajs=1crNjwjrpjr-vsr2-aj-1/m-1(18)

遵循GFCM算法的快速迭代思想,将(18)式中的表达成如下形式:

aj=α·minrNjwjrpjr-vlr2|l{1,,c}(19)

对比(19)式和(4)式可以看出,aj的表达形式由于不同的算法而有不同的表达形式。

3.3 GFCM_WP算法执行步骤

算法的执行流程可以总结如下。

输入:聚类数c,最大迭代次数T,终止条件ε,设置控制收敛速度系数α,图像块大小q,模糊指数m

1) 根据(5)式求解每个图像块的权重wjr;

2) 随机初始化V(0)=[ v1r(0), v2r(0),…, vcr(0)];

3) For t=1 to T;

4) 将V(t-1)代入(19)式得到相应的 aj(t);

5) 将V(t-1)aj(t)代入(18)式得到隶属度U(t);

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所示。

表 1. 算法参数

Table 1. Parameters of algorithms

AlgorithmParameter setting
mαβλαλsλgTεq
GFCM20.930010-4
KGFCM_S120.9430010-43
KGFCM_S220.9430010-43
EnFCM2430010-43
FGFCM23330010-43
NDFCM213330010-43
WIPFCM230010-43
GFCM_WP20.930010-43

查看所有表

分割的评价指标如下。

1)当图像具有标准分割结果时,采用分割准确率(SA)[10]和调整兰德指数(ARI)[11]来评价,这两个指标均越大越好,表达式为

NSA=i=1c(AiCi)/j=1cCj,(20)

式中:Ai为算法得到的第i类的集合;Ci为标准分割结果中第i类的集合。

NARI=2(ab-cd)(a+d)(d+b)(a+c)(c+b)(21)

RT分别为算法得到的分割集合和标准分割结果,那么abcd分别为RTR-T-RT-R-T

2)当图像没有标准分割结果时,采用一种基于熵信息的评价指标[14],其表达式为

E=Hr(I)+Hl(I),(22)

式中:Hr(I)和Hl(I)分别为预期区域熵和布局熵。该有效性指标的思想是:分割应使每个分割区域内像素的均匀性最大化,并使区域间的均匀性最小,E越小,表明算法分割效果越好。

4.2 合成图像分割实验

人工合成了如图2(a)所示的方块图,该图像包含4类,每一类灰度值分别为20、90、160和230,分别对该图像添加高斯白噪声(WGN)、椒盐噪声(SPN)和二者的混合噪声。采用表1的算法对含噪图像进行图像分割实验,所得的SA和ARI指标分别如图2、3和表2所示。由表2可以看出,在不同噪声的干扰下,GFCM_WP算法取得了较高的SA和ARI指标值。由图2、3的分割结果也可以看出,GFCM_WP算法得到了较好的视觉直观分割效果,尤其是在不同区域的交接处分割效果更好。

图 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

Noise levelIndexGFCMKGFCM_S1KGFCM_S2EnFCMFGFCMNDFCMWIPFCMGFCM_WP
WGN(0,0.01)SA0.87770.95610.98120.95580.97900.97030.97810.9903
ARI0.83690.94150.97350.94110.97200.96040.97090.9871
WGN(0,0.015)SA0.80440.95170.97340.95180.97310.96490.97050.9810
ARI0.73920.93550.96430.93580.96420.95320.96070.9747
WGN(0,0.02)SA0.75280.94210.96330.94210.96470.96030.95700.9659
ARI0.67040.92280.95110.92290.95300.94710.94260.9546
SPN (0.1)SA0.92310.88420.97010.88030.96490.96240.98410.9961
ARI0.89790.84560.96020.84040.95320.94990.97880.9948
SPN (0.2)SA0.85070.81370.93870.82110.91960.92720.97450.9890
ARI0.80090.75150.91830.76150.89280.90290.96600.9854
WGN(0,0.01) &SPN(0.1)SA0.81560.87010.94610.86960.95370.95600.94240.9838
ARI0.75420.82670.92810.82620.93820.94130.92320.9784

查看所有表

4.3 实际图像分割实验

采用四幅图像进行分割评价,如图4所示。图4(b)为图4(a)的标准分割结果,采用SA和ARI指标评价各算法;图4(c)~(e)没有标准分割结果,用上文提到的基于信息熵的指标E评价各算法。对图4中的图像分别添加不同程度的噪声进行分割测试,将图4(a)分割为2类、图4(c)分割为3类、图4(d)分割为3类、图4(e)分割为4类,其分割指标结果见表3,部分可视化分割结果见图5~8。

图 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

下载图片 查看所有图片

表3可以看出,在不同的噪声水平干扰下,对于Bird图,所提出的GFCM_WP算法得到的SA和ARI指标均为最大;对于House图和Coins图,GFCM_WP算法得到的E指标均最小,而对于Rocks图,GFCM_WP算法除了在WGN(0,0.005)情形下得到的E指标不是最小值外,其他情况下得到的E指标均最小。从视觉上看,如图5~8所示, GFCM_WP算法分割的区域更平滑,能够消除大多数噪声,保留更多的细节。

图 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

ImageNoise levelIndexGFCMKGFCM_S1KGFCM_S2EnFCMFGFCMNDFCMWIPFCMGFCM_WP
WGN(0,0.02)SA0.90070.95150.95350.96160.95600.95570.95560.9558
ARI0.80130.90300.90710.90310.91210.91130.91120.9114
WGN(0,0.03)SA0.85210.94960.95000.94970.95460.95500.95390.9550
ARI0.70420.89920.90000.89930.90850.90990.90780.9099
BirdSPN (0.1)SA0.91460.93370.95920.93430.95990.95930.95970.9626
ARI0.82920.83730.91850.86860.91980.91860.91940.9251
SPN (0.2)SA0.86720.89320.95480.89320.94830.95170.95760.9612
ARI0.73450.78630.90960.78630.89650.90350.91450.9225
WGN(0,0.02) &SPN(0.1)SA0.87680.92430.94610.92440.95070.95240.94600.9526
ARI0.75360.84860.89210.84880.90150.90470.89200.9053
Hr(I)5.09775.08775.08745.08984.23924.26514.23284.2191
WGN(0,0.002)Hl(I)1.04951.03591.04191.03951.03751.03991.03961.0388
E6.14736.12366.12936.12935.27675.30505.27235.2579
Hr(I)5.24115.22925.22895.23054.50524.49594.50284.4828
WGN(0,0.005)Hl(I)1.04261.03991.04421.04301.04251.04121.04201.0416
E6.28376.26916.27316.27355.54775.53715.54485.5244
Hr(I)4.46904.46194.46394.46883.61273.65673.59083.5832
HouseSPN (0.05)Hl(I)1.04511.03501.03621.04331.03721.04731.03641.0351
E5.51415.49695.50015.51214.64994.70404.62714.6183
Hr(I)4.41194.40044.40344.41023.61613.67383.60393.5713
SPN (0.1)Hl(I)1.05441.03911.04241.05451.03951.05521.04291.0345
E5.46635.43955.44585.46474.65554.72904.64684.6059
Hr(I)5.06795.06375.06495.06764.28944.29954.25644.2557
WGN(0,0.002) &SPN(0.05)Hl(I)1.04951.03641.04451.04931.04291.04861.03801.0383
E6.11746.10016.10956.11695.33235.34815.29455.2940
Hr(I)4.92224.82114.82104.82184.26694.26424.26654.2529
WGN(0,0.005)Hl(I)0.74560.71640.70680.73960.71670.73660.70930.7157
E5.66785.53745.52785.56134.98365.00084.97584.9686
Hr(I)4.83214.82964.82954.82984.32134.31224.32294.3089
WGN(0,0.01)Hl(I)0.75170.72110.71510.74280.72040.74130.71680.7203
E5.58385.55075.54465.57265.04175.05455.03965.0292
Hr(I)3.43573.43443.43363.43532.86422.84202.85162.8602
ImageNoise levelIndexGFCMKGFCM_S1KGFCM_S2EnFCMFGFCMNDFCMWIPFCMGFCM_WP
CoinsSPN (0.05)Hl(I)0.80710.71900.73200.80600.71810.80450.71250.7015
E4.24284.15344.16554.24133.58243.64653.56413.5617
Hr(I)3.44153.43073.43073.43212.89642.86652.89332.8884
SPN (0.1)Hl(I)0.89010.72730.77540.88250.72160.88370.71910.7033
E4.33164.15804.20614.31463.61813.75023.61303.5916
Hr(I)4.80114.73684.76404.76444.25284.20874.24454.2357

WGN(0,0.005) &
SPN(0.05)
Hl(I)0.81820.71800.74090.80740.71960.79820.70870.7161
E5.61935.48185.50495.57184.97245.00704.95324.9519
Hr(I)5.20015.19715.19225.19854.20274.20224.18054.1879
WGN(0,0.005)Hl(I)1.31571.31321.31351.31561.31441.31361.31341.3139
E6.51586.51036.51076.51415.51715.51585.49395.5018
Hr(I)5.25665.24375.24375.24524.38744.36614.38394.3400
WGN(0,0.01)Hl(I)1.31931.31511.31811.31941.31721.31691.31781.3166
E6.57596.55886.56186.56465.70465.68305.70175.6566
Hr(I)4.65504.64794.64944.65453.61133.65143.62883.6086
RocksSPN (0.05)Hl(I)1.31931.31221.31231.31911.31201.31771.31241.3104
E5.97445.96015.96175.97364.92334.96914.94124.9190
Hr(I)4.55524.54274.54524.55443.65933.59893.60823.5792
SPN (0.1)Hl(I)1.32721.31351.31441.32711.32341.31511.31271.3103
E5.88245.85625.85975.88144.98284.91404.92094.8895
Hr(I)5.12395.10695.10805.11114.20344.23414.18844.1966
WGN(0,0.005) &
SPN(0.05)
Hl(I)1.32611.31471.31761.32591.31381.32201.32341.3149
E6.45006.42156.42576.43705.51725.55615.51185.5115

查看所有表

5 结论

GFCM_WP算法直接以图像块为单位构建目标函数进而将邻域信息直接引入迭代进程。图像块内各像素权重由邻域像素和中心像素的空间关系以及图像块内各像素灰度关系综合确定,权重的确定过程体现了自适应性。利用新的目标函数得到的聚类中心和隶属度设计了算法的执行步骤。结合合成图像和实际图像对GFCM_WP算法和对比算法进行了分割性能测试,分割效果证明了所提算法的有效性。

参考文献

[1] Manikandan S, Ramar K, Iruthayarajan M W, et al. Multilevel thresholding for segmentation of medical brain images using real coded genetic algorithm[J]. Measurement, 2014, 47: 558-568.

[2] 聂方彦, 李建奇, 张平凤, 等. 一种基于Tsallis相对熵的图像分割阈值选取方法[J]. 激光与光电子学进展, 2017, 54(7): 071002.

    Nie F Y, Li J Q, Zhang P F, et al. A threshold selection method for image segmentation based on Tsallis relative entropy[J]. Laser & Optoelectronics Progress, 2017, 54(7): 071002.

[3] Ahmed M N, Yamany S M, Mohamed N, et al. A modified fuzzy c-means algorithm for bias field estimation and segmentation of MRI data[J]. IEEE Transactions on Medical Imaging, 2002, 21(3): 193-199.

[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.

    Jia H, Zheng C J, Li C B, et al. Retinal blood vessel segmentation based on fuzzy C-means clustering according to the local line structural constraints[J]. Acta Optica Sinica, 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.

    Huang H, Jin Y Y, Li Z Y, et al. Fluorescent microsphere segmentation and classification based on watershed and semi-supervised minor reconstruction error[J]. Chinese Journal of Lasers, 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.

    Zhu Z L, Wang J F. Image segmentation based on adaptive fuzzy C-means and post processing correction[J]. Laser & Optoelectronics Progress, 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.

[13] Zhang H, Wang Q M, Shi W Z, et al. A novel adaptive fuzzy local information-means clustering algorithm for remotely sensed imagery classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(9): 5057-5068.

[14] Gong M G, Liang Y, Shi J, et al. Fuzzy C-means clustering with local information and kernel metric for image segmentation[J]. IEEE Transactions on Image Processing, 2013, 22(2): 573-584.

[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.

[17] Zhao F, Jiao L C, Liu H Q. Kernel generalized fuzzy c-means clustering with spatial information for image segmentation[J]. Digital Signal Processing, 2013, 23(1): 184-199.

[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.

[19] Ji Z X, Xia Y, Chen Q, et al. Fuzzy c-means clustering with weighted image patch for image segmentation[J]. Applied Soft Computing, 2012, 12(6): 1659-1667.

[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.

朱占龙, 董建彬, 李明亮, 郑一博, 王远. 基于自适应加权图像块的广义模糊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.

本文已被 2 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!