基于全散度的自適應(yīng)魯棒圖形模糊聚類算法

吳成茂， 孫佳美

(西安郵電大學(xué) 電子工程學(xué)院， 陜西 西安 710121)

0 引言

圖像分割是將具有相似性特征的像素進(jìn)行分類或聚類的過程，它是數(shù)據(jù)挖掘、機(jī)器視覺和模式識(shí)別等領(lǐng)域的重要內(nèi)容，受到國(guó)內(nèi)外學(xué)者的高度重視。由于人的視覺特性和圖像本身所具有的模糊性，使得圖像分割問題成為典型的結(jié)構(gòu)不良問題，而模糊集理論具有描述這類問題的能力，將其應(yīng)用于圖像分割已成為一種必然趨勢(shì)，特別是模糊聚類分割方法[1-2]越來越受到廣大學(xué)者的關(guān)注。

傳統(tǒng)模糊c-均值聚類(FCM)[3]算法應(yīng)用于圖像分割，因僅考慮像素本身的聚類問題，忽視了當(dāng)前像素與其鄰域像素之間的相關(guān)性，導(dǎo)致對(duì)灰度分布不均勻的復(fù)雜醫(yī)學(xué)和遙感影像等難以獲得滿意的分割效果。為了解決FCM僅用隸屬度表達(dá)樣本分類信息的不足，Chaira[4]將FCM算法推廣至直覺模糊集并提出直覺模糊c-均值聚類(IFCM)算法，該算法相比FCM算法更有利于醫(yī)學(xué)影像分割需要；但是，該算法對(duì)猶豫度定義不夠明確，且算法分類仍然依賴于模糊隸屬度，不是嚴(yán)格意義上的IFCM算法。近年來，Cuong[5]提出了具有4個(gè)屬性參數(shù)的圖形模糊集，相比直覺模糊集能更有效地表達(dá)信息的不確定性和未知性，并將直覺模糊集視為其特例。隨后Son[6]將圖形模糊集引入傳統(tǒng)FCM，提出了一種基于圖形模糊集的模糊聚類 (PFCM)算法，它將樣本圖形模糊正隸屬度、中立度和拒分度相結(jié)合，來刻畫樣本分類信息的不確定性和未知性，但該算法存在嚴(yán)重的聚類不適定性，未得到相關(guān)學(xué)者的高度重視。近年來，Thong等[7]又提出了一種新的正則化圖形模糊聚類(PFCA)方法，解決了PFCM有關(guān)樣本圖形模糊中立度和拒分度存在的不適定性，對(duì)醫(yī)學(xué)影像和衛(wèi)星云圖等[8-9]分割能獲得滿意的效果。另外，Thong等[10]采用分而治之原理對(duì)PFCA的收斂性進(jìn)行了合理證明，給其應(yīng)用推廣奠定了堅(jiān)實(shí)的數(shù)學(xué)理論基礎(chǔ)。由于圖形模糊聚類是從FCM演化而來的計(jì)算智能方法，它與傳統(tǒng)FCM一樣對(duì)噪聲或奇異數(shù)據(jù)缺乏魯棒性。為此，基于圖形模糊聚類的改進(jìn)算法[11-13]也不斷涌現(xiàn)。吳成茂等[11]提出了一種對(duì)稱正則化圖形模糊聚類最優(yōu)化模型，將現(xiàn)有的FCM魯棒分割(FCM_s)算法[14]思想引入該模型，并得到了相應(yīng)的魯棒圖形模糊聚類分割方法，增強(qiáng)了現(xiàn)有的PFCA分割法對(duì)噪聲干擾的抑制能力。但是，這種利用平方歐氏距離導(dǎo)出的圖形模糊聚類分割法不僅難以滿足灰度空間分布不均勻的醫(yī)學(xué)和遙感等影像分割需要，而且對(duì)噪聲的抑制能力非常弱。

