一種改進的2維Otsu紅外圖像分割法研究

張建波，楊恢先，周彤彤，譚正華，李 淼

（1.湘潭大學信息工程學院，湘潭411105；2.湘潭大學材料與光電物理學院，湘潭411105）

一種改進的2維Otsu紅外圖像分割法研究

張建波1，楊恢先2*，周彤彤2，譚正華1，李 淼2

（1.湘潭大學信息工程學院，湘潭411105；2.湘潭大學材料與光電物理學院，湘潭411105）

為了在2維直方圖上用Otsu方法更好地分割紅外圖像、提高抗噪能力，提出了一種改進的方法。首先分析在2維灰度-鄰域均值直方圖上的分割存在不準確性，采用2維灰度-梯度直方圖，且改進對鄰域均值的求取算法；然后對Otsu法的閾值函數進行研究，引入類內的分離信息改進閾值函數，并簡化該閾值函數以降低運算復雜度，通過實驗給出了相應的實驗對比。結果表明，改進的方法能更好地分割目標，運行時間較少、抗噪性更強。

圖像處理；Otsu法；鄰域均值；閾值函數

引 言

圖像分割是把圖像分成互不重疊的不同區域的過程，同一區域內具有特性相似性，不同區域間具有特性相異性。圖像分割是圖像分析和理解的基礎，閾值分割因簡單有效，在圖像識別與計算機視覺等領域得到廣泛的應用。常用的閾值法有最大類間方差法（Otsu法）、最大熵法和最小交叉熵法等［1］。最大類間方差法是一種以目標和背景之間的類間方差作為分離性準則的閾值分割方法，它依據1維直方圖選取閾值，由于沒有利用圖像的局部空間信息，當存在噪聲干擾或灰度分布交叉時分割效果不佳，因此，LIU等人［2］將其拓展到了2維。2維Otsu法除了考慮像素點的灰度信息外，還考慮了像素點與其鄰域的空間相關信息，在由灰度值和鄰域均值灰度構建的2維直方圖上通過一個閾值點將圖像劃分成4個矩形區域。但是計算閾值時僅僅考慮兩個沿對角線的矩形區域，忽略了位于閾值向量附近且像素灰度級與其鄰域平均灰度級相近的一半區域，并簡單地認為該區域出現的概率為0，因此建立在該2維直方圖上的區域劃分所得到的分割結果并不準確［3］。同時Otsu分割方法只考慮前景類和背景類之間的方差，當圖像在2維直方圖上沒有明顯的雙峰（如只出現一個尖峰）情況時，該方法的分割性能迅速下降，甚至不能正確地分割出目標［4］，因此作者提出一種改進的2維Otsu分割方法，并應用于紅外圖像分割。

1 2維Otsu分割方法

1.1 2維直方圖

圖1a所示為1幅紅外圖像的原始圖像，參考文獻［2］中引入鄰域空間信息將分割直方圖由1維推廣到2維，在由像素灰度值和鄰域均值灰度值構建的2維直方圖上（如圖1b所示），通過一個閾值點（s，t）將直方圖劃分成4個區域，如圖1c所示。假設圖像的亮像素視為背景，暗像素視為目標，則區域0代表目標，區域2代表背景，而區域1和區域3代表邊緣和噪聲。計算閾值時僅考慮區域0和區域2，而假設區域1和區域3概率近似為0，其存在的缺點有：（1）區域0和區域2中可能存在邊緣點或噪聲點被視為目標或背景內點的情況；（2）忽略了位于閾值向量附近且靠近對角線的概率分布，會將一部分目標和背景內點劃分成噪聲和邊緣點。因此，在由灰度值和鄰域均值灰度值構建的2維直方圖上對圖像進行區域劃分會存在分割的不準確性。為此采用基于像素灰度和像素梯度的2維直方圖，該直方圖橫坐標采用像素灰度值，縱坐標采用像素梯度值即灰度值與其鄰域均值灰度值之差的絕對值［3］。因為目標內點和背景內點像素灰度級比較均勻、梯度較小、所占比例大；邊緣點和噪聲點梯度較大，在整幅圖像像素中所占比例相對較小，2維灰度-梯度直方圖上的區域劃分的優點是不僅盡可能地包含了所有目標點和背景點，同時也進一步剔除了2維灰度-鄰域均值直方圖區域劃分方法中目標和背景區域可能存在的噪聲點的影響?；诨叶扰c灰度梯度的2維直方圖如圖1d所示。

