郭海濤1,2,趙紅葉3,徐　雷3，侯一民3，焦圣喜3

(1.內(nèi)蒙古大學(xué)電子信息工程學(xué)院,內(nèi)蒙古 呼和浩特 010021;2.東北電力大學(xué)電氣工程學(xué)院,

吉林 吉林 132012;3.東北電力大學(xué)自動(dòng)化工程學(xué)院,吉林 吉林 132012)

?

多分辨分析和K均值聚類改進(jìn)FCM圖像分割

郭海濤1,2,趙紅葉3,徐雷3，侯一民3，焦圣喜3

(1.內(nèi)蒙古大學(xué)電子信息工程學(xué)院,內(nèi)蒙古 呼和浩特 010021;2.東北電力大學(xué)電氣工程學(xué)院,

吉林 吉林 132012;3.東北電力大學(xué)自動(dòng)化工程學(xué)院,吉林 吉林 132012)

摘要：模糊 C 均值(Fuzzy C-Means,FCM)聚類廣泛應(yīng)用于圖像分割,但FCM聚類中隨機(jī)確定初始聚類中心容易導(dǎo)致圖像的錯(cuò)誤分割.為了避免這個(gè)缺點(diǎn),提出一種用于圖像分割的FCM聚類初始聚類中心的選取方法.該方法利用圖像灰度-鄰域均值二維直方圖的峰值的個(gè)數(shù)確定圖像聚類數(shù)目，然后對(duì)圖像的低頻子帶圖像利用K均值聚類得到FCM聚類初始聚類中心.實(shí)測(cè)圖像的分割實(shí)驗(yàn)表明該方法具可行性.

關(guān)鍵詞：二維直方圖;多分辨分析;K均值聚類;FCM聚類;圖像分割

模糊C均值(Fuzzy C-Means，F(xiàn)CM)聚類[1]因具有無需訓(xùn)練樣本、無監(jiān)督性、能夠?qū)崿F(xiàn)多區(qū)域分割等特點(diǎn),被廣泛應(yīng)用于圖像分割領(lǐng)域[2-7].FCM聚類中初始聚類中心是隨機(jī)選取的.隨機(jī)選取初始聚類中心有時(shí)會(huì)導(dǎo)致錯(cuò)誤的聚類,因而導(dǎo)致錯(cuò)誤的圖像分割.筆者提出一種用于圖像分割的FCM聚類初始聚類中心的選取方法.該方法的步驟是:(ⅰ)利用圖像的灰度-鄰域均值二維直方圖(簡(jiǎn)稱二維直方圖)[3]的峰值的個(gè)數(shù)確定圖像聚類的數(shù)目(分割的區(qū)域數(shù)目);(ⅱ)對(duì)圖像進(jìn)行小波多分辨分析得到低頻子帶圖像;(ⅲ)對(duì)低頻子帶圖像利用K均值聚類(用步驟(ⅰ)中得到的聚類數(shù)目)得到FCM聚類的初始聚類中心.

利用上述初始聚類中心選取方法獲得圖像的FCM聚類初始聚類中心后,再利用FCM聚類(用步驟(ⅰ)中得到的聚類數(shù)目)重新對(duì)原圖像聚類得到分割后的圖像.下文闡述FCM聚類初始聚類中心的選取方法和完整的FCM聚類分割過程,并以實(shí)際的聲吶圖像為例給出文中方法的實(shí)驗(yàn)結(jié)果,同時(shí)將文中方法與其他常用方法的實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比.

1圖像FCM聚類初始聚類中心選取方法

利用圖像二維直方圖[3]的峰值的個(gè)數(shù)確定圖像聚類數(shù)目，然后對(duì)圖像的低頻子帶圖像利用K均值聚類得到FCM聚類初始聚類中心.下面介紹用于獲得聚類數(shù)目的二維直方圖和用于獲得低頻子帶圖像的小波多分辨分析[8].關(guān)于K均值聚類可參考相關(guān)的文獻(xiàn).

1.1 二維直方圖

