数字图像处理上机作业五

数字图像处理上机作业五

ID:5579970

大小:1.69 MB

页数:12页

时间:2017-12-19

数字图像处理上机作业五_第1页
数字图像处理上机作业五_第2页
数字图像处理上机作业五_第3页
数字图像处理上机作业五_第4页
数字图像处理上机作业五_第5页
资源描述:

《数字图像处理上机作业五》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、数字图像第四讲作业数字图像第四讲作业1.设计一个程序对受到高斯白噪声及椒盐噪声干扰的图像进行3x3,5x5邻域的平均平滑以及中值滤波.(添加噪声参看imnoise函数,空域卷积可用imfilter2函数实现)。分析:1.邻域平均平滑可以采用imfilter函数,选择正确的卷积核就可以进行相应的邻域平均平滑操作了。3x3的卷积核为:H1=1/8*[111101111];5x5的卷积核为:H2=1/24*[1111111111110111111111111];2.中值平滑可以先编写中值平滑子函数zhongzhi

2、(),然后在主函数中调用即可。以3*3中值平滑为例来分析其操作过程,3*3中值平滑就是将以各项素为中心的9个像素值的中间值作为平滑后的新的像素值赋给该像素。因此可以通过I(i-1:i+1,j-1:j+1)得到对应于I(i,j)点的九个像素值,然后在由median函数可求出这九个值的中值,赋给新矩阵的(i,j)点即可。注意I(i-1:i+1,j-1:j+1)操作可能会有i-1=0,j-1=0或i+1、j+1大于矩阵最大行列数的情况,从而出现错误。在这里我的处理是在I矩阵的外围补上一圈0,即出现上述情况时像素值

3、以0来代替。具体代码为:I0=zeros(m+2,n+2);fori=2:m+1forj=2:n+1I0(i,j)=I(i-1,j-1);endend同理,5*5的中值平滑也可以同样操作,只不过是在外围补上两圈零而已。Page12of12数字图像第四讲作业代码及注释如下:主函数:clearI=imread('Lenna.bmp');J=imnoise(I,'gaussian');K=imnoise(I,'salt&pepper');%H1为3*3邻域平滑的卷积核,H2为5*5邻域平滑的卷积核H1=1/8*[

4、111101111];H2=1/24*[1111111111110111111111111];J1=imfilter(J,H1);%高斯白噪声的3*3邻域平滑J2=imfilter(J,H2);%高斯白噪声的5*5邻域平滑K1=imfilter(K,H1);%椒盐噪声的3*3邻域平滑K2=imfilter(K,H2);%椒盐噪声的5*5邻域平滑J3=zhongzhi(J,3);%高斯白噪声的3*3中值平滑J4=zhongzhi(J,5);%高斯白噪声的5*5中值平滑K3=zhongzhi(K,3);%椒盐噪

5、声的3*3中值平滑K4=zhongzhi(K,5);%椒盐噪声的5*5中值平滑subplot(131);imshow(J);title('高斯白噪声');subplot(132);imshow(J1);title('高斯白噪声的3*3邻域平滑');subplot(133);imshow(J2);title('高斯白噪声的5*5邻域平滑');figuresubplot(131);imshow(J);title('高斯白噪声');subplot(132);imshow(J3);title('高斯白噪声的3*3中

6、值平滑');subplot(133);imshow(J4);title('高斯白噪声的5*5中值平滑');figuresubplot(131);imshow(K);title('椒盐噪声');subplot(132);imshow(K1);title('椒盐噪声的3*3邻域平滑');subplot(133);imshow(K2);title('椒盐噪声的5*5邻域平滑');figuresubplot(131);imshow(K);title('椒盐噪声');subplot(132);imshow(K3);t

7、itle('椒盐噪声的3*3中值平滑');subplot(133);imshow(K4);title('椒盐噪声的5*5中值平滑');中值平滑子函数zhongzhi()如下:Page12of12数字图像第四讲作业functionJ=zhongzhi(I,k)[m,n]=size(I);ifk==3%3*3的中值平滑I0=zeros(m+2,n+2);fori=2:m+1forj=2:n+1I0(i,j)=I(i-1,j-1);%将到操作的图像矩阵I外围不上0endendfori=2:m+1forj=2:n+

8、1a=I0(i-1:i+1,j-1:j+1);b=a(1:9);%将3*3的矩阵化成1*9的矩阵,便于median操作J(i-1,j-1)=median(b);%取中值,保存为平滑后矩阵J的i-1行、j-1列endendelsek==5%5*5的中值平滑I0=zeros(m+4,n+4);fori=3:m+2forj=3:n+2I0(i,j)=I(i-2,j-2);endendfori=3:m+2forj=3:n

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。