李志暉

（1.河南省有色金屬地質(zhì)勘查總院，河南 鄭州 450052；2.河南省有色金屬深部找礦勘查技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，河南 鄭州 450052）

獨(dú)立坐標(biāo)系在建立境外礦區(qū)GPS控制網(wǎng)中的應(yīng)用

李志暉1,2

（1.河南省有色金屬地質(zhì)勘查總院，河南 鄭州 450052；2.河南省有色金屬深部找礦勘查技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，河南 鄭州 450052）

總結(jié)了獨(dú)立坐標(biāo)系在老撾巴松縣鋁土礦勘探測(cè)量中的應(yīng)用，以期提供一套行之有效的解決類(lèi)似問(wèn)題的方法。主要論述了坐標(biāo)系的選擇，獨(dú)立坐標(biāo)系的建立及實(shí)現(xiàn)，數(shù)據(jù)處理及平差計(jì)算。在對(duì)GPS控制網(wǎng)進(jìn)行約束平差及精度評(píng)定后，得出“該GPS控制網(wǎng)優(yōu)于E級(jí)GPS控制網(wǎng)的精度，完全可以作為該礦區(qū)的平面首級(jí)控制”的結(jié)論，印證了獨(dú)立坐標(biāo)系建立境外礦區(qū)GPS控制網(wǎng)的可靠性。

獨(dú)立坐標(biāo)系；GPS控制網(wǎng)；測(cè)量方法；鋁土礦區(qū)；老撾巴松縣

1 測(cè)區(qū)概況

老撾巴松鋁土礦區(qū)位于老撾南部占巴色省巴松縣境內(nèi)南東部，礦區(qū)周邊地勢(shì)較高，最高海拔915．0 m，最低海拔416．0 m，相對(duì)高差約為500 m，具有“波羅芬高原”的地貌特征；屬熱帶季風(fēng)型—熱帶雨林氣候，一年兩季（雨季、旱季），年平均降雨量為1 400～2 000 mm，進(jìn)駐工區(qū)時(shí)正值雨季，給測(cè)量工作帶來(lái)不利影響。礦區(qū)距巴色市東約90 km，距巴松縣城東約40 km。13號(hào)公路沿湄公河?xùn)|岸南下可直達(dá)占巴色省巴色市，從巴色市有省際公路經(jīng)過(guò)巴松縣，巴松縣到礦區(qū)有公路經(jīng)過(guò)，礦區(qū)北部和南部有大車(chē)路經(jīng)過(guò)，交通較方便。測(cè)區(qū)內(nèi)居民地不多，呈零星分布；農(nóng)業(yè)生產(chǎn)以種植咖啡為主，區(qū)內(nèi)咖啡園較多；植被均為原始森林，地表被森林完全覆蓋，除零星居民地和幾條土路外基本見(jiàn)不到地面，開(kāi)展測(cè)量工作難度較大。測(cè)區(qū)內(nèi)沒(méi)有任何測(cè)繪資料，要進(jìn)行測(cè)繪工作，首先需要建立控制網(wǎng)。由于GPS網(wǎng)具有高精度、高效、無(wú)需通視等優(yōu)點(diǎn)，因此在該測(cè)區(qū)布設(shè)獨(dú)立GPS控制網(wǎng)。

2 基礎(chǔ)資料與坐標(biāo)系選擇

該礦區(qū)及其周邊地區(qū)未經(jīng)工業(yè)開(kāi)采和資源普查，測(cè)量基礎(chǔ)資料缺乏，導(dǎo)致本次礦區(qū)測(cè)量工作開(kāi)展前后一直未收集到礦區(qū)平面和高程控制點(diǎn)。由于收集的礦區(qū)主要地質(zhì)資料是以1∶25 000地質(zhì)地形圖為基礎(chǔ)的，與1∶50 000地形圖坐標(biāo)系一致（采用UTM投影）[1]，因此，采用與1∶50 000地形圖接近的獨(dú)立坐標(biāo)系——1∶50 000地形圖的UTM投影方式和印度1960坐標(biāo)系（Everest1830橢球參數(shù)）。

3 獨(dú)立坐標(biāo)系的建立及其實(shí)現(xiàn)

3.1 建立地方獨(dú)立坐標(biāo)系的主要元素

