第七讲 图像的傅里叶变换.pptVIP

假设f和h分别由A个和B个点组成。对两个函数同时添加零，以使它们具有相同的周期，表示为P。这个过程产生扩展的或延拓的函数，如下所示。卷积的独立周期将会混叠，这种现象通常归于缠绕误差。 （a） 在[0 N-1]周期中有两个背靠背半周期 （b） 同一区间内有一个完整的周期 这就意味着，坐标原点移到了频谱图像的中间位置，这一点十分重要，尤其是对以后的图像显示和滤波处理。 例题：利用(-1)x对f(x)曲线进行调制，达到平移频域坐标原点至屏幕正中央的目的。 %以一维情况为例，说明空域调制对应着频域坐标原点移位。 f(1:512)=0; f(251:260)=1; %产生宽度为10的窗口函数 subplot(221),plot(f),title(宽度为10 的窗口函数) F=fft(f,512); %进行快速傅里叶变换，延拓周期周期为512 subplot(222) plot(abs(F)) %绘幅度频谱（频谱坐标原点在左边界处） title(幅度谱（频谱坐标原点在左边界处）) x=251:260; f(251:260)=(-1).^x; %把曲线f(x)乘以(-1)^x，可以把频谱 %坐标原点移至屏幕正中央 subplot(223),plot(f),title(宽度为10 的调制窗口函数) F=fft(f,512); %进行快速傅里叶变换 subplot(224); plot(abs(F)) %直接显示幅度频谱（频谱坐标原点在正中央） title(幅度谱（频谱坐标原点在中央）) figure f(1:512)=0; f(251:270)=1; %产生宽度为20的窗口函数 subplot(221),plot(f),title(宽度为20 的窗口函数) F=fft(f,512); %进行快速傅里叶变换，延拓周期周期为512 subplot(222) plot(abs(F)) %绘幅度频谱（频谱坐标原点在左边界处） title(幅度谱（频谱坐标原点在左边界处）) x=251:270; f(251:270)=(-1).^x; %把曲线f(x)乘以(-1)^x，可以把频谱坐标原点移至屏幕正中央 subplot(223),plot(f),title(宽度为20 的调制窗口函数) F=fft(f,512); %进行快速傅里叶变换 subplot(224); plot(abs(F)) %直接显示幅度频谱（频谱坐标原点在正中央） title(幅度谱（频谱坐标原点在中央）) （2）空域移位： 周期性和共轭对称性 周期性： 共轭对称性： Periodicity and conjugate symmetry The discrete Fourier transform has the following periodicity properties: F(u, v) = F(u+M, v) = F(u, v+N) = F(u+M, v+N) The inverse transform also is periodic: f(x, y) = f(x+M, y) = f(x, y+N) = f(x+M, y+N) The idea of conjugate symmetry was introduced in previous section, and is repeated here for convenience: F(u, v) = F*(-u, -v) The spectrum also is symmetric about the origin: Conjugate symmetry 证明： （1）周期性： (2) 共轭对称性： 旋转不变性 证明： 注：为看清问题的实质、简化旋转不变性的证明，以上用二维连续傅里叶变换进行证明。实际上，由连续积分公式进行离散化处理，即可得到离散公式，证明可参照连续情况进行。 f=zeros(512,512); f(246:266,230:276)=1; subplot(221); imshow(f,[]) title(原图) F=fftshift(fft2(f)); subplot(222); imshow(log(1+abs(F)),[]) title(原图的频谱) f=imrotate(f,45,bilinear,crop); subplot(223) imshow(f,[]) title(旋转45^0图) Fc=fftshift(fft2(f)); subplot(224); imsho