数字图像处理04课件

数字图像处理第四章频域处理4/2/20231本章主要内容4/2/20232漳州师范学院计算机科学与工程系4.1二维离散傅里叶变换令f(x,y)表示一幅大小为M×N的图像,其中x=0,1,2,…,M-1和y=0,1,2,…,N-1.的二维傅里叶变换可表示为傅里叶逆变换可表示为式中,F(u,v)的值称为傅里叶系数。空间域是由f(x,y)所张成的坐标系.频域系统是由F(u,v)所张成的坐标系.4/2/20233漳州师范学院计算机科学与工程系4.2在MATLAB中计算并可视化二维DFT快速傅里叶变换(函数fft2)格式1：F=fft2(f)格式2：F=fft2(f,P,Q)快速傅里叶逆变换(函数ifft2)格式：f=ifft2(F)将变换的原点移动到频率矩形中心(函数fftshift)格式：Fc=fftshift(F)将频率矩形中心移动到变换的原点(函数ifftshift)格式：F=ifftshift(Fc)上述函数使用演示：图4.3（f43.m）P,Q的含义？4/2/20235漳州师范学院计算机科学与工程系4.3频域滤波4.3.1基本概念空间域和频域线性滤波的基础都是卷积定理。f(x,y)*h(x,y)H(u,v)F(u,v)f(x,y)h(x,y)H(u,v)*F(u,v)滤波传递函数H(u,v)低通滤波器—衰减高频分量而保持低频分量不变折叠误差的干扰—对周期函数执行卷积导致的相邻周期之间的干扰.可通过使用0来填充函数的方法来避免.图像的填充及用于图像填充的函数paddedsize例4.1使用填充和不使用填充的滤波效果(e4_1.m)4/2/20236漳州师范学院计算机科学与工程系4.3频域滤波4/2/20237漳州师范学院计算机科学与工程系4.3频域滤波4.3.2DFT滤波的基本步骤f(x,y)输入图像g(x,y)滤波后的图像傅里叶变换滤波函数H(u,v)逆傅里叶变换F(u,v)H(u,v)F(u,v)4/2/20239漳州师范学院计算机科学与工程系4.3频域滤波4.3.3用于频域滤波的M函数Functiong=dftfilt(f,H)F=fft2(f,size(H,1),size(H,2))G=H.*Fg=real(ifft(G))g=g(1:size(f,1),1:size(f,2))4/2/202310漳州师范学院计算机科学与工程系4.4从空间滤波器获得频域滤波器IPT提供了生成空间滤波器的函数（3.5）fspecial对给定的空间滤波器h，对应的频域滤波器H为：H=fft2(h,PQ(1),PQ(2))函数freqz2:计算FIR滤波器的频率响应H=freqz2(h,R,C)演示4/2/202311漳州师范学院计算机科学与工程系4.5在频域中直接生成滤波器4.5.2低通频域滤波器理想低通滤波器（ILPF）n阶巴特沃兹低通滤波器（BLPF）高斯低通滤波器（GLPF）4/2/202313漳州师范学院计算机科学与工程系4.5在频域中直接生成滤波器4.5.3线框图与表面图函数mesh函数gridon/off函数axison/off函数view4/2/202314漳州师范学院计算机科学与工程系4.6锐化频域滤波器4.6.1基本的高通滤波