數(shù)字圖像處理上機作業(yè)

《數(shù)字圖像處理上機作業(yè)》由會員上傳分享，免費在線閱讀，更多相關內容在應用文檔-天天文庫。

1、數(shù)字圖像第四講作業(yè)數(shù)字圖像第四講作業(yè)1.設計一個程序對受到高斯白噪聲及椒鹽噪聲干擾的圖像進行3x3,5x5鄰域的平均平滑以及中值濾波.(添加噪聲參看imnoise函數(shù),空域卷積可用imfilter2函數(shù)實現(xiàn))。分析：1.鄰域平均平滑可以采用imfilter函數(shù)，選擇正確的卷積核就可以進行相應的鄰域平均平滑操作了。3x3的卷積核為：H1=1/8*[111101111];5x5的卷積核為：H2=1/24*[1111111111110111111111111];2.中值平滑可以先編寫中值平滑子函數(shù)zhongzhi()，然后在主函數(shù)中調用即可。以3*3中值平滑為例來分析其操

2、作過程，3*3中值平滑就是將以各項素為中心的9個像素值的中間值作為平滑后的新的像素值賦給該像素。因此可以通過I(i-1:i+1,j-1:j+1)得到對應于I(i,j)點的九個像素值，然后在由median函數(shù)可求出這九個值的中值，賦給新矩陣的(i,j)點即可。注意I(i-1:i+1,j-1:j+1)操作可能會有i-1=0，j-1=0或i+1、j+1大于矩陣最大行列數(shù)的情況，從而出現(xiàn)錯誤。在這里我的處理是在I矩陣的外圍補上一圈0，即出現(xiàn)上述情況時像素值以0來代替。具體代碼為：I0=zeros(m+2,n+2);fori=2:m+1forj=2:n+1I0(i,j)=I(

3、i-1,j-1);endend同理，5*5的中值平滑也可以同樣操作，只不過是在外圍補上兩圈零而已。Page12of12數(shù)字圖像第四講作業(yè)代碼及注釋如下：主函數(shù)：clearI=imread('Lenna.bmp');J=imnoise(I,'gaussian');K=imnoise(I,'salt&pepper');%H1為3*3鄰域平滑的卷積核，H2為5*5鄰域平滑的卷積核H1=1/8*[111101111];H2=1/24*[1111111111110111111111111];J1=imfilter(J,H1);%高斯白噪聲的3*3鄰域平滑J2=imfilter

4、(J,H2);%高斯白噪聲的5*5鄰域平滑K1=imfilter(K,H1);%椒鹽噪聲的3*3鄰域平滑K2=imfilter(K,H2);%椒鹽噪聲的5*5鄰域平滑J3=zhongzhi(J,3);%高斯白噪聲的3*3中值平滑J4=zhongzhi(J,5);%高斯白噪聲的5*5中值平滑K3=zhongzhi(K,3);%椒鹽噪聲的3*3中值平滑K4=zhongzhi(K,5);%椒鹽噪聲的5*5中值平滑subplot(131);imshow(J);title('高斯白噪聲');subplot(132);imshow(J1);title('高斯白噪聲的3*3鄰域平

5、滑');subplot(133);imshow(J2);title('高斯白噪聲的5*5鄰域平滑');figuresubplot(131);imshow(J);title('高斯白噪聲');subplot(132);imshow(J3);title('高斯白噪聲的3*3中值平滑');subplot(133);imshow(J4);title('高斯白噪聲的5*5中值平滑');figuresubplot(131);imshow(K);title('椒鹽噪聲');subplot(132);imshow(K1);title('椒鹽噪聲的3*3鄰域平滑');subplot(

6、133);imshow(K2);title('椒鹽噪聲的5*5鄰域平滑');figuresubplot(131);imshow(K);title('椒鹽噪聲');subplot(132);imshow(K3);title('椒鹽噪聲的3*3中值平滑');subplot(133);imshow(K4);title('椒鹽噪聲的5*5中值平滑');中值平滑子函數(shù)zhongzhi()如下：Page12of12數(shù)字圖像第四講作業(yè)functionJ=zhongzhi(I,k)[m,n]=size(I);ifk==3%3*3的中值平滑I0=zeros(m+2,n+2);fori

7、=2:m+1forj=2:n+1I0(i,j)=I(i-1,j-1);%將到操作的圖像矩陣I外圍不上0endendfori=2:m+1forj=2:n+1a=I0(i-1:i+1,j-1:j+1);b=a(1:9);%將3*3的矩陣化成1*9的矩陣，便于median操作J(i-1,j-1)=median(b);%取中值，保存為平滑后矩陣J的i-1行、j-1列endendelsek==5%5*5的中值平滑I0=zeros(m+4,n+4);fori=3:m+2forj=3:n+2I0(i,j)=I(i-2,j-2);endendfori=3:m+2forj=3:n