(廣東科學(xué)技術(shù)職業(yè)學(xué)院 軟件學(xué)院，廣東 珠海 519090)

0 引言

隨著電子技術(shù)和計算機技術(shù)的飛速發(fā)展，圖像采集技術(shù)得到了極大提高，圖像處理技術(shù)的應(yīng)用逐漸成為圖像分析領(lǐng)域研究的熱點，已經(jīng)滲透至人們?nèi)粘Ｉ畹母鱾€方面，在視頻會議、遠程醫(yī)療、農(nóng)業(yè)生產(chǎn)自動化、以及工業(yè)產(chǎn)品質(zhì)量監(jiān)測方面有著重要應(yīng)用[1-2]。圖像分割作為圖像分析的最關(guān)鍵步驟，無論是圖像處理還是圖像理解與識別，其基本工作均建立在圖像分割的基礎(chǔ)上，圖像分割的目的就是將圖像劃分成各具特征的區(qū)域并提取目標區(qū)域的細節(jié)信息的過程[3-4]。在圖像處理過程中只要涉及到目標檢測，信息提取、圖像診斷時，首先就需要對圖像進行分割，在圖像處理過程中十分關(guān)鍵[5]。圖像分割技術(shù)發(fā)展至今，研究人員提出了很多分割方法，但仍然存在很多不足，例如分割精度不高、分割后的圖像分辨率低、分割計算效率低等眾多問題[6]。在這種情況下，尋求一種可擴展性強、計算效率高、抗噪性強的圖像分割方法具有重要的現(xiàn)實意義[7]。

目前，可擴展性較強的圖像分割方法有：文獻[8]提出一種基于馬爾科夫隨機場的圖像智能分割方法。該方法先采用分數(shù)階微分運算對圖像細節(jié)信息的敏感性進行分析，并將其引入圖像特征提取過程中，以得到全面完整的圖像紋理信息，此外，為智能分割圖像不同的紋理區(qū)域，采用模糊熵準則對圖像分割結(jié)果進行平滑去噪，降低圖像噪聲干擾。該方法可以有效提高圖像智能分割準確率，但分割后的圖像視覺效果較差。文獻[9]提出一種基于蟻群優(yōu)化的圖像智能分割方法。該方法先表征圖像的近似度圖，并將其定義為蟻群的棲息環(huán)境，對其進行歸一化分割指導(dǎo)，依據(jù)螞蟻的覓食行為將近似的頂點匯聚在一起，以圖像多層的形式完成圖像分割。該方法分割的準確率較高，但圖像分割速度較慢。文獻[10]提出一種基于形狀共享初始化的圖像智能分割方法。該方法先提取測試圖像的局部形狀，再找出測試樣本圖像庫中與其局部形狀近似的形狀集，依據(jù)測試圖像與樣本圖像中局部形狀的相對位置以及大小進行全局形狀映射，最后依據(jù)圖像全局形狀的覆蓋率進行分組，組合成一系列初始形狀，將圖像的初始輪廓作為主動輪廓模型的初始迭代函數(shù)，對測試圖像進行區(qū)域分割。該方法收斂速度較快，但分割精度較低。

本文充分利用Contourlet變換在圖像方向信息提取和圖像稀疏性表達方面的優(yōu)勢，通過采用小波變換替換拉普拉斯變換進行圖像多個方向子帶分解，并依據(jù)圖像最大似然函數(shù)估計準則獲得圖像的初始分割，此過程達到了抑制圖像噪聲，保護圖像邊緣輪廓的目的。在此基礎(chǔ)上采用自適應(yīng)上下文結(jié)構(gòu)從圖像粗尺度的分割結(jié)果融合至最細尺度，使得最終獲得的分割結(jié)果清晰度較高。

1 基于人工蜂群的圖像智能分割方法

在進行圖像智能分割過程中，首先對采集的圖像進行灰度形態(tài)學(xué)中的閉操作預(yù)處理，將圖像閾值定義為人工蜂群算法中的蜜蜂，采用二維Otsu法設(shè)計人工蜂群算法的適應(yīng)度值，依據(jù)人工蜂群不同工種蜜蜂的分工協(xié)作和信息共享逐步逼近最優(yōu)分割閾值，具體過程如下所述：

圖像智能分割過程中，假設(shè)f(x,y)代表采集的圖像，b(x,y)代表圖像結(jié)構(gòu)元素，其本身為一幅子圖像，用結(jié)構(gòu)元素b對輸入圖像f進行灰度膨脹：

