陳志華 徐吉松 張克銘

摘 要：闡述直接導(dǎo)出、同名點坐標(biāo)轉(zhuǎn)換、投影變形處理等基于DWG圖形文件制作復(fù)雜地形區(qū)域kml文件的制作方法，并對三種方法進(jìn)行實例介紹；對直接導(dǎo)出法與投影變形處理法進(jìn)行精度分析比較，對同名點坐標(biāo)轉(zhuǎn)換、投影變形處理三種方法進(jìn)行精度分析比較；提出投影變形處理法對基于DWG圖形文件制作復(fù)雜地形區(qū)域kml文件精度最高的結(jié)論。

關(guān)鍵詞：kml文件；坐標(biāo)轉(zhuǎn)換；投影變形

中圖分類號：TP311 文獻(xiàn)標(biāo)志碼：A 文章編號：2095-2945（2019）02-0131-02

Abstract： This paper expounds the methods of making KML files of complex terrain area based on DWG graphic files， such as direct export， coordinate transformation of same-name points， and processing of projection deformation， and introduces three methods with examples. The precision of the direct derivation method and the projection deformation processing method are analyzed and compared， and the accuracy analysis and comparison of the coordinate transformation of the same name point and the projection deformation processing method are carried out. This paper puts forward the conclusion that the precision of kml file in complex terrain area based on DWG drawing file is the highest by the method of projection deformation processing.

Keywords： kml file； coordinate transformation； projection deformation

1 概述

隨著測繪技術(shù)的飛速發(fā)展，各行各業(yè)對測繪新技術(shù)的應(yīng)用越來越廣泛。kml文件用途廣泛，如可用于三維專題圖表示[1]，也可與GIS結(jié)合用于農(nóng)村新能源規(guī)劃與評價[2]，而基于DWG圖形文件所明確的范圍制作的kml文件可用于無人機(jī)航線規(guī)劃、野外地質(zhì)礦產(chǎn)資源勘探、地質(zhì)勘察、環(huán)境保護(hù)、城鄉(xiāng)規(guī)劃與設(shè)計、施工等。

一方面，DWG圖形文件中所使用的坐標(biāo)系屬平面直角坐標(biāo)系，而且存在多種測量坐標(biāo)系的異同；而GoogleEarth所使用的是WGS-84坐標(biāo)系，屬橢球坐標(biāo)系[3]；顯然，將平面直角坐標(biāo)系的DWG圖形文件表達(dá)于GoogleEarth球面之上將存在著投影變形誤差。另一方面，復(fù)雜地形區(qū)域的DWG圖形文件存在著高差變化較大的情況，當(dāng)選擇的投影面為低海拔面時，則高海拔區(qū)域比低海拔的區(qū)域投影變形大，反正亦然[4]；這時，將兩者表達(dá)于GoogleEarth之上也會存在誤差大小的差異。第三方面，離中央子午線越遠(yuǎn)其投影變形越大；這時，將DWG圖形文件表達(dá)于GoogleEarth之上也會存在誤差大小的差異。

2 方法介紹

基于DWG圖形文件制作kml文件的方法可分為直接導(dǎo)出法、同名點坐標(biāo)轉(zhuǎn)換法、投影變形處理法。

2.1 直接導(dǎo)出法

對于中央子午線和坐標(biāo)系均未改變過的標(biāo)準(zhǔn)DWG圖形文件，可使用南方CASS9.1軟件直接導(dǎo)出kml文件。

2.2 同名點坐標(biāo)轉(zhuǎn)換法

同名點坐標(biāo)轉(zhuǎn)換即將同名點的地心大地坐標(biāo)系成果（BLH）轉(zhuǎn)換成地心空間直角坐標(biāo)系成果。

