myself.m

给出一组标定点的二维坐标即可得到摄象机的内外参数．

clc
clear
x=load('pyt\M.txt');      %model,ie pattern (world axis)
X=load('pyt\Dscfl3.txt'); %image pixels
Xf=X(:,1);
Yf=X(:,2);
xw=x(:,1);
yw=x(:,2);
[M,N]=size(x);
zw=zeros(M,1);
Ncx=1;
Nfx=1;
Cx=320/2;  %center of image,ie u0 and v0
Cy=240/2;
dx=6.45e-4; %6.45um
dy=6.45e-4;
sx=1;
plot(X(:,1),X(:,2),'*'),title('X:image')
figure
plot(x(:,1),x(:,2),'.'),title('x:axis xy')
[R, T, f, k1] = Tsai(Xf, Yf, xw, yw, zw, Ncx, Nfx, dx, dy, Cx, Cy, sx)