网上下的SUSAN程序总是运行不起来，需要修改，求大家帮助...

网上下了很多的susan程序，但是发现都是同一个，程序直接运行就开始报错，我对matlab不是很了解，所以很多小白错误都改不出来，现在希望懂的朋友能帮助下，帮我把程序改的可以正常运行，程序如下：
function image_out = susan(im,threshold)
close all
clc
im = imread('rice.png');
% check to see if the image is a color image...
d = length(size(im));
if d==3
    image=double(rgb2gray(im));
elseif d==2
    image=double(im);
end
% mask for selecting the pixels within the circular region (37 pixels, as
% used in the SUSAN algorithm
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]);  

% the output image indicating found edges
R=zeros(size(image));

% define the USAN area
nmax = 3*37/4;
% padding the image
[a b]=size(image);
new=zeros(a+7,b+7);
[c d]=size(new);
new(4:c-4,4:d-4)=image;
  
for i=4:c-4
   
    for j=4:d-4
        
        current_image = new(i-3:i+3,j-3:j+3);
        current_masked_image = mask.*current_image;
   
%   Uncomment here to implement binary thresholding
%         current_masked_image(find(abs(current_masked_image-current_masked_image(4,4))>threshold))=0;         
%         current_masked_image(find(abs(current_masked_image-current_masked_image(4,4))<=threshold))=1;

%   This thresholding is more stable
                              
        current_thresholded = susan_threshold(current_masked_image,threshold);
        g=sum(current_thresholded(:));
        
        if nmax             R(i,j) = g-nmax;
        else
            R(i,j) = 0;
        end
    end
end
image_out=R(4:c-4,4:d-4);
程序首先缺少个图像输入，我加了一句：im = imread('rice.png');
然后是提示阈值threshold没有设定，我就直接有加了一句   threshold=27.7不知道这样对不对，接下来的报错是：
??? Undefined command/function 'susan_threshold'.
Error in ==> susanedge at 91
        current_thresholded = susan_threshold(current_masked_image,threshold);
这个我就不知道怎么改了，是说函数定义不对么？应该怎么改，求大家帮忙！！

0