单片空间后方交会实习报告

单片空间后方交会

程序设计

实习报告

班级： 学号： 姓名：

实习任务：用C或VC++语言实现单片后方交汇的计算。 实习目的：

1、深入理解单片空间后方交会的原理，

2、在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过

程。

3．通过上机调试程序加强动手能力的培养

实习原理

以单幅影象为基础，从该影象所覆盖地面范围内若干控制点的已知地面坐标和相应点的像 坐标量测值出发，根据共线条件方程，求解该影象在航空摄影时刻的像片外方位元素Xs， Ys，Zs,ф,ω,κ.共线条件方程如下：

x-x0=-f*[a1(X-Xs)+b1(Y-Ys)+c1(Z-Zs)]/[a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)] y-y0=-f*[a2(X-Xs)+b2(Y-Ys)+c2(Z-Zs)]/[a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)] x,y为像点的像平面坐标；x0,y0,f为影像的外方位元素；

Xs,Ys,Zs为摄站点的物方空间坐标；X,Y,Z为物方点的物方空间坐标；

. 实习主要仪器：电脑（创天中文VC++）

程序框图：

输入原始数据 归算像点坐标x,,y 计算和确定初值Xs0, Ys0, Zs0, φ0,ω0,κ0 组成旋转矩阵R 计算（x），（y）和lx，ly 逐点组成误差方程并法化 迭代次数小于限差否？ 所有点完否？ 解法方程，求未知数改正数 计算改正后外方位元素 未知数改正数<限差否？ 整理并输出计算结果 输出中间结果和出错信息 正常输出 非正常结束 实验数据：航摄仪的内方位元素：fk=153.24mm，x0=y0=0.0mm，摄影比例

尺为1:50000；

4个地面控制点的地面坐标及其对应像点的像片坐标：

结果输出：

已知条件: 像点坐标x,y: -53.4 82.21 -14.78 -76.63 10.46 64.43 153.24 -86.15 地点坐标Xa,Ya,Za:

37631.1 31324.5 728.69 39101 24935 2386.5 40426.5 30319.8 757.31 -68.99 36589.4 25273.3 内方位元素:

x0=y0=0 f=153.24mm

计算结果: 旋转矩阵：

0.997709 0.0675334 0.00398399 -0.0675254 0.997715 -0.00211178 -0.0041175 0.00183792 0.99999

像点坐标位:(单位：mm) -86.15 -68.99 -53.41 82.21 -14.78 -76.63 10.47 64.43

单位权中误差：7.2602632e-006 外方位元素：

Xs=39795.435 精度为:1.1254813 Ys=27476.479 精度为:1.2437625 Zs=7572.6929 精度为:0.48380521 q=-0.0039840098 精度为:0.00018182003 w=0.0021117837 精度为:0.00015959235 k=-0.067576934 精度为:7.2440432e-005 迭代次数：4

Press any key to continue

实习总结：

通过这次实习，在最初的程序中遇到了很多地困难，但最终还是克服了它。掌握运用程序设计实现单片后方交汇的计算，了解在多余观测情况下，用最小二乘平差方法计算相片外方位元素，及其相关精度。

本文来源：https://www.bwwdw.com/article/am7v.html