(完整word版)2018秋西安电子科技大学数字图像处理上机实验题答案及matlab程序

西安电子科技大学

数字图像处理上机实验题

李XX 2018-12-16

数字图像处理上机实验题

1. 产生右图所示图像 f1(m,n)，其中图像大小为256×256，中间亮条为128×32，暗处=0，亮处=100。对其进行FFT：(matlab程序附于文末附录)

① 同屏显示原图f1(m,n)和FFT(f1)的幅度谱图；

② 若令f2(m,n)=(-1)m+n f1(m,n)，重复以上过程，比较二者幅度谱的异同，简述理由；

结论：

相同点：图像f1,f2的幅度谱的实质没有改变，幅度等都没有发生变化。

不同点：f2的频谱是对f1频谱的移位，它是f1的频谱从原点(0,0)移动到中心点(64,64)得到的频谱。

原因：f2图像是对f1图像做频谱中心化变换得来的，这一过程不改变图像的幅度等特性。

③ 若将f2(m,n)顺时针旋转90 度得到f3(m,n)，试显示FFT(f3)的幅度谱，并与FFT(f2)的幅度谱进行比较；

结论：

由图可知，FFT(f3)的幅度谱是由FFT(f2)的幅度谱顺时针旋转90度得到的。

④ 若将f1(m,n) 顺时针旋转90 度得到f4(m,n)，令f5(m,n)=f1(m,n)+f4(m,n)，试显示FFT(f5)的幅度谱，并指出其与 FFT(f1)和FFT(f4)的关系；

结论：

FFT(f5)的幅度谱是FFT(f1)和FFT(f4)幅度谱叠加的结果。

1

⑤ 若令f6(m,n)=f2(m,n)+f3(m,n)，试显示FFT(f6)的幅度谱，并指出其与 FFT(f2)和FFT(f3)

的关系，比较FFT(f6)和FFT(f5)的幅度谱。

结论：

FFT(f6)的幅度谱是FFT(f2)和FFT(f3)幅度谱叠加的结果。FFT(f6)是FFT(f5)频谱中心化的结果。

2. 产生教材104 页题图4.18（右图）所示的二值图像（白为1，黑为0），编程实现习题4.18 所要求的处理（3×3 的平均滤波和中值滤波）功能（图像四周边界不考虑，处理结果按四舍五入仍取0 或1），显示处理前后的图像，比较其异同。(matlab程序附于文末附录)

结论：

图像fa的平均滤波图像和平均加权图像与原图像相同；fb的平均滤波图像的处理点在白、黑块组成的4块交界处时，由于窗内1值点等于0值点个数，0值点的值变为1，形成如图所示图像，加权平均滤波图像与原图像相同。

3. 产生教材104 页题图4.16 所示的灰度图像（白为255，黑为0），分别加入高斯白噪声和椒盐噪声，再分别进行3×3 的平均滤波和中值滤波，显示原图像、加噪图像和滤波结果图像，并比较四种滤波结果。(matlab程序附于文末附录)

2

结论：

由图可知，对于加高斯白噪声后的图像，平均滤波后对噪声的滤除总体上效果好，但在局部边缘处存在误将噪声当作图像保留和误将图像细节当作噪声滤除的现象；中值滤波则对图像的边缘信息保留的较好，但总体滤除效果不佳。对于加椒盐噪声后的图像，平均滤波后在图像边缘上的损失较多，相对的，中值滤波后边缘信息损失较少。

4. 对某一灰度图像，进行如下处理：(matlab程序附于文末附录)

（1）分别利用Roberts、Prewitt 和Sobel 边缘检测算子进行边缘检测；

本题选取此图片进行处理。

3