(f⊕b)(s,t)=

max{f(s-x,t-y)+b(x,y)|(s-x),}

(1)

其中：Df和Db分別代表f和b的定義域，x和y分別代表圖像空間中的兩個點，(s,t)代表圖像最優(yōu)分割閾值。(s-x)和(t-y)必須在f的定義域內(nèi)，x和y需在b的定義域內(nèi)，對圖像進行灰度膨脹的結(jié)果是，比圖像背景亮的區(qū)域得到擴張，比圖像背景暗的區(qū)域受到收縮，選用結(jié)構(gòu)元素b對待分割的圖像進行灰度腐蝕處理，即：

f°b=(fΘb)⊕b

(2)

利用下式表示用結(jié)構(gòu)元素b對待分割圖像f進行開操作運算，可得：

f·b=(fΘb)⊕b

(3)

通過對待分割的圖像進行形態(tài)學(xué)處理，進一步抑制了圖像噪聲，在此基礎(chǔ)上將人工蜂群算法應(yīng)用于圖像智能分割過程中，各個蜜源位置對應(yīng)一個分割閾值，蜂群內(nèi)部蜜源的花蜜量等同于適應(yīng)度函數(shù)值，首先，采蜜蜜蜂依據(jù)記憶在鄰近區(qū)域內(nèi)進行搜索產(chǎn)生一個新的位置，當新的采蜜位置優(yōu)于蜂群搜索記憶中的最優(yōu)位置時，用新位置取代記憶最優(yōu)位置，蜂群內(nèi)部的觀察峰依據(jù)所獲得的信息按照與花蜜量相關(guān)的概率選取一個蜜源位置，并在該蜜源位置鄰近區(qū)域進行搜索，具體步驟如下：

1)蜂群X={xij|i=1,2,…,N}的初始化，N代表人工蜂群的規(guī)模大小，xi={xij|i=1,2}代表第i只蜜蜂所表征的圖像分割閾值。

2)計算蜂群內(nèi)部各個采蜜蜂的適應(yīng)度函數(shù)值，將當前蜂群中的最大花蜜量及其蜜源位置發(fā)布在公告板內(nèi)。

3)種群內(nèi)部采蜜峰依據(jù)記憶的信息在其鄰域內(nèi)進行搜索，產(chǎn)生新的位置vi={vij|j=1,2}：

vij=xi+(xi-xkj)×φ

(4)

式中，xkj代表第k只蜜蜂所表征的圖像分割閾值，φ代表取值區(qū)間[-1,1]內(nèi)的隨機數(shù)，計算采蜜蜂新位置vi的適應(yīng)度函數(shù)值，依據(jù)貪婪準則，假設(shè)vi優(yōu)于xi，則更新fi；依據(jù)xi的適應(yīng)度函數(shù)值fi，計算概率Pi：

(5)

(6)

式中，abs(fi)代表花蜜量最大的位置，n代表蜜蜂總數(shù)。

5)依據(jù)從蜂群那里得到的信息，觀測峰依據(jù)概率Pi選取一個蜜源位置，在該位置鄰近區(qū)域進行搜索，從而產(chǎn)生一個新的位置，并計算該位置下適應(yīng)度函數(shù)值，對比采蜜峰搜索的新的位置和種群內(nèi)部記憶中位置的花蜜量。

6)判定蜂群搜索過程中是否有應(yīng)放棄的位置，假設(shè)有，則用偵查峰隨機選取的一個新位置來取代它：

xij=min+(max-min)×φ

(7)

其中：min和max代表圖像智能分割問題域的上、下界，定義分割閾值的取值區(qū)間在[0,255]之內(nèi)。

7)重復(fù)步驟(3)～(6)，直到達到最大迭代次數(shù)，輸出最優(yōu)蜂群蜜源位置，即圖像最優(yōu)分割閾值(s,t)。

采用人工蜂群算法進行圖像分割過程中，相對于對比度高的圖像能夠起到對噪聲的抑制作用，且分割速度較快，但對于實際對比度低的圖像，采用上述方法難以避免地會帶來噪聲增強過度以及細節(jié)特征增強不足的現(xiàn)象。

2 基于Contourlet變換的圖像智能分割方法

2.1 基于小波變換的圖像增強

在進行圖像智能分割過程中，采用方向濾波器組實現(xiàn)圖像各個方向紋理分離，采用小波變換替換拉普拉斯變換進行圖像多個方向子帶分解，具體過程如下所述：

(8)

