白雪冰，陳 凱，郭景秋，祝 賀，張庭亮

(東北林業(yè)大學(xué) 機(jī)電工程學(xué)院，哈爾濱 150040)

自20世紀(jì)開始，計(jì)算機(jī)技術(shù)得到了長足的進(jìn)步，數(shù)字圖像處理也隨之快速的發(fā)展起來。數(shù)字圖像處理就是利用計(jì)算機(jī)或其他數(shù)字硬件，對(duì)從圖像信息轉(zhuǎn)換而得的電信號(hào)進(jìn)行某些數(shù)學(xué)運(yùn)算，以提高圖像的實(shí)用性[1-3]。目前數(shù)字圖像處理所處理的圖像信息大多數(shù)是二維信息，所以信息量很大，加大了圖像處理的難度，因此也對(duì)計(jì)算機(jī)會(huì)有較高的要求，如存儲(chǔ)容量和計(jì)算速度等。由于數(shù)字圖像中各個(gè)像素點(diǎn)之間具有較大的相關(guān)性與不獨(dú)立性，所以，圖像處理中能夠具有較快的處理速度和精度與圖像信息的壓縮還有很大的發(fā)展空間[4-6]。

本文結(jié)合目前的彩色圖像處理算法提出了基于k-mean聚類與灰度-梯度最大熵的樹木圖像分割算法，該算法為了更好地區(qū)分目標(biāo)圖像與背景，先采用k-mean聚類初分割，再利用灰度-梯度空間能夠明確的描述圖像中各像素點(diǎn)灰度與梯度的分布規(guī)律及圖像中的所需目標(biāo)圖像與背景之間的邊緣情況[7]，基于灰度-梯度最大熵進(jìn)行精分割，因此，圖像分割所得到的結(jié)果將更為理想，精度得到一定提高，處理速度也有所提高。

1 k-mean聚類初分割

k-means算法是一種經(jīng)典的基于距離的算法，用距離的遠(yuǎn)近來評(píng)價(jià)其相似性，即認(rèn)為兩個(gè)對(duì)象距離越遠(yuǎn)，相似性就越小，而距離越近，其相似度就越大，因此該算法最終得到緊湊且獨(dú)立的簇。

算法過程如下：

(1)從整體N個(gè)類別中任意選出K個(gè)類別作為聚類中心。

(2)對(duì)其他的每個(gè)類別計(jì)算其到每個(gè)聚類中心的距離，并把該類別歸為最近的聚類中心的類。

(3)再重新計(jì)算已經(jīng)分出的每一個(gè)類的聚類中心。

(4)迭代(2)～(3)步直至新的中心與原中心相等或小于指定閾值，算法結(jié)束。

K-mean聚類算法將對(duì)象分為K個(gè)類別C={ck，i=1，2，…，K}每類含ck有一個(gè)聚類中心μi，應(yīng)用歐式距離公式計(jì)算類內(nèi)各點(diǎn)到聚類中心μi的距離平方和：

。

(1)

圖1 K-mean聚類初分割

從圖中看出，圖像被分別聚成了兩類，三類和四類。聚成三類或者四類的時(shí)候，出現(xiàn)了目標(biāo)和背景區(qū)域混淆的情況，效果較差。圖1中原圖是一副背景簡(jiǎn)單的樹木圖像，聚成兩類比較合適。

2 灰度-梯度最大熵算法

2.1 灰度-梯度共生矩陣

描述圖像的過程中，將灰度分布與梯度信息融合得灰度梯度共生矩陣[2](GGCM)，設(shè)一幅M×N灰度圖像f(m,n);m=0,1…M-1,n=0,1,…N-1，則灰度歸一化圖像和梯度歸一化圖像分別為：

(2)

(3)

式中：Lf為最大灰度級(jí)；Lg為最大梯度級(jí)；fmax為灰度矩陣最大值；gmax為梯度矩陣最大值[4]。

GGCM定義為H=(hxy),x=0,…,Lf-1,y=0,…,Lg-1，其中元素H(x,y)定為F(m,n)與G(m,n)中具有灰度級(jí)x和梯度級(jí)y的像點(diǎn)數(shù)，歸一化圖像為：

