saltandpeppermedian.asv

来自「Image enhancement in frequency domain us」· ASV 代码 · 共 18 行

clear all
clc
I=imread('Fig0427(a)(woman).tif');
J=imnoise(I,'salt & pepper',0.02);
a=double(J);
b=a;
[row col]=size(a);
for x=2:1:row-1;
    for y=2:1:col-1;
        a1=[a(x-1,y-1) a(x-1,y) a(x-1,y+1) a(x,y-1) a(x,y) a(x,y+1) a(x+1,y-1) a(x+1,y) a(x+1,y+1)];
        a2=sort(a1);
        med=a2(5);
        b(x,y)=med;
    end
end
figure(1),imshow(uint8(I)),title('Original image');
figure(2),imshow(uint8(J)),title('Original image with salt and pepper Noise');
figure(3),imshow(uint8(b)),title('Noise Reduction Median Filtering of the Original Image');