2022年gamit学习总结解读,推荐文档

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

2、mit/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 的连接情况rm otl.grid 移除已有的连接ln -s /opt/o

3、tlmap/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_got00.grid 45.6 mb otl_nao99b.grid 87.1 mb atm

4、osphere 模型也许年更新：atml files 185 mb/year 3.sites.defaults 准备介绍sites.defaults 文件是站处理控制文件。格式如下：all_sites expt xstinfo brus_gps expt ftprnx ftpraw localrx xstinfo xsite glrepu glreps glts 其中ftprnx ftpraw localrx xstinfo xsite为 gamit 选项，glrepu glreps glts 为 globk 选项，格式具体含义如下：expt 为 4 个字符的项目名称（工程名称）ftprnx

5、表示该站点的rinex 数据从 ftp 站点下载ftpraw 表示该站点的raw 数据从 ftp 站点下载localrx 表示站点的rinex 数据位本地数据xstinfo 表示该站点不自动更新 表中的信息xsite 该站点的所有天或指定日期的数据不处理glrepu 表示站点用于计算重复解（默认全部）glreps 用于定义站点的在glred 处理中的参考框架（默认为igs 站点列表）glts 在 glred 处理中以时间序列绘制站点图（默认全部）一般通过自己编制的sites.defaults 文件可以设置所有站点不自动更新 表中的信息。示例如下：

6、all_sites test xstinfo （test 为工程名）lhas_gps test localrx glrepu glreps 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 14 页 - - - - - - - - -shao_gps test localrx glrepu 4. 文件准备 文件可以采用sh_upd_stnfo 命令

7、生成也可以人工编辑得到。对于部分不规范的rinex 文件建议采用人编辑的方法，并参照gamit 的 rcvant（天线与接收）文件进行编辑，确保后期处理的顺利。在更新前应将sites.defaults 文件中的all_sites expt xstinfo 注释掉，否则容易报错，无法更新。采用 sh_upd_stnfo 命令生成 文件方法：将 tables 中的 文件拷贝至rinex 中，在 rinex 文件夹下运行sh_upd_stnfo -l sd 提取 igs 站点的信息， 得到.new 然后见 station.in

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

9、格式： fixed format file with station information ? *trck sitestation 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 ? yeardoy sn start stop ? 2000 110 0 00 00 00 24 00 00 ? ( = line break for this display) this file must be so

10、rted by trck, then epoch trck 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 mea

11、sured 原始测量高度 critical measurement! 至关重要? 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 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共

12、 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共 14 页 - - - - - - - - - occupied point w.r.t. benchmark rcvr 6-character codes for the receiver type not critical except for some old receivers or buggy firmware versions ? vers version number of firmware in the receiver 接收机固件

13、版本号 should be correct to handle snr coding properly antcod and htcod are critical 关键 sets phase center model 涉及相位中心模型 conversion of antenna height measurement 天线高度的解译 slbgp - slant bottom ground plane 倾斜的底部大地平面? year and doy specify year and day of year when thevalues apply ? times set part of day (

14、if it changes during utc day) ? this repres entation will be augmented in future version of gamit 5.sittbl. 编辑打开 sittble 文件删除无关站点，增加工程站点设置各站点的约束属性。对于 igs 框架点或你需要约束的点设置为0.050 0.050 0.050，其他点100. 100. 100. 一般来说约束点为igs 站点且不少于3 个。示例：site fix -coord.constr.- bjfs bjfs_gps nnn 100. 100. 100. guao guao_gps

15、 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 nnn 0.050 0.050 0.05 下载 sittbl.refined 文件 ,里面有最近的各站的可靠性的一个评测。要想点位置精度高则设置igs 站的约束小一些,即认为igs 站已经很准 ,强约束 ,若目的是求基线后面平差则可以将igs 站的约束设置大些

16、。6.process.defaults 文件可修改部分的解释process.default 文件中需要修改的有采样间隔，最大历元数，开始时间等，示例如下：set sint = 30 采样间隔依据需要和rinex 的情况可修改为15 30 60 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 14 页 - - - - - - - - -set nepc = 2880 最大历元数依据观测时长计算调

17、整。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.文件依据实际需要进行修改。示例如下：

