图像匹准

%图像匹准
unregistered=imread('westconcordaerial.png');
load westconcordpoints;
t_concord=cp2tform(input_points,base_points,'projective');
info=imfinfo('westconcordorthophoto.png');
registered=imtransform(unregistered,t_concord,...
        'XData',[1 info.Width],'YData',[1 info.Height]);
subplot(3,1,1);imshow(unregistered)
subplot(3,1,2);imshow('westconcordorthophoto.png')
subplot(3,1,3);imshow(registered)