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)~=113 H=H+1;14 end

15 if mod(W,2)~=116 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+HH32 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 end36 end37

38 figure;

39 re=mat2gray(re(HH+1:m+HH,WW+1:n+WW));40 imshow(1-re);