Slide1第4章图像增强Slide2内容提要图像增强是一类对图像降质进行改善的方法,它们有选择性地突出图像的边缘、轮廓、对比度等特征,以便于显示、观察或进一步分析与处理。图像增强的目的在于采用一系列技术改善图像的视觉效果,提高图像的可懂度;或者将图像转化成一种更适合于人或机器进行分析处理的形式。Slide3内容提要单点增强:灰度级校正、灰度变换、灰度直方图变换。区域增强的平滑方法:邻域平均法、中值滤波和各种边界保持类滤波方法。区域增强的锐化方法:梯度锐化法、拉普拉斯算子、高通滤波及其他常用的锐化算子。Slide44.1概述4.1.1图像增强的目的首要目标:处理图像,使其比原始图像更适合于特定应用。增强的方法是因应用不同而不同的。图像增强方法只能有选择地使用。增强的结果一般靠人的主观感觉加以评价。Slide54.1.2图像增强技术的分类图像增强技术大致分为空间域增强和频率域增强两类。1.空间域增强法在空间域直接对像素灰度值进行运算。f(x,y)是待增强的原始图像,g(x,y)是已增强的图像,h(x,y)是空间运算函数。Slide6空间域增强模型对点操作(如灰度变换、直方图变换等)有g(x,y)=f(x,y)·h(x,y)(4.1)对于区域操作(如平滑、锐化等)有g(x,y)=f(x,y)*h(x,y)(4.2)图4.1空间域增强模型Slide72.频率域增强法在频率域利用二维滤波器H(u,v)对f(x,y)进行滤波,得到新的频谱G(u,v),即G(u,v)=F(u,v)·H(u,v)(4.3)图4.2频率域增强模型Slide8H(u,v)的性质可能是低通,起平滑作用;也可能是高通,起锐化作用。实际的图像增强方案可能综合上述两种技术。如同态滤波增强包含了空间域灰度的非线性运算,也有高频增强环节。Slide94.2灰度修正图像的退化会引起灰度级的变化。通过简单和有效的点运算,可以改善图像的显示效果,达到灰度修正的目的。点运算:一幅输出图像上每个像素的灰度值仅由相应输入像素的灰度值决定,而与像素点所在的位置无关,与相邻的像素之间也没有运算关系。Slide10点运算:指原始图像的像素灰度值通过运算后产生新图像的对应的灰度值。像素值通过运算改变之后,可以改善图像的显示效果。这是一种像素的逐点运算。是旧图像与新图像之间的映射关系。典型的点运算:对比度增强、对比度拉伸或灰度变换。Slide11灰度修正的方法一般有三种方法:(1)灰度级校正解决成像不均匀问题。(2)对比度增强解决图像曝光不足问题。(3)直方图修正以突出所需要的图像特征。Slide124.2.1灰度级校正在成像过程中,如光照的强弱、感光部件的灵敏度、光学系统的不均匀性、元器件特性的不稳定等均可引起图像亮度分布的不均匀。灰度级校正在图像采集系统中对图像像素进行逐点修正,使得整幅图像能够均匀成像。Slide13使理想图像发生畸变的比例运算设理想真实的图像为f(x,y),实际获得的含噪声的图像为g(x,y),则有g(x,y)=e(x,y)f(x,y)(4.1)e(x,y)是使理想图像发生畸变的比例因子。知道了e(x,y),就可以求出不失真图像。采用一幅灰度级为常数C的图像成像,若经成像系统的实际输出为gC(x,y),则有gC(x,y)=Cf(x,y)(4.2)Slide14标定系统失真系数的方法可得比例因子:可得实际图像g(i,j)经校正后所恢复的原始图像。注意:乘了一个系数C/gc(i,j),校正后可能出现“溢出”现象灰度级值可能超过某些记录器件或显示设备输入信号的动态可范围需再作适当的灰度变换,最后对变换后的图像进行量化。1,,ceijgijCSlide15图4.3非均匀光照的校正(a)校正前的图像(b)校正后的图像Slide164.2.2灰度变换为了将图像灰度级的整个范围或一段范围扩展或压缩到记录或显示设备的动态范围内,可使图像动态范围增大,图像对比度扩展使图像变得清晰/图像上的特征变得明显。环境光源太暗,使灰度值偏小,就会使图像太暗看不清。如果环境光源太亮,又使图像泛白。通过灰度变换,就可以将灰度值调整到合适的程度。灰度变换可分为线性变换、分段线性变换和非线性变换几种方法。Slide171.线性变换灰度g与灰度f之间的关系为(1)变换使得图像灰度范围增大,即对比度增大,图像会变得清晰;(2)变换使得图像灰度范围缩小,即对比度减小。[]bagafaba图4.4线性变换Slide18【例4.1】采用线性变换进行图像增强。应用函数imadjust将图像在0.3×255~0.7×255灰度之间的值通过线性变换映射到0~255之间。【解】实现的程序如下:I=imread('pout.tif');imshow(I);figure,imhist(I);%显示原始图像的直方图J=imadjust(I,[0.30.7],[]);%使用imadjust函数进行灰度的线性变换figure,imshow(J);figure,imhist(J)%显示变换后图像的直方图Slide19图4.5图像线性变换Slide202.分段线性变换对整个灰度区间进行分段,采用分段线性函数进行变换。这种变换突出了感兴趣的目标或灰度区间,相对抑制那些不感兴趣的灰度区间。常用的是三段线性变换。图4.6三段线性变换Slide21对灰度区间[a,b]进行了线性拉伸,而灰度区间[0,a]和[b,fmax]则被压缩。调整折线拐点的位置及控制分段直线的斜率,对图像的任一灰度区间进行拉伸或压缩。在遥感图像分类中,感兴趣的地貌特征可能有明显的灰度变化,而那些过黑或过白的像素往往对应于玄武岩、水、冰等。Slide22图4.7三段线性变换实例(a)原始图像(b)增强效果Slide233.非线性灰度变换当用某些非线性函数如对数、指数函数等作为映射函数时,可实现灰度的非线性变换。对数变换的一般表达式为:g=a+clg(f+1)对数变换可以增强低灰度级的像素,压制高灰度级的像素,使灰度分布与视觉特性相匹配。Slide244.8对数变换后的Couple图像Slide254.2.3灰度直方图变换1.直方图的概念对于连续图像,其灰度分布的统计特性用概率密度函数(PDF)刻画。离散图像直方图:指图像中各种不同灰度级像素出现的相对频率。在数字图像处理中,灰度直方图是简单且实用的工具,对图像的采集、处理和分析都可以有效地利用直方图。Slide26归一化的直方图直方图的横坐标是图像的各灰度级,纵坐标是各个灰度出现的像素个数。归一化的直方图(histogram)定义为灰度级出现的相对频率。即(4.13)式中,N表示像素的总数;nk表示灰度级为k的像素的数目。()/rkPknNSlide27【例4.2】求图像的归一化直方图。【解】lena图像是彩色图像,进行格式转换。I=imread('lena.jpg');J=rgb2gray(I);%将彩色图像转换为灰度图像imshow(J);N=numel(J);%求图像像素的总数Pr=imhist(J)/N;%显示原始图像的直方图k=0:255;figure,stem(k,Pr)Slide28图4.9直方图(a)lena图像(b)lena图像的直方图灰度直方图描述了图像的概貌。Slide292.直方图的修正直方图变换后可使图像的灰度间距拉开或使灰度分布均匀,从而增大对比度,使图像细节清晰,达到增强的目的。直方图变换有两类直方图均衡化;直方图规定化。Slide30直方图均衡化通过对原图像进行某种变换,使得图像的直方图变为均匀分布的直方图。灰度级连续的灰度图像:当变换函数是原图像直方图累积分布函数时,能达到直方图均衡化的目的。对于离散的图像,用频率来代替概率。Slide31【例4.3】对图像进行直方图均衡化。假定有一幅总像素为n=64×64的图像,灰度级数为8,各灰度级分布列于表4.1中。(1)按式(4.14)求变换函数Sk’(2)计算Sk’’(3)Sk的确定(4)计算对应每个sk的nsk(5)计算ps(sk)Slide32表4.1一幅图像的灰度级分布k01234567rk01/72/73/74/75/76/71nk790102385065632924512281pr(rk)0.190.250.210.160.080.060.030.02Sk’0.190.440.650.810.890.950.981Sk’’1/73/75/76/76/7111Sk1/73/75/76/71nsk7901023850985448pr(sk)0.190.250.210.240.11Slide33【例4.4】直方图均衡对图像进行增强。在MATLAB环境中,待处理图像为tire.tif。I=imread('tire.tif');J=histeq(I);%完成直方图均衡化imshow(I);%显示直方图均衡化前的图像figure,imhist(I);%均衡化前的直方图figure,imshow(J);%显示直方图均衡化后的图像figure,imhist(J);%均衡化后的直方图Slide34图4.10直方图均衡Slide354.3同态增晰4.3.1问题的由来物体受到不均匀的照度一类图像的灰度级动态范围很大,黑与白形成强烈反差,而感兴趣目标的灰度级范围却很小,分不清目标的灰度层次和细节,图像上对应照度暗的部分,其细节较难辨别。目的消除不均匀照度的影响而又不损失图像细节。同态滤波可以消除不均匀照度的影响,增强图像细节同时也是图像复原的一种方法。Slide364.3.2增晰原理图像同态增晰系统采用合适的滤波特性函数,可以既使图像灰度动态范围压缩,又能让感兴趣的物体图像灰度级扩展,从而使图像清晰。同态系统是服从广义叠加原理的各类非线性系统。对同态系统信号处理,特别适合处理乘法组合信号和卷积组合信号。(4.21)(,)(,)(,)fxyixyrxySlide37图4.11同态增晰原理按照高通滤波器设计,压缩低频分量,提升高频分量。照明函数虽然频率变化缓慢,但幅度变化大,数字化占用很多比特数,所以要压缩;反射函数描述的是人们感兴趣的景物区频率变化较快,但灰度变化很小,层次不清,细节不明,应该扩展。Slide384.3平滑(Smoothing)区域增强算法包括平滑算法和锐化算法。从频率域看,低通滤波可以对图像进行平滑去噪处理高通滤波可对图像进行边缘锐化处理。本节将介绍邻域平均法中值滤波边界保持类滤波等。Slide394.4.1图像噪声数字图像往往要经过采集、处理、存储、传输等一系列加工变换,而由电气系统和外界引入的图像噪声也将在这些过程中随之引入,可能严重影响图像的质量。这些过程将使得图像噪声的精确分析变得十分复杂。图像噪声消除或减低在图像预处理中的地位显得十分重要Slide401.图像噪声的分类(1)按其产生的原因,可以分为外部噪声和内部噪声。(2)按统计特性是否随时间变化,可以分为平稳噪声和非平稳噪声。(3)按噪声幅度随时间分布形状来定义,有高斯噪声、瑞利噪声、泊松噪声等。(4)按噪声频谱形状来分类,有白噪声、1/f噪声、三角噪声等。(5)按噪声和信号之间的关系,可分为加性噪声和乘性噪声。Slide412.MATLAB为图像加噪声的函数函数形式为J=imnoise(I,type,parameters)其中,I为原图像的灰度矩阵,J为加噪声后的灰度矩阵。type为噪声种类,parameters是允许修改的参数,可以默认。type可以有五种。Slide424.4.2邻域平均法1.空间域分析大部分的噪声都可以看作是随机信号,对图像的影响可以看作是孤立的。某一像素,如果