摄影测量学后方交会matlab实习报告.docVIP

PAGE PAGE 1 摄影测量原理 单张影像后方交会实习  目录 529_WPSOffice_Level1 一 实习目的 3 26104_WPSOffice_Level1 二 实习原理 3 26104_WPSOffice_Level2 1. 间接平差 3 28026_WPSOffice_Level2 2. 共线方程 3 2047_WPSOffice_Level2 3. 单向空间后方交会 4 28026_WPSOffice_Level1 三 计算流程 4 4418_WPSOffice_Level2 1. 求解步骤 4 26420_WPSOffice_Level2 2.计算机框图 4 2047_WPSOffice_Level1 四 程序实现 5 4418_WPSOffice_Level1 五 结果分析 6 23452_WPSOffice_Level2 1.外方位元素 6 31072_WPSOffice_Level2 2.误差 6 22028_WPSOffice_Level2 3.旋转矩阵R 7 26420_WPSOffice_Level1 六 实习体会 7 26956_WPSOffice_Level2 1. 平台的选择 7 4655_WPSOffice_Level2 2.问题的解决 7 5202_WPSOffice_Level2 3.心得体会 8 23452_WPSOffice_Level1 七 代码展示 8  一 实习目的 为了增强同学们对后方交会公式的理解，培养同学们对迭代循环编程的熟悉感，本次摄影测量课间实习内容定为用C语言或其他程序编写单片空间后方交会程序，最终输出像点坐标、地面坐标、单位权中误差、外方位元素及其精度。 已知四对点的影像坐标和地面坐标如下。内方位元素fk=153.24mm，x0=y0=0。比例尺系数m=50000。 影像坐标 地面坐标 x(mm) y(mm) X(m) Y(m) Z(m) 1 -86.15 -68.99 36589.41 25273.32 2195.17 2 -53.40 82.21 37631.08 31324.51 728.69 3 -14.78 -76.63 39100.97 24934.98 2386.50 4 10.46 64.43 40426.54 30319.81 757.31 试计算近似垂直摄影情况下空间后方交会的解。 本次实习，我使用了matlab2014进行后方交会程序实现。结果与参考答案一致，精度良好。 二 实习原理 题干中有四个控制点在地面摄影测量坐标系中的坐标和对应的像点坐标，由此可列出8个误差方程，存在2个多余观测（n=2）。故可利用间接平差的最小二乘法则求解。 由于共线方程是非线性函数模型，为了方便计算，需要将其“线性化”。但如果仅取泰勒级数展开式的一次项，未知数的近似值改正是不精确的。因此必须采用迭代趋近法计算，直到外方位元素的改正值小于限差。 间接平差 间接平差为平差计算最常用的方法。在确定多个未知量的最或然值时，选择它们之间不存在任何条件关系的独立量作为未知量组成用未知量表达测量的函数关系、列出误差方程式，按最小二乘法原理求得未知量的最或然值的平差方法。 在一个间接平差问题中，当所选的独立参数X个数与必要观测值t个数相等时，可将每个观测值表达成这t个参数的函数，组成观测方程。 函数模型为：L = BX + d。 共线方程 共线方程是中心投影构像的数学基础，也是各种摄影测量处理方法的重要理论基础。 式中： x，y 为像点的像平面坐标； x0，y0，f 为影像的内方位元素； XS，YS，ZS 为摄站点的物方空间坐标； XA，YA，ZA 为物方点的物方空间坐标； ai，bi，ci （i = 1,2,3）为影像的 3 个外方位角元素组成的 9 个方向余弦。 单向空间后方交会 利用至少三个已知地面控制点的坐标A、B、Z与其影像上对应的三个像点的影像坐标a、b、c，根据共线方程，反求该像点的外方位元素。这种解算方法以单张相片为基础，亦称单像空间后方交会。 三 计算流程 求解步骤 （1）获取已知数据。 （2）量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。 （3）确定未知数的初始值。在竖直航空摄影且地面控制点大体对称分布的情况下，可按如下方法确定初始值： ZS0 = H = m*f XS0 = 1/n*(X1+X2+…+Xn) YS0 = 1/n*(Y1+Y2+…+Yn) q0 = w0 = 0 k0可在航迹图上找出，或根据控制点坐标通过坐标正反变换求出。 利用角元素近似值计算方向余弦值，组成旋转矩阵R。 利用未知数的近似值按共线方程条件式计算控制点像点坐标的近似值x0，y