設圖中任意一點（x，y）的灰度值為f（x，y），其灰度級為L，則該點的m×n鄰域均值灰度值g（x，y）的求取如下式所示：

式中，m，n一般取奇數，??為取下整數?？紤]m，n取得過大，算法消耗的時間更長，大部分文獻里面采用3×3模板來求取鄰域均值，如參考文獻［2］中將像素點周圍所有的鄰域信息都考慮在內；參考文獻［5］中只考慮水平和垂直方向的鄰域信息；參考文獻［6］中只考慮對角方向上的鄰域信息。灰度共生矩陣中的相關特征統計量［7-8］反映了局部灰度之間相關性的程度，如果像素周圍的鄰域信息和灰度值越接近，則像素點分布越集中，在直方圖上的峰值越明顯。顯然，將像素周圍鄰域信息都考慮在內和只考慮某方向上的鄰域信息算法會削弱分割的效果。因此，需要對鄰域均值灰度值的求取算法進行改進。同樣是在3×3鄰域內將每個鄰域灰度值與灰度值在某一范圍內進行比較，實驗表明，這一范圍取為10時，鄰域信息能盡可能地被包含在內，得到的鄰域均值灰度值跟灰度值越接近，2維直方圖上峰值越明顯。因此鄰域均值灰度值的求取算法為：（1）在3×3鄰域范圍內，比較鄰域灰度值和灰度值差值的絕對值是否在10范圍內，若滿足條件則把該鄰域灰度值計算在內，否則不計算；（2）記錄鄰域范圍內滿足條件的個數；（3）由前兩步得到的所有滿足條件的鄰域灰度值除以滿足的個數，即可得到像素點的鄰域均值灰度值?？紤]到圖像難免會受噪聲的影響，當鄰域信息有受到污染時，像素的鄰域均值灰度信息就會失去其在含噪圖像分割中積極的指導作用，而對于圖像中的每一個像素而言，都可以找到與它具有相似鄰域結構的很多像素，因此，可利用參考文獻［9］中的方法找到具有相似鄰域結構的像素，對其進行加權平均得到該像素的鄰域信息，再通過改進的算法求鄰域均值。

1.2 Otsu方法

在改進的2維灰度-梯度直方圖上通過閾值點（s，t）將直方圖重新劃分為4個區域，如圖1f所示，則區域0代表目標，區域2代表背景，而區域1和區域3代表噪聲和邊緣。因此可以求得目標和背景的概率w0，w1及均值矢量μ0，μ1分別為：

式中，μ0i，μ0j分別表示目標類在灰度值和灰度梯度值兩個方向上的均值矢量。

式中，μ1i，μ1j分別表示背景類在灰度值和灰度梯度值兩個方向上的均值矢量。

總體均值矢量μt為：

基于類間方差的準則［10］采用離散矩陣Sb的跡Sb，tr作為閾值函數，則：

Otsu法是使得離散矩陣的跡Sb，tr最大時取得最佳閾值點（s，t）從而將圖像分割開來。

2 改進的Otsu紅外圖像分割法

圖2是圖1所示紅外圖像［11-12］的1維直方圖，左邊較平坦的代表飛機目標，右邊的尖峰代表天空背景。

由（2）式～（5）式可知：

