王　義

(中鐵二院工程集團有限責(zé)任公司，四川成都　610031)

Research on Autocad Based on Topographic Map Coordinate System Conversion

WANG Yi

基于Autocad地形圖坐標轉(zhuǎn)換系統(tǒng)研究

王義

(中鐵二院工程集團有限責(zé)任公司，四川成都610031)

Research on Autocad Based on Topographic Map Coordinate System Conversion

WANG Yi

摘要介紹同橢球和不同橢球的地形圖坐標轉(zhuǎn)換方法，并基于Autocad平臺，利用ObjectARX工具開發(fā)地形圖坐標轉(zhuǎn)換系統(tǒng)，應(yīng)用于多個高鐵項目，取得了較好的效果。

關(guān)鍵詞高斯投影四參數(shù)七參數(shù)ObjectARX坐標轉(zhuǎn)換

1概述

隨著我國高鐵技術(shù)的迅速發(fā)展，鐵路勘察設(shè)計標準更加規(guī)范，對鐵路設(shè)計各環(huán)節(jié)的精度要求不斷提高。為了減少長距離鐵路全線制圖帶來的投影變形誤差，鐵路勘察設(shè)計中必須針對鐵路全線進行多個投影分帶，以便更好的控制投影變形，提高鐵路勘察設(shè)計精度。伴隨投影分帶的不斷細化，帶來了大量的地形圖坐標轉(zhuǎn)換工作。

詳細介紹兩種地形圖坐標轉(zhuǎn)換方法，一是同橢球下高斯投影正反算坐標換帶法，二是不同坐標系下通過公共點解算四參數(shù)、七參數(shù)轉(zhuǎn)換法。在此基礎(chǔ)上，首次將高斯投影正反算、四參數(shù)、七參數(shù)公式引入AutoCAD平臺，利用ObjectARX工具開發(fā)了地形圖坐標轉(zhuǎn)換軟件，保證了換帶精度，提高了作業(yè)效率。

2數(shù)字地形圖坐標轉(zhuǎn)換方法

常用的地形圖坐標轉(zhuǎn)換方法有兩種，同橢球高斯投影正反算法及不同橢球四參數(shù)、七參數(shù)轉(zhuǎn)換法。

2.1　同橢球高斯投影正反算法

為了限制長度變形，以中央子午線進行分帶，把投影范圍限制在中央子午線東、西兩側(cè)一定的范圍內(nèi)，使得統(tǒng)一的坐標系分割成各帶的獨立坐標系。在工程應(yīng)用中，往往要用到相鄰帶中的點坐標，這樣便涉及到換帶計算。

地形圖轉(zhuǎn)換前后的坐標系如果是基于同橢球的，并且投影方式相同，可直接采用投影換帶計算，投影方式較常用的有高斯投影、UTM投影及蘭勃特投影等。高斯投影正算公式為

(1)

高斯投影反算公式為

(2)

式中ηf、tf分別為按Bf值計算的相應(yīng)量。

a為橢球長半軸，帶高程抵償面時：長半軸=長半軸理論值+高程抵償面高程值。為消除長度變形誤差，通常做法是選擇合適的基準投影面以使綜合變形在容許變形范圍之內(nèi)，這個基準投影面我們稱之為抵償高程面。

b為橢球短半軸，f為橢球扁率，X為子午線弧長，Bf為底點緯度，帶有高程抵償面的解算，在具有高程抵償面進行高斯變換時，應(yīng)該將緯度B的差值代入

ΔB=sin2B(1-e2sin2B)Z0f/((a+Z0)(1-e2))，

式中，Z0為高程抵償面的高程。

2.2　不同橢球四參數(shù)、七參數(shù)轉(zhuǎn)換法

地形圖轉(zhuǎn)換前后的坐標系如果是基于不同的橢球，需要利用公共點在兩個坐標系下不同的坐標值計算四參數(shù)、七參數(shù)，通過這些參數(shù)計算每個實體轉(zhuǎn)換后的坐標，修改實體坐標為轉(zhuǎn)換后坐標，完成不同橢球下地形圖的轉(zhuǎn)換工作。

四參數(shù)轉(zhuǎn)換法公式為

(3)

式(3)中，x2，y2為轉(zhuǎn)換后的平面直角坐標，x1，y1為原坐標系下平面直角坐標，坐標單位為m。x0，y0為平移參數(shù)，α為旋轉(zhuǎn)參數(shù)，m為尺度參數(shù)。解算這4個參數(shù)，需要兩個以上公共點。

七參數(shù)轉(zhuǎn)換法公式為

