Fast/Robust Template Matching

Template Matching with SSD Block Matching and Normalized CC, (color) 2D and 3D

Description

TEMPLATE_MATCHING is a CPU efficient function which calculates matching score images between template and (color) 2D image or 3D image volume.

It calculates:
- The sum of squared difference (SSD Block Matching), robust template matching.
- The normalized cross correlation (NCC), independent of illumination, only dependent on texture

The user can combine the two images, to get template matching which works robust with his application.

Both matching methods are implemented using FFT based correlation.

Try the examples.

Please report: bugs, successes and other comment.

• MATLAB release: MATLAB 7.8 (R2009a)

[Masoud]

Perfect. It is fast and robust.
Thank you

[Tony]

The function 'normxcorr2' available in matlab is faster than the code you submitted. Please profile it and see for yourself. But thanks for sharing the code.

[Dirk-Jan Kroon]

*Tony
The submitted function supports 3D, SSD and color which is not in the function 'normxcorr2'. Also normxcorr2 is only a little bit faster when the template is only a few pixels in size...

[Steffen]

I tried only the 3D normalized cross-correlation. It works very well. Thanks for sharing.

[Shahbaba2011]

Thanx for sharing your code. How can I detect same template apping more than once in image.

[Marcos Belmonte]

Hi, thanks for sharing your code, it's great. I have a 3D image with a sphere, and another sphere as a template. To locate the maximum cross correlation I write:
[I_SSD,I_NCC]=template_matching(T,I);
[x,y,z]=ind2sub(size(I_NCC),find(I_NCC==max(I_NCC(:))));
daspect([1 1 1]);
figure,isosurface(I,0); hold on; plot3(y,x,z,'r*'); title('MAX NCCC3D')

Is it correct?
I don't know if the real coordinates are displayed well because I don't understant very well the normalized cross-correlation in 3D.

Thanks very much.
Regards

[Marcos Belmonte]

Hi, I have a problem with the NCCC three-dimensional, with the code I have written above, it plots a * in all the slices not only in the slices that match the template.
What should I do? Thanks in advance.

[Marcos Belmonte]

The problem is that z is logical and I want it to be the number of the slice.

[Carlos]

Hi, You used any paper to develop this code? Can you give me the name of this paper?

Thanks a lot