GAMIT批处理流程

一、处理前准备

1、在主文件夹内新建test项目文件，项目内新建brdc、igs和rinex三个文件夹，分别存放广播星历，精密星历几观测值文件，所用的命令分别为sh_get_nav、sh_get_orbits和sh_get_rinex（若文件为.Z,用gunzip命令解压，若仍为d,用命令sh_crx2rnx -f 命令解压为o文件）

2、进入test项目文件夹，链接tables,运行 sh_setup -yr 2015 -apr itrf08_com.apr 3、生成station.info文件。 将test/tables下的station.info文件拷贝到rinex文件夹下，打开并编辑，仅保留以#或*开头的前几行，保存并关闭。 打开终端进入rinex文件夹，运行 sh_upd_stnfo -files *.150 4、建立lfile.文件（先验坐标文件）。（非IGS站时，才需此步） 打开终端进入rinex文件夹，用批处理的方式生成lfile.文件，分为三步： (1) 提取观测值文件的先验坐标： grep POSITION *.15o >lfile.rnx

(2) 将.rnx文件转化为.apr文件 rx2apr lfile.rnx 2015 001 (3) 由.apr文件生成lfile.文件

gapr_to_llfile.rnx.aprlfile. \

将station.info和lfile.两个文件复制到test/tables下，覆盖原文件。 5、sestbl.文件的配置。 首先配置sestbl.这里列出主要需要配置的几项内容: Choice of Experiment = BASELINE

如果是解算基线的话建议选择 BASELINE,即固定轨道 Choice of Observable = LC_AUTCLN

选择观测值类型,如果是双频接收机并且是大于 5KM 的基线话建 议选取 LC_AUTCLN,如果是小于 5KM 的基线也可以选择 LC_AUTCLN, 其中 LC_HELP 和 LC_AUTCLN 的差别不大,LC_AUTCLN 略好些。 Zenith Delay Estimation = Y Interval zen = 2

对流层参数的设置,如果是短基线(大概小于 40KM)可以设置

成 N,即不估计对流层天定延迟,采用双差的方法来消除对流层延迟。

如果是大于 40KM 或者 50KM 的基线要设置成 Y。第二项是如果估计对 流层延迟的话,设置多长时间估计一个对流层延迟参数,单位是小时。 Tides applied = 31 Use otl.list = Y Use otl.grid = Y

设置潮汐改正参数,如果为 31 则考虑海潮改正,为 23 则不加海 潮改正,如果考虑海潮改正的话,需要 tables 下有两个海潮文件, otl.list和otl.grid,其中前者是可有可无的,当两者都采用 Y 时, 程序会根据测站位置首先利用otl.list(只包括 IGS 站的海潮模型 参数,以及位于该 IGS 站方圆 10KM 范围内都是一套参数)中的海潮 参数,如果该测站没有出现在otl.list中,则会根据该测站的位置

和otl.grid文件内差出该测站的海潮改正。 以上为sestbl.中主要的设置项,其他项对基线解算的影响很小, 可以采用默认值。 6、sittbl.文件的配置。 下面设置sittbl.文件,即设置各个站坐标的约束量,对其算点 坐标采取紧约束,对非起算点的站采用松弛约束。 执行 visittbl. SITE

打开sittbl.文件。 FIX

--COORD.CONSTR.--

<> ALL NNN30. 30. 30.

<< IGS core stations >>

BJFS BJFS_GPS NNN 0.050 0.050 0.05 LHAZ LHAZ_GPS NNN 0.050 0.050 0.05 TWTF TWTF_GPS NNN 0.050 0.050 0.05 USUD USUD_GPS NNN 0.050 0.050 0.05

对起算点的 XYZ 坐标分量约束为 5cm,非起算点的站采用松弛约束, 约束量为 30 米。

7、sites.defaults文件的配置。 用来控制需要参与解算的测站，可定义不参与结算的测站。 8、process.default文件的配置。 用来控制处理过程中的细节。 二、理由GAMIT解算基线。 用终端进入test项目文件夹，输入批处理命令： sh_gamit -expt test -s 2015 001 008 -orbit IGSF -noftp 参数说明： -d:处理数据的指定日期，-d 2015 001 008 ,处理001和008天的数据 -s:处理数据的指定日期，-s 2015 001 008 ,处理001到008天的数据

结果文件：

qtesta.001:解算记录

otesta.001:解算记录的精简版，一般用此文件 htesta.15001:协方差矩阵、参数平均值

postfitnrms优于0.3时为最佳，大于1时有问题。

本文来源：https://www.bwwdw.com/article/gxto.html