GAMIT/GLOBK在高精度长基线像控点解算中的应用

杨德芳 丁开华 许伟

摘 要：青海省藏区大比例尺基础地理信息都兰测区数据采集项目，航摄结束后摄区像控点布设南北基线跨度160 km，东西基线跨度135 km，基线跨度很大，如果采用传统RTK方法，受基线距离的限制，测区必须重新建设多个基准站，外业工作效率将大幅度降低，并且常规坐标解算方法很难保障解算精度。为了解决长基线像控点高精度解算问题，利用GAMIT/GLOBK解算方式，引入了9个IGS基准站、精密星历及精密钟差对测区106个静态观测的像控点进行坐标解算，用另外一组以摄区高精度CORS（连续运行基准站）站作为控制点，区域网平差后cm级的解算成果检验GAMIT/GLOBK的解算精度。通过分析发现，采用GAMIT/GLOBK解算模式，在不依赖摄区CORS站及已知控制点的情况下，其基线解算所有单天的NRMS标准化均方根值精度值均小于0.25，基线重复率相对重复性精度达到了10-3 m、绝对重复性精度达到了10-8 m，通过TBC成果对GLOBK坐标平差精度进行评定，得出N，E，U等3方向的中误差分别为：11.9，10.4，4.2 cm，通过利用距测区150 km范围内的3个已知CORS站点对上述N，E方向中的系统性共模误差进行改正，改正后中误差为-0.728 4，-1.995 cm.解算精度满足1∶1万外业像控点要求，说明利用GAMIT/GLOBK解算方式，可以为青海省大部分无基准、无控制的地区开展GPS像控点数据测量及解算工作提供实践精度参考。

關键词：长基线像控点解算;GAMIT/GLOBK;CORS区域网平差;精度评定

中图分类号：P 228.4           文献标志码：A

文章编号：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组合平差解算方式，通过利用全球分布的IGS台站观测文件、GPS卫星精密星历、精密钟差、对流层、电离层等各类改正模型可以获取测站高精度的三维坐标信息。其解算过程中采用IGS站作为已知控制站点，因此待解算测站点可以不考虑基线距离的限制、区域网站点间共模误差的影响，可以直接实现待解算测站精密单点解算，大大解决了常规GPS测量依赖公共站点、短基线解算的缺点。

目前对于GAMIT/GLOBK解算方式的研究，纪海源通过GAMIT/GLOBK软件对汾渭断陷带及邻近地区2011—2013年的36个连续站和227个区域站GPS观测数据进行处理[1]，其2次解算的GPS连续站水平方向平均点位误差在2 mm左右，垂直分量平均点位误差在5 mm左右，一般点位误差的限差为5～10 mm.林超才利用GAMIT/GLOBK对四川地震监测网CORS站及国内部分IGS站2008年5月12日地震前后一段时间的观测数据按照地震前及地震后分别进行单日的基线解算[2]，并对基线解算结果进行精度评定和原因分析。刘宗强等对进行陆态网解算时只需要加入北半球的IGS站就可使基线得到较好的U方向解算精度，陆态网选择8个左右的IGS站较为合适[3]。高旺等实验结果表明对于中长基线，使用IGR或IGU星历进行GAMIT基线解算与使用IGS最终星历解算的坐标分量较差优于10-9数量级[4]。杨登科等从L1_ONLY和L1L2_INDEPENDENT仅适用5 km以下基线向量解算，LC_HELP和LC_AUTCLN的解算效果相当，LC_AUTCLN相对较好[5]。罗全等通过GAMIT软件解算不同截止高度角条件下的基线结果，得出当截止高度角取10°～15°时，基线解算结果最佳[6]。

由于上述优势，摄区像控点布设及施测南北基线跨度160 km，东西基线跨度135 km，点与点间基线跨度大，采用GAMIT/GLOBK组合平差模式长基线解算時必须考虑GPS长基线解算时对流层、电离层及地球曲率等因素的影响，通常情况下可以考虑2种解算模式。

第1种：测区待测点与覆盖测区范围CORS基准站联测，采用商用软件进行平差解算;

第2种：采用测区周边均匀分布的全IGS站作为基准站，利用GAMIT/GLOBK进行长基线平差解算。

第1种模式必须依赖高精度CORS站作为控制点，如果CORS站无法覆盖整个测区或是测区内没有有效的CORS站，则需联测周边已知控制点重新布设控制网进行测量，这样就增加了外业工作量及工作效率，第2种模式不受测区CORS站及已知控制点的制约，灵活采用全球分布的IGS站作为控制点，人员到达测区可直接开始静态施测，无需考虑点间基线长短，点、边连接等测量方法，就可以实现cm级的精度，文中测区站点及IGS测站分布情况及测站日观测情况分别如图1，2所示。

1 GAMIT/GLOBK解算策略

1.1 GAMIT/GLOBK解算意义

