f=imread('lena256.bmp')        %讀原圖
f=im2double(f);   %改變數(shù)值類型（unit-double)
imshow(f)            %顯示原來(lái)圖象
L=length(f);  %矩陣長(zhǎng)度
a=[-1,-2,-1;0,0,0;1,2,1];%水平算子
b=[-1,0,1;-2,0,2;-1,0,1];%垂直算子
cc=conv2(a,f);           %水平算子卷積圖象
dd=conv2(b,f);           %垂直算子卷積圖象
for i=1:1:L,             %截去多于象素
    for j=1:1:L
        c(i,j)=cc(i,j);
        d(i,j)=dd(i,j);
    end
end
for i=1:1:L,       %求得c,d矩陣每個(gè)值的平方
    for j=1:1:L
        c(i,j)=c(i,j)*c(i,j);
        d(i,j)=d(i,j)*d(i,j);
    end
end
t=c+d;
for i=1:1:L,
    for j=1:1:L
        t(i,j)=sqrt(t(i,j)) ; %開(kāi)方
    end
end
n=input('輸入threshold=');    %輸入閥值
for i=1:1:L,        %大于閥值1 小于閥值0
    for j=1:1:L
        if t(i,j)>n
            t(i,j)=1;
        else t(i,j)=0;
        end
    end
end
imshow(t)%顯示檢測(cè)后圖象