%-----------------------------------------------------------------
% 归一化相关算法 天外游侠
%-----------------------------------------------------------------
clear all ;
T = imread('T.bmp'); %打开模板图像，存储在矩阵T
M = imread('M.bmp'); %打开原图像，存储在矩阵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(T,2));(ypeak-size(T,1))];
%size(T,2)为T的第二维的大小即T的列下标的个数
%确定模板图像的左上角点在源图像的坐标
X = offset(1)+1; %图像匹配的X坐标
Y = offset(2)+1; %图像匹配的Y坐标
time = toc %结束即时，显示时间
%-----------------------结束计时-----------------------------
figure,imshow(M); %显示原图
title('Original Image')
figure,imshow(T); %显示模板图像
title('Template Image')
figure,imshow(M); %显示图像
title('Tracking Result')
%用红色矩形框出目标物体
%rectangle('Position',[坐标X,坐标Y,矩形宽度,矩形高度],'EdgeColor','r');
rectangle('Position',[X,Y,size(T,2),size(T,1)],'EdgeColor','r');
