%程序L5_2.m: 平滑濾波器%
%Image: 原圖像%
%Image_noisy: 受污染的圖像%
%Image_low: 執行低通濾波器后的圖像%  
%Image_med: 執行中值濾波器后的圖像%
Image=imread('L5_2.bmp');                   %讀取圖像
%將胡椒鹽似的噪聲加入原圖像中，其中0.02是噪聲密度%
Image_noisy= imnoise (Image ,'salt & pepper',0.06); 
Image2_noisy=double(Image_noisy)/255;       %變為 double格式
h=[1/9 1/9 1/9;1/9 1/9 1/9 ; 1/9 1/9 1/9];  %濾波器的子模板
Image_low=filter2(h, Image2_noisy);         %執行低通濾波
Image_med=medfilt2(Image_noisy,[3 3]);      %中值濾波采用3×3矩陣
%顯示圖像%
imshow(Image)                     
figure,imshow(Image_noisy)
figure,imshow(Image_low)
figure,imshow(Image_med)