18、hoice of experiment = baseline （计算高精度基线）baseline/relax./orbit默认为 baseline ，但一般修改为relax.(basline为固定轨道）（如只需要基线文件可以选择 baseline ）choice of observable = lc_autcln 观测值类型选择lc_autcln (default), lc_help (codeless l2), l1_onl y (l1 soln from dual freq),l2_onl y (l2 soln from dual freq),l1,l2_independent (l1 +

19、 l2 from dual freq), l1&l2 (same as l1,l2_independent but with ion constraint);l1_receiver (must add l1only in autcln.cmd) 对于 15km 以下的短基线应选择l1,l2_independent 或 l1_onl y type of analysis =0-iter 迭代多少次 计算结果作为初始值再次计算一般设置1-iter zenith delay =n 天顶延迟对于短基线选n interval zen = 2 ;2 hrs = 13 knots/day (defa

20、ult is 1 zd per day) 天顶延迟修改为1h，默认 2h update t/l files = l_only ;t_and_l (default), t_onl y, l_onl y, 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 z

21、hd inertial frame = j2000 ; b1950/j2000 (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

22、 ? sh_gamit -d 1998 240 -orbit siof -expt emed eops usno ? process gps data from day 240 of 1998 (aug 28) ? use any data in the local rinex directory精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 4 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 4 页，共 14 页 - - - - - - - - -? f

23、tp data for sites specified by the user? optionally scan multiple directories to find data? use final orbit estimate from scripps institution of oceanography (sio) 轨道估计数据来源? use eop data from us naval observatory (usno) iers bulletin a eop 数据来源? name the run “ emed工程名 identify files select data list

24、 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 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 op

25、tions? -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 if not the same as processing day ? -orbt 4 char id for orbit files. default 4 char extracted from the orbit selected ? -eo

26、ps 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 search rinex file names this many days into the future for current-day data 寻找往前推迟n天的 rinex 数据? -rx_doy_minus n search rinex fi

27、le names this many days into the past for current-day data 寻找往后 n天的 rinex 数据? -xsite list list of sites to be excluded from the processing. 解算排除的站? -netext add network suffix to day directories. i.e. 035r 日期文件夹添加网的后缀? -yrext ? add year prefix to day directories. ie 1997_053日期文件夹添加年前缀? -pres ? plot g

28、amit residuals as skyplots . yes, no, elev. default n 绘制天空残差图? -nogifs ? do not create skyplot gif graphics files不产生天空 gif图片精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 5 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 5 页，共 14 页 - - - - - - - - -copt file types compress rem

29、aining files for each dayserve when d default x k autcln.out 压缩指定类型的每天文件? -aopt file types copy the specified files for each day to the directory specified by archivepth 将指定每天文件复制到archivepth 文件夹。? -dopt file types delete specified files after processing each day. default c 删除指定类型文件 -mailto name of u

30、ser and machine 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 ch

31、ar 1-character 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 a

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

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

34、解算成功与否判断：postfit nrms: *e+00 是否在 0.25 左右，最差不能大于0.5，基线分量的误差大小是否小于预设值。至此基线解算全部结束。基线解算结果：sh_gamit_doy.summary 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 6 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 6 页，共 14 页 - - - - - - - - -基线解算结果：qexpta.doy,oexpta.doy,hexpta.yrdoy 10.

35、若软件限制了最大测站数据，最大卫星数目，最大对流层延迟参数，最大解算历元数目。更改于makefile.config文件。然后重新编。精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 7 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 7 页，共 14 页 - - - - - - - - -11.tables 表介绍。12.格式问题o 文件名需要小写，并且内容格式符合uinx 文件名改小写：sh_casefold -dir d -files * 格式转 uni

36、x: sh_dos2unix * 13.gamit 模块以及结构精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 8 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 8 页，共 14 页 - - - - - - - - -14.data types? l1 and l2 phase used for parameter estimation ? p1 and p2 range measurements used to help clean phase dat

37、a ? used to estimate receiver clock offset with tolerance of 1microsecond (= 300 m) ? lc phase proportional 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

38、 of l1, l2, p1 and p2 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 9 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 9 页，共 14 页 - - - - - - - - - melbourne-wubbena should be constant (+/- noise)15. aim of gamit basic aim of gamit is to generate positions/orbits for one day o

39、f 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 directory provides mapping from rinex header information

40、 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 igs rinex files ? default lines allowed, but be caref

41、ul17.模型说明：1.otl 潮汐改正2.vmfl gmf 投影函数3.atml 大气荷载模型，对高程影响较大，可消除周跳波动，可靠性需要进一步证实4.atl 大气抄袭荷载模型和met 气象模型下载地址： 1) 海洋潮汐。例如otl_fes2004.grid 放在软件 talbels 目录下。链接到otl.grid 。2) 大气负荷。例如atmldisp_cm.2006 ，每年更新一次。连接到atmldisp_grid.2006 3) vmfl 投影函数栅格，例如vmflgrd2006 ，连接到map.grid.2006。每年更新一次。18. 解算所

