莫 懦,韋廖軍,鄧 娟,林仁超

(南寧市勘察測繪地理信息院,廣西 南寧 530023)

0 引 言

NNCORS系統(tǒng)由5個連續(xù)運行參考站點組成,維持著南寧地區(qū)的測量框架和基準(zhǔn),隨著CGCS2000大地坐標(biāo)系的推廣,南寧地區(qū)正逐漸從以前坐標(biāo)系過渡到CGCS2000大地坐標(biāo)系,但NNCORS站點坐標(biāo)是基于ITRF2005框架,參考?xì)v元為瞬時歷元;而CGCS2000坐標(biāo)對應(yīng)著ITRF1997框架,參考?xì)v元為2000.0歷元;如何正確求解NNCORS站點的CGCS2000坐標(biāo),是南寧市測量工作者關(guān)心的一個問題。目前求解方法有兩種[1]∶1)基線平差法,即聯(lián)測南寧地區(qū)已知的CGCS2000坐標(biāo)成果的高級控制點,例如國家的A、B級點;2)框架轉(zhuǎn)換法,利用目前國際上高精度GPS后處理軟件解算出當(dāng)前最新ITRF框架下的坐標(biāo),再利用SOPAC網(wǎng)站公布的框架之間的轉(zhuǎn)換參數(shù)進行坐標(biāo)轉(zhuǎn)換,得出CGCS2000坐標(biāo)。本文分別論述了兩種方法,分析兩種方法的解算精度,并對兩種解算結(jié)果進行比較。

1 NNCORS概況

NNCORS系統(tǒng)共有5個連續(xù)運行參考站,布設(shè)在南寧市及周邊地區(qū),設(shè)為那樓站、那陳站、七塘站、金陵站和南寧站,站間最長距離69 km,最短距離30 km,平均間距48.6 km,數(shù)據(jù)處理中心設(shè)在南寧市。系統(tǒng)完全覆蓋南寧市六城區(qū),全天后向南寧六城區(qū)提供高精度、高時空分辨率、高效率、高覆蓋率的全球?qū)Ш叫l(wèi)星系統(tǒng)綜合信息服務(wù),維護著南寧地區(qū)的測量框架和基準(zhǔn)。

2 Gamit軟件介紹

GAMIT/GLOBK是美國麻省理工學(xué)院(MIT)與斯克里普斯海洋研究所(SIO)共同研發(fā)的高精度GPS定位定軌軟件,并得到美國哈佛大學(xué)和美國國家科學(xué)基金會的支持,目前已廣泛用于大尺度、高精度GPS定位和地殼形變分析。Gamit解算參數(shù)估計采用最小二乘算法,進行雙差觀測量處理,完全消除了接收機和衛(wèi)星鐘差的影響,削弱了軌道,大氣折射等系統(tǒng)性誤差,其對于長基線解算精度可達到10-7～10-9量級[2-3]。

3 NNCORS站點CGCS2000坐標(biāo)的確定

為了分析兩種求解結(jié)果差異性,本文選取了同一時間段的NNCORS站點數(shù)據(jù),采用相同的GAMIT數(shù)據(jù)處理策略進行基線處理。

南寧地區(qū)距離北部灣海岸300多公里,數(shù)據(jù)處理策略中不考慮海潮模型改正,并從極潮中移去平均值;考慮固體潮改正、依賴于頻率的改正、極潮改正,所有數(shù)據(jù)采用間隔均為30″,并用IGS的精密星歷進行逐日處理。

3.1 基線平差法

基線平差法基本思想是將NNCORS站點聯(lián)測周邊具有CGCS2000精確坐標(biāo),采用GAMIT軟件進行基線處理,再提取GAMIT處理的基線文件,用武漢大學(xué)研制的COSAGPS進行平差, 三維約束平差采用具有CGCS2000精確坐標(biāo)的高級點進行,進而獲得NNCORS站點的精確CGCS2000坐標(biāo)。根據(jù)NNCORS站點分布情況,本文選取了8個南寧周邊精度指標(biāo)較好的高級點,與NNCORS站點進行24 h以上的聯(lián)測。將聯(lián)測數(shù)據(jù)采用GAMIT進行基線處理,采用表1的數(shù)據(jù)處理策略對數(shù)據(jù)進行處理,數(shù)據(jù)處理后,查看處理文件中的Q文件中的均方根殘差(nrms),nrms是衡量基線解算質(zhì)量的重要指標(biāo)之一,nrms值要求小于0.5,如果nrms值大于0.5,則認(rèn)為數(shù)據(jù)中有大周跳沒有剔除,需要進行重新處理或者使用CVIEW進行人工剔除，本例中nrms=0.187,滿足要求。

表1 數(shù)據(jù)處理策略

提取GAMIT處理結(jié)果文件中的Q-FILE或者0-FILE基線解算結(jié)果文件,GAMIT對基線進行了兩次迭代處理,第二次迭代處理的基線數(shù)據(jù)是平差采用的數(shù)據(jù),利用COSAGPS軟件對GAMIT處理的基線結(jié)果進行網(wǎng)平差,三維約束網(wǎng)平差控制點采用高級點的CGCS2000坐標(biāo)進行約束,得出NNCORS站點的CGCW2000坐標(biāo)。網(wǎng)平差解算結(jié)果如表2所示。

表2 網(wǎng)平差邊長相對精度

表3 NNCORS站點CGCS2000坐標(biāo)及精度

限于篇幅,表3只示出了部分邊長及最弱邊相對精度,其中JZ01-NANN為最弱邊,相對精度達到了2.32×10-7,其他邊長相對精度均小于此值,坐標(biāo)量中誤差均在毫米級,達到了國家GPS B級網(wǎng)的精度要求。

