zuobiaozhuanhuan.txt

图像内定像中的坐标转换程序代码

fid = fopen('zuobiao.txt','r');
[X1,count] =fscanf(fid, '%f,', [16,inf]);
[X2,count]=fscanf(fid, '\n%f,', [2,inf]);
[X3,count]=fscanf(fid, '\n%f,', [2,inf]);
fclose(fid);
[m,n]=size(X1);
for i=1:2:m
    X1(i,1)=X1(i,1)+0.004;
    X1(i+1,1)=X1(i+1,1)+0.008;
 end   
for i=1:2:m
A(i,1)=1;
A(i+1,2)=1;
A(i+1,1)=0;
A(i,2)=0;
A(i+1,3)=0;
A(i,4)=0;
A(i+1,5)=0;
A(i,6)=0;
A(i+1,7)=0;
A(i,8)=0;
end
j=1;
for i=1:2:m
A(i,3)=X2(1,j);
A(i+1,4)=X2(1,j);
A(i,5)=X2(2,j);
A(i+1,6)=X2(2,j);
A(i,7)=X2(1,j)*X2(2,j);
A(i+1,8)=X2(1,j)*X2(2,j);
j=j+1;
end
t=inv(A'*A)*A'*X1;
[m1,n1]=size(X3);
for i=1:2:n1*2
B(i,1)=1;
B(i+1,2)=1;
B(i+1,1)=0;
B(i,2)=0;
B(i+1,3)=0;
B(i,4)=0;
B(i+1,5)=0;
B(i,6)=0;
B(i+1,7)=0;
B(i,8)=0;
end
k=1;
for i=1:2:n1*2
B(i,3)=X3(1,k);
B(i+1,4)=X3(1,k);
B(i,5)=X3(2,k);
B(i+1,6)=X3(2,k);
B(i,7)=X3(1,k)*X3(2,k);
B(i+1,8)=X3(1,k)*X3(2,k);
k=k+1;
end
X=B*t;
fid = fopen('zhuanhuan.txt', 'wt');
fprintf(fid, '%f\n',t);
fprintf(fid,'%f,%f\n',X);
fclose(fid);