数字图像处理噪声处理与均值滤波与中值虑波.docVIP

数字图像处理噪声处理和均值滤波和中值虑波 一、第一题程序： I=imread(01.jpg); gray=rgb2gray(I); subplot(231);imshow(gray,[]),title(输入图像), [ra,ca]=size(gray); r=0:255; lc=gray(:); %将a按列拉长 p=hist(lc,r); %统计各个灰度值出现的次数 pr=p/length(lc); subplot(234);bar(pr);title(输入图像直方图); gaussian_noise=imnoise(gray,gaussian,0,0.001); subplot(232);imshow(gaussian_noise,[]);title(高斯噪声); bb=gaussian_noise(:); pg=hist(bb,r); pgr=pg/length(bb); subplot(235);bar(pgr);title(高斯噪声直方图); salt_noise=imnoise(gray,salt pepper,0.1); subplot(233);imshow(salt_noise,[]);title(椒盐噪声); cc=salt_noise(:); ps=hist(cc,r); psr=ps/length(cc); subplot(236);bar(psr);title(椒盐噪声直方图); 二、结果： 三、第二题程序： %演示算术均值滤波器 a=imread(111.jpg);a=a(:,:,3); [ra,ca]=size(a); a1=double(imnoise(a,gaussian,0,0.01));%叠加严重的高斯噪声 %选定滤波器的大小为st×st，st可调 st=3, h=ones(st,st); g=filter2(h, a1)./sum(sum(h)); subplot(2,3,2),imshow(a1,[]),title(有高斯噪声的输入图像), subplot(2,3,5),imshow(g,[]),title(算术均值滤波后的图像，用3×3掩模) a=imread(111.jpg);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(2,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=sort