式中，0≤k<2l-1，2l-1≤k<2l包含了圖像中接近于垂直方向、水平方向的紋理信息，l代表圖像分解級數(shù)，將圖像頻域分解為2l個子帶。

(9)

設(shè)定，圖像灰色區(qū)域?qū)儆诙S小波空間Wj,HL、Wj,LH、Wj,HH，j代表空間尺度，Ψj,HL(n)、Ψj,LH(n)、Ψj,HH(n)代表其相應(yīng)的基函數(shù)，則可利用下式描述圖像尺度空間Vj與小波空間Wj之間的關(guān)系：

(10)

(11)

(12)

(13)

(14)

(15)

式中，M代表圖像原始數(shù)據(jù)總量，J代表圖像處理過程中小波轉(zhuǎn)換層數(shù)，Ld代表圖像第d層小波子帶進行方向濾波的級數(shù)。

2.2 基于自適應(yīng)上下文尺度間融合的圖像智能分割

在進行圖像智能分割過程中，以3.1節(jié)通過Contourlet變換獲得的圖像各個尺度上的初始分割為基礎(chǔ)，計算圖像多尺度似然函數(shù)，并依據(jù)圖像最大似然函數(shù)估計準則獲得圖像的初始分割，在此基礎(chǔ)上采用自適應(yīng)上下文結(jié)構(gòu)從圖像粗尺度的分割結(jié)果融合至最細尺度，獲得最終的圖像智能分割結(jié)果。具體過程如下所述：

(16)

(17)

(18)

(19)

其中：p(ci|di,vi)代表圖像塊di從屬于類別ci的概率，依據(jù)最大后驗概率融合圖像上下文信息，獲得圖像各個像素點的類別標簽由下式確定：

(20)

從圖像的某一穩(wěn)健尺度進行分析，通過選取圖像數(shù)據(jù)塊大小8×8，逐一向圖像細尺度融合，結(jié)合基于上下文的尺度間融合方法依次獲得圖像分辨率為4×4的數(shù)據(jù)塊、像素級的圖像精細分割結(jié)果。

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

為了驗證本文提出的基于Contourlet變換的圖像智能分割方法的有效性，選取多種類型灰度不均的醫(yī)學(xué)圖像作為實驗樣本，并與文獻[8]方法和文獻[10]方法進行了比較，實驗采用Matlab2016a對本文方法進行了驗證，實驗計算機配置為Windows7操作系統(tǒng)下CPU為Intel 2.2 GHz，內(nèi)存為4 G。

對于以下實驗選取相同實驗參數(shù)，設(shè)定時間步長為Δt=0.2，當實驗測試圖像灰度分布均勻時，控制參數(shù)α′=1和β′=1，當測試圖像灰度不均勻時，參數(shù)設(shè)定為α=0.1和β′=1，由于要分割出目標區(qū)域，長度參數(shù)需要設(shè)定較小的值，將該值設(shè)定為μ=0.01×2552。

3.1 參數(shù)設(shè)定及評價方法

為了定量分析和評價本文方法和文獻[8]方法以及文獻[10]方法的分割質(zhì)量，選用Jaccard相似度系數(shù)來衡量3種不同分割方法的精度。Jaccard相似度系數(shù)定義如下：

JS(A,B)=|A∩B|/|A∪B|

(21)

式中，A代表標準的圖像分割結(jié)果，B代表待評價的圖像分割結(jié)果，Jaccard相似度系數(shù)越高說明圖像分割方法的性能越好。

清晰度指標可以準確描述圖像分割后的細節(jié)反差和紋理變化特征，利用下式進行定義：

(22)

式中，ΔIx和ΔIy分別x和y方向上的差分，圖像的平均梯度越大，說明分割后的圖像清晰度越好。

3.2 圖像分割結(jié)果

實驗選取兩組灰度不均的圖像為例驗證本文提出的基于Contourlet的圖像智能分割方法的有效性。選取一幅醫(yī)學(xué)圖像進行分割，本文方法分割灰度不均勻的醫(yī)學(xué)圖像的效果分割結(jié)果用圖1進行表示，圖2為Montreal神經(jīng)研究所提供的醫(yī)學(xué)腦部圖像。

圖1 本文方法灰度不均醫(yī)學(xué)圖像分割結(jié)果

圖2 腦圖像中心提供的白質(zhì)、灰質(zhì)、腦脊液標準分割

