数字图像处理第4章直方图和灰度变换

1、第四章 图像增强与平滑 第四章 图像增强与平滑 4.1 背景知识背景知识4.2 直方图直方图4.3 灰度变换灰度变换4.4 图像噪声图像噪声4.5 去除噪声去除噪声4.6 图像锐化图像锐化4.7 图像的伪彩色处理图像的伪彩色处理第四章 图像增强与平滑 第四章 图像增强与平滑 图像处理中最具有吸引力的领域之一图像处理中最具有吸引力的领域之一 q增强的增强的首要目标首要目标是处理图像是处理图像, ,使其比原始图像更适合于使其比原始图像更适合于特定特定应用。其技术主要包括直方图修改处理、图像平应用。其技术主要包括直方图修改处理、图像平滑、图像锐化及彩色处理等。滑、图像锐化及彩色处理等。q图像增强的图

2、像增强的方法方法分为两大类：空间域方法和频域方法。分为两大类：空间域方法和频域方法。q当图像为视觉解释而进行处理时，由观察者最后判断当图像为视觉解释而进行处理时，由观察者最后判断特定方法的效果。图像质量的特定方法的效果。图像质量的视觉评价视觉评价是一种高度主是一种高度主观的过程，因此，定义一个观的过程，因此，定义一个“理想图像理想图像”标准，通过标准，通过这个标准去比较算法的性能。这个标准去比较算法的性能。第四章 图像增强与平滑 4.1 4.1 背景知识背景知识 “空间域增强空间域增强”是指增强构成图像的像素。是指增强构成图像的像素。空间域方法是直接对这些像素操作的过程。空间域方法是直接对这些

3、像素操作的过程。空间域处理可由下式定义：空间域处理可由下式定义：第四章 图像增强与平滑 4.2 4.2 直直 方方 图图 4.2.1 基本概念基本概念概率密度函数（概率密度函数（probability density function (pdf) ）：）：设设r表示图像中像素灰度级，它可看作是一个随机变量。表示图像中像素灰度级，它可看作是一个随机变量。作归一化处理后，被限定在作归一化处理后，被限定在0, 1之内。假定对每一之内。假定对每一瞬间它是连续的随机变量，那么就可以用瞬间它是连续的随机变量，那么就可以用pr(r)来表示来表示原始图像的灰度分布。原始图像的灰度分布。第四章 图像增强与平滑

4、概率密度曲线：概率密度曲线：用直角坐标系的横轴代表用直角坐标系的横轴代表灰度级灰度级r，纵轴代表灰度级的概率密度函数，纵轴代表灰度级的概率密度函数pr(r)，可作出一条曲线。这条曲线在概率，可作出一条曲线。这条曲线在概率论中就是概率密度曲线。论中就是概率密度曲线。第四章 图像增强与平滑 图像灰度分布的概率密度函数图像灰度分布的概率密度函数 pr(r)r10pr(r)10r(a)(b)第四章 图像增强与平滑 直方图直方图 (histogram)： 灰度直方图是灰度级的函数，它表示图像中具有某灰度直方图是灰度级的函数，它表示图像中具有某种灰度级的像素的个数，反映了图像中每种灰度出现种灰度级的像素的

5、个数，反映了图像中每种灰度出现的频率。灰度直方图的横坐标是灰度级，纵坐标是该的频率。灰度直方图的横坐标是灰度级，纵坐标是该灰度级出现的频度，它是图像最基本的统计特征。灰度级出现的频度，它是图像最基本的统计特征。 第四章 图像增强与平滑 图像灰度直方图图像灰度直方图 6646313266416665436646611223466543211426545654321第四章 图像增强与平滑 4.2.2 直方图的性质直方图的性质 （1）直方图只包含了图像中某一灰度值的像素出现）直方图只包含了图像中某一灰度值的像素出现的概率信息，而丢失了其所在位置的信息。的概率信息，而丢失了其所在位置的信息。 （2）图

6、像与直方图之间是多对一的映射关系。）图像与直方图之间是多对一的映射关系。第四章 图像增强与平滑 图像与直方图间的多对一关系图像与直方图间的多对一关系 (a)(b)第四章 图像增强与平滑 （3） 由于直方图是对具有相同灰度值的像素统计得到的，由于直方图是对具有相同灰度值的像素统计得到的， 因此，一幅图像各子区的直方图之和就等于该图像全图的直方因此，一幅图像各子区的直方图之和就等于该图像全图的直方图。图。 直方图的分解直方图的分解 (a)(b)(c)第四章 图像增强与平滑 4.2.3 直方图的计算直方图的计算 在离散形式下，在离散形式下， 用用rk代表离散灰度级，用代表离散灰度级，用pr(rk)代

