空间后方交会编程 单张像片解析 摄影测量教研室.pptVIP

主要内容 一、数据 二、误差方程 误差方程系数 误差方程系数 误差方程系数 误差方程系数 外方位元素结构定义 旋转矩阵 共线条件方程(误差方程)程序 共线条件方程(误差方程)程序 共线条件方程(误差方程)程序 法方程 法方程 法方程解算 精度评定 计算过程 计算结果 计算结果 课间实习 * * 空间后方交会编程 摄影测量教研室 赵双明 smzhao2000@ 第三章 单张像片解析 一、数据 二、误差方程 三、法方程 四、解算法方程 五、精度评定 航摄仪内方位元素：f = 153.24mm，x0 = y0 = 0.0mm 摄影比例尺1：m = 1:50000 757.31 30319.81 40426.54 64.43 10.46 4 2386.50 24934.98 39100.97 -76.63 -14.78 3 728.69 31324.51 37631.08 82.21 -53.40 2 2195.17 25273.32 36589.41 -68.99 -86.15 1 Z (m) Y (m) X (m) y (mm) x (mm) 控制点坐标 像片坐标 点号 以Y轴为主轴的?-?-?转角系统 typedef struct{ double Xs; double Ys; double Zs; double omega; double phi; double kappa; }ORISTRUCT; void RotationMatrix(double omega, double phi, double kappa, double *matrix) { matrix[0]= cos(phi)*cos(kappa)-sin(phi)*sin(omega)*sin(kappa); matrix[1]= -cos(phi)*sin(kappa)-sin(phi)*sin(omega)*cos(kappa); matrix[2]= -sin(phi)*cos(omega); matrix[3]= cos(omega)*sin(kappa); matrix[4]= cos(omega)*cos(kappa); matrix[5]= -sin(omega); matrix[6]= sin(phi)*cos(kappa)+cos(phi)*sin(omega)*sin(kappa); matrix[7]= -sin(phi)*sin(kappa)+cos(phi)*sin(omega)*cos(kappa); matrix[8]= cos(phi)*cos(omega); } void ErrorEquation(ORISTRUCT ori, double f, double x, double y, double X, double Y, double Z, double *Bx, double *By, double *lx, double *ly) { //参数初值 double Xs = ori.Xs, Ys = ori.Ys, Zs = ori.Zs; double omega = ori.omega, phi = ori.phi, kappa = ori.kappa; //旋转矩阵的计算 double R[9]; RotationMatrix(omega, phi, kappa, R); //辅助参数计算 double U,V,W; U = R[0]*(X-Xs)+R[3]*(Y-Ys)+R[6]*(Z-Zs); V = R[1]*(X-Xs)+R[4]*(Y-Ys)+R[7]*(Z-Zs); W = R[2]*(X-Xs)+R[5]*(Y-Ys)+R[8]*(Z-Zs); //x: 偏导数计算 Bx[0]=(R[0]*f+R[2]*x)/W;//Xs Bx[1]=(R[3]*f+R[5]*x)/W;//Ys Bx[2]=(R[6]*f+R[8]*x)/W;//Zs Bx[3]= y*sin(omega)-(x*(x*cos(kappa)- y*sin(kappa))/f+f*cos(kappa))*cos(omega);//phi Bx[4]=-f*sin(kappa)-x*(x*sin(kappa)+y*cos(kappa))/f;//omega Bx[5]= y;//kappa