【正文】 N, Br, Lr, H, X, Y, Z As Double a = WGSa e2 = WGSe2 H = Val() Br = Rad() Lr = Rad() N = a / Sqr(1 e2 * (Sin(Br)) ^ 2) X = (N + H) * Cos(Br) * Cos(Lr) Y = (N + H) * Cos(Br) * Sin(Lr) Z = (N * (1 e2) + H) * Sin(Br) = Format(X, ) = Format(Y, ) = Format(Z, ) End Sub Private Sub Command6_Click() Dim a, e2, e12, b, q, N, Br, Ba, La, Lr, H, X, Y, Z As Double a = WGSa e2 = WGSe2 b = WGSb e12 = WGSe12 X = Val() Y = Val() Z = Val() q = Atn((Z * a) / (Sqr(X ^ 2 + Y ^ 2) * b)) If X 0 Then Lr = Atn(Y / X) ElseIf Y 0 Then Lr = Atn(Y / X) + PI Else Lr = Atn(Y / X) PI End If Br = Atn((Z + e12 * b * (Sin(q)) ^ 3) / (Sqr(X ^ 2 + Y ^ 2) e2 * a * (Cos(q)) ^ 3)) N = a / Sqr(1 e2 * (Sin(Br)) ^ 2) H = Sqr(X ^ 2 + Y ^ 2) / Cos(Br) N Ba = Txt(Br) La = Txt(Lr) = Format(Ba, ) = Format(La, ) = Format(H, ) End Sub Private Sub Command7_Click() = = = = = = End Sub Private Sub Command8_Click() End End Sub Private Sub Form_Load() PI = XAa = 6378140: XAb = : XAc = : XAe2 = : XAe12 = BJa = 6378245: BJb = : BJc = : BJe2 = : BJe12 = WGSa = 6378137: WGSb = : WGSc = : WGSe2 = : WGSe12 = End Sub Function Rad(Txt As String) As Double s = Val(Txt) du = Fix(s) fe = Fix((s du) * 100) mi = Round((s du fe * ) * 10000, 4) Rad = (du + fe / 60 + mi / 3600) * PI / 180 End Function Function Txt(ra) As Double d = ra * 180 / PI du = Fix(d) fe = Fix((d du) * 60) mi = Round(((d du) * 60 fe) * 60) Txt = du + fe * + mi * End Function三維坐標轉換程序設計代碼Private Sub Command1_Click() Dim k, Ex, Ey, Ez, dx, dy, dz Dim x, z, y, Yy, Xx, Zz k = Val() Ex = Val() Ey = Val() Ez = Val() dx = Val() dy = Val() dz = Val() x = Val() y = Val() z = Val() Xx = (k + 1) * (x * Cos(Ey) * Cos(Ez) + y * Cos(Ey) * Sin(Ez) z * Sin(Ey)) + dx Yy = (k + 1) * (x * (Cos(Ex) * Sin(Ez) + Sin(Ex) * Sin(Ey) * Cos(Ez)) + y * (Cos(Ex) * Cos(Ez) + Sin(Ex) * Sin(Ey) * Sin(Ez)) + z * Sin(Ex) * Cos(Ey)) + dy Zz = (k + 1) * (x * (Sin(Ex) * Sin(Ez) + Cos(Ex) * Sin(Ey) * Cos(Ez)) + y * (Sin(Ex) * Cos(Ez) + Cos(Ex) * Sin(Ez) * Sin(Ey)) + z * (Cos(Ex) * Cos(Ey))) + dz = Format(Xx, ) = Format(Yy, ) = Format(Zz, ) End SubPrivate Sub Command2_Click() Dim k, Ex, Ey, Ez, dx, dy, dz Dim x, z, y, Yy, Xx, Zz k = Val() Ex = Val() Ey = Val() Ez = Val() dx = Val() dy = Val() dz = Val() Xx = Val() Yy = Val() Zz = Val() x = ((Xx dx) * Cos(Ey) * Cos(Ez) + (Yy dy) * (Cos(Ex) * Sin(Ez) + Sin(Ex) * Sin(Ey) * Cos(Ez)) + (Zz dz) * (Sin(Ex) * Sin(Ez) + Cos(Ex) * Sin(Ey) * Cos(Ez))) / (k + 1) y = ((Xx dx) * Cos(Ey) * Sin(Ez) + (Yy dy) * (Sin(Ex) * Sin(Ey) * Sin(Ez) + Cos(Ex) * Cos(Ez)) + (Zz dz) * (Sin(Ex) * Cos(Ez) + Cos(Ex) * Sin(Ez) * Sin(Ey))) / (k + 1) z = ((Xx dx) * (Sin(Ey)) + (Yy dy) * Sin(Ex) * Cos(Ey) + (Zz dz) * (Cos(Ex) * Cos(Ey))) / (k + 1) = Format(x, ) = Format(y, ) = Format(z, ) End SubPrivate Sub Command3_Click() = = = = = = = = = = = = = End SubPrivate Sub Command4_Click()EndEnd Sub致謝 本文在王愛生老師的細心指導下完成，從選題和大方向上到修改、定稿都有老師足跡。老師是治學嚴謹，謙虛，不斷追求完善的的精神給我了深刻的教誨，老師的平易近人和穩(wěn)定的心態(tài)讓我明白做事情不要指望一蹴而就，要腳踏實地，務實才是最重要的，這些財富必將影像我以后的日常工作和學習。由衷的感謝老師。同時，在此次畢業(yè)設計過程中我也學到了許多了關于坐標轉換方面的知識，編程技能也有了很大的提高。另外，我還要特別感謝老師對我論文寫作的指導，他們?yōu)槲彝瓿蛇@篇論文提供了巨大的幫助。最后我要向我的父母致謝，感謝他們對我的理解和支持。掩卷而思，唯有感激，唯有自省，唯有不懈地追尋！36