资源描述:
《数字图像实验报告一图像的加噪处理与几何变换》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验一图像的加噪处理与几何变换一、实验目的1•给Lena图像加高斯噪声,椒盐噪声,均匀分布噪声,观察图像;用平滑滤波器(均值滤波器,中值滤波器)分析效果。2•对lena图像作裁剪、放大、缩小、旋转、平移等几何变换。二、实验内容1•采用中值滤波、均值滤波对受椒盐噪声干扰的图像滤波;2•采用中值滤波、均值滤波对受高斯噪声干扰的图像滤波;3•采用中值滤波、均值滤波对受均匀噪声干扰的图像滤波;4•将图像lena.bmp裁剪成200X200大小;5.制作动画,将一幅图像逐渐向左上角平移移岀图像区域,空白的地
2、方用白色填充;6•利用剪切图像函数制作动画;7•将图像分别放大1.5倍和缩小0.8倍,插值方法使用双线性插值法,分别显示图像;9•将图像分别进行水平方向30度错切,垂直方向45度错切,分别显示结果。8•将图像水平镜像,再顺时针旋转45度,显示旋转后的图;三、实验步骤(一)采用中值滤波,均值滤波对受椒盐噪声干扰的图像滤波a=imread(实验报告数字图像处理实验报告lena.jpg*);%读取图像b=rgb2gray(a);%转化为灰度图像%给图像加入噪声I=imnoise(b,'salt&pep
3、per*)%椒盐噪声%扩展矩阵,生成待处理矩阵n二3;%模板阶数m=(n-1)/2;[p,q]=size(l);PI=zeros(p+2*m,q+2*m);%待处理矩阵fori=1:pforj=1:qPl(i+m,j+m)=l(i,j);endendfori=1:pforii=1:mPl(i+m,ii)=l(i,1);Pl(i+m,q+m+ii)=l(i,q);endendforj=1:qforjj=1:mPl(jj,j+m)=l(1,j);Pl(p+m+jj,j+m)=l(p,j);endend
4、forii=1:mforjj=1:mPl(ii,jj)=1(1,1);Pl(q+m+ii,jj)=l(p,1);Pl(ii,p+m+jj)=1(1,q);Pl(q+m+ii,p+m+jj)=l(p,q);endend%中值滤波&均值滤波derta=zeros(n,n);%n阶模板矩阵Plmid=PI;Plmean=PI;fori=m+1:p+mforj=m+1:q+mfork=1:mforh=1:mderta(k,h)=Pl(i+k-m-1,j+h-m-1);derta(k,m+1)=Pl(i+k
5、-m-1,j);derta(k,n-h+1)=Pl(i+k-m-1,j+m+1-h);derta(m+1,h)=Pl(i,j+h-m-1);derta(m+1,m+1)=Pl(i,j);derta(m+1,n-h+1)=Pl(i,j+m+1-h);derta(n-k+1,h)=Pl(i+m+1-k,j+h-m-1);derta(n-k+1,m+1)=Pl(i+m+1-k,j);derta(n-k+1,n-h+1)=Pl(i+m+1-k,j+m+1-h);Plmid(ij)=median(media
6、n(derta));%中栓滤波Plmean(ij)=round(mean(mean(derta)));%均值滤波endendendend%输出结果Imid=zeros(p,q);Imean=zeros(p,q);fori=1:pforj=1:qImid(ij)=Plmid(i+mj+m);Imean(ij)=Plmean(i+mJ+m);endend%显示结果figure()subplot(2,2,1);imshow(b);title('原图像');subplot(2,2,2);imshow(l)
7、;title('加入椒盐噪声的图像*);subplot(2,2,3);imshow(lmid,[0,255]);title('中值滤波处理后的图像*);subplot(2,2,4);imshow(lmean,[0,255]);title('均值滤波处理后的图像*);生成图像如下:Figure1原團像加入椒盐嗓声的图傑中值腌波处理后的图像均值淤液处理后的图僕由图可见,对于椒盐噪声,中值滤波效果更好。(二)采用中值滤波,均值滤波对受高斯噪声干扰的图像滤波FileEditViewInsertToolsD
8、esktopWindowHelpd3*>w^4?□0■Qa=imread(实验报告数字图像处理实验报告lena.jpg*);%读取图像b=rgb2gray(a);%转化为灰度图像%给图像加入噪声l=imnoise(b,'gaussian',0.02);%高斯噪声%扩展矩阵,生成待处理矩阵n=3;%模板阶数m=(n-1)/2;[p,q]=size⑴;PI=zeros(p+2*m,q+2*m);%待处理矩阵fori=1:pforj=1:qPl(i+m,j+m)=l(i,j);ende