為此，本文引入一種具有旋轉(zhuǎn)不變性的全散度，用其替代PFCA中樣本與聚類中心之間的平方歐氏距離，獲得一種全散度圖形模糊聚類算法。為了進(jìn)一步增強(qiáng)該算法對(duì)噪聲的抑制能力，將像素空間鄰域信息嵌入該聚類算法并得到魯棒全散度圖形模糊聚類分割法；利用當(dāng)前聚類像素點(diǎn)與其鄰域像素均值之間的偏差自動(dòng)調(diào)節(jié)魯棒全散度圖形模糊聚類正則化參數(shù)。最終獲得基于全散度的自適應(yīng)魯棒圖形模糊聚類分割算法，以期對(duì)噪聲干擾及灰度分布不均勻圖像獲得更加滿意的分割效果。

1 PFCA模型

針對(duì)傳統(tǒng)FCM僅用模糊隸屬度表達(dá)樣本分類信息的不足，文獻(xiàn)[7]提出了一種中立度和拒分度相結(jié)合的正則化圖形模糊聚類最優(yōu)化模型，可描述為

(1)

s.t.

1) 0≤μik+ηik+ξik≤1, 0≤μik,ηik,ξik≤1,

i=1,2,…,n,k=1,2,…,c；

式中：Jm(μ,η,ξ,V)為正則化圖形聚類模型的目標(biāo)函數(shù)，Jm為模糊聚類函數(shù)映射關(guān)系表示符號(hào)，μ為樣本隸屬度矩陣，η為樣本中立度矩陣，ξ為樣本拒分度矩陣，V為聚類中心矩陣；μik表示樣本xi屬于第k類的圖形正隸屬度；ξik表示樣本xi屬于第k類的圖形模糊拒分度；n為樣本數(shù)；c為聚類數(shù)；m為模糊指數(shù)且常取值為2；xi為數(shù)據(jù)集X中第i個(gè)樣本數(shù)據(jù)；vk表示第k類的聚類中心；ηik表示樣本xi屬于第k類的圖形模糊中立度。

針對(duì)該聚類最優(yōu)化模型，采用拉格朗日乘子法可獲得樣本聚類的圖形模糊正隸屬度μik、中立度ηik和聚類中心vk的迭代表達(dá)式，而圖形模糊拒分度ξik則根據(jù)約束條件1并借助Yager補(bǔ)算子構(gòu)造所得。

2 基于全散度的PFCA法

2.1 Bregman散度與全散度

定義1假設(shè)可微函數(shù)Φ為嚴(yán)格凸函數(shù)，x和y是函數(shù)Φ上兩個(gè)不同的點(diǎn)，則函數(shù)Φ上x和y之間的Bregman散度DΦ(x,y)[15]為

(2)

定義2假設(shè)可微函數(shù)Φ為嚴(yán)格凸函數(shù)，x和y是函數(shù)Φ上兩個(gè)不同的點(diǎn)，則函數(shù)Φ的全Bregman散度δΦ(x,y)[16]為

(3)

文獻(xiàn)[16]指出，Bregman散度DΦ(x,y)會(huì)隨著坐標(biāo)軸的變化而發(fā)生改變，而全Bregman散度δΦ(x,y)則對(duì)坐標(biāo)軸旋轉(zhuǎn)變化具有不變性。因此，全Bregman散度相比Bregman散度具有一定的潛在優(yōu)勢(shì)，且已在醫(yī)學(xué)影像分析、圖像檢索和目標(biāo)檢測(cè)等[16-18]方面得到了成功應(yīng)用。

2.2 基于全散度的PFCA分割法

(4)

(5)

3 自適應(yīng)魯棒TPFCA

3.1 算法模型創(chuàng)建

為了進(jìn)一步增強(qiáng)TPFCA分割法的魯棒性，將文獻(xiàn)[14]中空間鄰域信息正則約束的魯棒聚類思想引入TPFCA，并提出具有一種噪聲抑制能力的TPFCA分割最優(yōu)化模型，其對(duì)應(yīng)的目標(biāo)函數(shù)為