1）中央子午線(xiàn)。中央子午線(xiàn)可以和國(guó)家坐標(biāo)系標(biāo)準(zhǔn)帶中央子午線(xiàn)重合，但當(dāng)測(cè)區(qū)離標(biāo)準(zhǔn)帶中央子午線(xiàn)較遠(yuǎn)時(shí)，可取過(guò)測(cè)區(qū)中心點(diǎn)或過(guò)某點(diǎn)的經(jīng)線(xiàn)作為中央子午線(xiàn)。如果僅移動(dòng)中央子午線(xiàn)能夠解決投影變形，那么將起算點(diǎn)坐標(biāo)進(jìn)行換帶就可建立地方獨(dú)立坐標(biāo)系。

2）投影面。若移動(dòng)中央子午線(xiàn)不能解決投影變形，就要考慮選擇適當(dāng)?shù)耐队懊??？蛇x擇測(cè)區(qū)的平均高程面，也可選擇抵償高程面作為投影面[2]。

3）地方坐標(biāo)系參考橢球。通過(guò)對(duì)地方坐標(biāo)系參考橢球幾何元素、定位及定向的確定，使得橢球面與投影面擬合最好，這樣投影變形可減到最小，同時(shí)要求便于與國(guó)家坐標(biāo)系統(tǒng)進(jìn)行換算[3-4]。

3.2 獨(dú)立坐標(biāo)系的建立方案和測(cè)量方法

首先確定獨(dú)立坐標(biāo)系的基點(diǎn)坐標(biāo)，再施測(cè)礦區(qū)GPS控制網(wǎng)。由于測(cè)區(qū)范圍內(nèi)森林覆蓋，且全部是山地，通視條件差，根據(jù)GPS布網(wǎng)原則[5]，在礦區(qū)布設(shè)1個(gè)坐標(biāo)原點(diǎn)（G011），這個(gè)坐標(biāo)原點(diǎn)和與它通視的方向點(diǎn)（G017）及其他34個(gè)點(diǎn)組成GPS控制網(wǎng) （見(jiàn)圖1），導(dǎo)線(xiàn)布置形式見(jiàn)圖2。該獨(dú)立坐標(biāo)系建立方案和測(cè)量方法的具體作法如下：

1）閉合導(dǎo)線(xiàn)觀測(cè)。根據(jù)踏勘情況，選擇工區(qū)北部較平坦的公路進(jìn)行獨(dú)立閉合環(huán)觀測(cè)，全站儀用三連架的方式從G011開(kāi)始經(jīng)過(guò)G017布設(shè)成直伸導(dǎo)線(xiàn)到達(dá)G018，水平距離全長(zhǎng)1 903．244 m；又從G018折回閉合在G011，形成閉合導(dǎo)線(xiàn)，共20條邊、18個(gè)折角。經(jīng)過(guò)各項(xiàng)改正計(jì)算，長(zhǎng)度閉合差為0．01 m，角度閉合差為7 s，高程閉合差為0．002 m，達(dá)到規(guī)范要求，并進(jìn)行了平差計(jì)算。具體數(shù)據(jù)見(jiàn)獨(dú)立坐標(biāo)系閉合導(dǎo)線(xiàn)計(jì)算表（表1、2）。

圖1 老撾巴松縣鋁土礦測(cè)區(qū)GPS控制網(wǎng)

圖2 老撾巴松縣鋁土礦測(cè)區(qū)導(dǎo)線(xiàn)布置示意圖

表1 平距平差計(jì)算表

2）平面坐標(biāo)推算。起算點(diǎn)確定：為了精確定位和定向，在坐標(biāo)原點(diǎn)（G011）和方位點(diǎn)（G018）用兩臺(tái)南方9600 GPS接收機(jī)進(jìn)行12 h的數(shù)據(jù)采集。為了使建立的獨(dú)立坐標(biāo)系盡可能接近1∶50 000地形圖的坐標(biāo)系，首先利用靜態(tài)GPS觀測(cè)結(jié)果（B：15°04'2896224"N，L： 106°32'3345180"E），采用Everest1830橢球參數(shù)，進(jìn)行UTM投影得到G011的平面坐標(biāo)（X=1 666 846.003，Y=666 223.587）。方向確定：用同樣的投影方式得到G018（B：15°04'3800656"N，L：106°33'2986746"E）的平面坐標(biāo)（X=1 667 135.816，Y=667 906.321），計(jì)算出G011到G018的方位，作為推算閉合導(dǎo)線(xiàn)上點(diǎn)平面坐標(biāo)的方位閉合條件。由此，利用G011的假定坐標(biāo)，最終推算閉合導(dǎo)線(xiàn)上各導(dǎo)線(xiàn)點(diǎn)的平面坐標(biāo)。

