數(shù)字圖像處理噪聲處理和均值濾波和中值慮波

1、數(shù)字圖像處理噪聲處理和均值濾波和中值慮波一、第一題程序：I=imread('01.jpg');gray=rgb2gray(I);subplot(231);imshow(gray,),title('輸入圖像'),ra,ca=size(gray);r=0:255;lc=gray(:); %將a按列拉長(zhǎng)p=hist(lc,r); %統(tǒng)計(jì)各個(gè)灰度值出現(xiàn)的次數(shù)pr=p/length(lc);subplot(234);bar(pr);title('輸入圖像直方圖');gaussian_noise=imnoise(gray,'gaussian'

2、,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

3、=hist(cc,r);psr=ps/length(cc);subplot(236);bar(psr);title('椒鹽噪聲直方圖');二、結(jié)果：三、第二題程序：%演示算術(shù)均值濾波器a=imread('111.jpg');a=a(:,:,3);ra,ca=size(a);a1=double(imnoise(a,'gaussian',0,0.01);%疊加嚴(yán)重的高斯噪聲 %選定濾波器的大小為st×st，st可調(diào)st=3,h=ones(st,st);g=filter2(h, a1)./sum(sum(h);subplot(2,3,2),i

4、mshow(a1,),title('有高斯噪聲的輸入圖像'),subplot(2,3,5),imshow(g,),title('算術(shù)均值濾波后的圖像，用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可調(diào),但取奇數(shù)st=3,Sxymin=zeros(st,st);%預(yù)先留出鄰域%選定第二濾波器的大小為（st+2

5、）×（st+2）Sxymed=zeros(st+2,st+2);%預(yù)先留出鄰域%選定第二濾波器的大小為（st+4）×（st+4）Sxymax=zeros(st+4,st+4);%預(yù)先留出鄰域 a=zeros(ra+2.*st,ca+2.*st);%預(yù)先留出一個(gè)圖像a(st+1:ra+st,st+1:ca+st)=a1;%將輸入圖像放大到（ra+2.*st）×（ca+2.*st）g=zeros(ra+2.*st,ca+2.*st);%預(yù)先留出輸出圖像subplot(2,3,1),imshow(a,),title('有嚴(yán)重椒鹽噪聲的輸入')for n=

6、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);%賦掩模內(nèi)的值 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=sor

7、tSxymin(1,1);%排序后的最小像素的值為灰度的最小值z(mì)minmin zmedmin=sortSxymed(1,1);%排序后的最小像素的值為灰度的最小值z(mì)medmin zmaxmin=sortSxymax(1,1);%排序后的最小像素的值為灰度的最小值z(mì)maxmin zminmed=sortSxymin(st+2);%排序后的第st+2個(gè)像素的值是中值z(mì)minmed zmedmed=sortSxymed(st.*4+1);%排序后的第st.*4+1個(gè)像素的值是中值z(mì)medmed zmaxmed=sortSxymax(st.*8+1);%排序后的第st.*8+1個(gè)像素的值是中值z(mì)max

8、med zminmax=sortSxymin(st.*st);%第st.*st個(gè)像素的值是最大值z(mì)minmax zmedmax=sortSxymed(st+2).*(st+2);%第(st+2).*(st+2)個(gè)像素的值是最大值z(mì)medmax zmaxmax=sortSxymax(st+4).*(st+4);%第(st+4).*(st+4)個(gè)像素的值是最大值z(mì)maxmaxif zminmin<zminmed&zminmed<zminmax if zminmin<zxy&zxy<zminmax g(n,m)=zxy; else g(n,m)=zminmed

9、; end else if zmedmin<zmedmed&zmedmed<zmedmax if zmedmin<zxy&zxy<zmedmax g(n,m)=zxy; else g(n,m)=zmedmed; end else if zmaxmin<zmaxmed&zmaxmed<zmaxmax if zmaxmin<zxy&zxy<zmaxmax g(n,m)=zxy; else g(n,m)=zmaxmed; end else g(n,m)=zxy; end end end endendsubplot(2,3

10、,4),imshow(g,),title('自適應(yīng)中值濾波器的處理結(jié)果')a=imread('20.jpg');subplot(2,3,3),imshow(a,);title('斜條紋')fxy=a(:,:,1);fxy=double(fxy);ra,ca=size(fxy); Fuv=fftshift(fft2(fxy);u=1:ra;v=1:ca;u,v=meshgrid(u,v);u0=129-99;%噪聲在頻譜面上的中心v0=129-99;%噪聲在頻譜面上的中心D0=30%陷波區(qū)域的半徑Duv1=(u-ra/2-u0).2+(v-ca/2-v0).2).0.5;Duv2=(u-ra/2+u0).2+(v-ca/2+v0).2).0.5;%下面使用高斯陷波濾波器Huv3=zeros(ra,ca