数字图像处理实验八图像复原.pptx
《数字图像处理实验八图像复原.pptx》由会员分享,可在线阅读,更多相关《数字图像处理实验八图像复原.pptx(23页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、一、噪声的类型和基本特点一、噪声的类型和基本特点噪声噪声随机噪声随机噪声周期噪声周期噪声高斯噪声高斯噪声瑞利噪声瑞利噪声伽马噪声伽马噪声指数噪声指数噪声均匀噪声均匀噪声脉冲噪声脉冲噪声椒盐噪声椒盐噪声盐噪声盐噪声胡椒噪声胡椒噪声第1页/共23页以高斯噪声为例灰度概率分布函数:=120,=30时:灰度概率分布频 谱第2页/共23页MATLAB中叠加噪声的函数:输出图像名=imnoise(输入图像名,噪声类型,噪声参数);加高斯噪声:加高斯噪声:fxynoise1=imnoise(fxy,gaussian,0,0.001);加椒盐噪声:加椒盐噪声:fxynoise2=imnoise(fxy,sal
2、t&pepper,0.4);加均匀噪声:加均匀噪声:fxynoise3=imnoise(fxy,speckle,0.03);第3页/共23页第4页/共23页第5页/共23页周期噪声实例第6页/共23页二、白噪声引起退化图像复原的基本方法w1w2w3w4w5w6w7w8w91.均值滤波:111111111%本例演示算术均值滤波器a=imread(circuit.bmp);a=a(:,:,3);subplot(1,3,1),imshow(a,),title(原始图像),ra,ca=size(a);a1=double(imnoise(a,gaussian,0,0.01);%叠加严重的高斯噪声%选定滤
3、波器的大小为stst,st可调st=3,h=ones(st,st);g=filter2(h,a1)./sum(sum(h);subplot(1,3,2),imshow(a1,),title(有高斯噪声的输入图像),subplot(1,3,3),imshow(g,),title(算术均值滤波后的图像,用33掩模)第7页/共23页第8页/共23页2.自适应中值滤波器:输入待处理图像给出Smax的大小设定St起始大小给出Sxy区域计算Zmin、Zmax、Zmed和A1=Zmed-Zmin,A2=Zmed-ZmaxA10A20B20B1=Zxy-Zmin,B2=Zxy-Zmax增大Sxy是否是是否否检
4、查Zmed是否是冲击噪声检查Zxy是否是冲击噪声第9页/共23页a=imread(circuit.bmp);a=a(:,:,3);ra,ca=size(a);a1=imnoise(a,salt&pepper,0.25);a1=double(a1);%选定最小滤波器的大小为stst,st可调,但取奇数st=3,Sxymin=zeros(st,st);%预先留出邻域%选定第二滤波器的大小为(st+2)(st+2)Sxymed=zeros(st+2,st+2);%预先留出邻域%选定第二滤波器的大小为(st+4)(st+4)Sxymax=zeros(st+4,st+4);%预先留出邻域 a=zeros
5、(ra+2.*st,ca+2.*st);%预先留出一个图像a(st+1:ra+st,st+1:ca+st)=a1;%将输入图像放大到(ra+2.*st)(ca+2.*st)g=zeros(ra+2.*st,ca+2.*st);%预先留出输出图像subplot(1,3,1),imshow(a,),title(有严重椒盐噪声的输入)第10页/共23页for n=st+1:ra+st for m=st+1:ca+st for nn=1:st+4 for mm=1:st+4 Sxymax(nn,mm)=a(n+nn-st-1,m+mm-st-1);%赋掩模内的值 end end zxy=a(n,m);
6、Sxymin=Sxymax(3:st+2,3:st+2);Sxymed=Sxymax(2:st+3,2:st+3);sortSxymin=sort(Sxymin(:);%将像素按升序排列 sortSxymed=sort(Sxymed(:);%将像素按升序排列 sortSxymax=sort(Sxymax(:);%将像素按升序排列 zminmin=sortSxymin(1,1);%排序后的最小像素的值为灰度的最小值zminmin zmedmin=sortSxymed(1,1);%排序后的最小像素的值为灰度的最小值zmedmin zmaxmin=sortSxymax(1,1);%排序后的最小像素的
7、值为灰度的最小值zmaxmin zminmed=sortSxymin(st+2);%排序后的第st+2个像素的值是中值zminmed zmedmed=sortSxymed(st.*4+1);%排序后的第st.*4+1个像素的值是中值zmedmed zmaxmed=sortSxymax(st.*8+1);%排序后的第st.*8+1个像素的值是中值zmaxmed zminmax=sortSxymin(st.*st);%第st.*st个像素的值是最大值zminmax zmedmax=sortSxymed(st+2).*(st+2);%第(st+2).*(st+2)个像素的值是最大值zmedmax z
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字图像 处理 实验 图像 复原
限制150内