GAMIT使用手册范本.docx

GAMIT/GLOBK软件使用手册一 软解介绍GAMIT软件最初由美国麻省理工学院研制 , 后与美国 SCRIPPS海洋研究所共同开发改进该软件是世界上最优秀的GSP定位和定轨软件之一,采用精密星历和高精度起算点时,其解算长基线的相对精度能达到10-9量级 ,解算短基线的精度能优于1mm,特点是运算速度快、版木更新周期短以及在精度许可范围内自动化处理程度高等, 因此应用相当广泛GAMIT软件由许多不同功能的模块组成 , 这些模块可以独立地运行 按其功能可分成两个部分 : 数据准备和数据处理此外 , 该软件还带有功能强大的 shell 程序目前，比较著名的 GPS数据处理软件主要有美国麻省理工学院 （ MIT）和海洋研究所（ SIO）联合研制的 GAMIT/GLOBK软件、瑞士伯尔尼大学研制的 BERNESE软件、美国喷气推进实验室（ JPL）研制的 GIPSY软件等 GAMIT/GLOBK和 BERNESE软件采用相位双差数据作为基本解算数据， GIPSY软件采用非差相位数据作为基本解算数据，在精度方面，三个软件没有明显的差异，都可得到厘米级的点位坐标精度 相比较而言， GIPSY软件为美国军方研制的软件，国内只能得到它的执行程序，在国内，它的用户并不多， BERNESE软件需要购买，它的用户稍微多一点， GAMIT/GLOBK软件接近于自由软件，在国内拥有大量用户。

GLOBK 软件核心思想是卡尔曼滤波（卡尔曼滤波理论是一种对动态系统进行数据处理的有效方法 , 它利用观测向量来估计随时间不断变化的状态向量） ，其主要目的是综合处理多元测量数据 GLOBK的主要输人是经 GAMIT处理后的 h-file 和近似坐标 , 当然 , 它亦己成功地应用于综合处理其它的 GPS软件（如 Bernese 和 GIPSY）产生的数据以及其它大地测量和SLR观测数据 GLOBK的主要输出有测站坐标的时间序列、测站平均坐标、测站速度和多时段轨道参数， GLOBK可以有效地检验不同约束条件下的影响 , 因为单时段分析使用了非常宽松的约束条件，所以在 GLOBK中就可以对任一参数强化约束GAMIT/GLOBK和 BERNESE采用双差作为数据分析的基本观测量， 它们的缺陷是不能直接解算钟差参数， 只能给出测站的基线结果， 除测站坐标参数之外， 这些软件还可以解算的参数有：卫星轨道参数、卫星天线偏差、光压参数、地球自转参数、地球质量中心变化、测站对流层延迟参数、 电离层改正参数等， 这使这些软件的应用从大地测量学已逐渐延伸到地球动力学、卫星动力学、气象学以及地球物理学等领域，并取得了很多成果。

GAMIT软件的运行平台是 UNIX 操作系统， 目前，它可在 Sun、HP、IBM/RISC、DEC、LINUX等基于 intel 处理器的工作站上运行软件可处理的最大测站和卫星数目可在编译时设定它的基本输出文件是 H-文件，可作为 GLOBK软件的输入文件，进而估计测站坐标与速度、卫星轨道参数和地球定向参数 数据处理前， 用户需准备所需要的文件， 如测站先验坐标文件（ L- 文件和 vg-in 文件）、广播星历文件、观测数据文件以及其他辅助文件等 GAMIT每个时段观测数据要求的周期最长为 1 个 UTC 天，即从 UTC 的 0 点到 24 点（北京时间8:00~24:00 ），原则上不要跨天作业GAMIT软件的组成结构见图 1.1 ，它由不同功能模块组成，主要包括数据准备、生成参考轨道、计算残差和偏导数、周跳检测与修复、最小二乘平差等模块， 这些模块即可以单独运行， 也可以用批处理命令联在一起运行，最大限度地减少人为操作，提高运算效率 软件的执行程序放在 /com、 /kf/bin和 /gamit/bin三个目录下binGen_utilhtoglbglinitcomglobkUpdate_pmuglorgghandlerskfblsumgloutplotincludeshelpglistutilsexampleglbak...globctrackGAMIT/GLOBKlibrariescomlibmatrixarcfixdrvtablescleanlibtemplatesficamakexpincludessolvegamitmakexcfmrgorbitsexamplemapsutilshibinmodelctoxtform图 1.1 GAMIT/GLOBK 软件组成结构图二 批处理实例确认 GAMIT/GLOBK软件安装成功之后，为了能正确处理数据和绘图，还需要做两件事情。