由于青海省独特的高原地理环境，目前建成的CORS站还不能覆盖该省全境，如图3所示，红色圆形区域代表每一个CORS站所覆盖的范围，若要在CORS无覆盖的测区范围内开展GPS测量工作，则存在重新布设基准控制点及联测问题[7]，此过程涉及大量人力物力支出问题，如果采用GAMIT/GLOBK解算方式，就可以解决测区无起算控制点的限制。

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

GAMIT/GLOBK在数据解算之前，首先要根据测站的观测时间对“tables表文件”进行配置，表文件主要包括系统自带的表文件和待更新的表文件2类[8-9]。由于文中观测数据是2017年，在处理之前分别对系统自带的表文件，系统定期更新的表文件及GAMIT/GLOBK控制表文件进行了更新，见表1，表2，表3.

1.2.1 系统自带的主要表文件

系统自带的主要表文件更新内容包括：天线高及相位中心偏移模式参数，大地水准面参数，海潮文件参数，周跳自动探测和修改命令参数，P1-C1，P1-P2码差等参数的更新。具体参数见表1.

1.2.2 GAMIT/GLOBK平差解算

分别采用分步处理与批处理方式进行解算，具体解算流程如图4所示。

3 解算参数配置

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

2）sestbl.中修改基线类型为LC_HELP，筛检因子Decimation Factor=1，Interval zen=2，卫星高度角改为15°，Use met.grid=N，Use otl.grid=N，Use atml.grid=N，Tides applied = 23;

3）station.info手动更新待解算点测站信息，天线高改化到天线相位中心LEPHC;

4）sittbl.中修改控制点约束平面0.02，高程005，待求点约束平面0.05，高程0.05;

5）lfile.中对于IGS站坐标从itrf08.apr中提取，对于待解算站坐标利用双差定位法提取;

6）GLOBK_comb.cmd中禁用eq_file，增加make_svs@.svs Z，改变prt_opt GDLF CMDS，改变org_opt PSUM CMDS GDLF，启用所有控制点参与平差use_site all;

7）glorg_comb.cmd修改stab_it网平差迭代次数为5次，禁用source.

4 解算精度评定

GAMIT/GLOBK基线解算及平差结束后生成结果[15-16]主要包括：H-file基线松弛解（h0526a.17146），O-file约束解（o0526a.146），Q-file过程记录文件（q0526a.146），GLOBK坐标平差文件（GLOBK_0524.org），针对上述文件进行精度分析。

4.1 标准化均方根误差NRMS精度统计

解算完后，以单天解给出了18个年积日的NRMS精度值，从表4可以看出，最大的为0.193 29E+00，最小的为0.180 08E+00，根据国内外025的数据处理经验[17]，所有单天的NRMS精度值均小于0.25，说明单时段解算出的基线值偏离其加权平均值的程度较小，处理过程中周跳得到了较好地修复，基线估算精度较高，解算质量指标符合要求。

4.2 天空图统计

LC phase Residual[18]表示接收机接受信号的好坏，由于测站较多，以测站JB70为例给出解释，如图5所示，该测站信号整体趋势在-20～+20 mm之间围绕蓝色收敛线波动，收敛趋势平顺，并且逐渐接近于0，测站RMS值为9.1 mm，说明解算过程中天线相位模型引用正确，信号没有受到较强的电磁干扰，收敛情况较好。

4.3 基线重复率

以单日解解算出了18个年积日的基线重复率，从表5，图6可以看出，绝对重复性最大值为0006 814 51 m，最小的为0.000 309 02 m，精度在mm级，相对重复性最大值为1.000E-08，最小的为-4.736E-10，相对精度达到了10-8，可见，基线解算结果精度较高。

4.4 GAMIT/GLOBK精度检验

对测区106个像控点采用GAMIT/GLOBK进行解算合格后的结果，直接采用魏子卿研究员[20]关于“中国大陆速度场”一文中给出的中国大陆3°×3°格网平均速度。将ITRF08框架下得到的空间直角坐标，换算到CGCS 2000坐标系。转换过程中分别利用公式（8）、公式（9）将2017.* * *历元转换到2 000.000历元，将ITRF 2008框架坐标转换到ITRF 97框架坐标，这样就得到了CGCS 2000坐標[19]。

为了验证GAMIT/GLOBK解算的最终点位三维坐标，用另外一组以摄区高精度CORS站作为控制点，平差后cm级的解算成果对文中的解算精度进行检验。由上图7可知，106个测站N方向残差值在0.12 m附近波动，E方向残差值在0.10 m附近波动，U方向残差值在0.04 m附近波动，经公式（6）计算得出RMS N=0.119 m，RMS E=0.104 m，RMS U=0.042 m，U方向解算精度达到了cm级，平面方向解算精度在dm级。

4.5 CORS站偏差改正

