GAMIT/GLOBK在高精度長(zhǎng)基線像控點(diǎn)解算中的應(yīng)用

楊德芳 丁開華 許偉

摘 要：青海省藏區(qū)大比例尺基礎(chǔ)地理信息都蘭測(cè)區(qū)數(shù)據(jù)采集項(xiàng)目，航攝結(jié)束后攝區(qū)像控點(diǎn)布設(shè)南北基線跨度160 km，東西基線跨度135 km，基線跨度很大，如果采用傳統(tǒng)RTK方法，受基線距離的限制，測(cè)區(qū)必須重新建設(shè)多個(gè)基準(zhǔn)站，外業(yè)工作效率將大幅度降低，并且常規(guī)坐標(biāo)解算方法很難保障解算精度。為了解決長(zhǎng)基線像控點(diǎn)高精度解算問(wèn)題，利用GAMIT/GLOBK解算方式，引入了9個(gè)IGS基準(zhǔn)站、精密星歷及精密鐘差對(duì)測(cè)區(qū)106個(gè)靜態(tài)觀測(cè)的像控點(diǎn)進(jìn)行坐標(biāo)解算，用另外一組以攝區(qū)高精度CORS（連續(xù)運(yùn)行基準(zhǔn)站）站作為控制點(diǎn)，區(qū)域網(wǎng)平差后cm級(jí)的解算成果檢驗(yàn)GAMIT/GLOBK的解算精度。通過(guò)分析發(fā)現(xiàn)，采用GAMIT/GLOBK解算模式，在不依賴攝區(qū)CORS站及已知控制點(diǎn)的情況下，其基線解算所有單天的NRMS標(biāo)準(zhǔn)化均方根值精度值均小于0.25，基線重復(fù)率相對(duì)重復(fù)性精度達(dá)到了10-3 m、絕對(duì)重復(fù)性精度達(dá)到了10-8 m，通過(guò)TBC成果對(duì)GLOBK坐標(biāo)平差精度進(jìn)行評(píng)定，得出N，E，U等3方向的中誤差分別為：11.9，10.4，4.2 cm，通過(guò)利用距測(cè)區(qū)150 km范圍內(nèi)的3個(gè)已知CORS站點(diǎn)對(duì)上述N，E方向中的系統(tǒng)性共模誤差進(jìn)行改正，改正后中誤差為-0.728 4，-1.995 cm.解算精度滿足1∶1萬(wàn)外業(yè)像控點(diǎn)要求，說(shuō)明利用GAMIT/GLOBK解算方式，可以為青海省大部分無(wú)基準(zhǔn)、無(wú)控制的地區(qū)開展GPS像控點(diǎn)數(shù)據(jù)測(cè)量及解算工作提供實(shí)踐精度參考。

關(guān)鍵詞：長(zhǎng)基線像控點(diǎn)解算;GAMIT/GLOBK;CORS區(qū)域網(wǎng)平差;精度評(píng)定

中圖分類號(hào)：P 228.4 ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1672-9315（2019）05-0849-10

Abstract：For the project of basic geographic information data acquisition in Dulan area of Qinghai Tibetan，the span of baselines formed by photo-control-point is quite large， with 160 km from north to south and 135 km from east to west.If the traditional RTK method is used，some reference stations must be built in the survey area to overcome the effects of baseline distance.Thus，the efficiency of field work will be greatly reduced，and the accuracy of coordinate calculation is difficult to be guaranteed.In order to solve the problem about the high-precision positioning for long baseline photo-control-point，this paper used the GAMIT/GLOBK mode，while nine IGS reference stations，precise ephemeris and precise clock bias were introduced，to calculate the photo-control-point for 106 static observations in the survey area.Using another set of high accuracy Continuously Operating Reference Stations（CORS） stations as control points，to check the accuracy of GAMIT/GLOBK mode by the cm level solution results after block adjustment.Based on the analysis，we found the baseline accuracy of all single day Normalized Root-Mean-Square value（NRMS） is less than 0.25 by using GAMIT/GLOBK mode without depending on the CORS station and the known control points.The relative repeatability accuracy of baseline repetition rate has reached 10-3 m，absolute repeatability accuracy has reached 10-8 m.By evaluating the accuracy of GLOBK coordinate adjustment through TBC results，we concluded the mean error of N，E and U in three directions is 0.119，0.104，0.042 m respectively.The systematic common-mode errors in the N and E directions are corrected by using three known CORS stations within the range of 150 km from the measuring area.The standard deviation is -0.728 4，-1.995 cm after correction.The accuracy of the solution meets the requirements of 1∶10 000 external photo-control-point，It means the GAMIT/GLOBK mode，can provide a reference for the GPS photo-control-point data measurement and calculation in most areas of Qinghai Province without datum and control.

