廖林峰，邱曉暉

(南京郵電大學 通信與信息工程學院，江蘇 南京 210003)

基于模糊C-均值聚類醫(yī)學圖像分割的優(yōu)化算法

廖林峰，邱曉暉

(南京郵電大學 通信與信息工程學院，江蘇 南京 210003)

模糊C-均值聚類(FCM)算法在分割模糊的醫(yī)學圖像中有很好的效果，通過設(shè)置初始聚類中心，根據(jù)每個像素的隸屬度來劃分屬于哪一類，采用迭代的方式來得到分割結(jié)果。針對FCM算法容易受到聚類中心初始值和噪聲的影響，采用遺傳算法和粒子群算法的結(jié)合算法來確定一組合適的初始聚類中心，通過遺傳算法和粒子群算法的結(jié)合算法加快了單純使用遺傳算法確定初始聚類中心的收斂速度；再通過引入像素的鄰域信息，重構(gòu)標準FCM算法中的目標函數(shù)，以提高鄰域像素和中心像素之間的相似程度，使得相鄰的像素更容易劃分到同一類別，克服了標準FCM算法只考慮像素間的灰度值而導致對噪聲和異常值的敏感問題。將該方法應(yīng)用到核磁共振成像(MRI)腦部圖像分割實驗中，相比標準的FCM分割算法和遺傳模糊聚類算法，分割效果更好。

模糊C-均值聚類；遺傳算法；粒子群算法；鄰域像素；核磁共振成像

0 引 言

由于醫(yī)學圖像含有的噪聲高，而且各種組織之間邊界比較模糊，且組織內(nèi)部的灰度很不均勻，對比度也相對較低，所以采用傳統(tǒng)的分割技術(shù)很難達到預(yù)期的效果[1]。針對這些問題，將FCM算法(Fuzzy C-Means)廣泛應(yīng)用在這類圖像的分割中，國內(nèi)外專家也對FCM算法應(yīng)用在醫(yī)學圖像分割領(lǐng)域做了大量研究，

并且取得了一定的基礎(chǔ)研究成果，為后續(xù)研究提供了一定的參考。

FCM算法是1974年由Dunn[2]在硬C均值算法的基礎(chǔ)上提出的，轟動一時，很多專家都參與了研究。在同一年，BEZEDEK[3]又進行了一般化的推廣，建立了一套比較完善的FCM算法的基礎(chǔ)理論。1980年，BEZEDEK進一步對FCM算法的收斂性進行了證明[4]。由于FCM算法對初始聚類中心比較敏感，所以文獻[5]提出將遺傳算法(GA)應(yīng)用于FCM算法中，通過GA優(yōu)秀的全局尋優(yōu)能力來確定FCM算法中的聚類數(shù)目并尋找到一組全局最優(yōu)的聚類中心，也就是GA-FCM算法[5]。鑒于其有較好的分割效果，GA-FCM算法應(yīng)用范圍很廣泛[6-8]。

由于FCM算法的目標函數(shù)僅僅考慮了像素的灰度信息，所以對異常值和噪聲比較敏感。劉曉妮等[9]采用FCM算法和離散正則化算法相結(jié)合的分割方法，首先用標準FCM算法對目標圖像進行初始分割,再將得到的分割圖像采用半監(jiān)督的離散正則化算法進行二次分割，但是算法的計算復(fù)雜度較大，分割效率也變低了。文獻[10-13]也對模糊聚類使用遺傳算法進行改進，但算法復(fù)雜度太高。

針對上述問題，在尋找最優(yōu)初始聚類中心時，文中在傳統(tǒng)的粒子群算法中同時加入遺傳算法中的兩個算子，交叉算子和變異算子，使得兩個算法的優(yōu)勢互補，克服粒子群算法容易陷入局部最優(yōu)的缺陷，為后續(xù)分割提供更精確的初始聚類中心；同時引入像素的標號場和特征場信息，重新構(gòu)造傳統(tǒng)FCM算法的目標函數(shù)，加強鄰域像素和聚類中心像素間的作用強度，克服僅僅考慮像素的灰度信息而導致對噪聲和異常值敏感的問題，從而得到更精準的分割結(jié)果。

1 標準FCM算法

FCM算法是一種模糊、無監(jiān)督的聚類算法，它的分類是柔性的，而且實現(xiàn)簡單，所以應(yīng)用廣泛。簡言之，F(xiàn)CM算法將像素的灰度值作為一個參數(shù)，把每個像素到聚類中心像素的歐氏距離作為相似性測度，然后使得歐氏距離的加權(quán)取最小值，從而進行圖像數(shù)據(jù)的聚類。FCM算法根據(jù)歐氏距離的加權(quán)方式，把目標函數(shù)定義為:

(1)

目標函數(shù)J的物理意義表示目標圖像中的每個像素到各聚類中心加權(quán)距離的平方和。當目標圖像中的每個像素點到某個聚類中心的歐氏距離加權(quán)值最小，而距離其他聚類中心的歐氏距離盡量大時,FCM算法的基本原則是尋找到一組合適的聚類中心和隸屬度矩陣，使得目標函數(shù)J取得最小值min(J)。

通過求目標函數(shù)的最小值可以得到uij(i=1,2,…,C),uij是第j個像素隸屬于不同聚類的隸屬度；最后將每個像素劃分到不同的類中，把第j個像素劃分到第k類的劃分公式如下[14]：

k=argmax{uij,i=1,2,…,C}

(2)

(3)

(4)

(5)

在計算目標函數(shù)的最小值時，需要根據(jù)式(4)和式(5)不斷更新隸屬度矩陣和聚類中心，直到求得最小值為止。在目標函數(shù)取得最小值時，保留此時的隸屬度，根據(jù)式(2)對目標圖像的每個像素進行分類，從而得到圖像的分割結(jié)果。在迭代過程中要選取合適的迭代終止條件，這里不能根據(jù)目標函數(shù)前后兩次迭代的差值小于某個閾值就停止迭代，這樣使得算法很容易陷入局部最優(yōu)，得不到理想的分割結(jié)果。標準的FCM算法一般是根據(jù)前后兩次迭代的像素的隸屬度的差值小于某個設(shè)定的閾值就停止迭代，因為隸屬度更新和聚類中心有關(guān)，當初始聚類中心選取不恰當時算法容易陷入局部最優(yōu)；也可以設(shè)置一個最大迭代次數(shù)，到達迭代次數(shù)時立即停止迭代。

通過上述的推導分析，可以將FCM算法分割圖像的過程總結(jié)為：

(2)根據(jù)式(5)計算像素的聚類中心vi(i=1,2,…,C)；

(5)根據(jù)式(2)對圖像中的每個像素進行分類，輸出分類后的圖像。

從式(1)可知,要使目標函數(shù)取得最小值，就需要保證目標圖像中的像素離它最近的聚類中心的隸屬度盡可能大，而屬于其他聚類中心的隸屬度盡可能小。

2 改進算法

在使用FCM算法進行分割前，利用遺傳-粒子群算法尋找一組最優(yōu)的初始聚類中心，在粒子群算法中引入遺傳算法的交叉算子和變異算子。交叉算子充當?shù)淖饔檬鞘钩蓪Φ牧Ｗ又g可以互相交流信息，提高粒子往新的方向搜索的能力，而變異算子是用來提高粒子群算法跳出局部最優(yōu)的能力，并保留粒子群算法中的種群分割策略來保持種群的多樣性，克服粒子群容易陷入局部最優(yōu)的缺陷。其中將FCM算法的目標函數(shù)的倒數(shù)作為遺傳-粒子群算法的適應(yīng)度函數(shù)的一個參數(shù)，適應(yīng)度函數(shù)為：

(6)

其中，J為標準FCM算法中的目標函數(shù)。由此可知，當遺傳-粒子群算法中的適應(yīng)度函數(shù)F的取值越大時，聚類效果就越好[16]。

通過引入鄰域像素標號場和特征場，重構(gòu)目標函數(shù)，加強鄰域像素和聚類中心像素之間的作用強度，有利于削弱異常值和噪聲對分割結(jié)果的影響，得到更準確的分類結(jié)果。文中通過馬爾可夫理論引入鄰域像素的標號場信息。從理論上講,相鄰像素之間有很大的可能具有相同的標號,用Nj表示第j個像素的鄰域,則像素j劃分為聚類i的先驗概率可以表示為[17]：

(7)

其中，lj表示第j個像素的標號；β表示鄰域像素的作用強度；W(β)表示將先驗概率轉(zhuǎn)化到0到1之間的歸一化項；U(lj=i|β)表示能量函數(shù)，為方便，下文將用πij表示。

為了更好地表示第j像素與聚類i之間的關(guān)系,采用式(8)作為新的非相似性測度：

dij=-logp(xj|lj=i)

(8)

根據(jù)上述的分析和推導，可以將改進算法的目標函數(shù)定義為：

