pipei.m

832658matlab编的特征匹配（控制点）,做这方面的朋友不妨参考。

orthophoto=imread('westconcordorthophoto.png');
subplot(1,2,1),imshow(orthophoto)
unregistered=imread('westconcordaerial.png');
subplot(1,2,2),imshow(unregistered)
cpselect(unregistered(:,:,1),orthophoto)
input_points=[120.7086,93.9772;
             319.2222,78.9202;
             127.9838,291.6312;
             352.0729,281.1445]
             input_points_corr=cpcorr(input_points,base_points,unregistered(:,:,1),orthophoto)
             input_points_corr=[120.7086,93.9772;
             319.2222,78.9202;
             127.1046,289.8935;
             352.0729,281.1445]
             mytform=cp2tform(input_points,base_points,'projective');
             registered=imtransform(unregistered,mytform)