完整版GAMITGLOBK软件操作

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

2、PS数据，本文选择 2012.07.01 （ DOY 183）天如下共15 个全球IGS跟踪站建立全球观测网。由于该网最初是用于评定北斗电离层模型的改正精度,故而在选站上更加偏重中国及周边地区。-12Q-600601M*.TELL60*O.Alic-6or-IMeoISO在完成跟踪网选择后需要下载相应的导航电文和精密星历数据。这些数据亦可以通过GAMIT中的sh_get_rinex, sh_get_navs, sh_get_orbits脚本根据指定的站点名称和时间直接从CDDIS , SOPAC等服务器上下载。这里需要注意的是，由于这些脚本均是基于 csh （或tcsh）解释器，故而在bash

3、环境中无法正确执行。总结准备数据的相关信息如下:IGS站ALIC BJFS GMSD HRAO INEG KEL Y KUNM LHAZ MOBS NRIL P ENC SHAOUSUD YAKT YELL （其中GMSD由于观测值类型多余 15种而被禁用）时间段2012 年 7 月 1 日（ DOY 183）星历导航电文，IGS后处理精密星历二建立工程将上述更新后的文件复制在/gg/tables/目录下，并在工程目录下运行sh_setup脚本将根据GAMIT软件处理要求，需要建立相关目录。一个GAMIT工程主要包括如下几个 工程目录:DOY: P rocess ing data, fi na

4、l soluti ons, etc.rin ex: observati on file in RINEX o format.igs: p recise orbit file from IGS in sp3 orbitP roject Namebrdc: broadcast file in RINEX N formattables: table files lin ked to / gg/tablesOther directories created duri ng p rocess ing最初建立工程只需要在主工程目录下建立相应的DOY, igs, igs, brdc四个目录，并在相以DOY应

5、的目录存放数据。DOY目录用于存放单天基线解算的中间数据文件和结果文件,命名；rinex目录则是存放工程所需要的所有rinex观测数据文件；igs目录存放工程所需要的所有sp3精密星历文件；tables目录用于存放表文件，从/gg/tables目录中链接而来，具体的做法将在下文中给出。其他的一些目录如 gfiles, archive, control, gsoln等是在GAMIT解 算过程中生成的。tables更新与链接表文件（tables）主要包括地球自转参数EOP （含极移、ut1）,跳秒，章动表，太阳和月亮星历表，卫星和接收机天线类型信息，硬件延迟dcb数据等。这些表文件是软件在运行 中

6、所必需的文件，且有部分表文件需要随时进行更新，所有表文件均可以在以下服务器下载。ftp: /gar /archive/gar ner/gamit/tables/每日更新：地球自转参数utl.usno、pole.usno;每年更新：章动表nutabl.lAUOO、太阳表soltab.2012、月亮表luntab.2012以及跳秒表 leap.sec；/gg/tables 目录链接到主工程目录。链接完后需要进行检查是否出现死链接，如有出现则需 要手动链接。四 预备文件制作上述步骤完成之后，需要制作相应的准备文件用于运行GAMIT 。和谐准备文件主要包括：Ifiles.、st

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

8、脚本即可生成新的 文件。 将新的 文件复制 到 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.。同样需要将

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

10、it -s 2012 183 183-ex pt chlg-eo ps us no - nogifs y执行该命令后可以进行自动的数据处理。基线处理的结果在 DOY 目录下，在本文中即是 183目录下。为了评价数据处理结果的精度好坏，GAMIT软件在Q文件中给出了精度指标门口5值（q*.a.DOY给出的是模糊度固。查看本文 qchlla.183最定后解算结果，q*.p.DOY给出的是模糊度是浮点解时的解算结果） 终处理的精度指标值 nrms如下:FROM qchlla.183 FILE:* Summary of biases-free solution *Total P arameters: 1088 live p arameters: 658Prefit nrms: 0.13216E+01 Po stfit nrms: 0.18730E+00 -Uncertainties not scaled by nrmsEnd of loose solution with LC observable and ambiguities freePerforming LC biases-fixed loose solution* Summary of biases-fixed solution *Total p arameters: 1088 live p aramete