tongjipaixu.m

来自「这是一些滤波的小程序,用MATLAB编的,没有用到MATLAB自身的一些函数」· M 代码 · 共 27 行

clear all;
map=imread('Fig3.37(a).jpg');
[N,M]=size(map);
B=zeros(9);
for n=2:N-1
    for m=2:M-1        
        t=1;
        for x=n-1:n+1
            for y=m-1:m+1
               B(t)=map(x,y);
               t=t+1;
            end
        end
        for t=1:8
            for p=1:9-t
                if B(p)<B(p+1)
                    q=B(p);
                    B(p)=B(p+1);
                    B(p+1)=q;
                end
            end
        end
     map(n,m)=B(5);
    end
end
imagesc(map);
colormap(gray)