基于贝叶斯算法的二值化算法

基于贝叶斯算法的二值化算法

白　洁1,杨耀权1,陈余梅2

(1.华北电力大学控制科学与工程学院,河北保定071003;

2.广东省湛江电力有限公司,广东湛江524000)

摘要:针对在图像二值化过程中动态选取阈值难的问题,在分析了全局阈值法和局部阈值法各自优缺点的基础上,提出了一种基于贝叶斯算法的全局阈值法和局部阈值法相结合的二值化方法。经实验证明,该方法既能够有效地消除光照不均匀对图像的影响,较好地保留目标图像的细节,又能够有效地消除伪影,提高处理速度。

关键词:二值化;全局阈值法;局部阈值法;贝叶斯算法

中图分类号:T P391.41 文献标识码:A 文章编号:1007-2691(2007)03-0065-03 Application of binarization based on Bayes algorithm

BAI Jie1,YANG Yao-quan1,CHEN Yu-mei2

(1.School of Co ntrol Science and Eng ineering,Nor th China Electric Pow er U niversity,Bao ding071003,China;

2.Zhanjiang Electric Power Co.Ltd.,Zhanjiang524000,China)

A bstract:Aiming at the difficulties of selecting threshold value dynamically in the process of imag e binarization,a bi-narization method based on Bay es algo rithm is derived,this method is a combination of g lobal threshold method and lo-cal threshold method.Ex periments prove that this method is able to eliminate the effects of unequal illumination and ar tifacts,retain target image's details and raise the processing speed.

Key words:binarization;global threshold;local threshold;Bayes algo rithm

0　引　言

图像的二值化在计算机图像处理技术中广泛应用于目标识别、字符识别、牌照识别等领域。而图像阈值自动选取方法的研究长期以来吸引着众多学者,寻找简单实用、自适应强的阈值自动选取方法是这些研究者们的共同目标。

常用的二值化方法有各自的优缺点,如全局阈值算法简单,对于目标和背景明显分离、直方图分布呈双峰的图像效果良好,但对噪声较大和阴影明显的图像二值化效果不佳。局部阈值法可以克服上述不足,但要比较图像中所有点,所以速度慢。本文在分析了全局阈值法和局部阈值法各自优缺点的基础上,提出了一种全局阈值法和局部阈值法相结合的二值化方法。经实验证明,该方法抗噪能力强、并能较好的保留图像细节。1　基于灰度图的二值化

假设用f(x,y)(0≤x

f(x,y)=

1　f(x,y)≥T

0　otherw ise

式中:T为阈值。经过二值化处理后,图像中的目标和背景就由黑白两种颜色分开,选择不同的阈值会得到不同的划分结果。常用的二值化方法主要可以分为全局法和局部法两种。

1.1　全局动态二值化和局部自适应二值化

全局阈值法根据图像的直方图或灰度空间分布确定一个阈值,以此实现灰度图像到二值图像的转化。典型的全局阈值方法包括Ostu方法[1]、

第34卷第3期2007年5月 华北电力大学学报

Journal of No rth China Electric P ower U niversity

Vol.34,No.3

M ay,2007

收稿日期:2006-11-24.

最大熵方法等。全局阈值算法简单,对于目标和背景明显分离、直方图分布呈双峰的图像效果良好,但其对于由于光照不均匀、噪声干扰较大等原因使直方图分布不呈双峰的图像,二值化效果明显变差。局部阈值算法通过定义考察点的邻域,并由邻域计算模板,实现考察点灰度与邻域点的比较。非均匀光照条件等情况虽然影响整体图像的灰度分布却不影响局部的图像性质,从而使局部阈值算法比全局阈值算法有更广泛的应用,Bemsen[2]算法是一种典型的局部阈值算法。但局部阈值算法也存在以下问题和缺点:如实现速度比全局阈值算法慢;不能保证字符笔划的连通性;容易出现伪影(g host)现象(背景区域受到噪声干扰出现笔划)等。动态阈值法是一种自适应得二值化算法,它利用了像素自身及其邻域灰度变化特征,充分考虑了每个像素邻域的特征,能够更好的突出背景和目标的边界,使得相距很近的两条线不会产生粘连现象,效果较好。

1.2　常用方法简介

(1)双峰法

首先,根据输入的数字图像按某一灰度级内的像素频数做出该图像的直方图。大量统计表明,对于某一类图像(如目标和背景有较强的对比)直方图中出现两个峰值,其中一个处于背景灰度区中,另一个处于目标灰度区。在这样的前提下,在两个峰值的中间波谷处取阈值。

(2)微分直方图法

此算法利用灰度变化率即微分直方图来确定阈值。如果图像中的目标和背景的边界处于灰度值急剧变化的部分,不直接利用其灰度值,而是利用其微分值作为阈值。微分值是灰度的变化率,可以有多种定义,比如某一像素和其周围邻域各像素的灰度差的最大值,或是各个差的绝对平均值,由此求出微分直方图。微分直方图中的峰值就是所求的阈值。此方法适用于图像的目标与背景的边界位于灰度值急剧发生变化的情况,它的微分直方图会有一个峰值。但对于边界附近灰度变化复杂的图像,此方法是不太奏效的。

(3)最大方差法

在图像的灰度直方图中,把图像的灰度值集合用阈值分为两组,通过基于两组的最佳分离来确定阈值,即根据两组的组间方差和组内方差之比为最大来确定阈值。使组间方差与组内方差之比达到最大的即为阈值[3]。此方法在直方图不存在峰值时也可以使用,但是不能反映图像的几何结构,有时判断标准与人的视觉不一致。

2　基于贝叶斯公式的全局法和局部法相结合的二值化算法

在分析了各种全局阈值法和局部阈值法各自优缺点的基础上,提出一种全局阈值法和局部阈值法相结合的二值化方法,该方法抗噪能力强、并能较好的保留图像细节。

2.1　全局阈值选取

这里采用一种基于贝叶斯公式和最大熵法则求全局最佳阈值的方法,该方法能较好的保持视觉效果,适当保留细节[4]。

设G为输入图像的灰度图,则其有L层灰度的直方图H={h0,h1,……,h L-1},灰度概率密度为P(i)=h(i)/N,i=0,1…,L-1,N为G像素总数。G被最优阈值T分为G b(背景)和G f (前景)两部分,G b=∪

L-1

i=0

G bi,G f=∪

L-1

i=0

G bi G fi(2)所以由全概率公式得:

P(G b)=∑

L-1

i=0

p i p b│i(1)

P(G f)=∑

L-1

i=0

p i p f│i(2)式中:

p b│i=

P(G bi)

P(G i)

,……,p f│i=

P(G f i)

P(G i)

(3)

因为:p b│i=1-p f│i,则

P(G b)=1-P(G f)(4)采用香农最大熵原理作为判别准则,则熵函数E为

E=-P(G b)lg P(G b)-P(G f)lg P(G f)(5)把式(4)代入,则式(5)可改写为

E=-P(G f)lg P(G f)-

(1-P(G f))lg(1-P(G f))(6)当P(G f)=0.5时熵函数E取最大。然而在实际应用中,由于直方图离散,P(G f)很少能完全等于0.5,可以转而求最小误差:

E min:=P(G f)-0.5)

=΢

L-1

i=0

P i P f│i-0.5(7)条件概率函数P f i应满足式(7)的性质,选用非线性凸抛物线形式。

66华北电力大学学报 2007年