在GoogleEarth中打開kml文件，選取特征點明顯的同名點提取其大地坐標(biāo)（BLH）并保存為txt文件，使用GPS工具箱5.0將其轉(zhuǎn)換成投影坐標(biāo)（xyz）并整理成南方CASS的格式展點于和DWG圖形文件統(tǒng)一窗口中，以所展點位為實際理想位置，將DWG圖形文件轉(zhuǎn)換；轉(zhuǎn)換步驟：南方CASS9.0-地物編輯-坐標(biāo)轉(zhuǎn)換-轉(zhuǎn)換前、轉(zhuǎn)換后公共點坐標(biāo)拾?。ǜ叱梯斎?）、添加-計算四參數(shù)-圖形-使用四參數(shù)轉(zhuǎn)換即可將圖形轉(zhuǎn)換并導(dǎo)出kml文件。

2.3 投影變形處理法

將地球上的各種地理事物與現(xiàn)象表示到平面上，需要采取一系列數(shù)學(xué)變換，才能使地面上的點與平面上的點建立起對應(yīng)關(guān)系。本文使用ArcMap10.1對DWG圖形文件進(jìn)行投影變形處理，再導(dǎo)出kml文件，方法如下：

第一步，在ArcMap10.1窗口打開dwg數(shù)據(jù)線層（若為范圍），dwg圖右鍵，數(shù)據(jù)-導(dǎo)出shp格式數(shù)據(jù)；定義投影坐標(biāo)系：系統(tǒng)工具箱-Data Managemet tools-投影和變換-定義投影-坐標(biāo)系選投影坐標(biāo)系-高斯投影-北京54/西安80-根據(jù)實際情況選擇中心經(jīng)度；圖層右鍵-屬性-坐標(biāo)系-地理坐標(biāo)系-World-WGS84-shp數(shù)據(jù)右鍵-導(dǎo)出數(shù)據(jù)（勾選數(shù)據(jù)框）-關(guān)閉ArcMap10.1。

第二步，準(zhǔn)備同名點數(shù)據(jù)文件：從Google地球上選擇同名點并保存為excl文件，格式如下：

重新打開ArcMap10.1，加載新導(dǎo)出的勾選數(shù)據(jù)框的shp文件，打開從Google地球上選擇的同名點excl文件，sheet表右鍵-顯示XY數(shù)據(jù)。

第三步，調(diào)出空間校正工具條-編輯器-開始編輯shp-新建位移鏈接工具-選擇轉(zhuǎn)換前和轉(zhuǎn)換后同名點-空間校正設(shè)置校正數(shù)據(jù)（選擇以下圖層中的所有數(shù)據(jù)）-校正方法（選擇仿射或相似等）-校正-保存編輯。

第四步，系統(tǒng)工具箱-Conversion Tools-轉(zhuǎn)為kml-圖層轉(zhuǎn)kml。

3 實例

3.1 直接導(dǎo)出法實例

本文引用實例投影轉(zhuǎn)換參數(shù)為：

投影橢球：北京54-克拉索夫斯基橢球

坐標(biāo)加常數(shù)：X（緯度方向）：0；Y（經(jīng)度方向）：500000

橢球半徑：a 6378245.000000；b 6356863.018773

中央子午線： 102°00′ 00″；帶寬：3度

即可直接導(dǎo)出kml文件。

3.2 同名點坐標(biāo)轉(zhuǎn)換法實例

本文引用實例共選取同名點6個，中心經(jīng)度（中央子午線）為99°00′ 00″， 源橢球系為WGS 84橢球， 目標(biāo)橢球系為WGS 84橢球。轉(zhuǎn)換時大地坐標(biāo)以度分秒輸入，如24°15'9.31"應(yīng)輸為24.150931。轉(zhuǎn)換成果數(shù)據(jù)見表1。

使用轉(zhuǎn)換后圖形即可重新導(dǎo)出kml 文件。

3.3 投影變形處理法實例

本文實例屬復(fù)雜地形區(qū)域，相對高差變化較大，投影變形不一致；因此，將兩者表達(dá)于GoogleEarth之上時存在誤差大小的差異；當(dāng)?shù)秃０蔚膮^(qū)域與GoogleEarth符合較好時，高海拔的區(qū)域與GoogleEarth符合相差很大，反正亦然；當(dāng)使用ArcMap10.1先對DWG圖形文件進(jìn)行投影變形處理再導(dǎo)出kml文件時，則兩者均能。