以每個(gè)像素周圍的8個(gè)像素作為一個(gè)鄰域,計(jì)算其鄰域灰度均值.設(shè)原始圖像有L個(gè)灰度等級(jí),nij為圖像中灰度值為i而鄰域灰度均值為j的像素點(diǎn)數(shù),M，N分別為圖像的行數(shù)和列數(shù),則pij(i,j=1,2,…,L)就是該圖像的灰度-鄰域均值二維直方圖[3]：

二維直方圖由于同時(shí)考慮了像素灰度值和像素鄰域灰度均值,因此它比灰度直方圖的抗干擾性強(qiáng).在二維直方圖中目標(biāo)和背景更容易區(qū)分[3].

1.2 圖像的小波多分辨分析

對(duì)二維圖像信號(hào)f(x,y)進(jìn)行小波變換后,可以將信號(hào)分解為1個(gè)低頻子帶和3個(gè)高頻子帶,其中低頻子帶包含了圖像的概貌，3個(gè)高頻子帶包含了圖像在水平方向、垂直方向和對(duì)角方向上的細(xì)節(jié)[8].在變換的每一層次,圖像都被分解為4個(gè)1/4大小的圖像,它們都是由原圖像與一個(gè)小波基的內(nèi)積后,再經(jīng)過在行和列方向進(jìn)行2倍的間隔抽樣而生成的.圖像尺寸變?yōu)樯弦粚拥?/4[9].若要對(duì)二維圖像信號(hào)f(x,y)進(jìn)行2層小波分解,則可以在完成1層分解后,在低頻子帶上進(jìn)行小波變換,得到2層分解.

2基于FCM聚類的圖像分割

2.1 FCM聚類算法

FCM聚類算法如下[1,10]：

設(shè)數(shù)據(jù)集X={x1,x2,…,xN}為待分類對(duì)象,c(c>1)為分類數(shù)目,vi為第i(i=1,2,…,c)類類心,uij為第j(j=1,2,…,N)個(gè)對(duì)象屬于第i類的隸屬度.定義FCM聚類的目標(biāo)函數(shù)為[1,10]

目標(biāo)函數(shù)Jm(U,V)的約束條件為

(1)

其中:d2(xj,vi)為xj與vi之間的歐氏距離;m為模糊加權(quán)指數(shù),m>1;V為聚類中心集,V={v1,v2,…,vc};U為模糊矩陣,U=(uij)c×N.

利用拉格朗日乘數(shù)法則,可以求得在(1)式的約束下使目標(biāo)函數(shù)Jm(U,V)取極小值的必要條件為[1,10]

(2)

(3)

目標(biāo)函數(shù)Jm(U,V)的極小值，是根據(jù)二維直方圖確定的聚類數(shù)目c和K均值聚類確定的初始聚類中心通過(2)和(3)式的迭代完成的.迭代時(shí)可以事先指定模糊矩陣U的一個(gè)范數(shù)或最大迭代次數(shù)作為迭代終止條件.具體迭代步驟見參考文獻(xiàn)[1,10].

2.2 FCM聚類圖像分割的步驟

利用FCM聚類分割圖像時(shí),初始聚類中心利用K均值聚類確定,聚類數(shù)目利用圖像二維直方圖[3]的峰值的個(gè)數(shù)確定(即與確定初始聚類中心的K均值聚類的聚類數(shù)目一致).FCM聚類圖像分割步驟如下:

(ⅰ)計(jì)算濾波后圖像的二維直方圖；

(ⅱ)對(duì)二維直方圖進(jìn)行3×3均值濾波；

(ⅲ)由均值濾波后的二維直方圖峰值的個(gè)數(shù)確定圖像K均值聚類和FCM聚類的聚類數(shù)目；

(ⅳ)對(duì)濾波后的圖像進(jìn)行小波多分辨分析,得到低頻子帶圖像；

(ⅴ)利用K均值聚類得到低頻子帶圖像的聚類中心，將其作為FCM聚類的初始聚類中心；

(ⅵ)設(shè)定模糊加權(quán)指數(shù)m；

(ⅶ)根據(jù)(2)式計(jì)算初始隸屬度uij；