3）高程推算。在6 339Ⅲ幅1∶50 000圖上查得G011位置的高程為791 m，利用閉合導(dǎo)線(xiàn)的高差結(jié)果求出G017、G018的高程。

經(jīng)過(guò)以上步驟，已在本礦區(qū)建立了比較可靠的獨(dú)立坐標(biāo)系。據(jù)此將推算的G011、G017、G018坐標(biāo)作為本礦區(qū)GPS控制網(wǎng)的起算點(diǎn)。

4 數(shù)據(jù)處理和平差計(jì)算

4.1 基線(xiàn)處理及檢驗(yàn)

用南方測(cè)繪 GPS 4．4000 數(shù)據(jù)處理軟件對(duì)外業(yè)觀測(cè)數(shù)據(jù)進(jìn)行基線(xiàn)解算，全部結(jié)果均為雙差固定解。

依據(jù)GB/T 18314-2009《規(guī)范》12．2．7的要求，E級(jí)GPS網(wǎng)獨(dú)立閉合環(huán)或附合路線(xiàn)閉合差Ws和各坐標(biāo)分量閉合差（Wx，Wy，Wz）應(yīng)滿(mǎn)足以下要求；復(fù)測(cè)基線(xiàn)的長(zhǎng)度較差應(yīng)符合。經(jīng)檢驗(yàn)，同步環(huán)、異步環(huán)、復(fù)測(cè)基線(xiàn)較差均滿(mǎn)足要求。

4.2 GPS網(wǎng)的平差計(jì)算及精度評(píng)定

表2 高程平差計(jì)算表

網(wǎng)平差處理使用南方測(cè)繪 GPS 4．4000 數(shù)據(jù)處理軟件，選定合格的基線(xiàn)解，構(gòu)成GPS控制網(wǎng)。先對(duì)基線(xiàn)進(jìn)行WGS84坐標(biāo)系下的自由網(wǎng)平差，基線(xiàn)向量的改正數(shù)絕對(duì)值應(yīng)滿(mǎn)足以下要求[6]：VΔX≤3σ，VΔY≤3σ，VΔZ≤3σ。三維無(wú)約束平差后，基線(xiàn)向量的改正數(shù)絕對(duì)值見(jiàn)表3，均滿(mǎn)足限差要求。

表3 基線(xiàn)向量的改正數(shù)絕對(duì)值/cm

在地方獨(dú)立坐標(biāo)系下，采用印度1960坐標(biāo)系Everest1830橢球參數(shù)：

長(zhǎng)半軸為6 377 276.345 m；短半軸為6 356 075.000 m；扁率為1/300.801 7；投影為橫軸墨卡托投影；中央子午線(xiàn)為105°；坐標(biāo)原點(diǎn)（G011）的數(shù)值為X=1 666 846.003，Y=666 223.587。以坐標(biāo)原點(diǎn)（G011）、方位點(diǎn)（G018）及方向點(diǎn)（G017）的坐標(biāo)作為已知點(diǎn)進(jìn)行二維約束平差，經(jīng)過(guò)二維約束平差后得到各精度指標(biāo)見(jiàn)表4。最佳相對(duì)中誤差為1/1 139 728，最差相對(duì)中誤差為1/76 388，平均相對(duì)中誤差為1/287 854。

表4 點(diǎn)位中誤差/cm

從精度分析可知，最弱點(diǎn)的點(diǎn)位中誤差為0．475 cm，最弱邊的相對(duì)中誤差為1/76 388，而E級(jí)GPS控制網(wǎng)最弱邊的相對(duì)中誤差為1/45 000，點(diǎn)位誤差為5．0 cm，因此該GPS控制網(wǎng)優(yōu)于E級(jí)GPS控制網(wǎng)的精度，完全可以作為該礦區(qū)的平面首級(jí)控制[7]。

5 結(jié) 語(yǔ)

獨(dú)立坐標(biāo)系在建立老撾巴松縣鋁土礦區(qū)GPS控制網(wǎng)中的應(yīng)用，體現(xiàn)了其不同于在國(guó)內(nèi)布網(wǎng)的特點(diǎn)：

①區(qū)內(nèi)無(wú)可供參考的測(cè)量基礎(chǔ)資料，無(wú)法向測(cè)區(qū)內(nèi)引入控制點(diǎn)；②在坐標(biāo)原點(diǎn)和方位點(diǎn)進(jìn)行12 h的觀測(cè)，其精度指標(biāo)完全可以作為礦區(qū)E級(jí)GPS控制網(wǎng)的起算數(shù)據(jù)；③經(jīng)本礦區(qū)實(shí)踐證明，建立的獨(dú)立GPS控制網(wǎng)優(yōu)于E級(jí)GPS控制網(wǎng)的精度，完全可以作為該礦區(qū)的平面首級(jí)控制。

