可以证明,用r的累积分布函数作为变换函数可产生一幅灰度级分布具有均匀概率密度的图像。
其结果扩展了像素取值的动态范围。
通常把为得到均匀直方图的图像增强技术叫做直方图均衡化处理或直方图线性化处理。
用离散形式表示累积分布函数为:
(2-11)
kns
sk=T(/-J=X—R=0,12…厶1
戶on
2.3.2直方图均衡化步骤
(1)列岀原始图像灰度级fj,j=01・・・k…丄
⑵统讣各灰度级的象素数目,nj,
⑶计算原始图像直方图Pf(fj)=nj/n,n为原始附像总的象素数目;
⑷计算累积分布函数c(f);
(5)应用转移函数,计算映射后的灰度级,gi=INT[(gmax-gmin)c(f)+gmin+0.5]
⑹统让映射后各灰度级的象素数目ni,匸0,1p-1;
(7)计算输出图像直方图Pg(gi)=nj/n,i=0,l,-,p-l;
⑻用fj和gi的映射关系,修改原始图像灰度级,获得直方图近似均匀分布的输出图像。
2.4直方图规定化
(2-13)
(2-14)
(2-15)
在某些情况下,人们并不一定需要具有均匀直方图的图像,有时需要具有特左的直方图的图像,以便能够增强图像中某些灰度级。
直方图规泄化方法就是针对上述思想提出来的。
直方图规定化是使原图像灰度直方图变成规定形状的直方图而对图像作修正的增强方法,英基本思想。
直方图规定化是在运用均衡化原理的基础上,通过建立原始图像和期望图像之间的关系,选择地控制直方图,使原始图像的直方图变成规定的形状,从而弥补了直方图均衡化不具备交互作用的特性人们希望增强后的图像,英灰度级的分布不是均匀的,而是具有规赵形状的直方图,这样可突出感兴趣的灰度范圉。
令Pr(r)和Pz(z)分别为原始图像和期望图像的灰度概率函数,对两者均作直方图均衡化处理,应有:
N=G~x(V)
V=G(Z)=f^(Z>/Z
式子表明可以由均衡化后的灰度变疑v获得期望图像的灰度变量z,这就总味着可以由原始图像均衡化后的图像灰度值来il•算期望图像的灰度值。
因为对原始图像和期望图像都进行了均衡化处理,所以Ps(s)和Pv(v)具有相同的槪率密度,直方图规定化处理后的新图像将具有事先规左的概率密度Pz(z),从而达到预期处理效果
由此可见,直方图规范化就是把图像直方图均衡化结果映射到期望的理想直方图上,使图像按照人们的意愿转换。
直方图规左化使图像增强的实质是••
(1)增大两个占有较多像素灰度之间的差距,一般来讲,背景和目标占有较多的象素,这样实际上加大了背景和目标的对比度,增大了反差。
(2)归并占有较少的像素,通常,目标和背景的过渡处的像素较少,由于归并,其或者变为背景点或者变为目标点,从而使边界变得陡峭,使图像细右淸晰,达到图像增强的目的.
(3)
图(C)、(c)是将图像(A)按图(b)的直方图进行规定化得到的结果及其直方图。
通过对比可以看出图(C)的对比度同图(B)接近一致,对应的直方图形状差异也不大。
这样有利于影像融合处理,保证融合影像光谱特性变化小。
2.5图像平滑
一幅原始图像在获取和传愉过程中会受到各种噪声的干扰,使图像质量下降,时分析图像不利。
这些噪声干扰使图像退化,质量下降。
表现为图像模糊,特征淹没,对图像分析不利为了抑制噪声、改善图像质童,要时图像进行平滑处理。
图像平滑处理的方法多种多样,有邻域平均法、掩膜平滑法,空间低通滤波、噪声门限法、中值滤波法,多幅图像平均法等。
2.5.1图像平滑方法原理
在空间域平滑滤波有很多种算法,苴中最常见的有线性平滑、非线性平滑、自适应平滑。
(1)线性平滑就是对每一个像素的灰度值用它的邻域值来代替,其邻域的大小为NxN,N—般取奇数。
经过线性平滑滤波,相当于图像经过了一个二维的低通滤波器,可是虽然是降低了噪声,但同时也模糊了图像边缘和细卩,这是这类滤波器存在的通病。
(2)非线性平滑是对线性平滑的一种改进,即不对所有像素都用它的邻域平均值来代替,而是取一个闭值,当像素灰度值与其邻域平均值之间的差值大于已知值时才以均值代替当像素灰度值与其邻域平均值之间的差值不大于闭值时取其本身的灰度值。
非线性平滑可消除一些孤立的噪声点,对图像的细节影响不大,但对物体的边缘会带来一立的失真。
(3)自适应平滑是一种根据当时、当地情况来尽量不模糊边缘轮廓为目标进行控制的方法,所以这种算法要有一个适应的目标。
根据目的的不同,可以有各种各样的自适应图像处理方法。
下而再分别简单介绍几种线性平滑、非线性平滑算法和其它一些算法。
2.5.2邻域平均法
邻域平均法是一种局部空间域处理的算法。
设一幅图像f(x,y)为NxN的阵列,平滑后的图像为f(x,y),它的每个像素的灰度级由包含在(x,y)的预怎邻域的几个像素的灰度级的平均值所决定,即用下
式得到平滑的图像。
g(忑y)=—■工W・j)(2-17)
Mijes
式中的x/y=l,l,2N-1,S是(x,y)点邻域中心点的坐标的集合[不包括点(x,y)],M是
S内坐标点的总数。
以上方法简单,计算速度快,但它的主要缺点是在降低噪声的同时使图像产生模糊,特别在边沿和细肖处,邻域越大,模糊越厉害。
为了减少这种效应,可以采用阈值法。
这样平滑后的图像会比邻域平均法模糊度减少。
当某些点的灰度值与并邻点灰度的均值差别较大时,它必然是噪声,则取英邻域平均值作为该点的灰度值,它的平滑效果仍然是很好的。
为了克服简单局部平均的弊病,目前己提出许多保边沿保细巧的局部平滑算法,它们讨论的课题都在如何选择邻域的大小、形状和方向,如何选择参加平均的点数以及邻域各点的权重系数等,它们有:
灰度最相近的及个邻点平均法,梯度倒数加校平滑,最大均匀性平滑,小斜而模型平滑等等.如果将受噪声干扰的图像看成是一个二维随机场,则可以运用统汁理论来分析受噪声干扰的图像平滑后的信噪比问题,一般的噪声属于加性噪声,在独立和分布的髙斯噪声的情况下,我们定义信噪比为含噪图像的均值与噪声方之比,则含噪图像经邻域平均法平滑之后,其信噪比将提高M1/2倍(M为邻域中包含的像素数目),可见邻域取得愈大,像点愈多,则信噪比提高愈大,平滑效果好.
2.5.3空间域低通滤波法
从信号的角度看,信号缓慢变化主要分布在频率域的低频部分,而信号迅速变化的部分主要集中在髙频部分。
对图像来说,它的边缘以及噪声干扰的频率分量都处于频率较高的部分,因此可以用低通滤波方法去除噪声。
而频率域滤波可以用空间域的卷积来实现,为此只要恰当地设计空间域系统冲激响应矩阵就可以达到滤波的效果。
设f(x,y)为带有噪声的原始图像(大小NxN),g(x,y)为经滤波后的输岀图像(大小MxM),h(x,y)为滤波系统的脉冲响应函数(大小LxL),则存在
g(x,y)=f(x,y)*h(x,y)(2-18)
苴中f(x,y)是含有噪声图像的博里叶变换,g(x,y)是平滑处理后的图像之傅里叶变换,g(x,y)是传递函数。
选择传递函数h(x,y),利用h(x,y)使f(x,y)的髙频分量得到衰减,得到g(x,y)后再经反傅里叶变换就可以得到所希望的平滑图像g(x,y)。
根据前而的分析,显然h(x,y)应该具有低通滤波特性,所以这种方法被称之为低通滤波法平滑化处理。
2.5.4多图像平均法
多幅图像平均法是利用对同一景物的多幅图像相加取平均来消除噪声产生的高频成分。
多幅图像取平均处理常用于摄像机的视频图像中,以减少电视摄像机光电摄像管或CCD器件所引起的噪声。
这时对同一景物连续摄取多幅图像并将其数字化,再对多幅图像求平均
一般选用幅图像取平均,这种方法在实际应用中的难点在于如何把多幅图像配准,以便使相应的像素能正确地对应排列。
设g(x,y)为有噪声图像,n(x,y)为噪声,f(x,y)为原始图像,可用下式表示:
g(x,y)=n(x,y)+f(x,y)(2-19)
多图像平均法是把一系列有噪声的图像{g(x,y)}迭加起来,然后再取平均值以达到平滑的目的.当作平均处理的噪声图像数目增加时,其统讣平均值就越接近原始无噪声图像。
这种方法在实际应用中的最大因难在于把多幅图像配准,以便使相应的像素能正确地对应排列。
2.5.5中值滤波法
中值滤波也是一种典型的空间域低通滤波器,它的目的是保护图像边缘的同时去除噪声。
所谓中值滤波,就是指把以某点(XxY)了为中心的小窗口内的所有像素的灰度按从大到小的顺序排列,将中间值作为少处的灰度值若窗口中有偶数个像素,则取两个中间值的平均。
例:
采用43窗口进行中值滤波
原图像为:
22621244424
处理后为:
22222244444
它对脉冲丁扰及椒盐噪声的抑制效果好,在抑制随机噪声的同时能有效保护边缘少受模糊。
但它对点、线等细节较多的图像却不太合适。
对中值滤波法来说,正确选择窗口尺寸的大小是很重要的环节。
一般很难事先确立最佳的窗口尺寸,需通过从小窗口到大窗口的中值滤波试验,再从中选取最佳的。
中值滤波容易去除孤立点、线的噪声,同时保持图像的边缘,它能很好地去除二值噪声,但对髙斯噪声无能为力。
要注意的是,当窗口内噪声点的个数大于窗口宽度一半时,中值滤波的效果不是太好。
2.5.6噪声门限法
噪声门限法是一种简单易行的消除噪声的方法,它对于因噪声传感器或者信道引起的呈现离散分布的单点噪声具有较好的效果,运用噪声门限法进行图像平滑时,首先设左门限值,然后顺序检测图像中的每一个像素,将该像素与其他像素进行比较判断,以确定是否为噪声点:
若为噪声点,则以英邻域内所有像素灰度平均值代替,否则,以原灰度值输出。
假设像素(i,j)岀的灰度为f(i,j),以给该像素为中心取一个NxN的窗口(N=3,5,7...),该窗口内的和汁数器像素点组成集合A,定义灰度差值门限T,误差计算器Cnt和计数器门限值Y。
对每个窗口,Cnt的初始值都是0。
对集合A中的每一个像素点(i/)的灰度f(i;T),若满足
e
(2-21)
(2-22)
则误差计算器加1,位置(i,j)的输出为
g(bj)=f(i,j)(Cntg(i,j)pf(i:
j')/(N2-l)(Cnt>Y)
当窗口顺序移过整幅图像,即可完成噪声平滑。
需要注意的是,该方法中门限值T的选择至关重要,T太大,则噪声平滑不够,T太小,平滑图像就会变得模糊,汁数器门限值的选择一般在窗口内像素的一半附近。
2.5.7掩膜平滑法
图像中存在这样一个基本事实:
同一区域内部的像素之间灰度变化平缓,起伏较小,统计方差小:
在区域边缘,像素之间灰度值得起伏变化大,统讣方差大。
掩膜平滑法的目的在于进行滤波操作的同时,尽可能不破坏区域边缘的细节。
掩膜平滑以一个5x5的窗口为基准,中心位宜为(j,k),在这个窗口中确左9种不同的掩膜模版。
在平滑时,首先计算各模版的均值和方差。
(2-23)
(2-24)
Ai=[》f(j+rrbk+n)]/Q
Bi=》{[f(j+m,k+n)]2-Ai2}
式中,i表示掩膜板编号,QM应掩膜模版中包含像素的个数,(m,n)为掩膜模版中像素相对于中心像素(j,k)的位移量。
也就是说,掩膜平滑的输出为具有最小方差的模版所对应的灰度均值。
当同样的方法作用于图中的每一个像素后,即可得到平滑的图像,平滑图像中相对很好的保留了图像区域边缘的细肖。
2.5.8图像平滑处理方法比较
图像平滑处理的方法多种多样,每种方法在不同的方而各有优点和缺点。
中值滤波比较容易去除椒盐噪声,同时能较好的保持图像的边缘,它还能很好地去除二值噪声,但对髙斯噪声却无能为力。
可是邻域平均法却在去处奇斯噪声方而有比较好的效果,但是用它处理椒盐噪声效果却并不理想,虽然降低了噪声可也使图像出现了模糊,而且这种模糊在边缘和细节岀特別严重。
空间域低通滤波法可以说是对邻域平均法的一种改进,如果能选择比较合适的单位冲激响应阵列,那么在达到图像平滑的同时还可以很好的保留图像细节,可是如何选择单位冲激响应阵列却是一个难点。
多幅图像平均法一般用于摄像机的视频图像,用以减少电视摄像机光电摄像管或CCD器件所引起的噪声,该算法有比较好的去噪效果,但是在实际使用时,该算法要用到多幅图像,所以占用的空间也比较大,尤其是把多幅图像配准也比较难。
梯度倒数加权算法能够在降低噪声的同时,较好地保持图像的边缘和细节信息,然而,对于被椒盐噪声或脉冲噪声污染的图像,梯度倒数加权平滑算滤波效果并不明显,而且该方法的汁算量也比前而的几种算法要大。
因此,在对一幅图像进行平滑处理前,
必须仔细分析其产生噪声的原因、噪声的特点与类型,并选择合适的平滑方法,才能既消除图像噪声,又不使图像边缘轮廉或线条变模糊。
经过这样的处理后,图像才能更符合人的视觉特性。
参考文献
[1]孚弼程,彭天强,彭波等编著•智能图像处理技术㈣]•北京:
电子匸业出版社,2004.⑵阮秋琦•数字图像处理学[M]•北京:
电子工业岀版社,2003.
⑶徐飞,施晓红.Matlab应用图像处理[M]•西安:
西安电子科技大学出版社,2002.
⑷陈传波,金先级•数字图像处埋[M]・北京:
机械12业出版社,2004.
⑸冯清枝•基于直方图修正的图像増强技术[J]•广东公安科技,2004.(9).
⑹陈书海,傅录祥•实用数字图像处理[M]・北京:
科学出版社,2005.
⑺沈志光数字图像平滑算法的探尤[J]•科学大众,2006.(4).
⑻李世进•数字图像的平滑处理[J]・湖南科技学院学报,2008.(12).
⑼杨乐•图像増强算法及其实现[J]・现代电子技术,2007.(16).
(注:
可编辑下载,若有不当之处,请指正,谢谢!
)