王建營,張志全,汪偉
(天津市測繪院,天津 300381)
隨著國民經(jīng)濟的快速增長,我國對海島資源的逐步開發(fā)、利用日益迫切。精確的島礁測繪數(shù)據(jù),對我國宣示領海主權(quán)、維護海洋權(quán)益、開發(fā)海洋資源具有重要意義。傳統(tǒng)的測量方法,只能測量距大陸架較近的島礁,隨著GPS技術(shù)的廣泛應用結(jié)合航空攝影測量等技術(shù),理論上可測量海洋上所有島嶼。
本文主要討論的就是如何利用GPS方法在海島上建立CGCS2000國家大地坐標系控制網(wǎng)。
經(jīng)國務院批準,自2008年7月1日起在全國啟用2000國家大地坐標系(CGCS2000)。為獲取海島測量控制網(wǎng)CGCS2000坐標,需要與周邊已知控制點進行聯(lián)測。由于我國海島主要分布在東海和南海,如果僅與我國內(nèi)陸的已知點聯(lián)測,控制網(wǎng)網(wǎng)形較差,所以最好的方法是合理選取周邊3個或以上國際IGS站進行聯(lián)測。同時求取IGS站點CGCS2000坐標,以其作為起算數(shù)據(jù),進行聯(lián)合數(shù)據(jù)處理,最終獲得海島GPS控制網(wǎng)國家2000坐標??刂凭W(wǎng)建立的流程如圖1所示。
圖1 海島測量控制網(wǎng)建立流程圖
根據(jù)《無居民海島使用測量規(guī)范》,界址點坐標的點位中誤差不超過±0.5 m??梢圆捎盟牡菺PS控制網(wǎng)作為海島測量的首級控制網(wǎng)。
GPS網(wǎng)布設、觀測均應參照《衛(wèi)星定位城市測量技術(shù)規(guī)范》CJJ/T 73-2010中的要求。
ITRF當前發(fā)布的是ITRF2000、2005和2008三個框架下觀測時段中心歷元坐標,而我國CGCS2000系為ITRF1997框架2000.0歷元,兩者之間存在坐標差。所以要獲得CGCS2000坐標,還需對IGS站數(shù)據(jù)進行板塊運動改正、歷元轉(zhuǎn)換、框架轉(zhuǎn)換等工作。
(1)獲取國際IGS站數(shù)據(jù)
IGS站數(shù)據(jù)可以在很多網(wǎng)站下載,如ftp://garner.ucsd.edu等,還需同時下載導航文件及天線信息等文件。
(2)板塊運動改正
即根據(jù)IGS站點速度場,將觀測時段中心歷元站點坐標轉(zhuǎn)換至所需歷元(2000.0歷元)時段站點坐標。IGS站速度可在ITRFxx_GPS_SSC文件中得到。轉(zhuǎn)換公式如下:
其中:Xt1、Yt1、Zt1為觀測時段站點坐標;
Xt2、Yt2、Zt2為所求時段(2000.0歷元)站點坐標;
VX、VY、VZ為站點速度場。
實際上,不同歷元下站點坐標也可直接在http://sopac.ucsd.edu網(wǎng)站中查詢。由此可以獲得ITRF2000框架2000.0歷元站點坐標。
(3)歷元歸算
由于ITRF發(fā)布的ITRF 2000與以前框架的轉(zhuǎn)換關系參數(shù)為歷元1997.0(如下表),因此需要將表中各參數(shù)歸算至歷元2000.0。
轉(zhuǎn)換公式:
mr=1/206265為秒與弧度轉(zhuǎn)換關系。
從ITRF2000轉(zhuǎn)換到ITRF1997框架的轉(zhuǎn)換參數(shù)(1997.0歷元),如表1所示。
從框架ITRF2000轉(zhuǎn)換到ITRF1997的轉(zhuǎn)換參數(shù)(歷元1997.0) 表1
代入公式,則 ITRF2000轉(zhuǎn)換到 ITRF1997框架(2000歷元)的轉(zhuǎn)換參數(shù)為:Tx=0.67,Ty=0.43,Tz=-2.27,s=1.58,εx=0,εy=0,εz=0.06/206265
(4)框架轉(zhuǎn)換(獲取CGCS2000坐標)
根據(jù)上述轉(zhuǎn)換關系對經(jīng)板塊運動改正得到的ITRF2000框架2000.0歷元站點坐標進行轉(zhuǎn)換,獲得IGS站ITRF1997框架2000歷元(即CGCS2000)坐標。轉(zhuǎn)換公式為:
(5)聯(lián)合基線處理、平差計算
將四等GPS靜態(tài)觀測數(shù)據(jù)與周邊選取IGS站數(shù)據(jù)聯(lián)合基線處理、三維約束平差,獲得海島測量控制網(wǎng)CGCS2000大地坐標。
浙江省象山縣旦門山島測量是我國首例對無居民海島實施測量,其首級控制網(wǎng)即采用上述方法建立而成。旦門山島總面積約為 0.94 km2,全島共布設四等GPS控制點4個,由3臺TrimbleR8雙頻GPS接收機,分兩個時段完成??刂泣c布設情況如圖2所示。
圖2 控制點布設示意圖
為了獲得2000國家大地坐標系成果,在網(wǎng)站ftp://garner.ucsd.edu 上,下 載了 SHAO、WUHN、TNML三個IGS站的同期觀測數(shù)據(jù),并應用上述方法對其進行了板塊運動改正、歷元轉(zhuǎn)換、框架轉(zhuǎn)換等處理,得到了3點的CGCS2000大地坐標。
通過與四等GPS控制網(wǎng)觀測數(shù)據(jù)進行聯(lián)合基線解算、平差計算,最終獲得了旦門山島GPS控制網(wǎng)的CGCS2000大地坐標??刂凭W(wǎng)的最弱邊相對中誤差為1/68000<1/45000,其他各項精度指標均滿足四等GPS控制網(wǎng)精度要求。
采用上述方法建立海島控制網(wǎng),不僅打破了傳統(tǒng)控制網(wǎng)布設方法的局限性,而且也避免了尋找、觀測測區(qū)周邊高等級控制點等工作。直接選取周邊IGS站點作為其起算數(shù)據(jù),方便、靈活、快捷,節(jié)省了大量的人力、物力,證明了采用該方法的高效性和可行性。
[1]陳俊勇.中國現(xiàn)代大地基準——中國大地坐標系統(tǒng)2000(CGCS2000)及其框架[J].測繪學報,2008,37(3):269 ~271.
[2]林曉靜,張小紅,郭斐.ITRF2005與CGCS2000坐標轉(zhuǎn)換方法與進度分析[J].大地測量與地球動力學,2010,30(2):117~119.
[3] 成英燕,程鵬飛,秘金鐘等.基于現(xiàn)框架下的省市級CORS站到CGCS2000的轉(zhuǎn)換[J].測繪通報,2011(7):1~3.
[4]程鵬飛,文漢江,成英燕等.2000國家大地坐標系橢球參數(shù)與GRS80和 WGS84的比較[J].測繪學報,2009,38(3):189~194.
[5]魏子卿.2000中國大地坐標系及其與WGS84的比較[J].大地測量與地球動力學,2008,28(5):1~5.
[6] 劉經(jīng)南,劉暉,鄒容等.建立全國CORS更新國家地心動態(tài)參考框架的幾點思考[J].武漢大學學報·信息科學版,2009,34(11):1261 ~1265.