7、表代表pr(r)，并，并且有下式成立：且有下式成立： nnrpkkr)(1, 2 , 1 , 010lkrk式中：式中：nk为图像中出现为图像中出现rk级灰度的像素数，级灰度的像素数，n是图像像素总数，而是图像像素总数，而nk/n即为频数。在直角坐标系中做出即为频数。在直角坐标系中做出rk与与pr(rk)的关系图形，即称的关系图形，即称为该图像的直方图。为该图像的直方图。第四章 图像增强与平滑 lena图像及直方图图像及直方图（a） lena图像；（图像；（b） lena图像的直方图图像的直方图 第四章 图像增强与平滑 钟楼图像及直方图钟楼图像及直方图（a）钟楼图像；（）钟楼图像；（b）钟楼图

8、像的直方图）钟楼图像的直方图 第四章 图像增强与平滑 4.2.4 直方图修改技术的基础直方图修改技术的基础 一幅给定图像的灰度级分布在一幅给定图像的灰度级分布在0r1范围内，可以对范围内，可以对0, 1区间内的任一个区间内的任一个r值进行如下变换：值进行如下变换： s=t(r) 通过上述变换，每个原始图像的像素灰度值通过上述变换，每个原始图像的像素灰度值r都对应产生一个都对应产生一个s值。变换函数值。变换函数t(r)应满足下列条件：应满足下列条件： （1） 在在0r1区间内，区间内，t(r)值单调增加；值单调增加； （2） 对于对于0r1， 有有 0t(r)1。 第四章 图像增强与平滑 单调递

9、增的灰度级变换函数单调递增的灰度级变换函数 s1skr1rko)(krt第四章 图像增强与平滑 由由概率论理论概率论理论可知，如果可知，如果pr(r) 和变换函数和变换函数t(r)已知，且已知，且t(r)满足条件（满足条件（1），则变换变量），则变换变量s的概率密度函数的概率密度函数ps (s) 可由以下公可由以下公式得到：式得到： 由上式可见，通过变换函数由上式可见，通过变换函数t(r)可以控制图像灰度级的概率可以控制图像灰度级的概率密度函数改变图像的灰度层次。这就是直方图修改技术的理论密度函数改变图像的灰度层次。这就是直方图修改技术的理论基础。基础。 )(11)()()()(strrrsd

10、sdrrpstdsdrpsp第四章 图像增强与平滑 4.2.5 直方图均衡化处理直方图均衡化处理 直方图均衡化处理是以累积分布函数变换法为基直方图均衡化处理是以累积分布函数变换法为基础的直方图修正法。假定变换函数为础的直方图修正法。假定变换函数为 dprtsrr)()(0 式中：式中：是积分变量，是积分变量， 是是r的累积分布的累积分布函数。函数。dprr0)(第四章 图像增强与平滑 直方图均衡化处理直方图均衡化处理 pr(r)21s t(r)pr(s)12rr0 0.2 0.4 0.6 0.8 1.00.20.40.60.81.0r0 0.2 0.4 0.6 0.8 1.00.20.40.6

11、0.81.0(a)(b)(c)0第四章 图像增强与平滑 当灰度级是离散值时，当灰度级是离散值时， 1, 1 , 010)(lkrnnrpkkkr式中：式中：l是灰度级的总数目，是灰度级的总数目，pr(rk)是取第是取第k级灰度值的概率，级灰度值的概率，nk是图像中出现第是图像中出现第k级灰度的次数，级灰度的次数，n是图像中像素总数。则：是图像中像素总数。则：1, 1 , 010)()(00lkrrpnnrtsjkjjrkjjkk其反变换式为其反变换式为 )(1kkstr第四章 图像增强与平滑 直方图均衡化处理举例：书直方图均衡化处理举例：书p58 第四章 图像增强与平滑 经直方图均衡化后的经直

12、方图均衡化后的lena图像及直方图图像及直方图（a） 经直方图均衡化后的经直方图均衡化后的lena图像；图像； （b） 均衡化后的均衡化后的lena图像的直方图图像的直方图 第四章 图像增强与平滑 第四章 图像增强与平滑 直方图均衡化处理的特征：直方图均衡化处理的特征： 自动化处理；自动化处理； 图像动态范围增加；图像动态范围增加； 灰度简并现象。灰度简并现象。第四章 图像增强与平滑 4.2.6 直方图规定化处理（直方图匹配）直方图规定化处理（直方图匹配）算法来源背景：算法来源背景： 直方图均衡化的缺陷：直方图均衡化的缺陷：不能用于交互方式不能用于交互方式的图象增强应用，因为直方图均衡化只能的

13、图象增强应用，因为直方图均衡化只能产生唯一结果。产生唯一结果。 希望通过一个指定的函数（如高斯函数）希望通过一个指定的函数（如高斯函数）或用交互图形产生一个或用交互图形产生一个特定特定的直方图。的直方图。第四章 图像增强与平滑 算法思想：算法思想：设：设：rk是原图象的灰度级是原图象的灰度级,zk是符合指定直方图结果图象的灰度级是符合指定直方图结果图象的灰度级我们的目标是：找到一个灰度级变换我们的目标是：找到一个灰度级变换h，有：，有： z = h(r)第四章 图像增强与平滑 算法步骤：算法步骤：1) 对对rk、zk分别做直方图均衡化分别做直方图均衡化 2) 求求g变换的逆变换变换的逆变换 z

