張海波，焦殿陽

（1.遼寧省第三地質(zhì)大隊，遼寧朝陽122000；2.遼寧省第九地質(zhì)大隊，遼寧鐵嶺112000）

地質(zhì)調(diào)查中手持GPS轉(zhuǎn)換參數(shù)的確定

張海波1，焦殿陽2

（1.遼寧省第三地質(zhì)大隊，遼寧朝陽122000；2.遼寧省第九地質(zhì)大隊，遼寧鐵嶺112000）

手持GPS在地質(zhì)礦產(chǎn)調(diào)查中常用來測量地質(zhì)觀察點和一般地質(zhì)工程點點位坐標(biāo).本文根據(jù)多年的實踐經(jīng)驗探討手持GPS的使用及坐標(biāo)轉(zhuǎn)換，獲取需要的各種測量成果，具有較強(qiáng)的實用性.

手持GPS；坐標(biāo)系統(tǒng)；轉(zhuǎn)換參數(shù)；地質(zhì)調(diào)查

GPS（Global Positioning System）是由美國在20世紀(jì)70年代開始研究，歷時20年，于1994年全面建成的，具有在海、陸、空進(jìn)行全方位實時三維導(dǎo)航與測量定位能力的新一代衛(wèi)星導(dǎo)航與定位系統(tǒng).隨著技術(shù)的不斷改進(jìn)，硬、軟件的不斷完善，GPS技術(shù)被越來越廣泛地應(yīng)用于各個領(lǐng)域，并開始逐步深入人們的日常生活［1］.

在GPS定位技術(shù)的應(yīng)用和發(fā)展過程中，根據(jù)不同的市場需求，生產(chǎn)廠家研制生產(chǎn)出了不同用途和型號的接收機(jī).無論是高精度的測量型GPS，還是導(dǎo)航型手持GPS，其所提供的坐標(biāo)都是以美國WGS84而建立的，而我們使用的地圖卻屬于54北京坐標(biāo)系或80西安坐標(biāo)系（簡稱地方坐標(biāo)系）.不同的坐標(biāo)系統(tǒng)之間按一定的數(shù)學(xué)模型可計算出其間的轉(zhuǎn)換關(guān)系，也就是我們常說的轉(zhuǎn)換參數(shù).這給我們使用GPS提供了便利.但是求解轉(zhuǎn)換參數(shù)必定是一項專業(yè)性很強(qiáng)的工作，需要收集或采集一些可靠的相關(guān)數(shù)據(jù)并借助計算軟件.盡管如此，導(dǎo)航型手持GPS仍以其獨到的體積小、能耗低、操作簡單、快捷存點、實時快速定位、Teack-Back功能、完善的坐標(biāo)系統(tǒng)、高精度單機(jī)定位諸多優(yōu)越特點而倍受地質(zhì)礦產(chǎn)調(diào)查者的青睞.下面結(jié)合手持GPS在地質(zhì)礦產(chǎn)調(diào)查工作中的實際應(yīng)用，探討如何解決坐標(biāo)轉(zhuǎn)換的問題.

1 確定自定義坐標(biāo)格式

確定自定義坐標(biāo)格式（User UTM Grid）是高斯-克呂格投影的需要，屬于地圖投影方面的研究范疇.測量工作的一個重要內(nèi)容就是將地球表面上的地物地貌按一定的要求用平面圖形表示出來，這就是測繪各種比例尺地形圖.這就需要將橢球（近似地說是圓球）上的圖形投影到平面上，在平面上表示坐標(biāo)的方法是平面直角坐標(biāo)系.尤其重要的是，如果不用平面直角坐標(biāo)系，則在橢球（或圓球）上的各點位置要用經(jīng)緯度來表示.這樣，計算經(jīng)緯度或已知經(jīng)緯度來反算距離及方位角，都將是一件很復(fù)雜的工作.如果有一個平面直角坐標(biāo)系，在平面上（而不是在球面上）處理測量數(shù)據(jù)，將使計算工作大大簡化.高斯-克呂格投影實現(xiàn)了將球面投影轉(zhuǎn)為平面的過程.但是高斯-克呂格投影只是保持了投影后角度不變，長度變形隨著遠(yuǎn)離中央子午線會逐漸變大，為此常采用分帶投影來限制長度變形誤差的影響［2］.確定自定義坐標(biāo)格式就是確定投影帶和坐標(biāo)原點的平移.值得一提的是，除通用的投影帶劃分以外，還有一個根據(jù)地域所處地理位置不同而并不使用坐標(biāo)的自然值的問題.例如，我國位于赤道以北，縱軸（或稱X軸）均為正值，而每一投影帶由赤道投影后形成的橫軸（或稱Y軸）由縱軸分成東西兩部分，東者為正，西者為負(fù)，為了使橫坐標(biāo)均為正值起見，在3°及6°帶均加入了常數(shù)500km，并冠以投影帶號來表示.

