%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?%絕對誤差和算法（SAD）??
%?clc;
%?clear;
%?
%?%原始圖
%?src=imread（‘src1.jpg‘）;??
%?[a?b?d1]=size（src）;??
%?if?d1==3??
%?????src=rgb2gray（src）;??
%?end??
%?
%?%模板
%?mask=imread（‘template1.jpg‘）;??
%?[m?n?d2]=size（mask）;??
%?if?d2==3??
%?????mask=rgb2gray（mask）;??
%?end??
%?N=n;%模板尺寸，默認模板為正方形??
%?M=a;%代搜索圖像尺寸，默認搜索圖像為正方形??
%?dst=zeros（M-NM-N）;??
%?for?i=1:M-N?????????%子圖選取，每次滑動一個像素??
%?????for?j=1:M-N??
%?????????temp=src（i:i+N-1j:j+N-1）;%當前子圖??
%?????????dst（ij）=dst（ij）+sum（sum（abs（temp-mask）））;??
%?????end??
%?end??
%?
%?abs_min=min（min（dst））;
%?[xy]=find（dst==abs_min）;%返回最小值在矩陣中的行列號??
%?figure;??
%?imshow（mask）;title（‘模板‘）;??
%?figure;
%?imshow（src）;??
%?hold?on;??
%?rectangle（‘position‘[yxN-1N-1]‘edgecolor‘‘r‘）;??
%?hold?off;title（‘搜索結果‘）;??

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%基于灰度的模板匹配算法??
%絕對誤差和算法（SAD）?
%誤差平方和算法（SSD）
%思想一致，只是相似度測度公式不一樣而已！！！
clc;
clear;

%原始圖
src=imread（‘E:\任務\圖像匹配\模板匹配\SAD\images\src1.jpg‘）;??
[m?n?d1]=size（src）;%m行n列?
if?d1==3??
????src=rgb2gray（src）;??
end??

%模板
mask=imread（‘E:\任務\圖像匹配\模板匹配\SAD\images\template1.jpg‘）;??
[M?N?d2]=size（mask）;%M行N列??
if?d2==3??
????mask=rgb2gray（mask）;??
end??
?
%存放各個子圖與模板的絕對誤差和
dst=zeros（m-M+1n-N+1）;??
for?i=1:m-M+1?????????%子圖選取，每次滑動一個像素??
????for?j=1:n-N+1??
????????temp=src（i:i+M-1j:j+N-1）;%當前子圖??
????????
????????%絕對誤差和算法（SAD）
????????dst（ij）=dst（ij）+sum（sum（abs（temp-mask）））;?
????????
????????%誤差平方和算法（SSD）
%???????dst（ij）=dst（ij）+sum（sum（（temp-mask）.^2））;%.^：矩陣中每個元素的平方??
????end??
end??

%?m=min（dst）;%找出矩陣dst中每列中的最小元素，構成行向量m
%?mm=min（m）;%進一步找出m中的最小元素mm，當然也就是矩陣dst中的最小元素
abs_min=min（min（dst））;
[rc]=find（dst==abs_min）;%返回最小值在dst（同樣也在src）中的行列號?r：行?c：列??
x=c;%x坐標對應的是列
y=r;%y坐標對應的是行
figure;??
imshow（mask）;title（‘模板‘）;??
figure;
imshow（src）;??
hold?on;??
rectangle（‘position‘[xyN-1M-1]‘edgecolor‘‘r‘）;?%[x坐標y坐標寬高]?
hold?off;title（‘匹配結果‘）;??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2324??2018-04-10?15:11??SAD\BM_SAD.m

?????文件??????50544??2017-12-25?10:54??SAD\src.jpg

?????文件???????6889??2017-12-25?11:37??SAD\src1.jpg

?????文件??????33528??2017-12-25?10:54??SAD\template.jpg

?????文件??????18182??2017-12-25?11:37??SAD\template1.jpg

?????目錄??????????0??2018-04-10?15:10??SAD

-----------?---------??----------?-----??----

???????????????111467????????????????????6