Key words：long baseline image control point calculation;GAMIT/GLOBK;CORS regional network adjustment;precision evaluation

0 引 言

GAMIT/GLOBK組合平差解算方式，通過(guò)利用全球分布的IGS臺(tái)站觀測(cè)文件、GPS衛(wèi)星精密星歷、精密鐘差、對(duì)流層、電離層等各類改正模型可以獲取測(cè)站高精度的三維坐標(biāo)信息。其解算過(guò)程中采用IGS站作為已知控制站點(diǎn)，因此待解算測(cè)站點(diǎn)可以不考慮基線距離的限制、區(qū)域網(wǎng)站點(diǎn)間共模誤差的影響，可以直接實(shí)現(xiàn)待解算測(cè)站精密單點(diǎn)解算，大大解決了常規(guī)GPS測(cè)量依賴公共站點(diǎn)、短基線解算的缺點(diǎn)。

目前對(duì)于GAMIT/GLOBK解算方式的研究，紀(jì)海源通過(guò)GAMIT/GLOBK軟件對(duì)汾渭斷陷帶及鄰近地區(qū)2011—2013年的36個(gè)連續(xù)站和227個(gè)區(qū)域站GPS觀測(cè)數(shù)據(jù)進(jìn)行處理[1]，其2次解算的GPS連續(xù)站水平方向平均點(diǎn)位誤差在2 mm左右，垂直分量平均點(diǎn)位誤差在5 mm左右，一般點(diǎn)位誤差的限差為5～10 mm.林超才利用GAMIT/GLOBK對(duì)四川地震監(jiān)測(cè)網(wǎng)CORS站及國(guó)內(nèi)部分IGS站2008年5月12日地震前后一段時(shí)間的觀測(cè)數(shù)據(jù)按照地震前及地震后分別進(jìn)行單日的基線解算[2]，并對(duì)基線解算結(jié)果進(jìn)行精度評(píng)定和原因分析。劉宗強(qiáng)等對(duì)進(jìn)行陸態(tài)網(wǎng)解算時(shí)只需要加入北半球的IGS站就可使基線得到較好的U方向解算精度，陸態(tài)網(wǎng)選擇8個(gè)左右的IGS站較為合適[3]。高旺等實(shí)驗(yàn)結(jié)果表明對(duì)于中長(zhǎng)基線，使用IGR或IGU星歷進(jìn)行GAMIT基線解算與使用IGS最終星歷解算的坐標(biāo)分量較差優(yōu)于10-9數(shù)量級(jí)[4]。楊登科等從L1_ONLY和L1L2_INDEPENDENT僅適用5 km以下基線向量解算，LC_HELP和LC_AUTCLN的解算效果相當(dāng)，LC_AUTCLN相對(duì)較好[5]。羅全等通過(guò)GAMIT軟件解算不同截止高度角條件下的基線結(jié)果，得出當(dāng)截止高度角取10°～15°時(shí)，基線解算結(jié)果最佳[6]。

由于上述優(yōu)勢(shì)，攝區(qū)像控點(diǎn)布設(shè)及施測(cè)南北基線跨度160 km，東西基線跨度135 km，點(diǎn)與點(diǎn)間基線跨度大，采用GAMIT/GLOBK組合平差模式長(zhǎng)基線解算時(shí)必須考慮GPS長(zhǎng)基線解算時(shí)對(duì)流層、電離層及地球曲率等因素的影響，通常情況下可以考慮2種解算模式。

第1種：測(cè)區(qū)待測(cè)點(diǎn)與覆蓋測(cè)區(qū)范圍CORS基準(zhǔn)站聯(lián)測(cè)，采用商用軟件進(jìn)行平差解算;

第2種：采用測(cè)區(qū)周邊均勻分布的全I(xiàn)GS站作為基準(zhǔn)站，利用GAMIT/GLOBK進(jìn)行長(zhǎng)基線平差解算。