圖1中，(a)代表待加權(quán)灰度醫(yī)學(xué)圖像，(b)、(c)、(d)分別代表本文方法分割的腦部圖像白質(zhì)區(qū)域、灰質(zhì)區(qū)域、腦脊液區(qū)域。將本文方法進行圖像智能分割的結(jié)果與標準分割結(jié)果進行對比可知，本文方法能夠?qū)︶t(yī)學(xué)腦部圖像進行準確分割，與標準分割結(jié)果輪廓劃分基本一致，并且從圖中可以看出本文方法分割后的圖像梯度較大，分割后的圖像清晰度較好，證明了本文方法在腦部圖像整體分割上的應(yīng)用是有效的。

表1給出了本文方法和文獻[8]方法以及文獻[10]方法進行圖像智能分割的分割精度(%)、迭代次數(shù)(次)以及運行時間(s)。

表1 不同方法圖像分割精度、迭代次數(shù)以及運行時間

分析表1可知，文獻[8]方法、文獻[9]方法和文獻[10]方法在每次迭代過程中進行了卷積運算，因此運行時間更長，而本文方法僅需在迭代之前進行卷積運算，相比文獻提及方法具有更高的分割計算效率，從圖像分割結(jié)果來看，本文方法相比文獻提及方法對灰度不均的醫(yī)學(xué)圖像具有更好的分割結(jié)果和更快的計算效率。

采用相同的收斂條件，對腦部圖像加入椒鹽噪聲測試，對比4種不同方法分割帶噪聲醫(yī)學(xué)圖像的能力，如表2所示。

分析表2可知，采用相同的收斂條件，在添加相同的椒鹽噪聲強度0.05 dB下，文獻[8]方法的醫(yī)學(xué)腦部圖像分割相對錯誤率為6.58%，文獻[9]方法的醫(yī)學(xué)腦部圖像分割相對錯誤率為6.91%，文獻[10]方法的醫(yī)學(xué)腦部圖像分割相對錯誤率為7.52%，而本文方法分割相對錯誤率遠低于文獻[8]方法、文獻[9]方法和文獻[10]方法，說明本文所提方法對于醫(yī)學(xué)腦部圖像中普遍存在的噪聲具有較好的抗噪性能。

表2 不同方法對帶噪聲的醫(yī)學(xué)圖像分割結(jié)果的影響

4 結(jié)論

近年來，圖像分割方法已經(jīng)成為圖像處理研究領(lǐng)域的的一個研究重點。由于Contourlet變換方法的優(yōu)點是對圖像高頻細節(jié)可進行任意方向的分解，提出一種基于Contourlet變換的圖像智能分割方法。實驗結(jié)果表明，所提方法能夠有效區(qū)分醫(yī)學(xué)腦部圖像不同的紋理區(qū)域，在腦部圖像邊緣細節(jié)以及方向信息的保持上相比其他分割方法具有優(yōu)勢，充分證明了所提方法的有效性。

[1]王 艷.圖像分割的自適應(yīng)快速水平集演化模型[J].計算機仿真,2016,33(2):269-272.

[2]張新明,涂 強,尹欣欣.混合遷移的高效BBO算法及其在圖像分割中的應(yīng)用[J].計算機科學(xué)與探索,2016,10(10):1459-1468.

[3]王靈矯,鐘益群,郭 華,等.基于高階條件隨機場模型的改進型圖像分割算法[J].計算機工程,2016,42(6):241-246.

[4]吳 實,熊淑華,李 杰,等.基于簡單線性迭代聚類算法和動態(tài)區(qū)域合并的巖屑圖像分割[J].科學(xué)技術(shù)與工程,2016,16(35):238-243.

[5]李竹林,王 靜.一種基于濾波的分段點乘圖像分割算法[J].電子設(shè)計工程,2016,24(23):1-3.

[6]龍群芳,張 東.一種基于KNN紋理分類的超聲圖像自動分割方法[J].科技通報,2016,32(7):138-142.

[7]張永梅,巴德凱,邢 闊.基于模糊閾值的自適應(yīng)圖像分割方法[J].計算機測量與控制,2016,24(4):126-128.

[8]曹家梓,宋愛國.基于馬爾科夫隨機場的紋理圖像分割方法研究[J].儀器儀表學(xué)報,2015,36(4):776-786.

[9]葛 亮,楊竣鐸.基于蟻群優(yōu)化多層圖劃分的彩色圖像分割方法[J].計算機應(yīng)用研究,2015,32(4):1265-1268.

[10]梅夢杰,徐 軍.基于形狀共享初始化主動輪廓模型的圖像分割方法[J].計算機應(yīng)用研究,2015,32(6):1902-1905.