張彥清+++侯偉+++李玲+++辛化梅

摘 要：為了得到更好的凝膠圖像分割效果，文章提出了一種基于馬爾科夫隨機(jī)場(chǎng)（MRF）模型分割算法的改進(jìn)算法，首先用非局部均值（NL-means）對(duì)圖像進(jìn)行濾波，然后通過(guò)模糊c均值算法（FCM）實(shí)現(xiàn)了圖像的初始聚類分割，進(jìn)一步通過(guò)二階邏輯模型（MLL）模糊聚類作為先驗(yàn)知識(shí)獲取其先驗(yàn)概率并因此得到后驗(yàn)概率，在此過(guò)程中引入了灰度點(diǎn)密度權(quán)值來(lái)實(shí)現(xiàn)聚類和均值方差的更新，最終實(shí)現(xiàn)凝膠圖像的最優(yōu)化分割。

關(guān)鍵詞：馬爾科夫隨機(jī)場(chǎng)；分割；灰度點(diǎn)密度；凝膠圖像

引言

圖像分割就是把一幅圖像分成若干個(gè)特定的區(qū)域并提出感興趣目標(biāo)的技術(shù)和過(guò)程，也是從圖像處理進(jìn)入圖像分析的關(guān)鍵步驟，屬于基本計(jì)算機(jī)視覺(jué)技術(shù)范疇。只有在圖像分割成功的基礎(chǔ)上才有可能實(shí)現(xiàn)對(duì)目標(biāo)特征的提取和參數(shù)的測(cè)量，才能使更高層圖像的分析和圖像理解成為可能。

20世紀(jì)70年代中期，二維凝膠電泳這個(gè)概念由O'Farrell與Klose等人建立[1]，但那時(shí)并未有蛋白質(zhì)組學(xué)這個(gè)概念，隨著二維凝膠技術(shù)的逐漸成熟發(fā)展，蛋白質(zhì)組學(xué)的發(fā)展也突飛猛進(jìn)，直至今天的輝煌成就。

二維凝膠圖像[2]可以被看作是一幅數(shù)字圖像，在這幅圖像上，蛋白質(zhì)點(diǎn)被掃描成不同形狀、不同大小和不同灰度的點(diǎn)，分割是其分析的關(guān)鍵步驟之一，它的目的是有效的將目標(biāo)蛋白質(zhì)點(diǎn)從凝膠圖像中分割出來(lái)。基于馬爾可夫隨機(jī)場(chǎng)的分割方法是一種利用圖像中像素間的空間相關(guān)性進(jìn)行分割的方法，能夠準(zhǔn)確地描述每個(gè)像素所屬類別與周圍像素類別之間的重要依賴關(guān)系。國(guó)內(nèi)外的研究者們已對(duì)MRF的圖像分割方法進(jìn)行了深入了研究，其最早是由S. German和D. German[3]提到的有關(guān)隨機(jī)場(chǎng)概念的問(wèn)題開(kāi)始的，而后在1994年，Bouman等再次[4]提出了具有因果特性的馬爾科夫先驗(yàn)?zāi)Ｐ?，文獻(xiàn)[5]利用馬爾科夫隨機(jī)場(chǎng)的MLL模型來(lái)描敘像素標(biāo)號(hào)場(chǎng)的先驗(yàn)概率，文獻(xiàn)[6]提出的HMRF-FCM圖像分割算法實(shí)現(xiàn)了用均值場(chǎng)來(lái)逼近馬爾科夫隨機(jī)場(chǎng)，但是他的算法復(fù)雜度較高，計(jì)算效率較低。

本文在綜合分析已有文獻(xiàn)算法的基礎(chǔ)上，結(jié)合凝膠圖像的特點(diǎn)，提出了一種改進(jìn)的MRF凝膠圖像分割方法。改進(jìn)算法首先利用非局部均值算法對(duì)圖像進(jìn)行預(yù)處理，根據(jù)像素點(diǎn)之間的相似度進(jìn)行去噪，使其在結(jié)構(gòu)和細(xì)節(jié)保護(hù)上具有一定優(yōu)勢(shì)，同時(shí)對(duì)各類灰度均值點(diǎn)密度聚類采用歸一化加權(quán)方式引入權(quán)重ω，使其后驗(yàn)概率的確定更加準(zhǔn)確從而實(shí)現(xiàn)凝膠圖像的分割。

1 算法介紹

