基于harris算法的角点检测matlab原代码

这是一个用matlab编写的harris角点检测源代码，实验结果很好，是图像配准中的前期工作。

function [rowcolmax_local] = findLocalMaximum（valradius）
    % Determine the local maximum of a given value
    %
    % Author :: Vincent Garcia
    % Date   :: 09/02/2007
    %
    % INPUT
    % =====
    % val    : the NxM matrix containing values
    % radius : the radius of the neighborhood
    %
    % OUTPUT
    % ======
    % row       : the row position of the local maxima
    % col       : the column position of the local maxima
    % max_local	: the NxM matrix containing values of val on unique local maximum
    %
    % EXAMPLE
    % =======
    % [lcm] = findLocalMaximum（imgradius）;
    



    % FIND LOCAL MAXIMA BY DILATION （FAST） /!\ NON UNIQUE /!\
    % mask = fspecial（‘disk‘radius）>0;
    % val2 = imdilate（valmask）;
    % index = val==val2;
    % [rowcol] = find（index==1）;
    % max_local = zeros（size（val））;
    % max_local（index） = val（index）;


    % FIND UNIQUE LOCAL MAXIMA USING FILTERING （FAST）
    mask  = fspecial（‘disk‘radius）>0;
    nb    = sum（mask（:））;
    highest          = ordfilt2（val nb mask）;
    second_highest   = ordfilt2（val nb-1 mask）;
    index            = highest==val & highest~=second_highest;
    max_local        = zeros（size（val））;
    max_local（index） = val（index）;
    [rowcol]        = find（index==1）;


    % FIND UNIQUE LOCAL MAXIMA （FAST）
    % val_height  = size（val1）;
    % val_width   = size（val2）;
    % max_local   = zeros（val_heightval_width）;
    % val_enlarge = zeros（val_height+2*radiusval_width+2*radius）;
    % val_mask    = zeros（val_height+2*radiusval_width+2*radius）;
    % val_enlarge（ （1:val_height）+radius  （1:val_width）+radius ） = val;
    % val_mask（    （1:val_height）+radius  （1:val_width）+radius ） = 1;
    % mask  = fspecial（‘disk‘radius）>0;
    % row = zeros（val_height*val_width1）;
    % col = zeros（val_height*val_width1）;
    % index = 0;
    % for l = 1:val_height
    %     for c = 1:val_width
    %         val_ref = val（lc）;
    %         neigh_val  = val_enlarge（l:l+2*radiusc:c+2*radius）;
    %         neigh_mask = val_mask（   l:l+2*radiusc:c+2*radius）.*mask;
    %         neigh_sort = sort（neigh_val（neigh_mask==1））;
    %         if val_ref==neigh_sort（end） && val_ref>neigh_sort（end-1）
    %             index          = index+1;
    %             row（index1）   = l;
    %             col（index1）   = c;
    %             max_local（lc） = val_ref;
    %         end
    %     end
    % end
    % row（index+1:end:） = [];
    % col（index+1:end:） = [];


end

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

     文件     304182  2008-08-25 16:45  keypointExtraction\keypointExtraction\1.bmp

     文件     102454  2009-04-15 09:56  keypointExtraction\keypointExtraction\1G.BMP

     文件      66182  2008-08-25 16:45  keypointExtraction\keypointExtraction\2.bmp

     文件      77878  2009-04-15 09:56  keypointExtraction\keypointExtraction\2G.BMP

     文件     170474  2007-08-01 13:21  keypointExtraction\keypointExtraction\door.jpg

     文件       2563  2007-12-07 10:25  keypointExtraction\keypointExtraction\findLocalMaximum.m

     文件       1133  2007-12-07 10:28  keypointExtraction\keypointExtraction\kp_gilles.m

     文件       2032  2009-04-15 10:36  keypointExtraction\keypointExtraction\kp_harris.m

     文件       2032  2009-04-15 10:36  keypointExtraction\keypointExtraction\kp_harris1.m

     文件       3512  2007-12-07 10:38  keypointExtraction\keypointExtraction\kp_harrislaplace.m

     文件       1962  2007-12-07 10:10  keypointExtraction\keypointExtraction\kp_log.m

     文件        631  2007-12-07 10:19  keypointExtraction\keypointExtraction\kp_susan.m

     文件      18357  2007-08-03 12:53  keypointExtraction\keypointExtraction\patrol.jpg

     文件       1969  2007-12-06 09:45  keypointExtraction\keypointExtraction\README.txt

     文件     139536  2007-08-01 07:15  keypointExtraction\keypointExtraction\sunflower.jpg

     文件      40960  2006-09-23 16:22  keypointExtraction\keypointExtraction\susan.dll

     文件       2113  2006-09-23 17:16  keypointExtraction\keypointExtraction\susan.m

     文件       1536  2009-04-15 16:11  keypointExtraction\keypointExtraction\test.m

     文件      65551  2006-12-14 15:46  keypointExtraction\keypointExtraction\test.pgm

     文件      26457  2009-04-15 09:52  keypointExtraction\keypointExtraction\测试结果\1H.jpg

     文件      26587  2009-04-15 10:05  keypointExtraction\keypointExtraction\测试结果\1s.jpg

     文件      21546  2009-04-15 09:54  keypointExtraction\keypointExtraction\测试结果\2H.jpg

     文件      19167  2009-04-15 10:02  keypointExtraction\keypointExtraction\测试结果\2S.jpg

    ..A.SH.     23040  2009-04-15 10:06  keypointExtraction\keypointExtraction\测试结果\Thumbs.db

     目录          0  2009-04-15 10:45  keypointExtraction\keypointExtraction\测试结果

     目录          0  2009-04-15 16:02  keypointExtraction\keypointExtraction

     目录          0  2009-04-15 16:02  keypointExtraction

----------- ---------  ---------- -----  ----

              1121854                    27



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