由于目標所占比例小于背景所占比例，即w0＜w1，又有w0+w1=1，所以w0＜1/2，即1-2w0＞0。另外目標為整幅圖像較暗的部分，可得關系式：

又由（3）式得：μ0i=，代入到（10）式中可得＜μt，i，即＜A+B。所以w0A+w0B-A＞0。

因此可知（8）式μt，i-μ0i-μt，i-μ1i＞0，同理可證明μt，j-μ0j-μt，j-μ1j＞0，說明圖像均值μt，i更接近背景區。而OTSU只考慮了類間的方差，采用各類中心到總體中心距離的加權和，當類間距離最大時就會使分割閾值偏向于峰大的背景一方［13］，從1維直方圖上分析，其分割的閾值就會偏大。因此引入類內信息，類內分離信息反映的是每一類像素中包含的分類信息［14］，當分割出來的兩類內聚性最小時，表明其分割效果較好。因此分別求出兩類的類內方差，并修改平方項以降低運算復雜度，得到目標類絕對差散布矩陣δ0和背景類絕對差散布矩陣δ1，如下所示：

則目標類和背景類的類內散布矩陣Sw為：

由（2）式～（5）式，將（6）式類間方差法的閾值函數化簡成如下形式：

修改閾值函數去掉平方項得到基于絕對差的形式，如下式所示：

式中，μi=w0μ0i，μj=w0μ0j。

新閾值分割函數采用類間和類內的離散矩陣跡的絕對值比值形式，即：

3 實驗結果及分析

分別對3幅人為加入均值為0、方差為0.01的高斯噪聲的圖像進行實驗對比。其中圖3a為普通圖像，其在2維直方圖上具有較好的峰值信息；圖4a和圖5a均為紅外圖像，其在2維直方圖上都表現為一個尖峰；用傳統Otsu方法對圖像進行分割的結果如圖3b、圖4b和圖5b所示；圖3c、圖4c和圖5c為用Otsu方法在灰度-梯度直方圖上的分割結果；圖3d、圖4d和圖5d為提出的改進方法的分割結果。所有實驗均在AMD Athlon（tm）64×2主頻為2.71G的CPU、內存為1GB DDR RAM的機器上用MATLAB進行實驗對比。

cameraman圖像在2維直方圖具有較好的峰值信息，用傳統Otsu方法也可以較好地分割圖中的目標，如人和相機部分，但是在傳統2維直方圖上的區域劃分存在分割的不準確性，因此相機下面部分的三角架不能被較好地分割。而建立在灰度-梯度直方圖的Otsu方法以及提出的改進方法都可以較好的分割三角架，同時提出的改進方法改進了對鄰域均值灰度值的求取算法，其去噪的效果更好。對于圖4a和圖5a兩幅紅外圖像，因為在直方圖上沒有表現為較明顯的峰值信息，而是呈現出一個尖峰的情況，用Otsu方法對圖像進行分割的效果很差，不能較好地分割目標，而通過引入類內的分離信息對閾值函數進行修正可以較好地分割圖像。表1為分割圖像的運行時間對比。

從表1可以看出，修改閾值函數平方項，降低運算的復雜度，可以提高分割的運行速度。

4 結 論

分析了傳統2維直方圖區域劃分分割不準確的情況，并對鄰域均值灰度值的求取算法進行改進，重新構建2維灰度-梯度直方圖；對Otsu紅外圖像分割方法進行研究，引入類內分離信息，提出一種改進的Otsu紅外圖像分割方法。實驗表明，改進的2維Otsu方法可以較好地分割紅外圖像，且抗噪性更強，具有一定實用價值。

［1］ SEZGIN M，SANKUR B.Survey over image thresholding techniques and quantitative performance evaluation［J］.Journal of Electronic Imaging，2004，13（1）：145-165.

［2］ LIU J Zh，LI W Q.The automatic thresholding of gray-level pictures via 2-D Otsu method［J］.Acta Automatica Sinica，1993，19（1）：101-105（in Chinese）.

