文档介绍：中南大学
本科生课程设计(实践)任务书、设计报告
(摄影测量与遥感概论)
题目
空间后方-前方交会
学生姓名
指导教师
邹峥嵘
学院
地球科学与信息物理学院
专业班级
测绘0902班
学生学号
实验目的
通过对数字影像空间后交前交的程序设计实验,要求我们进一步理解和掌握影像外方位元素的有关理论、原理和方法。利用计算机程序设计语言编写摄影测量空间交会软件进行快速确定影像的外方位元素及其精度,然后通过求得的外方位元素求解未知点的地面摄影测量坐标,达到通过摄影测量量测地面地理数据的目的。
实验要求
用C、VB或者Matlab编写空间后方交会-前方交会计算机程序。
提交实验报告:程序框图,程序源代码、计算结果及体会。
计算结果:地面点坐标、外方位元素及精度。
完成时间:2011年11月17日。
实验数据
点号
左片
右片
地面摄影测量坐标
x
y
x
y
X
Y
Z
GCP1
-
GCP2
-
GCP3
-
-
-
GCP4
-
-
-
1
-
2
-
-
3
-
-
-
4
-
-
-
5
-
-
-
f=,x0=0,y0=0
实验思路
利用后方交会得出两张像片各自的外方位元素
获取已知数据:从摄影资料中插曲像片比例尺、平均航高、内方位元素以及控制点的地面摄影测量坐标及对应的像点坐标。
确定未知数的初始值:在竖直摄影的情况下,胶原素的初始值为0,线元素其中Zs=m*f+,Xs=,Ys=。
计算旋转矩阵R。
逐点计算像点坐标的近似值:利用共线方程。
组成误差方程并法化。
解求外方位元素。
检查计算是否收敛。
利用解求出的外方位元素进行前方交会
用各自像片的角元素计算出左右像片的旋转矩阵R1和R2。
根据左右像片的外方位元素计算摄影基线分量Bx,By,Bz。
逐点计算像点的空间辅助坐标。
计算投影系数。
计算未知点的地面摄影测量坐标。
重复以上步骤完成所有点的地面坐标的计算。
实验过程
程序流程框图
后方交会函数
确定已知数据比例尺m
确定各外方位元素初始值
计算旋转矩阵
逐点计算像点坐标近似值
不满足限差则重复计算
逐点计算误差方程系数项,组成误差系数矩阵A
利用矩阵运算求解外方位元素
检查是否满足限差
若满足则输出外方位元素
将整个过程作为一个函数继续进行右片的外方位元素求解
求解各外方位元素精度
空间前方交会
利用已求得的角元素计算2张像片各自的旋转矩阵
利用已求得的线元素Xs1,Ys1,Zs1,p0,,w01,k01;Xs2,Ys2,Zs2,p02,w02,k02,计算基线分量:
Bx=Xs2-Xs1;
By=Ys2-Ys1;
Bz=Zs2-Zs1;
输入像片坐标,利用旋转矩阵求解想空间辅助坐标
计算点投影系数:
N1=(Bx*Z2-Bz*X2)/(X1*Z2-X2*Z1);
N2=(Bx*Z1-Bz*X1)/(X1*Z2-X2*Z1);
计算地面摄影测量坐标
Xt=((N1*X1+Xs1)+(N2*X2+Xs2))/2;Yt=((N1*Y1+Ys1)+(N2*Y2+Ys2))/2;Zt=((N1*Z1+Zs1)+(N2*Z2+Zs2)) /2 ;
结束程序
程序中的主要函数设计
子函数(矩阵求积multiply,计算函数Resection,矩阵转置transpose,矩阵求逆inMerse1,输出函数shuchu,左片的外方位元素求解函数zuobian。右片的外方位元素求解函数youbian。)
程序源