(6)

(7)

其約束條件與模型(1)式完全相同。

3.2 算法最優(yōu)解推導(dǎo)

利用拉格朗日乘子法構(gòu)造(7)式的無約束最優(yōu)化目標(biāo)函數(shù)，具體表達(dá)式為

(8)

式中：λ為第2個(gè)約束條件的懲罰因子矩陣；β為第3個(gè)約束條件的懲罰因子矩陣；λi為第i個(gè)樣本滿足第2個(gè)約束的懲罰因子；βi為第i個(gè)樣本滿足第3個(gè)約束的懲罰因子。

目標(biāo)函數(shù)L(μ,η,ξ,V,λ,β)存在極小值的必要條件為

(9)

當(dāng)(8)式獲得最優(yōu)解時(shí)，樣本分類的隸屬度μik使得(9)式中的第1個(gè)方程式和第5個(gè)方程式同時(shí)成立，即

(10)

(11)

聯(lián)立(10)式和(11)式，獲得正隸屬度μik表達(dá)式為

(12)

當(dāng)(8)式獲得最優(yōu)解時(shí)，樣本分類的中立度ηik滿足1階必要條件(9)式中的第2個(gè)方程式和第5個(gè)方程式同時(shí)成立，即

(13)

(14)

聯(lián)立(13)式和(14)式，獲得中立度ηik表達(dá)式為

(15)

當(dāng)(8)式獲得最優(yōu)解時(shí)，則樣本分類的聚類中心vk滿足1階必要條件(9)式中的第3個(gè)方程式成立，即

(16)

由(16)式獲得聚類中心vk表達(dá)式為

(17)

另外，考慮到模型(7)式的第1個(gè)約束條件，利用Yager補(bǔ)算子獲得猶豫度ξik的表達(dá)式為

ξik=1-(μik+ηik)-(1-(μik+ηik)α)1/α.

(18)

由于最優(yōu)化模型(1)式所對(duì)應(yīng)的PFCA是收斂的，從而保證了嵌入像素空間鄰域信息的最優(yōu)化模型(7)式所對(duì)應(yīng)的基于ATPFAC也是收斂的。

3.3 算法實(shí)現(xiàn)步驟

本文針對(duì)最優(yōu)化模型(7)式所對(duì)應(yīng)的迭代求解表達(dá)式(12)式、(15)式、(17)式以及(18)式，設(shè)計(jì)了ATPFAC圖像分割算法。具體迭代步驟詳細(xì)描述如下：

步驟3選取圖像聚類分割區(qū)域數(shù)c，設(shè)置算法迭代終止誤差ε=10-4，算法迭代次數(shù)初值t=0，最大迭代次數(shù)tmax=103，Yager補(bǔ)算子參數(shù)α常選取區(qū)間為[0.6,1.0)。

步驟5利用(17)式計(jì)算不同類的聚類中心，即

式中：n表示圖像總像素?cái)?shù)，即n=MN.

步驟6利用(12)式計(jì)算樣本新的圖形模糊正隸屬度，即

步驟7利用(15)式計(jì)算樣本新的圖形模糊中立度，即

步驟8利用(18)式計(jì)算樣本新的圖形模糊拒分度，即

步驟9更新算法迭代次數(shù)t，即t←t+1.

步驟10重復(fù)上述步驟5～步驟9，直至t>tmax或max {‖ρ(t)-ρ(t-1)‖1}≤ε，ρ∈{μ，η，ξ}，終止算法迭代過程。

最后，利用本文聚類算法得到的樣本圖形模糊劃分信息，根據(jù)分類最大隸屬度原則對(duì)樣本像素進(jìn)行歸類標(biāo)記處理，完成圖像分割任務(wù)。

4 ATPFAC收斂性分析

