双像解析摄影测量习题.docxVIP

PAGE PAGE 10 / 11 双像解析摄影测量作业 报 告 姓名：*** 学号：*** 班级：*** 作业： 已知一个立体像对： 左像片的外方位元素为: (14922.9778, 11532.4077, 3230.3383)m, (-0.0566, -0.1215, -0.6456)°， 右像片的外方位元素为 : (16296.9601, 11524.2256, 3239.6070)m, (0.0429, -0.5148, -0.1538)°， 相机的内方位元素（x0, y0, f）为（-0.004，-0.008, 152.72）mm。现量测了 6 个像点的框标坐标，如下（单位mm）： 左 片 :1(85.7255,69.6561), 2(0.3206,59.1841),3(83.7954, -1.9617),4(-2.2603 ， -15.1519), 5(88.2371,-73.7021), 6(2.4792,-67.5320) 右片: 1 ( 0.6313, 70.1653), 2 ( -84.4506, 60.4553), 3( -3.3174, -1.3205) 4( -97.0211, -13.7425), 5(2.7516, -72.5590), 6( -82.9161, -65.7818) 请采用点投影系数法和共线方程严密法分别计算上述 6 个点的地面坐标。 计算原理 1、点投影系数法 设地面点 A 在地面摄影测量坐标系中的坐标为（,对应像点在各自的像空间坐标系中的坐标为（）和（）,在像空间辅助坐标系中的坐标分别为（和（。若已知两张像片的外方位元素，就可以由像点的像空间坐标系计算出该点的像空间辅助坐标，即 =R,=R （3-1-1） 式中：R、R 为由已知的外方位元素算得的左右像片的旋转矩阵。右测站S 在S-中的坐标，即摄影基线B 的三个分量、、可由外方位直线元素算得 (3-1-2) 因左、右像空间辅助坐标系与地面摄影测量坐标系相互平行，且摄站 点、像点、地面点三点共线在，则由式(3-1-1)可得 （3-1-3） 式中：、为左右像点的点投影系数。一般情况下，不同的点有不同的点投影系数值。根据式（3-1-3）可以得到前方交会法计算地面点坐标的公式，即 （3-1-4） 式（3-1-4）中、仍然未知，为此，结合（3-1-1）有 （3-1-5） 由式(3-1-5)中的一、三两式联立解得 （3-1-6） 式（3-1-4）和式（3-1-6）就是利用立体像对，在已知外方位元素的前提下，由像点坐标计算对应地面空间坐标的前方交会公式。 2、共线方程严密法 一个地面点分别对应左右两张像片两个点，每个像点对应两个共线方程，即一个地面未知点对应四个共线方程。再将四个共线方程分别对X、Y、Z 进行线性化得到关于 X、Y、Z 三个未知数的四个误差线性方程，即 （3-2-1） 进一步表示为 式中（i=1·· 12）,、、分别表示为 其中a1,a2,a3,b1,b2,b3,c1,c2,c3为旋转矩阵R中的元素. 可表示为：V=AX-L 解得： 计算过程： 1、 点投影系数法： 程序源代码： private void button1_Click(object sender, EventArgs e) { double Xs1 = 14922.9778, Ys1 = 11532.4077, Zs1 = 3230.3383,φ1 = -0.0566 * Math.PI / 180, ω1 = -0.1215 * Math.PI / 180, κ1 = -0.6456 * Math.PI / 180; double Xs2 = 16296.9601, Ys2 = 11524.2256, Zs2 = 3239.6070,φ2 = 0.0429 * Math.PI / 180, ω2 = -0.5148 * Math.PI / 180, κ2 = -0.1538 * Math.PI; double f = 0.15272, X, Y, Z,Bx,By,Bz,N1,N2,X1,X2,Y1,Y2,Z1,Z2; double x0 = -0.004, y0 = -0.008; double[,] R1 = new double[3, 3]; double[,] R2 = new double[3, 3]; double x1 = double.Parse(textBox1.Text); double y1 = double.Parse(textBox2.Text); double x2 = double.Parse(textBox3.Text); double y2 = double.Parse(textBox7.Text); double[,]