ArcGIS坐標轉換方法及其精度評估

趙慧慧，葛 瑩，肖勝昌，王 沖，楊林波

（1.河海大學 地球科學與工程學院，江蘇 南京 210098；2.中國電建集團昆明勘測設計研究院有限公司 測繪地理信息分院，云南 昆明 650051）

ArcGIS坐標轉換方法及其精度評估

趙慧慧1，葛 瑩1，肖勝昌2，王 沖2，楊林波2

（1.河海大學 地球科學與工程學院，江蘇 南京 210098；2.中國電建集團昆明勘測設計研究院有限公司 測繪地理信息分院，云南 昆明 650051）

ArcGIS提供了靜態轉換、動態轉換和即時轉換3種坐標轉換方法。基于我國1954北京坐標系、1980西安坐標系和2000國家大地坐標系，選定等級較高、分布均勻的坐標成果點，利用靜態轉換、動態轉換和即時轉換進行坐標轉換方法精度分析。

坐標轉換；精度評估；ArcGIS軟件；地理數據庫

在地理信息系統建設與應用中，經常需要進行空間坐標轉換[1-3]。在我國現行的測繪成果中，仍有大量數據采用1980西安坐標系，甚至是1954北京坐標系[4-6]，按照國務院要求，我國將在2016年前完成現行國家大地坐標系向2000國家大地坐標系的過渡[7]。ArcGIS作為主流的地理信息系統平臺，廣泛應用于我國的地理信息數據生產、建庫和應用系統的開發中，形成了大批基于ArcGIS軟件的矢量數據[8,9]。當前關于ArcGIS軟件的坐標轉換研究，主要集中在坐標系討論和坐標轉換程序開發方面，坐標轉換方式對空間數據精度的影響評估卻不多見。本文將針對該問題進行深入探討，在綜述ArcGIS坐標轉換方式的基礎上，針對我國常用的坐標系，進行點位坐標轉換的精度評估與檢核。

1 ArcGIS的坐標轉換

ArcGIS地理參照處理策略是將空間數據和坐標系分離存儲[10]，所以坐標轉換時，不僅要定義地理參照系，還要考慮空間數據坐標處理方式。一般來說，ArcGIS軟件包含2套坐標系統：地理坐標系和投影坐標系[11,12]。前者是用經度和緯度定義球或橢球面上點位的參照系[13]，后者是為二維或三維點、線、面要素的位置定位的（x，y，z）參照系[14]。ArcGIS軟件預置了全世界上百種地理參照系，其中我國常用的地理坐標系有1954北京坐標系、1980西安坐標系和2000國家大地坐標系。

在此基礎上，ArcGIS軟件提供了靜態轉換、動態轉換和即時轉換3種空間坐標轉換方法，見表1。

表1 ArcGIS 3種坐標轉換方式的比較

由此可知，對于ArcGIS坐標轉換，無論是地理參照系定義還是空間坐標轉換方式，都存在較明顯的差異。如果不能深入理解ArcGIS空間坐標轉換的實現機理，空間坐標轉換方式的誤用則不可避免。

2 實驗的設計與組織

2.1 數據準備

本文選定某區域12個均勻分布的坐標成果點，覆蓋范圍約為1 209 km2。按照研究目的，選擇了3套地理坐標系（1954北京坐標系、1980西安坐標系和2000國家大地坐標系），以及基于這3套地理坐標系的高斯-克呂格投影坐標。

2.2 實驗設計和步驟

2.2.1 實驗設計

為了分析ArcGIS靜態轉換、動態轉換和即時轉換3種坐標轉換方法，針對1954北京坐標系、1980西安坐標系和2000國家大地坐標系，本文設計了3組點位坐標轉換實驗：①1954北京坐標系向1980西安坐標系轉換（簡稱Beijing1954 to Xi'an1980）；②1954北京坐標系向2000國家大地坐標系轉換（簡稱Beijing1954 to CGCS2000）；③1980西安坐標系向2000國家大地坐標系轉換（簡稱Xi'an1980 to CGCS2000）。

2.2.2 實驗流程

實驗前，以“Beijing1954 to CGCS2000”為例，說明本文的實驗流程：

1）選擇點位坐標轉換模型。選擇二維七參數轉換模型，利用ArcGIS的Data Management Tools/Projection and Transformation/Create Custom Geographic Transformation工具，定義“Beijing1954 to CGCS2000”坐標轉換模型參數。

2）生成shapefile文件。在ArcGIS中，生成基于Beijing1954的 shapefile文件，將12個坐標成果展繪到該文件。

3）實現靜態轉換。利用ArcGIS的Data Management Tools/Projection and Transformation/Feature/Projection工具，實現靜態轉換。

