function [J]=medi(I,W)
[m,n]=size(I);%原來圖像的高度和寬度
[mm,nn]=size(W);%模板的高度和寬度
J=zeros(m,n);%處理后的圖像
%以下對(duì)圖像中每個(gè)點(diǎn)進(jìn)行處理
for i=1:m
    for j=1:n
        %如果是邊緣點(diǎn)，則不處理，賦予原來圖像的值
        if i<=(mm-1)/2||i>=m-(mm-1)/2+1||j<=(nn-1)/2||j>=n-(nn-1)/2+1
            J(i,j)=I(i,j);
        else
            %對(duì)于不是邊緣點(diǎn)，則按模板的大小進(jìn)行計(jì)算，把模板對(duì)應(yīng)原圖像上的點(diǎn)放到一個(gè)數(shù)組里
            k=0;
            temp=zeros(1,mm*nn);
            for a=-(mm-1)/2:(mm-1)/2
                for b=-(nn-1)/2:(nn-1)/2
                    temp(k+1)=I(i+a,j+b);
                    k=k+1;
                end
            end
            %對(duì)于模板上對(duì)應(yīng)原圖像上的點(diǎn)進(jìn)行排序，采用冒泡排序
            for s=1:k-2
               for l=s+1:k
                   if temp(s)>temp(l)
                       tem=temp(s);
                       temp(s)=temp(l);
                       temp(l)=tem;
                   end
               end
            end
            J(i,j)=temp((k+1)/2);%把排序后的值賦給處理后的矩陣
        end
    end
end