數字圖像處理——均值濾波

數字圖像處理——均值濾波

%讀取圖像，並簡單處理
Image= double(imread('2-HistogramEqualization.jpg'));
IMAGE1= Image(:,:,1);
[lines,samples,bands] = size(Image);
IMG1 = reshape(IMAGE1,1,[]);
MAX = max(IMG1);
MIN = min(IMG1);

%均值濾波
Meanfilter = [1/8 1/8 1/8;1/8 0 1/8;1/8 1/8 1/8];%定義模板
IM4 = IMAGE1;
for i= 2:lines -1
    for j = 2:samples - 1
        IM4(i,j) = round(sum(sum(Meanfilter.*[IMAGE1(i-1,j-1) IMAGE1(i,j-1) IMAGE1(i+1,j-1);
                                    IMAGE1(i-1,j) IMAGE1(i,j) IMAGE1(i+1,j);
                                    IMAGE1(i-1,j+1) IMAGE1(i,j+1) IMAGE1(i+1,j+1)])));
    end
end
figure(4)
subplot(2,2,1);
imshow(uint8(IMAGE1));
title('原始圖像');
subplot(2,2,3);
hist(IMG1,MAX);
title('原始圖像直方圖');
subplot(2,2,2);
imshow(uint8(IM4));
title('均值濾波後的圖像');
subplot(2,2,4);
hist(reshape(IM4,1,[]),255);
title('均值濾波後的直方圖');