adjust.bas

来自「测绘工程中的前方交会,主要是摄影测量中运用的较多」· BAS 代码 · 共 44 行

Attribute VB_Name = "Adjust"
Option Base 1
Sub calc()
 Dim i As Integer
 Dim R(2, 2) As Double
 Dim A As Double
 Dim g_X(14) As Double
 Dim g_Y(14) As Double
 Dim g_Z(14) As Double
 Dim L As Double
    g_X_B = Sqr(g_S_AB ^ 2 - g_H_AB ^ 2)
    For i = 1 To g_PotNum
      L = g_X_B * Sin(g_B_H(i)) / Sin(PI - g_A_H(i) - g_B_H(i))

      
      g_X(i) = L * Cos(g_A_H(i))
      g_Y(i) = L * Sin(g_A_H(i))
      g_Z(i) = L * Tan(g_A_V(i)) + g_HA
      Next i
      A = Atn(g_Y(6) - g_Y(1)) / (g_X(6) - g_X(1))
      R(1, 1) = Cos(A)
      R(1, 2) = Sin(A)
      R(2, 1) = -R(1, 2)
      R(2, 2) = R(1, 1)
    For i = 1 To g_PotNum
      g_XC(i) = g_X(i) * R(1, 1) + g_Y(i) * R(1, 2)
      g_YC(i) = g_X(i) * R(2, 1) + g_Y(i) * R(2, 2)
      g_ZC(i) = g_Z(i)
    Next i
'
End Sub
Sub d_r(g_Dir() As Double)
  Dim A As Double
  Dim b As Double
  Dim c As Double
   For i = 1 To g_PotNum
       
        A = Int(g_Dir(i))
        b = Int((g_Dir(i) - A) * 100) / 60#
        c = (100 * g_Dir(i) - Int(100 * g_Dir(i))) * 100 / 3600#
        g_Dir(i) = (A + b + c) * PI / 180#
   Next i
End Sub