(9)

(10)

(11)

(12)

其中，eij定義為：

(13)

其中，γ表示鄰域像素的作用強度；t(li=lj)表示一個指數(shù)函數(shù)。

綜上，文中改進算法的步驟為：

(1)設(shè)置算法中的初始值，包括聚類中心數(shù)目C、模糊因子m、先驗概率中的β、最大迭代次數(shù)K、特征場作用強度γ和迭代停止條件θ;

(2)根據(jù)GA-PSO確定一組最優(yōu)的初始聚類中心，并根據(jù)式(4)初始化隸屬度矩陣U(0)；

(3)根據(jù)式(7)、(11)、(12)分別計算每個聚類中的先驗概率、均值u(k)和協(xié)方差∑(k);

(4)根據(jù)式(10)更新隸屬度矩陣U(k)；

(5)如果{U(k+1)-U(k)}<θ或者到了最大迭代次數(shù)，則停止，否則令k=k+1，轉(zhuǎn)步驟(3)繼續(xù)進行迭代運算。

3 實驗結(jié)果及分析

將改進算法應(yīng)用到比較有代表性的MRI腦部圖像分割實驗中，目標圖像是一幅大小為256×256的灰度圖像。因為人腦的MRI圖像主要包括腦脊液、腦白質(zhì)、腦灰質(zhì)和背景四個部分，所以將聚類數(shù)目取4，模糊因子取2，分別用標準的FCM算法、遺傳模糊聚類算法和文中算法對其進行分割，結(jié)果如圖1所示。

從圖1可以看出，傳統(tǒng)FCM算法和遺傳模糊聚類算法受噪聲影響，分割結(jié)果不是很理想，有很多細小和模糊的部分，存在明顯的噪聲干擾；而遺傳模糊算法相比傳統(tǒng)算法分割效果較好；文中算法因為加入了鄰域像素作用，使得算法抗噪性能得到了較大提升，分割結(jié)果也相對更加精確。

分割參數(shù)如表1所示。其中Vpc為BEZDEK劃分系數(shù)，是一個衡量FCM算法聚類效果的參數(shù)，主要目的是對分割后像素的隸屬度進行刻畫。該指標的定義如下：

圖1 各算法分割結(jié)果

(14)

表1 分割參數(shù)對比

一般而言，對于目標圖像中的每個像素點來說，一個聚類效果好的算法會使它屬于其中一類的隸屬度盡可能大而屬于其他類的隸屬度盡可能小。因此，一個聚類算法的Vpc越大，說明這個聚類算法的效果越好。表1中，傳統(tǒng)的FCM算法的Vpc比遺傳聚類算法和文中算法都要小，而遺傳聚類算法和文中算法的Vpc相近，但是文中算法的迭代次數(shù)相比遺傳聚類算法少了很多，比傳統(tǒng)FCM算法的迭代次數(shù)也相對較少。因此，文中算法可以有效解決傳統(tǒng)FCM算法對聚類中心初始值、圖像中的噪聲以及異常值敏感的缺陷，分割結(jié)果改善明顯。

4 結(jié)束語

在標準的FCM算法及其改進算法的基礎(chǔ)上，通過結(jié)合遺傳算法和粒子群算法來尋找最優(yōu)的初始聚類中心，減少了算法的迭代次數(shù)；同時引入了特征場和標號場，重新構(gòu)造FCM算法中的目標函數(shù)對目標圖像進行分割，使得算法的抗噪性能得到了一定的提升。但是文中算法還有一些參數(shù)需要依賴人工經(jīng)驗選取，比如模糊因子、標號場以及特征場作用強度參數(shù)，下一步的研究中將尋找一種合適的方法來自動選取最優(yōu)參數(shù)，從而將圖像的半自動化分割轉(zhuǎn)為完全自動化分割。

[1] 張 翡，范 虹.基于模糊C均值聚類的醫(yī)學圖像分割研究[J].計算機工程與應(yīng)用，2014，50(4)：144-151.

[2] DUNN J C.A fuzzy relative of the ISODATA process and its use in detecting compact well-separated clusters[J].Journal of Cybernetics，1973，3(3)：32-57.

[3] BEZEDEK J C.Pattern recognition with fuzzy objective function algorithm[M].New York：Plenum Press，1981.

[4] CANNON R L,DAVE J V,BEZDEK J C.Efficient implementation of the fuzzy C-means clustering algorithms[J].IEEE Transactions on Pattern Analysis and Machine Intelligence，1986，8(2)：248-255.