4）添加點位坐標屬性。利用Data Management Tools/Features下的Add XY Coordinates工具，將坐標字段添加到屬性表，得到轉換后坐標。

5）實現動態轉換。在ArcCatalog環境下，創建Geodatabase數據庫，先以CGCS2000為地理坐標系建立要素數據集，再用ArcGIS Import功能將建立的要素類導入要素數據集，完成動態轉換。

6）實現即時轉換。在ArcMap環境下，先添加CGCS2000控制點坐標文件，再添加Beijing1954的shapefile文件。此時，12個控制點坐標自動從Beijing1954坐標系配準至CGCS2000坐標系，實現坐標系的即時轉換。

7）評估坐標轉換精度。轉換后坐標分別與CGCS2000真實坐標進行比較，對實驗結果進行統計性描述，繪制誤差曲線，分析坐標轉換精度，得出相關結論。

3 實驗結果與分析

本文選取點位中誤差作為精度評價指標，計算公式如下：

3.1 統計描述

本文完成了3組實驗，分別是“Beijing1954 to CGCS2000”、“Xi'an1980 to CGCS2000”以及“Beijing1954 to Xi'an1980”，基本上涵蓋了目前國內測繪生產中常見的坐標轉換。各組實驗結果及其統計描述性質列于表2～4。

表2 Beijing1954 to CGCS2000實驗結果的描述性統計

1）實驗1 ：Beijing1954 to CGCS2000。由表2可知，ArcGIS坐標轉換方式對點位坐標精度會產生較大影響。在ArcGIS軟件中，“Beijing1954 to CGCS2000”靜態轉換最好，動態轉換和即時轉換方式次之，且后兩種方式之間的點位坐標精度差異較小，而靜態轉換與動態轉換、即時轉換之間相差卻較大。因此，在空間坐標轉換時，盡量考慮靜態轉換。

2）實驗2：Xi'an1980 to CGCS2000。由表3可知，從坐標轉換方式上來看，實驗2與實驗1的研究結論相似。但是，對于靜態轉換，“Xi'an1980 to CGCS2000”坐標轉換精度稍高于“Beijing1954 to CGCS2000”；對于動態轉換和即時轉換，“Beijing1954 to CGCS2000”明顯高于“Xi'an1980 to CGCS2000”坐標轉換精度。

3）實驗3：Beijing1954 to Xi'an1980。由表4可知，從不同坐標系間的轉換上來看，實驗3與實驗1有相似結論，即在ArcGIS環境下，“Beijing1954 to Xi'an1980”靜態轉換最好，動態轉換和即時轉換次之，且后兩種方式之間的空間數據精度差異較小，靜態轉換與動態轉換、即時轉換之間相差很大。

表3 Xi'an1980 to CGCS2000實驗結果的描述性統計

表4 Beijing1954 to Xi'an1980實驗結果的描述性統計

總的來說，對于靜態轉換，“Beijing1954 to Xi'an1980”坐標轉換精度最高，“Beijing1954 to CGCS2000”坐標轉換精度最低；而對于動態轉換和即時轉換，“Beijing1954 to CGCS2000”坐標轉換精度最高，“Xi'an1980 to CGCS2000”坐標轉換精度最低。

3.2 實驗精度的分析

從點位中誤差的角度，運用ArcGIS靜態轉換、動態轉換和即時轉換3種方式，對“Beijing1954 to CGCS2000”、“Xi'an1980 to CGCS2000”、“Beijing1954 to Xi'an1980”3種坐標轉換的空間坐標轉換精度進行分析（見表5）。

表5 3種坐標轉換方式平均點位中誤差的比較

由表5可知，就平均中誤差來說，靜態轉換最好，動態轉換和即時轉換次之，且靜態轉換與動態轉換或即時轉換的精度差距較大。對于靜態轉換而言，“Beijing1954 to Xi'an1980”數據精度最高（0.199），其次是“Xi'an1980 to CGCS2000”（1.260），精度最低的是“Beijing1954 to CGCS2000”（1.480）。

進一步地繪制3種坐標轉換的點位中誤差曲線如圖1、圖2所示。

圖1 靜態轉換的點位中誤差曲線圖

由圖1可知，針對靜態轉換方式，“Beijing1954 to Xi'an1980”的精度比“Beijing1954 to CGCS2000”、“Xi'an1980 to CGCS2000”的精度高得多。另外，12個點位中誤差曲線形狀相似，表明誤差變化趨勢一致。

圖2 點位中誤差曲線圖

