数字图像处理_旋转与幅度谱(含MATLAB代码) .pdfVIP

数字图像处理实验一 15 生医 一、实验内容 产生右图所示图像 f1(m,n) ，其中图像大小为 256 ×256，中间亮条为 128×32，暗处 =0，亮处 =100。 对其进行 FFT： ① 同屏显示原图 f1(m,n) 和 FFT(f1) 的幅度谱图； ② 若令 f2(m,n)= （-1 ）^ （m+n）f1(m,n) ，重复 以上过程，比较二者幅度谱的异同，简述理由； ③ 若将 f2(m,n) 顺时针旋转 90 度得到 f3(m,n) ，试显示 FFT(f3) 的 幅度谱，并与 FFT(f2) 的幅度谱进行比较； ④ 若将 f1(m,n) 顺时针旋转 90 度得到 f4(m,n) ，令 f5(m,n) = f1(m,n) + f4(m,n) ，试显示 FFT(f5) 的幅度谱，指出其与 FFT(f1) 和 FFT(f4) 的关系； ⑤ 若令 f6(m,n)=f2(m,n)+f3(m,n) ，试显示 FFT(f6) 的幅度谱， 并指 出其与 FFT(f2) 和 FFT(f3) 的关系，比较 FFT(f6) 和 FFT(f5) 的幅度 谱。 二、运行环境 MATLAB R2014a 三、运行结果及分析 1. 同屏显示原图 f1(m,n) 和 FFT(f1) 的幅度谱图： 50 50 100 100 150 150 200 200 250 250 50 100 150 200 250 50 100 150 200 250 (a)亮 块 图像 f1(m,n) (b) f1(m,n)的 频 谱 2. 令 f2(m,n)= （-1 ）^ （m+n）f1(m,n) ，对其进行 FFT，比较 f2 与 f1 幅度谱的异同，简述理由： 50 50 100 100 150 150 200 200 250 250 50 100 150 200 250 50 100 150 200 250 (a)亮 块 图像 f2(m,n) (b) f2(m,n)的 频 谱 异同及理由：①空域： f2 由于前边乘了系数（ -1 ）^ （m+n），导致灰 度值有正有负，而在 MATLAB的 imshow 函数中默认把负值变为 0 （有 些情况是取反），所以形成了如左图所示的黑白花纹。 ②频域：FFT(2) 为 FFT(1) 中心化后的图像。空域进行乘以（