張紅梅，張超，馬祥

(蕪湖市勘察測(cè)繪設(shè)計(jì)研究院有限責(zé)任公司，安徽蕪湖 241000)

1 引言

經(jīng)國(guó)務(wù)院批準(zhǔn)，根據(jù)《中華人民共和國(guó)測(cè)繪法》，我國(guó)于2008年7月1日，啟用2000國(guó)家大地坐標(biāo)系(CGCS2000)。CGCS2000與現(xiàn)行國(guó)家大地坐標(biāo)系(1954年北京坐標(biāo)系，1980西安坐標(biāo)系)轉(zhuǎn)換、銜接的過(guò)渡期為8年～10年。2008年7月1日后新生產(chǎn)的各類測(cè)繪成果應(yīng)采用2000國(guó)家大地坐標(biāo)系?，F(xiàn)有地理信息系統(tǒng)，在過(guò)渡期內(nèi)應(yīng)逐步轉(zhuǎn)換到2000國(guó)家大地坐標(biāo)系;2008年7月1日后新建設(shè)的地理信息系統(tǒng)應(yīng)采用2000國(guó)家大地坐標(biāo)系。

針對(duì)上述形勢(shì)，通過(guò)對(duì)蕪湖市現(xiàn)有測(cè)繪成果的綜合分析，結(jié)合數(shù)據(jù)內(nèi)容、格式等，對(duì)關(guān)鍵技術(shù)進(jìn)行研究，建立完整的工藝，研究適合蕪湖測(cè)繪成果轉(zhuǎn)換方法，并利用EPS地理信息工作站坐標(biāo)系統(tǒng)轉(zhuǎn)換參數(shù)計(jì)算以及腳本處理模塊實(shí)現(xiàn)基礎(chǔ)測(cè)繪成果數(shù)據(jù)坐標(biāo)系的整體轉(zhuǎn)換。

2 數(shù)據(jù)轉(zhuǎn)換流程

2.1 中央子午線選取

根據(jù)城市區(qū)域范圍選取適合的中央子午線，保證每千米長(zhǎng)度投影變形不超過(guò) 2.5 cm。

2.2 重合點(diǎn)選取

當(dāng)兩種不同坐標(biāo)系進(jìn)行轉(zhuǎn)換時(shí)，坐標(biāo)轉(zhuǎn)換的精度除取決于坐標(biāo)轉(zhuǎn)換的數(shù)學(xué)模型和求解轉(zhuǎn)換參數(shù)的公共點(diǎn)(重合點(diǎn))坐標(biāo)精度外，還和公共點(diǎn)(重合點(diǎn))的多少、幾何形狀結(jié)構(gòu)有關(guān)。

重合點(diǎn)選取的原則是盡量選取足夠的高等級(jí)、高精度且分布均勻的點(diǎn)作為坐標(biāo)轉(zhuǎn)換的重合點(diǎn)。重合點(diǎn)的選取主要是根據(jù)外業(yè)GPS網(wǎng)施測(cè)技術(shù)總結(jié)、點(diǎn)之記與同坐標(biāo)差比較綜合確定重合點(diǎn)。為提高坐標(biāo)轉(zhuǎn)換精度，必須使重合點(diǎn)的數(shù)量足夠且分布均勻、合理，同時(shí)還要經(jīng)過(guò)大量的試算與分析剔除一些變化大的重合點(diǎn)(粗差點(diǎn))，從而確定出最終用于求取轉(zhuǎn)換參數(shù)的重合點(diǎn)。

2.3 轉(zhuǎn)換參數(shù)計(jì)算

常用的坐標(biāo)轉(zhuǎn)換模型有二維七參數(shù)轉(zhuǎn)換模型、平面四參數(shù)轉(zhuǎn)換模型、三維七參數(shù)坐標(biāo)轉(zhuǎn)換模型等，按照《現(xiàn)有測(cè)繪成果轉(zhuǎn)換到2000國(guó)家大地坐標(biāo)系技術(shù)指南》模型選擇中指出，省級(jí)以下的坐標(biāo)轉(zhuǎn)換可選擇三維四參數(shù)模型或平面四參數(shù)模型。對(duì)于相對(duì)獨(dú)立的平面坐標(biāo)系統(tǒng)與2000國(guó)家大地坐標(biāo)系的聯(lián)系可采用平面四參數(shù)模型:

其中，2個(gè)平移參數(shù)x0、y0，1個(gè)旋轉(zhuǎn)參數(shù)α和1個(gè)尺度因子m。

利用平面四參數(shù)轉(zhuǎn)換模型求解轉(zhuǎn)換參數(shù)。

2.4 數(shù)據(jù)成果轉(zhuǎn)換

利用EPS地理信息工作站進(jìn)行數(shù)據(jù)成果轉(zhuǎn)換，EPS地理信息工作站是以數(shù)據(jù)庫(kù)為核心的一體化解決數(shù)據(jù)采集、數(shù)據(jù)編輯、數(shù)據(jù)監(jiān)理、跨平臺(tái)數(shù)據(jù)轉(zhuǎn)換、整合共享與多格式數(shù)據(jù)分發(fā)、數(shù)據(jù)庫(kù)建設(shè)、數(shù)據(jù)更新到GIS分析應(yīng)用等諸多問(wèn)題的綜合信息系統(tǒng)體系。

3 具體做法

3.1 蕪湖市現(xiàn)有數(shù)據(jù)情況

目前有采用蕪湖獨(dú)立坐標(biāo)系制作的1∶500地形圖、1∶1 000地形圖，其中2005年的數(shù)據(jù)為*.DWG格式，2005年后直接利用EPS地理信息工作站地形圖的數(shù)據(jù)格式為*.EDB。*.EDB格式的數(shù)據(jù)可直接進(jìn)行轉(zhuǎn)換，*.DWG數(shù)據(jù)可通過(guò)EPS地理信息工作站調(diào)入轉(zhuǎn)換為*.EDB后進(jìn)行轉(zhuǎn)換。

3.2 蕪湖獨(dú)立坐標(biāo)系統(tǒng)現(xiàn)狀

蕪湖獨(dú)立坐標(biāo)系是以1954年北京坐標(biāo)系為基礎(chǔ)建立的獨(dú)立坐標(biāo)系，投影面高度為橢球面，采用高斯－克呂格投影，該系統(tǒng)不利于GPS測(cè)量快速、方便獲取高精度成果，不利于提升城市坐標(biāo)成果的精度，不利于國(guó)家與地方資料共享。另外，國(guó)家規(guī)定CGCS 2000啟用8年～10年后，1954年北京坐標(biāo)系將停止使用。當(dāng)國(guó)家參心坐標(biāo)系統(tǒng)的消失，該系統(tǒng)也隨之與國(guó)家坐標(biāo)系失去直接、嚴(yán)密的聯(lián)系。

3.3 蕪湖獨(dú)立坐標(biāo)轉(zhuǎn)CGCS 2000重合點(diǎn)的選取

根據(jù)蕪湖市似大地水準(zhǔn)面精化GPS點(diǎn)位的分布情況及已有資料，經(jīng)過(guò)綜合分析，選取重合點(diǎn)11點(diǎn)用于坐標(biāo)轉(zhuǎn)換參數(shù)的計(jì)算，重合點(diǎn)CGCS 2000成果采用2011年蕪湖市似大地水準(zhǔn)面精化GPS成果，蕪湖獨(dú)立坐標(biāo)系成果采用蕪湖市勘察測(cè)繪設(shè)計(jì)研究院有限責(zé)任公司提供的2007年蕪湖市三等GPS擴(kuò)網(wǎng)平差成果及CORS站成果。

在進(jìn)行蕪湖獨(dú)立坐標(biāo)系到CGCS 2000轉(zhuǎn)換時(shí)，經(jīng)過(guò)反復(fù)試算，最終采用重合點(diǎn)11點(diǎn)計(jì)算坐標(biāo)轉(zhuǎn)換參數(shù)。重合點(diǎn)分布圖如圖1所示:

圖1 重合點(diǎn)及GPS點(diǎn)點(diǎn)位分布圖

3.4 利用EPS地理信息工作站計(jì)算轉(zhuǎn)換參數(shù)