4.4中平面方向出现dm级的解算精度初步估计，是由于受测区速度场模型精度影响。图8应用GMT[21]软件，结合青藏高原东北缘水平运动GPS速度场[22-23]，给出了测区基于欧亚板块的青海省GPS速度场，图中红色矩形框为文中测区范围，从图可以看出，在测区范围内只有一个已知速度变化的测站，缺少起算点数据，无法应用局域欧拉矢量法内插速度，所以直接引用魏子卿速度场模型。测区所有站点都采用统一的速度场模型势必会引入系统性的共模误差[24]的影响，即在坐标分量上表现出整体性偏移，为了消除这一误差的影响，采用距离测区150 km范围内的3个省级CORS站进行共模误差的提取及分析。

根据图3结合测区情况，选取DGLE，YANH，GERM等3个CORS测站，利用同一Gamit/GLOBK参数解算模型，对CORS站点进行解算，将解算转换后的最终结果与CORS已知测站成果进行对比分析，得出CORS站残差为RMS N=0.126 m，RMS E=0.124 m（上述解算中高程精度已满足要求，不再进行改正分析），故将这一值，确定为本测区共模误差改正值，将这一改正值改正到测区所有像控点坐标值中。

从图7也可以看出，所有像控点N，E方向表现出整体性偏移，而且该偏移值在0.119 m与0104 m之间浮动，浮动范围±0.002 m，如果将CORS站残差值改正至上述测站整体偏移中，那么所有测站N，E方向精度将大幅度提高。图9是根据CORS站残差改正后的测站残差，从图10可以看出，经过改正后，所有测站N，E方向残差值均有所减小，N方向残差在-0.5 cm至2 cm之间浮动，E方向残差在-3 cm至-0.6 cm之间浮动，经过改正后的测站RMS N=-0.7284 cm，RMS E=-1995 cm，RMS U=0.042 m，说明经过CORS站残差改正后水平方向精度达到了亚厘米级，对于1∶1万的像控精度，上述解算精度已满足要求。

5 结 论

1）对解算策略及解算精度指标进行的分析发现：GAMIT/GLOBK基线解算所有单天的NRMS精度值均小于0.25，基线重复率相对重复性精度达到了mm级、绝对重复性精度达到了10-8;通过TBC成果对GLOBK坐标平差精度进行评定得出，解算结果RMS N=0.119 m，RMS E=0.104 m，RMS U=0042 m，N，E方向所有点位均表现出整体性偏移，并且该偏移值均值分别为11.9，10.4 cm.

2）为了解决这一偏移的影响，通过采用距离测区150 km范围内的3个省级CORS站进行该整体性偏移量的提取，得出CORS站整体性偏移改正量残差为RMS N=0.126 m，RMS E=0.124 m，将这一改正量改正到测区所有像控点坐标值中。经过改正后的所有测站RMS N=-0.728 4 cm，RMS E=-1.995 cm，RMS U=0.042 m，说明经过CORS站残差改正后水平方向精度达到了厘米级，高程精度达到了厘米级。解算结果满足文中1∶1万外业像控点要求，为青海大部门无基准、无控制的地区进行GPS数据解算工作的可行性提供了实践参考，达到了在大量减少外业工作量的同时控制成本的目的。

参考文献（References）：

[1] 纪海源.基于GAMIT/GLOBK与QOCA的汾渭断陷带地壳运動研究[D].西安：西安科技大学，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的汶川震区地壳形变分析与研究[D].成都：西南交通大学，2012.

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

[3] 刘宗强，党亚民，何 涛，等.IGS站的选取对陆态网解算精度的影响[J].测绘工程，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] 高 旺，高成发，潘树国，等.基于快速星历的GAMIT高精度基线解算研究[J].测绘科学，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].测绘地理信息，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]罗 权.不同截止高度角对GAMIT基线解算的影响分析[J].测绘地理信息，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在青海海东地区CORS基准站组网中的应用[J].测绘通报，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]徐 杰，孟 黎，唐诗华.GAMIT/GLOBK批处理在高精度海量数据处理中的应用[J].测绘科学，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，Fang 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]占 伟，武艳强，刘志广.GAMIT/GLOBK数据自处理程序的应用[J].测绘科学，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基线质量评价[J].测绘技术装备，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]张 祥.区域CORS站稳定性监测分析[D].西安：西安科技大学，2016.

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

[13]楊登科，安向东，汤 勰.GAMIT数据处理中基线解算模式的对比分析[J].测绘地理信息，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]于 强，易长荣，占 惠.ITRF2000转换到CGCS2000框架的分析[J].全球定位系统，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的中国地壳水平运动速度场模型研究[J].测绘学报，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中国大地坐标系——中国大陆速度场[J].测绘学报，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速度场的应用[J].华北地震科学，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]武永彩，保长燕，唐红涛.基于GPS速度场的青藏块体东北缘应变场分析[J].大地测量与地球动力学，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.