测量平差程序设计5.pptVIP

测量平差程序设计 Case 1 解算近似坐标 d = cds d是测站数 For i = 1 To nl(cds) 依次访问所有的方向值 p = 0 设识别变量 For j = 1 To d 依次访问所有测站 If dm(j) = lb(i) Then p = 1 查看目标点是否设过测站，是则对识别变量p赋值1。 Next j 测量平差程序设计 If p = 0 Then d = d + 1: redim preserve dm(d):dm(d) = lb(i) :ReDim Preserve ns(d): ns(d) = ns(d - 1):如 p=0，表明目标点未设过测站,将该点点名加入点名数组 Next i zds = d 将总点数存入模块级变量zds ReDim x(zds), y(zds) 重新定义坐标数组 x(1) = 100000: y(1) = 100000 为推算近似坐标,对第一个点赋假设坐标值 k = 1 测量平差程序设计 ss = sid(1, k) 调出第一点到未知点方向的边长，参数是测站点序号，照准方向号 If ss=0 then k=k+1:ss=sid(1,k) h = seqn(lb(k)) 查k方向照准点的计算序号 x(h) = x(1) + ss * Cos(0): y(h) = y(1) + ss * Sin(0) 计算第一点上第k方向值的目标点假设坐标 测量平差程序设计 For i = 1 To nl(cds) 遍访所有方向值，将其由角度值转换为弧度值. If l(i) 0.001 Then l(i) = radian(l(i)) 零方向值不参加转换 Next I Do For i = 1 To cds 按测站循环 If x(i) 1 Then 在该测站假设坐标已计算出的前提下，求照准点假设坐标 p1 = 0 测量平差程序设计 For j = nl(i - 1) + 1 To nl(i) 遍访i测站上所有方向值 h = seqn(lb(j)) 查目标点对应的序号 If x(h) 1 Then 目标点坐标已解出 t = azimuth(x(i), y(i), x(h), y(h)): p1 = j 反算方位角并记下方向号 Exit For 退出For－Next循环 End If Next j 测量平差程序设计 For j = nl(i - 1) + 1 To nl(i) 再次遍访i测站上的所有方向值 h = seqn(lb(j)) 查目标点对应的序号 If x(h) 1 Then 该照准方向点坐标未求出 测量平差程序设计 tt = t + l(j) - l(p1) 计算j方向的方位角 ss = sid(i, j) 根据测站号、方向号去查找边长 If ss = 0 Then ss = side(i, h) x(h) = x(i) + ss * Cos(tt): y(h) = y(i) + ss * Sin(tt) 计算近似坐标 End If ‘IF x(h)1 Next j End If’ If x(i) 1 Next i 测量平差程序设计 p = 0 For k = 1 To zds If x(k) 1 Then p = 1 查看是否还有没解算出坐标的点,有则对p赋值，使进入循环再次搜索计算 Next k Loop Until p = 0 测量平差程序设计 For i = 1 To yds 查点名数组中那些点是已知点,用m()数组存其序号 For j = 1 To zds If ym(i) = dm(j) Then m(i) = j Next j Next i 测量平差程序设计 t1 = azimuth(x(m(1)), y(m(1)), x(m(2)), y(m(2))) 利用一对公共点求旋转角 t2 = azimuth(xo(1), yo(1), xo(2), yo(2)) dt = t2 - t1: x1 = x(m(1)): y1 = y(m(1)) 测量平差程序设计 For i = 1 To zds 将假设坐标转换到实际坐标 dx = x(i) - x1: dy = y(i) - y1 x(i) = xo(1