張書真

(吉首大學信息科學與工程學院，湖南 吉首 41 6000)

基于三維直方圖修正和灰度熵分解的圖像分割

張書真

(吉首大學信息科學與工程學院，湖南 吉首 41 6000)

圖像噪聲容易引起圖像誤分割，而常用閾值選取方法僅依賴于圖像直方圖的概率信息，未直接考慮圖像中類內(nèi)灰度分布的均勻性。為此，提出一種修正三維直方圖和分解處理灰度熵的圖像分割算法。分析圖像的噪聲對其鄰域灰度造成的影響，通過修正三維直方圖來減弱噪聲干擾，給出三維灰度熵閾值的選取公式，并將三維灰度熵分解至一維進行處理，使計算復雜度由O(L3)降為O(L)。實驗結(jié)果表明，與二維最大熵斜分法、二維交叉熵遞推法、降維三維Otsu法相比，該算法抗噪性能更強、分割效果更好，同時能使運算時間縮短10%以上。

圖像分割；閾值選?。粓D像噪聲；三維直方圖；分解處理；灰度熵

1 概述

圖像分割技術(shù)在計算機視覺、模式識別、醫(yī)學圖像處理等領域應用十分廣泛，閾值法則因其簡單且性能穩(wěn)定而成為圖像分割的基本技術(shù)之一，其中基于熵概念的閾值法頗受關(guān)注，常見的有Shannon熵[1-2]、交叉熵[3-4]、Renyi熵[5-6]等算法。但是依據(jù)一維直方圖選擇閾值，由于沒有利用圖像的局部空間信息，當存在噪聲干擾或灰度分布交叉時分割效果不佳。為此，人們推導出相應的二維算法[7-9]和三維算法[10-11]，從而獲得了更好的分割效果。在這些方法中，閾值的選擇僅依賴于二維或三維直方圖的概率信息，而沒有直接考慮圖像中目標和背景類內(nèi)灰度分布的均勻性，因此，對有些圖像分割的效果不理想。此外，算法將搜索空間擴大到二維或三維后，計算復雜度大幅增加。目前，基于三維直方圖的快速算法[12-13]也相應提出，可是算法因存在區(qū)域誤分問題，抗噪性能較弱。文獻[14]提出利用三維直方圖重建來解決區(qū)域誤分問題，但在直方圖重建過程中計算量較大。

針對上述問題，本文提出一種修正三維直方圖和分解灰度熵的圖像分割算法。該算法分析噪聲灰度值對像素鄰域造成的影響，通過對比像素灰度值、鄰域均值、鄰域中值3個參數(shù)，通過簡單計算克服噪聲點干擾。此外，算法中采用的灰度熵不僅利用了直方圖中的概率信息，而且能直接反映目標和背景類內(nèi)灰度值的差異，據(jù)此來選取閾值可改善分割效果。為解決運算速度問題，還針對三維灰度熵提出分解處理的方法，以降低算法的計算復雜度。

2 三維直方圖的修正

無噪Lena圖像及其三維直方圖如圖1所示，圖中3個坐標軸分別表示圖像像素的灰度值x、鄰域均值y和鄰域中值z，若(s, t, q)被選取為閾值點，則三維直方圖被分成目標類Co和背景類Cb2個區(qū)域。通常，圖像目標和背景內(nèi)像素的灰度值、鄰域均值和鄰域中值十分接近，因此，目標和背景內(nèi)的像素點主要分布在直方圖的對角線附近，而邊緣點和噪聲點將分布直方圖中遠離對角線的區(qū)域。由于噪聲點和邊界點與整幅圖像的像素點相比數(shù)量很少，因此可假設這部分像素點的頻數(shù)和為0。

圖1 Len a原圖及其三維直方圖

(a)Lena原圖 (b)三維直方圖

包含混合噪聲的Lena圖像及其三維直方圖如圖2(a)和圖2(b)所示。

圖2 含噪Lena圖及其三維直方圖

實際上，三維直方圖閾值法雖然考慮了像素空間鄰域均值和鄰域中值信息，但并沒有從根本上達到去噪的目的，在三維直方圖上遠離對角線的分布點并不能被忽略掉。文獻[13]提出了一種三維直方圖重建去噪的方法，該方法可將遠離對角線的分布點重新賦值以聚集到對角線附近。但由于該方法對于直方圖每個分布點都要計算灰度值、均值和中值的兩兩絕對差，再通過比較這些絕對差重新賦值，因此計算量較大。為此，本文提出一種快速修正三維直方圖的方法?？紤]到在閾值分割中，相對于邊緣點噪聲點對分割效果影響要大很多，它會導致二值化后的圖像出現(xiàn)大量的“白點”或“黑點”現(xiàn)象，因此，在三維直方圖修正中重點針對噪聲分布點進行修正。在含噪圖像中，無論是目標區(qū)域還是背景區(qū)域，噪聲點總體表現(xiàn)為暗區(qū)域出現(xiàn)亮噪聲或亮區(qū)域出現(xiàn)暗噪聲。而在三維直方圖中，噪聲點分布于遠離對角線的區(qū)域，其中亮噪聲滿足x＞y＞z，暗噪聲滿足x＜y＜z。若要修正噪聲點，只需將噪聲點的灰度值矯正到與其鄰域內(nèi)像素的灰度值接近即可，由于鄰域均值受到噪聲灰度值的影響會比較大，因此修正中將鄰域中值z賦給x，即令x=z。修正后的直方圖如圖3所示，可見通過本文算法，大量噪聲點的灰度值被有效修正，像素點匯聚到三維直方圖的對角線附近。