42、需文件说明类型文件名说明：数据文件bjfs0900.06o bjfs 为站点代码 ， 090 为年积日 ， 0为时段号 ， 06 为2006 年，o 为观测文件auto0900 06n 广播星历文件 ，前4 个字节改为站点代码则为站点广播星历文件igs13685 sp3 igs 为事后精密星历 ，1368 为gps 周，5为星期五igr13685 sp3 igr 代表快速精密星历igu13685 sp3 igu 代表预报精密星历配置文件lgigg6 090 l 代表 l 文件， gigg 代表工程名 ，6代表 2006 年， 090 为年积日station oct 测站信息文件 ，记录各测站接

43、收机和天线类型session info 时段信息文件 ，可以在数据处理过程中自动创建sestbl 解算参数控制文件，设置与解算有关的参数，达到不同的解算目的精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 10 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 10 页，共 14 页 - - - - - - - - -sittbl设置测站坐标约束、大气模型等参数表文件ut1 国际时间系统表pole 极移表leap sec 跳秒表nutabl 章动表luntal

44、 月亮表soltab 太阳表gdetic dat 地球形状参数表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 no

45、t given. - : output to stdout -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 operat

46、ion of all data are initialized 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)

47、 20某个人的参数设置：表2 解算控制参数的选择控制选项控制参数处理模式 relax 迭代方案 1-iter 观测量 lc_autcln 天顶延迟参数 25 number zen = 25 ; number of zenith-delay parameters (default 1) 梯度参数 3 -number gradients 海潮模型 otl.list 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 11 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - -

48、 第 11 页，共 14 页 - - - - - - - - -光压模型 berne 惯性框架 j2000 0 对流层误差模型 saas 截止高度角/( ) 10 - elevation cutoff 测站约束 ( bjfs) /m 0 05 /0 05 /0 1 测站约束 ( wuhn) /m 0 025 /0 01 /0 1 测站约束 ( kunm) /m 9 999 /9 999 /9 999 9 其他默认某篇论文写到 ; 不同截止高度角对基线精度的影响表4 不同截止高度角对基线(bjfs_kunm) 精度的影响精度指标高度角 /( )5 10 15 20 25 nrms 0 241 0

49、4 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.90 10 e-9 4 90 10 e-9 5 49 10 e-9 6 53 10 e-9 7 67 10 e-9 不过高度角太低时路径延迟对水汽变化非常敏感，因此在一般的解算情况下选择10的截止高度角比较适宜 。不同天顶延迟参数对基线精度

50、的影响表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 y /m 0 040 8 0 040 2 0 039 8 0 040 1 z /m 0 045 9 0 044 7 0 044 2 0 044 5 l /l 4 95 10 e-9 5 05 10 e-9 4 90 10 e-9 4 90 10 e-9 一般情况下 ，天顶延迟参数取为13，即每 2 h 产生一个天顶延迟估值是比较合适

51、的。 天顶延迟参数变化对长基线的精度影响不是非常明显。选取不同观测量对基线精度的影响表6 选取不同观测量对基线(bjfs_kunm) 精度的影响精度指标观测量lc_autcln l1l2 lc_help lc_only nrms 0 241 04 0 239 28 0 239 27 0 239 27 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

52、-9 6 73 10 e-9 lc_autcln( 观测值是无电离层的线性组合)不同测站坐标约束对基线(bjfs_kunm) 精度的影响精度指标测站约束 /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 1/1/1 精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 12 页，共 14 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 12 页，共 14 页 - - - - - - - - -nrms 0 239 38 0 241 18 0 241 07 0 240 64 0 240 59 x /m 0 013 3 0 037 2 0 047 4 0 119 4 0 223 9 y /m 0 019 9 0 034 8 0 040 1 0 081 1 0 146 2 z /m 0 015 3 0 036 0 0 046 1 0 119 9 0 225 9 l /l 3 66 10 e-9 4 41 10 e-9 4 75 10 e-9 8 02 10 e-9 9 11 10 e-9 测站约