第1種模式必須依賴高精度CORS站作為控制點(diǎn)，如果CORS站無(wú)法覆蓋整個(gè)測(cè)區(qū)或是測(cè)區(qū)內(nèi)沒(méi)有有效的CORS站，則需聯(lián)測(cè)周邊已知控制點(diǎn)重新布設(shè)控制網(wǎng)進(jìn)行測(cè)量，這樣就增加了外業(yè)工作量及工作效率，第2種模式不受測(cè)區(qū)CORS站及已知控制點(diǎn)的制約，靈活采用全球分布的IGS站作為控制點(diǎn)，人員到達(dá)測(cè)區(qū)可直接開始靜態(tài)施測(cè)，無(wú)需考慮點(diǎn)間基線長(zhǎng)短，點(diǎn)、邊連接等測(cè)量方法，就可以實(shí)現(xiàn)cm級(jí)的精度，文中測(cè)區(qū)站點(diǎn)及IGS測(cè)站分布情況及測(cè)站日觀測(cè)情況分別如圖1，2所示。

1 GAMIT/GLOBK解算策略

1.1 GAMIT/GLOBK解算意義

由于青海省獨(dú)特的高原地理環(huán)境，目前建成的CORS站還不能覆蓋該省全境，如圖3所示，紅色圓形區(qū)域代表每一個(gè)CORS站所覆蓋的范圍，若要在CORS無(wú)覆蓋的測(cè)區(qū)范圍內(nèi)開展GPS測(cè)量工作，則存在重新布設(shè)基準(zhǔn)控制點(diǎn)及聯(lián)測(cè)問(wèn)題[7]，此過(guò)程涉及大量人力物力支出問(wèn)題，如果采用GAMIT/GLOBK解算方式，就可以解決測(cè)區(qū)無(wú)起算控制點(diǎn)的限制。

1.2 GAMIT/GLOBK表文件及平差解算

GAMIT/GLOBK在數(shù)據(jù)解算之前，首先要根據(jù)測(cè)站的觀測(cè)時(shí)間對(duì)“tables表文件”進(jìn)行配置，表文件主要包括系統(tǒng)自帶的表文件和待更新的表文件2類[8-9]。由于文中觀測(cè)數(shù)據(jù)是2017年，在處理之前分別對(duì)系統(tǒng)自帶的表文件，系統(tǒng)定期更新的表文件及GAMIT/GLOBK控制表文件進(jìn)行了更新，見表1，表2，表3.

1.2.1 系統(tǒng)自帶的主要表文件

系統(tǒng)自帶的主要表文件更新內(nèi)容包括：天線高及相位中心偏移模式參數(shù)，大地水準(zhǔn)面參數(shù)，海潮文件參數(shù)，周跳自動(dòng)探測(cè)和修改命令參數(shù)，P1-C1，P1-P2碼差等參數(shù)的更新。具體參數(shù)見表1.

1.2.2 GAMIT/GLOBK平差解算

分別采用分步處理與批處理方式進(jìn)行解算，具體解算流程如圖4所示。

3 解算參數(shù)配置

1）process.defaults中修改“set sint”為5 s，“set nepc”為2 160 s，set.aprf=control.apr;

2）sestbl.中修改基線類型為L(zhǎng)C_HELP，篩檢因子Decimation Factor=1，Interval zen=2，衛(wèi)星高度角改為15°，Use met.grid=N，Use otl.grid=N，Use atml.grid=N，Tides applied = 23;

3）station.info手動(dòng)更新待解算點(diǎn)測(cè)站信息，天線高改化到天線相位中心LEPHC;

4）sittbl.中修改控制點(diǎn)約束平面0.02，高程005，待求點(diǎn)約束平面0.05，高程0.05;

5）lfile.中對(duì)于IGS站坐標(biāo)從itrf08.apr中提取，對(duì)于待解算站坐標(biāo)利用雙差定位法提取;

6）GLOBK_comb.cmd中禁用eq_file，增加make_svs@.svs Z，改變prt_opt GDLF CMDS，改變org_opt PSUM CMDS GDLF，啟用所有控制點(diǎn)參與平差use_site all;

7）glorg_comb.cmd修改stab_it網(wǎng)平差迭代次數(shù)為5次，禁用source.

