数字图像处理实验八图像复原.pptxVIP

会计学;以高斯噪声为例;MATLAB中叠加噪声的函数：;第4页/共23页;第5页/共23页;周期噪声实例;二、白噪声引起退化图像复原的基本方法;第8页/共23页;2.自适应中值滤波器：;a=imread(circuit.bmp);a=a(:,:,3); [ra,ca]=size(a); a1=imnoise(a,salt pepper,0.25); a1=double(a1); %选定最小滤波器的大小为st×st，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(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(有严重椒盐噪声的输入);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); 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);%排序后的最小像素的值为灰度的最小值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 zmaxmax=sortSxymax((st+4).*(st+4));%第(st+4).*(st+4)个像素的值是最大值zmaxmax; if zminminzminmedzminmedzminmax if zminminzxyzxyzminmax g(n,m)=zxy; else g(n,m)=zminmed; end else if zmedminzmedmedzmedmedzmedmax if zmedminzxyzxyzmedmax g(n,m)=zxy; else g(n,m)=zmedmed; end else