NCC匹配算法

function  [correlate_peakimg peak_coordination] = Matchimg（ imgimg0）
%UNtitleD3 Summary of this function goes here
%   Detailed explanation goes here

% img : origin image
% img0: template image
% 

% tic;

% img=imread（inputargs）;
% img0=imread（inputargs0）;
% %Convert RGB image to gray image
% 
% if size（img3） > 1
%         img = rgb2gray（img）;
% end
% if size（img03） > 1
%         img0 = rgb2gray（img0）;
% end

% % t1=isgray（img）;
% t2=isgray（img0）;
% if ~t1
%    img=rgb2gray（img）;
% end
% if ~t2
%     img0=rgb2gray（img0）;
% end

[M N]=size（img）;
[M0 N0]=size（img0）;
new=rand（M-M0N-N0）;
% to calculate the R‘s average
img0=im2double（img0）;
i=0;
j=0;
Ssum=0;
for j=1:M0
    for i=1:N0
        SGrayVlue=img0（ji）;
        Ssum=Ssum+SGrayVlue;
    end
end
Saverage=Ssum/（N0*M0）;
% --**
% NCCO
img=im2double（img）;
MAX=0;
u=0;
v=0;
for j=1:M-M0
    for i=1:N-N0
        cropimg=imcrop（img[i j N0-1 M0-1]）;
        % to calculate the S‘s average
        cropimg=im2double（cropimg）;
        Rsum=0;
        for j1=1:M0
            for i1=1:N0
                RGrayVlue=cropimg（j1i1）;
                Rsum=Rsum+RGrayVlue;