GAMIT使用手册范本

1、范文.范例.参考GAMIT/GLOB软件使用手册一软解介绍GAMIT软件最初由美国麻省理工学院研制，后与美国SCRIPPS海洋研究所共同开发改进。 该软件是世界上最优秀的GSP定位和定轨软件之一，采用精密星历和高精度起算点时，其解算长基线的相对精度能达到10-9量级，解算短基线的精度能优于 1mm,特点是运算速度快、版木更新周期短以及在精度许可范围内自动化处理程度高等，因此应用相当广泛。GAMIT软件由许多不同功能的模块组成，这些模块可以独立地运行。按其功能可分成两 个部分：数据准备和数据处理。此外，该软件还带有功能强大的 shell程序。目前，比较著名的GPS数据处理软件主要有美国麻省理工学

2、院（MIT）和海洋研究所（SIO）联合研制的GAMIT/GLOBK件、瑞士伯尔尼大学研制的 BERNES软件、美国喷气推进实验室（JPL）研制的GIPS丫软件等。GAMIT/GLOBK和BERNES软件采用相位双差数据作为基本解算数据，GIPSY软件采用非差相位数据作为基本解算数据，在精度方面，三个软件没有明显 的差异，都可得到厘米级的点位坐标精度。相比较而言，GIPSY软件为美国军方研制的软件， 国内只能得到它的执行程序，在国内，它的用户并不多，BERNES软件需要购买，它的用户稍微多一点，GAMIT/GLOBK件接近于自由软件，在国内拥有大量用户。GLOBK软件核心思想是卡尔曼滤波（卡尔曼

3、滤波理论是一种对动态系统进行数据处理的有效方法，它利用观测向量来估计随时间不断变化的状态向量），其主要目的是综合处理多元测量数据。GLOBK勺主要输人是经GAMIT处理后的h-file 和近似坐标，当然，它亦己成功 地应用于综合处理其它的GPS软件（如Bernese和GIPSY产生的数据以及其它大地测量和SLR观测数据。GLOBK勺主要输出有测站坐标的时间序列、测站平均坐标、测站速度和多时 段轨道参数，GLOBK可以有效地检验不同约束条件下的影响，因为单时段分析使用了非常宽松的约束条件，所以在 GLOBK中就可以对任一参数强化约束。GAMIT/GLOBK和 BERNES采用双差作为数据分析的基

4、本观测量，它们的缺陷是不能直接解算钟差参数，只能给出测站的基线结果，除测站坐标参数之外， 这些软件还可以解算的参数有：卫星轨道参数、卫星天线偏差、光压参数、地球自转参数、地球质量中心变化、测站 对流层延迟参数、电离层改正参数等，这使这些软件的应用从大地测量学已逐渐延伸到地球 动力学、卫星动力学、气象学以及地球物理学等领域，并取得了很多成果。GAMIT软件的运行平台是 UNIX操作系统，目前，它可在Sun、HF、IBM/RISC、DEC LINUX 等基于in tel处理器的工作站上运行。软件可处理的最大测站和卫星数目可在编译时设定。它的基本输出文件是 H-文件，可作为 GLOBK软件的输入文件

5、，进而估计测站坐标与速度、 卫星轨道参数和地球定向参数。数据处理前，用户需准备所需要的文件，如测站先验坐标文件（L-文件和vg-in文件）、广播星历文件、观测数据文件以及其他辅助文件等。GAMIT每个时段观测数据要求的周期最长为1个UTC天，即从 UTC的0点到24点（北京时间8:0024:00），原则上不要跨天作业。GAMIT软件的组成结构见图1.1，它由不同功能模块组成，主要包括数据准备、生成参考轨道、计算残差和偏导数、周跳检测与修复、最小二乘平差等模块，这些模块即可以单独运行，也可以用批处理命令联在一起运行，最大限度地减少人为操作，提高运算效率。软件的执行程序放在/com、/kf/bin