針對(duì)ATPFAC目標(biāo)函數(shù)(7)式，將其分解為兩部分組成，即

其中

(19)

(19)式和最優(yōu)化模型(1)式的第2個(gè)約束條件相結(jié)合，得到ωik的表達(dá)式為

(20)

(21)

A=(aik,i1k1)(nc)×(nc)，B=(bij)c×c，

式中：

很顯然，矩陣A和B是正定的。

其次，針對(duì)(15)式，利用函數(shù)g(x)=(1-x)ex是單調(diào)遞減的，可證明(15)式所對(duì)應(yīng)的ηik是單調(diào)遞減的。再利用函數(shù)h(x)=1-x-(1-xα)1/α的單調(diào)性，可證明：當(dāng)μik+ηik<(0.5)1/α?xí)r，ξik是單調(diào)遞增的。再利用(15)式可證明下列3種情況成立：

綜合上述兩個(gè)方面的分析討論，可知本文ATPFCA算法是收斂的。

5 測(cè)試結(jié)果及分析

為了驗(yàn)證本文算法的有效性和魯棒性，首先采用標(biāo)準(zhǔn)灰度圖像和CT圖像進(jìn)行算法有效性測(cè)試，其次對(duì)噪聲干擾的合成圖像和遙感圖像進(jìn)行算法魯棒性測(cè)試，并與FCM、PFCA、TPFCA、FCM_s、模糊局部信息c-均值聚類(FLICM)[19]算法、PFCA_s (PFCA算法與文獻(xiàn)[14]魯棒算法思想相結(jié)合)以及文獻(xiàn)[11]算法進(jìn)行分割性能比較。實(shí)驗(yàn)測(cè)試環(huán)境為MATLAB R2014a，不同算法中所含參數(shù)及取值情況如表1所示，各算法迭代終止參數(shù)：迭代終止誤差ε取10-4、最大迭代次數(shù)tmax取103.

針對(duì)模糊聚類族算法，模糊指數(shù)m取1時(shí)，聚類算法屬于硬聚類，聚類效果較差；m取1.5～2.5時(shí)，算法聚類效果相對(duì)較優(yōu)，且對(duì)算法性能影響不顯著，m一般取2.

表1 不同算法所含參數(shù)情況

本文迭代終止誤差ε和最大迭代次數(shù)tmax所取數(shù)值均不會(huì)對(duì)各算法的收斂分割結(jié)果產(chǎn)生影響。對(duì)于采用平方歐氏距離測(cè)度的圖形模糊聚類算法，一般情況下，拒分度參數(shù)α的取值不會(huì)影響算法的分割結(jié)果，僅對(duì)算法收斂迭代次數(shù)和耗時(shí)有影響，且無明顯規(guī)律，通常取參數(shù)α在[0.5, 0.8]區(qū)間；對(duì)于采用改進(jìn)全散度距離測(cè)度的圖形模糊聚類算法，拒分度參數(shù)α的取值越大，分割結(jié)果的誤分率越大，峰值信噪比、迭代收斂次數(shù)以及迭代收斂耗時(shí)越小，參數(shù)α宜取[0.6, 1.0)。綜合兩類距離測(cè)度聚類算法，參數(shù)α取[0.6, 0.8]，本文選取α為0.68. 正則因子ω和鄰域窗口大小的取值均會(huì)影響算法分割結(jié)果，且對(duì)模糊聚類同族算法影響相同。二者取值越大，算法魯棒性越好，同時(shí)分割結(jié)果細(xì)節(jié)越模糊，因此正則因子取3，考慮到FLICM算法分割性能，鄰域窗口大小均選取5×5.

