matlab.rar_Normalization_damage detection_normalization code_归一化

%----------------------------------------------------------------- % 归一化相关算法 天外游侠 %----------------------------------------------------------------- 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');