资源描述：
《前方后方空间交会实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、中南大学本科生课程设计(实践)任务书、设计报告(摄影测量与遥感概论)题目空间后方-前方交会学生姓名指导教师邹峥嵘学院地球科学与信息物理学院专业班级测绘0902班学生学号一、实验目的通过对数字影像空间后交前交的程序设计实验，要求我们进一步理解和掌握影像外方位元素的有关理论、原理和方法。利用计算机程序设计语言编写摄影测量空间交会软件进行快速确定影像的外方位元素及其精度，然后通过求得的外方位元素求解未知点的地面摄影测量坐标，达到通过摄影测量量测地面地理数据的目的。二、实验要求Ø用C、VB或者Matlab编写空间后方交会-前方交会计
2、算机程序。Ø提交实验报告：程序框图，程序源代码、计算结果及体会。Ø计算结果：地面点坐标、外方位元素及精度。Ø完成时间：2011年11月17日。三、实验数据点号左片右片地面摄影测量坐标xyxyXYZGCP116.01279.963-73.9378.7065083.2055852.099527.925GCP288.5681.134-5.25278.1845780.025906.365571.549GCP313.362-79.37-79.122-78.8795210.8794258.446461.81GCP482.24-80.02
3、7-9.887-80.0895909.2644314.283455.484151.75880.555-39.95378.463 214.618-0.231-76.0060.036 349.88-0.782-42.201-1.022 486.14-1.346-7.706-2.112 548.035-79.962-44.438-79.736 f=150.000mm，x0=0，y0=0一、实验思路Ø利用后方交会得出两张像片各自的外方位元素1)获取已知数据：从摄影资料中插曲像片比例尺、平均航高、内方位元素以及控制
4、点的地面摄影测量坐标及对应的像点坐标。2)确定未知数的初始值：在竖直摄影的情况下，胶原素的初始值为0，线元素其中Zs=m*f+，Xs=，Ys=。3)计算旋转矩阵R。4)逐点计算像点坐标的近似值：利用共线方程。5)组成误差方程并法化。6)解求外方位元素。7)检查计算是否收敛。Ø利用解求出的外方位元素进行前方交会1)用各自像片的角元素计算出左右像片的旋转矩阵R1和R2。2)根据左右像片的外方位元素计算摄影基线分量Bx,By，Bz。3)逐点计算像点的空间辅助坐标。1)计算投影系数。2)计算未知点的地面摄影测量坐标。3)重复以上步骤
5、完成所有点的地面坐标的计算。一、实验过程Ø程序流程框图后方交会函数确定已知数据比例尺m确定各外方位元素初始值计算旋转矩阵逐点计算像点坐标近似值不满足限差则重复计算逐点计算误差方程系数项，组成误差系数矩阵A利用矩阵运算求解外方位元素检查是否满足限差若满足则输出外方位元素将整个过程作为一个函数继续进行右片的外方位元素求解求解各外方位元素精度空间前方交会利用已求得的角元素计算2张像片各自的旋转矩阵利用已求得的线元素Xs1，Ys1，Zs1，p0，,w01，k01；Xs2，Ys2，Zs2，p02，w02，k02，计算基线分量：Bx=X
6、s2-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，计算函数Resec
7、tion，矩阵转置transpose，矩阵求逆inMerse1，输出函数shuchu，左片的外方位元素求解函数zuobian。右片的外方位元素求解函数youbian。）Ø程序源代码#include"stdio.h"#include"math.h"doubleXs1,Xs2,Ys1,Ys2,Zs1,Zs2,p01,p02,w01,w02,k01,k02;//求矩阵a的转置矩阵b，a为m行、n列voidtranspose(double*a,double*b,intm,intn);//矩阵a乘以矩阵b，结果存储在c中，a为m×n大
8、小，b为n×l大小voidmultiply(double*a,double*b,double*c,intm,intn,intl);//求矩阵a的逆intinMerse1(double*a,intn);//输出m行、n列的矩阵avoidshuchu(double*a,intm,intn);/
