- 1、本文档共26页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
一、编写程序完成不同滤波器的图像频域降噪和边缘增强的算法并进行比较,得出结论。
1、不同滤波器的频域降噪
1.1 理想低通滤波器(ILPF)
I1=imread(eight.tif); %读取图像
I2=im2double(I1);
I3=imnoise(I2,gaussian,0.01);
I4=imnoise(I3,salt pepper,0.01);
figure,subplot(1,3,1);
imshow(I2) %显示灰度图像
title(原始图像); %为图像添加标题
subplot(1,3,2);
imshow(I4) %加入混合躁声后显示图像
title(加噪后的图像);
s=fftshift(fft2(I4)); %将灰度图像的二维不连续Fourier 变换的零频率成分移到频谱的中心
[M,N]=size(s); %分别返回s的行数到M中,列数到N中
n1=floor(M/2); %对M/2进行取整
n2=floor(N/2); %对N/2进行取整
d0=40; %初始化d0
for i=1:M
for j=1:N
d=sqrt((i-n1)^2+(j-n2)^2); %点(i,j)到傅立叶变换中心的距离
if d=d0 %点(i,j)在通带内的情况
h=1; %通带变换函数
else %点(i,j)在阻带内的情况
h=0; %阻带变换函数
end
s(i,j)=h*s(i,j); %ILPF滤波后的频域表示
end
end
s=ifftshift(s); %对s进行反FFT移动
s=im2uint8(real(ifft2(s))); %对s进行二维反离散的Fourier变换后,取复
数的实部转化为无符号8位整数
subplot(1,3,3); %创建图形图像对象
imshow(s); %显示ILPF滤波后的图像
title(ILPF滤波后的图像(d=40));
运行结果:
1.2 二阶巴特沃斯低通滤波器(BLPF)
I1=imread(eight.tif); %读取图像
I2=im2double(I1);
I3=imnoise(I2,gaussian,0.01);
I4=imnoise(I3,salt pepper,0.01);
figure,subplot(1,3,1);
imshow(I2) %显示灰度图像
title(原始图像); %为图像添加标题
subplot(1,3,2);
imshow(I4) %加入混合躁声后显示图像
title(加噪后的图像);
s=fftshift(fft2(I4));%将灰度图像的二维不连续Fourier 变换的零频率成分移到频谱的中心
[M,N]=size(s); %分别返回s的行数到M中,列数到N中
n=2; %对n赋初值
n1=floor(M/2); %对M/2进行取整
n2=floor(N/2); %对N/2进行取整
d0=40; %初始化d0
for i=1:M
for j=1:N
d=sqrt((i-n1)^2+(j-n2)^2); %点(i,
文档评论(0)