由圖2可知，無論是動態轉換還是即時轉換，“Beijing1954 to CGCS2000”的精度最高，“Xi'an1980 to CGCS2000”的精度最低。

推測導致ArcGIS空間坐標轉換誤差的原因是：① 地理坐標系的定義完全不同，1954北京坐標系、1980西安坐標系是參心坐標系，而2000國家大地坐標系是地心坐標系；②ArcGIS軟件的空間坐標轉換方法機理不同，但ArcGIS并沒有對外公布這些轉換機理；③由于地理坐標轉換參數具有時空限制性，所以轉換參數解算誤差也是空間坐標轉換誤差來源；④ 對于動態轉換和即時轉換，其空間坐標轉換主要目的是單純滿足制圖可視化需要，對其坐標精度沒有做出明確要求。

4 結 語

ArcGIS空間坐標轉換方式不同，對點位坐標轉換精度的影響也不同；另一方面，在同一種空間坐標轉換方式下，不同坐標系之間的轉換精度也不相同，具體地說：

1）對于靜態轉換、動態轉換和即時轉換而言，無論是“Beijing1954 to Xi'an1980”，還是“Beijing1954 to CGCS2000”、“Xi'an1980 to CGCS2000”，靜態轉換的點位精度較高，動態轉換和即時轉換次之，且后兩者之間的精度相差不大。

2）在靜態轉換中，“Beijing1954 to Xi'an1980”的精度最高，“Beijing1954 to CGCS2000”的精度最低；在動態轉換和即時轉換中，“Beijing1954 to CGCS2000”的精度最高，“Xi'an1980 to CGCS2000”的精度最低。

[1] 李兆雄.大范圍地理空間信息點對點坐標轉換設計與實現[J].地理空間信息,2010,8(2):21-25

[2] 顧秀梅,楊斌,高德政. ArcGIS空間數據投影坐標轉換方法研究[J].測繪科學, 2012,37(4):164-166

[3] 何駿,呂琲,廖明.基于ArcGIS的數據建庫與更新中的坐標轉換技術探討[J].江西測繪,2011(3):15-17

[4] 李平,盧麗. ArcGIS中幾種坐標系轉換方法的應用研究[J].城市勘測,2012(1):87-88

[5] 王寶山,武繼軍,羅建元.大地坐標轉換模型與精度研究[J].河南理工大學學報,2006,25(2):1-5

[6] 廖永生,陳瑞波.從地方坐標系到2000國家大地坐標系的轉換方法[J].海洋測繪,2010,30(5):6-9

[7] 魏子卿.2000中國大地坐標系[J].大地測量與地球動力學,2008(6):1-5

[8] 賈艷紅,劉秀君,郝志強.基于ArcGIS平臺的數據處理與建庫思路探討[J].數字技術與應用,2012(11):209-210

[9] 王小華. ArcGIS平臺矢量數據高精度坐標轉換實現方法[J].測繪與空間地理信息,2012,35(9):90-94

[10] 邢超,李斌.ArcGIS學習指南——ArcToolBox[M].北京:科學出版社,2010

[11] 周衛娟,張成剛. GIS軟件坐標系統剖析[J].現代測繪,2007,30(2):13-16

[12] 黃爽兵,劉昌榮,張成功,等.地質調查中手持GPS和GIS坐標系統應用常見問題[J].四川地質學報,2009(增刊):307-309

[13] 韓晨,程軍,楊秀英.地理坐標系淺析[J].價值工程,2012(32): 318-319

[14] 高德章.大地坐標系與投影坐標系[J].物探化探計算技術, 2011,33(1):51-54

[15] 李恒民,陳永波,李娜,等. 基于ArcGIS Engine的地圖投影與坐標轉換的研究與實現[J].城市勘測,2012(5):105-106

[16] 葉建平. ArcGIS自定義坐標系在城市基礎地理信息系統建庫及地圖融合中的應用[J].測繪通報,2012(12):70-74

[17] 喬金海,許哲,賈士軍. ArcGIS中的坐標系統及投影變換應用研究[J].城市道橋與防洪,2011(12):162-164

[18] 牛亮,杜艷琴.采用ArcGIS進行坐標轉換的兩種方法對比研究[J].科技視野,2014,11(1):349-350

[19] 唐勇.GPS坐標向地方坐標轉換的兩種方法精度比較[J].山西建筑,2009,35(5):353-355

P208

B

1672-4623（2016）03-0073-04

10.3969/j.issn.1672-4623.2016.03.023

趙慧慧,碩士, 研究方向為GIS工程及空間分析應用。

2015-05-07。

項目來源：國家自然科學基金資助項目（41071347）；云南省重大科技專項資助項目（2013ZB006）。