【正文】 60) / 100S = (Angular D M * 100 / 60) * 3600 / 10000Angular_DMS = D + M + SEnd Function20計算計算卯酉圈半徑 N 如下：Public Function GetN(a As Double, b As Double, RB As Double) As Double 39。計算卯酉圈半徑 NDim e2 As Doublee2 = GetE2(a, b)GetN = a / Sqr(1 e2 * (Sin(RB) ^ 2))End Function計算 公式如下：2ePublic Function GetE2(a As Double, b As Double) As Double 39。計算e^2GetE2 = (a ^ 2 b ^ 2) / a ^ 2End Function在此，最重要的就是坐標轉換的公式了，代碼如下：空間大地坐標系轉換空間直角坐標系(a 為地球長半軸,b 為地球短半軸)Public Sub BLH_XYZ(RBLH As BLH, a As Double, b As Double) 39?？臻g大地坐標系轉換空間直角坐標系(a 為地球長半軸,b 為地球短半軸)Dim e2 As Double, N As Double39。Dim BtoX As New XYZe2 = GetE2(a, b)N = GetN(a, b, ) = (N + ) * Cos() * Cos() = (N + ) * Cos() * Sin() = (N * (1 e2) + ) * Sin()End Sub關于空間直角坐標系轉換空間大地坐標系，我們才用的是迭代法，代碼如下：Public Sub XYZ_BLH(IXYZ As XYZ, a As Double, b As Double) 39。空間直角坐標系轉換空間大地坐標系(迭代法)Dim e2 As Double, N As Double21Dim tB As Double, tH As Double, Lb As Double, lH As Doublee2 = GetE2(a, b)If = 0 Then = PI / 2Else = Atn( / )End IfIf = 0 And = 0 Then = PI / 2N = GetN(a, b, ) = / Sin() N * (1 e2)ElsetB = Atn( / Sqr( ^ 2 + ^ 2))mGoTo:N = GetN(a, b, tB)tH = Sqr( ^ 2 + ^ 2) / Cos(tB) NLb = Atn(( + N * e2 * Sin(tB)) / Sqr( ^ 2 + ^ 2))If Abs(Lb tB) 1E99 ThentB = LbGoTo mGoToElse = LbN = GetN(a, b, ) = Sqr( ^ 2 + ^ 2) / Cos() NEnd IfEnd IfEnd Sub根據上述的設計語句，我們可以通過調用函數的方法實現坐標轉換的操作，代22碼如下：Public BtoX As New XYZPublic XtoB As New BLHPublic pi As DoublePublic Sub Command1_Click()Dim a As Double, b As Doublea = 6378245b = 6354950pi = Dim pBLH As New BLHDim pXYZ As New XYZIf = True Then = Val() = Val() = Val()Call BLH_XYZ(pBLH, a, b) = = = = Str() = Str() = Str()End IfIf = True Then = Val() = Val() = Val()Call XYZ_BLH(BtoX, a, b)Dim w As Double23 = Radian_Angular() = Angular_DMS() = Radian_Angular() = Angular_DMS() = = = = Format(, .00000000 ) = Format(, .00000000 ) = Format(, .00000000 ) = Str() = Str() = Str()End IfEnd Sub程序設計好之后就是要檢驗坐標轉換是否成功了，我們找來了一個已知點。在空間直角坐標系統下的坐標值：X=Y=Z=200在空間大地坐標系統的坐標值為：B=000:L=009:H=在文本框輸入已知點的坐標，檢測該點在另外一個坐標系的坐標是否正確：首先從空間直角坐標系統轉換到空間大地坐標系統，結果如下（圖 43）：24圖 43基本上滿足該點在空間大地坐標系統中的坐標值。25第五章 結論工程施工過程中，常常會遇到不同坐標系統間，坐標轉換的問題。目前國內常見的轉換有以下幾種：1，大地坐標（BLH）對平面直角坐標（XYZ）；2，北京 54 全國 80 及 WGS84 坐標系的相互轉換；3，任意兩空間坐標系的轉換。其中第 2 類可歸入第三類中。所謂坐標轉換的過程就是轉換參數的求解過程。常用的方法有三參數法、四參數法和七參數法。以下對上述三種情況作詳細描述如下： 大地坐標（BLH）對平面直角坐標（XYZ）常規(guī)的轉換應先確定轉換參數，即橢球參數、分帶標準（3 度，6 度）和中央子午線的經度。橢球參數就是指平面直角坐標系采用什么樣的橢球基準，對應有不同的長短軸及扁率。一般的工程中 3 度帶應用較為廣泛。對于中央子午線的確定有兩種方法，一是取平面直角坐標系中 Y 坐標的前兩位*3，即可得到對應的中央子午線的經度。如 x=3250212m，y=395121123m，則中央子午線的經度=39*3=117 度。另一種方法是根據大地坐標經度，如果經度是在 ~度之間，那么對應的中央子午線的經度=（+）/2=117 度，其他情況可以據此 3 度類推。另外一些工程采用自身特殊的分帶標準，則對應的參數確定不在上述之列。確定參數之后，可以用軟件進行轉換，以下提供坐標轉換的程序下載。 北京 54 全國 80 及 WGS84 坐標系的相互轉換這三個坐標系統是當前國內較為常用的，它們均采用不同的橢球基準。其中北京 54 坐標系，屬三心坐標系，大地原點在蘇聯的普而科沃，長軸6378245m，短軸 6356863，扁率 1/；西安 80 坐標系，屬三心坐標系，大地原點在陜西省徑陽縣永樂鎮(zhèn)，長軸 6378140m，短軸 6356755，扁率261/；WGS84 坐標系，長軸 ，短軸 ，扁率 1/。由于采用的橢球基準不一樣，并且由于投影的局限性，使的全國各地并不存在一至的轉換參數。對于這種轉換由于量較大，有條件的話，一般都采用 GPS 聯測已知點，應用 GPS 軟件自動完成坐標的轉換。當然若條件不許可，且有足夠的重合點，也可以進行人工解算。詳細方法見第三類。 任意兩空間坐標系的轉換由于測量坐標系和施工坐標系采用不同的標準，要進行精確轉換，必須知道至少 3 個重合點（即為在兩坐標系中坐標均為已知的點。采用布爾莎模型進行求解。27參考文獻［1］張坤宜，2022..交通土木工程測量（第三版） ，武漢：華中科技大學出版社［2］王勇， 測量技術，北京：中國電力出版社［3］王儂、過靜珺，（第二版） ，北京：清華大學出版社［4］陳英， Basic ，大連：大連理工大學出版社28致謝大學三年的生活、學習即將結束，在這里對在這三年中關心、幫助過我的老師、同學致以衷心的感謝。首先，感謝系里的各位領導、各位老師，他們?yōu)榱宋覀冞@些來自全國各地的同學在學校里安心學習、生活，做了大量的工作，正是由于他們面面俱到的工作，才使我順利完成學業(yè)，走向社會。還要感謝尊敬老師，本文是在老師三個月來的指導、幫助下完成的。感謝老師在百忙之中一絲不茍、不厭其煩的孜孜教誨，感謝老師在我懈怠之時給予我的鞭策和激勵。老師優(yōu)秀的道德品質、嚴謹的治學態(tài)度、一絲不茍的科研精神、淵博的知識儲備將是我終生學習的榜樣。感謝測繪工程專業(yè)同屆的所有同學，他們曾給過我許許多多的幫助，和他們相處的日子令我終生難忘。感謝善良而偉大的父親、母親給予我的一切。在將來的生活工作中，我將會以他們?yōu)榘駱优ぷ?，誠實做人。