1.1 馬爾可夫隨機(jī)場(chǎng)

馬爾可夫隨機(jī)場(chǎng)方法[7]是建立在MRF模型與貝葉斯（Bayes）理論的基礎(chǔ)上，提供了先驗(yàn)知識(shí)與不確定性描述聯(lián)系的紐帶，并根據(jù)觀測(cè)圖像的特點(diǎn)，利用估計(jì)理論和統(tǒng)計(jì)決策中某些最優(yōu)準(zhǔn)則確定分割問(wèn)題的目標(biāo)函數(shù)，以此來(lái)求解滿足這些條件的最大可能分布，從而將分割問(wèn)題轉(zhuǎn)化為最優(yōu)化問(wèn)題?；贛RF 的分割方法是一種利用圖像像素的空間相關(guān)性進(jìn)行分割的方法，它先將圖像的灰度信息進(jìn)行分類，并且能夠準(zhǔn)確的描述每個(gè)像素所屬分類與周圍像素類別之間的相互依賴關(guān)系。

1.2 基于改進(jìn)MRF的凝膠圖像分割問(wèn)題實(shí)現(xiàn)

在馬爾科夫隨機(jī)場(chǎng)中，常用兩個(gè)隨機(jī)場(chǎng)來(lái)描述帶分割的圖像，假設(shè)X和Y是二維平面上的隨機(jī)場(chǎng)，則令X是像素的標(biāo)號(hào)場(chǎng)，一般用先驗(yàn)分布描述標(biāo)號(hào)場(chǎng)的局部相關(guān)性；令Y為觀測(cè)圖像場(chǎng)，常以標(biāo)號(hào)場(chǎng)為條件用分布函數(shù)來(lái)描述觀測(cè)數(shù)據(jù)的分布。然后根據(jù)貝葉斯定理[8～12]實(shí)現(xiàn)圖像的分割，但是若依據(jù)傳統(tǒng)公式進(jìn)行計(jì)算，計(jì)算量很大，很難得到理想效果，所以我們假設(shè)觀測(cè)場(chǎng)Y=y每一個(gè)成分是相對(duì)獨(dú)立于給定標(biāo)記場(chǎng)X=x的，并且標(biāo)記場(chǎng)中共有C個(gè)標(biāo)記，這樣的話，凝膠圖像的分割問(wèn)題可表示為：

（1）

其中P（X=x|Y=y）是標(biāo)記場(chǎng)X=x在已知觀測(cè)數(shù)據(jù)Y=y條件下的后驗(yàn)概率，P（yc|X=x）指的是在X=x的條件下觀測(cè)場(chǎng)成分yc的概率分布，P（Y=y）則是已知觀測(cè)數(shù)據(jù)Y=y的概率，在這里由于觀測(cè)場(chǎng)Y=y指的是觀測(cè)圖像，是作為一個(gè)已知量存在的，因此在分析圖像分割的最大后驗(yàn)概率過(guò)程可以忽略，P（Y=y）的影響。P（X=x）是X=x的先驗(yàn)概率，本文通過(guò)基于MRF分割模型MLL來(lái)獲取標(biāo)記場(chǎng)的先驗(yàn)概率，這里僅考慮二階MLL模型，定義其能量函數(shù)如下：

（2）

其中V（xi，xj）是勢(shì)函數(shù)，當(dāng)xi=xj，V（xi，xj）=-1；當(dāng)xi≠xj，V（xi，xj）=1；Ni是第i個(gè)位置的領(lǐng)域；？茁指的是空間位置函數(shù)[8]，它的取值與隸屬度函數(shù)有關(guān)，定義為？茁=1-0.5？滋i（xi），所以當(dāng)隸屬度？滋i（xi）值越小，說(shuō)明其屬于類c的程度越小，反之越大，也就是說(shuō)在這里該點(diǎn)的分類很大程度上取決于它的領(lǐng)域狀況。這樣就得到了其先驗(yàn)概率：

（3）

通過(guò)前面的推導(dǎo)，僅有P（yc|X=x）作為未知量存在，由于一個(gè)類的特征數(shù)據(jù)通常被認(rèn)為是正態(tài)分布的，在此我們將做進(jìn)一步假設(shè)，令圖像中屬于同一類的像素服從高斯分布，那么在計(jì)算高斯分布的時(shí)候引入一個(gè)灰度密度的加權(quán)系數(shù)來(lái)實(shí)現(xiàn)高斯分布灰度均值與方差的獲取與更新，其公式如下：

