gamit学习总结

1、gamit 学习总结 打开 LINUX 系统，登陆。（选择安装语言包，关键更新 ） 重启。语言包的设置在 system settings-language supportinstall languages 选 择 chinese ( simplified ) 点 击 apply changes 等待系统自动安装， 重启，重复步骤进入 support 选 择系统和菜单语言 汉语 ok（文件夹尽量保持英文状态，以 便在终端中输入）； 2. GAMIT 安装包中没有包含海洋潮汐模型 otl.grid 和 map vmf1 函数模型需要自己下载连接。 otl.grid 在 ftp:/garner.uc

2、/archive/garner/gamit/tables/ 中下载 （约 696M ） vmf1 在 /pub/GRIDS/ 中下载（约 220M ） 如没有下载 otl.grid 模型可以通过设置禁用该模型。 模型的禁用： 打开 /opt/gamit/tables/sestbl. 将 use otl.grid 的之改为 N，表 示禁用海洋潮汐模型，同时将 Tides spplied 改为 23（原始为 31）。 独立下载模型的连接： 进入 /opt/gamit/tables 运行终端 运行： ls -al 进行连接 查看 otl 与 map

3、 的连接情况 rm otl.grid 移除已有的连接 ln -s /opt/otlmap/otl_FES2004.grid otl.grid 进行连接 ( /opt/otlmap/otl_FES2004.grid 为模型的路径) rm map.grid ln -s /opt/otlmap/vmf1grd.2012_063 map.grid （map.grid 目 前有 vmf 提供模型） ls -al 连接并查看是否已连接成功。 连接结束。（连接完成后安装更新时会自动连接） 其它 otl 模型： otl_FES99.grid 182.7 Mb otl_CSR4.grid 45.6 Mb otl

4、_GOT00.grid 45.6 Mb otl_NAO99b.grid 87.1 Mb Atmosphere 模型也许年更新： atml files 185 Mb/year 3. sites.defaults 准备介绍 sites.defaults 文件是站处理控制文件。格式如下： all_sites expt xstinfo glrepu brus_gps expt ftprnx ftpraw localrx xstinfo xsite glreps glts 其中 ftprnx ftpraw localrx xstinfo xsite 为 gamit 选项， glrepu glreps g

5、lts 为 GLOBK 选项，格式具体含义如下： expt 为 4 个字符的项目名称（工程名称） ftprnx 表示该站点的 rinex 数据从 FTP 站点下载 ftpraw 表示该站点的 raw 数据从 FTP 站点下载 localrx 表示站点的 rinex 数据位本地数据 xstinfo 表示该站点不自动更新 表中的信息 xsite 该站点的所有天或指定日期的数据不处理 glrepu 表示站点用于计算重复解（默认全部） glreps 用于定义站点的在 GLRED 处理中的参考框架（默认 为 IGS 站点列表） glts 在 GLRED 处理中以时间序列绘制站点

6、图（默认全部） 一般通过自己编制的 sites.defaults 文件可以设置所有站点不 自动更新 表中的信息。示例如下： all_sites test xstinfo （ test 为工程名） LHAS_GPS test localrx glrepu glreps SHAO_GPS test localrx glrepu 4. 文件准备 文件可以采用 sh_upd_stnfo 命令生成也可以人 工编辑得到。 对于部分不规范的 rinex 文件建议采用人编辑的方法，并参 照 gamit 的 rcvant （天线与接收

7、）文件进行编辑，确保后期 处理的顺利。 在更新前应将 sites.defaults 文件中的 all_sites expt xstinfo 注 释掉，否则容易报错，无法更新。 采用 sh_upd_stnfo 命令生成 文件方法：将 tables 中的 文件拷贝至 rinex 中，在 rinex 文 件夹下运行 sh_upd_stnfo -l sd 提取 IGS 站点的信息，得到 .new 然 后 见 .new 改 名 为 ， 在 运 行 sh_upd_stnfo -fi

8、les *12o -expt test 即可。 采用手工编辑的办法： 会自 将 tables 中的 文件拷贝至 rinex 中，在 rinex 文 件夹下运行 sh_upd_stnfo -l sd 提取 IGS 站点的信息。 动生成 .new 文件。 息。 打开 .new 删除无关条目，增加其他站点的信 此时应打开 rinex 文件进行核对，确保无误后将其保存 。应注意将天线高参考设为 DHARP （天线参考 中心一般为天线底部） 格式： Fixed format file with stati

