駱宗萬,曹良中,余芳濱,張雨霏

(九江學院 旅游與地理學院,江西 九江 332005)

1 引言

隨著我國與國外數(shù)據(jù)、技術(shù)交流愈來愈密切，國外的很多地理數(shù)據(jù)產(chǎn)品包括遙感影像產(chǎn)品、氣象數(shù)據(jù)以及再分析數(shù)據(jù)正愈來愈多的被國內(nèi)科技人員所采用，而國內(nèi)的一些數(shù)據(jù)產(chǎn)品如氣象站點觀測數(shù)據(jù)、風云系列數(shù)據(jù)也正慢慢走出國門。因國外地理數(shù)據(jù)產(chǎn)品，尤其是歐美國家的數(shù)據(jù)產(chǎn)品坐標系多采用WGS-84坐標系統(tǒng)，而國內(nèi)實際工作中所使用的坐標系統(tǒng)主要有1954北京坐標系、西安80坐標系、CGCS-2000坐標系等[1,2],其中尤以西安80坐標系的歷史數(shù)據(jù)最多，因此研究如何實現(xiàn)WGS-84坐標系與西安80坐標系的坐標轉(zhuǎn)換具有重要的現(xiàn)實意義[3]。

針對該問題，多位學者進行了研究。高寶華等借助三參數(shù)法對站點數(shù)據(jù)進行了坐標轉(zhuǎn)換[4]，趙飛燕等對四參數(shù)坐標轉(zhuǎn)換方法進行了介紹和應用[5]，姚朝龍等將六參數(shù)模型在平面坐標轉(zhuǎn)換中進行了應用[6]，以上所述幾種方法對于坐標轉(zhuǎn)換來說具有使用簡便的特點，但因考慮的轉(zhuǎn)換因子更少，在精度上有所欠缺。陳宇等從坐標系原點向聯(lián)測點重心位置平移的角度出發(fā)，提出了一種改進的布爾莎轉(zhuǎn)換模型[7]，張秋昭等為提高轉(zhuǎn)換參數(shù)精度，削弱方程的病態(tài)問題，提出了站心坐標轉(zhuǎn)換改進模型[8]，這兩種新的模型雖然精度有所提高，但因使用條件和情況的限制，整個求取轉(zhuǎn)換參數(shù)，實現(xiàn)兩種坐標系之間相互轉(zhuǎn)換的過程對于使用者而言復雜程度更大，不具有普適性。

前人的研究主要集中在對各種坐標轉(zhuǎn)換模型精度和簡便性的單一考量上，然而對二者綜合考量與實踐應用方面的研究較少。本文基于python這門簡單易學且不受開發(fā)平臺局限的跨平臺開源性語言[9～11]，以編程的方式實現(xiàn)了布爾莎七參數(shù)轉(zhuǎn)換模型，其過程借助該代碼完成了WGS-84坐標系統(tǒng)與西安80坐標系統(tǒng)轉(zhuǎn)換參數(shù)的計算，并通過ArcGIS軟件，將該參數(shù)在柵格數(shù)據(jù)類型中進行應用。最后通過4個站點數(shù)據(jù)驗證，該轉(zhuǎn)換系數(shù)的數(shù)據(jù)精度滿足工程需要。

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)

WGS-84坐標系是一種地心坐標系，其坐標原點位于地球的質(zhì)心[12]，西安80坐標系屬于參心大地坐標系,它們之間的坐標轉(zhuǎn)換，歸根到底是不同的空間直角坐標系之間的換算[13]，是WGS84橢球體與IAG75橢球體的轉(zhuǎn)換問題,兩種橢球體的坐標參數(shù)如表1所示。

表1 WGS-84與西安 80坐標參數(shù)

本文選用位于新疆北部布爾津縣8個同時具有西安80坐標系和WGS-84坐標系的坐標點。其中SJD2、SJE6、SJE10、SJE11四個站點用于計算坐標轉(zhuǎn)換參數(shù)，另外SJD5、SJE13、SJE14、SJE16四個站點用于檢核七參數(shù)的精度。其分布情況如圖1所示，其坐標信息如表2所示。

圖1 WGS-84坐標系坐標點分布

2.2 方法

布爾莎-沃爾夫模型，也稱默特模型，因在模型使用過程中七個參數(shù)Tx、Ty、Tz、Wx、Wy、Wz、m的部分參數(shù)的求解情況，它可被分為七參數(shù)轉(zhuǎn)換模型、六參數(shù)轉(zhuǎn)換模型、四參數(shù)轉(zhuǎn)換模型和三參數(shù)轉(zhuǎn)換模型[7,14,15]。本文主要研究的就是用以實現(xiàn)WGS-84坐標系到西安80坐標系轉(zhuǎn)換的布爾莎模型七參數(shù)法。兩個不同空間直角坐標的轉(zhuǎn)換模型如式(1)所示：

