应用GAMIT对加入长基线网的GPS数据解算精度分析

密级：页数：xxxxxx毕业设计（论文）题目：应用GAMIT对加入长基线网的GPS 数据解算精度分析学员姓名xxx学号xxxxxxxx所在单位xxxx指导教师xxx完成日期2012.6.15摘要GPS技术已在多学科各个领域有着很重要的应用，GPS相对定位精度己经达到10-9。为了得到GPS数据较好的解算，GPS数据处理越来越受到重视。本文从数据处理的整个流程出发，探讨GAMIT软件的使用，并应用GAMIT软件处理中国区7个IGS站、3个网络基准站和3个新测点的三天的数据，研究GPS数据在没有已知点的情况下加入IGS基准站构成长基线网前后解算精度的变化，说明加入IGS基准站构成长基线网后对未知测站坐标解算的精度有所提高。索引关键词:GAMIT，GPS数据处理，长基线网，精度ABSTRACT GPS technology has been used in lots of fields，and GPS relative positioning accuracy has been achieved the level of 10-9. In order to get the GPS data better solution to calculate, GPS data processing more and more attention to . This article from the starting from the data processing of the entire process to explore the use of GAMIT software, and application GAMIT software to handle the seven IGS stations, three network base station and three new points-day data, research and GPS data does not have add IGS reference station to join the case that a point constitute a long baseline solver before and after changes in the accuracy of the IGS reference stations constitute a long baseline unknown station coordinates solver accuracy improved.Keywords: GAMIT,GPS data processing, Long baseline nets, accuracy目录目录摘要ABSTRACT第一章 绪论11.1选题背景11.2 GAMIT/GLOBK软件介绍1第二章 GAMIT软件的安装32.1 GAMIT/GLOBK软件安装步骤32.2绘图工具的安装4第三章 数据解算63.1观测数据文件和星历文件预准备63.2参数表文件预准备63.3数据准备73.3数据处理7第四章 实验结果及分析114.1处理方案114.2解算结果及精度评定12第五章 总结20后记21参考文献22III第一章 绪论1.1 选题背景随着GPS技术的高速发展，地面点的坐标精度，由最初的分米量级，提高到目前的厘米甚至毫米级，同时大地测量、工程测量以及地球动力学等领域对GPS数据的精度要求越来越高，随机附带的基线解算软件，由于种种原因，已越来越不能满足这么高精度的要求。在此环境下，国内外研制出了很多优秀的数据处理软件,其中包括瑞士伯尔尼大学天文研究所研制的BERNESE软件, 美国宇航局(NASA)喷气推进实验室( JPL )研制的GIPSY /OASIS 软件,以及本文将要用到的GAMIT/GLOBK软件。在精度方面，这些软件，都可得到厘米级的点位坐标精度。但是由于在实际测量工作中可能没有精确的已知点来解算未知测站坐标，在此情况下如何解算未知测站坐标，如何提高解算未知测站坐标的精度，值得探讨。众所周知全球有很多IGS连续运行基准站，这些基准站具有较高的精度和良好的可靠性。因此，使用GAMIT/GLOBK软件研究GPS数据在没有已知点的情况下加入IGS基准站构成长基线网对未知测站坐标解算精度有多大提高，这便于我们进一步提高在没有已知测站信息情况下数据解算精度。1.2 GAMIT/GLOBK软件介绍GAMIT/GLOBK软件是美国麻省理工学院（MIT）和美国加利福尼亚大学斯克里普斯海洋研究所（SIO）合作开发,并得到美国哈佛大学和美国国家科学基金会的支持。它是一套数据分析处理软件包，用于GPS高精度数据的后处理和分析, 用来分析和估计陆地测站及人造卫星轨道的三维相关位置。一、 GAMIT介绍 GAMIT 软件包主要由以下几个程序构成: ARC ( 轨道积分程序) 、MODEL (观测值模型程序) 、SINCLN (单差自动修复周跳程序) 、DBCLN (双差自动修复周跳程序) 、CVIEW(人工交互式修复周跳程序) 、SOLVE(最小二乘解算程序) 、FXDRV (生成数据处理的程序) 、GLOBK(运用卡尔曼滤波进行网平差模块) 等。这些模块即可以单独运行，也可以用批处理命令联在一起运行，最大限度地减少人为操作，提高运算效率。软件的执行程序放在/com、/kf/bin和/gamit/bin三个目录下。GAMIT软件在解算时考虑了潮汐改正、大气负荷改正、岁差章动改正、极移改正等多项改正并将所有参数一起平差，可以解算卫星轨道、测站坐标、大气延迟、整周模糊度等可得到更可靠的基线数据。同时GAMIT 软件是开放源代码,可以安装在免费开放源代码的LINUX 操作系统上,用户可以根据需要对源程序修改,便于科研工作,因此在国内外广泛应用于长距离、高精度、长时间的GPS数据的处理。二、 GLOBK介绍GLOBK是一个卡尔曼虑波器，主要的目的是综合空间大地测量和经典大地测量的初步处理结果完成数据的后处理。因此它的输入一般是一些准观测量如测站坐标、地球自转参数、卫星轨道及它们的方差-协方差。实际上GLOBK是一个基线解的网平差软件，它将GAMIT的H-文件作为输入文件，采用卡尔曼滤波算法估计测站坐标、速度，卫星轨道参数以及EOP参数。GLOBK三个主要应用是：结合一个观测作业期内不同时段的初步处理结果，获取该观测作业期的测站坐标最佳估值。对GPS分析，轨道参数可做随即参数处理；结合不同年份获取的测站坐标结果估计测站的速度；将测站坐标作为随机参数，生成每个时段或每个观测作业期的坐标结果以评估观测质量。GLOBK软件的主要功能模块包括htoglb、glred、globk、glorg，各自的功能如下： htoglb：将GAMIT的H-文件转为GLOBK认可的二进制的h文件。glred：计算单天解的重复性。globk：估计测站坐标和速度。glorg：加约束的联合解。第二章 GAMIT软件的安装2.1 GAMIT/GLOBK软件安装步骤：一、在用户中建一个目录（我建的为/home/soft/gamit）, 并将source 文件夹(用于存放GAM IT安装文件)中的文件复制到该目录中。install_software 安装批处理文件install_updates 升级脚本com.10.4.tar.gz 组件压缩包tables.10.4.tar.gz 表格文件包gamit.10.4.tar.gz gamit软件压缩包help.10.4.tar.gz 帮助系统压缩包kf.10.4.tar.gz kf软件压缩包libraries.10.4.tar.gz 库文件压缩包maps.10.1.tar.gz 地图数据压缩包incremental_updates 增量升级包etopo5.grd.10.0.tar.gz 全球地形数据包example.10.4.tar.gz 示例包relnote.10.4 版本说明其中前八个为必选软件包，后五个是可选软件包, 但建议一并拷入, 完全安装。GAMIT/GLOBK软件的源代码可以从：/simon/gtgk/down.htm，上申请账户下载。二、正式安装，在当前目录下运行install-software（命令： ./ install-software），安装过程中根据提示给以回应即可顺利完成。三、设置工作路径。安装完后并不能直接工作，需要给linux指引程序运行的路径。由于本人的当前工作的为Tc shell,需要在系统中修改.tcshrc文件,并且要符合Tc shell语法要求:set path = (./opt/GMT/bin $HOME/gg/gamit/bin $HOME/gg/kf/bin $HOME/gg/com $path)setenv GMTHOME /opt/GMTsetenv HELP_DIR $HOME/gg/help/setenv INSTITUTE “chxy”修改完后就可以进行运行GAMIT软件了。2.2 绘图工具的安装GAMIT/GLOBK软件在处理数据时可能要绘图，这就要安装绘图工具netcdf。一、在安装文件夹中拷入以下文件：GMT4.4.0_doc.tar.bz2GMT4.4.0_share.tar.bz2GMT4.4.0_src.tar.bz2GMT4.4.0_suppl.tar.bz2GMT4.4.0_triangle.tar.bz2GSHHS_high.tar.bz2GSHHS1.10_coast.tar.bz2GSHHS1.10_full.tar.bz2GSHHS1.10_high.tar.bz2netcdf.tar.gz二、进入安装文件夹，执行以下命令：tar xvfz netcdf.tar.gzfor i in *.bz2do tar jxvf $idone这时候生成如下的文件夹netcdf4.0GMT4.4.0Share三、首先安装netcdf，执行以下命令：cd netcdf4.0./configure -prefix=/opt/netcdfmakesudo make install四、安装GMT，执行以下命令：cd GMT4.4.0./configure -prefix=/opt/GMT -enable-netcdf=/opt/netcdfmakesudo make install五、进入暗转文件夹，拷贝数据文件sudo cp -r GMT4.4.0/share/ /opt/GMT/sudo cp -r share/coast /opt/GMT/share这样GAMIT/GLOBK就可以正常绘图了。第三章 数据解算 3.1 观测数据文件和星历文件预准备GAMIT采用标准的RINEX格式观测数据O文件以及导航N文件，应利用GPS 接收机的随机软件将GPS 原始观测数据转换成RINEX 标准数据格式,以供GAMIT软件计算使用，下载并准备好所需的星历文件。N文件和IGS连续站O文件可以在3个IGS全球数据库下载：（1）/pub/rinex/YYYY/DOY/ N文件是以auto开头的文件，下载后改为brdc开头，O文件是下载d文件，然后用crx2rnx解压；（2）/pub/gps/data/YYYY/brdc/ 下载 N文件，以brdc开头的文件， /pub/gps/data/YYYY/DOY/YYo/ 下载O文件；（3）ftp:/igs.ensg.ign.fr/pub/igs/data/YYYY/DOY/ N文件是以brdc开头的文件，O文件是下载d文件，然后用crx2rnx解压精密星历（sp3文件）下载地址：/pub/products/WWWW//pub/gps/products/WWWW/ftp:/igs.ensg.ign.fr/pub/igs/products/WWWW/ 3.2 参数表文件预准备各种参数表可以通过IGS网站下载，也可通过下载GAMIT软件的网站下载表文件(tables)包括有：测段对应周的参数文件,包括:pole.（极移参数）、ut1. UT1（时间表）；测段对应年的参数文件, 包括:luntab. （月球星历表）、soltab.（太阳星历表）、leap.sec（跳秒表）、nutab. （章动历表）； 其他最新文件,包括: gdetic.dat（大地水准面参数表）、antmod.dat (天线高及相位中心偏移模式参数表)、svnav.dat (卫星数目编号信息)、rcvant.dat （接收机及天线类型信息）、autcln.cmd（模块执行命令数据）、hi.dat (天线型号及量高方式改正数参数)、svs_exclude.dat (卫星的质量情况数据)、dbc.dat (频率相位中心参数改正文件)等文件。下载地址：/pub/gamit/tables/ /pub/gamit/setup/ 3.3 数据准备一、在目录下建立一个工程目录,用于存放本次处理的所有数据,以下建立的所有目录均在此目录下。二、建立子目录rinex,用以存放处理好的RINEX格式的数据的观测值文件O文件。三、建立星历目录igs,用以存放从网上下载对应时段的精密星历sp3文件。四、建立子目录brdc,用来存放所有的N文件。五、建立tables目录,用于存放数据处理过程中所需的各种表格。此目录下的文件可以用sh_setup命令从GAMIT自带的文件中连接过来，如果有必要也可以把事先从网上下载的拷贝到此目录中。此外测站的概略坐标文件（vg_in）也存放在目录下，测站概略坐标文件存放着测站的先检验坐标和精度，一般情况下先检验坐标误差应该小于10m,这个可以由单点定位或者导航解得到。文件格式为：四字符测站名及四字符注释（_GPS）、测站XYZ坐标（米为单位）、测站XYZ速度（米/年为单位）、上述坐标历元、测站XYZ约束六、准备和文件文件存放着测站天线高、接收机代码、天线代码、天线高量测方式及观测时间范围等。文件可以使用sh_upd_stnfo命令添加，如果需要手工添加时，需要严格根据模板进行修改,不能出现多余的空格和字符。文件格式为：四字符测站名、四字符测站别名、测站全名、天线高（大地坐标系中HNE）、六字符接收机代码、六字符天线代码、六字符天线高量测方式代码、接收机软件版本号、观测年及年积日、观测时段起止时分秒及注释文字。文件也可以根据模板进行修改, 或使用makexp 自动创建。它包含测站先验约束、对流层模型、对流层解算时段、计算截至高度角及拟合钟差多项式阶数等十多项。七、其他文件的准备其他文件例如J文件、K文件、C文件以及非常重要的X文件的准备,可以通过运行gamit程序自动创建。以上工作完成后就可以进行数据解算和处理分析了。 3.4 数据处理准备好文件之后,就开始进入工程目录进行处理,具体操作如下:一、基线解算(1) 在工作目录下运行sh_setup命令，把表文件连接起来，命令： sh_setup yr 2010（2）用cd命令进入tables目录下，命令： cd tables 在tables目录下，用vi命令编辑sites.defaults文件，把该文件里没有的测站添加进去，并保存退出，命令： vi sites.defaults在当前目录下，用vi命令编辑process.defaults文件,把文件中的itrf05.apr修改成自己所需的坐标框架，命令： vi process.defaults在tables目录下，修改文件。运行sh_upd_stnfo命令，把文件中所有的包含sites.defaults文件中记录的测站的信息提取出来，会生成一个新的文件.new，命令： sh_upd_stnfo l sd然后把.new文件名修改为，原来的删除掉，可以用mv命令实现，命令：mv .new 最后要把中的没有的测站，用sh_upd_stnfo命令把未知测站的信息通过rinex目录下的O文件，写进文件中，命令：sh_upd_stnfo files ./rinex/*.00o（3）在tables目录下，还可以编辑sestbl.文件，根据需要对sestbl.文件进行编辑修改。sestbl.文件编辑修改时应该主要关注文件中以下内容：站及卫星约束控制利用这些设置可约束指定卫星或测站，其中：Station Number = * ，代表所有测站需要先验约束；Station Constraint = Y ，表示测站约束可用；Satellite Number = * ，表示全部卫星；Satellite Constraint = Y ，表示卫星约束可用。卫星可逐个约束，也可利用下面形式对所有卫星使用相同约束 all a e i node arg per M rad1 rad2 rad3 rad4 rad5 rad6 rad7 rad8 rad9 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01上述约束值对六个轨道根数单位为110-6，而对九个光压参数单位为初值的百分比。数据解算模式的设定可以利用设置“Type of Analysis”和“Choice of Experiment”参数值，指定解算模式。其中“Type of Analysis”参数有六个备选值，分别表示参数解算叠代次数及自动剔除周跳叠代次数。而“Choice of Experiment”共有三个备选值分别表示求解形式即Baseline（仅求基线解）、Relax（同时解算轨道及基线）和Orbit（仅解算轨道）。涉及数据量选用的参数本类参数确定数据处理时可选的数据量，其中：参数“Choice of Observable”共有六个备选值，分别表示数据处理时采用的观测值为单频、双频、无电离层观测量及是否使用伪距观测量。参数“Elevation cutoff”指定数据处理时选用卫星的截至高度角。参数“Decimation Factor”及“Quick-predecimation factor”分别指定解算时数据筛选因子。使用的潮汐模型制参数“Tides applied”确定数据处理时使用的潮汐改正，包含固体潮、极潮和海潮。（4）用cd命令进入工程目录，在 工程目录下执行sh_gami命令，进行基线解算，命令为：sh_gamit -expt chxy -d 2010 111 112 113 -pres ELEV -orbit IGSF -copt x k p -dopts c ao noftp注：chxy是本机注册名；-pres ELEV的功能是画天空图；因为本人的计算机不能上网下载数据，所以采用-noftp选项；如果没有-noftp，会首先下载需要的数据，以及表文件，然后才是解算。（5）使用sh_glred命令进行平差，查看每个测站的3天重复性，命令为：sh_glred -s 2010 111 2010 113 -expt scal -opt H G E &! sh_glred.log （6）使用h_cleanup命令，删除没用的文件，命令为：sh_cleanup -d 2010 111 112 113 -dopts p x k 二、 网平差通过GLOBK平差,即得在设定的框架坐标系下的坐标。解算步骤： 1）执行glred或globk之前，需将二进制的h-文件列入一个扩展名为.gdl的文件中，命令： ls ./glbf/h*.glx chxy.gdl 2）运行glred进行重复性精度评价，命令：glred std out print file log file list file command file（如：glred 6 globk_comb.prt globk_comb.log chxy.gdl globk_comb.cmd） grep POS STAT globk_如果安装有GMT绘图软件，则可以对输出结果进行图形显示，按照上面的输出文件，绘图过程为：sh_plotcrd -f globk_ -s long -res -o 1 -vert -col 1 -x 2000.0 2005.03）运行globk、glorg进行网平差，命令：globk std out print file log file list file command file（如：globk 6 globk_vel.prt globk_vel.log chxy.gdl globk_vel.cmd）如果安装有GMT绘图软件，则可以对输出结果进行图形显示，按照上面的输出文件，绘图过程为：sh_plotvel -ps chxy -f globk_ R73/135/04/54 -factor 0.5 -arrow_value 10 -page L第四章 实验结果及分析4.1 处理方案本文选取了2010年4月21日至2010年4月23日的三天的5080、5089、5093 三个新测点的实测数据（注：数据由教员提供）。三个新测点使用了三台TRIMBLE R7双频接收机/TRIMBLE GNSS CHOCK RING (59800-00C) 天线，采用同步观测3个时段，时段长度23小时40分，采样间隔30秒，截止高度角10度，进行观测。点位分布情况：未加入IGS站前的点位分布：图4-1未加入IGS站后的点位分布：图4-2处理方法：先使用3个网络工程基准站（YANC、TAIN、CHUN）作控制，应用GAMIT/GLOOBK软件进行数据处理；然后再在以上基础上添加7个连续运行IGS基准站（BJFS、CHAN、KUMN、LHAZ、SHAO、URUM、WUHN），构成长基线网，再应用GAMIT/GLOOBK软件进行数据处理。最后，对两组数据的精度进行分析。观测时间与GPS周对照表表4-1观测时间年积日GPS周星期2010.4.211111580三2010.4.221121580四2010.4.231131580五控制参数的选择如下:处理模式 RELAX.观测值的选择 LC_HELP批处理迭代方案 12ITER卫星的截止高度角 15.对流层延迟改正投影函数 GMF光压模型 BERNE数据筛选 AUTCLN惯性框架 CGCS2000测站坐标约束 10m10m10m轨道约束 100m100m 100m4.2解算结果及精度评定 经过GAMIT/GLOOBK软件进行数据处理之后,可以从生成的Q文件或O文件以及.org文件中查看结果并进行数据分析。一、 基线结果分析（1）均方根残差( nrms)比较均方根残差计算模型： （4-1）上式中为观测时段，为第时段基线s的中误差，为第时段基线解算结果，为总的时段数。未加入IGS站前的基线解算结果中nrms为:表4-2年积日nrms1110.165831120.164611130.16624加入IGS站后的基线解算结果中nrms为：表4-3年积日nrms1110.154001120.155151130.16339均方根残差( nrms)是衡量基线解算质量的最主要的指标之一。要求nrms的值原则上最好小于0.3,如果nrms值大于0.5,则说明本次处理中仍有大周跳未被清除,还需要查看自动清除周跳的总结文件autcln.sum或使用CV IEW 进一步人工剔除。通过表4-2和表4-3可看出加入IGS站构成长基线网后对数据的基线解算的均方根残差普遍减少了0.01左右，基线质量得到提高，这有助于提高未知测站坐标解算精度。（2） 不同时段较差及较差分布比较 未加入IGS站前的基线解算中时段较差及较差分布：不同时段较差表4-4基线名称X分量（mm）Y分量（mm）Z分量（mm）L(mm)基线长度5080-50894.515.041.613.231925080-50932.013.091.810.231145080-CHUN3.7613.333.6111.168805080-TAIN3.527.725.412.064245080-YANC8.923.082.178.677925089-50933.672.491.880.951165089-CHUN5.248.292.969.556895089-TAIN3.367.084.445.155375089-YANC6.796.870.948.929715093-CHUN1.7610.242.568.177815093-TAIN4.414.986.334.214235093-YANC0.984.932.8311.26904CHUN-TAIN5.989.865.5712.481098CHUN-YANC11.3715.173.3518.371651TAIN-YANC5.396.473.517.08878段基线总数：15较差分布情况表4-5分段分量0,2mm）（个）2,4mm）（个）4,6mm）（个）6,8mm）（个）（8mm）（个）最大值（mm）最小值（mm）X2551211.370.98Y0334515.172.49Z473106.330.67L2231818.370.23加入IGS站后的基线解算中时段较差及较差分布： 不同时段较差表4-6基线名称X分量（mm）Y分量（mm）Z分量（mm）L(mm)基线长度5080-50894.143.661.941.641925080-50931.451145080-BJFS3.6612.248.942.11445080-CHAN4.769.986.952.658805080-CHUN2.0510.414.964.558805080-KUNM11.0113.096.276.6820635080-LHAZ6.4513.586.83.5825315080-SHAO6.5112.447.513.410835080-TAIN4.268.744.541.864245080-UROM13.8910.187.3415.0523805080-WUHN5.7717.8712.291.4910575080-YANC4.045.363.431.25795089-50933.051165089-BJFS1.189.9373.932335089-CHAN0.87.685.023.546895089-CHUN2.098.13.046.146895089-KUNM6.8711.736.665.1922495089-LHAZ3.2612.096.490.7527105089-SHAO2.388.775.972.5911365089-TAIN2.056.432.62.155375089-UROM14.847.876.9616.0224865089-WUHN4.4314.2210.741.612045089-YANC2.921.751.892.49715093-BJFS3.559.368.113.31425093-CHAN3.3578175093-CHUN1.357.524.153.6178165093-KUNM9.5912.847.677.6121435093-LHAZ6.36726395093-SHAO92.8610445093-TAIN4.165.853.721.924235093-UROM15.217.296.2916.224845093-WUHN5.6714.7412.871.3810885093-YANC3.942.234.012.92904BJFS-CHAN1.992.251.992.4919BJFS-CHUN2.211.843.971.57919BJFS-KUNM7.5817.8211.396.592020BJFS-LHAZ2.7918.1911.233.072499BJFS-SHAO3.17.588.385.441058BJFS-TAIN0.863.624.411.93391BJFS-UROM13.652.074.0513.742374BJFS-WUHN3.2517.3115.151.971015BJFS-YANC1.738.195.571.91762CHAN-CHUN2.710.731.991.410CHAN-KUNM6.2515.6510.567.62899CHAN-LHAZ4.0216.0110.392.473372CHAN-SHAO1.985.47.532.351454CHAN-TAIN2.843.542.423.871098CHAN-UROM15.641.423.0715.42990CHAN-WUHN5.2215.1414.310.41761CHAN-YANC3.725.934.352.281651CHUN-KUNM8.9616.389.0110.352899CHUN-LHAZ5.016.748.845.493372CHUN-SHAO4.466.135.980.831453CHUN-TAIN2.813.240.673.421098CHUN-UROM14.490.934.6214.092990CHUN-WUHN4.3215.8712.761.891761CHUN-YANC2.596.352.83.61651KUNM-LHAZ7.534.21.037.091263KUNM-SHAO4.4913.314.68.91920KUNM-TAIN7.8517.958.349.161840KUNM-UROM19.1617.0813.6213.642483KUNM-WUHN8.7496.35101289KUNM-YANC8.9411.786.212.721477LHAZ-CHAU3.7910.613.5712864LHAZ-TAIN2.1918.318.173.992515LHAZ-UROM11.6217.4413.452.451597LHAZ-WUHN93.032228LHAZ-YANC2.4112.146.032.261752CHAU-TAIN3.358.565.324.52681CHAU-UROM15.46.8310.617.973225CHAU-WUHN4.999.736.781.7657CHAU-YANC4.447.074.084.631460TAIN-UROM12.792.315.2813.692622TAIN-WUHN2.3817.4412.11.22680TAIN-YANC1.03878UROM-WUHN10.4116.5617.3814.272755UROM-YANC11.926.127.4212.831791WUHN-YANC1.7312.519.963.241024长基线总数：68条较差分布情况表4-7分段分量0,2mm）（个）2,4mm）（个）4,6mm）（个）6,8mm）（个）（8mm）（个）最大值（mm）最小值（mm）X2328661519.160.8Y696134418.310.73Z71117192417.380.67L2328661517.970.4通过对比分析表4-5和表4-7，可见加入IGS站构成长基线网后，较差的最小值变小，但最大值也变大了，最小值所占的比例相对的变大了，说明加入IGS站构成长基线网后对数据解算的基线质量有一定的提高。二、 网平差结果分析 （1）点位中误差比较 未加入IGS站前的网平差结果中点位中误差：表4-8 点号X (m)Y (m)Z (m)50800.048240.032480.0269850890.092900.040400.0376850930.065820.030490.02740 加入IGS站后的网平差结果中点位中误差：表4-9点号X (m)Y (m)Z (m)50800.023440.013120.0104650890.044680.017190.0146850930.031920.014150.01132对比表4-8和表4-9中的网平差结果的中点位中误差，明显看出点位坐标的X，Y，Z三个分量上中误差显著减小14个厘米左右，说明加入IGS站构成长基线网后对数据解算的坐标精度有明显提高。（2）基线重复率(WRMS)比较基线重复性计算模型： （4-2）上式中为观测时段，为基线s的重复性统计值，为第时段基线s的中误差，为第时段基线解算结果，为总的时段数，为基线结果的加权平均值。有 （4-3）未加入IGS站前的网平差结果中WRMS：表4-105080508950931111.03.52.5 1128.713.49.41137.718.515.3加入IGS站后的网平差结果中WRMS：表4-11508050895093 1 1 1评定基线的精度通常用基线重复率（WRMS）来衡量。基线重复率是相对几个