(4)

2.2 灰度-梯度最大熵

共生矩陣是Lf×Lg維的矩陣的圖像，設(shè)(Lf-1)×(Lg-1)大小的區(qū)域?yàn)閳D像的二維直方圖，如圖2所示，共生矩陣坐標(biāo)系的原點(diǎn)在左上角，其中橫坐標(biāo)為圖像的梯度值且橫向遞增，縱坐標(biāo)為其灰度值且縱向遞增。定義閾值在(s,t)，由于閥值處灰度值與梯度值區(qū)分明顯，目標(biāo)的灰度值較低而背景的灰度值較高，所以共生矩陣就被分為A、B、C、D四個(gè)象限。本算法是以灰度-梯度共生矩陣和最大熵基礎(chǔ)，充分結(jié)合了圖像的灰度和梯度信息。

圖2 共生矩陣四個(gè)區(qū)域

物體和背景的灰度值變化不大，所以灰度值比較?。篈(0≤i≤s，0≤j≤t)表示物體，D(s+1≤i≤Lf-1,0≤j≤t)表示背景。隨著t值不斷變大，其對(duì)應(yīng)像素點(diǎn)是目標(biāo)與背景邊界的概率增加，因此共生矩陣在B(0≤i≤s，t+1≤j≤Lg-1)中的元素kij為灰度i屬于物體和梯度j屬于邊緣轉(zhuǎn)移數(shù)目，在C(s+1≤i≤Lf-1,t+1≤i≤Lg-1)中元素的kij為灰度i屬于背景和梯度j屬于邊緣數(shù)目，所以4個(gè)象限像素和為：

(5)

將pij正規(guī)化為：

(6)

由此可知其條件熵為：

(7)

(8)

所以圖像的條件熵為：

(9)

依據(jù)最大熵理論，H獲得最大值的(s*，t*)即是最佳閾值[6]：

(10)

最后對(duì)灰度-梯度最大熵的分割結(jié)果進(jìn)行形態(tài)學(xué)后處理，通過標(biāo)注連通對(duì)象，計(jì)算標(biāo)注的面積，去除面積比較小的部分，可以實(shí)現(xiàn)對(duì)背景非目標(biāo)區(qū)域信息比較理想的清除，再通過膨脹操作將圖像邊緣中的薄弱斷開處進(jìn)行連接，通過腐蝕操作將粘連的目標(biāo)部分?jǐn)嚅_。從分割結(jié)果看出，獲得了完整的閉合曲線輪廓。

3 實(shí)驗(yàn)結(jié)果

基于以上理論基礎(chǔ)，應(yīng)用計(jì)算機(jī)，Matla7.0集成環(huán)境[11-13]對(duì)樹木圖像進(jìn)行分割，分割過程及其結(jié)果分析如圖3所示。

圖3 實(shí)驗(yàn)結(jié)果

對(duì)經(jīng)過灰度-梯度最大熵方法的分割后結(jié)果進(jìn)行形態(tài)學(xué)后處理操作，從分割結(jié)果看出，獲得了完整的閉合曲線輪廓，獲得更多的樹木邊緣信息，細(xì)節(jié)處分割效果較好，對(duì)于樹木中的明顯空洞也分割的較為清晰[6]，分割時(shí)間為t=102 407.000us。

應(yīng)用灰度-梯度最大熵對(duì)不同圖像分割結(jié)果如圖4所示。

圖4 不同背景下樹木圖像分割

從分割結(jié)果看出，獲得了完整的閉合曲線輪廓，細(xì)節(jié)處分割效果較好，對(duì)于樹木中的明顯空洞也分割的較為清晰，整個(gè)分割過程運(yùn)行時(shí)間可以接受，但是在聚類過程可能會(huì)嘗試幾次確定合適的類數(shù)，導(dǎo)致分割具有一定程度的復(fù)雜性[5]。

二維最大熵與灰度-梯度最大熵分割結(jié)果對(duì)比如圖5～圖7所示。

