基于角点检测的图像匹配

基于harris角点算法的图像匹配，使用matlab语言实现

%-----------------------------------------------------------------
%                  归一化相关算法
%-----------------------------------------------------------------
%%
clc;clear all ;
[FileName2PathName2] = uigetfile（{‘*.bmp‘‘All Image Files‘;...
          ‘*.*‘‘All Files‘ }）;
M= imread（FileName2）;
[FileNamePathName] = uigetfile（{‘*.bmp‘‘All Image Files‘;...
          ‘*.*‘‘All Files‘ }）;
T=imread（FileName）;
%或者这样
%T = imread（‘T.bmp‘）;                 %打开模板图像，存储在矩阵T
%M = imread（‘M.bmp‘）;                 %打开原图像，存储在矩阵M
[AB]=size（M）;
%%
%-----------------------开始计时-----------------------------          
tic                                            %开始计时        
c = normxcorr2（T（::1）M（::1））;             %归一化相关
%求相关系数
%原图像M必须要比模板图像T大，返回值在-1.0与1.0之间

[max_c imax] = max（abs（c（:）））;    %求归一化相关系数的最大值点
%max_c列向量c（:）的最大值 imax为列向量c（:）最大值的下标

[ypeak xpeak] = ind2sub（size（c）imax（1））;    
%imax（1）为矩阵c列向量化后元素的下标，ypeak为元素在矩阵c的行下标 xpeak为列下标

offset = [（xpeak-size（T2））;（ypeak-size（T1））];
%size（T2）为T的第二维的大小即T的列下标的个数
%确定模板图像的左上角点在源图像的坐标

X = offset（1）+1;       %图像匹配的X坐标
Y = offset（2）+1;       %图像匹配的Y坐标
%%
figure（‘Numbertitle‘‘off‘‘Name‘‘归一化算法实现图像匹配‘）;
subplot（221）;imshow（M）;title（‘原图‘）      %显示原图

%%
subplot（222）;imshow（T）;title（‘模板图像‘）  %显示模板图像

h = waitbar（0‘1‘‘Name‘‘进度‘）;
set（findobj（h‘type‘‘patch‘）‘edgecolor‘‘w‘‘facecolor‘‘g‘）;

    for i=1:100
    s=sprintf（‘正在运行:%d‘ceil（i/10））;
    waitbar（i/100h[s ‘%‘]）;
    end
    waitbar（1h‘完成‘）;
    pause（0.7）;
    close（h）;
%%
subplot（223）;
i=1:A;j=1:B;abc=double（M）;dd=abc（ij）;mesh（dd）;
%surf（dd）;shading flat;
title（‘图像彩色分布图‘）;
subplot（224）imshow（M）;      %显示图像
title（‘匹配的结果‘）
toc;             %结束即时，显示时间
%-----------------------结束计时----------------------------- 
%%
                %用红色矩形框出目标物体
                %rectangle（‘Position‘[坐标X坐标Y矩形宽度矩形高度]‘EdgeColor‘‘r‘）; 
rectangle（‘Position‘[XYsize（T2）size（T1）]‘EdgeColor‘‘r‘）; 
%%
                      

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-09-06 20:23  毕设最终稿\
     文件      401153  2012-03-18 16:44  毕设最终稿\1.png
     文件       38826  2013-06-10 19:53  毕设最终稿\123.jpg
     文件      436630  2012-03-18 16:45  毕设最终稿\2.png
     文件        1306  1999-11-12 13:19  毕设最终稿\butterfly.GIF
     文件       70108  2013-06-10 19:53  毕设最终稿\CIMG0878.jpg
     文件        2292  2013-06-12 21:29  毕设最终稿\Correlation.m
     文件        2008  2013-06-08 09:36  毕设最终稿\Correlation.m.bak
     文件        8828  2013-06-10 15:27  毕设最终稿\descriptFPoints.m
     文件       98245  2013-06-10 19:53  毕设最终稿\DSC_0091.jpg
     文件       45365  2013-06-10 19:54  毕设最终稿\DSC_0209.jpg
     文件       67386  2013-06-10 19:54  毕设最终稿\DSC_0563.jpg
     文件        3225  2009-09-17 13:18  毕设最终稿\Findlclmxm.m
     文件      114380  2013-05-28 10:53  毕设最终稿\fruit.jpg
     文件      189519  2013-05-28 11:04  毕设最终稿\fruit_duibidu.jpg
     文件       88474  2013-06-10 21:06  毕设最终稿\fruit_duishu.jpg
     文件      143186  2013-06-10 21:18  毕设最终稿\fruit_gaosibai.jpg
     文件      125759  2013-06-10 20:49  毕设最终稿\fruit_xie.jpg
     文件       50131  2013-06-10 20:56  毕设最终稿\fruit_xie2.jpg
     文件        7964  2013-06-08 19:46  毕设最终稿\GUI_datatransmission.fig
     文件       13057  2013-06-13 09:21  毕设最终稿\GUI_datatransmission.m
     文件        7473  2013-06-23 23:11  毕设最终稿\harrislpls.m
     文件       14060  2003-02-17 17:45  毕设最终稿\image068.JPG
     文件       13579  2003-02-17 17:45  毕设最终稿\image069.JPG
     文件      578897  2013-05-22 19:48  毕设最终稿\img1.png
     文件      574557  2013-05-22 19:48  毕设最终稿\img2.png
     文件      464954  2006-04-05 06:45  毕设最终稿\img3.jpg
     文件      467808  2006-04-05 06:45  毕设最终稿\img5.jpg
     文件        6592  2010-10-11 21:12  毕设最终稿\immatch.m
     文件        3319  2011-12-18 14:11  毕设最终稿\invmoments.m
     文件        1382  2009-08-18 03:08  毕设最终稿\kNearestNeighbors.m
............此处省略20个文件信息