clear
clc
Background_Source=imread（‘Scene.jpg‘）;
Template_Source=imread（‘Template_1.jpg‘）;

if?size（Background_Source3）==3
????Template=rgb2gray（Background_Source）;
????Template=double（Template）;
else
????Template=double（Template_Source）;
end
if?size（Template_Source3）==3
????Background=rgb2gray（Template_Source）;
????Background=double（Background）;
else
????Background=double（Background_Source）;
end

[Bk_r?Bk_c]=size（Background）;
[Te_r?Te_c]=size（Template）;
%實現公式運算
for?i=1:（Bk_r-Te_r+1）
????for?j=1:（Bk_c-Te_c+1）
????????A=Background（i:i+Te_r-1j:j+Te_c-1）.*Template;
????????B=sum（A（:））;
????????C=Background（i:i+Te_r-1j:j+Te_c-1）.*Background（i:i+Te_r-1j:j+Te_c-1）;
????????D=sum（C（:））;
????????E=Template.*Template;
????????F=sum（E（:））;
????????Result（ij）=B/sqrt（D*F）;
????end
end

[C_maximax]=max（Result（:））;
[y?x]=ind2sub（size（Result）imax）;%函數可以用來把矩陣元素的index轉換成對應的下標，并輸出坐標
y1=Bk_r-Te_r-y+1
x1=x-1
figure;
subplot（132）;imshow（Background_Source）;
title（‘原始圖像‘）;
subplot（131）;imshow（Template_Source）;
title（‘模板圖像‘）;
subplot（133）;imshow（Background_Source）;
title（‘匹配結果‘）;
rectangle（‘Position‘[xysize（Template2）size（Template1）]‘EdgeColor‘‘r‘）;?