圖3 含噪Lena圖修正后的三維直方圖

3 三維灰度熵分割算法

3.1 一維灰度熵分割

現(xiàn)有常用的閾值分割方法，僅用到直方圖的概率分布信息，而沒有直接反映圖像中目標和背景內(nèi)類灰度值的差異，因此，本文考慮采用灰度熵進行分割。設一幅圖像總的灰度級數(shù)目為L，圖像中灰度為k的像素數(shù)目為g(k)，k=0,1,…,L-1，用閾值t將圖像劃分為目標類Co和背景類Cb，為討論方便，設低灰度區(qū)為目標類，高灰度區(qū)為背景類，則可令：

圖像灰度熵越大，表明類內(nèi)的像素灰度差異越小，當總灰度熵達到最大時，表明目標類和背景類的灰度趨于均勻，此時對應的t就是最佳閾值*t，即：

3.2 三維灰度熵分割

設一幅圖像總的灰度級數(shù)目為L，其三維直方圖的3個坐標軸分別表示圖像像素的灰度值x、均值y和中值z。若(s, t, q)被選取為閾值點，忽略遠離直方圖對角線上的分布點，則可得到目標類灰度熵向量為：

3.3 三維灰度熵的分解處理

三維閾值法在改善分割效果的同時，計算復雜度也以指數(shù)級方式大幅攀升。文獻[12]通過對三維Otsu進行分解，使計算復雜度降低到O(L)?；陬愃扑枷肟赏茖С鋈S灰度熵的分解公式。

根據(jù)三維直方圖的三元組(x, y, z)，令像素灰度值、鄰域均值、鄰域中值的邊緣頻數(shù)為：

則Fx，Vy，Rz分別表示對應原始圖像、鄰域均值圖像、鄰域中值圖像的灰度直方圖。

設遠離對角線分量近似為0，則有：

由式(19)和式(20)可知，通過分解處理可使最佳三維閾值的計算復雜度降至O(L)。

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

為驗證本文算法的有效性和可行性，實驗分別采用二維最大熵斜分法[8]、二維交叉熵遞推法[9]、重建三維直方圖與降維Otsu方法[14]以及本文算法進行圖像分割，并對分割結(jié)果進行比較。實驗在AMD Athlon Ⅱ X4 640、3.01 GHz CPU和內(nèi)存3.25 G B的微處理器上進行，編程環(huán)境為Matlab7.1。選取其中3幅圖像加以說明，圖像分別是含有高斯噪聲與椒鹽噪聲的Lena、Flower和Barbara圖像。4種算法的分割結(jié)果分別如圖4～圖6所示，閾值與時間比較如表1所示。

從圖4可以看出，相對于其他3種算法，本文算法分割出來的Lena的面部和肩部都比較干凈，視覺效果最好，說明算法中修正直方圖的措施相當有效，在去除噪聲干擾方面能力突出。

從圖5和圖6可以看出，由于其他算法中閾值的選取僅依賴于圖像直方圖，而沒有直接考慮圖像中目標和背景內(nèi)灰度的均勻性，因此對有些圖像分割的效果不理想。例如在Flower圖像中，花托、花柄和葉片部分，其他算法皆不能有效提取為目標區(qū)域；在Barbara圖像中，Barbara的眼部細節(jié)和臉部輪廓的切分，本文算法和文獻[14]算法的效果較好，但本文算法又具有明顯的計算優(yōu)勢。

從表1可以看出，除本文算法，其他3種算法中所需時間最少的是二維最大熵斜分法，而與該算法相比，本文算法運算時間也減少了10%以上。

圖4 含噪Lena圖的分割效果比較

圖5 含噪Flower圖的分割效果比較

表1 不同算法獲得的閾值以及運行時間比較

5 結(jié)束語

