MATLAB图像邻域运算去除噪声模板运算中值滤波均值滤波.docx
- 文档编号:23317095
- 上传时间:2023-05-16
- 格式:DOCX
- 页数:16
- 大小:862.47KB
MATLAB图像邻域运算去除噪声模板运算中值滤波均值滤波.docx
《MATLAB图像邻域运算去除噪声模板运算中值滤波均值滤波.docx》由会员分享,可在线阅读,更多相关《MATLAB图像邻域运算去除噪声模板运算中值滤波均值滤波.docx(16页珍藏版)》请在冰豆网上搜索。
MATLAB图像邻域运算去除噪声模板运算中值滤波均值滤波
1.列出常用的卷积模板
2.基于3×3的模板,编写均值滤波的处理程序,处理含有加性高斯噪声和椒盐噪声的图像,观察处理结果
3.编写中值滤波程序,处理相同的图像与均值滤波进行比较;改变模板尺寸观察处理结果
4.编程实现利用一阶微分算子和二阶拉普拉斯算子进行图像锐化的程序
5.对比不同的邻域运算结果,体会图像锐化与图像平滑的区别
均值滤波处理含有椒盐噪声图像程序代码
(1):
模板运算
f=imread('lenasp.bmp');
f=double(f);
[row,col]=size(f);
r=1;
fori=2:
row-1
forj=2:
col-1
g(i,j)=(f(i-1,j-1)+f(i-1,j)+f(i-1,j+1)+f(i,j-1)+f(i,j)+f(i,j+1)+f(i+1,j-1)+f(i+1,j)+f(i+1,j+1))/9;
g=uint8(g);
end
end
subplot(121);
imshow(uint8(f));
title('椒盐噪声图像');
subplot(122);
imshow(g);
title('均值滤波处理后的图像')
(2):
mean函数均值滤波
f=imread('lenasp.bmp');
f=double(f);
[row,col]=size(f);
g=size(f);
fori=2:
row-1
forj=2:
col-1
t=f(i-1:
i+1,j-1:
j+1);
t=double(t);
g(i,j)=mean(mean(t));
g(i,j)=uint8(g(i,j));
end
end
subplot(121);
imshow(uint8(f));
title('椒盐噪声图像');
subplot(122);
imshow(uint8(g));
title('均值滤波处理后的图像')
(3)可变模板处理
f=imread('lenagauss.bmp');
f=double(f);
[row,col]=size(f);
r=1;
g=size(f);
fori=2:
row-1
forj=2:
col-1
s=0;
s=double(s);
form=-r:
r
forn=-r:
r
s=s+f(i+m,j+n);
end
end
g(i,j)=s/power((r+2),2);
end
end
subplot(121);
imshow(uint8(f));
title('椒盐噪声图像');
subplot(122);
imshow(uint8(g));
title('均值滤波处理后的图像')
均值滤波处理含有加性高斯噪声图像程序代码
(1):
模板运算
f=imread('lenagauss.bmp');
f=double(f);
[row,col]=size(f);
r=1;
fori=2:
row-1
forj=2:
col-1
g(i,j)=(f(i-1,j-1)+f(i-1,j)+f(i-1,j+1)+f(i,j-1)+f(i,j)+f(i,j+1)+f(i+1,j-1)+f(i+1,j)+f(i+1,j+1))/9;
g=uint8(g);
end
end
subplot(121);
imshow(uint8(f));
title('加性高斯噪声图像');
subplot(122);
imshow(g);
title('均值滤波处理后的图像')
(2):
mean函数均值滤波
f=imread('lenagauss.bmp');
f=double(f);
[row,col]=size(f);
r=1;
g=size(f);
fori=2:
row-1
forj=2:
col-1
t=f(i-1:
i+1,j-1:
j+1);
t=double(t);
g(i,j)=mean(mean(t));
g(i,j)=uint8(g(i,j));
end
end
subplot(121);
imshow(uint8(f));
title('加性高斯噪声图像');
subplot(122);
imshow(uint8(g));
title('均值滤波处理后的图像')
(3):
可变模板处理
f=imread('lenagauss.bmp');
f=double(f);
[row,col]=size(f);
r=1;
g=size(f);
fori=2:
row-1
forj=2:
col-1
s=0;
s=double(s);
form=-r:
r
forn=-r:
r
s=s+f(i+m,j+n);
end
end
g(i,j)=s/power((r+2),2);
end
end
subplot(121);
imshow(uint8(f));
title('加性高斯噪声图像');
subplot(122);
imshow(uint8(g));
title('均值滤波处理后的图像')
中值滤波(椒盐噪声图像处理)
f=imread('lenasp.bmp');
f=double(f);
[row,col]=size(f);
g=size(f);
r=1;
fori=r+1:
row-r
forj=r+1:
col-r
t=f(i-r:
i+r,j-r:
j+r);
id=0;
fortm=1:
5
m=0;
fork=1:
9
if(m m=t(k); id=k; end end t(id)=0; end g(i,j)=m; end end subplot(121); imshow(uint8(f)); title('原图像'); subplot(122); imshow(uint8(g)); title('中值滤波处理后的图像') R=1时的图像 R=5 中值滤波(加性高斯噪声) f=imread('lenagauss.bmp'); f=double(f); [row,col]=size(f); g=size(f); r=1; fori=r+1: row-r forj=r+1: col-r t=f(i-r: i+r,j-r: j+r); id=0; fortm=1: 5 m=0; fork=1: 9 if(m m=t(k); id=k; end end t(id)=0; end g(i,j)=m; end end subplot(121); imshow(uint8(f)); title('原图像'); subplot(122); imshow(uint8(g)); title('中值滤波处理后的图像') R=5 不同领域运算结果比较 r=2 f=imread('lenasp.bmp'); f=double(f); [row,col]=size(f); r=2; g=size(f); fori=r+1: row-r forj=r+1: col-r s=0; s=double(s); form=-r: r forn=-r: r s=s+f(i+m,j+n); end end g(i,j)=s/power((r+2),2); end end subplot(121); imshow(uint8(f)); title('椒盐噪声图像'); subplot(122); imshow(uint8(g)); title('均值滤波处理后的图像') R=3 R=8
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 图像 邻域 运算 去除 噪声 模板 中值 滤波 均值