亚像素边缘提取程序

通过求高斯梯度的方法进行亚像素级别的图像边缘提取，效果不错，附有运行结果图。

% Function “d2dgauss.m“: 
% This function returns a 2D edge detector （first order derivative 
% of 2D Gaussian function） with size n1*n2; theta is the angle that 
% the detector rotated counter clockwise; and sigma1 and sigma2 are the 
% standard deviation of the gaussian functions. 
function h = d2dgauss（n1sigma1n2sigma2theta） 
r=[cos（theta） -sin（theta）; 
sin（theta） cos（theta）]; 
for i = 1 : n2 
    for j = 1 : n1 
        u = r * [j-（n1）/2 i-（n2）/2]‘; 
        h（ij） = gauss（u（1）sigma1）*dgauss（u（2）sigma2）; 
    end 
end 
h = h / sqrt（sum（sum（abs（h）.*abs（h））））; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-09-17 18:58  gauss-lena\
     文件     1178678  2011-09-19 15:49  gauss-lena\4-lena.bmp
     文件         591  2011-09-17 17:29  gauss-lena\d2dgauss.m
     文件         129  2011-09-17 17:29  gauss-lena\dgauss.m
     文件          98  2011-09-17 17:29  gauss-lena\gauss.m
     文件        2040  2011-09-19 15:37  gauss-lena\main-gauss-lena.m