4 解算精度評(píng)定

GAMIT/GLOBK基線解算及平差結(jié)束后生成結(jié)果[15-16]主要包括：H-file基線松弛解（h0526a.17146），O-file約束解（o0526a.146），Q-file過(guò)程記錄文件（q0526a.146），GLOBK坐標(biāo)平差文件（GLOBK_0524.org），針對(duì)上述文件進(jìn)行精度分析。

4.1 標(biāo)準(zhǔn)化均方根誤差NRMS精度統(tǒng)計(jì)

解算完后，以單天解給出了18個(gè)年積日的NRMS精度值，從表4可以看出，最大的為0.193 29E+00，最小的為0.180 08E+00，根據(jù)國(guó)內(nèi)外025的數(shù)據(jù)處理經(jīng)驗(yàn)[17]，所有單天的NRMS精度值均小于0.25，說(shuō)明單時(shí)段解算出的基線值偏離其加權(quán)平均值的程度較小，處理過(guò)程中周跳得到了較好地修復(fù)，基線估算精度較高，解算質(zhì)量指標(biāo)符合要求。

4.2 天空?qǐng)D統(tǒng)計(jì)

LC phase Residual[18]表示接收機(jī)接受信號(hào)的好壞，由于測(cè)站較多，以測(cè)站JB70為例給出解釋，如圖5所示，該測(cè)站信號(hào)整體趨勢(shì)在-20～+20 mm之間圍繞藍(lán)色收斂線波動(dòng)，收斂趨勢(shì)平順，并且逐漸接近于0，測(cè)站RMS值為9.1 mm，說(shuō)明解算過(guò)程中天線相位模型引用正確，信號(hào)沒(méi)有受到較強(qiáng)的電磁干擾，收斂情況較好。

4.3 基線重復(fù)率

以單日解解算出了18個(gè)年積日的基線重復(fù)率，從表5，圖6可以看出，絕對(duì)重復(fù)性最大值為0006 814 51 m，最小的為0.000 309 02 m，精度在mm級(jí)，相對(duì)重復(fù)性最大值為1.000E-08，最小的為-4.736E-10，相對(duì)精度達(dá)到了10-8，可見，基線解算結(jié)果精度較高。

4.4 GAMIT/GLOBK精度檢驗(yàn)

對(duì)測(cè)區(qū)106個(gè)像控點(diǎn)采用GAMIT/GLOBK進(jìn)行解算合格后的結(jié)果，直接采用魏子卿研究員[20]關(guān)于“中國(guó)大陸速度場(chǎng)”一文中給出的中國(guó)大陸3°×3°格網(wǎng)平均速度。將ITRF08框架下得到的空間直角坐標(biāo)，換算到CGCS 2000坐標(biāo)系。轉(zhuǎn)換過(guò)程中分別利用公式（8）、公式（9）將2017.* * *歷元轉(zhuǎn)換到2 000.000歷元，將ITRF 2008框架坐標(biāo)轉(zhuǎn)換到ITRF 97框架坐標(biāo)，這樣就得到了CGCS 2000坐標(biāo)[19]。

為了驗(yàn)證GAMIT/GLOBK解算的最終點(diǎn)位三維坐標(biāo)，用另外一組以攝區(qū)高精度CORS站作為控制點(diǎn)，平差后cm級(jí)的解算成果對(duì)文中的解算精度進(jìn)行檢驗(yàn)。由上圖7可知，106個(gè)測(cè)站N方向殘差值在0.12 m附近波動(dòng)，E方向殘差值在0.10 m附近波動(dòng)，U方向殘差值在0.04 m附近波動(dòng)，經(jīng)公式（6）計(jì)算得出RMS N=0.119 m，RMS E=0.104 m，RMS U=0.042 m，U方向解算精度達(dá)到了cm級(jí)，平面方向解算精度在dm級(jí)。

4.5 CORS站偏差改正

