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?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-12-05?14:28??ROA綆楁硶浠ｇ爜\
?????文件????????2873??2011-07-05?15:59??ROA綆楁硶浠ｇ爜\roa.asv
?????文件????????2854??2011-07-07?09:29??ROA綆楁硶浠ｇ爜\roa.m
?????文件???????66614??2011-07-05?15:58??ROA綆楁硶浠ｇ爜\seashore.bmp
?????文件???????66616??2011-07-07?09:29??ROA綆楁硶浠ｇ爜\棣栧皵L娉㈡.bmp
?????文件???????66614??2011-06-29?03:14??ROA綆楁硶浠ｇ爜\棣栧皵L娉㈡9綾沖垎杈ㄧ巼.bmp
?????目錄???????????0??2011-12-05?14:28??ROEWA綆楀瓙\
?????文件????????2165??2011-07-11?14:24??ROEWA綆楀瓙\roewa.m
?????文件???????66614??2011-07-05?15:49??ROEWA綆楀瓙\seashore.bmp
?????文件??????263224??2011-07-11?14:24??ROEWA綆楀瓙\鍙板寳.bmp
?????文件???????66616??2011-07-07?09:30??ROEWA綆楀瓙\棣栧皵L娉㈡.bmp
?????文件???????66614??2011-06-23?04:40??ROEWA綆楀瓙\棣栧皵L娉㈡9綾沖垎杈ㄧ巼.bmp