2022年gamit学习总结解读

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

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

3、/opt/otlmap/otl_fes2004.grid otl.grid进行连接（ /opt/otlmap/otl_fes2004.grid为模型的路径）rm map.gridln -s /opt/otlmap/vmf1grd.2021_063 map.grid map.grid目前有 vmf 供应模型）ls -al 连接并查看为否已连接胜利；连接终止；（连接完成后安装更新时会自动连接）其它 otl 模型：otl_fes99.grid 182.7 mbotl_csr4.grid 45.6 mbotl_got00.grid 45.6 mbotl_nao99b.grid 87.1 mbatmos

4、phere 模型或许年更新：atml files185 mb/year3. sites.defaults 预备介绍sites.defaults 文件为站处理掌握文件；格式如下：all_sites expt xstinfobrus_gps expt ftprnx ftpraw localrx xstinfo xsite glrepu glreps glts其中ftprnx ftpraw localrx xstinfo xsite为 gamit 选项,glrepu glreps glts 为 globk选项,格式详细含义如下：expt 为 4 个字符的项目名称（工程名称）ftprnx表示该站点的r

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

6、est xstinfo（ test 为工程名）lhas_gps test localrx glrepu glreps精品学习资料精选学习资料 - - - 欢迎下载shao_gps test localrx glrepu4. 文件预备 文件可以采纳sh_upd_stnfo 命令生成也可以人工编辑得到；对于部分不规范的rinex 文件建议采纳人编辑的方法,并参照gamit 的 rcvant（天线与接收）文件进行编辑,确保后期处理的顺当；在更新前应将sites.defaults 文件中的all_sites expt xstinfo 注释掉,否就简单报

7、错,无法更新；采纳 sh_upd_stnfo 命令生成 文件方法：将 tables 中的 文件拷贝至rinex 中,在 rinex 文件夹下运行sh_upd_stnfo -l sd提取 igs 站点的信息, 得到.new然后见 .new 改名为 ,在运行sh_upd_stnfo -files *12o -expt test即可；采纳手工编辑的方法：将 tables 中的 文件拷贝至rinex 中,在 rinex 文件夹下运行sh_upd_stnf

8、o -l sd提取 igs 站点的信息； 会自动生成 .new文件；打开 .new删除无关条目,增加其他站点的信息；此时应打开rinex 文件进行核对,确保无误后将其储存为 ；应留意将天线高参考设为dharp （天线参考中心一般为天线底部）格式：fixed format file with station information. *trck sitestation nameant htant nant e . 7odm 7odmseven oaks dam 0.0083 0.0000 0.0000 .rcvrantco

9、dhtcodvers . ashz12 atdmrb dhpab9.10 . yeardoysnstart stop. 20001100 00 00 0024 00 00. = line break for this displaythis file must be sorted by trck、 then epoch trck and site4-character codematch the 4-char code in rinex file nametrck = site for static analysis在静态分析中. change for kinematic analysis对于

10、动态分析要转变file must be sorted by this field、 then epoch. station_nameis a long 12-char descriptionanthtraw antenna height as measured原始测量高度critical measurement.至关重要. htcoddefines geometry of antht measurement定义几何类型dhpab is rinex standard is vertical height to antenna reference point arp（ dhpab 指相对于站点的垂

11、直高度）. n and e offsets基准偏移e.g.、 uncalibrated tribrachrarely used精品学习资料精选学习资料 - - - 欢迎下载occupied point w.r.t. benchmarkrcvr6-character codes for the receiver typenot critical except for some old receiversor "buggy" firmware versions. versversion number of firmware in the receiver接收机固件版本号shou

12、ld be correct to handle snr coding properlyantcod 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 the values apply. times set part of dayif it changes duri

13、ng utc day. this repres entation will be augmented in future version of gamit5. 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 nnn 100. 100. 100.lha

14、z lhaz_gps nnn 100. 100. 100.shao shao_gps nnn 100. 100. 100.urum urum_gps nnn 0.050 0.050 0.05wuhn wuhn_gps nnn 0.050 0.050 0.05kunm kunm_gps nnn 0.050 0.050 0.05下载 sittbl.refined文件 、里面有最近的各站的牢靠性的一个评测；要想点位置精度高就设置igs 站的约束小一些、即认为igs 站已经很准 、强约束 、如目的为求基线后面平差就可以将igs 站的约束设置大些；6.process.defaults 文件可修改部分的说