為了客觀評(píng)價(jià)不同算法的分割性能和噪聲抑制能力，本文采用文獻(xiàn)[20-21]提出的改進(jìn)峰值信噪比和誤分率對(duì)分割結(jié)果進(jìn)行定量評(píng)價(jià)。其中，改進(jìn)峰值信噪比是理想分割結(jié)果和算法分割結(jié)果的差異性度量，是一種與誤分率相反的刻畫圖像分割算法性能新的評(píng)價(jià)指標(biāo)，下文統(tǒng)稱為峰值信噪比。峰值信噪比(誤分率)的值越大(小)，表示分割結(jié)果越接近理想分割結(jié)果，算法分割性能越好；反之，算法分割性能越差。

5.1 無噪聲干擾圖像的測(cè)試與分析

針對(duì)Lena圖像、腦CT圖像和華盛頓遙感圖像，采用FCM、FCA、TPFCA、FCM_s、FLICM、PFCA_s、文獻(xiàn)[11]算法，以及ATPFCA對(duì)其進(jìn)行分割測(cè)試，結(jié)果如表2和表3所示。

從表2中可見，Lena圖采用FCM、PFCA、FCM_s、PFCA_s、FLICM以及文獻(xiàn)[11]算法所得結(jié)果均出現(xiàn)了嚴(yán)重的誤分現(xiàn)象，而本文TPFCA和ATPFCA算法在一定程度上保留了圖像輪廓以及細(xì)節(jié)部分，其中ATPFCA所得分割結(jié)果與理想分割結(jié)果相差不大。針對(duì)腦CT圖，采用FCM、PFCA、FCM_s、PFCA_s、FLICM以及文獻(xiàn)[11]算法所得結(jié)果中病變組織的完整性較差且存在一定噪聲，而本文TPFCA和ATPFCA相對(duì)其他6種算法所獲得病變部位區(qū)域更加清晰、完整，其中ATPFCA相對(duì)TPFCA算法噪聲點(diǎn)數(shù)明顯減少，病變部位也相對(duì)更完整，與實(shí)際的病變組織區(qū)域基本吻合(見圖1)，表明本文方法對(duì)復(fù)雜醫(yī)學(xué)影像分割具有一定的有效性和準(zhǔn)確性。針對(duì)華盛頓遙感圖像，采用TPFCA算法所得分割結(jié)果相比FCM和PFCA算法所得結(jié)果誤分點(diǎn)明顯減少，F(xiàn)CM_s、PFCA_s、FLICM以及文獻(xiàn)[11]算法所得分割結(jié)果的噪聲點(diǎn)數(shù)相對(duì)未加入鄰域信息的3種算法明顯減少，而本文算法ATPFCA所得分割結(jié)果相比前4種加入鄰域信息算法背景噪聲點(diǎn)顯著減少，目標(biāo)河流區(qū)域分割較為準(zhǔn)確，輪廓也更清晰。此外，從表3中不同算法分割結(jié)果所對(duì)應(yīng)的誤分率和峰值信噪比值可見，TPFCA算法對(duì)應(yīng)的誤分率(峰值信噪比)相比FCM和PFCA算法明顯減小(增大)，而嵌入鄰域信息正則化約束的ATPFCA算法誤分率或峰值信噪比值均小于或大于TPFCA、FCM_s、PFCA_s、FLICM和文獻(xiàn)[11]算法。結(jié)合表2、圖1和表3分析可得，本文算法ATPFCA更適合于分割灰度分布不均勻圖像，而且具有較強(qiáng)的抑制奇異數(shù)據(jù)對(duì)分割性能影響的能力。

表2 原始圖像及其分割結(jié)果

Tab.2 Original images and their segmentation results

表3 無噪聲圖像分割結(jié)果的誤分率及峰值信噪比