（4）

Dij指的是像素之間的歐式距離，兩像素值越接近？籽i值將會(huì)越大。之后再進(jìn)行一個(gè)歸一化處理得到其加權(quán)系數(shù)，用于表示像素 對(duì)于分類的影響程度：

（5）

將加權(quán)系數(shù)帶入廣義模糊集的質(zhì)心公式，用質(zhì)心來(lái)替代類的均值，即：

（6）

？滋ij是初始分割模糊聚類算法中的隸屬度。

（7）

由公式（6）和（7）便得到了均值與方差，那么其高斯分布函數(shù)就表示了像素的觀測(cè)場(chǎng)同一類型區(qū)域的分布，即：

（8）

則凝膠圖像分割問(wèn)題可以表示為：

（9）

這里，P（Xm=c|X■）表示標(biāo)記場(chǎng)的局部概率，Nm是位置m的鄰域位置集合。

2分割算法實(shí)現(xiàn)步驟與實(shí)驗(yàn)仿真分析

2.1 分割算法實(shí)現(xiàn)步驟

（1）凝膠圖像的預(yù)分割：利用非局部均值算法（NL-means）對(duì)圖像進(jìn)行濾波并且結(jié)合模糊C均值算法實(shí)現(xiàn)圖像的初始分割，將圖像分為3類；

（2）利用公式（3）獲取圖像的先驗(yàn)概率；

（3）依據(jù)已獲得的先驗(yàn)知識(shí)的先驗(yàn)概率利用公式（9）得到圖像的后驗(yàn)概率；

（4）利用權(quán)重值來(lái)更新類的均值與方差；

（5）判斷收斂與否，若不收斂，重復(fù)步驟（2）～（5），反之算法結(jié)束。

2.2實(shí)驗(yàn)仿真分析

2.2.1 實(shí)驗(yàn)仿真

本文采用凝膠模擬圖像和凝膠真實(shí)圖像進(jìn)行實(shí)驗(yàn)測(cè)試，且選取的圖像大小均為256 256。模擬圖像的分割結(jié)果如圖1所示：圖（a）為凝膠模擬圖像，圖（b）為加入均值為零，方差為0.02的高斯噪聲模擬圖像，圖（c）為原始方法對(duì)于模擬圖像的分割結(jié)果，圖（d）為改進(jìn)算法的分割結(jié)果。真實(shí)圖像的分割結(jié)果如圖2所示：圖（a）為真實(shí)凝膠圖像的原始圖像，圖（b）為原始圖像預(yù)處理結(jié)果，圖（c）為原始算法對(duì)于真實(shí)圖像的分割結(jié)果，圖（d）為改進(jìn)算法對(duì)真實(shí)圖像的分割結(jié)果。

2.2.2 算法的性能描述

首先從人的主觀觀察上可以發(fā)現(xiàn)，原算法與改進(jìn)算法無(wú)論對(duì)模擬圖像還是真實(shí)圖像均能在一定程度上實(shí)現(xiàn)有效分割，但是本文所提出的改進(jìn)算法相對(duì)原始算法來(lái)說(shuō)，明顯提高了對(duì)噪聲的魯棒性，并且對(duì)微弱蛋白點(diǎn)的分割有一定的保障作用，增強(qiáng)了分割效果。

表1是對(duì)原始算法和改進(jìn)算法分割結(jié)果的客觀分析，這里通過(guò)對(duì)總體分割精度（錯(cuò)分率）、kappa系數(shù)[11]以及算法效率進(jìn)行對(duì)比研究。從表1中我們可以看到，本文采用的改進(jìn)算法不僅提高了凝膠圖像的分割精度，而且還減少了計(jì)算量，提高了計(jì)算效率，除此之外，kappa系數(shù)的值越大，說(shuō)明其分割效果越好。

3 結(jié)束語(yǔ)

本文在綜合分析已有文獻(xiàn)算法的基礎(chǔ)上，結(jié)合凝膠圖像的特點(diǎn)，提出了一種改進(jìn)的MRF凝膠圖像分割方法，首先用非局部均值對(duì)圖像進(jìn)行濾波，為凝膠圖像的有效分割提供了一個(gè)良好的基礎(chǔ)，而后通過(guò)模糊c均值算法實(shí)現(xiàn)了圖像的初始聚類分割，將圖像分為3類，再通過(guò)MLL模型模糊聚類作為先驗(yàn)知識(shí)獲取其先驗(yàn)概率并因此得到后驗(yàn)概率，在此過(guò)程中引入了灰度點(diǎn)密度權(quán)值來(lái)實(shí)現(xiàn)聚類均值方差的更新，最終實(shí)現(xiàn)凝膠圖像的最優(yōu)化分割。仿真結(jié)果表明，本文提出的算法具有較高的分割精度和可行性，提高了對(duì)噪聲的消除能力和算法計(jì)算效率。