EPS地理信息工作站集信息化測(cè)繪生產(chǎn)技術(shù)體系、工藝流程、生產(chǎn)工具、數(shù)據(jù)管理于一體，是GIS技術(shù)與CAD技術(shù)高度融合的產(chǎn)品，能較好的實(shí)現(xiàn)數(shù)據(jù)轉(zhuǎn)換、圖屬關(guān)聯(lián)、數(shù)據(jù)處理、GIS建庫(kù)、動(dòng)態(tài)更新與成果輸出一體化。

利用EPS地理信息工作站下的坐標(biāo)轉(zhuǎn)換軟件計(jì)算坐標(biāo)轉(zhuǎn)換所需要的四參數(shù)如圖2所示:

圖2 坐標(biāo)轉(zhuǎn)換參數(shù)計(jì)算界面

坐標(biāo)轉(zhuǎn)換的精度是通過(guò)求取轉(zhuǎn)換參數(shù)的重合點(diǎn)的殘差中誤差體現(xiàn)的。

轉(zhuǎn)換時(shí)利用的重合點(diǎn)數(shù)及坐標(biāo)轉(zhuǎn)換殘差中誤差如表1所示。

坐標(biāo)轉(zhuǎn)換殘差中誤差統(tǒng)計(jì)表 表1

坐標(biāo)轉(zhuǎn)換精度估計(jì)依據(jù)下式計(jì)算。

①V(殘差)=重合點(diǎn)轉(zhuǎn)換坐標(biāo)－重合點(diǎn)已知坐標(biāo)

3.5 利用EPS地理信息工作站進(jìn)行數(shù)據(jù)成果轉(zhuǎn)換

(1)編寫轉(zhuǎn)換腳本

EPS腳本語(yǔ)言遵循VBScript和JScript語(yǔ)言標(biāo)準(zhǔn)，支持用戶在EPS平臺(tái)上進(jìn)行二次開(kāi)發(fā)，主要實(shí)現(xiàn)EPS平臺(tái)工程管理、選擇集操作、數(shù)據(jù)編輯、數(shù)據(jù)整理、數(shù)據(jù)轉(zhuǎn)換和系統(tǒng)設(shè)置等處理。

EPS腳本以*.VBS或*.JS為擴(kuò)展名，采用文本文件格式存儲(chǔ)，無(wú)需編譯即可在EPS平臺(tái)中以菜單、工具按鈕或命令行方式執(zhí)行。

腳本分為有模腳本和無(wú)模腳本:

有模腳本

以O(shè)nClick為入口函數(shù)，End Sub結(jié)束，在一個(gè)有模腳本程序中必須定義OnClick函數(shù)。有模腳本不支持與視圖窗口進(jìn)行交互操作。

腳本結(jié)構(gòu)

無(wú)模腳本

以O(shè)nInitScript為入口函數(shù)，在一個(gè)無(wú)模腳本程序中必須定義 OnInitScript、OnExitScript、OnOK和 On-Cancel函數(shù)。有模腳本允許與視圖窗口進(jìn)行交互操作，可以響應(yīng)視圖窗口的鼠標(biāo)消息、按鍵消息和繪圖消息。同時(shí)無(wú)模腳本支持通用參數(shù)錄入對(duì)話框，通過(guò)該對(duì)話框可實(shí)現(xiàn)用戶參數(shù)錄入。

腳本結(jié)構(gòu)

在本次轉(zhuǎn)換過(guò)程中由于轉(zhuǎn)換參數(shù)已定所以采用有模腳本。

坐標(biāo)轉(zhuǎn)換腳本實(shí)現(xiàn):

(2)填寫參數(shù)進(jìn)行轉(zhuǎn)換

在腳本編輯器中調(diào)用腳本，選中需要轉(zhuǎn)換的數(shù)據(jù)點(diǎn)擊執(zhí)行即可完成坐標(biāo)的轉(zhuǎn)換，如圖3所示。

圖3 腳本編輯器界面

4 成果檢核

各項(xiàng)資料坐標(biāo)轉(zhuǎn)換完成后，采用內(nèi)外業(yè)相結(jié)合的方法進(jìn)行檢查驗(yàn)收。