15、明process.default 文件中需要修改的有采样间隔,最大历元数,开头时间等,示例如下：set sint = '30' 采样间隔依据需要和rinex 的情形可修改为15 30 60精品学习资料精选学习资料 - - - 欢迎下载set nepc = '2880' 最大历元数依据观测时长运算调整；set stime = '0 0' 开头时间（小时分）依据rinex 文件修改,一般要求时间与rinex 时间差不大于0.1h 否就报错；set use_rxc = "n"设置为 n 不依据 rinex 更新 lfile 文件；s

16、et aprf = itrf08.apr采纳 08 框架数据（默认itrf2000 ） 如不采纳globk平差可不设置；set minxf = '300' xfile文件大小的最小值,小于该值的将被抛弃；依据历元数设置；历元数少时应减小该值；set mailto = '*'设置为有效邮箱以便接受gamit 结算结果；7. sestbl.文件依据实际需要进行修改；示例如下：hoice of experiment = baseline（运算高精度基线）baseline/relax./orbit默认为 baseline ,但一般修改为relax.basline为固定轨

17、道）（如只需要基线文件可以挑选 baseline ）choice of observable = lc_autcln观测值类型挑选lc_autclndefault、 lc_helpcodeless l2、l1_onl y l1soln fromdual freq、l2_onl y l2soln from dual freq、l1、l2_independent l1 + l2 from dual freq、l1&l2 same as l1、l2_independent but with ion constraint;l1_receiver must add 'l1only'

18、; in autcln.cmd对于 15km 以下的短基线应挑选l1、l2_independent或 l1_onl ytype of analysis =0-iter迭代多少次 运算结果作为初始值再次运算一般设置1-iter zenith delay =n天顶推迟对于短基线选ninterval zen = 2;2 hrs = 13 knots/day default is 1 zd per day天顶推迟修改为1h,默认 2hupdate t/lfiles = l_only;t_and_ldefault、 t_onl y、 l_onl y、 none人工编辑时可设为none# dmap= gm

19、f; gmfdefault/nmfh/vmf1假如已下载vmf1这可设为vmf1-重新连接；# wmap = gmf; gmfdefault/nmfw/vmf1假如已下载vmf1这可设为vmf1 use map.grid = y;vmf1 grid file with mapping functions and zhdinertial frame = j2000; b1950/j2000 default = j20008. sh_gamit 批处懂得释> & . sh_gamit_test.log生成过程详解文件sh_gamit_test.log当采纳快速星历处理时需要设置-or

20、bit igsu/igsr#要处理最近的数据需要将gamit/gamit/fixdrv文件夹下,打开armake.f 文件,将2021 改为 2100.然后再 gamit 文件夹下运行 ./install_software. 即重新编译安装； （未验证）-mailto lyzlrz联网后发送-basic example. sh_gamit -d 1998 240 -orbit siof -expt emedeops usno. process gps data from day 240 of 1998 aug 28. use any data in the local rinex direct

21、ory精品学习资料精选学习资料 - - - 欢迎下载. ftp data for sites specified by the user. optionally scan multiple directories to find data. use final orbit estimate from scripps institution ofoceanography sio轨道估量数据来源. use eop data from us naval observatory usno iersbulletin aeop 数据来源. name the run“ emed" 工程名ident

22、ify filesselect datalist optionssh_gamit -help. select days of data to process. -d yr dayswhere yr and days are the data to be processed e.g.、 -d 1997 153 156 179. -s yr d1 d2where d1 is start day and d2 is stop day to be processed e.g. 1997 153 178. -r daywhere day is number of days before current

23、date.sh_gamit options. -expt4 char name of experiment/solution being run default expt. -orbittype of orbit to use. default siof. optional second argument gives the day-of-year for the orbit if not the same as processing day. -orbt4 char id for orbit files. default 4 char extracted from the orbit sel

24、ected. -eopsname of eop series to use. bull_b、 bull_a、 etc. default usno-sessinfosession - sampling interval、 number of epochs、 start time default 30 2880 0 0. -rx_doy_plus nsearch rinex file names this many days into the future for current-day data查找往前推迟n天的 rinex 数据. -rx_doy_minus nsearch rinex fil

25、e names this many days into the past for current-day data查找往后 n天的 rinex 数据. -xsite listlist of sites to be excluded from the processing.解算排除的站. -netext <char>add network suffix <char> to day directories. i.e. 035r日期文件夹添加网的后缀. -yrext. add year prefix to day directories. ie 1997_053日期文件夹添加

26、年前缀. -pres. plot gamit residuals as skyplots <y/n/e>. yes、 no、elev. default n绘制天空残差图. -nogifs. do not create skyplot gif graphics files不产生天空 gif 图片精品学习资料精选学习资料 - - - 欢迎下载-copt file typescompress remaining files for each dayserve when ddefault x k autcln.out压缩指定类型的每天文件. -aopt file typescopy the