參考文獻(xiàn)

[1]宋革，姜勇.二維凝膠電泳的新技術(shù)及其應(yīng)用[J].中國(guó)微循環(huán)，2005，9（1）：62-65.

[2]Marc R.Wilkins， Christian Pasquali， et al. From proteins to proteomes large scale protein identification by two-dimensional electrophoresis and arnino acid analysis[J]. Bio/Technology， 1996， 14（1）： 61-65.

[3] LEIT， SEWCHAND W. Statistical approach to X-ray CT imaging and its applications in image analysis-Part 2： A new stochastic model-based image segmentation technique for X-ray CT images [J]. IEEE Transactions on Medical Image， 1992， 11（1）： 62-69.

[4] Charles. A Bouman， Michael Shapiro. A multiscal random field model for Bayesian image segmentation[J]. IEEE Transactions on image processing， 1994， 3（2）： 162-177.

[5] Yanqiu Feng， Wufan Chen. Brain MR image segmentation using fuzzy clustering with spatial constraints based on markov random field theory[J]. Lecture Notes in Computer Science， 2004： 188-195.

[6] Sotirios P Chatzis， Theodora A Varigou. A fuzzy clustering approach toward hidden markov random field models for enhanced spatially constrained image segmentation[J]. IEEE Transactions on Fuzzy Systems， 2008， 16（5）： 1351-1361.

[7] YANG Qinghai， LU Bo， YAN Ziye， et al. Touched string segmentation algorithm based on Markov random field[J]. Computer Engineering， 2013， 39（4）： 258-262.

[8] Huawu Deng？， David A. Clausi. Unsupervised image segmentation using a simple MRF model with a new implementation scheme[J]. Pattern Recognition， 2004， 2（8）： 691 - 694.

[9] Yan Gang， Chen Wufan， Feng Yanqiu. Generalized fuzzy gibbs random filed and research on algorithm for MR image segmentation[J].Journal of Image and Graphics， 2005， 10（9）： 1082-1088.

[10]李旭超. 小波域馬爾可夫隨機(jī)場(chǎng)在圖像處理中的應(yīng)用[M]. 北京：電子工業(yè)出版社.

[11]C. S. Won， H. Derin. Unsupervised segmentation of noisy and textured images using Markov random ？elds[J]， CVGIP： Graphical Models Image Process， 1992， 54 （4） ： 308-328.

[12]劉國(guó)英. 基于Markov隨機(jī)場(chǎng)的小波域圖像建模與分割：Matlab環(huán)境[M]. 北京：北京出版社， 2008.

作者簡(jiǎn)介：張彥清（1988）女，山東泰安人，碩士研究生。主要研究方向?yàn)橹悄苄盘?hào)與信息處理。

侯偉，女，山東濟(jì)南人，副研究員。主要研究方向?yàn)樾盘?hào)與信息處理。

李玲，女，山東濟(jì)南人，副教授。主要研究方向?yàn)樾盘?hào)與信息處理。

辛化梅，女，山東青島人，博士，副教授。主要研究方向?yàn)樾盘?hào)與信息處理。

2分割算法實(shí)現(xiàn)步驟與實(shí)驗(yàn)仿真分析

2.1 分割算法實(shí)現(xiàn)步驟

（1）凝膠圖像的預(yù)分割：利用非局部均值算法（NL-means）對(duì)圖像進(jìn)行濾波并且結(jié)合模糊C均值算法實(shí)現(xiàn)圖像的初始分割，將圖像分為3類；

（2）利用公式（3）獲取圖像的先驗(yàn)概率；

（3）依據(jù)已獲得的先驗(yàn)知識(shí)的先驗(yàn)概率利用公式（9）得到圖像的后驗(yàn)概率；

（4）利用權(quán)重值來(lái)更新類的均值與方差；

（5）判斷收斂與否，若不收斂，重復(fù)步驟（2）～（5），反之算法結(jié)束。