[5] VIJAYACHITR S,TAMILARSI A,KASTHURI N.Multiple input single output (MISO) process optimization using[C]//International conference on education technology and computer.Singapore:IEEE,2009:248-252.

[6] LIU Suhua,HOU Huifang.A combination of mixture genetic algorithm and fuzzy c-means clustering[C]//International symposium on IT in medicine & education.Ji’nan，China:IEEE，2009:254-258.

[7] BELAHBIB F Z B,SOUAMI F.Genetic algorithm clustering for color image quantization[C]//3rd European workshop on visual information processing.Paris,French:[s.n.],2011:83-87.

[8] BEZDEK J C, EHRLICH R, FULL W. FCM:the fuzzy c-means clustering algorithm[J].Computer and Geosciences，1984,10(2):191-203.

[9] 劉曉妮,盧奕南,雷 玲.基于FCM和離散正則化的多目標圖像分割[J].計算機輔助設(shè)計與圖形學學報,2015,27(1):142-146.

[10] 張永庫，尹靈雪，孫勁光.基于改進的遺傳算法的模糊聚類算法[J].智能系統(tǒng)學報，2015,10(4)：627-635.

[11] 楊 凱，蔣華偉.模糊C均值聚類圖像分割的改進遺傳算法研究[J].計算機工程與應(yīng)用，2009,45(33)：179-182.

[12] 邱雙雙.基于核模糊C-均值聚類與閾值分割的SAR影像分割算法[J].科技創(chuàng)新與應(yīng)用，2014(35)：15.

[13] 張佳骕，蔣亦樟，王士同.基于特征選擇聚類方法的稀疏TSK模糊系統(tǒng)[J].智能系統(tǒng)學報，2015,10(4)：583-591.

[14] 張自嘉，岳邦珊，潘 琦，等.基于蟻群和自適應(yīng)濾波的模糊聚類圖像分割[J].電子技術(shù)應(yīng)用，2015,41(4)：144-147.

[15] 宋孌孌.一種基于圖像濾波的加權(quán)FCM圖像分割算法[J].商丘師范學院學報，2014,30(12)：10-14.

[16] 張紅旗，王春光，李海軍.基于遺傳算法的草莓圖像FCM分割方法研究[J].農(nóng)機化研究,2015(4)：55-57.

[17] 趙雪梅，李 玉，趙泉華.結(jié)合高斯回歸模型和隱馬爾可夫隨機場的模糊聚類圖像分割[J].電子與信息學報，2014，36(11):2730-2736.

AnOptimalAlgorithmforMedicalImageSegmentationBasedonFuzzyC-MeansClustering

LIAO Lin-feng，QIU Xiao-hui

(School of Communication and Information Engineering,Nanjing University of Posts and Telecommunications,Nanjing 210003,China)

Fuzzy C-Means Clustering (FCM) has a good effect in the segmentation of fuzzy medical images.By setting the initial cluster center,the division is carried out according to the membership degree of each pixel,and segmenting results are obtained by means of iteration.Aiming at the problem that FCM is susceptible to initial value of cluster center and noise,the genetic algorithm and particle swarm algorithm are combined to determine a set of suitable initial clustering centers.After combination,the convergence rate of initial clustering center is accelerated than that of only use of genetic algorithm.Then the objective function of standard FCM is reconstructed by introducing neighborhood information of pixels so as to improve the similarity between the neighborhood and the center pixel,which makes the adjacent pixels more divided into the same class and overcomes the problem that the standard FCM only considers the gray value between the pixels and causes the sensitivity to noise and outliers.The proposed method is applied into the MRI brain image segmentation experiment which shows that it is superior to the standard FCM and genetic algorithm on segmentation effect.

Fuzzy C-Means clustering；genetic algorithm；particle swarm optimization；neighborhood pixels；nuclear magnetic resonance imaging

TP301.6

A

1673-629X(2017)12-0081-04

10.3969/j.issn.1673-629X.2017.12.018

2016-12-11

2017-04-13 < class="emphasis_bold">網(wǎng)絡(luò)出版時間

時間：2017-09-27

江蘇省自然科學基金(BK2011789)；東南大學毫米波國家重點實驗室開放課題(K201318)

廖林峰(1992-)，男，碩士研究生，研究方向為圖像處理；邱曉暉，博士，教授，研究方向為圖像處理。

http://kns.cnki.net/kcms/detail/61.1450.TP.20170927.0957.020.html