近景摄影测量实验报告.docx

课间实验报告 2015 年——2016 年第 2 学期 实验课程： 近 景 摄 影 测 量实验班级： 学生姓名： 学 号 ： 指导教师： 重庆交通大学测量与空间数据处理实验室 目 录 TOC \o 1-1 \h \z \u 实验一 近景单张像片空间后方交会 3 实验二 近景立体像对空间前方交会 8 实验三 直接线性变换求 L 系数精确解 13 实验四 直接线性变换求空间坐标精确值 17 实验一 近景单张像片空间后方交会 一、实验目的 通过单张像片空间后方交会算法编程，掌握空间后方交会的基本原理和基本步骤，为近景摄影测量解析处理方法打下理论基础。 二、实验内容 利用 VB 编程平台，通过给定的单片像点的像点坐标、内方位元素和控制点物方空间坐标，计算出像片的外方位元素。 三、实验步骤： 开始获取已知数据： 开始 获取已知数据： 比例尺，地面点坐标，像点坐标和内方位元素 定义旋转矩阵R 组成矩阵A和L 计算内方位元素改正值 否 判断内方位元素改 正值是否收敛 是 输出外方位元素 结束 2、程序代码 定义内方位元素变量Dim x0, yo, f As Double 定义外方位元素变量Dim Xs, Ys, Zs As Double Dim Q, w, k As Double 定义旋转变量元素 Dim a1, a2, a3 As Double Dim b1, b2, b3 As Double Dim c1, c2, c3 As Double 定义 A 矩阵 Dim a11, a12, a13, a14, a15, a16 As Double Dim a21, a22, a23, a24, a25, a26 As Double 定义像片共线方程变量 Dim X8 As Double, Y8 As Double, Z8 As Double 定义像点坐标变量Dim x, y, z As Double 定义像片 L 变量 Dim L1 As Double, L2 As Double 定义外方位角元素三角函数值Dim sk, sw, ck, cw As Double Private Sub Command1_Click() x0 = 48.6401 Txtx0 = x0 y0 = 11.4765 Txty0 = y0 f = 4548.5949 Txtf = f Xs = 3391.658 TxtXs = Xs Ys = -135.645 TxtYs = Ys Zs = 97.25 TxtZs = Zs TxtQ = -0°35′9.00″ Txtw = 3°33′13.00″ Txtk = -5°4′59.00″ End Sub Private Sub Command3_Click() 计算外方位角元素的三角函数值sw = 0.061982457 cw = 0.998077239 sk = -0.088599727 ck = 0.996067311 输入像点坐标x = 599.331 y = 1894.198 z = -4929.805 输入旋转矩阵 a1 = 0.995959: a2 = 0.089231: a3 = 0.010205 b1 = -0.088434: b2 = 0.994152: b3 = -0.061983 c1 = -0.015676: c2 = 0.06083: c3 = 0.998025 计算像片共线方程系数 X8 = a1 * (x - Xs) + b1 * (y - Ys) + c1 * (z - Zs) Y8 = a2 * (x - Xs) + b2 * (y - Ys) + c2 * (z - Zs) Z8 = a3 * (x - Xs) + b3 * (y - Ys) + c3 * (z - Zs) 计算矩阵 A a11 = (a1 * f + a3 * (x - x0)) / Z8 a12 = (b1 * f + b3 * (x - x0)) / Z8 a13 = (c1 * f + c3 * (x - x0)) / Z8 a14 = (y - y0) * sw - ((x - x0) * ((x - x0) * ck - (y - y0) * sk) / f + f * ck) * cw a15 = -f * sk - (x - x0) * ((x - x0) * sk + (y - y0) * ck) / f a16 = y - y0 a21 = (a2 * f + a3 * (y - y0)) / Z8 a22 = (b2 * f + b3 * (y - y0)) / Z8 a23 = (c2 * f + c3 * (y - y0)) / Z8 a24 = -(x - x0) * sw