原圖像評(píng)價(jià)指標(biāo)FCM算法PFCA算法TPFCA算法FCM_s算法FLICM算法PFCA_s算法文獻(xiàn)[11]算法ATPFCA算法Lena圖誤分率/%信噪比/dB42.729.690 442.709.692 127.1611.653 243.669.536 544.349.486 043.659.536 943.349.622 826.1311.739 8腦CT圖誤分率/%信噪比/dB18.2713.397 517.4113.605 19.1616.382 617.4213.600 117.1713.662 917.3813.602 017.2413.648 17.6217.142 4華盛頓遙感圖誤分率/%信噪比/dB11.429.422 311.519.452 75.6212.504 04.4513.516 93.1215.064 24.3613.619 46.0412.192 22.6315.804 5

圖1 腦CT圖局部區(qū)域及分割結(jié)果Fig.1 Local region of brain CT image and its segmentation results

5.2 高斯噪聲干擾圖像的測(cè)試與分析

對(duì)圖2所示的3幅圖像分別添加均值為0且均方差為44、51和26(歸一化方差分別為0.03、0.04和0.01)的高斯噪聲，采用FCM、PFCA、TPFCA、FCM_s、FLICM、PFCA_s、文獻(xiàn)[11]算法，以及ATPFCA算法對(duì)其進(jìn)行分割測(cè)試，結(jié)果如表4和表5所示。

圖2 原始圖像Fig.2 Original images

表5 高斯噪聲圖像分割結(jié)果的誤分率及峰值信噪比

圖像名稱評(píng)價(jià)指標(biāo)FCM算法PFCA算法TPFCA算法FCM_s算法FLICM算法PFCA_s算法文獻(xiàn)[11]算法ATPFCA算法人工合成圖誤分率/%信噪比/dB45.8814.213 846.0914.206 146.6614.067 627.8823.271 730.7120.522 227.5023.342 526.7824.670 125.4027.191 7機(jī)場(chǎng)遙感圖誤分率/%信噪比/dB33.974.688 633.544.744 729.255.338 613.998.542 512.828.921 513.878.578 210.869.641 87.5911.196 9工業(yè)遙感圖誤分率/%信噪比/dB34.4913.919 434.0913.966 731.0314.463 021.2316.260 521.3516.233 620.6216.387 519.9916.520 915.6617.569 2

由表4所示的8種算法所得分割結(jié)果可見：嵌入鄰域信息正則約束的FCM_s、PFCA_s、FLICM、文獻(xiàn)[11]算法和本文算法ATPFCA所得分割結(jié)果圖，相比未嵌入鄰域信息的FCM、PFCA和TPFCA算法噪聲點(diǎn)數(shù)明顯減少，且區(qū)域分割效果也相對(duì)較優(yōu)；而嵌入鄰域信息的5種算法，其中本文算法ATPFCA相比FCM_s、PFCA_s、FLICM和文獻(xiàn)[11]算法，所得分割結(jié)果噪聲點(diǎn)數(shù)均顯著減少，表明其具有更強(qiáng)的抑制噪聲干擾能力。另外，從圖3所示的工業(yè)遙感局部區(qū)域放大圖來看，采用FCM、PFCA和TPFCA算法得到的分割結(jié)果中均含有大量的噪聲點(diǎn)且各區(qū)域無明顯分界，采用FCM_s、FLICM、PFCA_s和文獻(xiàn)[11]算法得到的分割結(jié)果中噪聲點(diǎn)數(shù)相對(duì)前4種算法較少，各區(qū)域邊界也相對(duì)較為清晰，而采用ATPFCA算法得到的分割結(jié)果中幾乎沒有噪聲點(diǎn)，且其區(qū)域邊界輪廓也更為清晰。此外，從表5所示的誤分率和峰值信噪比值可見，本文算法ATPFCA相比FCM、PFCA、TPFCA、FCM_s、FLICM、PFCA_s和文獻(xiàn)[11]算法，能獲得較小的誤分率和較大的峰值信噪比，進(jìn)一步表明本文算法ATPFCA的聚類分割性能更好、抑制高斯噪聲干擾的能力更強(qiáng)。

5.3 椒鹽噪聲干擾圖像的測(cè)試與分析

