SUSAN算法的MATLAB实现

通过MATLAB验证的SUSAN算法实现

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