27、 specified files for each day to the directoryspecified byarchivepth将指定每天文件复制到archivepth 文件夹；. -dopt file typesdelete specified files after processing each day. default c删除指定类型文件-mailtoname of user and machine to send results summary. e.g.、 -snrrescalculate theoretical phase resi

28、duals from snr. requires raw data or rinex with valid snr values. -noaccdo not exclude satellites base on accuracy codes in the .sp3 file. default exclude不排除 .sp3 中基于精度码方式排除卫星；. -xver char1-character x-file version 6th character of x-file. default is single-digit year. -noftpuse this option if no ft

29、p 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 available default c. n - don't remake x-files under any circumstances如tables 中 消逝, sh_gamit 将会使用 rinex 的

30、头文件信息；guess_rcvant.dat来将rinex 的头文件信息翻译成 gamit可阅读形式；antenna height errors在 globk 可降低,但类型识别错误无法排除；-9.打开日期目录查看结算结果 判定解算质量主要有两个方面：一为数据利用率（即为否有足够多的可用于进行合理估量的数据参加运算）； 二为数据的噪声水平为否与选用的模型匹配；（这一项主要有两个指标,基线重量误差的大小和nrms（标准化的均方根误差,对于基线重量误差的判定主要依 据预先设定的站点坐标和轨道参数的先验约束,假如大于先验可以通过查看q 文件或者autcln.sum 文件,为否有大量被 autcln

31、丢弃的数据；对于nrms 一般应在0.25 左右；假如大于0.5 就意味着周跳没有排除,或显现了模型错 误（固定站点坐标错误或者未模型化的卫星）等；解算结果文件autcln.sum.post 中的 rms 残差为显示掌握网站点的相对质量；打开sh_gamit_doy.summary , otesta.doy、 qtesta.doy 文件查看解算总结和基线解算情形,解算胜利与否判定：postfit nrms: *e+00为否在 0.25 左右,最差不能大于0.5,基线重量的误差大小为否小于预设值；至此基线解算全部终止；基线解算结果：sh_gamit_doy.summary精品学习资料精选学习资料

32、 - - - 欢迎下载基线解算结果：qexpta.doy、oexpta.doy、hexpta.yrdoy10.如软件限制了最大测站数据,最大卫星数目,最大对流层推迟参数,最大解算历元数目；更换于makefile.config文件；然后重新编；精品学习资料精选学习资料 - - - 欢迎下载11.tables 表介绍；12.格式问题o 文件名需要小写,并且内容格式符合uinx 文件名改小写：sh_casefold -dir d-files* 格式转 unix:sh_dos2unix*13. gamit 模块以及结构精品学习资料精选学习资料 - - - 欢迎下载14. data types. l1

33、and l2 phaseused for parameter estimation. p1 and p2 range measurementsused to help clean phase data. used to estimate receiver clock offset with tolerance of 1 microsecond = 300 m. lc phaseproportional to geometric rangesensitive insensitive to ionospheric delaylinear combination of l1 and l2. lg p

34、hasemeasures ionospheric delay onlyinsensitive to geomery. widelane wllinear combination of l1、 l2、 p1 and p2精品学习资料精选学习资料 - - - 欢迎下载melbourne-wubbenashould be constant +/- noise15. aim of gamitbasic aim of gamit is to generate positions/orbits for one day of data. standard product of gamit is a mini

35、mally constrained solution16. gg/tables/rcvant.datfile 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 fromrinex header information to 6-char codes. example linesrcv ashtechz-xii a

36、shz12rcv z-xii3 ashz12swv 5f 1.00. explanation given at top of file. used a represent space. sometimes need to be careful that pattern correctlygets correct type. should not be a problem for igs rinex files. default lines allowed、 but be careful17. 模型说明：1. otl潮汐改正2. vmfl gmf投影函数3. atml 大气荷载模型,对高程影响较

37、大,可排除周跳波动,牢靠性需要进一步证明4. atl 大气抄袭荷载模型和met 气象模型下载地址： 1) 海洋潮汐；例如otl_fes2004.grid放在软件 talbels 目录下；链接到otl.grid ；2) 大气负荷；例如atmldisp_cm.2006 ,每年更新一次；连接到atmldisp_grid.20063) vmfl 投影函数栅格,例如vmflgrd2006 ,连接到map.grid.2006 ；每年更新一次；18. 解算所需文件说明类型文件名说明：数据文件bjfs0900.06objfs为站点代码 , 090为年积日 , 0为时段号

38、 , 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测站信息文件 ,记录各测站接收机和天线类型session info时段信息文件 ,可以在数据处理过程中自动创建sestbl 解算参数掌握文件,设置与解算有