第一，更新表文件在ftp服务器上下载最新的tables表 文 件（），然后复制到安装目录下的tables中，遇到相同文件选择覆盖即可；第二，安装并配置好绘图工具GMT（ The GenericMapping Tools ）下面的实例中， 处理的数据为 2010 年第 56 到 60 共 5 天的 3 个 IGS 站（ BJFS、和 WUHN、KUNM）和 2 个待求站（ chdu、pixi ）的数据一、处理前的准备1、在主文件夹内新建 test 项目文件夹，项目内新建 brdc 、igs 、 rinex 三个文件夹，分别存放当天的广播星历、精密星历和观测值文件注意：如果文件采用 .z 压缩格式，则应当使用 gunzip 命令进行解压；如果解压出的观测值文件仍采用 .d 的压缩格式，则应使用 crx2rnx 命令将其转化为 .o 的标准 RINEX 格式（ sh_crx2rnx -f *.*d可用于 .d 到 .o 文件的批量转换） 2、终端进入 test 项目文件夹内链接 tables ，运行：sh_setup -yr 20103、生成 station.info 文件。

将 test/tables 下的 station.info 文件拷贝到 rinex 文件夹下，打开并编辑，仅保留以 #或 * 开头的前几行，保存并关闭打开终端并进入 rinex 文件夹，运行sh_upd_stnfo -files *.10o运行成功以后 station.info 文件便存放了此项目内的开始、结束时间，站名和接收机、天线类型等信息4、建立 lfile. 文件 lfile. 是测站的先验坐标文件打开终端并进入 rinex 文件夹，用批处理的方式生成 lfile. 文件可以分为三步：a）提取观测值 .o 文件的先验 XYZ坐标grep POSITION *.10o > lfile.rnxb）将 .rnx 文件转化为 .apr 文件rx2apr lfile.rnx 2010 056c）由 .apr 文件生成 lfile. 文件将 station.info 和 lfile. 两个文件拷贝到 test/tables 文件夹，覆盖原文件5、 sestbl. 的配置 test/tables 下的 sestbl. 文件是测段分析策略文件，该文件内部有详细的说明一般来说采用默认配置即可，通常需要修改的三个地方是：Choiceof Experiment 选择处理方式Choice of Observable 选择观测值类型和模糊度解算Use otl.grid 选择是否使用潮汐文件对于 Choice of Experiment ，选择 BASELINE时将固定轨道并在 GAMIT处理中和输出 h-文件时忽略轨道参数；选择 RELAX时将采用松弛解，合并全球 IGS h- 文件时需要。

要想点位置精度高用 RELAX；若目的是求基线后面平差则用 BASELINE在此实例中采用默认的BASELINE对于 Choice of Observable ，选择 LC_AUTCLN为采用宽巷模糊度值并用伪距在 autcln中解算；对于小于几公里的基线，用 L1 和 L2 独立载波相位观测值（ L1,L2_INDEPENDENT）或者仅用 LI （ L1_ONLY），相比用无电离层组合（ LC_HELP）可以减少噪声水平对于 Use otl.grid ，由于这里已在 ftp 上更新 tables ，有了最新的 otl.grid 文件，所以这里选择 Y6、 sittbl. 的配置 test/tables 下的 sittbl. 文件对各个测站的先验坐标（或钟差、大气模型等） 进行约束 对高精度的已知坐标采取强约束， 而对待求点采用松弛约束 如 IGS站的坐标分量约束在较小的 1～75px，对未知点的约束可以到 5～ 10m7、 sites.defaults 和 process.defaults 的配置sites.defaults 文件用来控制需要参与解算的测站在文档末尾可以根据提示编辑，来给定那些不参与解算的测站或是测站的某些天。

process.defaults 文件用来控制处理过程中的很多细节，比如 sampling interval,number of epochs, start time for processing ， default globk .apr file 等等，根据需求和提示进行编辑在这里均使用它们的默认值二、利用 GAMIT解算基线用终端进入 test 的项目文件夹，输入批处理命令进行解算：sh_gamit -expt test -s 2010 056 060 -noftp -dopt D ao c x >& sh_gamit.log参数说明：-expt ：指定四个字符的项目名称-d ：指定需要处理的指定日期，例如 -d 2010 56 60，指的是处理 2010 年第 56 和。