BS=uint8（zeros（288352））;
B=imread（‘1-HLS.bmp‘）;
B=medfilt2（B）;
Count_Susan_plot=0;
AREA=0;
tic

for?i=5:285?????????????%??行?????????????%i=1：100是天空
????for?j=5:348????????%??列?????
????????????for?m=i-1:i+1
????????????????for?n=j-3:j+3
????????????????????if??abs（B（mn）-B（ij））<4
???????????????????????AREA=AREA+1;
????????????????????end????
????????????????end
????????????end??
????????????for?n=j-1:j+1
????????????????if??abs（B（i-3n）-B（ij））<4
????????????????????AREA=AREA+1;
????????????????end????
????????????end
????????????for?n=j-2:j+2
????????????????if??abs（B（i-2n）-B（ij））<4
????????????????????AREA=AREA+1;
????????????????end????
????????????end
????????????for?n=j-1:j+1
????????????????if??abs（B（i+3n）-B（ij））<4
????????????????????AREA=AREA+1;
????????????????end????
????????????end
????????????for?n=j-2:j+2
????????????????if??abs（B（i+2n）-B（ij））<4
????????????????????AREA=AREA+1;
????????????????end????
????????????end
????????????
????????????
????????????
?????????if?AREA<14
?????????????BS（ij）=255;?
?????????????Count_Susan_plot=Count_Susan_plot+1;
?????????end???
?????????AREA=0;
???????
????end
end

Count_Susan_plot
toc
subplot（121）;
imshow（B）;
subplot（122）;
imshow（BS）;
imwrite（BS‘image0012-1-1.bmp‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1313??2008-01-24?16:24??SUSAN37\SUSAN.m

?????文件??????12465??2007-12-17?17:32??SUSAN37\SUSAN37.m

?????目錄??????????0??2009-06-26?22:23??SUSAN37

-----------?---------??----------?-----??----

????????????????13778????????????????????3