6第六章频率域图像增强资料.ppt

在谱中叠加的圆周分别有5，15，30，80，230像素的半径。 这些圆周包围的图像功率的百分比分别为92.0%，94.6%，96.4%，98%，99.5%。 理想低通滤波器 5 15 30 80 230 频域： 空间域： §6.4 .3高斯低通滤波器（GLPF） 有更加平滑的过渡带，平滑后的图像没有振铃现象 与BLPF相比，衰减更快，经过GLPF滤波的图像比BLPF处理的图象更模糊一些 高斯低通滤波器 5 15 30 80 230 §6.4 .4 图像锐化——高通滤波器 高通滤波器—空域通过用带有正有负的系数模版实现高通滤波； 频率域滤波器越宽，滤除的频率成份越多，在空间域意味着滤波器越窄，模板越小，检测边缘越少； 频域越窄，空域越宽，模板越大，检测边缘越多。 高通滤波与低通滤波的作用相反，它使高频分量顺利通过，而使低频分量受到削弱。 与低通滤波器相对应，频率域内常用的高通滤波器有2种： 理想高通滤波器 高斯高通滤波器 理想高通滤波器 高斯高通滤波器 滤波器的空间表示： §6.4 .5理想高通滤波器(IHPF) D0=15 D0=80 D0=30 §6.4 .5理想高通滤波器(IHPF) 截止频率距原点的距离为D0的GHPF的传递函数： §6.4 .6高斯高通滤波器(IHPF) D0=15 D0=80 D0=30 §6.4 .6高斯高通滤波器(IHPF) * * * 。数学表达式： 设周期信号为f(t)，其重复周期为T1，基波角频率 为 , 当满足狄利赫利条件时：可有： (3.1) (3.1)式表明：满足以上条件的周期信号可以分解成直流分量 和许多正、余弦分量；这些分量的频率必定是基频的整倍数；叫谐波 分量：二次谐波、三次谐波.....； 系数 a0、an、bn 代表各次分量的幅度，它们 的大小决定了周期信号的波形。  因为三角函数集是正交函数集，所以  * 各正交函数分量若 且 或 (3.5)这里,把(3.5)代入(3.2)、(3.3)、(3.4) 得到直流分量 (3.6)余弦分量的幅度 (3.7)正弦分量的幅度 (3.8)其中 n=1,2,3,.....通常取 积分区间为（0,T1 )或 (- T1 /2,+T1 /2)若同频率项合并，则另一种表示为 (3.9) (3.10)比较两种表示，得 (3.11) (3.12) (3.13) (3.14) (3.15) (3.16) 以上n=1,2,3,... * * 优点：概念清楚，通阻分明； 1 缺点：产生模糊和振铃现象， D（u,v） D0越小，模糊越厉害； 特性曲线 * 优点：具有比较平滑的过渡带，较快的衰减特性，无振铃现象； 缺点：实用效果较巴特沃思稍差些，稍微模糊一点； FT 移中的变换： 能量分布于四角 能量集中于中心 移中FT 原图像f（x，y） 移中FT Matlab显示频谱 I=imread(lena2.png); I=im2double(I); fcoef=fft2(I); spectrum=fftshift(fcoef); temp=log(1+abs(spectrum)); subplot(1,2,1); imshow(temp,[]); subplot(1,2,2); imshow(I); §6.2.4直接计算DFT的问题及改进途径 设x(n)为N点有限长序列，其DFT为 一般来说，x(n)和 都是复数，X(k)也是复数， 因此每计算一个X(k)值，需要N次复数乘法以及 （N-1)次复数加法。而X(k)一共有N个点，所以 完成整个DFT运算总共需要 次复数乘法及 N(N-1)次复数加法。 时 共需16次乘法，12次加法。 1）直接计算DFT的问题及改进途径 由于一次复数乘法需用四次实数加法；一次复数 加法则需二次实数加法。因此每运算一个X(k)需 要4N次实数乘法及2N+2(N-1)=2(2N-1)次实数加 法。所以整个DFT运算总共需要 次实数乘法和 次加法。