焦点检测Harris算子 SUSAN算子 文档 代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%   Harris角点提取算法和改进的Harris角点提取算法                              % 
%   2005-11-25                                                                % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear; 
FileInfo = imfinfo（‘E:Image1MagImg1.bmp‘）; 
[ImageDatamap] = imread（‘E:Image1
4.bmp‘）; 
%  
% if（strcmp（‘truecolor‘FileInfo.ColorType） == 1） 
%    ImageData = im2uint8（rgb2gray（ImageData））; 
% elseif（strcmp（‘indexed‘FileInfo.ColorType） == 1） 
%    ImageData = im2uint8（ind2gray（ImageDatamap））;   
% end     
ori_im=rgb2gray（ImageData）; 
%%ori_im=histeq（ori_im）;
%ori_im = imread（‘arroyo-r.tiff‘）;     % 读取图像 
 
% fx = [5 0 -5;8 0 -8;5 0 -5];          % 高斯函数一阶微分，x方向（用于改进的Harris角点提取算法） 
fx = [-2 -1 0 1 2];                 % x方向梯度算子（用于Harris角点提取算法） 
Ix = filter2（fxori_im）;              % x方向滤波 
% fy = [5 8 5;0 0 0;-5 -8 -5];          % 高斯函数一阶微分，y方向（用于改进的Harris角点提取算法） 
fy = [-2;-1;0;1;2];                 % y方向梯度算子（用于Harris角点提取算法） 
Iy = filter2（fyori_im）;              % y方向滤波 
Ix2 = Ix.^2; 
Iy2 = Iy.^2; 
Ixy = Ix.*Iy; 
clear Ix; 
clear Iy; 
 
h= fspecial（‘gaussian‘[7 7]2）;      % 产生7*7的高斯窗函数，sigma=2 
 
Ix2 = filter2（hIx2）; 
Iy2 = filter2（hIy2）; 
Ixy = filter2（hIxy）; 
 
height = size（ori_im1）; 
width = size（ori_im2）; 
result = zeros（heightwidth）;         % 纪录角点位置，角点处值为1 
 
R = zeros（heightwidth）; 
 
Rmax = 0;                              % 图像中最大的R值 
for i = 1:height 
    for j = 1:width 
        M = [Ix2（ij） Ixy（ij）;Ixy（ij） Iy2（ij）];             % auto correlation matrix 
        R（ij） = det（M）-0.06*（trace（M））^2;                     % 计算R 
        if R（ij） > Rmax 
            Rmax = R（ij）; 
        end; 
    end; 
end; 
 
cnt = 0; 
for i = 2:height-1 
    for j = 2:width-1 
        % 进行非极大抑制，窗口大小3*3 
        if R（ij） > 0.01*Rmax && R（ij） > R（i-1j-1） && R（ij） > R（i-1j） && R（ij） > R（i-1j+1） && R（ij） > R（ij-1） && R（ij） > R（ij+1） && R（ij） > R（i+1j-1） && R（ij） > R（i+1j） && R（ij） > R（i+1j+1） 
            result（ij） = 1; 
            cnt = cnt+1; 
        end; 
    end; 
end; 
 
[posc posr] = find（result == 1）; 
cnt                                      % 角点个数 
imshow（ori_im）; 
hold on; 
plot（posrposc‘r+‘）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     631648  2010-06-30 16:11  角点检测2_2_CSS角点检测.pdf

     文件       2139  2000-01-22 12:00  角点检测Cornera.gif

     文件       1707  2000-01-22 12:00  角点检测Corner.gif

     文件       1745  2000-01-22 12:00  角点检测Cornerc.gif

     文件      10900  2010-07-09 10:29  角点检测Cornercorner.jpg

     文件      13922  2006-01-24 14:24  角点检测Cornercorner.m

     文件        338  2010-07-09 10:29  角点检测Cornercout.mat

     文件       1426  2000-01-22 12:00  角点检测Cornerd.gif

     文件       1759  2000-01-22 12:00  角点检测Cornere.gif

     文件       9603  2010-07-09 10:29  角点检测Corneredge.jpg

     文件       1879  2000-01-22 12:00  角点检测Cornerf.gif

     文件       1916  2000-01-22 12:00  角点检测Cornerg.gif

     文件      19024  2005-05-10 16:46  角点检测Cornergui_corner.fig

     文件       8565  2005-05-10 16:54  角点检测Cornergui_corner.m

     文件       2182  2006-01-17 17:45  角点检测Cornerh.gif

     文件      66568  2003-10-17 05:43  角点检测Cornerhouse.bmp

     文件      21406  2010-07-07 16:21  角点检测FART关键点C和MATLab程序fast-C-src-1.2.zip

     文件     174649  2010-07-07 16:21  角点检测FART关键点C和MATLab程序fast-matlab-src.zip

     文件      23782  2004-03-07 17:05  角点检测Forstner 角点检测 自编Forstner 角点检测 自编dibapi.cpp

     文件       1876  2000-02-03 07:29  角点检测Forstner 角点检测 自编Forstner 角点检测 自编dibapi.h

     文件      33304  2010-07-08 11:38  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.aps

     文件       4096  2010-07-08 15:18  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.clw

     文件       4230  2005-09-22 21:51  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.cpp

     文件       4999  2007-09-11 21:15  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.dsp

     文件        541  2005-09-22 21:51  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.dsw

     文件       1389  2005-09-22 21:51  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.h

     文件     181248  2010-07-08 15:18  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.ncb

     文件     386048  2010-07-08 15:18  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.opt

     文件       1587  2007-09-27 15:09  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.plg

     文件      14718  2007-09-11 14:54  角点检测Forstner 角点检测 自编Forstner 角点检测 自编imagetest.rc

............此处省略91个文件信息