投影仪标定程序

投影仪标定程序，反投影精度在1个象素之内

img=load（‘image2d1.txt‘）;
wrd=load（‘word3d1.txt‘）;
[ab]=size（img）;%获得矩阵的大小
p=1;
for i =2: 2: b-2
   Chy（1:8p）=img（1:8i+2）-img（1:8i）;%保存y方向两个圆心象素差值的结果
   p=p+1;%循环变量
end
save（‘Chy‘）;
k=1;
for i =1: 2: b-3
  Chx（1:8k）=img（1:8i+2）-img（1:8i）;%保存x方向两个圆心象素差值的结果
  k=k+1;%循环变量
end
save（‘Chx‘）;

[cd]=size（Line）;
%[ef]=size（Line）;
t=1;
for i=1:6:c
    Lin（1:6t）=Line（i:i+51）;
    t=t+1;
end
save（‘Lin‘）;

%for i=1:6
   %Q= Lin（img（iimg（:i）））-img（:i）;
%end
for i=1 :6
    Avg（i）=mean（Lin（i1:1020））;
end
for i=1:6
    AP（1:8i）=Avg（1i）*ones（81）-img（1:82.*i）;
end

lmda=AP./Chy;

%q=1;
%for i=1:2:11
%Xa（1:8q）=img（1:8i）;
%q=q+1;
%end

%s=1;
%for i=2:2:12
%Ya（1:8s）=img（1:8i）;
%s=s+1;
%end

q=1;
for i=1:2:11
Xa（1:8q）=wrd（1:8i）;
q=q+1;
end
%Xa=XA‘;

s=1;
for i=2:2:12
Ya（1:8s）=wrd（1:8i）;
s=s+1;
end

Xp=Xa+lmda.*Chx;
Yp=Ya+lmda.*Chy;
save（‘Xp‘）;
save（‘Yp‘）;
Zp=load（‘axz.txt‘）;

%XPZ=Xp‘;
%YPZ=Yp‘;
%ZPZ=Zp‘;
one=ones（18）;
zero=one*0;
%B1=[Xp（1:）;Yp（1:）;Zp（1:）;one;zero;zero;zero;zero];
%B2=[zero;zero;zero;zero;Xp（1:）;Yp（1:）;Zp（1:）;one];
%B=[B1‘; B2‘];
A1=[Xp（:1）‘;Yp（:1）‘;Zp（:1）‘;one;zero;zero;zero;zero];
A2=[zero;zero;zero;zero;Xp（:1）‘;Yp（:1）‘;Zp（:1）‘;one];
A3=[Xp（:2）‘;Yp（:2）‘;Zp（:2）‘;one;zero;zero;zero;zero];
A4=[zero;zero;zero;zero;Xp（:2）‘;Yp（:2）‘;Zp（:2）‘;one];
A5=[Xp（:3）‘;Yp（:3）‘;Zp（:3）‘;one;zero;zero;zero;zero];
A6=[zero;zero;zero;zero;Xp（:3）‘;Yp（:3）‘;Zp（:3）‘;one];
A7=[Xp（:4）‘;Yp（:4）‘;Zp（:4）‘;one;zero;zero;zero;zero];
A8=[zero;zero;zero;zero;Xp（:4）‘;Yp（:4）‘;Zp（:4）‘;one];
A9=[Xp（:5）‘;Yp（:5）‘;Zp（:5）‘;one;zero;zero;zero;zero];
A10=[zero;zero;zero;zero;Xp（:5）‘;Yp（:5）‘;Zp（:5）‘;one];
A11=[Xp（:6）‘;Yp（:6）‘;Zp（:6）‘;one;zero;zero;zero;zero];
A12=[zero;zero;zero;zero;Xp（:6）‘;Yp（:6）‘;Zp（:6）‘;one];
A1=A1‘;
A2=A2‘;
A3=A3‘;
A4=A4‘;
A5=A5‘;
A6=A6‘;
A7=A7‘;
A8=A8‘;
A9=A9‘;
A10=A10‘;
A11=A11‘;
A12=A12‘;
%A=[A1‘;A2‘];
A=[A1（1:）;A2（1:）;A1（2:）;A2（2:）;A1（3:）;A2（3:）;A1（4:）;A2（4:）;A1（5:）;A2（5:）;A1（6:）;A2（6:）;A1（7:）;A2（7:）;A1（8:）;A2（8:）;A3（1:）;A4（1:）;A3（2:）;A4（2:）;A3（3:）;A4（3:）;A3（4:）;A4（4:）;A3（5:）;A4（5:）;A3（6:）;A4（6:）;A3（7:）;A4（7:）;A3（8:）;A4（8:）;A5（1:）;A6（1:）;A5（2:）;A6（2:）;A5（3:）;A6（3:）;A5（4:）;A6（4:）;A5（5:）;A6（5:）;A5（6:）;A6（6:）;A5（7:）;A6（7:）;A5（8:）;A6（8:）;A7（1:）;A8（1:）;A7（2:）;A8（2:）;A7（3:）;A8（3:）;A7（4:）;A8（4:）;A7（5:）;A8（5:）;A7（6:）;A8（6:）;A7（7:）;A8（7:）;A7（8:）;A8（8:）;A9（1:）;A10（1:）;A9（2:）;A10（2:）;A9（3:）;A10（3:）;A9（4:）;A10（4:）;A9（5:）;A10（5:）;A9（6:）;A10（6:）;A9（7:）;A10（7:）;A9（8:）;A10（8:）;A11（1:）;A12（1:）;A11（2:）;A12（2:）;A11（3:）;A12（3:）;A11（4:）;A12（4:）;A11（5:）;A12（5:）;A11（6:）;A12（6:）;A11（7:）;A12（7:）;A11（8:）;A12（8:）];

