后方交会实验文档.docxVIP

PAGE \* MERGEFORMAT13 单向空间后方交会 实验报告 潘志富 2015221108130008 实验目的 以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像在航空摄影时刻的六个外方位元素，确定被摄物体与航摄影像的关系。 实验数据 摄影机主距f=153.24mm，x0=0.01mm，y0=0.02mm, 像片比例尺为1:40000，有四对点的像点坐标与相应的地面坐标如下表。 点号像点坐标地面坐标　x(mm)y(mm)X(m)Y(m)Z(m)1-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.4664.4340426.5430319.81757.31 算法流程 （1）获取已知数据。从航摄资料中差取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄影坐标。 （2）量测控制点的像点坐标并作系统误差改正。 （3）确定未知数的初始值。在竖直摄影且地面控制点大体对称分布的情况下，按如下方法确定初始值，即 ，， 用三个角元素的初始值按下式，计算各个方向余弦值，组成旋转矩阵R 逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面坐标；带入共线方程式，逐点近似像点坐标的近似值（x）、（y）。 逐点计算误差方程式的系数和常数项，组成误差方程式。 计算法方程的系数矩阵和常数项，组成法方程式。 解法方程，求得外方位元素的改正数。 用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素的新值。 将求得的外方位元素改正数与规定的限差比较，若小于限差，则迭代结束。负责用新的近似值重复（4）-（9），直到满足要求为止。 用共线方程进行空间后方交会的程序框如图所示。 源程序 #include stdio.h #include MatrixOperation.h #define N 4 //控制点数 #define PI 3.1415926 #define MAXITERATION 8 //最大允许迭代次数 struct EOEO // elements of exterior orientation { double Xs; double Ys; double Zs; double omega; double phi; double kappa; }; struct SourceData //保存原始数据 { double x; double y; double X; double Y; double Z; }; //函数声明 void InitInterface(); void InitData(double PhotographicScale, double focus, char* filename); bool CheckPrecision(double* deta); void Iterator(EOEO *eoeo, SourceData* sd, double PhotographicScale, double focus); void OutputResult(EOEO *eoeo, double* RotationMatrix, double* Pricision, double m0); ////////////////////////////////////////////////////////////////////////// // 函数功能：初始化坐标数据 // 参数说明：sd：保存原始数据的结构体数组，PhotographicScale：摄影比例尺，focus：摄影机主距 // filename：保存坐标数据的文件名 ////////////////////////////////////////////////////////////////////////// void InitData(SourceData* sd, char* filename) { //打开文件以备读取数据 cout 开始读取数据...endl; ifstream datafile(filename,ios::in | ios::nocreate); if( !datafile) { cout打开文件失败, 文件可能不存在endl; system(pause); _exit(1); } memset(sd, 0