单像空间后方交会实习报告.docxVIP

摄影测量学 单像空间后方交会实习报告 1 一、实习目的 1. 掌握空间后方交会的定义和实现算法 定义：空间后方交会是以单幅影像为基础，从该影像所覆盖地面范围内若干控 制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像 在航空摄影时刻的外方位元素 Xs，Ys，Zs，φ，ω，κ。 算法：由于每一对像方和物方共轭点可列出 2 个方程，因此若有 3 个已知地面 坐标的控制点，则可列出 6 个方程，解求 6 个外方位元素的改正数△Xs，△Ys，△ Zs，△φ，△ω，△κ。实际应用中为了提高解算精度，常有多余观测方程，通常是 在影像的四个角上选取 4 个或均匀地选择更多的地面控制点，因而要用最小二乘平差 方法进行计算。 2. 了解摄影测量平差的基本过程 获取已知数据。从摄影资料中查取影像比例尺 1/m，平均摄影距离（航空摄影的 航高）、内方位元素 x0，y0，f；获取控制点的空间坐标 X，Y，Z。 量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。 确定未知数的初始值。单像空间后方交会必须给出待定参数的初始值，在竖直航 空摄影且地面控制点大体对称分布的情况下，Xs0 和 Ys0 为均值，Zs0 为航高， φ、ω、κ 的初值都设为 0。或者 κ 的初值可在航迹图上找出或根据控制点坐 标通过坐标正反变换求出。 计算旋转矩阵 R。利用角元素近似值计算方向余弦值，组成 R 阵。 逐点计算像点坐标的近似值。利用未知数的近似值按共线条件式计算控制点像点 坐标的近似值（x），(y)。 逐点计算误差方程式的系数和常数项，组成误差方程式。 计算法方程的系数矩阵 ATA 与常数项 ATL，组成法方程式。 解求外方位元素。根据法方程，解求外方位元素改正数，并与相应的近似值求和， 得到外方位元素新的近似值。 检查计算是否收敛。将所求得的外方位元素的改正数与规定的限差比较，通常对 φ，ω，κ 的改正数△φ，△ω，△κ 给予限差，通常为 0.000001 弧度，当 3 个改正数均小于 0.000001 弧度时，迭代结束。否则用新的近似值重复(4)～（8） 步骤的计算，直到满足要求为止。 通过对提供的试验数据进行计算，输出像片的外方位元素并评定精度。深入理解 单片空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程 实现解求影像外方位元素的过程。通过上机调试程序加强动手能力的培养，通过 对实验结果的分析，增强综合运用所学知识解决实际问题的能力。 实习过程：4.1 学习单张像片空间后方交会的基本理论，掌握其基本思想。如果我 们知道每幅影像的 6 个外方位元素，就能确定被摄物体与航摄影像的关系。而单 像空间后方交会就是用于测定像片的外方位元素的，它的基本思想是：以单幅影 像为基础，从影像所覆盖的地面范围内若干控制点的已知地面坐标和相应点的像 坐标量测值出发，根据共线方程，解求该影像在航空摄影时刻的外方位元素 Xs,Ys,Zs,p,w,k.由于空间后方交会所采用的数学模型共线方程是非线性函数，为了 便于外方位元素的解求，首先将其线性化。4.2 在纸上绘出空间后方交会的计算机 程序框图。为了能够在宏观上指导我们编写程序，我们需要在草稿纸上绘出程序 2 框图。 3 输入原始数据归算像点坐标 x，y计 输入原始数据 归算像点坐标 x，y 计算和确定初值 Xs0，Ys0，Zs0，φ0，ω0，κ0 所有点完否？ 解法方程，求未知数改正数 未知数改正数限差否？ 整理并输出计算结果 正常结束 非正常结束 计算改正后的外方位元素 组成旋转矩阵 R 计算(x),(y)和 lx，ly 逐点组成误差方程式并法化 代 次 数 小 于 限 差 输出中间结果 和出错信息 是 否 迭 否？ 否  否 4 二、源代码 using System; using System.Collections.Generic; using System.Linq; using System.Text; namespace 单像空间的后方交会 { public class calculate { private double j, k, l, Xs, Ys, Zs;//六个外方位元素 private double f = 28.1539;//主距 struct point //像点和地面点坐标 { public double x, y, X, Y, Z; } private point[] p = new point[4];//存取控制点的坐标 private double[] R = new double[9];//旋转矩阵 private double[] a = new double[8];//近似值坐标 private double[] L = new d