MATLAB圖像鄰域運算,去除噪聲,模板運算中值濾波,均值濾波

1、數(shù)字圖像處理實驗報告 #領(lǐng)域運算 領(lǐng)域運算列出常用的卷積模板基于3X3的模板，編寫均值濾波的處理程序，處理含有加性高斯噪聲和椒鹽噪聲的圖像，觀察處理結(jié)果編寫中值濾波程序，處理相同的圖像與均值濾波進行比較；改變模板尺寸觀察處理結(jié)果編程實現(xiàn)利用一階微分算子和二階拉普拉斯算子進行圖像銳化的程序?qū)Ρ炔煌泥徲蜻\算結(jié)果，體會圖像銳化與圖像平滑的區(qū)別均值濾波處理含有椒鹽噪聲圖像程序代碼：模板運算f=imread(lenasp.bmp);f=double(f);row,col=size(f);r=1;fori=2:row-1forj=2:col-1g(i,j)=(f(i_1,j_1)+f(i_1,j)+f(

2、i_1,j+1)+f(i,j_1)+f(i,j)+f(i,j+1)+f(i+1,j-1)+f(i+1,j)+f(i+1,j+1)/9;g=uint8(g);endendsubplot(121);imshow(uint8(f);title(椒鹽噪聲圖像);subplot(122);imshow(g);title(均值濾波處理后的圖像)數(shù)字圖像處理實驗報告椒鹽噪聲圖像均值濾波處理后的圖像數(shù)字圖像處理實驗報告椒鹽噪聲圖像均值濾波處理后的圖像 領(lǐng)域運算 #領(lǐng)域運算：mean函數(shù)均值濾波f=imread(lenasp.bmp);f=double(f);row,col=size(f);g=size(f);

3、fori=2:row-1forj=2:col-1t=f(i-1:i+1,j-1:j+1);t=double(t);g(i,j)=mean(mean(t);g(i,j)=uint8(g(i,j);endendsubplot(121);imshow(uint8(f);title(椒鹽噪聲圖像);subplot(122);imshow(uint8(g);title(均值濾波處理后的圖像)數(shù)字圖像處理實驗報告均值濾波處理后的圖像數(shù)字圖像處理實驗報告均值濾波處理后的圖像 #領(lǐng)域運算 領(lǐng)域運算可變模板處理f=imread(lenagauss.bmp);f=double(f);row,col=size(f)

4、;r=1;g=size(f);fori=2:row-1forj=2:col-1s=0;s=double(s);form=-r:rforn=-r:rs=s+f(i+m,j+n);endendg(i,j)=s/power(r+2),2);endendsubplot(121);imshow(uint8(f);title(椒鹽噪聲圖像);subplot(122);imshow(uint8(g);title(均值濾波處理后的圖像)數(shù)字圖像處理實驗報告椒鹽噪聲圖像均值濾波處理后的圖像數(shù)字圖像處理實驗報告椒鹽噪聲圖像均值濾波處理后的圖像 領(lǐng)域運算 #領(lǐng)域運算均值濾波處理含有加性高斯噪聲圖像程序代碼:模板運算

5、f=imread(lenagauss.bmp);f=double(f);row,col=size(f);r=1;fori=2:row-1forj=2:col-1g(i,j)=(f(i_1,j_1)+f(i_1,j)+f(i_1,j+1)+f(i,j_1)+f(i,j)+f(i,j+1)+f(i+1,j-1)+f(i+1,j)+f(i+1,j+1)/9;g=uint8(g);endendsubplot(121);imshow(uint8(f);title(加性高斯噪聲圖像);subplot(122);imshow(g);title(均值濾波處理后的圖像)數(shù)字圖像處理實驗報告加性高斯噪聲圖像均值濾

6、波處理后的圖像數(shù)字圖像處理實驗報告加性高斯噪聲圖像均值濾波處理后的圖像 #領(lǐng)域運算 領(lǐng)域運算：mean函數(shù)均值濾波f=imread(lenagauss.bmp);f=double(f);row,col=size(f);r=1;g=size(f);fori=2:row-1forj=2:col-1t=f(i-1:i+1,j-1:j+1);t=double(t);g(i,j)=mean(mean(t);g(i,j)=uint8(g(i,j);endendsubplot(121);imshow(uint8(f);title(加性高斯噪聲圖像);subplot(122);imshow(uint8(g);

7、title(均值濾波處理后的圖像)數(shù)字圖像處理實驗報告均值濾波處理后的圖像數(shù)字圖像處理實驗報告均值濾波處理后的圖像 領(lǐng)域運算 #領(lǐng)域運算:可變模板處理f=imread(lenagauss.bmp);f=double(f);row,col=size(f);r=1;g=size(f);fori=2:row-1forj=2:col-1s=0;s=double(s);form=-r:rforn=-r:rs=s+f(i+m,j+n);endendg(i,j)=s/power(r+2),endendsubplot(121);imshow(uint8(f);title(加性高斯噪聲圖像);subplot(1

8、22);imshow(uint8(g);數(shù)字圖像處理實驗報告數(shù)字圖像處理實驗報告 領(lǐng)域運算 領(lǐng)域運算title（均值濾波處理后的圖像）椒鹽噪聲圖像均值濾波處理后的圖像中值濾波(椒鹽噪聲圖像處理)f=imread(lenasp.bmp);f=double(f);row,col=size(f);g=size(f);r=1;fori=r+1:row-rforj=r+1:col-rt=f(i-r:i+r,j-r:j+r);id=0;fortm=1:5m=0;fork=1:9if(mt(k)m=t(k);id=k;endendt(id)=0;endg(i,j)=m;endendsubplot(121);

9、imshow(uint8(f);title(原圖像);subplot(122);imshow(uint8(g);title(中值濾波處理后的圖像)R=1時的圖像原圖像中值濾波處理后的圖像R=5原圖像中值濾波處理后的圖像中值濾波(加性高斯噪聲)f=imread(lenagauss.bmp);f=double(f);row,col=size(f);g=size(f);r=1;fori=r+1:row-rforj=r+1:col-rt=f(i-r:i+r,j-r:j+r);id=0;fortm=1:5m=0;fork=1:9if(mt(k)m=t(k);id=k;endendt(id)=0;endg

10、(i,j)=m;endendsubplot(121);imshow(uint8(f)title(原圖像)subplot(122)imshow(uint8(g)title(中值濾波處理后的圖像)原圖像中值濾波處理后的圖像R=5原圖像中值濾波處理后的圖像數(shù)字圖像處理實驗報告均值濾波處理后的圖像數(shù)字圖像處理實驗報告 領(lǐng)域運算 領(lǐng)域運算不同領(lǐng)域運算結(jié)果比較r=2f=imread(lenasp.bmp);f=double(f);row,col=size(f);r=2;g=size(f);fori=r+1:row-rforj=r+1:col-rs=0;s=double(s);form=-r:rforn=-r:rs=s+f(i+m,j+n);end