MATLAB 模板匹配

这里只是简单的实现一下，这个和直方图反向投影很像。

　　书上说了，“优势是理论上的，很难发展处一种满意的模板匹配技术”。

　　代码：

 1 clear all;
 2 close all;
 3 clc;
 4 
 5 img=imread('lena.jpg');
 6 imshow(img);
 7 img=double(img);
 8 mask=double(imcrop());
 9 [m n]=size(img);
10 [H W]=size(mask);
11 
12 if mod(H,2)~=1
13     H=H+1;
14 end
15 if mod(W,2)~=1
16     W=W+1;
17 end
18 mask=imresize(mask,[H W]);
19 
20 HH=floor(H/2);
21 WW=floor(W/2);
22 
23 imgn=zeros(m+2*HH+1,n+2*WW+1);
24 imgn(HH+1:m+HH,WW+1:n+WW)=img;
25 imgn(1:HH,WW+1:n+WW)=img(1:HH,1:n); 
26 imgn(1:m+HH,n+WW+1:n+2*WW+1)=imgn(1:m+HH,n:n+WW);
27 imgn(m+HH+1:m+2*HH+1,WW+1:n+2*WW+1)=imgn(m:m+HH,WW+1:n+2*WW+1);
28 imgn(1:m+2*HH+1,1:WW)=imgn(1:m+2*HH+1,WW+1:2*WW);
29 
30 re=imgn;
31 for i=HH+1:m+HH
32     for j=WW+1:n+WW
33         tmp=imgn(i-HH:i+HH,j-WW:j+WW);
34         re(i,j)=sum(sum((tmp-mask).^2));       %最小平方差 
35     end
36 end
37 
38 figure;
39 re=mat2gray(re(HH+1:m+HH,WW+1:n+WW));
40 imshow(1-re);