光束法区域网空中三角测量作业教程解读.docVIP

光束法区域网空中三角测量求像片外方位元素和加密点的地面坐标 一、原理：首先给出的原始数据像素坐标转换成像平面直角坐标；然后利用后方交会原理并采用控制点的向平面直角坐标和地面坐标迭代求解每张像片的外方位元素；再通过前方交会求加密点的地面坐标；最后进行精度评定。 二、计算流程： 求出四个控制点及所有待求点的像点坐标（x1,y1）与(x2,y2)； 空间后方交会计算两像片的外方位元素：根据编好的计算机程序读取原始文件中的数据控制点的地面坐标和相应的像点像素坐标（第一步已经转换成像平面直角坐标），对两张像片各自进行空间后方交会，计算各自的6个外方位元素Xs1,Ys1,Zs1,ψ1,w1,k1和Xs2,Ys2,Zs2,ψ2,w2,k2； 空间前方交会计算待定点的地面坐标：用各片的外方位角元素计算左右片的方向余弦值，组成旋转矩阵R1与R2，逐点计算像点的像空间爱你辅助坐标(u1,v1,w1)及(u2,v2,w2)；根据外方位元素计算基线分量（Bu,Bv,Bw）;计算投影系数N1,N2；计算待定点在各自的像空间辅助坐标系中的坐标(U1,V1,W1)及(U2,V2,W2)；最后计算待定点的地面摄影测量坐标。 进行精度评定：利用前方交会求出地面控制点的坐标，再与地面实测坐标求差，将两者差值视为真差值，由这些真差值计算出点位坐标精度。 三、基本公式： 1、像素坐标（I,J）转换成像平面直角坐标(x,y)： x = (J- 2136) * 5.19 / 1000 / 1000; y = (1424 - I) * 5.19 / 1000 / 1000; 2、后方交会： 共线条件方程： 外方位元素系数矩阵： Zp=a3*(L_X-Xs)+b3 *(L_Y-Ys)+c3 *(L_Z-Zs); a11= 1/Zp*(a1*f+a3*x); a12 =1/Zp*(b1*f+b3*x); a13= 1/Zp*(c1*f+c3*x); a14= y* Sin (w)-(x /f*(x*Cos(k)-y* Sin (k))+f* Cos (k))* Cos (w ); a15= -f* Sin (k)-x /f*(x*Sin (k)+y *Cos (k)); a16 = y; a21 = 1/Zp*(a2*f+a3*y); a22 = 1/Zp*(b2*f+b3*y); a23 = 1/Zp*(c2*f+c3*y); a24 = -x* Sin (w )-(x/f*(x* Cos(k)-y* Sin (k))-f* Sin (k))* Cos (w ); a25 = -f* Cos (k)-y/f*(x * Sin (k)+y* Cos (k)); a26= -x; 3、前方交会： 计算旋转矩阵： a = Cos(ψ) * Cos(k) - Sin(w) * Sin(ψ) * Sin(k); a2 = Cos(ψ) * Sin(k) - Sin(w) *Sin(ψ) * Cos(k); a3 = - Sin(ψ) *Cos(w); b1 = Cos(w) * Sin(k) ; b2 = Cos(w) * Math.Cos(k); b3 = Sin(w); c1 = =Sin(ψ) *.Cos(k) + Cos(ψ) * Sin(w) *Sin(k); c2 =- Sin(ψ) * Sin(k) + Cos(ψ) * Sin(w) * Cos(k); c3 =Cos(ψ) * Cos(w); 基线分量： Bu = Xs2 - Xs1; Bv = Ys2 - Ys1; Bw = Zs2 - Zs1; 投影系数： 计算地面点坐标： 精度评定： 四、实例： 使用数码相机拍摄了两张航测相片Img261和Img262，量测了6对点的影像坐标，其中前4个点为地面控制点（GCP）, 后两个点为加密点。如下表所示。 相片参数：大小为4272（列）*2848（行），像主点位于影像中心，像素大小为5.19micron，相机的焦距为24mm。. Img261的外方位元素初始值：X=397510.760，Y=3445853.978, Z=1455.153，角元素为0； Img262的外方位元素初始值：X=397513.320，Y=3445979.811, Z=1453.685，角元素为0； 请采用光束法区域网空中三角测量方法，解求每张相片的外方位元素和加密点的地面坐标。 点名 Img261 Img262 地面坐标 x（pixel） y（pixel） x（pixel） y（pixel） X(m) Y(m) H(m) 1 3156.625 95.125 3064.875 903.625 397534.502 3446082.7