sar-edge-detection sar 图像边缘检测算子

sar 图像边缘检测算子，代码简洁，能够运行成功。

a=imread（‘首尔L波段.bmp‘）;
% a=imread（‘seashore.bmp‘）;
X=1+double（a）/255*254;
[MN]=size（X）;
T1=0.2;
for i=1:M
    for j=1:N
        output（ij）=0;
        direction（ij）=0;
        r=[];
        if  i>4 & i4 &j            %angle 0
            u1=mean（X（i-4:i-2j））;
            u2=mean（X（i-1:i+1j））; 
            u3=mean（X（i+2:i+4j））;
            r12=1-min（u1/u2u2/u1）;
            r13=1-min（u2/u3u3/u2）;
            r=[rmin（r12r13）];
            %angle 22.5
            u1=mean（[X（i-4j-2）X（i-3j-2）X（i-2j-1）]）;
            u2=mean（[X（i-1j-1）X（ij）X（i+1j+1）]）;
            u3=mean（[X（i+2j+1）X（i+3j+2）X（i+4j+2）]）;         
            r12=1-min（u1/u2u2/u1）;
            r13=1-min（u2/u3u3/u2）;
            r=[rmin（r12r13）];
            %angle 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2854  2011-07-07 09:29  sar edge detection\roa.m

     文件       2165  2011-07-11 14:24  sar edge detection\roewa.m

     目录          0  2012-10-12 10:20  sar edge detection

----------- ---------  ---------- -----  ----

                 5019                    3