《数字图像处理与分析第3章.ppt》由会员分享，可在线阅读，更多相关《数字图像处理与分析第3章.ppt（100页珍藏版）》请在三一办公上搜索。
1、3.1 概述 3.2 空域增强 3.3 频域增强 3.4 图像的锐化 3.5 彩色图像增强 3.6 小结,第3章 图像增强,问题的引入,看两个图例，分析画面效果不好的原因,亮暗差别不是很大,解决问题的思路,提高对比度，增加清晰度,对比度的概念,对比度通俗地讲，就是亮暗的对比程度。对比度通常表现了图像画质的清晰程度。,问题1：灰度分布不合理,没有充分利用灰度动态范围,典型场合:曝光不足、曝光过度、对比过于强烈,问题2：噪声干扰,原因：强噪声成像通道,问题3：图像模糊,影响图像细节分辨,原因：成像通道分辨率不足、景物移动等,3.1 概述:,定义:图像增强是通过某种技术有选择地突出对某一具体应用有用
2、的信息，削弱或抑制一些无用的信息。,结果：改善后的图像不一定逼近原图像,前提：不考虑图像降质的原因,3、图像增强处理最大的困难增强后图像质量的好坏主要依靠人的主观视觉来评定，也就是说，难以定量描述,注意：,1、图像增强处理并不能增加原始图像的信息，其结果只能增强对某种信息的辨别能力，而这种处理肯定会损失一些其它信息,2、强调根据具体应用而言，更“好”，更“有用”的视觉效果图像,3.1 图像增强的两类方法,空域增强:直接在图像所在的二维空间进行处理，即直接对每一像素的灰度值进行处理。频域增强:首先经过傅里叶变换将图像从空间域变换到频率域，然后在频率域对频谱进行操作和处理，再将其反变换到空间域，从
3、而得到增强后的图像。,图象增强方法总结,3.1 概述空域增强,按技术不同可分为灰度变换和空间滤波灰度变换基于点操作，将每一个像素的灰度值按照一定的数学变换公式转换为一个新的灰度值。常用的有:对比度增强、直方图变换等方法。空域滤波基于邻域处理，应用某一模板对每个像素及其周围邻域的所有像素进行某种数学运算，得到该像素的新的灰度值。图像平滑与锐化技术就属于空域滤波。,3.2.1 灰度变换法增强,灰度变换可使图像对比度扩展，图像清晰，特征明显。它是图像增强的重要手段。,1.线性灰度变换。,2.分段线性灰度变换。,3.非线性灰度变换。,当图像成象时曝光不足或过度,或由于成象设备的非线性和图像记录设备动态
4、范围太窄等因素。都会产生对比度不足的弊病，使图像中的细节分辨不清。这时可将灰度范围线性扩展。设f(x,y)灰度范围为a,b，g(x,y)灰度范围为c,d,(1)线性灰度变换,线性拉伸前:图像灰度集中在a,b之间.线性拉伸后:图像灰度集中在c,d之间.,图像灰度变换前后效果对比图:变换前 变换后,函数imadjust介绍,Imadjust(f,low_in high_in,low_out high_out,gamma)对图像进行灰度变换所有输入输出都指定在0和1之间gamma指定曲线的形状 Gamma=1是线性变换 Gamma1高灰度区扩展，低灰度区压缩。,使用imadjust函数的灰度变换步骤
5、,(1)绘制直方图，观察灰度范围(2)将灰度范围转换为0-1之间的数，使得灰度范围low_in,high_in传递给imadjust函数,线性拉伸将原始输入图像中的灰度值不加区别地扩展。,(2)分段线性变换,分段线性拉伸是将某一段感兴趣的灰度范围线性扩展，相对抑制不感兴趣的灰度区域。,0,f(x,y),g(x,y),a,b,a,b,M,M,分段线性变换,分段线性变换,例：,原始图像 分段线性变换函数图 线性变换后的图像,分段线性变换 灰级窗,当256个灰度级所表示的亮暗范围内的信息量太大，没办法很好地表述时，通过开窗的方式，每次只把窗内的灰度级展宽，而把窗外的灰度级完全抑制掉。例如，CT图像的
6、原始数据为12bit（或是16bit），要将其显示出来，则只能转换为8bit，于是有了16归并为1（或256归并为1）的需要。这时，开骨窗、肌肉窗、组织窗就可分别清晰地显示相应的内容。,灰级窗效果示例,原图,肺窗,肌肉窗,骨窗,分段线性变换 灰级窗的实现方法,如图所示，灰级窗实际上是线性对比度展宽的一种特殊形式。,线性对比度展宽,灰级窗,(3)非线性灰度变换,对数变换幂次变换,对数变换,对数变换,低灰度区扩展，高灰度区压缩。,一个简单的例子,image=imread(moon.tif);subplot(1,2,1);imshow(image);image=double(image);log_i
7、mage=log(image);subplot(1,2,2);imshow(log_image),原图像,对数变换后,指数变换,高灰度区扩展，低灰度区压缩。,实 例,原始图像,灰度倒置底片效果,原始图像,非线性灰度变换对数效应,原始图像,非线性灰度变换指数效应,3.2.2 直方图变换增强,灰度直方图是灰度值的函数，它描述了图像中各灰度值的像素个数。通常用横坐标表示像素的灰度级别，纵坐标表示对应的灰度级出现的频率（像素的个数）。,数字图像的灰度直方图 计算例,灰度直方图,直方图变换增强灰度直方图,直方图的性质只能反映图像的灰度分布情况，而不能反映图像像素的位置,。一幅图像对应唯一的灰度直方图，反
8、之不成立。,不同的图像具有相同直方图,常用的直方图是规格化和离散化的，即纵坐标用相对值表示。设图像总像素为N，某一级灰度像素数为nk，则直方图表示为：p(rk)=nk/N,灰度直方图反映了一幅图像的灰度分布情况。(a)大多数像素灰度值取在较暗区域，图像会较暗.一般在摄影过程中曝光过弱就会造成这种结果。(b)大多数像素灰度值集中在亮区,图像将偏亮.一般在摄影中曝光太强将导致这种结果。(c)图像的像素窄而集中，对比度低。从三幅图像的灰度分布来看图像的质量均不理想。,注意高对比度的图像有更平坦的直方图。一幅图像应该利用全部或几乎全部可能的灰度级,灰度直方图的又一应用 分割阈值选取,假设某图像的灰度直
9、方图具有 二峰性，则表明这个图像较亮的区域和较暗的区域可以较好地分离。取二峰间的谷点为阈值点，可以得到好的二值处理的效果。,具有二峰性的灰度图的二值化,(2)直方图均衡化,直方图均衡方法的基本思想是，对在图像中像素个数多的灰度级进行展宽，而对像素个数少的灰度级进行缩减。从而达到清晰图像的目的。,直方图均衡化,s=T(r)r代表原始图像的灰度级，s为变换后的灰度级。通过上述变换，每个原始图像的像素灰度值r都对应产生一个s值。,连续灰度的直方图非均匀分布,连续灰度的直方图均匀分布,直方图均衡化目标,直方图均衡化,直方图均衡化,要找到一种变换 S=T(r)使直方图变平直，为使变换后的灰度仍保持从黑到
10、白的单一变化顺序，且变换范围与原先一致，以避免整体变亮或变暗。必须规定：（1）在0r1中，T(r)是单调递增函数，且0T(r)1；（2）反变换r=T-1(s),T-1(s)也为单调递增函数，0s1。,直方图均衡化,考虑到灰度变换不影响象素的位置分布，也不会增减像素数目。所以有,直方图均衡化,应用到离散灰度级，设一幅图像的象素总数为n，分L个灰度级。nk:第k个灰度级出现的频数。第k个灰度级出现的概率 P(rk)=nk/n 其中0rk1，k=0,1,2,.,L-1 形式为：,直方图均衡化的步骤,1、计算每个灰度级的像素个数在整个图像中所占的概率(百分比),2、计算图像各灰度级的累计分布概率,直方
11、图均衡化的步骤,3、根据 的值判断变换后的灰度值 设图像的灰度级只有8级，因此需用1/7为量化单位进行舍入运算，得到如下结果：,的值落到的哪个区间,则对应变换到该灰度值,例,例：设图像有64*64=4096个象素，有8个灰度级，灰度分布如表所示。进行直方图均衡化。,rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1,nk 790102385065632924512281,p(rk),例,例：设图像有64*64=4096个象素，有8个灰度级，灰度分布如表所示。进行直方图均衡化。,rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6
12、=6/7r7=1,nk 790102385065632924512281,p(rk)0.190.250.210.160.080.060.030.02,1.由（2-2）式计算sk。,rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1,nk 790102385065632924512281,p(rk)0.190.250.210.160.080.060.030.02,sk计算 0.190.440.650.810.890.950.981.00,例,步骤：,rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1,nk 790
13、102385065632924512281,p(rk)0.190.250.210.160.080.060.030.02,sk计算 0.190.440.650.810.890.950.981.00,sk舍入 1/73/75/76/76/7111,2.把计算的sk就近安排到8个灰度级中。,例,rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1,nk 790102385065632924512281,p(rk)0.190.250.210.160.080.060.030.02,sk计算 0.190.440.650.810.890.950.981.00,sk舍入
14、 1/73/75/76/76/7111,sk s0s1s2s3s4,nsk 7901023850985448,p(sk)0.190.250.210.240.11,3.重新命名sk，归并相同灰度级的象素数。,例,直方图均衡化,均衡化前后直方图比较,例,均衡化,直方图均衡化效果示例,直方图规定化直方图均衡化存在的问题,直方图均衡化的优点是得到近似均匀分布的直方图。但由于变换函数采用累积分布函数，只能产生近似均匀的直方图的结果,实际应用中，有时需要具有特定直方图的图像，以便能够有目的地对图像中的某些灰度级分布范围内的图像加以增强。,例：图像均衡化效果,原图像,均衡化后的图像,直方图规定化的思想,直方
15、图规定化方法是使原图像灰度直方图变成规定形状的直方图而对图像作修正的增强方法。,Matlab函数,imhist（f，n）函数：计算和显示图像的直方图。n为指定的灰度级数目，缺省值为256。g=histeq(f,nlev)%f为输入图像，nlev是输出图像的灰度级数，默认值为64，通常我们设置为256。g=histeq(f,hspec)%f为输入图像，hspec为指定的直方图(一个由指定值构成的行向量)。,3.2.3 空间平滑滤波增强,背景 图像在传输过程中，由于传输信道、采样系统质量较差，或受各种干扰等影响，会造成图像毛糙，此时，就需对图像进行平滑处理。,图像平滑的作用类似剃须刀,平滑可以去除
16、毛糙，噪声，但也使图像变得模糊。,1.邻域平均法:线性滤波(均值滤波)2.中值滤波：非线性滤波,假设图像由许多灰度恒定的小块组成，相邻像素间存在很高的空间相关性，而噪声则相对独立。可以将一个像素及其邻域内的所有像素的平均灰度值赋给平滑图像中对应的像素，从而达到平滑的目的。（常用的邻域有4-邻域和8-邻域）作用：减噪，去除不相干的细节，对灰度级不足引起的伪轮廓进行平滑等等。,1.邻域平均法,1.邻域平均法-非加权邻域平均,最简单的邻域平均法为非加权邻域平均：一幅图像大小为NN的图像f(x,y)，用邻域平均法得到的平滑图像为g(x,y)，则 x,y=0,1,N-1;s为(x,y)邻域中像素坐标的集
17、合，其中不包括(x,y);M表示集合s内像素的总数。,3.2.3 空间平滑滤波增强非加权邻域平均,像素相邻：四连接：当前像素为黑，其四个近邻像素中至少有一个为黑；八连接：当前像素为黑，其八个近邻像素中至少有一个为黑。,3.2.3 空间平滑滤波增强非加权邻域平均,在图像上，对待处理的像素给定一个模板，该模板包括了其周围的邻近像素。将模板中的全体像素的均值来替代原来的像素值的方法。,1.邻域平均法-非加权邻域平均,非加权邻域平均法可以用模板求得，即在待处理图像中逐点地移动模板，求模板系数与图像中相应像素的乘积之和，模板数为1。下图是非加权邻域平均33模板。,非加权邻域平均33模板：,3,4,4,4
18、,5,6,6,7,8,原图像,处理后的图像,非加权邻域均值滤波器 处理方法,待处理像素,3.2.3 空间平滑滤波增强非加权邻域平均,边界处理：,3.2.3 空间平滑滤波增强非加权邻域平均,边界处理：,（a）为含有随机噪声的灰度图像(b)（c）(d）是分别用33、55、77模板得到的平滑图像。,均值滤波器的缺点是，会使图像变的模糊，原因是它对所有的点都是同等对待，在将噪声点分摊的同时，将景物的边界点也分摊了。为了改善效果，就可采用加权平均的方式来构造滤波器。,1.邻域平均法-加权邻域平均,3.2.3 空间平滑滤波增强加权邻域平均,所有模板系数可以有不同的权值(p61),3.2.3 空间平滑滤波增
19、强加权邻域平均,下面几个典型的加权平均滤波器。,邻域平均法虽然可以平滑图像，但在消除噪声的同时，会使图像中的一些细节变得模糊。中值滤波则在消除噪声的同时还能保持图像中的细节部分，防止边缘模糊。,2.中值滤波,中值滤波方法对脉冲干扰和椒盐噪声的抑制效果好，在抑制随机噪声的同时能够保持边缘减少模糊。,中值滤波器 设计思想,因为噪声（如椒盐噪声）的出现，使该点像素比周围的像素亮（暗）许多。如果在某个模板中，对像素进行由小到大排列的重新排列，那么最亮的或者是最暗的点一定被排在两侧。对于图像中的每个像素,先确定一个奇数像素窗口W，窗口内各像素按灰度值从小到大排序后，用中间位置灰度值代替原灰度值，就可以达
20、到滤除噪声的目的。,中值滤波器 例题,2,3,4,5,6,6,6,7,8,（a）为含有随机噪声的灰度图像(b)（c）(d）是分别用33、55、77模板得到的平滑图像。,可以看出，中值滤波的效果优于均值滤波的效果，图像中的边缘轮廓比较清晰,中值滤波器与均值滤波器的比较,对于椒盐噪声，中值滤波效果比均值滤波效果好。,中值滤波器与均值滤波器的比较,原因：椒盐噪声是幅值近似相等但随机分布在不同位置上，图像中有干净点也有污染点。中值滤波是选择适当的点来替代污染点的值，所以处理效果好。因为噪声的均值不为0，所以均值滤波不能很好地去除噪声点。,中值滤波器与均值滤波器的比较,对于高斯噪声，均值滤波效果比中值滤波效果好。,中值滤波器与均值滤波器的比较,原因：高斯噪声是幅值近似正态分布，但分布在每点像素上。因为图像中的每点都是污染点，所以中值滤波选不到合适的干净点。因为正态分布的均值为0，所以均值滤波可以消除噪声。（注意：实际上只能减弱，不能消除。）,线性滤波函数,g=imfilter(f,w,filtering_mode,boundary_options,size_options)其中：f为输入图像，w为滤波掩模。w可通过fspecial函数生成。,非线性滤波器,g=medfilt2(f,m n)中值滤波器,椒盐噪声示例,高斯噪声示例,
