作业报告
空间后方交会
专业：测绘工程
班级： 2008级(1)班姓名：陈闻亚
指导教师：陈强
2010 年 4 月16 日
1 作业任务------------------------------------------------------------------------------------ 3
2 作业思想--------------------------------------------------------------------------------------- 3
3 作业条件及数据-------------------------------------------------------------------- 3
4 作业过程--------------------------------------------------------------------------- 3
5 源程序----------------------------------------------------------------------------- 4
6 计算结果--------------------------------------------------------------------------- 1 7 7心得体会与建议----------------------------------------------------------------------------- 17
1 作业任务
计算近似垂直摄影情况下后方交会解。即利用摄影测量空间后方交会的方法，获取相片的6个外方位元素。限差为0.1。 2作业思想
利用摄影测量空间后方交会的方法求解。该方法的基本思想是利用至少三个一直地面控制点的坐标A（X A，Y A，Z A）、B（X B，Y B，Z B）C（X C，Y C，Z C），与其影像上对应的三个像点的影像坐标a（x a，y a）、b（x b，y b）、c（x c，y c），根据共线方程，反求该相片的外方位元素X S、Y S、Z S、φ、ω、κ。 3作业条件及数据
已知摄影机主距f=153.24mm，四对点的像点坐标与相应的地面坐标列入下表：
表1
4作业过程
4．1 获取已知数据
相片比例尺1/m=1：10000，内方位元素f=153.24mm，x0，y0；获取控制点的地面测量坐标X t、Y t、Z t。
4．2 量测控制点的像点坐标：
本次作业中为已知。见表1。
4．3 确定未知数的初始值：
在近似垂直摄影情况下，胶原素的初始值为0，即φ0 = ω0 = κ0=0；线元素中，Z S0=H=mf=1532.4m，X S0、Y S0的取值可用四个控制点坐标的平均值，即：
X S0= 1
4
4
i
i1
X
=
∑=38437.00
Y S0= 1
4
4
i
接地母排
i1
Y
=
∑=89106.62
4．4 计算旋转矩阵R：
利用胶原素的近似值计算方向余弦值，组成R阵。
4．5 逐点计算像点坐标的近似值：
利用未知数的近似值按共线方程式计算控制点像点坐标的近似值（x）（y）。4．6 组成误差方程：
逐点计算误差方程式的系数和常数项。
移动隔断墙
4．7 组成法方程式：
计算法方程的系数矩阵A T A与常数项A T L。
4．8 求解外方位元素：
牙刷消毒器根据法方程，由式X=（A t A）-1 A T L解求外方位元素改正数，并与相应的近似值求和，得到外方位元素新的近似值。
4．9 求解外方位元素：平衡梁
将求得的外方位元素的改正数与规定的限差（0.1）比较，小于限差则计算终止，否则用新的近似值重复第4.4至4.8步骤的计算，知道满足要求为止。
5 源程序
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
const double PRECISION=1e-5;
typedef double DOUBLE[5];
int InputData(int &Num, DOUBLE *&Data,double &m,double &f);
int Resection(const int &Num,const DOUBLE *&Data,const double &m,const double &f);
int InverseMatrix(double *matrix,const int &row);
int main(int argc, char* argv[])
{
DOUBLE *Data=NULL;
int Num;
double f(0),m(0);
if(InputData(Num,Data,m,f))
{
if (Data!=NULL)
{
delete []Data;
udiab}
return 1;
}
if(Resection(Num,Data,m,f))
{
if (Data!=NULL)
{
delete []Data;
}
return 1;
}
if (Data!=NULL)
{
delete []Data;
}
printf("解算完毕...\n");
ap劫do{
printf("计算结果保存于\"结果.txt\"文件中\n"
"请选择操作（输入P打开结果数据，R打开原始数据，其它退出程序）：");
fflush(stdin); //刷新输入流
char order=getchar();
if ('P'==order || 'p'==order)
{
system("结果.txt");
}
else if ('R'==order || 'r'==order)
{
system("");
}
else
break;
system("cls");
}while(1);
system("PAUSE");
return 0;
}
/**********************************************
*函数名：InputData
*函数介绍：从文件（）中读取数据,
*文件格式如下：
*点数m(未知写作0)
* 内方位元素(f x0 y0)
*编号x y X Y Z
*实例：
4 0
153.24 0 0
1 -86.15 -68.99 36589.41 25273.3
2 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
*参数：(in/out)Num(点数)，
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议