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