vb坐標(biāo)轉(zhuǎn)換程序設(shè)計(jì)(已修改)
2025-10-18 08:05 本頁(yè)面
【正文】 Option Explicit Dim k2, e2, dX2, dY2 Dim x2, Xx2, y2, Yy2 Dim k3, Ex, Ey, Ez, dX3, dY3, dZ3 Dim X3, Y3, Z3, Xx3, Yy3, Zz3 Const PI = Private Sub Check1_Click() If = 1 Then = 5175 ElseIf = 0 Then = 4440 End If End Sub Private Sub cmdBrowFile_Click() = 控制點(diǎn)文件 (*.gcp)|*.gcp|所有文件 (*.*)|*.* = 1 = End Sub Private Sub cmdCalc_Click() Dim s As String, iPos%, i%, iCent! Dim n%, x1(), y1(), x2(), y2() Dim A() As Double, L() As Double, x(1 To 4) As Double Dim At(), Naa(), W() Open For Input As 1 Line Input 1, s n = Val(s) ReDim x1(n), y1(n), x2(n), y2(n) For i = 1 To n Line Input 1, s iPos = InStr(s, ,) x1(i) = Val(Left(s, iPos 1)) s = Mid(s, iPos + 1) iPos = InStr(s, ,) y1(i) = Val(Left(s, iPos 1)) s = Mid(s, iPos + 1) iPos = InStr(s, ,) x2(i) = Val(Left(s, iPos 1)) s = Mid(s, iPos + 1) y2(i) = Val(s) Next i Close 1 39。計(jì)算轉(zhuǎn)換參數(shù) ReDim A(1 To 2 * n, 1 To 4) As Double, L(1 To 2 * n) As Double ReDim At(1 To 4, 1 To 2 * n), Naa(1 To 4, 1 To 4), W(1 To 4) 系數(shù)矩陣 A： For i = 1 To n A(2 * i 1, 1) = 1: A(2 * i 1, 2) = 0: A(2 * i 1, 3) = x1(i): A(2 * i 1, 4) = y1(i) A(2 * i 1, 1), A(2 * i 1, 2), A(2 * i 1, 3), A(2 * i 1, 4) A(2 * i, 1) = 0: A(2 * i, 2) = 1: A(2 * i, 3) = y1(i): A(2 * i, 4) = x1(i) A(2 * i, 1), A(2 * i, 2), A(2 * i, 3), A(2 * i, 4) L(2 * i 1) = x2(i): L(2 * i) = y2(i) Next i 常數(shù)向量 L： For i = 1 To 2 * n L(i) Next i MatrixTrans A, At A的轉(zhuǎn)置矩陣： ShowMatrix At Matrix_Multy Naa, At, A Naa： ShowMatrix Naa Matrix_Multy W, At, L W： For i = 1 To 4 W(i) Next i MajorInColGuass Naa, W, x X For i = 1 To 4 x(i) Next i 39。分離旋轉(zhuǎn)和尺度參數(shù) If Abs(x(3)) Then If x(4) 0 Then e2 = PI / 2 Else e2 = PI * 3 / 2 End If Else e2 = Atn(x(4) / x(3)) 39。得到的是弧度 If x(3) 0 And x(4) 0 Then e2 = PI e2 ElseIf x(3) 0 And x(4) 0 Then e2 = PI + e2 ElseIf x(3) 0 And x(4) 0 Then e2 = PI * 2 + e2 End If End If k2 = x(3) / Cos(e2) 39。將轉(zhuǎn)換參數(shù)寫(xiě)入相應(yīng)文本框 txtK2 = Str(k2 1) e2 = e2 * 180 / PI Dim du%, fen% du = Int(e2): e2 = (e2 du) * 60 fen = Int(e2): e2 = (e2 fen) * 60 e2 = Val(Format(e2, )) e2 = du + fen / 100 + e2 / 10000 txtE2 = Str(e2) = Str(x(1)) = Str(x(2)) End Sub Private Sub cmdCalc2_Click() k2 = Val() e2 = Val() e2 = DoToHu(e2) dX2 = Val() dY2 = Val() x2 = Val() y2 = Val() Xx2 = (k2 + 1) * (x2 * Cos(e2) + y2 * Sin(e2)) + dX2 Yy2 = (k2 + 1) * (y2 * Cos(e2) x2 * Sin(e2)) + dY2 = Format(Xx2, ) = Format(Yy2, ) End Sub Private Sub cmdCalc3_Click() k3 = Val() Ex = Val() Ex = DoToHu(Ex) Ey = Val() Ey = DoToHu(Ey) Ez = Val() Ez = DoToHu(Ez) dX3 = Val() dY3 = Val() dZ3 = Val() X3 = Val() Y3 = Val() Z3 = Val() Xx3 = (k3 + 1) * (X3 * Cos(Ey) * Cos(Ez) + Y3 * Cos(Ey) * Sin(Ez) Z3 * Sin(Ey)) + dX3 Yy3 = (k3 + 1) * (X3 * (Cos(Ex) * Sin(Ez) + Sin(Ex) * Sin(Ey) * Cos(Ez)) + Y3 * (Cos(Ex) * Cos(Ez) + Sin(Ex) * Sin(Ey) * Sin(Ez)) + Z3 * (Sin(Ex) * Cos(Ey))) + dY3 Zz3 = (k3 + 1) * (X3 * (Sin(Ex) * Sin(Ez) + Cos(Ex) * Sin(Ey) * Cos(Ez)) + Y3 * (Sin(Ex) * Cos(Ez) + Cos(Ex) * Sin(Ey) * Sin(Ez)) + Z3 * (Cos(Ex) * Cos(Ey))) + dZ3 = Format(Xx3, ) = Format(Yy3, ) = Format(Zz3, ) End Sub Private Sub cmdClear2_Click() = = = = End Sub Private Sub cmdClear3_Click() = = = = = = End Sub Private Sub cmdconCalc2_Click() k2 = Val() e2 = Val() e2 = DoToHu(e2) dX2 = Val() dY2 = Val() Xx2 = Val() Yy2 = Val() x2 = ((Xx2 dX2) * Cos(e2) (Yy2 dY2) * Sin(e2)) / (k2 + 1) y2 = ((Yy2
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1