6、 和/gamit/bin 三个目录下。图1.1 GAMIT/GLOBK软件组成结构图二批处理实例确认GAMIT/GLOB软件安装成功之后，为了能正确处理数据和绘图，还需要做两件事情。 第一，更新表文件。在ftp 服务器上下载最新的tables 表文件（/pub/gamit/tables/），然后复制到安装目录下的tables 中，遇到相同文件选择覆盖即可；第二，安装并配置好绘图工具GMTThe Generic Mapping Tools ）。下面的实例中，处理的数据为 2010年第56到60共5天的3个IGS站（BJFS和WUHN KUNM和2个待求站（

7、chdu、pixi ）的数据。一、处理前的准备1、 在主文件夹内新建test项目文件夹，项目内新建 brdc、igs、rinex三个文件夹， 分别存放当天的广播星历、精密星历和观测值文件。注意：如果文件采用.z压缩格式，则应当使用gun zip命令进行解压；如果解压出的观测值文件仍采用.d的压缩格式，则应使用crx2rnx命令将其转化为.0的标准RINEX格式（sh_crx2rnx -f *.*d可用于.d到.o文件的批量转换）。2、终端进入test项目文件夹内链接tables，运行：sh_setup -yr 20103、 生成 文件。将 test/tables 下的

8、文件拷贝到 rinex 文件夹下，打开并编辑，仅保留以#或*开头的前几行，保存并关闭。打开终端并进入rinex文件夹，运行sh_upd_st nfo -files *.10o运行成功以后stati on .i nfo文件便存放了此项目内的开始、结束时间，站名和接收机、天线类型等信息。4、 建立lfile. 文件。lfile.是测站的先验坐标文件。打开终端并进入rinex文件夹，用批处理的方式生成lfile.文件可以分为三步：a） 提取观测值.o文件的先验XYZ坐标grep POSITION *.10o lfile.rnxb）将.rnx文件转化为.apr文件rx2apr

9、lfile.rnx 2010 056c）由.apr文件生成lfile. 文件gapr_to_l lfile.rnx.apr lfile. 2010 056将和lfile.两个文件拷贝到test/tables文件夹，覆盖原文件。5、 sestbl.的配置。test/tables 下的sestbl.文件是测段分析策略文件，该文件内部 有详细的说明。一般来说采用默认配置即可，通常需要修改的三个地方是：Choiceof Experime nt选择处理方式Choice of Observable选择观测值类型和模糊度解算Use otl.grid选择是否使用潮汐文件对于Choice

10、 of Experiment，选择BASELINE时将固定轨道并在 GAMIT处理中和输出 h- 文件时忽略轨道参数；选择RELAX时将采用松弛解，合并全球IGS h-文件时需要。要想点位置精度高用 RELAX若目的是求基线后面平差则用BASELINE在此实例中采用默认的BASELINE对于Choice of Observable ，选择LC_AUTCLF为采用宽巷模糊度值并用伪距在 autcln 中解算；对于小于几公里的基线，用L1和L2独立载波相位观测值（L1,L2_INDEPENDENT或者仅用LI（ L1_ONLY，相比用无电离层组合（LC_HELP可以减少噪声水平。对于Use otl

11、.grid，由于这里已在ftp上更新tables，有了最新的otl.grid 文件，所以这里选择Y。6、 sittbl.的配置。test/tables 下的sittbl.文件对各个测站的先验坐标（或钟差、大气模型等）进行约束。对高精度的已知坐标采取强约束，而对待求点采用松弛约束。如IGS站的坐标分量约束在较小的175px，对未知点的约束可以到510m=7、 sites.defaults禾口 process.defaults 的配置。sites.defaults文件用来控制需要参与解算的测站。在文档末尾可以根据提示编辑，来给定那些不参与解算的测站或是测站的某些天。process.defaults

12、文件用来控制处理过程中的很多细节，比如sampli ng in terval,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

13、60 ，指的是处理 2010年第56和60天。-s :指定需要处理的时间序列，例如-s 2010 56 60 ，指的是处理 2010年第56到60天。-orbit:卫星轨道类型。-yrext:给日目录前添加年前缀，例如2010_006。-noftp:处理过程中不连接ftp下载数据。-copt :数据处理完成后待压缩的文件类型，例如-copt o q m k x-dopt :数据处理完成后待删除的文件类型，例如-dopt D ao c x结果文件将存放在名称为年积日的文件夹内，此例中，可供参考的结果文件为：sh_gamit_2010_006.summary解算总结qtesta.006解算记录ot

