空间后方交会报告.docVIP

空间后方交会报告 空间后方交会报告 空间后方交会报告 任务已知f=153.24mm，m=10000，限差0.1’各点坐标 点号 像点坐标 地面坐标 x（mm） y（mm） X（m） Y（m） Z（m） 1 -86.15 -68.99 36589.41 25273.32 2195.17 2 -53.40 82.21 37631.08 31324.51 728.69 3 -14.78 -76.63 39100.97 24934.98 2386.50 4 10.46 64.43 40426.54 30319.81 757.31 求近似垂直摄影情况下后方交会解 设计任务 确定未知数的初始值： Φ0 =ω0 =К0 = 0 , 内方位元素，，f=153.24mm。 ； = 38437m ； = 27963.16mbpurA6p。 计算旋转矩阵R 利用角元素的近似值计算方向余弦值，组成R阵 根据《摄影测量学》P32中的公式（3-9），初步计算R阵 R[0][0]=cos(Φ)*cos(K)-sin(Φ)*sin(W)*sin(K); R[0][1]=-cos(Φ)*sin(K)-sin(Φ)*sin(W)*cos(K); R[0][2]=-sin(Φ)*cos(W); R[1][0]=cos(W)*sin(K); R[1][1]=cos(W)*cos(K); R[1][2]=-sin(W); R[2][0]=sin(Φ)*cos(K)+cos(Φ)*sin(W)*sin(K); R[2][1]=-sin(Φ)*sin(K)+cos(Φ)*sin(W)*cos(K); R[2][2]=cos(Φ)*cos(W); 得初始R阵 逐点计算近似值（x），（y）： 带入《摄影测量学》P61的公式（5-1）；得 xUkuqcw。 组成误差方程式：按（5-8）；（5-9b）、（5-4）式逐点计算误差方程式的系数和常数项 根据Lx=x-(x)；Ly=y-(y)得 mJ26USm。 解得A阵为 pHlzBq6。 根据公式解除X阵，检验Φ1 、ω1 、К1 并不在限差0.1’之内;利用解得的Φ1 、ω1 、К1， 带入2中循环；知道解出合格结果！R8KoBFC。 结果，通过C++程序解得结果 源程序代码 #include iostream #include fstream #include cmath #include string #include iomanip using namespace std; const int n=6; void inverse (double c[n][n]); templatetypename T1,typename T2void transpose (T1*mat1,T2*mat2,int a,int b);sTcr9T6。 templatetypename T1,typename T2void multi(T1*mat1,T2 * mat2,T2 * result,int a,int b,int c);ghnVNB5。 templatetypename Tvoid input (T*mat,int a,int b);X0eOlwE。 templatetypename Tvoid output(T*mat,char*s,int a,int b);JLD0K61。 int main() { ofstream outFile; cout.precision(5); double x0=0.0, y0=0.0; double fk=0.15324; //内方位元素k9sHRJS。 double m=10000; //估算比例尺 double B[4][5]={0.0},R[3][3],XG[6][1],AT[6][8],ATA[6][6],ATL[6][1];JzIoqKg。 input (B,4,5); //从文件中读取控制点的影像坐标和地面坐标，存入数组B double Xs=0.0, Ys=0.0, Zs=0.0,Q=0.0,W=0.0,K=0.0; double X,Y,Z,L[8][1],A[8][6]; //确定未知数的出始值 for(int i=0;i4;i++) {Xs=Xs+B[i][2]; Ys=Ys+B[i][3]; Zs=Zs+B[i][4]; } Xs=Xs/4; Ys=Ys/4; Zs=Zs/4+m*fk; // coutXs Ys Zsendl; // int f=0; do//迭代计算 {f++; //组成旋转矩阵 R[0][0]=cos(Q)*cos(K)-sin(Q)*