最新图像增强精品课件.ppt

《最新图像增强精品课件.ppt》由会员分享，可在线阅读，更多相关《最新图像增强精品课件.ppt（71页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、图像增强图像增强4.1 4.1 图像增强概述图像增强概述 4.1.1 4.1.1 图像增强的定义图像增强的定义u对图像的某些特征，如对图像的某些特征，如u边缘边缘u轮廓轮廓u对比度等对比度等u进行强调或锐化，以便于显示、观察或进一进行强调或锐化，以便于显示、观察或进一步分析与处理。步分析与处理。l 可得比例因子：可得比例因子： l 可得实际图像可得实际图像g(i,j)经校正后所恢复的原始图像经校正后所恢复的原始图像 l 乘了一个系数乘了一个系数C/ gc(i,j) ，校正后可能出现，校正后可能出现“溢出溢出”现现象象l 灰度级值可能超过某些记录器件或显示设备输入信号的动态灰度级值可能超过某些记

2、录器件或显示设备输入信号的动态可范围可范围l 需再作适当的灰度变换，最后对变换后的图像进行量化。需再作适当的灰度变换，最后对变换后的图像进行量化。 1,ce ijgij C（4.3）,cgijfijCgij（4.4）4.2.2 灰度变换灰度变换n 灰度变换灰度变换n可使图像动态范围增大，图像对比度扩展可使图像动态范围增大，图像对比度扩展n从而使图像变得清晰以及图像上的特征变得明显。从而使图像变得清晰以及图像上的特征变得明显。1线性变换线性变换n 令原图像令原图像f (i, j)的灰度范围为的灰度范围为a，bn 线性变换后图像线性变换后图像g(i, j)的范围为的范围为a, b。f(i,j)g(

3、i,j)a bba图图4.2 线性变换线性变换g(i, j)与与f(i, j)之间的关系为：之间的关系为： 分段线性变换的目的分段线性变换的目的突出感兴趣的目标或灰度区间，相对抑突出感兴趣的目标或灰度区间，相对抑制那些不感兴趣的灰度区间。制那些不感兴趣的灰度区间。常用的是三段线性变换。常用的是三段线性变换。 ,bagi jafi jaba（4.5）abfcdggMabfM图图4.3 三段线性三段线性l 对灰度区间对灰度区间a,b进行了线性拉伸，而灰度区进行了线性拉伸，而灰度区间间0,a和和b,Mf则被压缩。则被压缩。l 仔细调整折线拐点的位置及控制分段直线的仔细调整折线拐点的位置及控制分段直线

4、的斜率，可以对图像的任一灰度区间进行拉伸斜率，可以对图像的任一灰度区间进行拉伸或压缩或压缩 。 /,/,/,gfcafijgijdcbafijacMdMbfijbd 0,ffi jaafi jbbfi jM 【例例4.1 】在在MATLAB环境中，采用图像线性变环境中，采用图像线性变换进行图像增强。应用换进行图像增强。应用MATLAB的函数的函数imadjust将图像将图像0.32550.7255灰度之间的灰度之间的值通过线性变换映射到值通过线性变换映射到0255之间。之间。 解：解：分别取：分别取：a0.3255，b0.7255， a=0，b=255。 A=imread(pout.tif);

5、 %读入图像读入图像 imshow(A); %显示图像显示图像 figure,imhist(A); %显示图像的直方图显示图像的直方图 J1=imadjust(A,0.3 0.7,); %函数将图像在函数将图像在0.3*2550.7*255灰度之间的值灰度之间的值通过线性变换映射到通过线性变换映射到0255之间之间 figure,imshow(J1); %输出图像效果图输出图像效果图 figure,imhist(J1) %输出图像的直方图输出图像的直方图 实现的程序：实现的程序： （a）原图）原图 （b）原图的直方图）原图的直方图 （c）输出图像）输出图像 （d）输出图像的直方图）输出图像的直