14、esta.006解算记录的简略版，一般关注此文件htesta.10006协防差矩阵、参数平差值在上述文件中，可根据描述来判断解算结果是否符合相应的需求。其中，基线解算结果（O文件，即这里的otesta.006 ）中的postfit_nrms项优于0.3左右时最佳；如果大于1.0,则表示此解存在问题。三、利用GLOB进行平差处理用终端进入test的项目文件夹内，运行：sh_glred -expt test -s 2010 050 2010 065 -yrext -opt H G E & sh_glred.log参数说明：H：运行htoglb，把文件转换为二进制文件G：运行glred，合并文件E:

15、绘图运行成功后，进入 gslon文件夹，以psbase开头的文件即为各个站三维坐标时间序列的图形表示。用终端进入gsoln目录下，依次运行:ls./glbf/h*glx test.gdlglred6 globk_comb.prt globk_comb .log test.gdl globk_comb.cmd 生成的globk_文件中就包含了解算点的三维坐标和相关参数。 至此，利用GAMIT-GLOB解算GPS基线并进行平差的实例叙述完毕。三分布处理实例1. 在桌面上新建一个test的文件夹，在test中新建名称为brdc （广播星历）、igs （存放精密星历）、rinex （下

16、载的观测值文件和自己要解算的观测值文件）。2. 用doy在test的终端里测试，比如输入：doy 2013 02 01查看要解算的是GPS第几周，第几天，以便下载数据用。 Applications Places SystemTerminalFile Ecfit View Terminal Help/root/桌宙丿址听耳 doy 2013 G2 01Qate 213/02/fil 百:盹 hrs, 30V32 ZD 2456324.51900 PJD 56324.SGGSGPS Week I7Z5 Day ofGPS Seconds 432G&0 Qay of Week rriDecimal

17、Year 2613 984932/root/Jg|3. 在终端输入gftp，下载相关的文件，一定要选对自己解算需要的数据。一般下载*.d.z文件（比较小，好下载）。下载完所需要的数据后，要解压：语句：gunzip*.Z4. 新建一个文件夹，比如：名字为032，把要解算的测站 o文件和n文件，以及精密星历文件放进来。把032和test进行如下链接。ab/reat/t54t/0321 links.day 2013 932 test/root/面#性砒/姑2利 In -s ./tables/soltab.2513.J20M5. 打开rinex文件，输入下面的语句，会产生一个lfile.rnx 文件:

18、/roovit/test/rinexv grep pcsition 130nx|；root/246血APPAWt POUTiON KfZllid9320.13Q：附汩2.IT沁W19274,75903139222.2700APPftttXPOSITIONXYZshdo832e.l3a：811733,5B504675M5,9560327539.413APPROXWJSITIOMXYZur砒32tL】ir19Jen 2854輛飾小5,轴却4393315.7337Apf*R&xPOSITIONXV2|6. 输入 rx2apr Ifile.rnx 2013 032,会多出来一个名为：lfile.rnx

19、.apr 的文件,接着输入：gapr_l Ifile.r nx.apr Ifile.2013 032生成一个用纬度和精度来表示点的坐标的文件(lfile.),并把它复制到tables文件夹。7. 在tables文件夹中新建一个文件sitelist，在里面输入参考站点的名称，每行的第一列是空格，从第二行开始输。DjtS lhaz rum shaoj8. 在tabie 文件夹启动一个终端，输入：sh_upd_stnfo -l sitelistgrep： STaiion.ififo: Nc such flic or dtrectory fron siation.iJifo antrl viih si

20、te Iietpm put file is nation.lnfa.ncstates ； 13口彌；眄52 3日阳 MsnNF/Hoglb/minfnecQc cr update siaU file 5： version 2.1 朗 10 的 呢F AL : 133730:952： 3E. 0 bETLNF/ htaq.t；Tst：nf: trrcir cpnincj Jtizn use 1ile si legist ERROR 25TQP FP.T3L Etd： stop fcoiT report stat Output file is staTion. in to.new p