當(dāng)我們使用一部新的GPS或到一個新的勘查區(qū)工作時，首先要做的就是確定自定義坐標(biāo)格式.確定自定義坐標(biāo)格式的關(guān)鍵是如何確定勘查區(qū)所在投影帶中央子午線的經(jīng)度.這是因為在使用統(tǒng)一的國家坐標(biāo)系或地方坐標(biāo)系時，因勘查區(qū)所處地理位置不同，中央子午線的經(jīng)度是需要相應(yīng)變更的.那么如何方便快捷地確定中央子午線的經(jīng)度呢？我們在進(jìn)入新的勘查區(qū)之前，要收集包括地形圖在內(nèi)的有關(guān)資料.國家基本比例尺地形圖中的中小比例尺地形圖橫軸坐標(biāo)前都冠以帶號，內(nèi)圖廓還標(biāo)有經(jīng)緯度，這就為我們確定勘查區(qū)中央子午線經(jīng)度提供了最基本的條件.

1.1 根據(jù)投影帶號確定勘查區(qū)中央子午線經(jīng)度

我國區(qū)域內(nèi)帶號均為兩位正整數(shù)，只要把確定的帶號代入下式（1）就可計算出中央子午線的經(jīng)度.根據(jù)投影帶號計算中央子午線經(jīng)度的公式：

式中L0代表中央子午線經(jīng)度，N3、N6分別代表3°帶帶號、6°帶帶號［3］.

1.2 根據(jù)大地坐標(biāo)經(jīng)度值確定勘查區(qū)中央子午線經(jīng)度

用經(jīng)度值（L）確定勘查區(qū)中央子午線經(jīng)度需按以下步驟進(jìn)行：首先要確定工作底圖是3°分帶還是6°分帶，然后利用勘查區(qū)內(nèi)任一點經(jīng)度值計算其所在投影帶帶號，再將帶號代入上式（1）即可算出中央子午線經(jīng)度.利用經(jīng)度值計算6°帶帶號的算法為，用L值整數(shù)位除以 6，取整數(shù)商加 1.例如，某地 L 為 120°18′28″，根據(jù)計算（120°÷6+1=21）得知其所在 6°帶帶號為 21.利用經(jīng)度值計算3°帶帶號的算法為，將L值換成以度為單位后除以3，按四舍五入取整數(shù)商即為帶號.例如：某地 L 為 120°18′28″，根據(jù)計算（120.31÷3=40.103，四舍五入取整數(shù)商為40），得知其所在3°分帶帶號為 40［4］.

3°和 6°分帶其東西跨度分別為 3°和 6°，被中央子午線分成東西對稱兩部分.

我國位于東半球，中央子午線設(shè)置均為東經(jīng)（E）并輸入相應(yīng)的中央子午線經(jīng)度值，投影比例設(shè)置為1，東西偏差設(shè)置為500 000，南北偏差設(shè)置為0，并設(shè)置單位為m.

2 確定自定義坐標(biāo)系統(tǒng)的轉(zhuǎn)換參數(shù)

GPS默認(rèn)坐標(biāo)系統(tǒng)為美國WGS84坐標(biāo)系統(tǒng).若想得到我們所用的地方坐標(biāo)，必須要有相應(yīng)的轉(zhuǎn)換參數(shù).這里所說的地方坐標(biāo)就是指我們常用的54北京坐標(biāo)系和80西安坐標(biāo)系以及改變了中央子午線或投影面高度的城市坐標(biāo)系的坐標(biāo)的統(tǒng)稱.這就要求在確定自定義坐標(biāo)系統(tǒng)（User）轉(zhuǎn)換參數(shù)之前，首先要確定需要什么坐標(biāo)系的成果.只有這樣才能按不同的橢球參數(shù)（確定橢球大小和形狀的參數(shù)）來計算與之相適應(yīng)的WGS84坐標(biāo)系與地方坐標(biāo)系的轉(zhuǎn)換參數(shù)［5］.

