- 1、本文档共4页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验五 二维DFT和频域滤波器实验报告
实验目的
实验目的
理解二维离散傅立叶变换的概念,掌握 Matlab中进行二维离散傅立叶变换并显示频谱图的方法。
理解频域滤波的基本概念,掌握频域滤波的基本方法。
理解理想低通滤波器、巴特沃斯低通滤波器的概念,并会利用其传递函数对已知图像进行滤波。
实验要求
实验要求
在Matlab中对给定图像进行二维离散傅立叶变换,并显示其居中的频谱图。
分别生成尺寸为 M x N (给定图像的尺寸),截止频率D0=50的理想低通滤波器和 1阶巴特沃斯低通滤 波器。并以图像的形式显示。
对给定图像分别利用上述滤波器进行滤波,显示滤波后图像的频谱以及经滤波后的图像,并与原始图 像及其频谱进行对比。
三、实验原理
根据卷积定理,两个空间函数的卷积可以通过计算两个傅立叶变换函数的乘积的逆变换得到,如果 f(x,
y)和h(x, y)分别代表图像与空间滤波器,F(u, v)和H(u, v)分别为响应的傅立叶变换 (H(u, v)又称为传递函数), 那么我们可以利用卷积定理来进行频域滤波。
理想低通滤波器具有传递函数
H(u,v)1 若 D(u,v) Do
H(u,v)
0 若 D(u,v) D0
其中D。为制定的非负数,D(u, v)为点(u, v)到滤波器中心的距离。 n阶巴特沃斯低通滤波器的传递函数为
H(u,v)
1
1 [D(u, v)/D°]2n
四、实验结果图像
理弊假询龙巧器茁舖谱图
匚贞询理匸恨诵总涯后的频谓
五、实验代码
%原始图片及其居中后的频谱:
J=imread( d:\lena.bmp );
J_gray=rgb2gray(J);
figure;
imshow(J_gray); %显示原始灰度图像
title(原始灰度图像’);
J2=fft2(J_gray); %二维离散傅立叶变换
J3=fftshift(J2); %居中
J4=log(abs(J3)); %对图像做二维DFT变换并居中
figure;
imshow(J4,[]); %显示居中后的频谱图
title(离散傅立叶变换居中频谱图 ’);
%理想低通滤波器:
[m,n ]=size(J_gray);
d0=50;
xO=rou nd(256);
yO=rou nd(256); %定义一个半径为 256的圆形区域
%截止频率d0=50的理想低通滤波器
for i=1:m
for j=1: n
d=sqrt((i-x0)A2+(j-y0)A2); %到原点中心的距离
if (dd0)
h(i,j)=1;
else
h(i,j)=0;
end
end
end
figure;
imshow(h); %画出理想低通滤波器的频谱图 title( 理想低通滤波器的频谱图 );
J5=J3.*h;
J6=log(abs(J5));
figure;
imshow(J6,[]) ; %画出滤波后的频谱
title( 原频谱理想低通滤波后的频谱 );
J7=ifftshift(J5);
J8=ifft2(J7);
J9=uint8(real(J8)); %对滤波后的频谱做反变换后居中
figure;
imshow(J9); %画出经过理想低通滤波器滤波后的图像 title( 经过理想低通滤波器滤波后的图像 );
%巴特沃斯滤波器:
[m,n]=size(J_gray); d0=50;n=1;
x0=round(256);
y0=round(256);
for i=1:m
for j=1:n
d=sqrt((i-x0)A2+(j-y0)A2); h(i,j)=1/(1+(d/dOF(2* n));
end
end
figure;
imshow(h); %画出巴特沃斯滤波器的频谱图 title( 巴特沃斯滤波器的频谱图 );
J10=J3.*h;
J11=log(abs(J10));
figure;
imshow(J11,[]); %画出滤波后的频谱
title( 巴特沃斯滤波器滤波后的频谱 );
J12=ifftshift(J10);
J13=ifft2(J12);
J14=uint8(real(J13)); %对巴特沃斯滤波后的图像进行反变换并居中
figure;
imshow(J14); %画出经过巴特沃斯滤波器滤波后的图像
title( 经过巴特沃斯滤波器滤波后的图像 );
您可能关注的文档
最近下载
- (2023)GB2894安全标志及使用导则标准大全.pdf
- 云-优加-菜鸟裹裹驿站1013出海0930在线岗前考试1011.docx VIP
- 芳香保健理论知识考试题库(含答案).pdf VIP
- 部编人教版小学道德与法治三年级上册:全册表格式教案.pdf
- 田家四季歌(第二课时) 逐字稿 二年级上册语文 智慧中小学.pdf
- 中级(监控类) 消防设施操作员理论考试题及答案.doc VIP
- 2024年秋季新湘教版七年级上册地理全册教案.docx
- 电缆生产企业安全风险分级管控(危险源辨识分析)资料.docx VIP
- 天天特卖售后结业售后场景考试 - 副本.doc VIP
- 铝镁锰金属屋面施工方案.docx
文档评论(0)