21、 /root/jftc2 . /labLei(9/test/tat h_upd_trtc -1EXELUTJhiG Ell upd fitflfabvina original statiorianfo to o,ona selecting froa entries with site Iist siiclist hurt put file is station .Lnfae*.STATUS ； 13即弼；活53：】臼 MSTlNF/hToglt/mstinf; Merge cr 如曲讥 statLQn.iMo rte s： Version 2.16

22、4 ft?S1ATUS rl3S73ti：*J553： H. 0TLNF/hTogIg/rstxnT: TTwrt 11L be 6 statlen, inf u f1 I p ard ft RTMFX fripariprs rrgw!S1ATUS : 130730:6953 l 14.9 hETIMF/lf DC lb； TStinf: TTitre will bt B 1G5 lug file, ard fnprgp!入:9.会生成一个 .new文件，把它给成statio 保存后，再输sh_upd_st nfo -files ./ri nex/*13o1

23、0. 在032文件中查看 的链接是否正确，然后做一个I文件的链接,打开终端输入：In -s ./tables/lfile ./到此准备工作完成。11.在032文件打开终端，输入：makexp 然后按照以下输入。Erl er -character project cedetestErt r Mh占2d orbit codErter j/carEr t r dy of the ya r 估的 to search all J932Er*e r &eLon number (的 to searcfi all)Eftcr Le . Tilc or 曲f f亠l(? rune ( i

24、 A h i、”*_. m a u 文件名称-:WORD格式整理版-file jauto3.03216. sh_check_sess - sess 32- type jfile检查j-文件/root/M a /test/0J2 sh check sess -sess i2 -type Jfile -fie jautoj.ud217. makex test.makex.batch生产x-文件1/root/桌面/test/e32t-6 makex test.maKex.batch18. fixdrv dtest3.032生成批处理文件19. csh btest3.bat执行批处理脚本，得到 q-文

25、件，即解算结果。ni n a sj _；jjt _j _j_e. u f_i r. _rI. .-r lj- i i wI_l ?. _j ar ui rrdL.cnnj/root/a面/test/殆2毛 csh btest3.bat四精度评估（1）GPS常用评价站坐标精度的指标是多时段基线重复性和多时段坐标重复性。 基线重复性计算模型：(Ri-R 2(5.1)上式中i为观测时段， J为基线s的重复性统计值，二i为第i时段基线s的中误差，Ri为第i时段基线解算结果，N为总的时段数，R为基线结果的加权平均值。有2i(5.2)基线结果的评价还以基线长度与误差的关系来衡量，有两种形式:打=a2 b2

26、L(5.3)(5.4)这里的a、b、c、d四个参数可以根据基线长度和它的重复性统计值 s拟合出来。也有以基线的三个分量分别检验它的重复性的，这三个分量是南北分量sn、东西分量se和高程分量二su。计算公式与基线重复性的公式（5.1 ）相同，只是用基线分量代替基线长度而已。坐标分量重复性计算模型也与（ 5.1 ）相同，只是用坐标分量和坐标分量中误差代替基线长 度和基线长度中误差计算而已。（2）通过天空图评定精度liBD11931BDn 4 hr* UTUTfi-IShraUT如上图：同一个地点在各个的时间段有高残差预示着存在多路径效应。I料血恃期舸和|1士 SmrSJCS Day 335 RMS

27、 6.3 mm如上图：同一个地点在特定的时间段有高残差预示着存在水波折射。（3）通过相位与高度角评定精度tiAfjy HMb=4.6inm error moaei 旷single町户乂t=o.4niniA1 幻J1ir1 A1dbJAx1 dI 4R PS*如上图，错误：整体趋势线在有的高角度离开蓝色线，表明错误引用了天线相位模型，或卫星信号受到了电磁干扰。如上图，正常情况下：红色整体趋势线平顺，在蓝色线上下均匀波动。(4) 查看 qexpta 文件：qexpta.doy-Adjust(units)!prioriAdjust (cd)Form且丄Fract1 *BJFSGEOC LATdms0,00140*00600 + 2 -2 * B-JESGEOC WNGdirts-0.00710x0062-1.13 BJFSRADIUS-o m0m00fi6一GEOC LATdm吕0+D0500+0075