圖5 樣本一

圖6 樣本二

圖7 樣本三

兩種分割方法運(yùn)行時(shí)間對(duì)比見表1。

表1 兩種分割方法運(yùn)行時(shí)間對(duì)比

將之前的分析結(jié)果結(jié)合表格1可以直觀的看出，對(duì)于樹木圖像分割，灰度-梯度最大熵分割結(jié)果的分割效果要好于二維最大熵分割結(jié)果[7]的分割效果，且運(yùn)行時(shí)間明顯小于其運(yùn)行時(shí)間。結(jié)合上述各組實(shí)驗(yàn)的結(jié)果圖像的分析可知，灰度-梯度最大熵分割結(jié)果可以分割出較為完整的閉合曲線，因?yàn)榻Y(jié)合了圖像的梯度信息，所以可以獲得較好的邊緣輪廓線。

4 結(jié)束語

該算法對(duì)彩色樹木圖像的分割具有較好的分割效果。對(duì)于背景比較復(fù)雜的樹木圖像，運(yùn)用基于k-mean聚類與灰度-梯度最大熵算法的分割算法對(duì)圖像進(jìn)行分割，可以獲得較完整的閉合曲線，效果較好，再通過形態(tài)學(xué)后處理可以消除背景非目標(biāo)區(qū)域的干擾，且程序運(yùn)行速度較快，精度較高。當(dāng)然這種方法也存在缺點(diǎn)，對(duì)于樹木中的空洞繁多，基于k-mean聚類與灰度梯度最大熵算法分割獲得的圖像邊緣和原圖像邊緣會(huì)有出入，這樣就可能丟失了部分的樹木信息。

【參 考 文 獻(xiàn)】

[1]王克奇，馬曉明，白雪冰.基于分形理論和數(shù)學(xué)形態(tài)學(xué)的木材表面缺陷識(shí)別的圖像處理[J].森林工程，2013，29(3)：48-50+53.

[2]梁 建.基于灰度-梯度共生矩陣的植被提取方法[J].海洋測(cè)繪，2013，33(1)：29-31.

[3]朱 峰，王海豐，任洪娥.基于Gabo變換的紋理圖像分割算法及應(yīng)用[J].森林工程，2013，29(5)：60-63.

[4]郭 烈，王榮本，金立生，等.基于二維最大熵閾值分割的坑識(shí)別方法[J].計(jì)算機(jī)工程與應(yīng)用，2006(21)：226-228.

[5]王曉松.復(fù)雜背景下樹木圖像提取研究[D].北京：北京林業(yè)大學(xué)，2010.

[6]蔡世捷.基于Matlab的樹木圖像分割方法研究[D].南京：南京林業(yè)大學(xué)，2005.

[7]洪繼光.灰度-梯度共生矩陣紋理分析方法[J].自動(dòng)化學(xué)報(bào)，1984，10(1)：22-25.

[8]周德龍，申石磊.基于灰度-梯度共生矩陣模型的最大熵閾值處理算法[J].小型微型計(jì)算機(jī)系統(tǒng)，2002，23(2)：136-138.

[9]陳 果，左洪福.圖像分割的二維最大熵遺傳算法明[J].計(jì)算機(jī)輔助與圖形算法，2006，14(6)：530-533.

[10]王再尚.Markov隨機(jī)場(chǎng)的樹木表面缺陷模式識(shí)別方法的研究[D].哈爾濱：東北林業(yè)大學(xué)，2012.

[11]蔣慶豐.K-Means聚類算法研究及圖形演示的實(shí)現(xiàn)[J].信息技術(shù)，2010，10(3)：22-25.

[12]闞江明，李文彬.基于數(shù)學(xué)形態(tài)學(xué)的樹木圖像分割方法[J].北京林業(yè)大學(xué)學(xué)報(bào)，2006，40(S2)：12-30.

[13]季 虎，孫即祥，邵曉芳，等.圖像邊緣提取方法及展[J].計(jì)算機(jī)工程與應(yīng)用，2004，40(14)：70-73.