14、 = g-1 (v) dprtsrr)()(0dpzgvzz)()(0)(rp为原始图像灰度分布的为原始图像灰度分布的pdf， 为希望得到的为希望得到的pdf。)(zp第四章 图像增强与平滑 3) 根据均衡化的概念，根据均衡化的概念，s，v都是常量都是常量 用用 s替代替代v有有 z = g-1 (s) 2) 求求g-1和和t的复合变换，有：的复合变换，有： z = g-1 (t(r) = g-1t(r) h = g-1t第四章 图像增强与平滑 算法实现：算法实现：1）求出灰度级变换）求出灰度级变换t 2）求出灰度级变换）求出灰度级变换g，同时求出逆变换，同时求出逆变换g-13）通过）通过t和

15、和g-1求出复合变换求出复合变换h4）用）用h对图象做灰度级变换对图象做灰度级变换第四章 图像增强与平滑 原始直方图原始直方图例：例：第四章 图像增强与平滑 规定化直方图规定化直方图第四章 图像增强与平滑 均衡化处理后的直方图数据均衡化处理后的直方图数据第四章 图像增强与平滑 第四章 图像增强与平滑 第四章 图像增强与平滑 结果直方图数据结果直方图数据第四章 图像增强与平滑 直直方方图图规规定定化化处处理理结结果果图图第四章 图像增强与平滑 第四章 图像增强与平滑 第四章 图像增强与平滑 第四章 图像增强与平滑 4.3 4.3 灰灰 度度 变变 换换 4.3.1 某些基本灰度变换某些基本灰度变

16、换 图像增强常用的三图像增强常用的三种类型函数：线性种类型函数：线性函数（正比函数（正比/反比）、反比）、对数函数（对数对数函数（对数/反反对数）、幂函数对数）、幂函数（n次幂次幂/n次方次方根）。根）。第四章 图像增强与平滑 1. 图像反转图像反转 该处理适用于该处理适用于增强嵌入图像增强嵌入图像暗色区域的白暗色区域的白色或灰色细节，色或灰色细节，特别是当黑色特别是当黑色面积占主导地面积占主导地位时。位时。第四章 图像增强与平滑 2. 对数变换对数变换 这种变换来扩展被这种变换来扩展被压缩的高值图像中压缩的高值图像中的暗像素。频谱值的暗像素。频谱值的范围从的范围从0 0到到10106 6或或

17、更高的情况是常见更高的情况是常见的。当计算机处理的。当计算机处理像这样的无误数字像这样的无误数字时，图像显示系统时，图像显示系统通常不能如实地再通常不能如实地再现如此大范围的强现如此大范围的强度值。度值。第四章 图像增强与平滑 3. 幂次变换幂次变换 随着随着值的变化将简值的变化将简单地得到一族变换单地得到一族变换曲线。如预期的一曲线。如预期的一样，我们看到图中样，我们看到图中11的值和的值和11的的值产生的曲线有相值产生的曲线有相反的效果。用于图反的效果。用于图像获取、打印和显像获取、打印和显示的各种装置根据示的各种装置根据幂次规律进行响应。幂次规律进行响应。第四章 图像增强与平滑 第四章

18、图像增强与平滑 用幂次变换进用幂次变换进行对比度增强。行对比度增强。相应的伽马值相应的伽马值分别为分别为0.6，0.4和和0.3 (c始终为始终为1) 人人的的脊脊椎椎骨骨折折的的核核磁磁共共振振图图像像第四章 图像增强与平滑 幂次变换幂次变换的另一例的另一例证证 令令= 3.0，4.0和和5.0的处理的处理结果。结果。 第四章 图像增强与平滑 4. 线性变换线性变换 假定原图像假定原图像f(x, y)的灰度范围为的灰度范围为a, b，希望变换后图像，希望变换后图像g(x, y)的灰度范围扩展至的灰度范围扩展至c, d，则线性变换可表示为，则线性变换可表示为 cayxfabcdyxg),(),

19、(第四章 图像增强与平滑 dcayxfabcdcyxg),(),(fmyxfbbyxfaayxf),(),(),(0 有时为了保持f(x, y)灰度低端和高端值不变，可以采用以下形式： ),(),(),(yxfcayxfabcdyxgbyxfa),(式中的a、 b、 c、 d这些分割点可根据用户的不同需要来确定。 第四章 图像增强与平滑 线性灰度变换线性灰度变换(a) 原始图像；原始图像； (b) 灰度变换后的图像灰度变换后的图像 第四章 图像增强与平滑 5. 分段线性变换分段线性变换 为了突出感兴趣的目标或灰度区间，相对抑制那些不感兴为了突出感兴趣的目标或灰度区间，相对抑制那些不感兴趣的灰度区域，可采用分段线性变换。常用的三段线性变换法趣的灰度区域