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

- - PAGE # - 摄影测量原理 单张影像后方交会实习 目录 TOC \o 1-5 \h \z \o Current Document 实习目的 3 \o Current Document 实习原理 3 间接平差 3 共线方程 3 单向空间后方交会 4... \o Current Document 计算流程 4 . 求解步骤 4 .计算机框图 4 \o Current Document 程序实现 5 \o Current Document 结果分析 6 .外方位元素 6 .误差 6 .旋转矩阵 R 7... 实习体会 7 . 平台的选择 7... .问题的解决 7 .心得体会 8 代码展示 8 - - PAGE # - 一实习目的 为了增强同学们对后方交会公式的理解，培养同学们对迭代循环编程的熟悉 感，本次摄影测量课间实习内容定为用 C语言或其他程序编写单片空间后方交会 程序，最终输出像点坐标、地面坐标、单位权中误差、外方位元素及其精度。 已知四对点的影像坐标和地面坐标如下。 内方位元素fk=153.24mm, x0=y0=0。 比例尺系数m=50000o 影像坐标 地面坐标 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 .共线方程 共线方程是中心投影构像的数学基础，也是各种摄影测量处理方法的重要理 论基础。 _ _ r —（孙 -x—）一瓦（匕-为）口（工力一 ） 0 fl3 （Xji - ） + 与（匕-）+ 白（z5- zs） 0 町区I + I 匕-匕）+ 向（2/-Zs） 式中： 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)确定未知数的初始值。在竖直航空摄影且地面控制点大体对称分布的情况下， 可按如下方法确定初始值： Zs = H = m*f Xs0 = 1/n*(Xi++…+M) Ys0 = 1/n*(Yi+Y2+…+询 q0 = w0 = 0 k0可在航迹图上找出，或根据控制点坐标通过坐标正反变换求出。 (4)利用角元素近似值计算方向余弦值，组成旋转矩阵 R。 (5)利用未知数的近似值按共线方程条件式计算控制点像点坐标的近似值 X0,y。 (6)逐点计算误差方程式的系数和常数项，组成误差方程式。 (7)解求外方位元素改正数，并与相应的近似值求和，得到外方位元素新的近似 值。 (8)规定限差限制迭代次数，检查计算是否收敛。 .计算机框图 riW*A?k^r迭代次粗小FM差管量出有程.束束以歌也正敷4所刷电宛者7本知我为正把父用胃舌1计W正好的外方位元?告善再井箱出计算站果 riW*A?k^r 迭代次