9、on information ? *TRCK SIT E Station Name Ant Ht Ant N Ant E ? 7ODM 7ODM Seven Oaks Dam 0.0083 0.0000 0.0000 ? Rcvr AntCod HtCod Vers ? ASHZ12 ATDMRB DHPAB 9.10 ? Year Doy SN Start Stop ? 2000 110 0 00 00 00 24 00 00 ? ( = line break for this display) This file must be sorted by TRCK, then epoch TRC

10、K and SITE 4-character code match the 4-char code in RINEX file name TRCK = SITE for static analysis在静态分析中 Change for kinematic analysis对于动态分析要改变 File must be sorted by this field, then epoch ? station_name is a long (12-char) description AntHt Raw antenna height as measured原始测量高度 Critical measureme

11、nt!至关重要 ? HtCod Defines geometry of AntHt measurement定义几何类型 DHPAB is RINEX standard is vertical height to antenna reference point (ARP) 度) 定义几何类型 DHPAB 指相对于站 点的垂直高 ? N and E offsets基准偏移 e.g., uncalibrated tribrach rarely used occupied point w.r.t. benchmark RCVR 6-character codes for the receiver ty

12、pe not critical except for some old receivers or buggy firmware versions version number of firmware in the receiver接收机固件版本号 Should be correct to handle SNR coding properly AntCod and HtCod are critical关键 Sets phase center model涉及相位中心模型 conversion of antenna height measurement天线高度的解译 SLBGP - Slant bo

13、ttom ground plane倾斜的底部大地平面 ? year and doy specify year and day of year when the values apply ? times set part of day (if it changes during UTC day) ? This representation will be augmented in future version of GAMIT 5. sittbl. 编辑 打开 sittble 文件删除无关站点，增加工程站点设置各站点的 约束属性。 对于 IGS 框架点或你需要约束的点设置为 0.050 0.05

14、0 0.050， 其他点 100. 100. 100. 一般来说约束点为 IGS 站点且不少于 3 个。 示例： SITE FIX -COORD.CONSTR.- BJFS BJFS_GPS NNN 100. 100. 100. GUAO GUAO_GPS NNN 100. 100. 100. LHAZ LHAZ_GPS NNN 100. 100. 100. SHAO SHAO_GPS NNN 100. 100. 100. URUM URUM_GPS NNN 0.050 0.050 0.05 WUHN WUHN_GPS NNN 0.050 0.050 0.05 KUNM KUNM_GPS NN

15、N 0.050 0.050 0.05 下载 sittbl.refined 文件 ,里面有最近的各站的可靠性的一个 评测。 要想点位置精度高则设置 IGS 站的约束小一些 ,即认为 IGS 站已经很准 ,强约束 ,若目的是求基线后面平差则可以将 IGS 站的约束设置大些。 6. process.defaults 文件可修改部分的解释 process.default 文件中需要修改的有采样间隔， 最大历元数 ， 开始时间 等，示例如下： set sint = 30 采样间隔依据需要和 rinex 的情况可修改为 15 30 60 set nepc = 2880 最大历元数 依据观测时长计算调整。

16、set stime = 0 0 开始时间（小时 分）依据 rinex 文件修改， 一般要求时间与 rinex 时间差不大于 0.1h 否则报错。 set use_rxc = N 设置为 N 不依据 rinex 更新 lfile 文件。 set aprf = itrf08.apr 采用 08 框架数据 （默认 itrf2000 ） 如不 采用 GLOBK 平差可不设置。 set minxf = 300 Xfile 文件大小的最小值， 小于该值的将被抛 弃。依据历元数设置。历元数少时应减小该值。 set mailto = * 设置为有效邮箱以便接受 gamit 7. sestbl.文件依据实际需要

17、进行修改。示例如下： hoice of Experiment = BASELINE (计算高精度基线) BASELINE/RELAX./ORBIT 默认为 BASELINE ，但一般修 改为 relax.(BASLINE 为固定轨道)(如只需要基线文件可以 选择 BASELINE ) Choice of Observable = LC_AUTCLN 观测值类型选择 LC_AUTCLN (default), LC_HELP (codeless L2), L1_ONLY (L1 soln from dual freq),L2_ONLY (L2 soln from dual freq),L1,L2_

18、INDEPENDENT (L1 + L2 from dual freq), L1L1_RECEIVER (must add L1only in autcln.cmd) 对于 15km 以下的短基线应选择 L1,L2_INDEPENDENT 或 L1_ONLY Type of analysis =0-iter 迭代多少次 计算结果作为初始值再 次计算一般设置 1-iter Zenith delay =N 天顶延迟 对于短基线选 N Interval zen = 2 ;2 hrs = 13 knots/day (default is 1 ZD per day) 天顶延迟修改为 1h，默 认 2h