6、方图 图图4.4 图像线性变换图像线性变换2 2非线性灰度变换非线性灰度变换l 当用某些非线性函数如对数、指数函数等当用某些非线性函数如对数、指数函数等作为映射函数时，可实现灰度的非线性变换。作为映射函数时，可实现灰度的非线性变换。l 对数变换的一般表达式为对数变换的一般表达式为：),(1log(),(jifcjig（4.7） l对数变换可以增强低灰度级的像素，压制高对数变换可以增强低灰度级的像素，压制高灰度级的像素灰度级的像素, ,使灰度分布与视觉特性相匹配。使灰度分布与视觉特性相匹配。 直方图（图图4.5 ）：）：指图像中各种不同指图像中各种不同灰度级像素出现的相对频率灰度级像素出现的相对

7、频率 。灰度级灰度级 相对频率相对频率4.2.3 灰度直方图变换灰度直方图变换l 灰度直方图描述了图像的概貌。灰度直方图描述了图像的概貌。l 直方图变换后可使图像的灰度间距拉开或使直方图变换后可使图像的灰度间距拉开或使灰度分布均匀，从而增大对比度，使图像细灰度分布均匀，从而增大对比度，使图像细节清晰，达到增强的目的。节清晰，达到增强的目的。l 直方图变换有直方图变换有l直方图均衡化及直方图规定化两类。直方图均衡化及直方图规定化两类。l 直方图均衡化直方图均衡化l通过对原图像进行某种变换，使得图像的直方图通过对原图像进行某种变换，使得图像的直方图变为均匀分布的直方图。变为均匀分布的直方图。l 灰

8、度级连续的灰度图像：当变换函数是原图灰度级连续的灰度图像：当变换函数是原图像直方图累积分布函数时，能达到直方图均像直方图累积分布函数时，能达到直方图均衡化的目的。衡化的目的。l 对于离散的图像，用频率来代替概率对于离散的图像，用频率来代替概率 。 【例例4.2】假定有一幅总像素为假定有一幅总像素为n6464的图的图像，灰度级数为像，灰度级数为8，各灰度级分布列于表，各灰度级分布列于表4.1中。试对其进行直方图均衡化。中。试对其进行直方图均衡化。 k01234567rk01/72/73/74/75/76/71nk790102385065632924512281pr(rk)0.190.250.21

9、0.160.080.060.030.020.190.440.650.810.890.950.9811/73/75/76/76/71111/73/75/76/717901023850985448pr(sk)0.190.250.210.240.11表表4.1 一幅图像的灰度级分布一幅图像的灰度级分布 ksksksksn 解解：（：（1 1）求变换函数）求变换函数 类似地计算出类似地计算出 ks000019. 0)()(jjrrprTs101144. 025. 019. 0)()(jjrrprTs1,98. 0,95. 0,89. 0,81. 0,65. 0765432ssssss（2）计算）计算l

10、 输出图像灰度是等间隔的，且与原图像灰度输出图像灰度是等间隔的，且与原图像灰度范围一样取范围一样取8个等级，即要求最终的值个等级，即要求最终的值sk=k/7, k=1,2,7。l 需要对需要对 进行重新量化后加以修正进行重新量化后加以修正： ksks1, 1, 1, 7/6,76,75,73, 7176543210ssssssss（3） 的确定的确定u 由由 可知，输出灰度级仅为可知，输出灰度级仅为5个级别：个级别：（4）计算对应每个）计算对应每个 的的 u 因为因为 映射到映射到 ，所以有，所以有790个像个像素在输出输出图像上变成素在输出输出图像上变成ksks1, 7/6, 7/5, 7/

11、3, 7/174210sssssksksn00r 7/10s7/10su 映射到映射到 ，所以有，所以有1023个像素取个像素取值值 。 映射到映射到 ，因此有，因此有850个像素取值个像素取值 。u因为因为 和和 都映射到都映射到 ，因此有，因此有656+329985个像素取值个像素取值 。同理有。同理有245+122+81488个像素变换个像素变换 。 （5）计算）计算 11/7r7/31s7/31s22/7r7/52s7/52s3r4r7/64s7/64s17s/kskspsnn【例例4.3】在在MATLAB环境中，采用直方环境中，采用直方图均衡的方法进行图像增强图均衡的方法进行图像增强

12、 。 解解：程序如下：程序如下 A=imread(p1.jpg);I=histeq(A); %调用函数完成直方图均衡化调用函数完成直方图均衡化subplot(1,2,1),imshow(A); %直方图均衡化前的图像效直方图均衡化前的图像效果果subplot(1,2,2),imshow(I); %直方图均衡化后的图像效果直方图均衡化后的图像效果figure,subplot(1,2,1),imhist(A); %均衡化前的直方图均衡化前的直方图subplot(1,2,2),imhist(I); %均衡化后的直方图均衡化后的直方图4.3 4.3 平滑平滑 l 一种区域增强的算法一种区域增强的算法

13、，平滑算法有：，平滑算法有：l邻域平均法，中值滤波和边界保持类滤波等。邻域平均法，中值滤波和边界保持类滤波等。 4.3.1 4.3.1 邻域平均法邻域平均法n 大部分的噪声都可以看作是随机信号，对图大部分的噪声都可以看作是随机信号，对图像的影响可以看作是孤立的。像的影响可以看作是孤立的。n 某一像素，如果它与周围像素点相比，有明某一像素，如果它与周围像素点相比，有明显的不同，则该点被噪声感染了。显的不同，则该点被噪声感染了。n 设当前待处理像素为设当前待处理像素为f (m,n) ，给出一个大，给出一个大小为小为33的处理模板。的处理模板。图图4.7 模板示意图模板示意图 l处理后的图像设为处理

14、后的图像设为 ，则处理过程，则处理过程可描述为可描述为 l其中其中Z=-1,0,1， 为门限，它可以根据为门限，它可以根据对误差容许的程度，选为图像灰度均方对误差容许的程度，选为图像灰度均方差的差的 若干倍，或者通过实验得到。若干倍，或者通过实验得到。 其他当),(),(91),(),(91),(nmfjnimfnmfjnimfnmgZiZjZiZj),(nmgfu也可以把平均处理看作是图像通过一个低也可以把平均处理看作是图像通过一个低通空间滤波器后的结果通空间滤波器后的结果u设该滤波器的冲激响应为设该滤波器的冲激响应为H(r,s) ，于是滤，于是滤波器输出的结果波器输出的结果g(m,n)表示

15、成卷积的形式，表示成卷积的形式，即即u k,l决定了所选邻域的大小决定了所选邻域的大小 ，为加权函数，为加权函数，又被称为掩模又被称为掩模(Mask)或模板或模板kkrllsNnmsrHsnrmfnmf1, 2 , 1 , 0,),(),(),(（4.8） 1111211111011H1212421211612H111101111813H0010041414141214H常用的模板： I1 = imread(blood1.tif);I=imnoise(I1,salt & pepper); %对图像加椒盐噪声对图像加椒盐噪声imshow(I);h1= 0.1 0.1 0.1; 0.1 0.2 0

16、.1; 0.1 0.1 0.1; %定义定义4种模板种模板h2=1/16.*1 2 1;2 4 2;1 2 1;h3=1/8.*1 1 1;1 0 1;1 1 1;h4=1/2.*0 1/4 0;1/4 1 1/4;0 1/4 0;I2=filter2(h1,I); %用用4种模板进行滤波处理种模板进行滤波处理I3=filter2(h2,I);I4=filter2(h3,I);I5=filter2(h4,I);figure,imshow(I2,) %显示处理结果显示处理结果figure,imshow(I3,)figure,imshow(I4,)figure,imshow(I5,)【例【例4.4

17、4.4】分别采用】分别采用4 4种模板对图像进行处理。种模板对图像进行处理。（a）有噪声的图像）有噪声的图像 （b）模板）模板1处理的结果图处理的结果图 （c）模板）模板2处理的结果图处理的结果图 （d）模板）模板3处理的结果图处理的结果图 （e）模板）模板4处理的结果图处理的结果图图图4.8 平滑处理的例子平滑处理的例子 l 邻域平均法：在去噪的同时也使边界变得邻域平均法：在去噪的同时也使边界变得模糊了。模糊了。l 中值滤波：非线性的处理方法，在去噪的中值滤波：非线性的处理方法，在去噪的同时可以兼顾到边界信息的保留。同时可以兼顾到边界信息的保留。l选一个含有奇数点的窗口选一个含有奇数点的窗口

18、W，将这个窗口在图，将这个窗口在图像上扫描，把该窗口中所含的像素点按灰度级像上扫描，把该窗口中所含的像素点按灰度级的升（或降）序排列，取位于中间的灰度值，的升（或降）序排列，取位于中间的灰度值，来代替该点的灰度值。来代替该点的灰度值。WlklnkmfMediannmg),(),(),(（4.9） 4.3.2 4.3.2 中值滤波中值滤波 例：选择滤波用的窗口例：选择滤波用的窗口W如图，是一个一维的如图，是一个一维的窗口，待处理像素的灰度取这个模板中灰度的窗口，待处理像素的灰度取这个模板中灰度的中值，滤波过程为：中值，滤波过程为： 图图4.9 一维窗口一维窗口 u 常用的窗口还有方形、十字形、圆

19、常用的窗口还有方形、十字形、圆形和环形等。形和环形等。图图4.10 中值滤波的常用窗口中值滤波的常用窗口u中值滤波对于消除孤立点和线段的干扰十中值滤波对于消除孤立点和线段的干扰十分有用。分有用。u特别是对于二进噪声尤为有效，对于消除高斯特别是对于二进噪声尤为有效，对于消除高斯噪声的影响效果不佳。噪声的影响效果不佳。u对于一些细节较多的复杂图像，还可以多次使对于一些细节较多的复杂图像，还可以多次使用不同的中值滤波，然后通过适当的方式综合用不同的中值滤波，然后通过适当的方式综合所得的结果作为输出，这样可以获得更好的平所得的结果作为输出，这样可以获得更好的平滑和保护边缘的效果。滑和保护边缘的效果。

20、I1 = imread(blood1.tif);I=imnoise(I1,salt & pepper,0.02);imshow(I);K = medfilt2(I); %中值滤波中值滤波figure,imshow(K); 【例【例4.5】选用】选用33的窗口进行中值滤波。的窗口进行中值滤波。 （a） 原图原图 （b）结果图）结果图图图4.11 中值滤波中值滤波 1. K近旁均值滤波器（近旁均值滤波器（KNNF）l 在在mm的窗口中，属于同一集合类的像素，它们的的窗口中，属于同一集合类的像素，它们的灰度值将高度相关。灰度值将高度相关。l 被处理的像素（对应于窗口中心的像素）可以用窗口被处理的像素

21、（对应于窗口中心的像素）可以用窗口内与中心像素灰度最接近的内与中心像素灰度最接近的k个邻近像素的平均灰度个邻近像素的平均灰度来代替。来代替。（1）作一个）作一个mm的作用模板。的作用模板。（2）在其中选择）在其中选择K个与待处理像素的灰度差为最小的像个与待处理像素的灰度差为最小的像素。素。（3）用这）用这K个像素的灰度均值替换掉原来的值。个像素的灰度均值替换掉原来的值。4.3.3 4.3.3 边界保持类滤波边界保持类滤波u模板为模板为33，k3的的K近旁均值滤波器。近旁均值滤波器。图图4.12 K近旁均值滤波器近旁均值滤波器u在在K近旁均值滤波器中，不选近旁均值滤波器中，不选k个邻近像素的个邻

22、近像素的平均灰度来代替，而选平均灰度来代替，而选k个邻近像素的中值灰度个邻近像素的中值灰度来代替。来代替。图图4.13 K近旁中值滤波器近旁中值滤波器 2. K近旁中值滤波器近旁中值滤波器(KNNMF)l对图像上待处理的像素（对图像上待处理的像素（m，n）选它）选它的的55邻域。邻域。l在此邻域中采用图在此邻域中采用图4.14所示的模板。所示的模板。l计算各个模板的均值和方差，按方差计算各个模板的均值和方差，按方差排序，最小方差所对应的模板的灰度排序，最小方差所对应的模板的灰度均值就是像素（均值就是像素（m，n）的输出值。）的输出值。 3. 最小均方差滤波器最小均方差滤波器图图4.14 最小均

23、方差滤波器模板最小均方差滤波器模板l计算步骤如下：计算步骤如下：（1）按图做出）按图做出9个模板，计算出各自的方差。个模板，计算出各自的方差。（2）选出方差为最小的模板。）选出方差为最小的模板。（3）用该模板的灰度均值代替原像素的灰度值。）用该模板的灰度均值代替原像素的灰度值。iifNf1 其中 是指对应的模板，N是模板中像素的数量。l以方差作为各个邻域灰度均匀性的测度。l若邻域含有尖锐的边缘，灰度方差必定很大，而不含边缘或灰度均匀的邻域，方差就很小，那么最小方差所对应的邻域就是灰度最均匀邻域。l通过这样的平滑既可以消除噪声，又能够不破坏邻域边界的细节。iffN22)(1（4.9）均值和方差公

24、式：4.4 4.4 锐化锐化l 基本思想：基本思想：l有时还需要加强图像中景物的边缘和轮廓。有时还需要加强图像中景物的边缘和轮廓。l 边缘和轮廓常常位于图像中灰度突变的地方，边缘和轮廓常常位于图像中灰度突变的地方，因而可以直观的想到用灰度的差分对边缘和因而可以直观的想到用灰度的差分对边缘和轮廓进行提取。轮廓进行提取。l 梯度向量的幅度：梯度向量的幅度： yfxfGGfyx（4.10）)( fmagf2122yxGG 2122)()(yfxf（4.11）4.4.1 4.4.1 梯度锐化法梯度锐化法l二元函数二元函数f (x,y)在坐标点在坐标点(x,y)处的梯度定义为处的梯度定义为n 数字微分将

25、用差分代替：数字微分将用差分代替：|yxGGf|)| |,max(|yxGGf （4.12）（4.13）),(), 1(jifjifGx),() 1,(jifjifGy（4.14）（4.15）l为了降低运算量，常用绝对值或最大值运为了降低运算量，常用绝对值或最大值运算代替平方与平方根运算近似求梯度的幅度：算代替平方与平方根运算近似求梯度的幅度：图图4.15 沿沿x和和y方向的一阶差分方向的一阶差分 图图4.16 罗伯茨差分罗伯茨差分 l 所有梯度值都和相邻像素之间的灰度差分成所有梯度值都和相邻像素之间的灰度差分成比例。比例。l 利用它来增强图像中景物的边界。利用它来增强图像中景物的边界。l 采

26、用梯度进行图像增强的方法有：采用梯度进行图像增强的方法有： 第一种方法第一种方法: :使其输出图像的各点等于该点使其输出图像的各点等于该点处的梯度。即处的梯度。即 ),(),(jifjig（4.16） u缺点缺点:输出的图像在灰度变化比较小的区域，输出的图像在灰度变化比较小的区域，g(i,j)很小，显示的是一片黑色。很小，显示的是一片黑色。 第二种方法第二种方法: 对梯度值超过某阈值对梯度值超过某阈值T的像素选用梯的像素选用梯度值，而小于度值，而小于T时选用原像素点值。时选用原像素点值。 u 适当的选取适当的选取T，可以有效地增强边界而不影响比，可以有效地增强边界而不影响比较平滑的背景。较平滑

27、的背景。其他),(),(),(),(jifTjifjifjig（4.17）l 第三种方法：对梯度值超过第三种方法：对梯度值超过T的像素选用固定的像素选用固定灰度灰度LG 代替，而小于代替，而小于T时仍选用原像素点值时仍选用原像素点值l 这种方法可以使边界清晰，同时又不损害灰这种方法可以使边界清晰，同时又不损害灰度变化比较平缓区域的图像特性。度变化比较平缓区域的图像特性。其他),(),(,),(jifTjifLjigG（4.18） 第四种方法：将梯度值超过第四种方法：将梯度值超过T的像素选用梯度的像素选用梯度值，而小于值，而小于T时选用固定的灰度时选用固定的灰度LB 。即。即 l 这种方法将背景

28、用一个固定的灰度级这种方法将背景用一个固定的灰度级LB来表来表示，便于研究边缘灰度的变化。示，便于研究边缘灰度的变化。（4.19）其他,),(),(),(BLTjifjifjig 第五种方法：是将梯度值超过某阈值第五种方法：是将梯度值超过某阈值T T的像素选用固的像素选用固定灰度定灰度L LG G ，而小于该阈值，而小于该阈值T T时选用固定的灰度时选用固定的灰度L LB B 。 u该法生成的是二值图，根据阈值将图像分成边缘和该法生成的是二值图，根据阈值将图像分成边缘和背景，便于研究边缘所在的位置。背景，便于研究边缘所在的位置。其他,),(,),(BGLTjifLjig（4.20）【例【例4.

29、54.5】利用罗伯茨梯度对图像进行锐化处理。】利用罗伯茨梯度对图像进行锐化处理。 解解: : 程序如下程序如下: : I = imread(rice.tif);imshow(I);BW1 = edge(I,roberts,0.1); figure, imshow(BW1); （a）原图像）原图像 （b）结果图）结果图图图4.17 罗伯茨梯度的锐化罗伯茨梯度的锐化 除一阶微分外，还可以选用二阶微分算子。除一阶微分外，还可以选用二阶微分算子。l 一个连续的二元函数一个连续的二元函数f (x,y)，其拉普拉斯运算，其拉普拉斯运算定义为定义为 l 对于数字图像，拉普拉斯算子可以简化为对于数字图像，拉普

30、拉斯算子可以简化为22222yfxff（4.21）4.4.2 4.4.2 拉普拉斯算子（拉普拉斯算子（LaplacianLaplacian）(, ) 4 (, )( 1 , )( 1 , )(,1 )(, 1 )gi jf i jf ijf ijf i jf i j （4.22） 也可以表示为卷积的形式，即也可以表示为卷积的形式，即其中，其中，k1，l1， 取下式取下式kkrllsNjisrHsjrifjig1, 2 , 1 , 0,),(),(),(（4.23）),(srH0101410101Hl拉普拉斯的增强算子：拉普拉斯的增强算子：其对应的模板为：其对应的模板为：) 1,() 1,(),

31、 1(), 1(),(5),(jifjifjifjifjifjig（4.24）0101510101H【例【例4.64.6】应用拉普拉斯算子进行图像锐化处】应用拉普拉斯算子进行图像锐化处理。理。解解: : 程序如下程序如下 I = imread(2.gif); imshow(I); h=0 -1 0;-1 4 -1;0 -1 0; I2=imfilter(I,h); figure,imshow(I2);（a）原图像）原图像 （b）结果图）结果图图图4.18 拉普拉斯算子的锐化拉普拉斯算子的锐化 4.4.3 4.4.3 高通滤波高通滤波常用的高通模板有：常用的高通模板有：4.4.4 4.4.4 其

32、他锐化算子其他锐化算子1. Sobel算子算子1111811112H1212421213H0101510104H2122yxddS（4.25）)1, 1(), 1(2) 1, 1()1, 1(), 1(2) 1, 1(jifjifjifjifjifjifdx)1, 1() 1, (2) 1, 1()1, 1() 1, (2) 1, 1(jifjifjifjifjifjifdy用模板来表示用模板来表示: : 101202101xd121000121yd2. Prewitt算子算子用模板表示用模板表示 :2122yxpddSyxdd ,101101101xd111000111yd3. Isotropic算子算子用模板表示 ：2122yxIddSyxdd ,101202101xd121000121yd【例【例4.7】利用】利用Sobel算子和算子和Prewitt算子对图像进算子对图像进行锐化处理。行锐化处理。 I = imread(ric.tif);imshow(I);hs=fspecial(sobel);S=imfilter(I,hs);hp=fspecial(prewitt)P=imfilter(I,hp);figure,imshow(S,);figure,imshow(P,);（a）Sobel算子算子 （b）Prewitt算子算子图图4.19 锐化结果锐化结果71 结束语结束语