中值滤波

摘 要：数字图像中，噪声主要来源于图像的获取和/或传输过程。当仅存在加性噪声的情况下，可以选择空间滤波方法。本实验主要编写了一个给图像中添加椒盐噪声的程序，程序的输入参数为两个噪声分量的概率值；除此之外还编写了一个3*3的中值滤波函数，对经过椒盐噪声污染的图像进行滤波。

---

function R=saltandpepper(ima,pa,pb)%功能：给图像添加椒盐噪声%输入参数：原图像ima,两个噪声分量的概率值pa，pb%输出参数：添加噪声后的图像    [m,n]=size(ima);%检查pa+pb是否小于1    if (pa+pb)>1        error('The sum Pa+Pb must not exceed 1.')    end    R=ima;    X=rand(m,n);    c=find(X<=pa);    R(c)=0;    c=find(X>pa&X<=(pa+pb));    R(c)=1;%------------------------------------------function noise    ima=imread('Fig5.07(a).jpg');    %添加椒盐噪声    imt=saltandpepper(ima,0.2,0.2);    imshow(ima);title('原图像');    figure;    imshow(imt);title('添加Pa=Pb=0.2的椒盐噪声后的图像');    %中值滤波    %一次中值滤波    [r,c]=size(imt);    imt=double(imt);    imta=zeros(r,c);    for i=2:r-1        for j=2:c-1            imta(i,j)=medain(imt,i,j);        end    end    imta=uint8(imta);    figure;imshow(imta);title('一次中值滤波之后的图像');    %二次中值滤波    imtb=zeros(r,c);    for i=2:r-1        for j=2:c-1            imtb(i,j)=medain(imta,i,j);        end    end    imtb=uint8(imtb);    figure;imshow(imtb);title('二次中值滤波之后的图像');    %三次中值滤波    imtc=zeros(r,c);    for i=2:r-1        for j=2:c-1            imtc(i,j)=medain(imtb,i,j);        end    end    imtc=uint8(imtc);    figure;imshow(imtc);title('三次中值滤波之后的图像');%---------------------------------------------function md=medain(ima,I,J)%中值滤波器%将一个3*3的矩阵转换为一行数组    b=zeros(1,9);    k=1;    for i=I-1:I+1        for j=J-1:J+1            b(1,k)=ima(i,j);            k=k+1;        end    end    %按从小到大排序    for i=1:8        for j=1:8            if (b(j)>b(j+1))                t=b(j);                b(j)=b(j+1);                b(j+1)=t;            end        end    end    %取中值返回    md=b(5);