4 精度分析

4.1 模型精度分析

由于修正參數(shù)法是在GoogleEarth軟件窗口中使用“顯示標(biāo)尺”按鈕在GoogleEarth模型上量取kml文件與GoogleEarth同名點之間的差值作為投影轉(zhuǎn)換參數(shù)進(jìn)行修正，同名點坐標(biāo)轉(zhuǎn)換法所需要的大地坐標(biāo)也是在GoogleEarth模型中獲取， 而GoogleEarth模型本身的精度為±30m，顯然，從模型精度分析得出投影變形處理法的精度高于修正參數(shù)法及同名點坐標(biāo)轉(zhuǎn)換法。

4.2 數(shù)學(xué)精度分析

4.2.1 直接導(dǎo)出法與投影變形處理法數(shù)學(xué)精度分析

本文實例選取分布于圖形區(qū)域內(nèi)高低海拔范圍內(nèi)同名點14個，以GoogleEarth上同名點位為基準(zhǔn)，統(tǒng)計直接導(dǎo)出法與投影變形處理法平均點位誤差，繪制數(shù)學(xué)精度柱狀圖，見圖1。

從圖中比較得出，直接導(dǎo)出法平均點位誤差47m，投影變形處理法平均點位誤差12m；因此，投影變形處理法精度高于直接導(dǎo)出法。

4.2.2 修正參數(shù)法、同名點坐標(biāo)轉(zhuǎn)換法、投影變形處理法數(shù)學(xué)精度分析

本文實例選取分布于圖形區(qū)域內(nèi)高低海拔范圍內(nèi)同名點15個，以GoogleEarth上同名點位為基準(zhǔn)，統(tǒng)計修改參數(shù)后點位誤差、同名點坐標(biāo)轉(zhuǎn)換后點位誤差、投影變形處理后點位誤差，繪制數(shù)學(xué)精度柱狀圖，見圖2。

從圖中比較得出，修正參數(shù)法平均點位誤差35m，同名點坐標(biāo)轉(zhuǎn)換法20m，投影變形處理法10m；因此，投影變形處理法精度最高，同名點坐標(biāo)轉(zhuǎn)換法精度次之，修正參數(shù)法精度最低；但是，投影變形處理法的方法步驟相對麻煩一些，同名點坐標(biāo)轉(zhuǎn)換法次之，修正參數(shù)法相對簡單一些。

5 結(jié)束語

基于DWG圖形文件制作kml文件的方法有多種，應(yīng)根據(jù)其坐標(biāo)系、地理位置、地形復(fù)雜情況，以及對kml文件使用的精度要求，采用與之相適應(yīng)的制作方法。在DWG圖形文件與其比例尺精度相一致的情況下，根據(jù)本文分析得出，對于復(fù)雜地形區(qū)域，先使用ArcMap10.1對DWG圖形文件進(jìn)行投影變形處理再導(dǎo)出kml文件的制作方法較為理想，表現(xiàn)為kml文件與GoogleEarth符合較好，精度最高。

參考文獻(xiàn)：

[1]閆鶴，張崢，楊森，等.基于KML的三維專題圖表示方法[J].測繪工程，2012，21（2）：50-53.

[2]杜英俊，于重重，劉杰.基于KML開發(fā)的GIS系統(tǒng)研究與應(yīng)用[J].計算機(jī)應(yīng)用與軟件，2010，27（10）：49.

[3]國家測繪局職業(yè)技能鑒定指導(dǎo)中心.測繪案例分析[M].北京：測繪出版社，2009：21-22.

[4]中華人民共和國國家標(biāo)準(zhǔn).GB 50026-2007.工程測量規(guī)范[S].北京：中國計劃出版社，2008：19.