(4)

3地形圖坐標轉(zhuǎn)轉(zhuǎn)換系統(tǒng)的設(shè)計與實現(xiàn)

為了提高坐標轉(zhuǎn)換自動化程度，同時減少人工轉(zhuǎn)換可能出現(xiàn)的差錯，根據(jù)實際情況，引用不同的轉(zhuǎn)換模型，利用ObjectARX工具，采用點對點轉(zhuǎn)換方法，即遍歷整個AutoCAD文件中的實體，對單個實體定位點坐標進行轉(zhuǎn)換計算，然后修改單個實體定位點的坐標為變換后的坐標，開發(fā)了數(shù)字地形圖坐標轉(zhuǎn)換軟件。該軟件基于AutoCAD平臺底層開發(fā)，程序運行穩(wěn)定，轉(zhuǎn)換速度快，可實現(xiàn)海量數(shù)據(jù)批量轉(zhuǎn)換。

基于相同橢球轉(zhuǎn)換時，可查閱1954北京橢球、1980西安橢球、WGS84橢球、國家2000橢球等通用橢球，或自定義橢球的長、短半軸，橢球扁率參數(shù)，代入公式(1)、(2)，可實現(xiàn)任意中央子午線、任意高程抵償面間的高精度轉(zhuǎn)換。數(shù)字地形圖坐標轉(zhuǎn)換軟件能實現(xiàn)自動轉(zhuǎn)換，開發(fā)環(huán)境采用ObjectARX2008與VS2005，同橢球坐標轉(zhuǎn)換流程如圖1所示。

圖1　同橢球坐標轉(zhuǎn)換流程

基于不同橢球或者任意坐標系的轉(zhuǎn)換，四、七參數(shù)轉(zhuǎn)換法都可以實現(xiàn)，但考慮地形圖轉(zhuǎn)換要保證高程值不改變，四參數(shù)轉(zhuǎn)換法更具有實際意義，帶有高程的坐標文本文件轉(zhuǎn)換可采用七參數(shù)轉(zhuǎn)換法。四參數(shù)轉(zhuǎn)換流程如圖2所示。

圖2　四參數(shù)坐標轉(zhuǎn)換流程

程序同時提供轉(zhuǎn)換精度報告，即每個公共點的改正量，改正量如異常，檢查公共點輸入格式及判斷粗差點，剔除粗差點，再進行參數(shù)計算，完成地形圖坐標系轉(zhuǎn)換。

對圖形文件進行轉(zhuǎn)換時，根據(jù)實體類型的不同，做出相應(yīng)的轉(zhuǎn)化。例如，針對圖案填充實體可轉(zhuǎn)換其邊界，刪除原有圖案填充，依據(jù)新邊界重新填充圖案；坐標范圍較大的塊(Block)及嵌套塊，如果直接轉(zhuǎn)換誤差會比較大，對于這種塊程序采取打散后轉(zhuǎn)換；文字注記、高程點及植被符號以插入點坐標進行轉(zhuǎn)換；圓弧曲線轉(zhuǎn)換起始點、終止點及中心點坐標，然后重新生成圓弧曲線；標注實體按標注線和標注文本分開轉(zhuǎn)換；通過以上手段解決了難以轉(zhuǎn)換的實體，保證換帶無遺漏。

在實際項目中，會同時涉及到同橢球和不同橢球的轉(zhuǎn)換，為了保證轉(zhuǎn)換精度，經(jīng)常綜合運用高斯正反算法及四參數(shù)法。

下面以一個實例來闡述地形圖換帶流程。某項目地形圖為北京54橢球，中央子午線經(jīng)度100.5°，高程抵償面大地高1 020 m；轉(zhuǎn)換到WGS84橢球，中央子午線經(jīng)度101.75°，高程抵償面大地高為1 060 m?？上壤酶咚狗此隳K將地形圖及公共點轉(zhuǎn)換到相同中央子午線，再利用四參數(shù)模塊進行轉(zhuǎn)換，這樣，公共點投影到相同中央子午線，四參數(shù)轉(zhuǎn)換殘差最小，精度最高。最后在同橢球下，利用高斯正算模塊將地形圖轉(zhuǎn)換到對應(yīng)的工程獨立坐標系。圖3為詳細的換帶流程。

圖3　地形圖換帶流程

4轉(zhuǎn)換精度檢核

4.1　同橢球換帶精度檢查

