高斯投影坐标正反算VB程序.pdfVIP

高斯投影坐标正反算VB程序 高斯投影坐标正反算 高斯投影坐标正反算 学 院： 班 级： 学 号： 姓 名： 课程名称 : 指导老师： - 1 - 高斯投影坐标正反算 实验目的： 1.了解高斯投影坐标正反算的基本思想； 2.学会编写高斯正反算程序，加深了解。 实验原理： 高斯投影正算公式中应满足的三个条件： 1. 中央子午线投影后为直线； 2. 中央子午线投影后长度不变； 3. 投影具有正形性质，即正形投影条件。 高斯投影反算公式中应满足的三个条件： 1. x 坐标轴投影成中央子午线，是投影的对称轴； 2. x 轴上的长度投影保持不变； 3. 正形投影条件，即高斯面上的角度投影到椭球面上后 角度没有变形，仍然相等。 操作工具： 计算机中的 VB6.0 - 2 - 高斯投影坐标正反算 - 1 - 高斯投影坐标正反算 - 1 - 高斯投影坐标正反算 a6 / 6 + Sin(8 * b_) * a8 / 8 t = Tan(b_) e1 = Sqr((a * a - b * b) / (b * b)) eta = Sqr(e1 * e1 * Cos(b) * Cos(b)) W = Sqr(1 - e * e * Sin(b_) * Sin(b_)) N = a / x = x_ + N * Sin(b_) * Cos(b_) * l * l / 2 + N * Sin(b_) * Cos(b_) ^ 3 * (5 - t * t + 9 * eta * eta + 4 * eta ^ 4) * l ^ 4 / 24 + N * Sin(b_) * Cos(b_) ^ 5 * (61 - 58 * t * t + t ^ 4) * l ^ 6 / 720 y = N * Cos(b_) * l + N * Cos(b_) ^ 3 * (1 - t * t + eta * eta) * l * l * l / 6 + N * Cos(b_) ^ 5 * (5 - 18 * t * t + t ^ 4 + 14 * eta * eta - 58 * eta * eta * t * t) * l ^ 5 / 120 Text7 = x If dh = 6 Then y_ = y + 500000 + 1000000 * ik1 Else If dh = 3 Then y_ = y + 500000 + 1000000 * ik2 Else MsgBox error, 48, error: Exit Sub End If End If - 1 - 高斯投影坐标正反算 Text8 = y_ End Sub Private Sub Command2_Click() Dim bf#, j%, Wf#, Vf#, Nf#, Mf#, c#, tf#, etaf#, dh%, ik% x = Val(Text7.Text) y_ = Val(Text8.Text) e = Sqr((a * a - b * b) / (a * a)) m0 = a * (1 - e * e) m2 = e * e * m0 * 3 / 2 m4 = e * e * m2 * 5 / 4 m6 = m4 * e * e * 7 / 6 m8 = e * e * m6 * 9 / 8