homework3.m

这是一个基于matlab的图像匹配算法

I1=imread('.\jjj.tif');
I2=imread('.\lena.tif');
subplot(1,2,1),imshow(I2);
subplot(1,2,2),imshow(I1);

sub_I1=imcrop(I1);

c=normxcorr2(sub_I1(:,:,1),I2(:,:,1));

figure,surf(c),shading flat;

[max_c,imax]=max(abs(c(:)));
[ypeak,xpeak]=ind2sub(size(c),imax(1));

[m,n]=size(sub_I1);
xbegin=xpeak-n+1;
ybegin=ypeak-m+1;
xend=xpeak;
yend=ypeak;

extracted_I1=I2(ybegin:yend,xbegin:xend,:);
disp(m);
disp(n);
figure,imshow(extracted_I1);
disp(xbegin);
disp(ybegin);

for i=1:m
    for j=1:n
    T(i,j)=extracted_I1(i,j)-sub_I1(i,j);
    end
end
disp(T);