数字摄影测量编程实习.pptVIP

LOGO LOGO 数字摄影测量编程实习 成员：许劲松 答辩：许劲松 Page ? * 实习过程 实习目的： 用Visual C++编写一个完整的双相片空间后方交会程序，通过对提供的试验数据 进行计算，输出像片的内、外方位元素并评定精度。深入理解单双相片空间后方 交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像 内、外方位元素的过程。通过上机调试程序加强动手能力的培养，通过对实验结 果的分析，增强综合运用所学知识解决实际问题的能力。 实习内容： 利用小组成员提供的地面控制点，相点坐标根据共线条件方程求解像片内、外方位元素。 Page ? * 实习过程 数据准备： 初始值设定： 已知航摄仪内方位元素f＝18mm，Xo＝Yo＝0.0m 数据准备： 左相片 右相片 ? 影像 坐标 地面X坐标 地面Y坐标 地面Z坐标 ? x（m） y（m） X(m) Y（m） Z(m) 1 0.009744596842 0.009744596842 -0.2448 3.5173 -0.6169 2 0.01008552183 0.005509022325 -0.1723 4.2091 -0.3997 3 0.01011899767 0.0005140748393 -0.1544 4.6903 -0.2094 4 0.01660890563 0.01177869451 0.594 3.484 -0.6363 5 0.01990756196 0.005615371482 0.9314 4.2254 -0.533 6 0.01805103297 0.001033135935 0.7529 4.9947 -0.465 ? 影像 坐标 ?地面X坐标 地面坐标 地面Z坐标 ? x（mm） y（mm） X(m) Y（m） Z(m) 1 0.004838340665 0.01120032335 -0.2448 3.5173 -0.6169 2 0.004569731795 0.005438111155 -0.1723 4.2091 -0.3997 3 0.004063790431 0.0005957923502 -0.1544 4.6903 -0.2094 4 0.01150940753 0.01163543292 0.594 3.484 -0.6363 5 0.01455333476 0.005465707957 0.9314 4.2254 -0.533 6 0.01252037037 0.0009159152499 0.7529 4.9947 -0.465 Page ? * 实习过程 实习流程： 输入原始数据 归算像点坐标x,,y 计算和确定初值Xs0, Ys0, Zs0, φ0,ω0,κ0 组成旋转矩阵R 计算（x），（y）和lx，ly 逐点组成误差方程并法化 所有点完否？ 解法方程，求未知数改正数 计算改正后外方位元素 未知数改正数限差否？ 整理并输出计算结果 正常输出 迭代次数小于限差否？ ?输出中间结果和出错信息 非正常结束 Page ? * 程序片段 确定未知数的出始值： for(int i=0;i6;i++) {Xs=Xs+B[i][2]; Ys=Ys+B[i][3]; Zs=Zs+B[i][4]; } Xs=Xs/6; Ys=Ys/6; Zs=Zs/6+m*fk; 旋转矩阵R： R[0][0]=cos(Q)*cos(K)-sin(Q)*sin(W)*sin(K); R[0][1]=-cos(Q)*sin(K)-sin(Q)*sin(W)*cos(K); R[0][2]=-sin(Q)*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(Q)*cos(K)+cos(Q)*sin(W)*sin(K); R[2][1]=-sin(Q)*sin(K)+cos(Q)*sin(W)*cos(K); R[2][2]=cos(Q)*cos(W); Page ? * 程序片段 共线条件方程系数矩阵： A[k][0]=(R[0][0]*fk+R[0][2]*(B[i][0]-x0))/Z; A[k][1]=(R[1][0]*fk+R[1][2]*(B[i][0]-x0))/Z; A[k][2]=(R[2][0]*fk+R[2][2]*(B[i][0]-x0))/Z; A[k][3]=(B[i][1]-y0)*sin(W)-((B[i][0]-x0)*((B[i][0]-x0)*cos(K)-(B[i][1]-y0)*sin(K))/fk+fk*cos(K))*