［3］ WU Y Q，ZHANG J K.Thresholding based on maximum entropic correcation of average gray level-gradient 2-D histogram［J］.Journal of Chinese Computer Systems，2009，30（8）：1675-1679（in Chinese）.

［4］ HOU Z，HU Q，NOWINSKI W L.On minimum variance thresholding［J］.Pattern Recognition Letters，2006，27（14）：1732-1743.

［5］ ZHANG X M，FENG Y Zh，YAN H L，et al.Improved two-dimensional minimum error image thresholding method［J］.Computer Science，2012，39（8）：259-262（in Chinese）.

［6］ ZHANG Y F，ZHANG Y.Automatic threshold of image segmentation using 2-D entropy［J］.Journal of Harbin Engineering University，2006，27（3）：353-356（in Chinese）.

［7］ YUAN L H，FU L，YANG Y，et al.Analysis of texture feature extracted by gray level co-occurrence matrix［J］.Journal of Computer Applications，2009，29（4）：1018-1021（in Chinese）.

［8］ WU G，TANG Zh M，CHENG Y，et al.Object tracking method based on gray level co-occurrence matrix texture characteristic［J］.Journal of Nanjing University of Science and Technology（Natural Science Edition），2010，34（4）：459-463（in Chinese）.

［9］ ZHAO F，FAN J L，PAN X Y，et al.Two-dimensional Otsu's curve thresholding segmentation method based on gray and non local spatial gray feature［J］.Application Research of Computers，2012，29（5）：1987-1989（in Chinese）.

［10］ JIANG Q Y，LI P，SUN L.Application of Otsu method in motion detection system［J］.Journal of Computer Applications，2011，31（1）：260-262（in Chinese）.

［11］ KANG L Zh，CHEN F Sh，WANG D Sh，et al.Detection method for infrared small target based on mathematical morphology［J］.Optoelectronic Engineering，2010，37（11）：26-31（in Chinese）.

［12］ WEI X，MA L H，LI Y X，et al.Infrared image enhancement algorithm based on image segmentation and platform histogram equalization［J］.Infrared Technology，2012，34（5）：271-275（in Chinese）.

［13］ WANG Zh B，GU Y，LI Zh Q.Threshold image segmentation based on maximum scatter difference discriminant criterion［J］.Journal of Applied Optics，2010，31（3）：403-407（in Chinese）.

［14］ ZHANG J K，WU Y Q.Image thresholding based on improved 2-D minimum within-cluster absolute difference method and its fast algorithm［J］.Signal Processing，2010，26（4）：552-557（in Chinese）.

Improved segmentation method of 2-D Otsu infrared image

ZHANG Jianbo1，YANG Huixian2，ZHOU Tongtong2，TAN Zhenghua1，LI Miao2
（1.College of Information Engineering，Xiangtan University，Xiangtan 411105，China；2.Faculty of Material and Photoelectronic Physics，Xiangtan University，Xiangtan 411105，China）

In order to gain better segmentation result of infrared images，and improve the ability to resist noise，an improved 2-D Otsu method was proposed.The inaccurate segmentation in the 2-D gray-neighborhood average histogram was analyzed，and 2-D gray-gradient histogram was adopted.A new algorithm to gain neighborhood average value was put forward.Information of within-cluster was added to amend threshold function，which was further simplified to reduce the calculation complex.Experiments show that the improved method can segment the target better，gain better noise resistance and cost less time.

image processing；Otsu method；neighborhood average value；threshold function

TN911.73

A

10.7510/jgjs.issn.1001-3806.2014.05.029

1001-3806（2014）05-0713-05

湖南省教育廳科研資助項目（10C1263）；湘潭大學科研資助項目（11QDZ11）

張建波（1987-），男，碩士研究生，主要研究方向為圖像處理、人工智能。

*通訊聯系人。E-mail：yanghx87@163.com

2013-09-03；

2013-10-16