图像频域剖析.pptVIP

第四章 频域处理 1、本章所有滤波都是通过傅里叶变换在频域中实现。 2、这一章重点关注是如何在matlab中实现频域滤波。 3、这一章中我们将举例说明在图像增强中频域滤波，包括低通滤波、基本高通滤波与高频强调滤波。 4、1 二维离散傅里叶变换 1、二维离散傅里叶变换 2、二维离散傅里叶逆变换 3、傅里叶频谱（直观剖析一个变换主要方法是计算它频谱即F(u，v)幅度 4、1 二维离散傅里叶变换 4、1 二维离散傅里叶变换 4、2 在MATLAB中计算可视化二维DFT 1、快速傅里叶变换 F=fft2（f）该函数将返回一个大小仍为M*N傅里叶变换，数据排列形式仍为原点在左上角，而四个四分之一周期交汇于频率矩形中心。 2、使用傅里叶变换进行滤波时，需要对输入数据进行零填充。F=fft2（f，P，Q）. 3、傅里叶频谱可以使用函数 S=abs（F） 在这里我们知道通过显示频谱图像来进行可视化剖析是频域处理一个重要方面。 4、2 在MATLAB中计算可视化二维DFT 4、将变换原点移动到频率矩形中心。Fc=fftshift（F） 5、iffshift用于颠倒这种居中。 6、ifft2（F）用于计算傅里叶逆变换。 f=imread(Fig0403(a)(image).tif); imshow(f) F=fft2(f); S=abs(F); imshow(S,[]) Fc=fftshift(F); imshow(abs(Fc),[]) S2=log(1+abs(Fc)); imshow(S2,[]) 4、3 频域滤波 将简单介绍频域滤波概念以及如何使用MATLAB来实现频域滤波。 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) 双箭头两边表达式组成了傅里叶变换对。两个空间函数卷积可以通过计算两个傅里叶变换函数乘积逆变换得到。相反地，两个空间函数卷积傅里叶变换恰好等于两个函数傅里叶变换乘积。 4、3 频域滤波 基本上，频域滤波目是：选择一个滤波器传递函数，以便按照指定方式修改 F（u，v）。 基于卷积理论，我们知道为了在空间域中得到相应滤波后图像，仅需要计算积H(u，v)F（u，v）。 4.3.2 DFT基本步骤 1、使用函数paddedsize获得填充参数： PQ=paddedsize（size（f））； 2、得到使用填充傅里叶变换： F=fft2（f，PQ（1），PQ(2)）; 3、生成一个大小为PQ（1）*PQ(2)滤波函数H; 4、将变换乘以滤波函数： G=H.*F; 5、获得G傅里叶逆变换实部： g=real（ifft2（G））； 6、将左上角矩形修剪为原始大小: g=g(1:size(f,1):size(f,1)) 4、4 从空间滤波器获得频域滤波器 函数freqz2： H=freqz2（h，R, C） 其中h是一个二维空间滤波器，H是相应频域滤波器。使用方法见下面例子： f=imread(Fig0409(a)(bld).tif); imshow(f) F=fft2(f); S=fftshift(log(1+abs(F))); S=gscale(S); imshow(S) h=fspecial(sobel) h =1 2 1 0 0 0 -1 -2 -1 freqz2(h) H=freqz2(h,PQ(1),PQ(2)); H1=ifftshift(H); imshow(abs(H),[]) figure,imshow(abs(H1),[]) gs=imfilter(double(f),h); gf=dftfilt(f,H1);Imshow（gs,[]） imshow(gs,[]) figure,imshow(abs(gs),[]) figure,imshow(abs(gf),[]) 4、5 在频域中直接生成滤波器 4、5、1 函数dftuv提供了距离计算及其类型应用所需要网格数组。 4、5、2低通频域滤波器 1、理想低通滤波器ILPF 2、n阶巴特沃兹低通滤波器 3、高斯低通滤波器