39、关的参数, 达到不同的解算目的精品学习资料精选学习资料 - - - 欢迎下载sittbl设置测站坐标约束.大气模型等参数表文件ut1 国际时间系统表pole 极移表leap sec跳秒表 nutabl 章动表luntal 月 亮表 soltab 太阳表gdetic dat地球外形参数表rcvant dat接收机及其天线型号对比表svnav dat星号对比表 svs_exclude dat坏卫星信息表 hi dat仪器高参数表dcb dat差分码偏差表autcln cmd 周跳自动探测和修改命令表antmod dat天线相位中心改正表ot1 list /grid全球格网或各测站的海洋潮汐参数表1

40、9 .crx2rnx压缩的 d文件转换成 o文件usage: crx2rnx file - -f -s -hstdin and stdout are used if input file name is not given.- : output to stdout-f : force overwrite of output file-s : skip strange epochs default:stop with errorthis option may be used for salvaging usable data when middle of the compact rinex fi

41、le is missing. the data after the missing part、 are、 however、 useless until the compression operation of all data are initialized at some epoch. combination with use of -e optionof rnx2crx ver.4.0 or after may be effective.caution : it is assumed that no change in # of data types happens in the lost

42、 part of the data.-h : display help messageexit code = 0 success= 1 error= 2 warning20某个人的参数设置： 表2 解算掌握参数的挑选掌握选项掌握参数处理模式 relax迭代方案 1-iter观测量 lc_autcln天顶推迟参数 25number zen = 25; number of zenith-delay parameters default 1梯度参数 3-number gradients海潮模型 otl.list精品学习资料精选学习资料 - - - 欢迎下载光压模型 berne惯性框架 j2000 0

43、对流层误差模型 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精 度的影响高度角 /°510152025nrms± 0 241 04± 0 241 04± 0 231 09± 0 222 52± 0 213 98x /m&#

44、177; 0 041 9± 0 041 9± 0 045 0± 0 048 5± 0 051 2y /m± 0 040 1± 0 040 1± 0 043 4± 0 049 1± 0 056 3z /m± 0 044 5± 0 044 5± 0 047 8± 0 051 1± 0 055 1l /l4.90 × 10 e-94 90× 10 e-95 49× 10 e-96 53 × 10 e-97 67×

45、 10 e-9不过高度角太低时路径推迟对水汽变化特别敏锐, 因此在一般的解算情形下挑选10° 的截止高度角比较相宜 ；精品学习资料精选学习资料 - - - 欢迎下载不同天顶推迟参数对基线精度的影响表5 不同天顶推迟参数对基线bjfs_kunm精度的影响精品学习资料精选学习资料 - - - 欢迎下载精度指标天顶推迟参数7131925nrms± 0 239 47± 0 242 77± 0 241 66± 0 241 04x /m± 0 042 5± 0 041 9± 0 041 8± 0 041 9y /m&

46、#177; 0 040 8± 0 040 2± 0 039 8± 0 040 1z /m± 0 045 9± 0 044 7± 0 044 2± 0 044 5l /l4 95× 10 e-9 5 05× 10 e-94 90× 10 e-94 90× 10 e-9一般情形下 , 天顶推迟参数取为13, 即每 2 h 产生一个天顶推迟估值为比较合适的； 天顶推迟参数变化对长基线的精度影响不为特别明显；精品学习资料精选学习资料 - - - 欢迎下载选取不同观测量对基线精度的影响表6 选取

47、不同观测量对基线bjfs_kunm精度的影响精品学习资料精选学习资料 - - - 欢迎下载精度指标观测量nrmslc_autcln± 0 241 04l1±l20 239 28lc_help± 0 239 27lc_only± 0 239 27x /m± 0 041 9±0 047 3± 0 047 9± 0 047 9y /m± 0 040 1±0 044 1± 0 044 6± 0 044 6z /m± 0 044 5±0 053 8± 0

48、055 3± 0 055 3l /l4 90× 10 e-96 53× 10 e-96 73× 10 e-96 73× 10 e-9lc_autcln观测值为无电离层的线性组合精品学习资料精选学习资料 - - - 欢迎下载不同测站坐标约束对基线bjfs_kunm精度的影响精品学习资料精选学习资料 - - - 欢迎下载精度指标测站约束 /m0.01/0.01/0.020.05/0.05/0.060.1/0.1/0.10.5/0.5/0.51/1/1精品学习资料精选学习资料 - - - 欢迎下载nrms± 0 239 38± 0 241 18± 0 241 07± 0 240 64± 0 240 59x /m± 0 013 3± 0 037 2± 0 047 4± 0 119 4± 0 223 9y /m± 0 019 9± 0 034 8± 0 040 1±