本文詳細分析了噪聲對圖像分割造成的不良影響，提出一種快速修正直方圖的方法，通過該方法能有效降低噪聲干擾，提高分割準確度。另外，與現(xiàn)有的閾值分割算法不同，算法中采用的灰度熵不僅考慮了直方圖中的概率信息，而且還能直接反映出圖像中目標和背景的灰度值差異，當總灰度熵達到最大時，各類內(nèi)灰度趨于均勻。而相對于三維灰度熵而言，經(jīng)過分解處理，其計算復雜度大大降低。通過對多幅圖像進行閾值分割，驗證了本文算法的可行性。但本文算法中三維直方圖的構(gòu)建僅采用了簡單的均值和中值濾波，對于被噪聲干擾非常嚴重的圖片，其分割效果是不理想的。因此，下一步將研究更有效的直方圖構(gòu)建方法。

[1] Kapur J N, Sah oo P K, Wang A K C. A New Method for Grey-level Pictur e T hresholding Using the Entropy of the Histogram[J]. Computer Vision, Graphics and Image Processing, 1985, 29(3): 273-285.

[2] 吳成茂. 基于加權(quán)香農(nóng)熵的圖像閾值法[J]. 計算機工程與應用, 2008, 44(18): 177-180.

[3] Li Ch unhung, Le e C K. Minimum Cross E ntropy Thresholding[J]. Pattern Recognition, 1993, 26(4): 617-625.

[4] 曹建農(nóng). 基于直方圖重構(gòu)的極大交叉熵圖像分割方法[J].計算機應用, 2011, 31(12): 3373-3377.

[5] Sahoo P K, Wilkins C, Y enger J. Threshold Selection Using Renyis’ Entropy[J]. Pattern Recognition, 1997, 30(1): 71-84.

[6] 雷 博, 范九倫. 一維Renyi熵閾值法中參數(shù)的自適應選取[J]. 光子學報, 2009, 38(9): 2439-2442.

[7] 吳一全, 潘 喆, 吳文怡. 二維直方圖區(qū)域斜分的最大熵閾值分割算法[J]. 模式識別與人工智能, 2009, 22(1): 162-168.

[8] 雷 博, 范九倫. 灰度圖像的二維交叉熵閾值分割法[J].光子學報, 2009, 38(6): 1572-1576.

[9] 卞紅雨, 劉 翠. 基于修正二維熵的水聲圖像分割[J].計算機工程, 2010, 36(14): 193-195.

[10] 孫文邦, 唐海燕. 基于三維直方圖聚類特性的非監(jiān)督圖像變化檢測方法[J]. 上海交通大學學報, 2012, 46(2): 245-249.

[11] 龍建武, 申鉉京, 魏 巍, 等. 一種結(jié)合紋理信息的三維Renyi熵閾值分割算法[J]. 小型微型計算機系統(tǒng), 201 1, 32(5): 947-952.

[12] 范九倫, 趙 鳳, 張雪峰. 三維Otsu閾值分割方法的遞推算法[J]. 電子學報, 2007, 35(7): 1398-1402.

[13] 龔 劬, 倪 麟, 唐萍峰, 等. 基于分解的三維Otsu圖像分割快速算法[J]. 計算機應用, 2012, 32(6): 1526-1528.

[14] 申鉉京, 龍建武, 陳海鵬, 等. 三維直方圖重建和降維的Otsu閾值分割算法[J]. 電子學報, 2011, 39(5): 1108-1114.

編輯 金胡考

Image Segmentation Based on Three-dimensional Histogram Correction and Gray Entropy Decomposition

ZHANG Shu-zhen

(School of Information Science and Engineering, Jishou University, Jishou 416000, China)

Aiming at the problem of inaccurating image segmentation caused by image noise and the common threshold selection methods which only rely on the probabilistic information from th e image histogram is without directly thinking of the uniformity of th e image inter-class gray distribution, a threshold selection algorithm based on a three-dimensional histogram correction and gray entropy decomposition is proposed. It analyzes the influence of image noise to the gray of pi xel’s neighborhood region, and reduces the noise interference by modifying the three-dimensional h istogram. A formula of threshold selection based on three-dimensional gray entropy is presented, and the dimension of gray entropy is decomposed to one dimension, which makes the computation complexity reduced fro m O(L3) to O(L). Experimental results show that, co mpared with two-dimensional maximum entropy a lgorithm based o n oblique segmentation, two-dimensional cross entropy algorithm based on recursi on and three-dimensional Otsu algorithm based on dim ension reduction, the presented algorithm has better anti-noise performance, visual quality and the operation time is reduced by about 10% at least.

image segmentation; threshold selection; image noise; three-dimensional histogram; decomposition processing; gray entropy

10.3969/j.issn.1000-3428.2014.05.048

國家自然科學基金資助項目(61262032)；湖南省教育廳科學研究基金資助項目(12C0314)。

張書真(1977－)，女，講師、碩士，主研方向：圖像處理，模式識別。

2013-03-26

2013-04-22E-mail：sunny_zsz@126.com

1000-3428(2014)05-0234-04

A

TP391