这里只是简单的实现一下，这个和直方图反向投影很像。
书上说了，“优势是理论上的，很难发展处一种满意的模板匹配技术”。
代码：
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);