Xd=[200 1 200 1 200 1 200 1 200 1 200 1 200 1 200 1 300 1 300 1 300 1 300 1 300 1 300 1 300 1 300 1 400 1 400 1 400 1 400 1 400 1 400 1 400 1 400 1 500 1 500 1 500 1 500 1 500 1 500 1 500 1 500 1 600 1 600 1 600 1 600 1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3168  2008-02-28 15:24  projectcalib\calibcomp.m

     文件       5452  2009-03-15 09:50  projectcalib\calib_data.mat

     文件       2353  2007-06-06 10:11  projectcalib\calib_gui.m

     文件       3768  2005-04-29 08:55  projectcalib\centerfinder.m

     文件        905  2008-05-05 17:02  projectcalib\cha.m

     文件        973  2002-09-30 14:32  projectcalib\check_active_images.m

     文件       4958  2001-11-27 12:26  projectcalib\check_directory.m

     文件        910  2005-04-27 08:17  projectcalib\check_extracted_images.m

     文件        105  2000-01-26 11:32  projectcalib\clear_windows.m

     文件       4728  2005-04-29 08:35  projectcalib\click_calib.m

     文件       8092  2005-04-29 08:58  projectcalib\click_ima_calib.m

     文件        466  2007-12-16 22:38  projectcalib\comp.m

     文件       3708  2000-05-15 17:03  projectcalib\compute_homography.m

     文件       1094  2006-09-06 10:51  projectcalib\Contents.m

     文件       2505  2007-06-06 10:49  projectcalib\data_calib.m

     文件        987  2008-03-22 15:47  projectcalib\data_process.m

     文件        780  2008-04-21 13:49  projectcalib\draw_1.m

     文件        447  2008-03-26 14:12  projectcalib\edgecolor.m

     文件        411  2008-03-22 15:13  projectcalib\eg_exa.m

     文件       1888  2005-04-29 08:25  projectcalib\export_calib_data.m

     文件        322  2008-03-10 14:27  projectcalib\ext_wu.m

     文件        668  2008-03-23 12:36  projectcalib\findcha_im.m

     文件        580  2001-10-22 11:49  projectcalib\fixvariable.m

     文件        399  2008-04-20 14:29  projectcalib\gauseidel.m

     文件       6312  2004-10-15 10:31  projectcalib\ginput2.m

     文件       6560  2005-04-29 10:51  projectcalib\ginput3.m

     文件       1864  2008-04-21 14:30  projectcalib\GRAY.m

     文件       1312  2008-04-21 14:33  projectcalib\hough.m

     文件       3282  2008-04-21 14:43  projectcalib\houghline.m

     文件    2044678  2008-03-06 20:34  projectcalib\Image000.bmp

............此处省略33个文件信息