表2 站點坐標信息

(1)

式(1)中的(Xold,Yold,Zold)表示W(wǎng)GS-84坐標系成果，(Xnew,Ynew,Znew)表示西安80坐標系成果，三個平移參數(shù)Tx、Ty、Tz因兩個坐標系的坐標原點不同產(chǎn)生，三個旋轉(zhuǎn)角參數(shù)Wx、Wy、Wz因各坐標系相對應的坐標軸不平行產(chǎn)生，而尺度參數(shù)m則是因兩個坐標系的尺度不一致所產(chǎn)生。式中共有Tx、Ty、Tz、Wx、Wy、Wz、m七個參數(shù)，被稱為布爾莎七參數(shù)公式。當Wx、Wy、Wz均為小角度時，布爾莎七參數(shù)模型可以簡化為式(2)：

(2)

列出誤差方程式(3)：

(3)

為了檢核坐標轉(zhuǎn)換參數(shù)估計值的可靠性，本文采用外符合精度檢測的方法得到外符合中誤差m外，對參數(shù)估計值進行精度評定[16]。

(4)

其中△為誤差，可由表2中的SJD5、SJE13、SJE14、SJE16四個站點數(shù)據(jù)和參數(shù)估計值代入誤差方程(3)得到，n為測量個數(shù)。

借助于布爾莎模型進行參數(shù)求解及驗證的流程如圖2所示。

圖2 工作流程

其中，用以求取參數(shù)估計值的部分代碼如圖3所示。

圖3 具體實現(xiàn)部分代碼

圖3在求取參數(shù)值的過程中，首先將SJD2、SJE6、SJE10、SJE11四個站點數(shù)據(jù)代入系數(shù)列矩陣中，得到坐標轉(zhuǎn)換的參數(shù)估計值，接著將坐標轉(zhuǎn)換參數(shù)估計值和SJD5、SJE13、SJE14、SJE16四個站點數(shù)據(jù)代入誤差方程中，得到用于計算外符合中誤差做精度驗證的各點中誤差。這些過程與圖2所示流程相一致。

在求取七參數(shù)并進行驗證之后，通過ArcGIS軟件實現(xiàn)WGS84坐標系下的柵格數(shù)據(jù)轉(zhuǎn)換到Xian80坐標系下，其具體流程如圖4所示。

圖4 參數(shù)應用流程

圖4中，首先，借助于ArcGIS中的自定義地理(坐標)變換工具，通過填入按照圖2工作流程得到的參數(shù)估計值，建立起WGS-84坐標系與西安80坐標系坐標轉(zhuǎn)換的轉(zhuǎn)換參數(shù)；其次，使用ArcGIS中的柵格投影工具，將建立的轉(zhuǎn)換參數(shù)寫入，即可完成柵格數(shù)據(jù)由WGS-84坐標系到西安80坐標系的轉(zhuǎn)換。

3 結(jié)果與討論

使用圖3中的代碼，求解得到研究區(qū)內(nèi)WGS-84坐標系到西安80坐標系的坐標轉(zhuǎn)換參數(shù)，其值如表3所示。

表3 坐標轉(zhuǎn)換參數(shù)

表3中的參數(shù)估計值分別為X平移，Y平移，Z平移，X旋轉(zhuǎn)，Y旋轉(zhuǎn)，Z旋轉(zhuǎn)，尺度變化m。

為了驗證七參數(shù)的精度，選取SJD5、SJE13、SJE14、SJE16 站點計算出四個站點的計算坐標和已知坐標的差值ΔX和ΔY，對求解得到的七參數(shù)進行精度驗證，其值如表4所示。結(jié)果表明：在X和在Y方向上的絕對誤差ΔX和ΔY均小于0.05 m,通過將SJD5、SJE13、SJE14、SJE16站點數(shù)據(jù)代入方程式(4)得到的外符合中誤差約為±0.071 m,精度滿足工程測量的精度要求，可以進行參數(shù)應用。

表4 SJE13等4個站點數(shù)據(jù)轉(zhuǎn)換成果

本文借助于ArcGIS軟件，使用表3的參數(shù)，成功將圖2中WGS-84坐標系下的柵格數(shù)據(jù)轉(zhuǎn)換成為圖5中西安80坐標系下的柵格數(shù)據(jù)。

圖5 WGS-84坐標系坐標點分布

4 結(jié)論

本文采用布爾莎七參數(shù)模型實現(xiàn)了局部區(qū)域內(nèi)WGS-84坐標系與西安80坐標系的轉(zhuǎn)換，并將該參數(shù)在柵格數(shù)據(jù)上進行了實例應用，整個求取轉(zhuǎn)換參數(shù)的過程不僅操作簡便，而且坐標轉(zhuǎn)換的結(jié)果經(jīng)精度驗證，誤差較小，可以滿足工程需求。