□ 譚榮杰 陳長歡
(內(nèi)蒙古自治區(qū)測繪院,內(nèi)蒙古 呼和浩特 010051)
經(jīng)國務(wù)院批準,2000國家大地坐標系于2008年7月1日,在中國開始全面啟用。2012年12月中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室公布的《中國北斗衛(wèi)星導(dǎo)航空間信號接口控制文件》公開服務(wù)信號B1Ⅰ(1.0版)中要求,北斗衛(wèi)星導(dǎo)航系統(tǒng)采用2000國家大地坐標系。由此,中國進入了以2000國家大地坐標系為基準的測量時代。
作為中國新一代地心坐標系,2000國家大地坐標系的定義與國際地球參考系(ITRS)一致:原點位于地球質(zhì)心,也是參考橢球的幾何中心;Z軸指向國際地球自轉(zhuǎn)服務(wù)組織(IERS)定義的參考極(IRP)方向;X軸為IERS定義的參考子午面(IRM)與通過原點且同Z軸正交的赤道面的交線。2000國家大地坐標系參考橢球定義的基本常數(shù)為:長半軸 a=6378137.0m;扁率:f=1/298.257222101。參考橢球定義常數(shù),除地心引力常數(shù)(GM)外,同GRS80橢球一致。
其具體實現(xiàn)方法是通過中國地殼運動觀測網(wǎng)絡(luò),國家GPSA、B級網(wǎng),全國GPS一、二級網(wǎng),地殼運動監(jiān)測網(wǎng)4個網(wǎng)內(nèi)2600多個站點數(shù)據(jù),統(tǒng)一在ITRF97框架下聯(lián)合處理,得到2000.00歷元下的三維坐標。2000中國大地坐標系實際上是對ITRF97框架2000.00歷元的。依據(jù)以上定義,本文以內(nèi)蒙古連續(xù)運行參考站(NMGCORS)系統(tǒng)中部分站點數(shù)據(jù)為例,探討2000國家大地坐標的獲取方法,希望能夠得出一些有益的結(jié)論。
1.1 不同ITRF框架之間的變化關(guān)系
國際參考框架(ITRF)對每年分布全球的IGS連續(xù)運行站坐標與速度場處理分析,通過ITRS以年報與備忘錄的形式陸續(xù)發(fā)布了ITRF88、89、90、91、92、93、94、96、97、2000、2005、2008等參考框架,2017年ITRF2014框架近期即將公布,目前使用的還是ITRF2008參考框架,而作為2000國家大地坐標系基礎(chǔ)框架的ITRF97參考框架,是我們關(guān)注的重點。
理論上,將ITRF2008參考框架轉(zhuǎn)換到ITRF97參考框架,通過速度場將歷元歸算到2000.00即可獲得CGCS2000坐標。
對于ITRF2008參考框架轉(zhuǎn)換到ITRF97參考框架,可通過相似性變化轉(zhuǎn)換得到,公式1為七參數(shù)轉(zhuǎn)換模型:
式中(X、Y、Z)97為 ITRF97框架下的三維坐標,(X、Y、Z)08為 ITRF08 下的三維坐標,(Tx、Ty、Tz)為平移量,(Rx、Ry、Rz)為旋轉(zhuǎn)量,D為尺度因子。
對于不同的框架之間轉(zhuǎn)換,隨時間變化參數(shù)也隨之變化,可表達為參數(shù)P(t)從基準歷元P(t0)到P(t)的變化量,公式2為:
ITRF2008與ITFR97參考框架之間七參數(shù)及其變化率可在ITRF網(wǎng)站上查詢,(如表1所示)從ITRF2008與ITFR97參考框架的轉(zhuǎn)換參數(shù)。
表1 從ITRF2008與ITFR97參考框架的轉(zhuǎn)換參數(shù)
1.2 不同歷元之間的變換方法
對于待求點的ITRF97框架坐標來說,最終結(jié)果需要轉(zhuǎn)換到2000.00歷元。雖然IGS站都有精確的速度場,但是對于其他待求站點,無法獲得精確的速度場。較為容易實現(xiàn)解決辦法有兩個,一是利用待定點周邊的4個以上IGS站已知速度場模型,利用反距離加權(quán)的方法,求取待定點的速度場,通過公式3對不同歷元坐標進行變換。
式中,Xt1、Yt1、Zt1為待求點 2000.00 歷元坐標;Xt0、Yt0、Zt0。為待求點當(dāng)前歷元坐標;Vx、Vy、Vz為獲得的待求點的速度場?;蛘咧苯硬樵兾鹤忧涞热搜芯繃?000大地坐標系后得出的中國區(qū)域內(nèi)3°×3°的格網(wǎng)平均速度值;二是利用IGS站公共點當(dāng)前歷元與2000.00歷元不同坐標之間的相似性變換,求取待定點坐標。
選擇內(nèi)蒙古自治區(qū)連續(xù)運行參考站內(nèi)DONS、HLGE、KBAS、TZQ_ 四站及中國周邊 IGS站CHAN、SHAO、BJFS、LHAZ、ULAB、URUM、WUHN七站2016年在102、103天數(shù)據(jù),利用Gamit軟件進行數(shù)據(jù)處理,利用Globk進行平差,得到DONS、HLGE、KBAS、TZQ_ 四站在 ITRF08 框架下的2016.279歷元下的坐標,通過表1的給定轉(zhuǎn)換參數(shù),求解出 DONS、HLGE、KBAS、TZQ_ 四站在ITRF97框架下的2016.279歷元下的坐標,以上兩種方法將2016.279歷元轉(zhuǎn)化為2000.00歷元。
2.1 利用周邊IGS站已知速度場通過反距離加權(quán)的方法求取 DONS、HLGE、KBAS、TZQ_ 四站的速度場。其中西安站(XIAN)速度場由于站點維護原因,與周圍站點差異很大,故未使用?;蛘咄ㄟ^查閱魏子卿等公布的中國區(qū)域內(nèi)3°×3°的格網(wǎng)平均速度值,獲得 DONS、HLGE、KBAS、TZQ_四站的速度場值。見表2兩種方法所獲得速度場。根據(jù)公式3計算等到2000.00歷元坐標。
表2 兩種方法所獲得速度場
2.2 通過中國周邊 IGS站 CHAN、SHAO、BJFS、LHAZ、ULAB、URUM、WUHN 七個公共點,剔除掉具有顯著性差異站LHAZ及URUM,求得在2016.279歷元及2000.00歷元下七參數(shù)轉(zhuǎn)換關(guān)系,從而等到 DONS、HLGE、KBAS、TZQ_ 四站在2000.00歷元下的坐標。
根據(jù)定義,以上兩種方法獲得的均為2000國家大地坐標系坐標。以武漢大學(xué)解算結(jié)果為真值,與其比較,比較結(jié)果見表3兩種方法獲得坐標與武漢大學(xué)計算真值比較。
3.1 以上兩種方法獲得的2000國家大地坐標具有高度的相近和相似性,而七參數(shù)方法獲取結(jié)果略高于速度場方法,分析原因認為主要為缺少區(qū)域內(nèi)精確速度場模型,導(dǎo)致精度略低。但總體認為,兩種速度場的方法及七參數(shù)方法,都具有可靠性。
表3 兩種方法獲得坐標與武漢大學(xué)計算真值比較
3.2 兩種方法比較而言,X、Z方向精度低于Y方向,其原因待查。Z方向最大中誤差為±0.045,完全可以滿足區(qū)域2000國家大地坐標系平差計算要求,且方法簡單,不依賴于高精度的起算點,即可以獲得良好的結(jié)果。
3.3 對于高精度的求解,建議獲得更高精度的區(qū)域速度場;或通過板塊模型(NNR-NUVEL1A)對基線矢量進行轉(zhuǎn)換,再直接約束2000.00歷元的已知數(shù)據(jù)。
3.4 本文著重從2000國家大地坐標系定義角度,探討2000坐標的獲取方法,但與實際2000國家大地坐標一致性存在差異,使用時,視具體精度要求取舍。