《(完整版)GAMITGLOBK软件操作.pdf》由会员分享,可在线阅读,更多相关《(完整版)GAMITGLOBK软件操作.pdf(5页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、GAMIT/GLOBK 数据处理报告键入文档副标题GAMIT/GLOBK 是一套高精度数据处理软件, 主要用于分析研究地壳变形、高精度GPS测量数据处理等领域。 它由美国麻省理工学院( MIT) 和斯克里普斯海洋研究所(SIO) 联合开发, 并得到美国哈佛大学和美国国家科学基金会的支持,是目前世界上应用最为广泛的高精度 GPS 数据处理软件之一。GAMIT/GLOBK 基于 UNIX(Linux)系统开发和运行。本文中所有数据处理工作均是基于Ubuntu9.0 与 csh SHELL 环境下完成的。一一 数据预备数据预备为了学习使用 GAMIT 处理 GPS 数据,本文选择 2012.07.0
2、1(DOY 183)天如下共 15个全球 IGS 跟踪站建立全球观测网。 由于该网最初是用于评定北斗电离层模型的改正精度,故而在选站上更加偏重中国及周边地区。在完成跟踪网选择后需要下载相应的导航电文和精密星历数据。这些数据亦可以通过GAMIT 中的 sh_get_rinex, sh_get_navs, sh_get_orbits脚本根据指定的站点名称和时间直接从CDDIS, SOPAC等服务器上下载。 这里需要注意的是, 由于这些脚本均是基于csh (或 tcsh)解释器,故而在 bash 环境中无法正确执行。总结准备数据的相关信息如下:IGS 站ALIC BJFS GMSD HRAO INE
3、G KELY KUNM LHAZ MOBS NRIL PENC SHAOUSUD YAKT YELL (其中 GMSD 由于观测值类型多余 15 种而被禁用)时间段星历2012 年 7 月 1 日(DOY 183)导航电文,IGS 后处理精密星历二二 建立工程建立工程根据 GAMIT 软件处理要求,需要建立相关目录。一个GAMIT 工程主要包括如下几个工程目录:DOY: processing data, final solutions, etc.rinex: observation file in RINEX o format.igs: precise orbit file from IGS
4、in sp3 orbitProject NameProject Namebrdc: broadcast file in RINEX N formattables: table files linked to / gg/tablesOther directories created during processing最初建立工程只需要在主工程目录下建立相应的DOY, igs, igs, brdc 四个目录,并在相应的目录存放数据。DOY 目录用于存放单天基线解算的中间数据文件和结果文件, 以 DOY命名;rinex 目录则是存放工程所需要的所有 rinex 观测数据文件;igs 目录存放工程所需
5、要的所有 sp3 精密星历文件;tables 目录用于存放表文件,从/gg/tables 目录中链接而来,具体的做法将在下文中给出。 其他的一些目录如 gfiles, archive, control, gsoln等是在 GAMIT 解算过程中生成的。三三 tablestables 更新与链接更新与链接表文件(tables)主要包括地球自转参数 EOP(含极移、ut1) ,跳秒,章动表,太阳和月亮星历表,卫星和接收机天线类型信息, 硬件延迟 dcb 数据等。这些表文件是软件在运行中所必需的文件, 且有部分表文件需要随时进行更新, 所有表文件均可以在以下服务器下载。ftp:/garner.ucs
6、d.edu/archive/garner/gamit/tables/每日更新:地球自转参数utl.usno、pole.usno;每年更新:章动表 nutabl.IAU00、太阳表 soltab.2012、月亮表 luntab.2012 以及跳秒表 leap.sec;可选更新:接收机天线信息rcvant.dat、卫星天线信息 svnav.dat。将上述更新后的文件复制在/gg/tables/目录下,并在工程目录下运行 sh_setup 脚本将/gg/tables 目录链接到主工程目录。链接完后需要进行检查是否出现死链接,如有出现则需要手动链接。四四 预备文件制作预备文件制作上述步骤完成之后,需要
7、制作相应的准备文件用于运行GAMIT。和谐准备文件主要包括: lfiles.、 station.info、 process.defaults、 sestbl.、 sittbl.、 sites.defaults。 其中, process.defaults、sestbl.、sittbl.、sites.defaults 用于定义处理策略,只需要在 tables 目录下对相应文件修改即可,lfiles.、station.info 是所处理测站的相关信息,需要根据rinex 观测文件生成。测站信息文件测站信息文件 station.infostation.info在 rinex 目录下,复制 tables
8、 目录下 station.info 文件并将其中测站信息删掉只留下文件头部分, 执行 sh_upd_stnfo 脚本即可生成新的 station.info 文件。 将新的 station.info 文件复制到 tables 目录下并覆盖原文件。近似坐标文件近似坐标文件 lfiles.lfiles.在 rinex 目录下,执行 grep POSITION *.12olfiles.rnx 命令将 o 文件中的近似坐标值获取到 lfiles.rnx 文件中,再采用命令 rx2apr lfiles.rnx 2012 006 可以得到 apr 格式的文件lfiles.rnx.apr, 最后采用命令 ga
9、pr_to_l lfiles.rnx.apr lfiles. 2012 006, 将 apr 文件转换为 lfiles.。同样需要将 lfiles.复制到 tables 目录下。在制作好上述两个文件后,一般情形下,目录控制文件process.defaults、测站列表文件sites.defaults、解算控制文件 sestbl.、精度控制文件 sittbl.均可以保持默认值。如数据处理中不考虑海洋潮汐的影响, 可以禁用海洋潮汐模型otl.grid 文件, 需要在 sestbl.文件中设置为:Tides applied = 23;Use otl.grid = N。本文在处理中禁用了海洋潮汐改正。
10、五五 数据处理及结果分析数据处理及结果分析GAMIT 软件数据处理整体流程比较复杂。本文直接采用了 sh_gamit 脚本进行批处理。具体命令:sh_gamit -s 2012 183 183-expt chlg-eops usno nogifs y执行该命令后可以进行自动的数据处理。基线处理的结果在 DOY 目录下,在本文中即是183 目录下。为了评价数据处理结果的精度好坏,GAMIT 软件在 Q 文件中给出了精度指标nrms 值(q*.a.DOY 给出的是模糊度固定后解算结果,q*.p.DOY 给出的是模糊度是浮点解时的解算结果) 。查看本文qchlla.183 最终处理的精度指标值 nr
11、ms 如下:FROMFROMqchlla.183qchlla.183FILE:FILE:- * Summary of biases-free solution *- Total parameters: 1088 live parameters: 658 Prefit nrms: 0.13216E+01 Postfit nrms: 0.18730E+00 - Uncertainties not scaled by nrmsEnd of loose solution with LC observable and ambiguities free-Performing LC biases-fixed loose solution- * Summary of biases-fixed solution *- Total parameters: 1088 live parameters: 343 Prefit nrms: 0.13174E+01 Postfit nrms: 0.19374E+00 - Uncertainties not scaled by nrmsEnd of loose solution with LC observable and ambiguities fixd-Normal stop in SOLVE
限制150内