单像空间后方交会实习报告.doc

单像空间后方交会实习报告 08012班 2008302590078 卢昊 程序流程图 是 是 否 是 否 读取数据 Xs, Ys, Zs, phi, omega, kappa 赋予初值 迭代次数是否超限 计算旋转矩阵R 计算x, y近似值和 Lx, Ly 计算误差方程系数并逐步法化 是否所有点都已计算 解法方程，求未知数改正数 是 改正数是否已达精度要求 计算改正后的外方位元素 输出或保存计算结果 否 正常退出 错误退出 部分代码 两个结构体，分别用来保存外方位元素和原始数据。 struct SourceData //保存原始数据{ struct SourceData //保存原始数据 { double x; double y; double X; double Y; double Z; }; struct EOEO // elements of exterior orientation { double Xs; double Ys; double Zs; double omega; double phi; double kappa; }; 几个主要函数 ////////////////////////////////////////////////////////////////////////// // 函数功能：初始化坐标数据 // 参数说明：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, N*sizeof(SourceData)); int i; for (i = 0; i != N; i ++) { datafile sd[i].x sd[i].y sd[i].X sd[i].Y sd[i].Z; } datafile.close(); cout 读取数据完毕...endl; } ////////////////////////////////////////////////////////////////////////// // 函数功能：检查改正数是否已达精度要求 // 参数说明：deta：保存改正数的数组 ////////////////////////////////////////////////////////////////////////// bool CheckPrecision(double* deta) { //2.9088820866572159615394846141477e-5 即 0.1′的弧度表示 bool ret; ret = (fabs(deta[0])0.000001 fabs(deta[1])0.000001 fabs(deta[2])0.000001 \ fabs(deta[3])2.90888208656e-5 fabs(deta[4])2.90888208656e-5 \ fabs(deta[5])2.90888208656e-5); return ret; } ////////////////////////////////////////////////////////////////////////// // 函数功能：迭代器，计算的主体部分 // 参数说明：sd：保存原始数据的结构体数组，PhotographicScale：摄影比例尺，focus：摄影机主距 ////////////////////////////////////////////////////////////////////////// void Iterator(SourceData sd[N], double PhotographicScale, double Focus) { double phi, omega, kappa, Xs, Ys, Zs; phi = omega = kappa =0.0;