WGS84坐標(biāo)系統(tǒng)為GPS默認(rèn)系統(tǒng)，與我們常用的地方坐標(biāo)系統(tǒng)是兩個不同橢球參數(shù)下定義的坐標(biāo)系，它們之間的轉(zhuǎn)換關(guān)系是通過特定的數(shù)學(xué)模型計算出來的，經(jīng)常用到的有布爾莎、莫洛登斯基、海爾曼特等幾種數(shù)學(xué)模型［6］.嚴(yán)密的轉(zhuǎn)換關(guān)系為七參數(shù)（Dx、Dy、Dz、Wx、Wy、Wz、K，它們分別是 X、Y、Z 三軸的平移與旋轉(zhuǎn)及比例系數(shù)）.由于手持GPS單機(jī)定位精度要求并不高，一般在±15m，因此在主版設(shè)計時選用了近似七參數(shù)，也就是我們常用的三參數(shù) Dx、Dy、Dz，忽略了旋轉(zhuǎn)參數(shù)Wx、Wy、Wz，實際上旋轉(zhuǎn)參數(shù)也非常?。槐壤齾?shù)近似為1.所以一般手持GPS自定義坐標(biāo)系統(tǒng)轉(zhuǎn)換參數(shù)設(shè)置的變量參數(shù)只有 3 個：Dx、Dy、Dz，至于 DA（橢球長軸之差）、DF（橢球扁率之差）對于不同橢球間是固定常數(shù).三參數(shù)Dx、Dy、Dz是隨著勘查區(qū)地理位置不同而變化的.

3 三參數(shù)的確定方法

3.1 利用已知坐標(biāo)計算

收集或采集勘查區(qū)及周邊1～3個或多個三等以上已知當(dāng)?shù)刈鴺?biāo)和WGS84坐標(biāo)的控制點成果.需要指出的是，1個點只能計算1組三參數(shù)，2個點計算2組有比較，3個點計算3組有檢核，3個（包括3個）以上可以比較取中數(shù)使用.一般情況下在省測繪局測繪科技資料處可收集到上述成果.現(xiàn)在GPS靜態(tài)測量也很普遍，這為我們獲取WGS84坐標(biāo)提供了便利.

三參數(shù)的計算程序可使用廣州中海達(dá)測繪儀器有限責(zé)任公司的靜態(tài)測量后處理軟件，采用莫洛登斯基數(shù)學(xué)模型來計算.具體做法是，進(jìn)入程序，打開工具到計算七參數(shù)界面，將WGS84坐標(biāo)作為源坐標(biāo)輸入，該點的當(dāng)?shù)刈鴺?biāo)作為目標(biāo)坐標(biāo)輸入，點擊計算.因為是單點，軟件會自動選用莫洛登斯基數(shù)學(xué)模型來計算并有對話框提示，記錄下計算值，如此輸入幾對點進(jìn)行計算，通過比較后取中數(shù)，并將計算出的Dx、Dy、Dz反號置入手持GPS就可以了.

下面就某勘查區(qū)的實例來說明計算三參數(shù)的具體做法：收集勘查區(qū)A、B、C三控制點當(dāng)?shù)刈鴺?biāo)成果，再做包括三點在內(nèi)的GPS靜態(tài)控制測量，取得WGS84坐標(biāo)，使用中海達(dá)靜態(tài)后處理軟件進(jìn)行計算（表 1、2）.

最終結(jié)果Dx=-103.8m，Dy=-52.2m，Dz=1.6m.在該勘查區(qū)取DA=-3，DF=0.設(shè)置此參數(shù)到手持GPS，并對已知點進(jìn)行測量比較，平面精度均在±3m，較接收機(jī)標(biāo)稱精度±15m有很大提高.高程測量值也比較接近已知成果.

事實上獲取WGS84橢球參數(shù)下約束平差的WGS84坐標(biāo)成果不是一件容易的事情，所以常采用自由網(wǎng)平差結(jié)果.那么采用手持GPS測量的WGS84坐標(biāo)來計算三參數(shù)可不可以呢？答案是肯定的.然而因手持GPS單點定位精度的限制，以其計算出的轉(zhuǎn)換參數(shù)來定位點，精度不高.

3.2 收集使用經(jīng)驗成果

在手持GPS的實際應(yīng)用中，會有大多數(shù)使用者不掌握三參數(shù)的計算程序和計算方法，有的即使掌握三參數(shù)的計算程序和計算方法，然而已知點的地方坐標(biāo)系和WGS84坐標(biāo)系成果（尤其是WGS84坐標(biāo)系成果）的收集或采集都是很困難的，那么只能靠收集使用經(jīng)驗成果，把收集來的三參數(shù)輸入GPS后到已知點進(jìn)行測量比較，只要滿足設(shè)計精度要求或達(dá)到接收機(jī)標(biāo)稱精度就可以采用.

表2 計算結(jié)果表Table 2 Calculated result

3.3 其他方法

從數(shù)學(xué)上來講，轉(zhuǎn)換參數(shù)是點的位置與橢球參數(shù)的矩陣函數(shù)，其數(shù)理關(guān)系復(fù)雜、計算過程較繁瑣，況且有些參數(shù)只能在專用數(shù)表中才能查得到.如果不是從事專門理論研究，而只限于手持GPS的使用方面，上述方法是方便實用的轉(zhuǎn)換參數(shù)計算方法.

