第六章1清华大学数字图像处理课件

第6章图像增强图像噪声图像增强方法分类基于灰度变换的图像增强方法基于直方图的图像增强方法空间域滤波增强技术频域滤波增强彩色增强第六章图像增强第6章图像增强图像噪声图像增强是数字图像处理的基本内容之一，其目的是根据应用需要突出图像中的某些“有用”信息，削弱或去除不需要的信息，以达到扩大图像中不同物体特征之间的差别。图像增强算法并不能增加原始图像的信息，而是通过某种技术手段有选择地突出对某一具体应用有价值的信息。第6章图像增强对于数字图像处理而言，噪声是指图像中的非本源信息。因此，噪声会影响人的感官对所接收的信源信息的准确理解。第6章图像增强1.图像噪声的产生在图像形成过程中，图像数字化设备、电气系统和外界影响将使得图像噪声的产生不可避免。第6章图像增强2.图像噪声分类图像噪声按其产生的原因可分为外部噪声和内部噪声。一般情况下，数字图像中常见的外部干扰主要包括如下几种：设备元器件及材料本身引起的噪声。如磁带、磁盘表面缺陷所产生的噪声;系统内部设备电路所引起的噪声，包括电源系统引入的交流噪声，偏转系统和箝位电路引起的噪声等;电器部件机械运动产生的噪声。如数字化设备的各种接头因抖动引起的电流变化所产生的噪声，磁头、磁带抖动引起的抖动噪声等;第6章图像增强3.图像噪声特点叠加性；分布和大小的不规则性；噪声与图像之间具有相关性；第6章图像增强图像增强方法分类图像增强的定义图像增强的目的图像增强的技术方法第6章图像增强图像增强的定义按我们的需要突出一幅图像中的某些“有用”信息，同时削弱另外一些“无用”信息的图像处理方法。返回第6章图像增强图像增强（突出边界）第6章图像增强图像增强（消除噪声）返回第6章图像增强图像增强的目的通过对图像的处理，使图像比处理前更适合一个特定的应用。也就是为了某种应用目的去改善图像的质量。返回1.改善图像的视觉效果。2.突出图像的特征，便于计算机处理。第6章图像增强图像增强的技术方法主要有空域处理法和频域处理法（1）空域处理法：直接在图像所在的二维空间进行处理，即直接对每一像元的灰度值进行处理。（2）频域处理法：将图像从空间域变换到频率域对图像进行处理。第6章图像增强空间域图像增强频率域灰度变换空域滤波直接灰度变换直方图修正法图像的代数运算直方图均衡化直方图规定化图像平滑图像锐化高通滤波低通滤波带通、带阻滤波图6.1图像增强的主要内容返回第6章图像增强基于灰度变换的图像增强方法基于空间域灰度变换的图像增强方法是一种点处理方法。点处理实际上是灰度到灰度的映射过程，设输入图像每个像元的灰度值为A(x,y)，输出图像的灰度值为B(x,y)，表示为：B(x,y)=f[A(x,y)]，显然点运算不会改变图像内像素点之间的空间关系。第6章图像增强基于空域灰度变换的图像增强方法直接灰度变换基于图像的直方图进行灰度变换对两个以上的图像进行代数运算实现灰度变换返回第6章图像增强直接灰度变换直接灰度变换非线性拉伸线性拉伸按比例线性拉伸分段线性拉伸对数扩展指数扩展第6章图像增强线性拉伸：是将输入图像的灰度值的动态范围按线性关系公式拉伸扩展到指定范围。包括按比例线性拉伸和分段线性拉伸。（1）按比例线性拉伸:将原始图像的灰度范围不加区别的扩展。fgaba’b’255218482550第6章图像增强第6章图像增强提高对比度举例第6章图像增强灰度动态范围较窄观察直方图分布第6章图像增强对比度拉伸灰度动态范围变宽第6章图像增强灰度动态范围变宽观察直方图分布第6章图像增强第6章图像增强25548255019621623第6章图像增强非线性拉伸：是在整个灰度值范围内按照统一的变换函数来实现对不同灰度值区间的扩展与压缩。常用的有按对数函数扩展和按指数函数扩展。按对数函数变换：低灰度区扩展，高灰度区压缩。第6章图像增强按指数函数变换：高灰度区扩展，低灰度区压缩。第6章图像增强原始图象第6章图像增强非线性灰度变换对数效应第6章图像增强返回非线性灰度变换指数效应第6章图像增强基于直方图的图像增强方法灰度直方图：数字图像中每一灰度级与它出现的频率之间的统计，可以理解为描述各个灰度级的像素出现多少的统计图示。若用横坐标表示灰度级，纵坐标表示频率，就可以看出图像中灰度的分布情况，例：水泥微观结构图与左图对应的直方图第6章图像增强第6章图像增强较暗图象的直方图P(rk)rk第6章图像增强较亮图象的直方图P(rk)rk第6章图像增强对比度较低图象的直方图P(rk)rk第6章图像增强基于灰度直方图的图像增强直方图均衡化直方图规定化目的原理及基本算法应用举例目的原理及处理步骤应用举例第6章图像增强直方图均衡化1)目的将原始图像的直方图变为均衡分布的形式，即将一已知灰度概率密度分布的图像，经过某种变换，变成一幅具有均匀灰度概率密度分布的新图像。图像均衡化处理后，图像的直方图是平直的，即各灰度级具有相同的出现频数，那么由于灰度级具有均匀的概率分布，图像看起来就更清晰了。第6章图像增强对比度较高图象的直方图P(rk)rk第6章图像增强第6章图像增强第6章图像增强第6章图像增强第6章图像增强第6章图像增强第6章图像增强直方图均衡化2.原理首先假定连续灰度级的情况，推导直方图均衡化变换公式，令r代表灰度级，P(r)为概率密度函数。r值已归一化，最大灰度值为1。第6章图像增强连续灰度的直方图非均匀分布第6章图像增强连续灰度的直方图均匀分布第6章图像增强直方图均衡化目标直方图均衡化第6章图像增强要找到一种变换S=T(r)使直方图变平直，为使变换后的灰度仍保持从黑到白的单一变化顺序，且变换范围与原先一致，以避免整体变亮或变暗。必须规定：（1）在0≤r≤1中，T(r)是单调递增函数，且0≤T(r)≤1；（2）反变换r=T-1(s),T-1(s)也为单调递增函数，0≤s≤1。第6章图像增强rjrj+rsjsj+s直方图均衡化变换公式推导图示第6章图像增强考虑到灰度变换不影响像素的位置分布，也不会增减像素数目。所以有)(1)()(000rTsdsdsspdrrprssrdrrprT0)12()()(第6章图像增强应用到离散灰度级，设一幅图像的像素总数为n，分L个灰度级。nk:第k个灰度级出现的个数。第k个灰度级出现的概率：P(rk)=nk/n其中0≤rk≤1，k=0,1,2,...,L-1形式为：)22()()(00kjkjjjkknnrprTs第6章图像增强3.基本步骤:(1)求出图像中所包含的灰度级rk,可以定为0~L-1,(2)统计各灰度级的像素数目nk(k=0,1,2,…L-1)(3)计算图像直方图(4)计算变换函数：(5)用变换函数计算映射后输出的灰度级Sk(6)统计映射后新的灰度级Sk的像素数目nk(7)计算输出图像的直方图第6章图像增强例例：设图象有64*64=4096个象素，有8个灰度级，灰度分布如表所示。进行直方图均衡化。rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1nk790102385065632924512281P(rk)0.190.250.210.160.080.060.030.02第6章图像增强步骤：例rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1nk790102385065632924512281P(rk)0.190.250.210.160.080.060.030.02第6章图像增强(1)由（2-2）式计算skSk计算0.190.440.650.810.890.950.981.00例rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1nk790102385065632924512281P(rk)0.190.250.210.160.080.060.030.02第6章图像增强Sk舍入1/73/75/76/76/7111(2)把计算的sk就近安排到8个灰度级中。例Sk计算0.190.440.650.810.890.950.981.00rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1nk790102385065632924512281P(rk)0.190.250.210.160.080.060.030.02第6章图像增强Sks0s1s2s3s4nsk7901023850985448P(sk)0.190.250.210.240.11(3)重新命名sk，归并相同灰度级的像素数。例Sk计算0.190.440.650.810.890.950.981.00rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1nk790102385065632924512281P(rk)0.190.250.210.160.080.060.030.02Sk舍入1/73/75/76/76/7111第6章图像增强直方图均衡化均衡化前后直方图比较例第6章图像增强直方图均衡化直方图均衡化实质上是减少图像的灰度级以换取对比度的加大。在均衡过程中，原来的直方图上频数较小的灰度级被归入很少几个或一个灰度级内，故得不到增强。若这些灰度级所构成的图像细节比较重要，则需采用局部区域直方图均衡。第6章图像增强第6章图像增强第6章图像增强第6章图像增强第6章图像增强第6章图像增强原图像及直方图均衡后的图像及直方图图像的反差大了，细节清楚了在直方图中的表现是直方图灰度范围窄且集中在低灰度值区域。现在直方图占据了整个图像灰度值的允许范围，增加了图像的动态范围。原图较暗且动态范围小第6章图像增强第6章图像增强(2)用变换函数计算映射后输出的灰度级:原图像的灰度只有8级,所以Sk需以1/7为量化单位进行舍入运算。（1/7=0.142/7=0.293/7=0.434/7=0.575/7=0.726/7=0.867/7=1）S0=0.19S0=1/7S1=0.44S1=3/7S2=0.65S2=5/7S3=0.81S3=6/7S4=0.89S4=6/7S5=0.95S5=1S6=0.98S6=1S7=1S7=1第6章图像增强(3)统计映射后各灰度级的像素数目ni：由上舍入结果可见，均衡化后的灰度级仅有5级，分别是S0=1/7；S1=3/7；S2=5/7；S3=6/7；S4=1对应的像素数目是S0=1/7n0=790S1=3/7n1=1023S2=5/7n2=850S3=6/7n3+n4=656+329=985S4=1n5+n6+n7=245+122+81=448第6章图像增强(4)计算输出图像的直方图:第6章图像增强直方图均衡化技术的缺点：（1）将原始函数的累积分布函数作为变换函数，只能产生近似均匀的直方图。（2）在某些情况下，并不一定需要具有均匀直方图的图像。返回第6章图像增强(课堂作业)一幅图像共有8个灰度级，每一灰度级概率分布如下表所示，要求对其进行直方图均衡化处理，并画出均衡化后的图像的直方图。rkr0r1r2r3r4r5r6r7Pr（rk）0.290.240.170.120.090.060.020.01第6章图像增强SkPs(Sk)1/72/73/74/75/76/7010.290.240.170.210.09第6章图像增强直方图规定化1.定义修改一幅图像的直方图，使它与另一幅图像的直方图匹配或具有一种预先规定的函数形状。目标：当需要具有特定的直方图的图像时，可按照预先设定的某个形状人为的调整图像的直方图。第6章图像增强连续灰度的直方图原图第6章图像增强连续灰度的直方图规定第6章图像增强2.直方图规定化的原理Pr(r)表示原始图像的灰度分布概率密度函数。Pz(z)表示目的图像的灰度分布概率密度函数。（1）首先对原始图像进行直方图均衡化处理，即求变换函数：（2）对目标图像也进行直方图均衡化处理，即：rdwwrTs0)Pr()(zdwwPzzGv0)()(其逆变换是即由均衡化后的灰度级v得到目标函数的灰