2.2實(shí)驗(yàn)仿真分析

2.2.1 實(shí)驗(yàn)仿真

本文采用凝膠模擬圖像和凝膠真實(shí)圖像進(jìn)行實(shí)驗(yàn)測(cè)試，且選取的圖像大小均為256 256。模擬圖像的分割結(jié)果如圖1所示：圖（a）為凝膠模擬圖像，圖（b）為加入均值為零，方差為0.02的高斯噪聲模擬圖像，圖（c）為原始方法對(duì)于模擬圖像的分割結(jié)果，圖（d）為改進(jìn)算法的分割結(jié)果。真實(shí)圖像的分割結(jié)果如圖2所示：圖（a）為真實(shí)凝膠圖像的原始圖像，圖（b）為原始圖像預(yù)處理結(jié)果，圖（c）為原始算法對(duì)于真實(shí)圖像的分割結(jié)果，圖（d）為改進(jìn)算法對(duì)真實(shí)圖像的分割結(jié)果。

2.2.2 算法的性能描述

首先從人的主觀觀察上可以發(fā)現(xiàn)，原算法與改進(jìn)算法無(wú)論對(duì)模擬圖像還是真實(shí)圖像均能在一定程度上實(shí)現(xiàn)有效分割，但是本文所提出的改進(jìn)算法相對(duì)原始算法來(lái)說(shuō)，明顯提高了對(duì)噪聲的魯棒性，并且對(duì)微弱蛋白點(diǎn)的分割有一定的保障作用，增強(qiáng)了分割效果。

表1是對(duì)原始算法和改進(jìn)算法分割結(jié)果的客觀分析，這里通過(guò)對(duì)總體分割精度（錯(cuò)分率）、kappa系數(shù)[11]以及算法效率進(jìn)行對(duì)比研究。從表1中我們可以看到，本文采用的改進(jìn)算法不僅提高了凝膠圖像的分割精度，而且還減少了計(jì)算量，提高了計(jì)算效率，除此之外，kappa系數(shù)的值越大，說(shuō)明其分割效果越好。

3 結(jié)束語(yǔ)

本文在綜合分析已有文獻(xiàn)算法的基礎(chǔ)上，結(jié)合凝膠圖像的特點(diǎn)，提出了一種改進(jìn)的MRF凝膠圖像分割方法，首先用非局部均值對(duì)圖像進(jìn)行濾波，為凝膠圖像的有效分割提供了一個(gè)良好的基礎(chǔ)，而后通過(guò)模糊c均值算法實(shí)現(xiàn)了圖像的初始聚類分割，將圖像分為3類，再通過(guò)MLL模型模糊聚類作為先驗(yàn)知識(shí)獲取其先驗(yàn)概率并因此得到后驗(yàn)概率，在此過(guò)程中引入了灰度點(diǎn)密度權(quán)值來(lái)實(shí)現(xiàn)聚類均值方差的更新，最終實(shí)現(xiàn)凝膠圖像的最優(yōu)化分割。仿真結(jié)果表明，本文提出的算法具有較高的分割精度和可行性，提高了對(duì)噪聲的消除能力和算法計(jì)算效率。

參考文獻(xiàn)

[1]宋革，姜勇.二維凝膠電泳的新技術(shù)及其應(yīng)用[J].中國(guó)微循環(huán)，2005，9（1）：62-65.

[2]Marc R.Wilkins， Christian Pasquali， et al. From proteins to proteomes large scale protein identification by two-dimensional electrophoresis and arnino acid analysis[J]. Bio/Technology， 1996， 14（1）： 61-65.

[3] LEIT， SEWCHAND W. Statistical approach to X-ray CT imaging and its applications in image analysis-Part 2： A new stochastic model-based image segmentation technique for X-ray CT images [J]. IEEE Transactions on Medical Image， 1992， 11（1）： 62-69.

[4] Charles. A Bouman， Michael Shapiro. A multiscal random field model for Bayesian image segmentation[J]. IEEE Transactions on image processing， 1994， 3（2）： 162-177.

[5] Yanqiu Feng， Wufan Chen. Brain MR image segmentation using fuzzy clustering with spatial constraints based on markov random field theory[J]. Lecture Notes in Computer Science， 2004： 188-195.

[6] Sotirios P Chatzis， Theodora A Varigou. A fuzzy clustering approach toward hidden markov random field models for enhanced spatially constrained image segmentation[J]. IEEE Transactions on Fuzzy Systems， 2008， 16（5）： 1351-1361.

