双向解析光束法

双向解析光束法双向解析光束法双向解析光束法双向解析光束法编制仅供参考审核批准生效日期地址：电话：传真:邮编:双向解析光束法光束法程序有问题，在Getelement这个函数里便出现索引超限，这个问题一直解决不了光束法的流程：根据同名像点对对相交理论求系数阵A，系数阵B和常数阵La11=(a1f+a3x)/Z;a12=(b1f+b3x)/Z;a13=(c1f+c3x)/Z;a14=ysin(omega)-[x/f(xcos(kappa)-ysin(kappa))+fcos(kappa)]cos(omega);a15=-fsin(kappa)-x/f(xsin(kappa)+ycos(kappa));a16=y;a21=(a2f+a3y)/Z;a22=(b2f+b3y)/Z;a23=(c2f+c3y);a24=-xsin(omega)-[x/f(xcos(kappa)-ysin(kappa))-fsin(kappa)]cos(omega)a25=-fcos(kappa)-y/f(xsin(kappa)+ycos(kappa));a26=-x;求方程的改化法方程求出外方位元素和物方坐标改正数判断改正数的值，如果小于限差则输出结果光束法是最严密的一种方法的原因：在一张相片中，待定点与控制点的像点与摄影中心及相应地面点均构成一条光束，该方法是以每张相片所组成的一束光线作为平差的基本单元，已共线条件方程作为平差的基础方程，通过各个光束在空间中的旋转和平移，使模型之间公共点的光线实现最佳交汇，并使整个区域纳入到已知的地面控制点坐标系中，所以要建立全区域统一的误差方程，整体解求全区域内每张相片的六个外方位元素及所有待定点坐标，光束法区域网平差是基于摄影时像点，物点和摄站点三点共线提出来的。由单张相片构成区域，其平差的数学模型是共线条件方程，平差单元是单个光束，像点坐标是观测值，未知数是每张相片的外方位元素及所有待定点坐标。误差方程直接由像点坐标的观测值列出，能对像点坐标进行系统误差改正。光束法的程序代码为：ranspose();etElement(2,i))/2);}ultiply(A);Class1N12=().Multiply(B);Class1N21=().Multiply(A);Class1N22=().Multiply(B);Class1u1=().Multiply(LL);Class1u2=().Multiply(LL);ultiply()));ultiply(u2));ultiply(N11).Multiply(N12));Class1lX=().Multiply(N11).Multiply(u1));//求物方坐标改正数();Class1X=(lX);count++;//改正外方位元素和物方坐标phil+=(0,0);omegal+=(1,0);kappal+=(2,0);Xsl+=(3,0);Ysl+=(4,0);Zsl+=(5,0);phir+=(6,0);omegar+=(7,0);kappar+=(8,0);Xsr+=(9,0);Ysr+=(10,0);Zsr+=(11,0);//物方坐标改正数for(intj=0;j<8;j++){(0,j,(0,j)+(0+3*j,0));(1,j,(1,j)+(1+3*j,0));(2,j,(2,j)+(2+3*j,0));}//将改正数带出作结束判断for(intj=0;j<12;j++){dett[j]=(j,0);}for(intj=0;j<24;j++){detxyz[j]=(j,0);}for(inti=0;i<12;i++){if(det1<dett[i])det1=dett[i];}for(inti=0;i<24;i++){if(det2<detxyz[i])