资源描述：
《摄影测量后方交会法求解外方位元素》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、摄影测量后方交会求外方位元素09地信一班肖明梅解题思路：定义要用到的变量并初始化，定义一个函数用于求解旋转矩阵R，系数矩阵A,近似坐标矩阵JSZB,常数矩阵L;矩阵的转置，逆，矩阵相乘，相减，求外方位元素的近似值初值以及结果输出都定义为相应的函数。最后定义一个用于循环求解的函数（程序中xhqiujie()）,在该函数中调用之前定义的函数，求出外方位元素近似值初值，改正数Dv[6,1],误差V[8,1]，用do…while语句进行循环，使精度达到20μm,限差低于20μm，然后调用结果输出函数用于
2、输出达到要求的结果。在主函数中创建对象的实例，引用该实例的方法即xhqiujie()函数，就可以求出外方位元素。代码：usingSystem;usingSystem.Collections.Generic;usingSystem.Linq;usingSystem.Text;namespace摄影测量后方交会求外方位元素{classProgram{doubleψ,ω,κ,Xs,Ys,Zs,m,f,Sx=0,Sy=0;doublen=206265/3600;//*角元素的单位从m到度的转换系数*//
3、double[,]zuobiao={{-86.15,-53.40,-14.78,10.46},{-68.99,82.21,-76.63,64.43},{36589.41,37631.08,39100.97,40426.54},{25273.32,31324.51,24934.98,30319.81},{2195.17,728.69,2386.50,757.31}};double[]a=newdouble[3];//*存储a1,a2.a3*//double[]b=newdouble[3];//*存
4、储b1,b2,b3*//double[]c=newdouble[3];//*存储c1,c2,c3*//double[]XX=newdouble[4];double[]YY=newdouble[4];double[]ZZ=newdouble[4];double[,]L=newdouble[8,1];double[,]JSZB=newdouble[2,4];//*像点近似坐标*//double[,]A=newdouble[8,6];//*系数阵*//double[,]AT=newdouble[6,8
5、];double[,]AB=newdouble[6,6];double[,]AN=newdouble[6,6];double[,]AL=newdouble[6,1];double[,]Dv=newdouble[6,1];double[,]AX=newdouble[8,1];//*系数阵与改正数矩阵的乘积*//double[,]V=newdouble[8,1];//*中误差矩阵*//privatevoidQchuzhi(){for(inti=0;i<4;i++){Sx+=zuobiao[2,i];
6、}for(inti=0;i<4;i++){Sy+=zuobiao[3,i];}for(inti=0;i<4;i++){zuobiao[0,i]=zuobiao[0,i]/1000;zuobiao[1,i]=zuobiao[1,i]/1000;}//*求外方位元素初始值*//ψ=ω=κ=0;m=50000;f=0.15324;Xs=Sx/4;Ys=Sy/4;Zs=m*f;}//*求旋转矩阵R,像点坐标近似值，常数项矩阵L,系数矩阵A*//privatevoidR(){//旋转矩阵Ra[0]=Mat
7、h.Cos(ψ)*Math.Cos(κ)-Math.Sin(ψ)*Math.Sin(ω)*Math.Sin(κ);a[1]=-Math.Cos(ψ)*Math.Sin(κ)-Math.Sin(ψ)*Math.Sin(ω)*Math.Cos(κ);a[2]=-Math.Sin(ψ)*Math.Cos(ω);b[0]=Math.Cos(ω)*Math.Sin(κ);b[1]=Math.Cos(ω)*Math.Cos(κ);b[2]=-Math.Sin(ω);c[0]=Math.Sin(ψ)*Math
8、.Cos(κ)+Math.Cos(ψ)*Math.Sin(ω)*Math.Sin(κ);c[1]=-Math.Sin(ψ)*Math.Sin(κ)+Math.Cos(ψ)*Math.Sin(ω)*Math.Cos(κ);c[2]=Math.Cos(ψ)*Math.Cos(ω);for(inti=0;i<4;i++){XX[i]=a[0]*(zuobiao[2,i]-Xs)+b[0]*(zuobiao[3,i]-Ys)+c[0]*(zuobiao[4,i]-Zs);YY[i]=a[1]*(zuob
