资源描述：
《双像解析摄影测量作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、双像解析摄影测量作业报告姓名：***学号：***班级：***作业：已知一个立体像对：左像片的外方位元素为:(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),
2、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在地面摄影测量坐标系中的坐标为（,对应像点在各自的像空间坐标系中的坐标为（）和（）,在像空间辅助
3、坐标系中的坐标分别为（和（。若已知两张像片的外方位元素，就可以由像点的像空间坐标系计算出该点的像空间辅助坐标，即=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）中、仍然未知，
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
5、=AX-L解得：计算过程：1、点投影系数法：程序源代码：privatevoidbutton1_Click(objectsender,EventArgse){doubleXs1=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;doubleXs2=16296.9601,Ys2=11524.2256,Zs2=3239.6070,φ2=0.0429*Math.PI/180,ω2=-0.5148*Math.PI/180,κ
6、2=-0.1538*Math.PI;doublef=0.15272,X,Y,Z,Bx,By,Bz,N1,N2,X1,X2,Y1,Y2,Z1,Z2;doublex0=-0.004,y0=-0.008;double[,]R1=newdouble[3,3];double[,]R2=newdouble[3,3];doublex1=double.Parse(textBox1.Text);doubley1=double.Parse(textBox2.Text);doublex2=double.Parse(textBox3.Text);doubley2=double.Parse(textBox
7、7.Text);double[,]A1=newdouble[3,1];double[,]A2=newdouble[3,1];double[,]B1=newdouble[3,1];double[,]B2=newdouble[3,1];x1=x1-x0;y1=y1-y0;R1[0,0]=Math.Cos(φ1)*Math.Cos(κ1)-Math.Sin(φ1)*Math.Sin(ω1)*Math.Sin(κ1);//a1R1[0,1]=-Math.Cos(φ1)*Math.Sin(