3.2 框架轉(zhuǎn)換法

ITRF國際地球參考框架是IERS的具體實現(xiàn),是一個地心四維坐標(biāo)參考框架,是基于VBVI、LLR、SSR、GPS和DORIS等空間技術(shù)所建立起來的,是全球IGS站和速度場的具體實現(xiàn)。目前最新的ITRF框架為ITRF2008參考框架,最新精密星歷計算的坐標(biāo)為ITRF2008框架下的坐標(biāo),而CGCS2000對應(yīng)的框架為ITRF97,Epoch=2 000.0,ITRF2008坐標(biāo)轉(zhuǎn)換到CGCS2000坐標(biāo),實質(zhì)就是框架之間的轉(zhuǎn)換。然而ITRF2008到ITRF97框架之間沒有直接的轉(zhuǎn)換參數(shù),根據(jù)ITRF所公布的框架轉(zhuǎn)換參數(shù),可根據(jù)ITRF2008與ITRF2000、ITRF2000與ITRF97的轉(zhuǎn)換關(guān)系得到。表4示出了ITRF2008與ITRF2000、ITRF2000與ITRF97的轉(zhuǎn)換參數(shù)及參數(shù)變率。

表4 ITRF2008至各框架轉(zhuǎn)換參數(shù)及其速率

表4所列數(shù)據(jù)為特定歷元下的轉(zhuǎn)換參數(shù),而當(dāng)前求解的ITRF2008框架下的坐標(biāo)與CGCS2000參考?xì)v元跨度達10年以上,需要考慮板塊漂移和地殼形變等因素。地球表面的觀測點會隨時間的推移進行漂移及升降,不同歷元基準(zhǔn)的坐標(biāo)應(yīng)進行歷元基準(zhǔn)的歸算,公式為[4]

(1)

公式中的IGS測站的速度為年變率。t1為框架歷元,t2為需轉(zhuǎn)換到的歷元,這里為CGCS2000所在的歷元,即2 000.0.為了獲取準(zhǔn)確的CGCS200精確坐標(biāo),必須獲得精度為毫米級的點位速度場,而南寧地區(qū)目前沒有CORS站維持本地區(qū)的速度場模型,所以實際操作中南寧地區(qū)速度場可以采用加權(quán)平均法[5]和NNR-NUVEL-1A板塊運動模型[6]所獲得,南寧地區(qū)上某一點的速度公式可以采用

(2)

選取了2013年,年積日為280,…,289十天的NNCORS站點數(shù)據(jù),選取CUSV、PIMO、LHAZ、SHAO、TWTF、BJFS、SUWN的ITRF2008框架下的坐標(biāo)做約束,采用IGS站精密星歷進行處理,處理后得到單天的松弛解,再用GLOBK軟件單日解合并,平差得到ITRF2008框架,Ecoph=2 013.781下的精確坐標(biāo),單日解水平精度5 mm,垂直精度10 mm.同時根據(jù)以上框架轉(zhuǎn)換和歷元基準(zhǔn)之間的轉(zhuǎn)換關(guān)系,編制了程序,計算出了NNCORS的CGCS2000坐標(biāo)，坐標(biāo)如表5所示。

表5 NNCORS CGCS2000坐標(biāo)

用基線平差法與框架轉(zhuǎn)換法計算出的NNCORS站點的CGCS2000坐標(biāo)作了比較,比較結(jié)果如表6所示,比較結(jié)果表明,各站點坐標(biāo)分量較差大部分都小于1 cm,達到毫米級。最大較差絕對值1.6 cm,最小較差絕對值為0.5 cm.

表6 CGCS2000坐標(biāo)較差

4 結(jié)束語

本文闡述的基線平差法和框架轉(zhuǎn)換法計算NNCORS的CGCS2000坐標(biāo)保持了高度的一致,但在計算過程中需注意以下幾個方面:

1)城市CORS站點CGCS2000精確坐標(biāo)的確定,在國家A、B級坐標(biāo)點密度分布不均勻的地區(qū)框架轉(zhuǎn)換法比基線平差法有優(yōu)勢。

2)IGS站的選取應(yīng)兼顧平衡性與精度一致性原則;

3)不同框架,不同歷元下的坐標(biāo)轉(zhuǎn)換,其點位速度場起關(guān)鍵性作用,因此需獲得轉(zhuǎn)換點的精確速度場,否則會影響轉(zhuǎn)換精度;

4)用基線平差法時,應(yīng)選擇精度高、穩(wěn)定性好的高級點作為聯(lián)測點。

[1]成英燕,2000國家大地坐標(biāo)框架轉(zhuǎn)換方法[R].2009.

[2]李 杰, 王曉強. 應(yīng)用GAMIT軟件進行高精度GPS定位解算[J].內(nèi)陸地震, 2001, 15 (4): 355-359.

[3]鄂棟臣, 詹必偉, 姜衛(wèi)平,等. 應(yīng)用GAMIT/GLOBK軟件進行高精度GPS數(shù)據(jù)處理[J].極地研究,2005, 17 (3): 173-182.

[4]劉 立,成英燕.坐標(biāo)框架轉(zhuǎn)換若干問題的研究[J].全球定位系統(tǒng),2010,35(1):20-25.

[5]林曉靜,張小紅,郭 斐.ITRF2005與CGCS2000坐標(biāo)轉(zhuǎn)換方法與精度分析[J].大地測量與地球動力學(xué),2010(2):117-120,124.

[6]榮 敏,周 魏.基于ITRF2005的全球版塊運動模型[J].海洋測繪,2009(1):1-3