19、Update T/L files = L_ONLY ;T_AND_L (default), T_ONLY, L_ONLY, NONE人工编辑时可设为 NONE # DMap = GMF ; GMF(default)/NMFH/VMF1 如果已下载 VMF1 这可设为 VMF1-重新连接。 # WMap = GMF; GMF(default)/NMFW/VMF1 如果已下载 VMF1 这可设为 VMF1 Use map.grid = Y;VMF1 grid file with mapping functions and ZHD Inertial frame = J2000 ; B1950/J20

20、00 (default = J2000) 8. sh_gamit 批处理解释 sh_gamit_test.log 件 sh_gamit_test.log 当采用快速星历处理时需要设置 生成过程详解文 -orbit IGSU/IGSR #要处理最近的数据需要将 gamit/gamit/fixdrv 文件夹下， 打开 armake.f 文件，将 2010 改为 2100.然后再 gamit 文件夹 下运行 ./install_software. 即重新编译安装。 （未验证） -mailto 联网后发送 Basic example ? sh_gamit -d 1998 240 -orbit SIOF

21、 -expt emedeops usno ? process GPS data from day 240 of 1998 (Aug 28) ? use any data in the local rinex directory ? ftp dat a for sites specified by the user ? optionally scan multiple directories to find data ? use Final orbit estimate from Scripps Institution ofOceanography (SIO) 轨道估计数据来源 ? use EO

22、P data from US Naval Observatory (USNO) IERSBulletin A EOP 数据来源 ? Name the run “ emed 工程名 identify files select data List options sh_gamit -help ? Select days of data to process ? -d yr days where yr and days are the data to be processed e.g., -d 1997 153 156 179 ? -s yr d1 d2 where d1 is start day

23、and d2 is stop day to be processed e.g. 1997 153 178 ? -r day where day is number of days before current date. sh_gamit options ? -expt 4 char name of experiment/solution being run Default expt ? -orbit Type of orbit to use. Default SIOF. Optional second argument gives the day-of-year for the orbit

24、if not the same as processing day ? -orbt 4 char ID for orbit files. Default 4 char extracted from the orbit selected ? -eops Name of EOP series to use. (bull_b, bull_a, etc.) Default usno -sessinfo Session - sampling interval, number of epochs, start time (default 30 2880 0 0) ? -rx_doy_plus n Sear

25、ch RINEX file names this many days into the future for current-day data ? -rx_doy_minus n Search RINEX file names this many days into the past for current-day data ? -xsite list List of sites to be excluded from the processing.解算排除的站 寻找往前推迟 n 天的 rinex 寻找往后 n 天的 rinex 数据 数据 ? -netext add network suff

26、ix to day directories. i.e. 035r ? -yrext ? add year prefix to day directories. Ie 1997_053 ? -pres ? Plot gamit residuals as skyplots . YES, NO, ? -nogifs ? Do NOT create skyplot GIF graphics files 日期文件夹添加网的后缀 日期文件夹添加年前缀 ELEV. Default N 绘制天空残差图 不产生天空 GIF 图片 -copt file types Compress remaining files

27、 for each dayserve when d ? -aopt file types Copy the specified files for each day to the directory archivepth 文件夹。 ? -dopt file types Delete specified files after processing each day. Default c -mailto Default x k autcln.out 压缩指定类型的每天文件 specified by archivepth将指定每天文件复制到 删除指定类型文件 Name of user and ma

28、chine to send results summary. (e.g., ) -snrres Calculate theoretical phase residuals from SNR. Requires raw data or RINEX with valid SNR values. ? -noacc Do not exclude satellites base on accuracy codes in the .sp3 file. Default exclude不排除 .sp3 中基于精度码 方式排除卫星。 ? -xver char 1-char

29、acter x-file version (6th character of x-file). Default is single-digit year ? -noftp Use this option if no ftp download connections are available.不联网 ? -remakex Y/C/N ? Y - Remake x-files without checking. ? C - Remake x-file if orbit and are incompatible and all RINEX files are availa

30、ble Default C ? N - Dont remake x-files under any circumstances 若tables 中 消失， sh_gamit 将会使用 rinex 的头文件信息。 guess_rcvant.dat 来将 rinex 的头文件信息翻译成 gamit 可阅读形式。 Antenna height errors 在 GLOBK 可降低，但类型识别错误无法消除。 9. 打开日期目录查看结算结果 判断解算质量主要有两个方面： 一是数据利用率 （即是否有足够多的可用于进行合理估计的数据参与计算） ； 二是数据的噪声水平是否与选用的模型匹