上述獲取轉(zhuǎn)換參數(shù)的方法都是在比較理想狀態(tài)下，地質(zhì)礦產(chǎn)調(diào)查往往會遇到意想不到的環(huán)境，如廣袤森林、戈壁灘、沼澤等人跡罕至的地方.這些地方或因為經(jīng)濟(jì)欠發(fā)達(dá)，國家基本控制網(wǎng)尚未覆蓋或是密度不夠，在很難收集到已知地方坐標(biāo)的情況下，只能采用相鄰?fù)队皫蛲鈹U(kuò)展，利用相鄰地區(qū)轉(zhuǎn)換參數(shù)進(jìn)行測量.測量過程中采集地形圖上能圖解坐標(biāo)的明顯地物點坐標(biāo)數(shù)據(jù)，再根據(jù)比較結(jié)果作固定差改正，最后根據(jù)需要對測量成果進(jìn)行換帶計算.分帶進(jìn)行投影的目的就是為了限制投影變形影響制圖的精度，特殊情況下采用相鄰?fù)队皫蛲鈹U(kuò)展勢必產(chǎn)生變形，但變形值的大小與手持GPS定位精度相比，完全可以忽略不計.

4 結(jié)論

利用已知坐標(biāo)計算三參數(shù)時，已知坐標(biāo)點幾何分布最好呈等邊三角形或正多邊形（中點多邊形最佳）.依上述實例做法，在不同地區(qū)地質(zhì)勘查項目中使用手持GPS，實踐表明，在已知坐標(biāo)點連成的多邊形內(nèi)及不超過已知坐標(biāo)點間距一半的外部范圍內(nèi)，使用這樣一組已知坐標(biāo)點計算的三參數(shù)，在定位點的坐標(biāo)時，平面精度明顯提高而且穩(wěn)定，只是在提高高程精度方面不太顯著.

表1 已知點成果表Table 1 Data of known points

三參數(shù)Dx、Dy、Dz是隨著勘查區(qū)地理位置不同而變化的，所以想謀求一組普遍適用所有地區(qū)的轉(zhuǎn)換參數(shù)是不現(xiàn)實的.

地質(zhì)礦產(chǎn)調(diào)查工作有其點多、面廣、流動性強(qiáng)的特點，一般勘查區(qū)工作周期只有幾個月，有遠(yuǎn)景的或一到幾年，且限于局域范圍內(nèi).除探槽、淺井、鉆孔等山地工程需要精確定位外，地質(zhì)觀察點、一般地質(zhì)工程點定位工作量相當(dāng)大，且精度要求并不高，特別是地質(zhì)界線點一般只要求平面位置準(zhǔn)確.盡管如此，合理的轉(zhuǎn)換參數(shù)才能提高和保證定位點的精度.比較原來靠地質(zhì)工作人員判讀并在地形圖上標(biāo)繪，手持GPS定位技術(shù)突顯優(yōu)越.

總之，使用靜態(tài)測量后處理軟件來計算局域范圍內(nèi)轉(zhuǎn)換參數(shù)是一種方便、快捷、實用、有效的途徑之一.

［1］GARMIN（Asia）Corp.GARMINGPS72說明書（中文版）［Z］.

［2］武漢測繪學(xué)院《測量學(xué)》編寫組.測量學(xué)［M］.北京：測繪出版社，1985.

［3］王德修.礦區(qū)地形測量［M］.北京：煤炭工業(yè)出版社，1987.

［4］孫洪梅.手持GPS在錦州排污口監(jiān)測中應(yīng)用［J］.遼寧測繪，2010（2）.

［5］GB50026-2007，工程測量規(guī)范［S］.

［6］中海達(dá)測繪儀器有限公司.手簿軟件說明書［Z］.

THE DETERMINATION OF CONVERSION PARAMETERS OF HANDHELD GPS IN GEOLOGICAL SURVEY

ZHANG Hai-bo1,JIAO Dian-yang2

（1.No.3 Geological Brigade of Liaoning Province,Chaoyang 122000,Liaoning Province,China;2.No.9 Geological brigade of Liaoning Province,Tieling 112000,Liaoning Province,China）

In geology and mineral resources survey,the handheld GPS can be applied to measure the coordinates of geological observation point and geological engineering point.Based on practical experiences,this paper discusses the application and coordinate conversion of handheld GPS to obtain surveying result.It is practically significant in geology and mineral resources survey.

handheld GPS;coordinate system;conversion parameter;geological survey

1671-1947（2011）03-0222-04

P628

A

2011－04－13；

2011-05-05.編輯：歐陽兆灼，張哲.

張海波（1962—），男，測繪工程師，從事地質(zhì)礦產(chǎn)勘查工作中的控制測量、地形測量、地質(zhì)工程測量等工作，通信地址遼寧省朝陽市長江路四段29號，E-mail//zhanghaibocy@163.com