(ⅷ)設(shè)定迭代運(yùn)算中止的最大迭代次數(shù)；

(ⅸ)按2.1節(jié)和參考文獻(xiàn)[1,10]中的方法進(jìn)行迭代運(yùn)算(即聚類)；

(ⅹ)去除聚類后圖像中的除孤立區(qū)和孔得到分割后的圖像.

3實(shí)例分析

圖1來源于文獻(xiàn)[11],是一幅海底目標(biāo)的合成孔徑聲吶圖像,圖像尺寸為152120.圖1大部分為海底混響區(qū),該部分構(gòu)成圖像的背景.圖1中較突出的部分為目標(biāo)區(qū).目標(biāo)區(qū)中較亮的部分是目標(biāo)亮區(qū)，與之相鄰(右側(cè))的較暗區(qū)域是目標(biāo)暗區(qū).圖2和3分別為中值濾波(3×3)和均值濾波(3×3)后的圖像.圖4為圖3的二維直方圖,圖4經(jīng)過均值濾波(3×3)后得到圖5.比較圖4和圖5可以看出，利用圖5確定圖像聚類數(shù)目可以避免虛假峰值的干擾.

根據(jù)圖5可知,聚類數(shù)目為3.對(duì)圖3利用db1小波進(jìn)行2層小波分解,得到圖6.對(duì)于圖6,利用K均值聚類得到FCM聚類的初始聚類中心;取m=2[10],最大迭代次數(shù)為100.聚類后得到圖7(未去除孤立區(qū)和孔,下文同).作為與文中方法的比較,這里給出最為常用的多閾值Otsu法[12]的閾值化結(jié)果(圖8)，閾值為50和200.

從圖7和圖8可以看出,雙閾值Otsu法分割出的目標(biāo)亮區(qū)比實(shí)際情況偏小,文中算法分割出的目標(biāo)亮區(qū)和目標(biāo)暗區(qū)的大小、形狀更接近實(shí)際情況,分割結(jié)果更準(zhǔn)確.

圖1　原始圖像　　　　　　　　圖2　圖1中值濾波后　　　　　　圖3　圖2均值濾波后

圖4　圖3的二維直方圖　　　　　　　　　　　　　　　　　圖5　圖4均值濾波后

圖6　圖3的2層小波分解　　　　　　圖7　文中方法聚類　　　　　　　　圖8　雙閾值Otsu法

4結(jié)語

針對(duì)圖像FCM聚類時(shí)隨機(jī)確定初始聚類中心容易導(dǎo)致圖像錯(cuò)誤分割的問題,提出一種用于圖像分割的FCM聚類初始聚類中心的選取方法,并進(jìn)而形成一種圖像FCM聚類分割的改進(jìn)方法.以某海底目標(biāo)的合成孔徑聲吶圖像為例,給出了詳盡的分割過程.實(shí)測(cè)圖像的分割實(shí)驗(yàn)表明文中方法具有可行性.該方法適用于目標(biāo)亮區(qū)、目標(biāo)暗區(qū)、海底混響區(qū)的面積相差不是很懸殊的情況.另外,該方法也存在一些需要進(jìn)一步改進(jìn)的地方,比如:利用二維直方圖確定聚類數(shù)目存在一定的誤差,有時(shí)會(huì)導(dǎo)致錯(cuò)誤的結(jié)果;用于確定FCM聚類初始聚類中心的K均值聚類有時(shí)會(huì)陷入局部最小,這有可能使得其確定的FCM聚類初始聚類中心不是最合適的.

參考文獻(xiàn)：

[1] BEZDEK J C,EHRLICH R,FULL W.FCM:The Fuzzy C-Means Clustering Algorithm[J].Computers & Geosciences,1984,10(2):191-203.

[2] 郭海濤,劉麗媛,趙亞鑫,等.基于MAR與FCM聚類的聲吶圖像分割[J].儀器儀表學(xué)報(bào),2013,34(10):2 322-2 327.

[3] 劉健莊.基于二維直方圖的圖象模糊聚類分割方法[J].電子學(xué)報(bào),1992,20(9):40-46.