31、配。 （这一项主要有两个指标， 基线分量误差的大小和 nrms（ 标 准化的均方根误差 ），对于基线分量误差的判断主要依据预先 设定的站点坐标和轨道参数的先验约束，如果大于先验可以 通过查看 Q 文件或者 autcln.sum 文件，是否有大量被 autcln 丢弃的数据。对于 nrms 一般应在 0.25 左右。如果大于 0.5就意味着周跳没有消除，或出现了模型错误（固定站点坐标 错误或者未模型化的卫星） 等。解算结果文件 autcln.sum.post 中的 rms 残差为显示控制网站点的相对质量。 打开 sh_gamit_doy.summary ，otesta.doy, qtesta.d

32、oy 文件查 看解算总结和基线解算情况，解算成功与否判断： Postfit nrms: *E+00 是否在 0.25 左右， 最差不能大于 0.5，基线分量的误差大小是否小于预设值。 至此基线解算全部结束。 基线解算结果： sh_gamit_doy.summary 基线解算结果： qexpta.doy,oexpta.doy,hexpta.yrdoy 10. 若软件限制了最大测站数据，最大卫星数目，最大对流层 延迟参数，最大解算历元数目。更改于 Makefile.config 文件。然后重新编。 11. tables 表介绍 12. 格式问题 o 文件名需要小写，并且内容格式符合 UINX 文件

33、名改小写： sh_casefold -dir d -files * 格式转 UNIX: sh_dos2unix * 13. gamit 模块以及结构 14. Data Types ? L1 and L2 phase used for parameter estimation ? P1 and P2 range measurements used to help clean phase data ? used to estimate receiver clock offset with tolerance of 1 microsecond (= 300 m) ? LC phase proport

34、ional to geometric range sensitive insensitive to ionospheric delay linear combination of L1 and L2 ? LG phase measures ionospheric delay only insensitive to geomery ? Widelane WL linear combination of L1, L2, P1 and P2 Melbourne-Wubbena Should be constant (+/- noise) 15. Aim of GAMIT Basic aim of G

35、AMIT is to generate positions/orbits for one day of data. Standard product of GAMIT is a minimally constrained solution 16. gg/tables/rcvant.dat File contains all receivers and antennas known to GAMIT ? Translates user -specified 6-character code to IGS SINEX standard. ? The file in templates direct

36、ory provides mapping fromRINEX header information to 6-char codes. ? Example lines rcv ASHTECHZ-XII ASHZ12 rcv Z-XII3 ASHZ12 swv 5F 1.00 ? Explanation given at top of file ? used a represent space ? Sometimes need to be careful that pattern correctly gets correct type. Should not be a problem for IG

37、S rinex files ? default lines allowed, but be careful 17. 模型说明： 1.otl 潮汐改正 2. vmfl GMF 投影函数 可消除周跳波动， 3. atml 大气荷载模型， 对高程影响较大， 可靠性需要进一步证实 4. atl 大气抄袭荷载模型和 met 气象模型 下载地址： 1) 海洋潮汐。例如 otl_FES2004.grid 放在软件 talbels目录下 链接到 otl.grid 。 2) 大气负荷。例如 atmldisp_cm.2006 ，每年更新一次。连接 到 atmldisp_g

38、rid.2006 3) vmfl 投 影 函 数 栅 格 ， 例 如 vmflgrd2006 ， 连 接 到 map.grid.2006 。每年更新一次。 18. 解算所需文件说明 类型文件名说明： 数据文件 bjfs0900.06o bjfs 为站点代码 ， 090 为年积日 ， 0为时段号 ， 06 为2006 年 ，o 为观测文件 auto0900 06n 广播星历文件 ，前 4 个字节改为站点代码则为站点广播星历文件 igs13685 sp3 igs 为事后精密星历 ，1368 为 GPS周 ，5为星期五 igr13685 sp3 igr 代表快速精密星历 igu13685 sp3 i

39、gu 代表预报精密星历 配置文件 lgigg6 090 l station Oct session Info sestbl sittbl 代表 L 文件 ， gigg 代表工程名 ， 6代表2006 年， 090 为年积日 测站信息文件 ， 记录各测站接收机和天线类型 时段信息文件 ， 可以在数据处理过程中自动创建 解算参数控制文件 ，设置与解算有关的参数 ， 达到不同的解算目的 设置测站坐标约束 、 大气模型等参数 表文件 ut1 国际时间系统表 pole 极移表 leap sec 跳秒表 nutabl 章动表 luntal 月亮表 soltab 太阳表 gdetic dat 地球形状参数表

