资源描述：
《摄影测量作业4》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、地球科学与环境工程学院摄影测量课后作业4课程名：《摄影测量测量学》学号：2014113248姓名：黄黎东指导老师：陈强日期：2016年5月21日15目录一、目的与要求3二、实验数学模型及算法分析3三、程序流程图5四、实验代码5五、实验结果截图13六、实验总结1415一．实验目的（1）掌握摄影测量立体像对的空间前方交会的计算方法及过程；（2）比较采用点投影系数法和共线方程严密法计算地面点坐标的差异。二．实验数据及实验设备计算机windows10系统，编程软件VS2015，两张影像各自的外方位元素和框标坐标及相机的内方位元素存在文件data.excel中。已知条件左像
2、片的外方位元素为:(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.151
3、9),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）点投影系数法计算过程：①根据两张影像量测得的框标坐标用旋转矩阵分别计算其对应的像空间辅助坐标，计算公式如下：②根据两张影像的外方位元素计算摄
4、影基线的三个分量，计算公式如下：③计算两张影像的投影系数N1,N2，计算公式如下：④计算两张影像各自测量得的框标坐标所对应的地面坐标，计算公式为：15（2）共线方程严密法计算过程：①根据共线条件方程推导基于共线严密方程的前方交会公式；共线条件方程为：根据共线条件方程推出共线条件严密方程形式如下：L1X+L2Y+L3Z-Lx=0L4X+L5Y+L6Z-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
5、-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=L1L2L3L4L5L6L=LxLy根据最小二乘间接平差原理推导出如下式子：(ATPA)XYZ=ATPL则XYZ=(ATA)-1ATL四．程序流程图。15开始读取文件点投影法or严密法？严密方程法点投影法点投影法数学模型计算处理严密方程法数学模型计算处理输出计
6、算结果结束五．程序的主要源代码如下所示：#include"stdafx.h"#include"HomeWork4.h"#include"HomeWork4Dlg.h"#include"afxdialogex.h"#include"Matrix.h"//矩阵类#ifdef_DEBUG#definenewDEBUG_NEW#endifCMatrixxy1,xy2,XYZqwk(2,6);boolButtonOpenFile=false;BOOLCHomeWork4Dlg::OnInitDialog(){CDialogEx::OnInitDialog();strx0=_
7、T("-0.004");stry0=_T("-0.008");strf=_T("152.72");UpdateData(false);DWORDdwStyle=m_list1.GetExtendedStyle();15dwStyle
8、=LVS_EX_FULLROWSELECT;//选中某行使整行高亮（只适用与report风格的listctrl）dwStyle
9、=LVS_EX_GRIDLINES;//网格线（只适用与report风格的listctrl）m_list1.SetExtendedStyle(dwStyle);//设置扩展风格m_list1.InsertCo
10、lumn(