[1] 熊忠招．淺談UTM 投影下獨(dú)立坐標(biāo)系統(tǒng)建立[J]．地理空間信息,2010,8(2)∶41-43

[2] 張勝利．地方獨(dú)立坐標(biāo)系的建立[J]．陜西水利水電技術(shù),2007(1)∶43-46

[3] 張述清,李永云．地方獨(dú)立坐標(biāo)系統(tǒng)的建立及其實(shí)現(xiàn)[J]．測(cè)繪工程,2007,16(4)∶22-29

[4] 陳士銀．建立地方獨(dú)立坐標(biāo)系的方法[J]．測(cè)繪通報(bào),1997(10)∶5-7

[5] 孔祥元,梅是義．控制測(cè)量學(xué)[M]．武漢∶武漢測(cè)繪科技大學(xué)出版社,1996

[6] 劉大杰,施一民,過(guò)靜珺．全球定位系統(tǒng)(GPS)的原理與數(shù)據(jù)處理[M]．上海∶同濟(jì)大學(xué)出版社,1996

[7] 朱照榮．赤道幾內(nèi)亞新城GPS獨(dú)立控制網(wǎng)的建立[J]．北京測(cè)繪,2004(3)∶27-29

圖5 地價(jià)水平值統(tǒng)計(jì)柱狀圖

圖6 地價(jià)水平狀況分析圖

4 結(jié) 語(yǔ)

該基準(zhǔn)地價(jià)查詢(xún)公示系統(tǒng)集成了基準(zhǔn)地價(jià)成果管理、顯示、查詢(xún)、統(tǒng)計(jì)、分析、網(wǎng)上發(fā)布等一體化功能，將其應(yīng)用于地價(jià)評(píng)估和日常管理工作，科學(xué)實(shí)用、方便快捷，為地價(jià)管理的信息化提供了依據(jù)，同時(shí)也為實(shí)現(xiàn)土地市場(chǎng)公開(kāi)、公平、公正的健康發(fā)展奠定了良好的基礎(chǔ)。

參考文獻(xiàn)

[1] 唐凌奇．WebGIS在地價(jià)信息發(fā)布中的應(yīng)用研究[D]．武漢∶武漢大學(xué),2004

[2] 成都市國(guó)土資源局．成都市中心城區(qū)基準(zhǔn)地價(jià)更新研究報(bào)告[R]．成都∶成都市國(guó)土資源局,2007

[3] 毛國(guó)梁．淺談GeoMedia WebMap的配置與開(kāi)發(fā)[J]．科技情報(bào)開(kāi)發(fā)與經(jīng)濟(jì),2005,15(22)∶214-215

[4] 曹志剛,李滿(mǎn)春．設(shè)計(jì)模式在GeoMedia WebMap地圖數(shù)據(jù)發(fā)布中的應(yīng)用[J]．計(jì)算機(jī)應(yīng)用研究,2005(4)∶84-86

[5] 方從剛,陳先偉,黃成毅．基于WebGIS的成都市基準(zhǔn)地價(jià)信息系統(tǒng)的構(gòu)建[J]．計(jì)算機(jī)與數(shù)字工程,2010(3)∶67-69

[6] 楊俊艷．淺談太原市基準(zhǔn)地價(jià)查詢(xún)系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J]．科技情報(bào)開(kāi)發(fā)與經(jīng)濟(jì),2011,21(31)∶129-131

[7] 徐京華,李何超,汪四文．四川省城鎮(zhèn)土地級(jí)別與基準(zhǔn)地價(jià)圖集及查詢(xún)系統(tǒng)研制 [J]．測(cè)繪,2010,33(5)∶195-198

第一作者簡(jiǎn)介：杜艷琴，碩士，工程師，從事土地利用及城市發(fā)展相關(guān)研究。

P258

B

1672-4623（2016）09-0079-03

10.3969/j.issn.1672-4623.2016.09.026

李志暉，碩士，工程師，主要從事地質(zhì)礦產(chǎn)勘查及測(cè)繪工作。

2015-02-11。

項(xiàng)目來(lái)源：中色國(guó)際礦業(yè)國(guó)際合作資助項(xiàng)目（L-2005-01）。