资源描述：
《双像解析摄影测量作业.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
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)m
2、m。现量测了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,
3、-65。7818）请采用点投影系数法和共线方程严密法分别计算上述6个点的地面坐标。计算原理1、点投影系数法设地面点A在地面摄影测量坐标系中的坐标为（,对应像点在各自的像空间坐标系中的坐标为(）和（),在像空间辅助坐标系中的坐标分别为（和(.若已知两张像片的外方位元素，就可以由像点的像空间坐标系计算出该点的像空间辅助坐标,即=R，=R（3-1-1）式中：R、R为由已知的外方位元素算得的左右像片的旋转矩阵。右测站S在S-中的坐标，即摄影基线B的三个分量、、个人收集整理勿做商业用途可由外方位直线元素算得(
4、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）就是利用立体像对，在已知外方位元素的前提下，由像点坐标计算
5、对应地面空间坐标的前方交会公式。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、点投影系数法：个人收集整理勿做商业用途程序源代码：
6、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*
7、Math。PI/180，κ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=dou
8、ble.Parse(textBox3。Text）；doubley2=double。Parse（textBox7。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）*M