[4] THAMARAICHELVI B,YAMUNA G,VANITHA U.An Improved Method of Image Segmentation Using Fuzzy C-Means[C].Proceedings of 2012 IEEE International Conference on Advances in Engineering,Science and Management,Tamil Nadu.New York:IEEE,2012:669-672.

[5] QIU Cunyong,XIAO Jian,YU Long,et al.A Modified Interval Type-2 Fuzzy C-Means Algorithm with Application in MR Image Segmentation[J].Pattern Recognition Letters,2013,34(12):1 329-1 338.

[6] FENG Jie,JIAO L C,ZHANG Xiangrong,et al.Robust Non-Local Fuzzy C-Means Algorithm with Edge Preservation for SAR Image Segmentation[J].Signal Processing,2013,93(2):487-499.

[7] 崔兆華,張萍,李洪軍,等.融合結(jié)構(gòu)特征的增強(qiáng)型FCM圖像分割算法[J].東北大學(xué)學(xué)報(bào):自然科學(xué)版,2013,34(7):922-926.

[8] 楊佐龍,王新超,賈明.基于小波域雙閾值分割的SAR圖像預(yù)處理方法[J].船舶電子對(duì)抗,2013,36(1):62-65.

[9] 袁蔚林,馬燕,劉圣偉,等.基于小波變換及Matlab的遙感圖像壓縮效果分析[J].地理信息世界,2013,20(3):39-44.

[10] 郭桂蓉.模糊模式識(shí)別[M].長(zhǎng)沙:國防科技大學(xué)出版社,1992:173-180.

[11] GROEN J,HANSEN R E,CALLOW H J,et al.Shadow Enhancement in Synthetic Aperture Sonar Using Fixed Focusing[J].IEEE Journal of Oceanic Engineering,2009,34(3):269-284.

[12] 茍中魁,張少軍,李忠富,等.一種Otsu閾值法的推廣——Otsu雙閾值法[J].機(jī)械,2004,31(7):12-14.

(責(zé)任編輯向陽潔)

Improved FCM Image Segmentation Based on Multi-Resolution

Analysis and K-Means Clustering

GUO Haitao1,2,ZHAO Hongye3,XU Lei3,HOU Yimin3,JIAO Shengxi3

(1.College of Electronic Information Engineering,Inner Mongolia University,Hohhot 010021,China;2.College of

Electrical Engineering,Northeast Dianli University,Jilin 132012,Jilin China;3.College of

Automation Engineering,Northeast Dianli University,Jilin 132012,Jilin China)

Abstract:The fuzzy C-Means (FCM) clustering is widely used in image segmentation,but the random determination of initial clustering centers of the FCM clustering is likely to generate incorrect segmentation of an image.To avoid the such deficiency,a method of choosing initial clustering centers in the FCM clustering for image segmentation is proposed.The method determines the number for the image clustering by means of the number of the peaks in the two-dimensional histogram of an image comprised of gray values of pixels and mean values of their neighborhoods.Then the K-means clustering is used to obtain the initial clustering centers of the FCM clustering for the low-frequency subband image of the original image.The image segmentation experiments show that the proposed method is feasible.

Key words:two-dimensional histogram;multi-resolution analysis;K-means clustering;FCM clustering;image segmentation

作者簡(jiǎn)介：郭海濤(1965—)，男，黑龍江安達(dá)人，內(nèi)蒙古大學(xué)電子信息工程學(xué)院教授，博士，主要從事圖像處理、模式識(shí)別等研究.

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目(41076060)；吉林省自然科學(xué)基金資助項(xiàng)目(20130101056JC)；內(nèi)蒙古自然科學(xué)基金資助項(xiàng)目(2014MS0601)；內(nèi)蒙古大學(xué)高層次人才引進(jìn)科研項(xiàng)目(135123)

收稿日期：2014-09-05

中圖分類號(hào)：TP391.4

文獻(xiàn)標(biāo)志碼：A

DOI:10.3969/j.issn.1007-2985.2015.02.007

文章編號(hào)：1007-2985(2015)02-0029-04