频域增强算法实现(共2页).doc

《频域增强算法实现(共2页).doc》由会员分享，可在线阅读，更多相关《频域增强算法实现(共2页).doc（2页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上南京信息工程大学 实验（实习）报告实验（实习）名称 频域增强算法实现 实验（实习）日期 * 得分 指导教师 * 学院 * 专业 * 年级 * 班次 * 姓名 * 学号 * 一、实验目的：1了解离散傅里叶变换的基本性质； 2熟练掌握图像傅里叶变换的方法及应用； 3通过实验了解二维频谱的分布特点； 4熟悉图像频域处理的意义和手段；5通过本实验掌握利用MATLAB的工具箱实现数字图像的频域处理。二、实验内容：1. 编程实现理想低通滤波器、二阶Butterworth低通滤波器、指数低通滤波器和梯形低通滤波器。要求：读入一幅图像，添加椒盐噪声，用上述四种滤波器进行滤波2. 对上

2、述加噪图像进行频域高斯低通和高斯高通滤波，并观察不同滤波半径 (5 100) 条件下滤波图像的变化（半径步长为10）；同屏显示原图、幅度谱图和低通、高通滤波的结果图。三、实验要求：1. 用Matlab语言进行编程，实现上述功能，并尽量使得程序具有通用性；2. 撰写实验报告并附上所用程序和结果。四、实验代码及实验结果1. 低通滤波-平滑clear;clc;I=imread(cameraman.tif);subplot(131),imshow(I);J1=imnoise(I,salt & pepper);subplot(132),imshow(J1);F = double(J1);G = fft2

3、(F);G= fftshift(G);M,N=size(G);nn = 2;d0 = 50;m = fix(M/2);n = fix(N/2);for i = 1 : M for j = 1 : N d = sqrt(i-m)2+(j-n)2); h = 1/(1+(d/d0)(2*nn); result(i,j) = h * G(i,j); end;end;result = ifftshift(result);J2 = ifft2(result);J3 = uint8(real(J2);subplot(133),imshow(J3);2. 高通滤波- 锐化clear;clc;I=imread

4、(cameraman.tif);subplot(121),imshow(I);F= double(I);G = fft2(F);G = fftshift(G);M,N=size(G);nn = 2;d0 = 5;m = fix(M/2);n = fix(N/2);for i = 1 : M for j = 1 : N d = sqrt(i-m)2+(j-n)2); if (d = 0) h = 0; else h=1/(1+(d0/d)(2*nn); end; result(i,j) = h * G(i,j); end;end;result = ifftshift(result);J2= ifft2(result);J3= uint8(real(J2);subplot(122),imshow(J3);五、实验总结通过这次实验，了解了离散傅里叶变换的基本性质，掌握利用MATLAB的工具箱实现数字图像的频域处理，并且实现了理想低通滤波器、二阶Butterworth低通滤波器、指数低通滤波器和梯形低通滤波器的算法实践。专心-专注-专业