4.4中平面方向出現(xiàn)dm級(jí)的解算精度初步估計(jì)，是由于受測(cè)區(qū)速度場(chǎng)模型精度影響。圖8應(yīng)用GMT[21]軟件，結(jié)合青藏高原東北緣水平運(yùn)動(dòng)GPS速度場(chǎng)[22-23]，給出了測(cè)區(qū)基于歐亞板塊的青海省GPS速度場(chǎng)，圖中紅色矩形框?yàn)槲闹袦y(cè)區(qū)范圍，從圖可以看出，在測(cè)區(qū)范圍內(nèi)只有一個(gè)已知速度變化的測(cè)站，缺少起算點(diǎn)數(shù)據(jù)，無(wú)法應(yīng)用局域歐拉矢量法內(nèi)插速度，所以直接引用魏子卿速度場(chǎng)模型。測(cè)區(qū)所有站點(diǎn)都采用統(tǒng)一的速度場(chǎng)模型勢(shì)必會(huì)引入系統(tǒng)性的共模誤差[24]的影響，即在坐標(biāo)分量上表現(xiàn)出整體性偏移，為了消除這一誤差的影響，采用距離測(cè)區(qū)150 km范圍內(nèi)的3個(gè)省級(jí)CORS站進(jìn)行共模誤差的提取及分析。

根據(jù)圖3結(jié)合測(cè)區(qū)情況，選取DGLE，YANH，GERM等3個(gè)CORS測(cè)站，利用同一Gamit/GLOBK參數(shù)解算模型，對(duì)CORS站點(diǎn)進(jìn)行解算，將解算轉(zhuǎn)換后的最終結(jié)果與CORS已知測(cè)站成果進(jìn)行對(duì)比分析，得出CORS站殘差為RMS N=0.126 m，RMS E=0.124 m（上述解算中高程精度已滿足要求，不再進(jìn)行改正分析），故將這一值，確定為本測(cè)區(qū)共模誤差改正值，將這一改正值改正到測(cè)區(qū)所有像控點(diǎn)坐標(biāo)值中。

從圖7也可以看出，所有像控點(diǎn)N，E方向表現(xiàn)出整體性偏移，而且該偏移值在0.119 m與0104 m之間浮動(dòng)，浮動(dòng)范圍±0.002 m，如果將CORS站殘差值改正至上述測(cè)站整體偏移中，那么所有測(cè)站N，E方向精度將大幅度提高。圖9是根據(jù)CORS站殘差改正后的測(cè)站殘差，從圖10可以看出，經(jīng)過(guò)改正后，所有測(cè)站N，E方向殘差值均有所減小，N方向殘差在-0.5 cm至2 cm之間浮動(dòng)，E方向殘差在-3 cm至-0.6 cm之間浮動(dòng)，經(jīng)過(guò)改正后的測(cè)站RMS N=-0.7284 cm，RMS E=-1995 cm，RMS U=0.042 m，說(shuō)明經(jīng)過(guò)CORS站殘差改正后水平方向精度達(dá)到了亞厘米級(jí)，對(duì)于1∶1萬(wàn)的像控精度，上述解算精度已滿足要求。

5 結(jié) 論

1）對(duì)解算策略及解算精度指標(biāo)進(jìn)行的分析發(fā)現(xiàn)：GAMIT/GLOBK基線解算所有單天的NRMS精度值均小于0.25，基線重復(fù)率相對(duì)重復(fù)性精度達(dá)到了mm級(jí)、絕對(duì)重復(fù)性精度達(dá)到了10-8;通過(guò)TBC成果對(duì)GLOBK坐標(biāo)平差精度進(jìn)行評(píng)定得出，解算結(jié)果RMS N=0.119 m，RMS E=0.104 m，RMS U=0042 m，N，E方向所有點(diǎn)位均表現(xiàn)出整體性偏移，并且該偏移值均值分別為11.9，10.4 cm.

2）為了解決這一偏移的影響，通過(guò)采用距離測(cè)區(qū)150 km范圍內(nèi)的3個(gè)省級(jí)CORS站進(jìn)行該整體性偏移量的提取，得出CORS站整體性偏移改正量殘差為RMS N=0.126 m，RMS E=0.124 m，將這一改正量改正到測(cè)區(qū)所有像控點(diǎn)坐標(biāo)值中。經(jīng)過(guò)改正后的所有測(cè)站RMS N=-0.728 4 cm，RMS E=-1.995 cm，RMS U=0.042 m，說(shuō)明經(jīng)過(guò)CORS站殘差改正后水平方向精度達(dá)到了厘米級(jí)，高程精度達(dá)到了厘米級(jí)。解算結(jié)果滿足文中1∶1萬(wàn)外業(yè)像控點(diǎn)要求，為青海大部門無(wú)基準(zhǔn)、無(wú)控制的地區(qū)進(jìn)行GPS數(shù)據(jù)解算工作的可行性提供了實(shí)踐參考，達(dá)到了在大量減少外業(yè)工作量的同時(shí)控制成本的目的。

