function?image_out=susan（imthresthold）
d=length（size（im））;
if?d==3
????image=double（rgb2gray（im））;
elseif?d==2
????image=double（im）;
end
mask=（[0?0?1?1?1?0?0;0?1?1?1?1?1?0;1?1?1?1?1?1?1;1?1?1?1?1?1?1;1?1?1?1?1?1?1;0?1?1?1?1?1?0;?0?0?1?1?1?0?0]）;
R=zeros（size（image））;
nmax=3*37/4;
[a?b]=size（image）;
new=zeros（a+7b+7）;
[cd]=size（new）;
new（4:c-44:d-4）=image;
for?i=4:c-4
???for?j=4:d-4
????????Current_image=new（i-3:i+3j-3:j+3）;
????????Current_masked_image=mask.*Current_image;???????
????????temp1?=?（Current_masked_image-Current_masked_image（44））/thresthold;
????????temp2?=?temp1.^6;
????????Current_thresholded?=?exp（-1*temp2）;
????????g=sum（Current_thresholded（:））;
????????if??nmax???????????R（ij）=0;
???????else
???????????R（ij）=nmax-g;
???????end
???end
end
image_out=R（4:c-44:d-4）;

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

?????文件????????831??2009-03-20?10:48??susan\susan.m

?????目錄??????????0??2009-03-20?10:55??susan

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

??????????????????831????????????????????2