《图像处理图像增强》PPT课件

《《图像处理图像增强》PPT课件》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

第4章图像增强4.1引言4.2直接灰度变换4.3直方图修正法4.4图像平滑4.5图像锐化4.6同态增晰4.7彩色增强 4.1引言根据所处理的空间不同：基于图像域的方法：直接在图像所在的空间进行处理基于变换域的方法：在图像变换域间接进行处理方法：空域频域图像增强处理策略：全局局部处理对象：灰度彩色空域处理：点处理、邻域处理、全图处理 4.2直接灰度变换4.2.1灰度线性变换4.2.2灰度非线性变换 灰度范围线性变换关系线性变换关系4.2.1灰度线性变换0abf(x,y)dcg(x,y)0abf(x,y)dcg(x,y) 4.2.1灰度线性变换1.全域线性变换设原始图像中所有像素灰度的最小值和最大值分别为f1（>0）和f2（<255）设结果图像中所有像素灰度的最小值和最大值分别为g1（≥0）和g2（≤255）线性变换表示公式： 4.2.1灰度线性变换1.全域线性变换若大部分像素的灰阶分布在[a,b]，小部分灰度级超出了此区域，为了改善增强效果，可以用如下所示的变换关系： 4.2.1灰度线性变换2．分段线性变换abcdMfMgf(x,y)g(x,y) 4.2.1灰度线性变换线性变换特例图像反转条件：原始图像（输入图像）：f(x,y)结果图像（输出图像）：g(x,y)灰度反转公式：f(x,y)=255-g(x,y) 线性变换特例255f(x,y)g(x,y)2552554.2.1灰度线性变换f(x,y)g(x,y)255 4.2.1灰度线性变换对比度扩展增强原图各部分的反差。也即增强原图里某两个灰度值间的动态范围来实现突出感兴趣的区间，相对抑制不感兴趣的灰度区域非线性变换往往以牺牲某些灰度范围的图像信息（灰度压缩），来换取其它灰度范围的图像信息的改善（灰度拉伸）。 4.2.2灰度非线性变换非线性变换：变换函数方程为例如：对数函数作为图像的映射函数 4.3直方图修正法4.3.1灰度直方图的定义4.3.2直方图的用途4.3.3直方图均衡化4.3.4直方图规定化 4.3.1灰度直方图的定义直方图直方图是图像的灰度——像素数统计图，即对于每个灰度值，统计在图像中具有该灰度值的像素个数，并绘制成图形，称为灰度直方图（简称直方图）。直方图模型表示图像中不同灰度级出现的相对频率Gray-levelhistogram 4.3.2直方图的用途P(k):具有该灰度级的像素的频数 4.3.2直方图的用途图像——直方图不可逆变换，多对一的变换直方图是多对一的映射结果，即多个图像可以生成相同的直方图，因此直方图作为一阶统计特征未反映相邻点之间的关系。但反映了图像的灰度散布范围等特征，在很多场合下，往往是重要特征。 六个像素0.24.3.2直方图的用途 4.3.2直方图的用途直方图性质1）2）如果一图像由两个不连接的区域组成，且每个区域的直方图已知，则整幅图像的直方图是该两个区域的直方图之和。 T双峰直方图4.3.2直方图的用途3）边界阈值的选择 4.3.2直方图的用途要点1．直方图表明在每一灰度级有多少个像素2．观察直方图可以看出不合适的数字化 动态范围宽了，对比度增强了4.3.2直方图的用途 4.3.3直方图均衡化直方图均衡把原始图的直方图变换为均匀分布的形式，增加像素灰度值的动态范围，提高图像对比度。Histogramequalization点处理增强：g=EH(f)表示。 4.3.3直方图均衡化g=EH(f)需满足：a）EH(f)在内单值递增，保证由黑到白b）保证动态范围一致原始图像的累计直方图满足上面两条件且能将f的分布转换为均匀分布 4.3.3直方图均衡化 DAnkPr(rk)f取成整数倍均衡后直方图07900.190.191/7(0.14)0.191/710230.250.443/7(0.428)0.252/78500.210.655/7(0.714)0.213/76560.160.816/7(0.857)0.16+0.08=0.244/73290.080.896/7(0.857)5/72450.060.957/7(1.00)0.06+0.03+0.02=0.116/71220.030.987/7(1.00)7/7810.021.007/7(1.00)8个灰级，总计64*64=4096点注意：离散均衡不可能拉平仅存5个灰级，宏观拉平，微观不可能平，层次减少，对比度提高。 原图f变换图0.11均衡后7/70.250.150.051/73/75/77/71/73/75/77/70.20.40.60.81.0DA1/73/75/70.190.250.210.240.250.05仅存5个灰级，层次减少，对比度提高。4.3.3直方图均衡化 4.3.3直方图均衡化直方图均衡化是一种非线性变换。直方图均衡的特点增加像素灰度值的动态范围，提高图像对比度。 4.3.4直方图规定化均衡化优点能自动增强整个图像的对比度，但具体的增强效果不易控制，处理的结果是全局均衡的直方图，实际中需特定形状的直方图，从而有选择的增强某个灰度值范围内的对比度。分别对原始直方图和规定化处理后的直方图进行均衡化处理 0.30.20.101/73/75/71rkPr(rk)0.30.20.101/73/75/71zkPz(zk)0.30.20.101/73/75/71zkPz(zk)1.00.80.60.40.201/73/75/71rks－T(rk) 0.115,6,7->771.000.151.000.02817步骤和结果0.240.210.250.19000变换后直方图93,4->62->51->40->3确定映射关系87766543映射|V2-V1|最小70.850.650.350.1500060.200.300.200.15000规定直方图P(z)规定累积直方图V250.980.950.890.810.650.440.19原始累积直方图V140.030.060.080.160.210.250.19原始直方图P(r)31222453296568501023790原始直方图各灰度级像素26543210原始图像灰度级1运算序号 a)原图b)规定化函数c)直方图规定化后的结果d)图c的直方图 4.4.1邻域平均法4.4.2中值滤波4.4.3多图像平均法4.4.4频域低通滤波法4.4图像平滑 4.4图像平滑空域滤波是在图像空间借助模板进行邻域操作完成线性、非线性运算功能1）平滑：低通滤波器。目的：在提取较大目标前去除太小的细节或将目标内的小间断连接起来消除噪声2）锐化：高通滤波器，增强被模糊的细节 滤波处理方法空域：取局部邻域（2M+1）×（2M+1）邻域的加权和局域处理K4K3K2K5K0K1K6K7K8R4.4图像平滑 4.4图像平滑目的：减少噪声1）加性噪声2）乘性噪声3）量化噪声4）盐和胡椒噪声噪声：独立同分布的高斯白噪声，均值为0，方差σ 4.4.1邻域平均法010101010111101111mask消除麻点噪声1/41/8 为减轻经平滑后产生的模糊效应，另一种阈值平均法例：4.4.1邻域平均法 4.4.1邻域平均法均值滤波实现图像平滑最常见的方法是在像素邻域内求局部均值，称为均值滤波。 4.4.1邻域平均法3x3均值滤波O(X,Y)=(I(X-1,Y-1)+I(X,Y-1)+I(X+1,Y-1)+I(X-1,Y)+I(X,Y)+I(X+1,Y)+I(X-1,Y+1)+I(X,Y+1)+I(X+1,Y+1))/9滤波核h(x，y)：111111111mask111121111121242121 4.4.1邻域平均法步骤:1）模板游走2）将mask下对应的灰度值相加，求平均值3）用均值代替f(x,y)h(x,y)矩阵的元素之和乘前面系数为1，h(x,y)矩阵中心的元素占的比例越小，越平滑，图像越模糊4）对图像的四周边缘：补0或者不处理边缘 [例1]设16x16点阵的假想图像如右图所示。00000000000000000000000000000000000000000000000000000000000000000000999999990000000099999999000000009999999900000000999999990000000099999999000000009999999900000000999999990000000099999999000000000000000000000000000000000000000000000000000000000000000000004.4.1邻域平均法 采用3x3均值滤波。滤波核为：结果如右图所示1111/9111111?????????????????00000000000000??00000000000000??00123333332100??00246666664200??00369999996300??00369999996300??00369999996300??00369999996300??00369999996300??00369999996300??00246666664200??00123333332100??00000000000000??00000000000000?????????????????4.4.1邻域平均法 [例2]设16x16点阵的假想图像如右图所示00000000000000000000000008000000000000000000000000000000000000000000999999990000000099999999000000009999999900000000999999990000000099991999000000009999999900000000999999990000000099999999000000000000000000000000000000000000008000000000000000000000000000004.4.1邻域平均法 采用3x3均值滤波滤波核为：结果如右图所示?????????????????00000001110000??00000001110000??00123333332100??00246666664200??00369999996300??00369999996300??00369988896300??00369988896300??00369988896300??00369999996300??00246666664200??00123333332100??11100000000000??11100000000000?????????????????4.4.1邻域平均法1111/9111111 5x5均值滤波滤波核：1/251/251/251/251/25111111/251/251/251/251/25111111/251/251/251/251/25或1/25111111/251/251/251/251/25111111/251/251/251/251/2511111模板越大，模糊作用越强4.4.1邻域平均法 均值滤波的特点局部求均值的运算或平均计算使数字信号变“平坦”，可以在图像中消除或抑制噪声。同时，图像中景物边缘也会不同程度地变得模糊。4.4.1邻域平均法 4.4.2中值滤波中值滤波中值滤波与均值滤波的区别仅限于：中值滤波是求局部中值而不是局部均值，即对参与计算的像素灰度值按大小排序，然后取位置居中的像素灰度值。Medianfiltering 目的：既要消除噪声又要保持图像的细节步骤:1）模板游走2）将mask下对应的灰度值（奇数）排序3）用中间值代替f(x,y),消除孤立的噪声点mask大小不一样，效果不一样，与叠加的噪声有关系窗口形状方形、十字形、圆形、圆环形4.4.2中值滤波 使用二维滤波需注意1）有尖顶角几何结构的图像，一般采用十字窗，大小不超过最小有效物体的尺寸2）有较多的点、线、尖顶角的细节结构，不用中值滤波用3×3方形，图像全为0用3×3十字，保留了线状细节，丢失了点状细节4.4.2中值滤波 [例3]设16x16点阵的图像如右图所示00000000000000000000000008000000000000000000000000000000000000000000999999990000000099999999000000009999999900000000999999990000000099991999000000009999999900000000999999990000000099999999000000000000000000000000000000000000008000000000000000000000000000004.4.2中值滤波 采用3x3中值滤波，结果如右图所示?????????????????00000000000000??00000000000000??00000000000000??00009999990000??00099999999000??00099999999000??00099999999000??00099999999000??00099999999000??00099999999000??00009999990000??00000000000000??00000000000000??00000000000000?????????????????4.4.2中值滤波 4.4.2中值滤波中值滤波的机理在图像局部中，杂散噪声点的灰度值较大或较小，求中值可以自动将其消除。但中值滤波不同于均值滤波，中值直接取自某个像素的灰度值，所以能较好地保持图像景物原状。中值滤波的特点可以消除杂散噪声点而不会或较小程度地造成边缘模糊。 4.4.3多图像平均法多图像平均法是利用对同一景物的多幅图像取平均来消除噪声产生的高频成分，在图像采集中常应用这种方法去除噪声。 4.4.4频域低通滤波法(1)理想低通滤波器（ILPF）H(u,v)=1当D(u,v)≤D00当D(u,v)>D0其中：D(u,v)＝(u2+v2)1/2是点(u,v)到频率平面原点的距离副作用：图象模糊，出现振铃效果H(u,v)D(u,v)D0 (2)巴特沃斯低通滤波器（BLPF）传递系数：4.4.4频域低通滤波法01H(u,v)D(u,v) （3）指数低通滤波器（ELPF）4.4.4频域低通滤波法01H(u,v)D(u,v) （4）梯形低通滤波器（TLPF）类别振铃程度图像模糊程度噪声平滑效果ILPFTLPFELPFBLPF严重较轻无无严重轻较轻很轻最好好一般一般4.4.4频域低通滤波法H(u,v)D0D101D(u,v) 4.5图像锐化4.5.1微分法4.5.2高通滤波法 边缘锐化（Sharpening）补偿图像的轮廓，突出图像中景物的边缘或纹理，使图像清晰——空域高通滤波（俗称勾边处理）。4.5图像锐化 4.5图像锐化为什么需要锐化图像传输变换（未聚焦好）、受到各种干扰而退化——图像模糊，而图像的判读和识别中，常需突出目标的轮廓或边缘信息。数学原理图像模糊的原因——图像被平均或积分，为实现图象的锐化，需反运算“微分”——增强高频分量，使图像边缘清晰，但同时也增强了噪声条件原图像有较高的SNR 4.5图像锐化边缘锐化的原理图像中景物的边缘或纹理是灰度值发生突变的地方，数字信号处理采用差分或梯度计算来检测边缘。若在原图上（或一定的灰度值上）叠加该检测结果有勾边效果。边缘锐化的缺点边缘突出的同时图像中噪声也会被突出。 图像模糊的两种理解1．平均或积分效应引起——用微分处理；2．模糊可看作图像上高频分量被削弱——高频增强，引入微分：4.5.1微分法注意噪声亦属高频分量，往往会随高频增强而突出。因此往往需要去噪后再锐化——微分 讨论微分算子一般情况：希望求微分算子是各向同性的，即其微分效果不随特征方向不同而不同。可证明：偏导数的平方和是各向同性的，梯度运算、拉普拉斯运算都符合上述条件。4.5.1微分法 x´y´yf(x,y)θ4.5.1微分法 1．梯度法4.5.1微分法 (i,j)(i+1,j)(i,j+1)(i,j)(i+1,j)(i,j+1)(i+1,j+1)4.5.1微分法1．梯度法 4.5.1微分法1．梯度法改进的梯度算法：Δ—阈值 2．拉普拉斯运算法4.5.1微分法 -1-14-1-1-1-1-18-1-1-1-1-1注意:1)图像中灰度变化小时，g(x,y)输出小或为负2）图像中灰度变化大时，g(x,y)的增强可能超出范围，需要变回0～L-1尺度中4.5.1微分法2．拉普拉斯运算法 2222222mask0-10-14-10-10例:灰度级为422241111mask0-10-14-10-104.5.1微分法2．拉普拉斯运算法 3．Roberts算子(i,j)(i+1,j)(i,j+1)(i+1,j+1)4.5.1微分法 4．高频增强滤波器：光学操作将聚焦的正像与散焦的负像在底版上叠加，散焦的负像相当于一个模糊掩模。-1-1-19-1-1-1-1-11-2-25-211-21-21161-2-21-24.5.1微分法经平滑后的图像（模糊）原图 5．定向滤波：沿特定方向增强，有塑像效果-1-c-11c1-1-cc1-11c11-1-1-c-111c-c-14.5.1微分法 频域：G(u,,v)=H(u,v)·F(u,v)其中：F(u,v)：原始图象傅立叶频谱G(u,v)：平滑后图象的傅立叶频谱H(u,v)：滤波转移函数H(u,v)函数的定义方法很多，针对具体情况选用不同方法。FFTH(u,v)IFFTf(x,y)F(u,v)G(u,v)g(x,y)4.5.2高通滤波法 图像中的均匀与不均匀反映了频率高低不同抑制低频（增强高频）—锐化抑制高频（增强低频）—平滑滤波传递函数g(x,y)=h(x,y)*f(x,y)G(u,v)=H(u,v)·F(u,v)以下讨论考虑对F(u,v)的实部、虚部影响完全相同的滤波转移函数——零相移滤波器4.5.2高通滤波法 (1)理想滤波器D04.5.2高通滤波法 (2)巴特沃斯滤波器4.5.2高通滤波法 (3)指数滤波器4.5.2高通滤波法 (4)梯形滤波器D1D0H(u,v)D(u,v)4.5.2高通滤波法 4.6同态增晰作用消除图像上照明不均的问题，增加暗区的图像细节，同时又不损失亮区的图像细节，它在频域中同时将图像亮度范围进行压缩和将图像对比度进行增强成像物理背景人眼对图象亮度响应具有类似于对数运算的非线性形式f(x,y)=I(x,,y)·R(x,,y)I(x,,y)：照射分量（低频）R(x,y)：反射分量（高频）(图象细节的不同在空间作快速变化) 分析关心反射信息，但室内外照射分量强度不同，图片明暗不均，能否消除照度不均，而增强反射部分比重?过程f(x,y)=i(x,y)r(x,y)i(x,y)——照射分量,低频区，r(x,y)——反射分量，反映图像的细节分量，处于高频区4.6同态增晰f(x,y)lnFFTH(u,v)高频增强FFT-1expg(x,y) 步骤：（1）z(x,y)=ln[f(x,y)]=ln[I(x,y)]+ln[R(x,y)]把频谱分开（2）Z(u,v)=I(u,v)+R(u,v)傅立叶变换（3）S(u,v)=H(u,v)·Z(u,v)同态滤波函数用H(u,v)处理Z(u,v)（4）s(x,y)=F-1[S(u,v)]=i(x,y)+r(x,y)（5）g(x,y)=exp[s(x,y)]=exp[i(x,y)+r(x,y)]=i0(x,y)·r0(x,y)4.6同态增晰 如果图像照明不均，图像上各部分的亮度会有起伏，压缩照度分量的灰度范围或频域上消弱照度分量的频谱分量。因反射分量反映细节，利用对比度增强这一分量的对比度或频域上加大反射频谱成分，使暗区细节增强，并保留亮区图像细节4.6同态增晰2.01.51.00.5H(u,v)D(u,v) a)同态滤波处理前b)同态滤波处理后（压缩图像的动态范围，增加了图像各部分之间的对比度）4.6同态增晰 4.7彩色增强4.7.1伪彩色增强4.7.2真彩色增强 4.7.1伪彩色增强(1)亮度切割x,yl1L红色蓝色亮度切割剖面示意 伪彩色效果图4.7.1伪彩色增强 （2）从灰度到彩色的变换红色转换器绿色转换器蓝色转换器IR(x,y)IG(x,y)IB(x,y)输入三个电子枪，可得到其颜色内容由3个变换函数调制的混合图像4.7.1伪彩色增强 4.7.1伪彩色增强a)红色变换特性b)绿色变换特性IRL0L/2Lf(x,y)L/2Lf(x,y)IGL0L/2Lf(x,y)IBL0L/2Lf(x,y)ILc)蓝色变换特性d)合成特性 4.7.2真彩色增强从彩色到彩色的一种转换目的1）引起人们的特别关注2）人眼对不同颜色的灵敏度不同 真彩色效果例图4.7.2真彩色增强 作业4-1试给出把灰度范围（0，10）拉伸为（0，15），把灰度范围（10，20）移到（15，25），并把灰度范围（20，30）压缩为（25，30）的变换方程。4-2试给出变换方程T(z),使其满足在10≤z≤100的范围内，T(z)是lgz的线性函数 0000.190.250.210.240.110007901023850985448012345670.140.220.260.170.090.060.040.0256092010467053562671707201234567nk/nnk/ng(x,y)nk/nnkf(x,y)表a表b作业4-5已知一幅64×64的数字图像，其灰度级有8个，各灰度级出现的频数如表a所示。试将此幅图像进行直方图变换，使其变换后的图像具有如表b所示的灰度级分布，并画出变换前后图像的直方图 4-7试设计一个程序实现n×n的中值滤波器。当模板中心移过图像中每个位置时，设计一种简便的中值更新方法。作业