對(duì)圖2所示的3幅圖像分別添加概率為30%、10%和20%的椒鹽噪聲，采用FCM、PFCA、TPFCA、FCM_s、FLICM、PFCA_s、文獻(xiàn)[11]算法以及ATPFCA算法進(jìn)行分割測(cè)試，所得測(cè)試結(jié)果如表6和表7所示。

表6 椒鹽噪聲干擾圖像及分割結(jié)果

Tab.6 Salt and pepper noise images and their segmentation results

從表6可見，F(xiàn)CM_s、FLICM、PFCA_s、文獻(xiàn)[11]算法和ATPFCA算法因加入了鄰域信息，使得其分割結(jié)果中的噪聲相對(duì)FCM、PFCA、TPFCA算法所得結(jié)果明顯減少，其中ATPFCA算法相比其他4種加入鄰域信息的算法具有更強(qiáng)的抗噪魯棒性。從圖4中工業(yè)遙感局部區(qū)域放大圖來看，F(xiàn)CM、PFCA和TPFCA算法所得結(jié)果含有大量的噪聲，F(xiàn)CM_s、FLICM、PFCA_s和文獻(xiàn)[11]算法所得結(jié)果噪聲相對(duì)較少，ATPFCA算法所得結(jié)果噪聲近乎為0，且獲取的目標(biāo)輪廓相對(duì)其他7種算法也更為清晰。此外，從表7可見，嵌入鄰域信息的FCM_s、PFCA_s、FLICM、文獻(xiàn)[11]算法和ATPFCA算法的峰值信噪比(誤分率)均大于(小于)FCM、PFCA和TPFCA算法，而ATPFCA算法的峰值信噪比(誤分率)相對(duì)FCM_s、PFCA_s、文獻(xiàn)[11]算法和FLICM算法更高(低)。上述結(jié)果進(jìn)一步表明，本文算法ATPFCA具有好的分割性能且抗椒鹽噪聲魯棒性更強(qiáng)。

圖4 遙感圖局部區(qū)域及分割結(jié)果Fig.4 Local region of SAR image and its segmentation results

6 結(jié)論

本文引入全散度距離度量，與像素鄰域信息相結(jié)合，得到改進(jìn)的全散度，并將像素空間鄰域信息嵌入算法目標(biāo)函數(shù)，提出了一種ATPFCA圖像分割算法。得出主要結(jié)論如下：

1)TPFCA算法能夠?qū)Σ煌愋突叶确植疾痪鶆驁D像進(jìn)行有效分割，其分割結(jié)果的誤分率(峰值信噪比)均比傳統(tǒng)FCM、PFCA減少約30%～55%(增加約19%～33%)。全散度的引入在一定程度上提高了算法對(duì)不均勻圖像的有效分割程度。

2)ATPFCA算法相比傳統(tǒng)FCM、PFCA、FCM_s、FLICM以及PFCA_s算法具有更強(qiáng)的分割灰度不均勻圖像的能力，其分割結(jié)果的誤分率(峰值信噪比)相比上述聚類算法減少約38%～76%(增加約21%～67%)，與TPFCA算法相比，對(duì)不均勻圖像的有效分割程度得到進(jìn)一步提升。針對(duì)噪聲干擾圖像，ATPFCA算法對(duì)應(yīng)分割結(jié)果的誤分率(峰值信噪比)相比FCM、PFCA、TPFCA、FCM_s、FLICM、文獻(xiàn)[11]算法和PFCA_s算法等減少約5%～77%(增加約6%～230%，特別是相對(duì)于FCM增加了高達(dá)230%)。自適應(yīng)像素空間鄰域信息的嵌入，不僅提升了算法對(duì)灰度不均勻圖像的有效分割程度，而且增強(qiáng)了算法的抗噪魯棒性。

綜上所述，本文算法能夠?qū)叶确植疾痪鶆驁D像得到滿意的分割結(jié)果，且具有良好的抑制噪聲干擾能力。