(完整版)(精心整理)图像的傅里叶变换.pptVIP

* . * . * . * . * . * . * . * . * . 主极大的值用Fmax表示，第一个旁瓣的峰值用Fmin表示 * . 例题：对一幅图像实施二维DFT，显示并观察其频谱。 解：源程序及运行结果如下： %对单缝进行快速傅里叶变换，以三种方式显示频谱， %即：直接显示（坐标原点在左上角）；把坐标原点平 %移至中心后显示；以对数方式显示。 f=zeros(512,512); f(246:266,230:276)=1; subplot(221),imshow(f,[]),title(单狭缝图像) F=fft2(f); %对图像进行快速傅里叶变换 S=abs(F); subplot(222) imshow(S,[]) %显示幅度谱 title(幅度谱（频谱坐标原点在坐上角）) * . Fc=fftshift(F); %把频谱坐标原点由左上角移至屏幕中央 subplot(223) Fd=abs(Fc); imshow(Fd,[]) ratio=max(Fd(:))/min(Fd(:)) %ratio = 2.3306e+007,动态范围太大，显示器无法正常显示 title(幅度谱（频谱坐标原点在屏幕中央）) S2=log(1+abs(Fc)); subplot(224) imshow(S2,[]) title(以对数方式显示频谱) 运行上面程序后，结果如下： * . * . 二维离散傅里叶变换的性质 线性性 证明： * . %imagelinear.m %该程序验证了二维DFT的线性性质 f=imread(D:\chenpc\data\thry\chpt4\Fig4.04(a).jpg); g=imread(D:\chenpc\data\thry\chpt4\Fig4.30(a).jpg); [m,n]=size(g); f(m,n)=0; f=im2double(f); g=im2double(g); subplot(221) imshow(f,[]) title(f) subplot(222) imshow(g,[]) title(g) * . F=fftshift(fft2(f)); G=fftshift(fft2(g)); subplot(223) imshow(log(abs(F+G)),[]) FG=fftshift(fft2(f+g)); title(DFT(f)+DFT(g)) subplot(224) imshow(log(abs(FG)),[]) title(DFT(f+g)) * . * . 可分离性 二维DFT可视为由沿x,y方向的两个一维DFT所构成。 其中： * . 例题：编程验证二维离散傅里叶变换可分离为两个一维离散傅里叶变换。 解： %myseparable.m %该程序验证了二维DFT的可分离性质 %该程序产生了冈萨雷斯《数字图像处理》（第二版） %P125 图4.4 * . f=imread(D:\chenpc\data\thry\chpt4\Fig4.04(a).jpg); subplot(211) imshow(f,[]) title(原图) F=fftshift(fft2(f)); subplot(223) imshow(log(1+abs(F)),[]) title(用fft2实现二维离散傅里叶变换) [m,n]=size(f); F=fft(f); %沿x方向求离散傅里叶变换 G=fft(F); %沿y方向求离散傅里叶变换 F=fftshift(G); subplot(224) imshow(log(1+abs(F)),[]) title(用fft实现二维离散傅里叶变换) * . * . 平移性 证明： （1）频域移位  * . 结论： 即如果需要将频域的坐标原点从显示屏起始点（0，0）移至显示屏的中心点只要将f(x,y)乘以(-1)x+y因子再进行傅里叶变换即可实现。 例题：利用(-1)x+y对单缝图像f(x,y)进行调制，实现把频谱坐标原点移至屏幕正中央的目标。 当 * . 解：完成本题的源程序为： %在傅里叶变换之前，把函数乘以(-1) x+y，相当于把频谱 %坐标原点移至屏幕窗口正中央。 f(512,512)=0; f=mat2gray(f); [Y,X]=meshgrid(1:512,1:512); f(246:266,230:276)=1; g=f.*(-1).^(X+Y); subplot(221),imshow(f,[]),title(原图像f(x,y)) subplot(222),imshow(g,[]),title(空域调制图像g(x,y)=f(x,y)*(-1)^{x+y})