參考文獻(xiàn)（References）：

[1] 紀(jì)海源.基于GAMIT/GLOBK與QOCA的汾渭斷陷帶地殼運(yùn)動(dòng)研究[D].西安：西安科技大學(xué)，2014.

JI Hai-yuan.Study of crustal movement in Fenwei rift zone based on GAMIT/GLOBK and QOCA[D].Xi’an：Xi’an University of Science and Technology，2014.

[2] 林超才.基于GAMIT/GLOBK與QOCA的汶川震區(qū)地殼形變分析與研究[D].成都：西南交通大學(xué)，2012.

LIN Chao-cai.Analysis and research of Wenchun earthquake zone based on GAMIT/GLOBK and QOCA[D].Chengdu：Southwest Jiaotong University，2012.

[3] 劉宗強(qiáng)，黨亞民，何 濤，等.IGS站的選取對(duì)陸態(tài)網(wǎng)解算精度的影響[J].測(cè)繪工程，2017，26（9）：28-31.

LIU Zong-qiang，DANU Ya-min，HE Tao，et al.Selection of IGS station in solution precision of CMONOC[J].Engineering of Surveying and Mapping，2017，26（9）：28-31.

[4] 高 旺，高成發(fā)，潘樹國(guó)，等.基于快速星歷的GAMIT高精度基線解算研究[J].測(cè)繪科學(xué)，2015，40（2）：22-25.

GAO Wang，GAO Cheng-fa，PAN Shu-guo，et al.Research on GAMIT high-precision baseline calculation based on rapid ephemeris[J].Surveying of Surveying and Mapping，2015，40（2）：22-25.

[5]楊登科，安向東.基于GAMIT的GPS基線解類型分析[J].測(cè)繪地理信息，2016，41（5）：25-28.

YANG Deng-ke，AN Xiang-dong.Analysis of GPS baseline solution type based on GAMIT[J].Journal of Geomatics，2016，41（5）：25-28.

[6]羅 權(quán).不同截止高度角對(duì)GAMIT基線解算的影響分析[J].測(cè)繪地理信息，2017，42（3）：14-17.

LUO Quan.Effect of different cut-off elevation angles on GAMIT baseline solution[J].Journal of Geomatics，2017，42（3）：14-17.

[7]楊鴻海，竇 超，邊雪清.D-TIN在青海海東地區(qū)CORS基準(zhǔn)站組網(wǎng)中的應(yīng)用[J].測(cè)繪通報(bào)，2012（11）：65-67.

YANG Hong-hai，DOU Chao，BIAN Xue-qing.Application of D-TIN to CORS fiducial networking of Qinghai Haidong region[J].Bulletin of Surveying and Mapping，2012（11）：65-67.

[8]徐 杰，孟 黎，唐詩(shī)華.GAMIT/GLOBK批處理在高精度海量數(shù)據(jù)處理中的應(yīng)用[J].測(cè)繪科學(xué)，2008，33（6）：187-188.

XU Jie，MENG Li，TANG Shi-hua.Application of GAMIT/GLOBK batch processing in high-precision mass data processing[J].Science of Surveying and Mapping，2008，33（6）：187-188.

[9]Li X，Xu L，F(xiàn)ang Y，et al.Estimation of the precipitable water vapor from ground-based GPS with GAMIT/GLOBK[C]//Lita International Conference on Geoscience & Remote Sensing，2010：210-214.

[10]占 偉，武艷強(qiáng)，劉志廣.GAMIT/GLOBK數(shù)據(jù)自處理程序的應(yīng)用[J].測(cè)繪科學(xué)，2009（S1）：171-173.

ZHAN Wei，WU Yan-qiang，LIU Zhi-guang.Application of GAMIT/GLOBK data self-processor[J].Science of Surveying and Mapping，2009（S1）：171-173.

[11]任秀波，尹偉言，陳小英，等.GPS基線質(zhì)量評(píng)價(jià)[J].測(cè)繪技術(shù)裝備，2017，19（2）：84-85.