[7] YANG Qinghai， LU Bo， YAN Ziye， et al. Touched string segmentation algorithm based on Markov random field[J]. Computer Engineering， 2013， 39（4）： 258-262.

[8] Huawu Deng？， David A. Clausi. Unsupervised image segmentation using a simple MRF model with a new implementation scheme[J]. Pattern Recognition， 2004， 2（8）： 691 - 694.

[9] Yan Gang， Chen Wufan， Feng Yanqiu. Generalized fuzzy gibbs random filed and research on algorithm for MR image segmentation[J].Journal of Image and Graphics， 2005， 10（9）： 1082-1088.

[10]李旭超. 小波域馬爾可夫隨機(jī)場(chǎng)在圖像處理中的應(yīng)用[M]. 北京：電子工業(yè)出版社.

[11]C. S. Won， H. Derin. Unsupervised segmentation of noisy and textured images using Markov random ？elds[J]， CVGIP： Graphical Models Image Process， 1992， 54 （4） ： 308-328.

[12]劉國(guó)英. 基于Markov隨機(jī)場(chǎng)的小波域圖像建模與分割：Matlab環(huán)境[M]. 北京：北京出版社， 2008.

作者簡(jiǎn)介：張彥清（1988）女，山東泰安人，碩士研究生。主要研究方向?yàn)橹悄苄盘?hào)與信息處理。

侯偉，女，山東濟(jì)南人，副研究員。主要研究方向?yàn)樾盘?hào)與信息處理。

李玲，女，山東濟(jì)南人，副教授。主要研究方向?yàn)樾盘?hào)與信息處理。

辛化梅，女，山東青島人，博士，副教授。主要研究方向?yàn)樾盘?hào)與信息處理。

2分割算法實(shí)現(xiàn)步驟與實(shí)驗(yàn)仿真分析

2.1 分割算法實(shí)現(xiàn)步驟

（1）凝膠圖像的預(yù)分割：利用非局部均值算法（NL-means）對(duì)圖像進(jìn)行濾波并且結(jié)合模糊C均值算法實(shí)現(xiàn)圖像的初始分割，將圖像分為3類；

（2）利用公式（3）獲取圖像的先驗(yàn)概率；

（3）依據(jù)已獲得的先驗(yàn)知識(shí)的先驗(yàn)概率利用公式（9）得到圖像的后驗(yàn)概率；

（4）利用權(quán)重值來(lái)更新類的均值與方差；

（5）判斷收斂與否，若不收斂，重復(fù)步驟（2）～（5），反之算法結(jié)束。

2.2實(shí)驗(yàn)仿真分析

2.2.1 實(shí)驗(yàn)仿真

本文采用凝膠模擬圖像和凝膠真實(shí)圖像進(jìn)行實(shí)驗(yàn)測(cè)試，且選取的圖像大小均為256 256。模擬圖像的分割結(jié)果如圖1所示：圖（a）為凝膠模擬圖像，圖（b）為加入均值為零，方差為0.02的高斯噪聲模擬圖像，圖（c）為原始方法對(duì)于模擬圖像的分割結(jié)果，圖（d）為改進(jìn)算法的分割結(jié)果。真實(shí)圖像的分割結(jié)果如圖2所示：圖（a）為真實(shí)凝膠圖像的原始圖像，圖（b）為原始圖像預(yù)處理結(jié)果，圖（c）為原始算法對(duì)于真實(shí)圖像的分割結(jié)果，圖（d）為改進(jìn)算法對(duì)真實(shí)圖像的分割結(jié)果。

2.2.2 算法的性能描述

首先從人的主觀觀察上可以發(fā)現(xiàn)，原算法與改進(jìn)算法無(wú)論對(duì)模擬圖像還是真實(shí)圖像均能在一定程度上實(shí)現(xiàn)有效分割，但是本文所提出的改進(jìn)算法相對(duì)原始算法來(lái)說(shuō)，明顯提高了對(duì)噪聲的魯棒性，并且對(duì)微弱蛋白點(diǎn)的分割有一定的保障作用，增強(qiáng)了分割效果。

表1是對(duì)原始算法和改進(jìn)算法分割結(jié)果的客觀分析，這里通過(guò)對(duì)總體分割精度（錯(cuò)分率）、kappa系數(shù)[11]以及算法效率進(jìn)行對(duì)比研究。從表1中我們可以看到，本文采用的改進(jìn)算法不僅提高了凝膠圖像的分割精度，而且還減少了計(jì)算量，提高了計(jì)算效率，除此之外，kappa系數(shù)的值越大，說(shuō)明其分割效果越好。

