摄影测量作业4.docxVIP

地球科学与环境工程学院摄影测量课后作业4课程名：《摄影测量测量学》学 号： 2014113248 姓 名： 黄黎东 指导老师： 陈 强 日 期： 2016年5月21日目录一、目的与要求3二、实验数学模型及算法分析3三、程序流程图5四、实验代码5五、实验结果截图13六、实验总结14实验目的掌握摄影测量立体像对的空间前方交会的计算方法及过程；比较采用点投影系数法和共线方程严密法计算地面点坐标的差异。二．实验数据及实验设备计算机windows10系统，编程软件VS2015，两张影像各自的外方位元素和框标坐标及相机的内方位元素存在文件data.excel中。已知条件左像片的外方位元素为: (14922.9778, 11532.4077, 3230.3383)m, (-0.0566, -0.1215, -0.6456)°，右像片的外方位元素为: (16296.9601, 11524.2256, 3239.6070)m, (0.0429, -0.5148, -0.1538)°，相机的内方位元素（x0, y0, f）为（-0.004，-0.008, 152.72）mm。现量测了6个像点的框标坐标，如下（单位mm）：左片: 1(85.7255, 69.6561), 2(0.3206, 59.1841), 3(83.7954, -1.9617), 4(-2.2603, -15.1519), 5(88.2371,-73.7021), 6(2.4792,-67.5320)右片:1 ( 0.6313, 70.1653), 2 ( -84.4506, 60.4553), 3( -3.3174, -1.3205)4( -97.0211, -13.7425), 5(2.7516, -72.5590), 6( -82.9161, -65.7818)采用点投影系数法和共线方程严密法，分别计算上述6个点的地面坐标。三．实验所用到的数学公式及程序计算步骤。（1）点投影系数法计算过程：根据两张影像量测得的框标坐标用旋转矩阵分别计算其对应的像空间辅助坐标，计算公式如下：根据两张影像的外方位元素计算摄影基线的三个分量，计算公式如下：计算两张影像的投影系数N1,N2，计算公式如下：计算两张影像各自测量得的框标坐标所对应的地面坐标，计算公式为：（2）共线方程严密法计算过程：根据共线条件方程推导基于共线严密方程的前方交会公式；共线条件方程为：根据共线条件方程推出共线条件严密方程形式如下：L1 X + L2 Y +L3 Z - Lx = 0L4 X + L5 Y +L6 Z - Ly = 0②求出共线严密方程的系数L1,L2, L3, L4, L5, L6, Lx, Ly.L1=fa1+(x-x0)a3L2=fb1+(x-x0)b3L3=fc1+(x-x0)c3L4=fa2+(y-y0)a3L5=fb2+(y-y0)b3L6=fc2+(y-y0)c3Lx=fa1Xs+fb1Ys+fc1Zs+(x-x0)a3Xs+(x-x0)b3Ys+(x-x0)c3ZsLy=fa2Xs+fb2Ys+fc2Zs+(y-y0)a3Xs+(y-y0)b3Ys+(y-y0)c3Zs③用平差方法计算出地面坐标根据严密方程得到系数阵为：A= L=根据最小二乘间接平差原理推导出如下式子：(ATPA)=ATP则=(ATA)-1ATL四．程序流程图。开始读取文件点投影法or严密法？严密方程法点投影法点投影法数学模型计算处理严密方程法数学模型计算处理输出计算结果结束五．程序的主要源代码如下所示：#include stdafx.h#include HomeWork4.h#include HomeWork4Dlg.h#include afxdialogex.h#include Matrix.h//矩阵类#ifdef _DEBUG#define new DEBUG_NEW#endifCMatrix xy1, xy2, XYZqwk(2,6);boolButtonOpenFile = false;BOOL CHomeWork4Dlg::OnInitDialog(){CDialogEx::OnInitDialog();strx0 = _T(-0.004);stry0 = _T(-0.008);strf = _T(152.72);UpdateData(false);DWORD dwStyle = m_list1.GetExtendedStyle();dwStyle |= LVS_EX_FULLROWSELECT;//选中某行使整行高亮（只适用与report风格的listctrl）dwStyle |= LVS_EX_GRIDLINES;/