選擇均勻分布蕪湖市區(qū)的11個(gè)沒(méi)參與參數(shù)計(jì)算的重合點(diǎn)進(jìn)行了外部檢核。

作為檢核點(diǎn)的外符合精度，外部檢核點(diǎn)誤差σ:

其中，m為檢核點(diǎn)個(gè)數(shù)，△為檢核點(diǎn)轉(zhuǎn)換坐標(biāo)與其已知值之差。

坐標(biāo)轉(zhuǎn)換精度要通過(guò)外部檢核點(diǎn)誤差與轉(zhuǎn)換(殘差)中誤差全面來(lái)衡量，當(dāng)轉(zhuǎn)換的檢核點(diǎn)誤差與轉(zhuǎn)換(殘差)中誤差數(shù)值接近時(shí)，可認(rèn)為是坐標(biāo)轉(zhuǎn)換精度。

進(jìn)行蕪湖獨(dú)立坐標(biāo)系到CGCS2000轉(zhuǎn)換時(shí)，坐標(biāo)轉(zhuǎn)換精度如表2所示:

外部檢核計(jì)算成果表 表2

外部檢核點(diǎn)誤差:

另外選擇了均勻分布蕪湖市區(qū)的15個(gè)點(diǎn)進(jìn)行蕪湖獨(dú)立坐標(biāo)系、CGCS2000的實(shí)地測(cè)量.

實(shí)際測(cè)量結(jié)果與轉(zhuǎn)換結(jié)果精度統(tǒng)計(jì)如表3所示:

外業(yè)實(shí)測(cè)點(diǎn)計(jì)算成果表 表3

外部實(shí)測(cè)點(diǎn)誤差:

5 小結(jié)

(1)EPS腳本語(yǔ)言遵循VBScript和JScript語(yǔ)言標(biāo)準(zhǔn)，支持用戶二次開(kāi)發(fā)，主要實(shí)現(xiàn)EPS平臺(tái)數(shù)據(jù)整理，數(shù)據(jù)轉(zhuǎn)換等處理。SSProcess是系統(tǒng)缺省定義的對(duì)象，是與EPS平臺(tái)交互的主要入口對(duì)象，SSProcess對(duì)象的功能包括工程管理、選擇集操作、圖形編輯操作、數(shù)據(jù)整理、數(shù)據(jù)轉(zhuǎn)換、坐標(biāo)體系轉(zhuǎn)換和系統(tǒng)設(shè)置等功能。

(2)利用腳本不僅可以實(shí)現(xiàn)不同坐標(biāo)系間二維四參數(shù)轉(zhuǎn)換模型的轉(zhuǎn)換，也同樣適用于二維七參數(shù)轉(zhuǎn)換模型、三維七參數(shù)坐標(biāo)轉(zhuǎn)換模型。

(3)只需計(jì)算出所選模型的參數(shù)寫入腳本執(zhí)行命令即可。利用坐標(biāo)轉(zhuǎn)換腳本轉(zhuǎn)換后的數(shù)據(jù)無(wú)丟失，轉(zhuǎn)換精度高。

［1］李征航，黃勁松.GPS測(cè)量與數(shù)據(jù)處理［M］.武漢大學(xué)出版社，2005.

［2］肖建華，王厚之，李江衛(wèi)等.武漢市現(xiàn)代測(cè)繪基準(zhǔn)體系建設(shè)［J］.城市勘測(cè)，2011(3).

［3］馬祥，李華平.蕪湖市現(xiàn)代測(cè)繪基準(zhǔn)體系的建立［J］.城市勘測(cè)，2013(4).

［4］國(guó)家測(cè)繪地理信息局.啟用2000國(guó)家大地坐標(biāo)系實(shí)施方案［R］.2010－09.

［5］國(guó)家測(cè)繪地理信息局.現(xiàn)有測(cè)繪成果轉(zhuǎn)換到2000國(guó)家大地坐標(biāo)系技術(shù)指南［R］.2010－09.

［6］國(guó)家測(cè)繪地理信息局.2000國(guó)家大地坐標(biāo)系推廣使用技術(shù)指南［R］.2013－6.