隨機抽取某幅地形圖內(nèi)26個實體坐標，轉(zhuǎn)換前坐標系采用WGS84橢球，中央子午線經(jīng)度為103.25°，高程抵償面大地高為1940 m；轉(zhuǎn)換后坐標系采用WGS84橢球，中央子午線經(jīng)度為106.5°，高程抵償面大地高為1 255 m。通過本軟件計算和西南交大的GPS2005軟件計算，比較計算結(jié)果，其差值如表1。

表1　同橢球換帶殘差比較

4.2　參數(shù)換帶精度檢查

實驗數(shù)據(jù)為機載Lidar點云數(shù)據(jù)，原始坐標系采用Clarke 1880(克拉克)橢球，中央子午線39°，高程抵償面大地高為0 m，投影方式為UTM投影；換帶后的坐標系采用WGS84橢球，中央子午線39°，高程抵償面大地高為0 m，投影方式為高斯投影；經(jīng)分析前后坐標系為不同橢球，只能利用公共點進行轉(zhuǎn)換，通過本軟件選擇四參數(shù)法計算和Cosawin(科傻)軟件選擇赫爾默特法計算，比較計算結(jié)果，其差值如表2。

表2　參數(shù)換帶殘差比較

從表2統(tǒng)計可見，同橢球高斯正反算轉(zhuǎn)換法精度優(yōu)于0.1 mm，不同橢球參數(shù)轉(zhuǎn)換法精度可控，與公共點精度相關(guān)。

通過多個項目驗證得知，高斯投影正反算法，轉(zhuǎn)換公式嚴密，換帶精度較高。精度不受長距離、中央子午線經(jīng)度跨度較大的影響。四參數(shù)、七參數(shù)轉(zhuǎn)換方法，通過公共點計算平移、旋轉(zhuǎn)以及縮放參數(shù)，近似實現(xiàn)兩個坐標系的轉(zhuǎn)換。該方法對于短距離、中央子午線經(jīng)度跨度不大的情況，轉(zhuǎn)換方法嚴密有效。但對于長距離、中央子午線經(jīng)度跨度較大時，轉(zhuǎn)換精度較差。

5結(jié)論

為了解決海量數(shù)字地形圖基于多個投影分帶的換帶問題，探索將大地測量學(xué)數(shù)學(xué)模型引入AutoCAD平臺，結(jié)合ObjectARX工具開發(fā)的地形圖轉(zhuǎn)換軟件，實現(xiàn)了高精度、無遺漏及快速批量轉(zhuǎn)換數(shù)字地形圖功能，經(jīng)過多個高鐵項目實踐驗證，轉(zhuǎn)換精度可靠，具有較強的實用性。

參考文獻

[1]周適.七參數(shù)坐標轉(zhuǎn)換研究及應(yīng)用[J].鐵道勘察，2013(5)

[2]楊鋒，林春峰，程昂.基于ObjectArx地形圖數(shù)據(jù)標準統(tǒng)一軟件的設(shè)計與實現(xiàn)[J].鐵道勘察，2014(5)

[3]梅熙.高速鐵路坐標轉(zhuǎn)換方法探討[J].高速鐵路技術(shù)，2012(8)

[4]管建華.基于ObjectARX的空間數(shù)據(jù)坐標變換及應(yīng)用[J].現(xiàn)代測繪，2007(1)

[5]鄭敏輝,林鐵.數(shù)字地形圖坐標轉(zhuǎn)換工具的設(shè)計與實現(xiàn)[J].測繪技術(shù)裝備，2008(1)

[6]劉秀軍,孫小英,袁月欣.省區(qū)矢量數(shù)字地圖WGS-84坐標系與1980西安坐標系的轉(zhuǎn)換[J].測繪，2011(6)

[7]孔祥元,郭際明,劉宗泉.大地測量學(xué)基礎(chǔ)[M].武漢：武漢大學(xué)出版社，2001

[8]呂忠剛,許世寧.關(guān)于抵償高程面與移動中央子午線最佳選取問題的研究[J].東北測繪，2002(2)

[9]王小華,程傳錄,范宏濤，等.基于AutoCAD.NET API的地形圖坐標轉(zhuǎn)換實現(xiàn)方法研究[J].測繪與空間地理信息，2013(11)

[10]楊本廷.數(shù)字地形圖坐標變換的探討[J].城市勘測，2010(1)

中圖分類號：P284

文獻標識碼：B

文章編號：1672-7479(2015)02-0001-04

作者簡介：王義(1975—)，男，2000年畢業(yè)于武漢測繪科技大學(xué)攝影測量與遙感專業(yè)，工學(xué)學(xué)士，高級工程師。

收稿日期：2014-12-31