用 VB 实现测量坐标转换系统

坐标转换系统（VB）

东华理工大学 Theory

北京 54 高斯坐标转换西安 80 高斯坐标转换系统 1.0 版

多点处理结果

核心源代码：

'高斯坐标转换成大地坐标过程

Public Sub GausReverse(a As Double, f As Double, x() As Double, y() As

Double, RB() As Double, RL() As Double, k As Integer)

beita2 = beita0 - 1

beita4 = (15 / 32) * e2 ^ 4 - (175 / 384) * e2 ^ 6 + (3675 / 8192) * e2 ^ 8

beita6 = (-35 / 96) * e2 ^ 6 + (735 / 2048) * e2 ^ 8

beita8 = (315 / 1024) * e2 ^ 8

For i = 1 To k

 B0 = x(i) / (C * beita0)

 a1 = a * Cos(B0) / Sqr(1 - e ^ 2 * (Sin(B0)) ^ 2)

 l0 = y(i) / a1

Do

 fxb = 0

 a3 = (1 / 6) * n * (Cos(B0)) ^ 3 * (1 - t ^ 2 + yita ^ 2)

 a4 = (1 / 24) * n * Sin(B0) * (Cos(B0)) ^ 3 * (5 - t ^ 2 + 9 * yita ^ 2 + 4 *

yita ^ 4)

 a5 = (1 / 120) * n * (Cos(B0)) ^ 5 * (5 - 18 * t ^ 2 + t ^ 4 + 14 * yita ^ 2 -

58 * yita ^ 2 * t ^ 2)

 a6 = (1 / 720) * n * Sin(B0) * (Cos(B0)) ^ 5 * (61 - 58 * t ^ 2 + t ^ 4)

 fxb = fxb + (C * beita6 + C * beita8 * (Cos(B0)) ^ 2) * (Cos(B0)) ^ 2

 fxb = (fxb + C * beita4) * (Cos(B0)) ^ 2