40、 rcvant dat 接收机及其天线型号对照表 svnav dat 星号对照表 svs_exclude Dat 坏卫星信息表 hi dat 仪器高参数表 dcb dat 差分码偏差表 autcln cmd 周跳自动探测和修改命令表 antmod dat 天线相位中心改正表 ot1 list /grid全球格网或各测站的海洋潮汐参数表 19 .crx2rnx压缩的 d文件转换成 o文件 Usage: crx2rnx file - -f -s -h stdin and stdout are used if input file name is not given. - : output to s

41、tdout -f : force overwrite of output file -s : skip strange epochs (default:stop with error) This option may be used for salvaging usable data when middle of the Compact RINEX file is missing. The data after the missing part, are, however, useless until the compression operation of all data are init

42、ialized at some epoch. Combination with use of -e option of RNX2CRX (ver.4.0 or after) may be effective. Caution : It is assumed that no change in # of data types happens in the lost part of the data. -h : display help message exit code = 0 (success) = 1 (error) = 2 (warning) 20某个人的参数设置： 表2 解算控制参数的选

43、择 控制选项控制参数 处理模式 RELAX 迭代方案 1-ITER 观测量 LC_AUTCLN 天顶延迟参数 25 Number Zen = 25; number of zenith-delay parameters (default 1) 梯度参数 3 Number gradients 海潮模型 otl.list 光压模型 BERNE 惯性框架 J2000 0 对流层误差模型 SAAS 截止高度角 /( ) 10 Elevation Cutoff 测站约束 ( BJFS) /m 0 05 /0 05 /0 1 测站约束 ( WUHN) /m 0 025 /0 01 /0 1 测站约束 ( K

44、UNM) /m 9 999 /9 999 /9 999 9 其他默认 某篇论文写到 ; 不同截止高度角对基线精度的影响 表4 不同截止高度角对基线 精度指标 (BJFS_KUNM精) 度的影响 高度角 /( nrms 0 241 04 0 241 04 0 231 09 0 222 52 0 213 98 X /m 0 041 9 0 041 9 0 045 0 0 048 5 0 051 2 Y /m 0 040 1 0 040 1 0 043 4 0 049 1 0 056 3 Z /m 0 044 5 0 044 5 0 047 8 0 051 1 0 055 1 L /L 4.9 0

45、10 E-9 4 90 10 E-9 5 49 10 E-9 6 53 10 E-9 7 67 10 ) 5 10 15 20 25 E-9 不过高度角太低时路径延迟对水汽变化非常敏感 因此在 般的解算情况下选择 10的截止高 度角比较适宜 。 不同天顶延迟参数对基线精度的影响 表5 不同天顶延迟参数对基线 (BJFS_KUNM精) 度的影响 精度指标 天顶延迟参数 7 13 19 25 Nrms 0 239 47 0 242 77 0 241 66 0 241 04 X /m 0 042 5 0 041 9 0 041 8 0 041 9 0 040 8 Z /m 0 045 9 L /L

46、4 95 10 E-9 5 05 10 E-94 90 10 E-9 4 90 10 E-9 一般情况下 ，天顶延迟参数取为 13，即每2 h 产生一个天顶 延迟估值是比较合适的 。 天顶延迟参数变化对长基线的精度 Y /m 0 040 2 0 044 7 0 039 8 0 044 2 0 040 1 0 044 5 影响不是非常明显 选取不同观测量对基线精度的影响 表6 选取不同观测量对基线 （BJFS_KUNM精） 度的影响 精度指标 观测量 LC _AUTCLN L1L2 LC_HELP LC_ONLY nrms 0 241 04 0 239 28 0 239 27 0 239 27

47、X /m 0 041 9 0 047 3 0 047 9 0 047 9 Y /m 0 040 1 0 044 1 0 044 6 0 044 6 Z /m 0 044 5 0 053 8 0 055 3 0 055 3 L /L 4 90 10 E-9 6 53 10 E-9 6 73 10 E-9 6 73 10 E-9 LC_AUTCLN（观 测值是无电离层的线性组合 ） 精度指标 测站约束 /m 0.01/0.01/0.02 0.05/0.05/0.06 0.1/0.1/0.1 0.5/0.5/0.5 nrms 0 239 38 0 241 18 0 241 07 0 240 64 X /m 0 013 3 0 037 2 0 047 4 0 119 4 Y /m 0 019 9 0 034 8 0 040 1 0 081 1 Z /m 0 015 3 0 036 0 0 046 1 0 119 9 L /L 3 66 10 E-9 4 41 10 E-9 4 75 10 E-9 8 02 10