(完整版)GAMITGLOBK软件操作.doc

《(完整版)GAMITGLOBK软件操作.doc》由会员分享，可在线阅读，更多相关《(完整版)GAMITGLOBK软件操作.doc（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、GAMIT/GLOBK 数据处理报告键入文档副标题GAMIT/GLOBK 是一套高精度数据处理软件， 主要用于分析研究地壳变形、高精度GPS 测量数据处理等领域。它由美国麻省理工学院( MIT) 和斯克里普斯海洋研究所(SIO) 联合开发， 并得到美国哈佛大学和美国国家科学基金会的支持，是目前世界上应用最为广泛的高精度 GPS 数据处理软件之一。GAMIT/GLOBK 基于 UNIX（Linux）系统开发和运行。本文中所有数据处理工作均是基于 Ubuntu9.0 与 csh SHELL 环境下完成的。一 数据预备为了学习使用 GAMIT 处理 GPS 数据，本文选择 2012.07.01（DO

2、Y 183）天如下共 15个全球 IGS 跟踪站建立全球观测网。由于该网最初是用于评定北斗电离层模型的改正精度，故而在选站上更加偏重中国及周边地区。在完成跟踪网选择后需要下载相应的导航电文和精密星历数据。这些数据亦可以通过GAMIT 中的 sh_get_rinex, sh_get_navs, sh_get_orbits 脚本根据指定的站点名称和时间直接从CDDIS，SOPAC等服务器上下载。这里需要注意的是，由于这些脚本均是基于 csh（或 tcsh）解释器，故而在 bash 环境中无法正确执行。总结准备数据的相关信息如下：IGS 站 ALIC BJFS GMSD HRAO INEG KELY

3、 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 in sp3 orb

4、itProject 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 目录存放工程所需要的所有 sp3 精密星历文件；tables

5、目录用于存放表文件，从/gg/tables 目录中链接而来，具体的做法将在下文中给出。其他的一些目录如 gfiles, archive, control, gsoln 等是在 GAMIT 解算过程中生成的。三 tables 更新与链接表文件（tables）主要包括地球自转参数 EOP（含极移、ut1），跳秒，章动表，太阳和月亮星历表，卫星和接收机天线类型信息，硬件延迟 dcb 数据等。这些表文件是软件在运行中所必需的文件，且有部分表文件需要随时进行更新，所有表文件均可以在以下服务器下载。ftp:/garner.ucsd.edu/archive/garner/gamit/tables/ 每日更新

6、：地球自转参数 utl.usno、pole.usno； 每年更新：章动表 nutabl.IAU00、太阳表 soltab.2012、月亮表 luntab.2012 以及跳秒表 leap.sec； 可选更新：接收机天线信息 rcvant.dat、卫星天线信息 svnav.dat。将上述更新后的文件复制在/gg/tables/目录下，并在工程目录下运行 sh_setup 脚本将/gg/tables 目录链接到主工程目录。链接完后需要进行检查是否出现死链接，如有出现则需要手动链接。四 预备文件制作上述步骤完成之后，需要制作相应的准备文件用于运行GAMIT。和谐准备文件主要包括 ：lfiles.、st

7、ation.info、process.defaults、sestbl.、sittbl.、sites.defaults。其中，process.defaults、sestbl.、sittbl.、sites.defaults 用于定义处理策略，只需要在 tables 目录下对相应文件修改即可，lfiles.、station.info 是所处理测站的相关信息，需要根据 rinex 观测文件生成。 测站信息文件 station.info在 rinex 目录下，复制 tables 目录下 station.info 文件并将其中测站信息删掉只留下文件头部分，执行 sh_upd_stnfo 脚本即可生成新的

8、station.info 文件。将新的 station.info 文件复制到 tables 目录下并覆盖原文件。 近似坐标文件 lfiles.在 rinex 目录下，执行 grep POSITION *.12olfiles.rnx 命令将 o 文件中的近似坐标值获取到 lfiles.rnx 文件中，再采用命令 rx2apr lfiles.rnx 2012 006 可以得到 apr 格式的文件lfiles.rnx.apr，最后采用命令 gapr_to_l lfiles.rnx.apr lfiles. 2012 006，将 apr 文件转换为 lfiles.。同样需要将 lfiles.复制到 ta

9、bles 目录下。在制作好上述两个文件后，一般情形下，目录控制文件process.defaults、测站列表文件sites.defaults、解算控制文件 sestbl.、精度控制文件 sittbl.均可以保持默认值。如数据处理中不考虑海洋潮汐的影响，可以禁用海洋潮汐模型 otl.grid 文件，需要在 sestbl.文件中设置为：Tides applied = 2 3；Use otl.grid = N。本文在处理中禁用了海洋潮汐改正。五 数据处理及结果分析GAMIT 软件数据处理整体流程比较复杂。本文直接采用了 sh_gamit 脚本进行批处理。具体命令：sh_gamit -s 2012 1

10、83 183-expt chlg-eops usno nogifs y执行该命令后可以进行自动的数据处理。基线处理的结果在 DOY 目录下，在本文中即是183 目录下。为了评价数据处理结果的精度好坏，GAMIT 软件在 Q 文件中给出了精度指标 nrms 值（q*.a.DOY 给出的是模糊度固定后解算结果，q*.p.DOY 给出的是模糊度是浮点解时的解算结果）。查看本文 qchlla.183 最终处理的精度指标值 nrms 如下：FROM qchlla.183 FILE:-* Summary of biases-free solution *-Total parameters: 1088 li

11、ve parameters: 658Prefit 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: 343Prefit 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