文档介绍：空间后交-前交程序设计实验报告
空间后交-前交程序设计实验报告
1 / 11
空间后交-前交程序设计实验报告
空间后交 - 前交程序设计
( 实验报告 )
姓名：
班级：
学号：
时间：
空间后交-前交程序设计实验报告
空间后交-前交程序设计实验报告
11 / 11
空间后交-前交程序设计实验报告
空间后交 - 前交程序设计
一、实验目的
用 C 、VB或 MATLAB语言编写空间后方交会 - 空间前方交会程序⑴提交实****报告：程序框图、程序源代码、计算结果、体会⑵计算结果：像点坐标、地面坐标、单位权中误差、外方位元素及其精度
二、实验数据
f=，x0=0，y0=0
三、实验思路
利用空间后方交会求左右像片的外方位元素
(1). 获取 m(于像片中选取两点，于地面摄影测量坐标系中选取同点，
分别计算距离，距离比值即为 m)，x,y,f,X,Y,Z
(2). 确定未知数初始值 Xs,Ys,Zs,q,w,k
(3). 计算旋转矩阵 R
(4). 逐点计算像点坐标的近似值（ x），（ y）
空间后交-前交程序设计实验报告
空间后交-前交程序设计实验报告
3 / 11
空间后交-前交程序设计实验报告
(5). 组成误差方程式
(6). 组成法方程式
(7). 解求外方位元素
(8). 检查是否收敛，即将求得的外方位元素的改正数与规定限差比
较，小于限差即终止；否则用新的近似值重复步骤（ 3）- （7）
利用求出的外方位元素进行空间前交，求出待定点地面坐标
(1). 用各自像片的角元素计算出左、右像片的方向余弦值，组成旋转
矩阵 R1,R2
(2). 根据左、右像片的外方位元素，计算摄影基线分量 Bx，By，Bz
(3). 计算像点的像空间辅助坐标（ X1,Y1,Z1 ）和（ X2,Y2,Z2 ）
(4). 计算点投影系数 N1和 N2
(5). 计算未知点的地面摄影测量坐标
四、实验过程
⑴程序框图
输入
testvar
spacehoujiao
AAndL
计算
spaceqianiao
输出 deg2dmsxy
ok
⑵程序代码
函数 AandL
空间后交-前交程序设计实验报告
空间后交-前交程序设计实验报告
4 / 11
空间后交-前交程序设计实验报告
%求间接平差时需要的系数
空间后交-前交程序设计实验报告
空间后交-前交程序设计实验报告
11 / 11
空间后交-前交程序设计实验报告
%%%已知
%a=像点坐标 x,b= 像点坐标 y，f 内方位元素主距
%φ=q, ψ=w,κ=k
%像空间坐标系 X,Y,Z
%地面摄影测量坐标系 Xs,Ys,Zs
function [A1,L1,A2,L2]=AandL(a,b,f,q,w,k,X,Y,Z,Xs,Ys,Zs)
%%%%%%%%%%%选择矩阵元素
a1=cos(q)*cos(k)-sin(q)*sin(w)*sin(k);
a2=-cos(q)*sin(k)-sin(q)*sin(w)*cos(k);
a3=-sin(q)*cos(w);
b1=cos(w)*sin(k);
b2=cos(w)*cos(k);
b3=-sin(w);
c1=sin(q)*cos(k)+cos(q)*sin(w)*sin(k);
c2=-sin(q)*sin(k)+cos(q)*sin(w)*cos(k);
c3=cos(q)*cos(w);
%%%%%%%共线方程的分子分母
X_=a1*(X-Xs)+b1*(Y-Ys)+c1*(Z-Zs);
Y_=a