文档介绍:中南大学本科生课程设计(实践)任务书、设计报告(摄影测量与遥感概论) 题目: 空间后方交会- 前交院系: 地球科学与信息物理学院班级: 测绘 1201 班指导老师: 邹峥嵘学号: 0103120506 姓名: 杨腊梅二零一四年十一月一、实验目的通过对数字影像空间后交前交的程序设计实验,要求我们进一步理解和掌握影像外方位元素的有关理论、原理和方法。利用计算机程序设计语言编写摄影测量空间交会软件进行快速确定影响的外方位元素及其精度,然后通过求得的外方位元素求解位置点的地面摄影测量坐标,达到通过摄影测量量测地面地理坐标的目的。二、实验要求 、 VB 、 C++ 或 MATLAB 语言编写空间后方交会- 空间前方交会程序 2. 提交实****报告:程序框图、程序源代码、计算结果、体会 3. 计算结果:像点坐标、地面坐标、单位权中误差、外方位元素及其精度 4. 完成时间: 11月 11 日前完成三、实验数据点号左片右片地面摄影测量坐标系 xyxyXYZ GCP1 - GCP2 - GCP3 - - - GCP4 - - - 1 - 2 - - 3 - - - 4 - - - 5 - - - f= , x0=0 , y0=0 四、实验思路?利用后方交会得出两张相片各自的外方位元素 1)获取已知数据影响比例尺 m,,内方位元素 x0、y0、f,控制点的地面摄影测量坐标 Xtp, Ytp, Ztp 2)量测控制点左片和右片的像点坐标 x,y 3 )确定未知数初值 Xs0, Ys0, Zs0, ω,φ,κ(线元素可用控制点均值代替,角元素可用 0初始化),即: ?? Xtp X4 1 s0 ,?? Ytp Y4 1 s0 ,fZ*ms?ω=0,φ=0,κ=0 4)计算旋转矩阵 R5)利用共线方程逐点计算像点坐标的近似值 6)组成误差方程式并法化 7)解求外方位元素改正数 8)检查迭代是否收敛(改正值是否小于某一特定常数) ?利用解出的外方位元素进行前方交会 1 )获取已知数据 x0,y0,f,XS1, YS1, ZS1, φ1,ω1,κ1,XS2, YS2, ZS2 ,φ2,ω2,κ22)量测像点坐标 x1,y1 ,x2,y2 3)由外方位线元素计算基线分量 BX,BY, BZ 4)由外方位角元素计算像空间辅助坐标 X1, Y1, Z1,X2, Y2, Z2 5)计算点投影系数 N1,N2 6)计算地面坐标 XA, YA, ZA 五、实验过程?程序流程图此过程完成空间后方交会求解像片的外方位元素,其中改正数小于限差( ,相当于 的角度值)为止。在这个过程中采用迭代的方法,是外方位元素逐渐收敛于理论值,每次迭代所得的改正数都应加到上一次的初始值之中。?程序中的类与其主要成员、主要函数设计 ControlPointData 类: 控制点、未知点的点名、地面摄影测量坐标、左右像片坐标 P hoto 类:左右片的内方位元素、外方位元素及其精度 Fileheader 类:相片对象、控制点个数、未知点个数、控制点数组、未知点数组 Adjustment 类: 头文件对象、控制点数组、未知点数组、控制点个数、未知点个数、 bool LoadData( const CString& strFileName,CString& strData); // 读取已知函数 CString* SplitString(CString str, char split, int & iSubStrs); // 分割函数 void adjustment1(CString &OutData); // 后方交会函数 void adjustment2(CString &OutData); // 前方交会函数 void SaveData(CString &OutData); // 保存结果?读取数据文件格式( .txt ) elements of interior orientation,0,0, Given