RENG Xiu-bo，YIN Wei-yan，CHEN Xiao-ying，et al.Baseline quality evaluation about GPS[J].Geomatics Technology and Equipment，2017，19（2）：84-85.

[12]張 祥.區(qū)域CORS站穩(wěn)定性監(jiān)測(cè)分析[D].西安：西安科技大學(xué)，2016.

ZHANG Xiang.Stability monitoring and analysis of regional CORS station[D].Xi’an：Xi’an University of Science and Technology，2016.

[13]楊登科，安向東，湯 勰.GAMIT數(shù)據(jù)處理中基線解算模式的對(duì)比分析[J].測(cè)繪地理信息，2016，41（2）：18-21.

YANG Deng-ke，AN Xiang-dong，TANG Xie.Comparative analysis of baseline solution based on GAMIT data processing[J].Journal of Geomatics，2016，41（2）：18-21.

[14]Zhao G，Ren L.Construction of GPS data processing platform based on GAMIT/GLOBK[J].Journal of Geomatics，2016.

[15]Jie X U，Meng L，Tang S H.Application of the GAMIT/GLOBK batch processing in high precision massive data processing[J].Science of Surveying & Mapping，2008，33（6）：187-189.

[16]于 強(qiáng)，易長(zhǎng)榮，占 惠.ITRF2000轉(zhuǎn)換到CGCS2000框架的分析[J].全球定位系統(tǒng)，2009，34（5）：49-51.

YU Qiang，YI Chang-rong，ZHAN Hui.Analysis from ITRF2000 to CGCS2000 framework[J].Gnss World of China，2009，34（5）：49-51.

[17]Herring T A，Kink R W，McClusky S C.GAMIT reference manual（Release 10.4）[R].Department of Earth，Atmospheric，and Planetary Sciences Massachusetts Institute of Technology，2010.

[18]Pohl M，Schenke H W.Precise baseline determination using GAMIT/GLOBK Software[J].Phytochemistry，1997，84（6）：47-55.

[19]蔣志浩，張 鵬，秘金鐘，等.基于CGCS2000的中國(guó)地殼水平運(yùn)動(dòng)速度場(chǎng)模型研究[J].測(cè)繪學(xué)報(bào)，2009，38（6）：471-476.

JIANG Zhi-hao，ZHANG Peng，BI Jin-zhong，et al.The model of crustal horizontal movement based on CGCS2000 frame[J].Acta Geodetica et Cartographical Sinica，2009，38（6）：471-476.

[20]魏子卿，劉光明，吳富梅.2000中國(guó)大地坐標(biāo)系——中國(guó)大陸速度場(chǎng)[J].測(cè)繪學(xué)報(bào)，2011，40（4）：403-410.

WEI Zi-qing，LIU Guang-ming，WU Fu-mei.China geodetic coordinate system 2000：velocity field in mainland China[J].Acta Geodaetica et Cartographica Sinica，2011，40（4）：403-410.

[21]占 偉，孟憲剛，劉志廣.GMT繪制GPS速度場(chǎng)的應(yīng)用[J].華北地震科學(xué)，2010，28（3）：61-64.

ZHAN Wei，MENG Xian-gang，LIU Zhi-guang.Application of GMT on the plotting of GPS velocity field[J].North China Earthquake Science，2010，28（3）：61-64.

[22]武永彩，保長(zhǎng)燕，唐紅濤.基于GPS速度場(chǎng)的青藏塊體東北緣應(yīng)變場(chǎng)分析[J].大地測(cè)量與地球動(dòng)力學(xué)，2016，36（11）：981-984.

WU Yong-cai，BAO Chang-yan，TANG Hong-tao.Strain field analysis of northeast margin of Qinghai-Tibet block based on GPS velocity fields[J].Journal of Geodesy and Geodynamics，2016，36（11）：981-984.

[23]Liang S，Gan W，Shen C，et al.Three-dimensional velocity field of present-day crustal motion of the Tibetan Plateau derived from GPS measurements[J].Journal of Geophysical Research Solid Earth，2013，118（10）：5722-5732.

[24]Xu C.Reconstruction of gappy GPS coordinate time series using empirical orthogonal functions[J].Journal of Geophysical Research Solid Earth，2016，121（12）：9020-9033.