3 結(jié)束語(yǔ)

本文在綜合分析已有文獻(xiàn)算法的基礎(chǔ)上，結(jié)合凝膠圖像的特點(diǎn)，提出了一種改進(jìn)的MRF凝膠圖像分割方法，首先用非局部均值對(duì)圖像進(jìn)行濾波，為凝膠圖像的有效分割提供了一個(gè)良好的基礎(chǔ)，而后通過(guò)模糊c均值算法實(shí)現(xiàn)了圖像的初始聚類分割，將圖像分為3類，再通過(guò)MLL模型模糊聚類作為先驗(yàn)知識(shí)獲取其先驗(yàn)概率并因此得到后驗(yàn)概率，在此過(guò)程中引入了灰度點(diǎn)密度權(quán)值來(lái)實(shí)現(xiàn)聚類均值方差的更新，最終實(shí)現(xiàn)凝膠圖像的最優(yōu)化分割。仿真結(jié)果表明，本文提出的算法具有較高的分割精度和可行性，提高了對(duì)噪聲的消除能力和算法計(jì)算效率。

參考文獻(xiàn)

[1]宋革，姜勇.二維凝膠電泳的新技術(shù)及其應(yīng)用[J].中國(guó)微循環(huán)，2005，9（1）：62-65.

[2]Marc R.Wilkins， Christian Pasquali， et al. From proteins to proteomes large scale protein identification by two-dimensional electrophoresis and arnino acid analysis[J]. Bio/Technology， 1996， 14（1）： 61-65.

[3] LEIT， SEWCHAND W. Statistical approach to X-ray CT imaging and its applications in image analysis-Part 2： A new stochastic model-based image segmentation technique for X-ray CT images [J]. IEEE Transactions on Medical Image， 1992， 11（1）： 62-69.

[4] Charles. A Bouman， Michael Shapiro. A multiscal random field model for Bayesian image segmentation[J]. IEEE Transactions on image processing， 1994， 3（2）： 162-177.

[5] Yanqiu Feng， Wufan Chen. Brain MR image segmentation using fuzzy clustering with spatial constraints based on markov random field theory[J]. Lecture Notes in Computer Science， 2004： 188-195.

[6] Sotirios P Chatzis， Theodora A Varigou. A fuzzy clustering approach toward hidden markov random field models for enhanced spatially constrained image segmentation[J]. IEEE Transactions on Fuzzy Systems， 2008， 16（5）： 1351-1361.

[7] YANG Qinghai， LU Bo， YAN Ziye， et al. Touched string segmentation algorithm based on Markov random field[J]. Computer Engineering， 2013， 39（4）： 258-262.

[8] Huawu Deng？， David A. Clausi. Unsupervised image segmentation using a simple MRF model with a new implementation scheme[J]. Pattern Recognition， 2004， 2（8）： 691 - 694.

[9] Yan Gang， Chen Wufan， Feng Yanqiu. Generalized fuzzy gibbs random filed and research on algorithm for MR image segmentation[J].Journal of Image and Graphics， 2005， 10（9）： 1082-1088.

[10]李旭超. 小波域馬爾可夫隨機(jī)場(chǎng)在圖像處理中的應(yīng)用[M]. 北京：電子工業(yè)出版社.

[11]C. S. Won， H. Derin. Unsupervised segmentation of noisy and textured images using Markov random ？elds[J]， CVGIP： Graphical Models Image Process， 1992， 54 （4） ： 308-328.

[12]劉國(guó)英. 基于Markov隨機(jī)場(chǎng)的小波域圖像建模與分割：Matlab環(huán)境[M]. 北京：北京出版社， 2008.

作者簡(jiǎn)介：張彥清（1988）女，山東泰安人，碩士研究生。主要研究方向?yàn)橹悄苄盘?hào)與信息處理。

侯偉，女，山東濟(jì)南人，副研究員。主要研究方向?yàn)樾盘?hào)與信息處理。

李玲，女，山東濟(jì)南人，副教授。主要研究方向?yàn)樾盘?hào)與信息處理。

辛化梅，女，山東青島人，博士，副教授。主要研究方向?yàn)樾盘?hào)與信息處理。