GEN_SUSAN - Generalized SUSAN 2-D filtering SUSAN filtering with filter kernel W scaled with generalized Gaussian of intensity difference. Different prefiltering functions can be selected as well as width and exponent of the intensity scaling. GEN_SUSAN can produce filtering with caracteristics similar to wiener2 and medfilt2 and "everything inbetween"
MATLAB 6.5 (R13)
Gen_susan can make use of wiener2 and medfilt2 (of the image processing toolbox) if available, but they are not required
I believe there is a bug in the line 81, should be
I_internal = I([ones(1,ceil((regsize(1)-1)/2)) 1:end size(I,1)*ones(1,(regsize(1)-1)/2)],...
[ones(1,ceil((regsize(2)-1)/2)) 1:end size(I,2)*ones(1,(regsize(2)-1)/2)]);