《测量坐标系转换代码》由会员分享，可在线阅读，更多相关《测量坐标系转换代码（2页珍藏版）》请在金锄头文库上搜索。
1、测量坐标变换方法1.VB/VC原始代码Const PI As Double = 3.14159265Private Sub Command1_Click()Rem 定义数据类型Dim x1 As Double, x2 As Double, y1 As Double, y2 As Double, z1 As Double, z2 As DoubleDim As Single, k As Single,彷 As SingleDim a1 As Single, a2 As Single, a3 As SingleDim b1 As Single, b2 As Single, b3 As Single
2、Dim c1 As Single, c2 As Single, c3 As SingleRem 从 text 中读取数据x1 = Val(txtx1.Text): y1 = Val(txty1.Text): z1 = Val(txtz1.Text)cd = deg(Val(Txt.Text): k = deg(Val(TxtK.Text):彷=deg(Val(Txt.TeXt) Debug.Pri nt 4 k,Rem 求解 a1, a2, a3 ，b1, b2, b3， c1, c2, c3al = Cos * PI / 180) * Cos(k * PI / 180) - Sin( * P
3、I / 180) * Sin(K * PI / 180) * Sin * PI / 180)a2 = - Cos * PI / 180) * Sin(K * PI / 180) - Sin仲 * PI / 180) * Sin( * PI / 180) * Cos(K * PI / 180)a3 = -Si n * PI / 180) * Cos( * PI / 180)b1 = Cos( * PI / 180) * Sin(K * PI / 180)b2 = Cos( * PI / 180) * Cos(K * PI / 180)b3 = -Si n( * PI / 180)c1 = Sin
4、(彷 * PI / 180) * Cos(K * PI / 180) + Sin( * PI / 180) * Cos * PI / 180) * Sin(K * PI / 180)c2 = - Sin(彷 * PI / 180) * Si n(K * PI / 180) + Sin ( * PI / 180) * Cos * PI / 180) * Cos(K * PI / 180)c3 = Cos( * PI / 180) * Cos * PI / 180)Debug.Print a1, a2, a3Debug.Print b1, b2, b3Debug.Print b1, b1, b3R
5、em 定义数组Dim a(3, 3) As Singlea(1, 1) = a1: a(1, 2) = a2: a(1, 3) = a3a(2, 1) = b1: a(2, 2) = b2: a(2, 3) = b3 a(3, 1) = c1: a(3, 2) = c2: a(3, 3) = c3 Dim x(3, 1) As Doublex(1, 1) = x1x(2, 1) = y1x(3, 1) = z1Rem 求解Dim y(3, 1) As DoubleFor i = 1 To 3y(i, 1) = a(i, 1) * x(1, 1) + a(i, 2) * x(2, 1) + a(
6、i, 3) * x(3, 1) Next ix2 = y(1, 1)y2 = y(2, 1) z2 = y(3, 1)Debug.Print x2, y2, z2Txtx2.Text = Format(x2, 0.000) Txty2.Text = Format(y2, 0.000)Txtz2.Text = Format(z2, 0.000)End SubRem 定义 deg 函数，即度分秒转换为度Private Function deg(a As Double) sign = Sgn(a)a = Abs(a) + 0.0000000001b = Int(a)c = Int(a - b) * 100)d = a - b - c / 100deg = sign * (b + c / 60 + d / 0.36)End Function2.程序运行结果图
