改進(jìn)的基于模糊C_均值聚類的圖像分割算法

1、第31卷第3期西南大學(xué)學(xué)報(自然科學(xué)版2009年3月Vol131No13Journal of Sout hwest University(Nat ural Science EditionMar12009文章編號:167329868(20090320169204改進(jìn)的基于模糊C2均值聚類的圖像分割算法王志兵,魯瑞華西南大學(xué)電子信息工程學(xué)院,重慶400715摘要:為了提高圖像分割算法的抗噪聲性能,提出了一種改進(jìn)的基于模糊C2均值聚類的圖像分割算法.該算法首先根據(jù)鄰域像素的隸屬度矩陣來計算出像素和聚類中心的空間距離,然后利用空間距離和歐氏距離來重新確定像素和聚類中心的距離,最后利用新提取的距離特征和

2、改進(jìn)的FCM聚類算法對圖像進(jìn)行分割.實驗結(jié)果表明,該算法能有效地提取目標(biāo)圖像,對噪聲具有較強(qiáng)的魯棒性,收斂速度快.關(guān)鍵詞:空間距離;歐氏距離;模糊C2均值聚類;圖像分割中圖分類號:TP391文獻(xiàn)標(biāo)識碼:A圖像分割的瓶頸問題是把圖像分成各具特性的區(qū)域并提出感興趣的目標(biāo).至今已提出多種不同的圖像分割方法,如直方圖閾值法、特征空間聚類法、基于區(qū)域的方法、模糊方法等122,其中模糊分割算法因為能比傳統(tǒng)的硬分割算法保留更多的原始圖像信息而引起了人們的關(guān)注3,特別是模糊c2均值(Fuzzy c2 mean,FCM算法3因描述簡潔,實現(xiàn)方便,已在圖像分析、醫(yī)療診斷、目標(biāo)識別和圖像分割等領(lǐng)域得到廣泛應(yīng)用3,5

3、26.傳統(tǒng)的FCM算法在特征提取中未考慮其空間信息,對于分割含噪聲的圖像不能得到滿意的結(jié)果.為了在圖像分割的過程中能有效抑制噪聲的影響,很多人對傳統(tǒng)的FCM算法做了改進(jìn)729,其中文獻(xiàn)9的基于空間模式聚類的圖像分割方法能夠有效地抑制噪聲,但對特征集的計算比較復(fù)雜,而且分割效果受權(quán)重因子的影響很大.因此本文提出了一種改進(jìn)的基于模糊c2均值聚類的圖像分割算法,根據(jù)鄰域像素的隸屬度矩陣來計算出像素和聚類中心的空間距離,利用空間距離和歐氏距離新確定的距離特征和改進(jìn)的FCM聚類算法對圖像進(jìn)行分割.實驗結(jié)果表明,該方法對特征集的計算比較簡單,受權(quán)重因子影響小,能夠有效地提取目標(biāo)圖像,具有較強(qiáng)的噪聲抑制能力

4、.1改進(jìn)的模糊C2均值聚類算法111基于空間信息的特征提取d(x k,v i=(d E ik/N+d S ik2(1收稿日期:2008207202基金項目:西南師范大學(xué)發(fā)展基金資助項目(SWNU F2004006;重慶市自然科學(xué)基金資助項目(2007BB2331.作者簡介:王志兵(19842,男,江蘇贛州人,碩士研究生,主要從事信號與信息處理的研究.通訊作者:魯瑞華,教授.其中N =256為灰度級,d E ik 為像素x k 和聚類中心v i 之間的歐氏距離,d S ik 為像素x k 和聚類中心v i 之間的空間距離.像素x k 和聚類中心v i 之間的歐氏距離d E ik 定義為:d E

5、ik =x k -v i 2(2其中2表示2范數(shù).像素x k 和聚類中心v i 之間的空間距離主要由其鄰域像素的隸屬度來確定,本文在計算像素和聚類中心的空間距離時,首先確定圖像中每個像素對聚類中心的隸屬度,然后根據(jù)隸屬度矩陣計算出像素x k 隸屬于聚類中心v i 的概率p S ik ,最后由概率p S ik 來確定像素x k 和聚類中心v i 之間的空間距離d S ik .對于任一像素點,p S ik 越大,像素x k 和聚類中心v i 的空間距離d S ik 就越小;反之,像素x k 和聚類中心v i 的空間距離d S ik 就越大.因此像素x k 和聚類中心v i 之間的空間距離d S i

6、k 計算如下:d S ik =1-p S ik (3對于任何一個像素x k ,如果它的所有鄰域像素對于聚類中心v i 的隸屬度之和越大,那么像素x k 隸屬于聚類中心v i 的概率p S ik 就越大.像素x k 在空間上隸屬于聚類中心v i 的概率可以描述為:如果像素x k 的所有鄰域像素都以概率1隸屬于聚類中心v i ,則p S ik =1;反之,如果該像素的所有鄰域都以概率0隸屬于聚類中心v i ,則p S ik =0.因此像素x k 在空間隸屬于聚類中心v i 的概率p S ik 定義為:p S ik =t k u itc c =1t k u it (4其中,k 為像素x k 的所有鄰

7、域的集合,u it 為像素x k 的某個鄰域x t 對于聚類中心v i 的隸屬度,0u it 1.112本文算法的實現(xiàn)模糊C 均值聚類算法,是用隸屬度確定每個數(shù)據(jù)點屬于某個聚類的程度的一種聚類算法3.根據(jù)式(1得出像素x k 到聚類中心v i 的距離d (x k ,v i ,則改進(jìn)的FCM 聚類算法的目標(biāo)函數(shù)為:Min J m (U ,V =n k =1ci =1(u ik m d 2(x k ,v i (5約束條件為:ck =1u ik =1,1k n (6u ik 0,1,1k n ,1i c(7這里,X =x 1,x 2,x n R p 為聚類樣本集合,n 是聚類空間的樣本數(shù),c 是聚類

8、的類別數(shù),2c n ,u ik 是第i 類中像素x k 的隸屬度,v i 是第i 類的聚類中心,d (x k ,v i 為像素x k 到聚類中心v i 的距離,U =u ik 是模糊分類矩陣,V =v 1,v 2,v c 為中心矩陣,m 1,是一個加權(quán)指數(shù).根據(jù)新的距離特征可以得出改進(jìn)的模糊隸屬度函數(shù):u ik =d (x k ,v i -2m-1cj =1d (x k ,v j -2m-1(8聚類中心的迭代函數(shù)為:v i =nk =1(uik m x k n k =1(u ik m(9本文算法的具體實現(xiàn)步驟為:Step1:設(shè)定閾值,聚類的類別數(shù)c 和加權(quán)指數(shù)m ,選擇初始聚類中心V 0.St

9、ep2:根據(jù)式(2計算像素x k 和聚類中心v i 之間的歐氏距離d E ik ,根據(jù)式(4計算像素x k 在空間上隸屬于聚類中心v i 的概率p S ik ,然后把p S ik 代入式(3計算像素x k 和聚類中心v i 之間的空間距離d S ik .Step3:根據(jù)式(1計算像素x k 到聚類中心v i 的距離d (x k ,v i .071西南大學(xué)學(xué)報(自然科學(xué)版投稿網(wǎng)址http :/xbgjxt 1swu 1cn 第31卷Step4:根據(jù)式(8更新模糊隸屬度矩陣u ik .Step5:根據(jù)式(9更新聚類中心v i .Step6:若V new -V old <,則算法停止;否則,返

10、回Step4.Step7:根據(jù)隸屬度矩陣,利用最大隸屬度原則,將軟分割結(jié)果轉(zhuǎn)化為硬分割結(jié)果,輸出分割結(jié)果.2實驗結(jié)果與分析為了驗證本文算法的有效性,選取傳統(tǒng)FCM 算法和本文算法進(jìn)行實驗.設(shè)算法中的參數(shù)=1e -6,m =2.圖1(a 為一幅128×128的圖像,由目標(biāo)T 和背景兩部分組成,灰度值分別對應(yīng)為255和0,故分割時類別數(shù)c =2,圖1(b 為添加5%高斯噪聲的圖像.利用傳統(tǒng)FCM 算法和本文算法分別對圖1(b 進(jìn)行分割可以看出,傳統(tǒng)FCM 算法僅根據(jù)當(dāng)前像素進(jìn)行分割,不能有效地抑制圖像中存在的噪聲.而本文算法由于考慮了像素點的空間信息,能很好地提取目標(biāo)T ,對噪聲的抑制能

11、力要比傳統(tǒng)FCM 算法強(qiáng) .圖1傳統(tǒng)FCM 算法和本文算法對噪聲圖像(b的分割結(jié)果圖2(a 是一幅含有多類目標(biāo)的顱腦圖像(MRI ,圖2(b 是添加9%高斯噪聲的圖像.對顱腦圖像分割時,通常將腦部組織分成白質(zhì)、灰質(zhì)、腦脊髓和背景4部分,所以算法中類別數(shù)c =4,對應(yīng)的灰度值分別為240,160,80,0.比較圖2中的(c 和(d 可以看出,傳統(tǒng)FCM 算法不能成功地將圖像分成四類,一些白質(zhì)、灰質(zhì)和背景被錯誤分割.而本文算法利用空間距離降低了噪聲隸屬度的權(quán)重,不僅能成功地分割出白質(zhì)、灰質(zhì)、腦脊髓和背景4類,而且去除了白質(zhì)和灰質(zhì)部分的許多誤分類點,白質(zhì)和灰質(zhì)的邊緣也變得更加清晰,背景中的噪聲也被去

12、除了 .圖2傳統(tǒng)FCM 算法和本文算法對噪聲圖像(b的分割結(jié)果為了更好地比較兩種算法的性能,本文分別對噪聲含量為1%,3%,5%,7%,9%的圖像進(jìn)行測試,表1列出了在不同噪聲影響下兩種方法的正確分割率(SA 11,SA 定義為:SA =正確分割的像素數(shù)所有的像素數(shù)×100%表1兩種算法對不同噪聲圖像的SA圖像算法1%噪聲3%噪聲5%噪聲7%噪聲9%噪聲圖1(b 從表1中可以看出,對于類別比較少的圖像,當(dāng)噪聲含量比較低的時候,兩種方法都能得到較高的SA ,隨著噪聲含量的增加,傳統(tǒng)FCM 算法的SA 迅速下降,而本文算法的SA 變化不大.當(dāng)對類別比較多,信噪比較低的顱腦圖像進(jìn)行分割時,

13、在不同噪聲影響下,本文算法的SA 都比傳統(tǒng)FCM 算法的SA 要高,表明在分割噪聲圖像時本文算法比傳統(tǒng)的FCM 算法抑制噪聲的能力強(qiáng),具有較強(qiáng)的魯棒性.171第3期王志兵,等:改進(jìn)的基于模糊C 2均值聚類的圖像分割算法271西南大學(xué)學(xué)報(自然科學(xué)版投稿網(wǎng)址http:/xbgjxt1swu1cn第31卷3結(jié)束語本文提出了一種改進(jìn)的模糊c2均值聚類算法,在計算像素和聚類中心的距離時考慮了像素的空間信息,通過像素在空間上的隸屬度來確定其空間距離,利用空間距離和歐氏距離新確定的距離特征和改進(jìn)的FCM聚類算法對圖像進(jìn)行分割,抑制了噪聲對聚類的影響,有效地糾正了被誤分的像素點,分割效果好,收斂速度快.參考

14、文獻(xiàn):1Cheng H D,Jiang X H,Sun Y.Color Image Segmentation:Advances and ProspectsJ.Pattern Recognition,2001,34:225922281.2陳婷婷,程小平.采用模糊形態(tài)學(xué)和形態(tài)學(xué)分水嶺算法的圖像分割J.西南大學(xué)學(xué)報,2008,30(3:1422145.3楊勇,鄭崇勛,林盤.基于改進(jìn)的模糊C均值聚類圖像分割新算法J.光電子激光,2005,16(9:111821122.4王力新.模糊系統(tǒng)與模糊控制教程M.王迎軍,譯.北京:清華大學(xué)出版社,2003:2762280.5田捷,韓博聞,王巖.模糊C2均值聚類法

15、在醫(yī)學(xué)圖像分析中的應(yīng)用J.軟件學(xué)報,2001,12(11:162321629.6李云松,李明.基于灰度空間特征的模糊C均值聚類圖像分割J.計算機(jī)工程與設(shè)計,2007,28(6:135821363.7劉華軍,任明武,楊靜宇.一種改進(jìn)的基于模糊聚類的圖像分割方法J.中國圖像圖形學(xué)報,2006,11(9:131221316.8Cai W L,Chen S C,Zhang D Q.Fast and Robust Fuzzy c2Means Clustering Algorithms Incorporating Local Informationfor Image SegmentationJ.Patte

16、rn Recognition,2007,40:8252838.9Xia Y,Feng D G,Wang T J.Image Segmentation by Clustering of Spatial PatternsJ.Pattern Recognition Letters,2007,28:154821555.10楊淑瑩.圖像模式識別VC+技術(shù)實現(xiàn)M.北京:清華大學(xué)出版社;北京交通大學(xué)出版社,2005:1612162.11Ahmed M N,Yamany S M,Mohamed N.A Modified Fuzzy c2Means Algorithm for Bias Field Estima

17、tion and Segmenta2tion of MRI DataJ.IEEE Trans on Medical Imaging,2002,21(3:1932199.An Improved Algorithm for Im age SegmentationB ased on Fuzzy C2Means ClusteringWAN G Zhi2bing,L U Rui2huaSchool of E lectronics and Information Engineering,Southwest University,Chongqing400715,ChinaAbstract:To enhanc

18、e t he noise immunity performance of t he image segmentation algorit hm,an improved algorit hm for image segmentation based on f uzzy C2means clustering is p roposed in t his paper.The spatial distance between a pixel and t he cluster center is calculated by t he membership mat rix of t he neighboring pixels,and a new distance is determined by t he spatial distance and t he Euclidean distance.This new dis2 tance feat ure and t he improved algorit hm based on f uzzy C2means clustering are used in image segmenta2 tion.T