registp.m

MATLAB基于点的图像配准程序

%function [R,t]= RegistP(Ref,Flo,Pref,Pflo,N )
Imref=imread('pout.bmp');
Imflo=imrotate(Imref,45);
figure(1)
subplot(121),imshow(Imref)
subplot(122),imshow(Imflo)
N=7;
[row1,col1]=size(Imref);
[row2,col2]=size(Imflo);
%cpselsct(Imflo,Imref);
%选取控制点,自己选取对应点，例如A对应A1
A=[1	1]';B=[1 236.53]';
C=[	291.68 239.51 ]';D=[288.7  1]';
E=[	103.86 50.193]';F=[	290.19 97.894]';
G=[	3.9814 111.31 ]';A1=[169.38   3.9802]';
B1=[5.4704    169.38]';C1=[209.62   373.53]';
D1=[372.04   206.64 ]';E1=[209.62   108.29]';
F1=[	303.49   276.67]';G1=[90.407  84.447 ]';
Pref=[A B C D E F G];
Pflo=[A1 B1 C1 D1 E1 F1 G1];
H=zeros(2);
xmean=[sum(Pref(1,:))/N;sum(Pref(2,:))/N];
ymean=[sum(Pflo(1,:))/N;sum(Pflo(2,:))/N];
Pref1=[Pref(1,:)-xmean(1);Pref(2,:)-xmean(2)];
Pflo1=[Pflo(1,:)-ymean(1);Pflo(2,:)-ymean(2)];
for i=1:N   
H=H+Pflo1(:,i)*(Pref1(:,i))';
end
[U,S,V]=svd(H);
r=det(U*V);
W=[1 0;0 r];
R=V*W*U';
t=xmean-R*ymean;
%%进行双线性插值
imm=im2double(Imflo);
ii=0;
 for i =1:row1%row
    ii=ii+1; jj=0;
    for j = 1:col1%col
        jj=jj+1;
        SS=inv(R)*([i j]'-t);
         if(((SS(1)<=1)||(SS(1)>=row2))||((SS(2)<=1)||(SS(2)>=col2)))
            S(ii,jj) = 0;
        else
             S(ii,jj) = BilinearInterplate(imm,SS);
         end
    end

 end

 S0= im2uint8(S);
 figure(2)
 imshow(S0)