直接线性变换数字摄影测量dlt

直接线性变换，从扫描坐标系到物方空间点直接的线性变换。有11参数l的初值计算，也有多个控制点的迭代计算，相机畸变参数也考虑进去。利用直接线性变换可以达到mm级精度，若是用DLT得到的初值，进行光束法平差，可以期待得到更高的精度。

function [fl1fl2fl3fl4fl5fl6fl7fl8fl9fl10fl11]=l__xishu__chuzhi（a1b1c1a2b2c2a3b3c3a4b4c4a5b5c5a6b6c6...
                                                                          xa1xb1xa2xb2xa3xb3xa4xb4xa5xb5xa6xb6）
% 计算l初值函数，左边输出11参数，右边，前6个对应地面空间坐标，后面x--像平面坐标
% 1 求B阵
B=[a1b1c110000xa1*a1xa1*b1xa1*c1;...
   0000a1b1c11xb1*a1xb1*b1xb1*c1;...
   a2b2c210000xa2*a2xa2*b2xa2*c2;...
   0000a2b2c21xb2*a2xb2*b2xb2*c2;...
   a3b3c310000xa3*a3xa3*b3xa3*c3;...
   0000a3b3c31xb3*a3xb3*b3xb3*c3;...
   a4b4c410000xa4*a4xa4*b4xa4*c4;...
   0000a4b4c41xb4*a4xb4*b4xb4*c4;...
   a5b5c510000xa5*a5xa5*b5xa5*c5;...
   0000a5b5c51xb5*a5xb5*b5xb5*c5;...
   a6b6c610000xa6*a6xa6*b6xa6*c6]; %小心，只需11个方程，舍去一个！
P=[xa1xb1xa2xb2xa3xb3xa4xb4xa5xb5xa6]‘;
ff=inv（B）*（-P）;
fl1=ff（1）;
fl2=ff（2）;
fl3=ff（3）;
fl4=ff（4）;
fl5=ff（5）;
fl6=ff（6）;
fl7=ff（7）;
fl8=ff（8）;
fl9=ff（9）;
fl10=ff（10）;
fl11=ff（11）;
kkk=0;
end










   
   
                           

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-12-03 16:09  DLT变换matlab\
     文件        1182  2012-12-02 21:03  DLT变换matlab\l__xishu__chuzhi.m
     文件        6065  2012-12-03 12:42  DLT变换matlab\l__xishu__diedai.m