高斯正反算VB实验代码.docx

1、高斯正反算VB实验代码高斯正反算实验报告一、实验目的：通过对高斯正反算的软件编程，增强学生对大地测量学基础课程的理解，使学生牢固掌握高斯正反算的基本原理和公式，熟悉高斯正反算的基本技能和计算方法。培养学生利用计算机编程的能力、综合运用知识的能力及理论联系实际的能力二、实验要求：要求学生综合运用测绘知识、大地测量知识、数学知识和计算机知识，设计数学模型和程序算法，编制程序实现数据的自动化处理。三、实验环境Visual studio 2008界面如下：输入数据计算: 代码如下：Public Class 高斯投影正反算 Private Sub 高斯投影正反算_Load(ByVal sender As

2、 System.Object, ByVal e As System.EventArgs) Handles MyBase.Load 自动化控件显示初始值 RadioButton2.Checked = True RadioButton3.Checked = True TextBox1.Focus() TextBox1.Text = 0 TextBox2.Text = 0 TextBox3.Text = 0 TextBox4.Text = 0 TextBox5.Text = 0 TextBox6.Text = 0 RichTextBox2.Text = 说明：输入的坐标需为按6带投影且采用克氏椭球参

3、数所得的国家统一坐标 End Sub Private Sub Button1_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button1.Click 获取输入数据 Dim bb1 As Double, bb2 As Double, bb3 As Double, ll1 As Double, ll2 As Double, ll3 As Double bb1 = TextBox1.Text bb2 = TextBox2.Text bb3 = TextBox3.Text ll1 = TextBox

4、4.Text ll2 = TextBox5.Text ll3 = TextBox6.Text Dim B As Double, L As Double, x As Double, y As Double, X1 As Double, Y1 As Double X1 Y1 Dim N As Integer Dim L0 As Double, l1 As Double Dim p As Double p = 206264.80625 检查输入格式的正确性 If (bb1 = 0 And bb1 = 0 And bb2 = 0 And bb3 = 0 And ll1 = 0 And ll2 = 0

5、And ll3 1 Then Y1 = N * i * 10 + y11 Y坐标加代号 i = i * 10 End If Dim tuoqiu As String tuoqiu = 采用克氏椭球参数, RichTextBox1.Text = 按经差 + tyd.ToString(0) + 进行投影分带， + tuoqiu + 其计算结果为： + vbCrLf + x= + x.ToString(0.00000000) + vbCrLf + y= + y.ToString(0.00000000) + vbCrLf + 国家统一坐标为： + vbCrLf + X= + X1.ToString(0

6、.00000000) + vbCrLf + Y= + N.ToString + Y1.ToString(0.00000000) End If 采用1975国际椭球参数的计算公式() If (ty = 2) Then Dim c As Double Dim c1 As Double Dim c2 As Double Dim l2 As Double Dim n1 As Double Dim a0 As Double Dim a3 As Double Dim a4 As Double Dim a5 As Double Dim a6 As Double Dim i As Integer i = 1

7、c = Math.Cos(b1) * Math.Cos(b1) c1 = Math.Sin(b1) * Math.Cos(b1) c2 = Math.Cos(b1) l2 = Math.Pow(l1, 2) n1 = 6399596.652 - (21565.045 - (108.996 - 0.603 * c) * c) * c a0 = 32144.5189 - (135.3646 - (0.7034 - 0.0041 * c) * c) * c a4 = (0.25 + 0.00253 * c) * c - 0.04167 a6 = (0.167 * c - 0.083) * c a3

8、= (0.3333333 + 0.001123 * c) * c - 0.1666667 a5 = 0.00878 - (0.1702 - 0.20382 * c) * c x = 6367452.1328 * b1 - (a0 - (0.5 + (a4 + a6 * l2) * l2) * l2 * n1) * c1 y = (1 + (a3 + a5 * l2) * l2) * l1 * n1 * c2 X1 = x Dim y11 As Double y11 = y + 500000.0 If y11 / i 1 Then Y1 = N * i * 10 + y11 Y坐标加代号 i =

9、 i * 10 End If Dim tuoqiu As String tuoqiu = 采用1975国际椭球参数, RichTextBox1.Text = 按经差 + tyd.ToString(0) + 进行投影分带， + tuoqiu + 其计算结果为： + vbCrLf + x= + x.ToString(0.00000000) + vbCrLf + y= + y.ToString(0.00000000) + vbCrLf + 国家统一坐标为： + vbCrLf + X= + X1.ToString(0.00000000) + vbCrLf + Y= + N.ToString + Y1.

10、ToString(0.00000000) End If End Sub Private Sub Button2_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button2.Click Dim y As Double, X1 As Double, Y1 As Double, B As Double, L As Double X1 = TextBox7.Text 获取输入数据 Y1 = TextBox8.Text Dim L0 As Double Dim i As Integer i = 1 D

11、o While Y1 / i = 10 y = Y1 - Int(Y1 / i) * i - 500000 L0 = 6 * Int(Y1 / i) - 3 i = i * 10 Loop 对Y坐标处理并求出中央子午线经度 按6带克氏椭球反算 Dim bt As Double, BT1 As Double, c3 As Double, c4 As Double, Bf As Double, c5 As Double, c6 As Double, Nf As Double, Z As Double Dim b2 As Double, b3 As Double, b4 As Double, b5

12、As Double, z2 As Double Dim p As Double p = 206264.80625 bt = X1 / 6367558.4969 * p BT1 = X1 / 6367558.4969 c3 = Math.Cos(BT1) * Math.Cos(BT1) c4 = Math.Sin(BT1) * Math.Cos(BT1) Bf = (bt + (50221746 + (293622 + (2350 + 22 * c3) * c3) * c3) * c4 * Math.Pow(10, -10) * p) / p c5 = Math.Pow(Math.Cos(Bf)

13、, 2) c6 = Math.Sin(Bf) * Math.Cos(Bf) Nf = 6399698.902 - (21562.267 - (108.973 - 0.612 * c5) * c5) * c5 Z = y / (Nf * Math.Cos(Bf) b2 = (0.5 + 0.003369 * c5) * c6 b3 = 0.333333 - (0.166667 - 0.001123 * c5) * c5 b4 = 0.25 + (0.16161 + 0.00562 * c5) * c5 b5 = 0.2 - (0.1667 - 0.0088 * c5) * c5 z2 = Mat

14、h.Pow(Z, 2) B = (Bf * p - (1 - (b4 - 0.12 * z2) * z2) * z2 * b2 * p) / 3600.0 L = L0 + (1 - (b3 - b5 * z2) * z2) * Z * p) / 3600.0 RichTextBox2.Text = 按6带克氏椭球反算后，结果为： + vbCrLf + B= + B.ToString(0.0000000000000) + vbCrLf + L= + L.ToString(0.0000000000000) End Sub Private Sub Button3_Click(ByVal sende

15、r As System.Object, ByVal e As System.EventArgs) Handles Button3.Click TextBox1.Focus() TextBox1.Text = 0 TextBox2.Text = 0 TextBox3.Text = 0 TextBox4.Text = 0 TextBox5.Text = 0 TextBox6.Text = 0 RichTextBox1.Text = End Sub Private Sub Button4_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button4.Click TextBox7.Text = TextBox8.Text = RichTextBox2.Text = 说明：输入的坐标需为按6带投影且采用克氏椭球参数所得的国家统一坐标 End Sub Private Sub Button5_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button5.Click End End SubEnd Class