Bácsatyai László: Magyarországi vetületek - NymE GEO portál

208
5.1.2.-1.
__________________________________
Private Sub GPS_WGS()
Static Fi As Double, Lambda As Double, H As Double
Static Y As Double, X As Double, Z As Double
Static n As Double, PE As Double, Teta As Double
Lambda = Atn(Y / X)
PE = Sqr(X ^ 2 + Y ^ 2)
Teta = Atn((Z * a) / (PE * b))
Fi = Atn((Z + E_V_Negyzet * b * Sin(Teta) ^ 3) / (PE - E_Negyzet * a * Cos(Teta) ^ 3))
n = Nagy_N(Fi)
H = PE / Cos(Fi) - n
End Sub
__________________________________
5.2.2.-1.
__________________________________
p = Pontok_Szama
For l = 1 To p
Xi(l) = X_Y_Z(l, 1): Rem Az X vektor
Yi(l) = X_Y_Z(l, 2)
Zi(l) = X_Y_Z(l, 3)
Xiv(l) = U_V_W(l, 1) : Rem Az X’ vektor
Yiv(l) = U_V_W(l, 2)
Ziv(l) = U_V_W(l, 3)
Next l
Rem ********
For k = 1 To 3 * p
For j = 1 To 7
Matrix(k, j) = 0
Next j
Next k
For j = 1 To p
i = 3 * (j - 1)
Matrix(i + 1, 1) = 1: Rem Az A együttható mátrix összeállítása
Matrix(i + 1, 4) = Xi(j)
Matrix(i + 1, 6) = -Zi(j)
Matrix(i + 1, 7) = Yi(j)
Matrix(i + 2, 2) = 1
Matrix(i + 2, 4) = Yi(j)
Matrix(i + 2, 5) = Zi(j)
Matrix(i + 2, 7) = -Xi(j)
Matrix(i + 3, 3) = 1
Matrix(i + 3, 4) = Zi(j)
Matrix(i + 3, 5) = -Yi(j)
Matrix(i + 3, 6) = Xi(j)
Tiszta(i + 1) = Xiv(j) - Xi(j): Rem Az l vektor összeállítása