Gamit 10.70 示例运行

2019-1-7 3:11:24 来源: GNSSer 发布人:

如果您是有经验的用户,您可能只想执行批处理文件'runtest'在大约90分钟内完成所有步骤而无需干预。事实上,2018年只运行一天或一个全球导航卫星系统确认安装正确。但是,如果你是新手,我们强烈建议您执行每个命令手动并仔细阅读结果,直到您了解它为止一步完成了。如安装README中所述,您将需要下载并链接otl_FES2004.grid或编辑sestbl。省略海洋潮汐负荷。
 
该示例设置为

1)从CDDIS下载RINEX 3文件2016年欧洲和周边地区的10个连续站点,为期3天,2017年和2018年(请注意,如果RINEX 3使用长名称文件,它们将被重命名为旧的短名称下载);

2)仅对GPS进行阶段处理(sh_gamit)或GPS,Glonass,Beidou和Galileo;

3)每天计算和绘图每年的重复性(sh_glred)并结合每日将h文件(.glx)转换为每年的单个H文件(.GLX);和

4)计算三年的可重复性和速度在一起(globk / glorg)。唯一提供的输入文件例如sites.defaults,它始终特定于网络你正在处理。所有其他人都被复制或链接gg / tables,适用于大多数网络。

该示例建立的结构具有三个GAMIT处理目录,按年命名(任意)(/ 2014,/ 2016,/ 2018),每个都有一个/ rinex,/ tables和/ gsoln目录具体到那一年。在顶层有一个处理(/ vsoln)目录和多年GLOBK解决方案的/ tables目录。
所描述的步骤假设您已经下载并链接了海洋潮汐装载网格(例如otl_FES2004.grid,链接到otl.grid在gg / tables中)并在处理时访问互联网;如果你这样做没有这些,见下面的注释3和4。在目录/ check_files中是每天保存的q文件和.org文件的副本,和多年可重复性的.org文件和postscript文件和速度。
                            
在开始之前,请确保已构建路径和安装README中描述的别名。可以运行该示例从您系统上的任何目录,最好是您想要的位置处理你自己的数据,而不是gg。

第1步:从2018年开始运行GAMIT三天。

在示例/ 2018目录类型中
  

  sh_setup -yr 2018

                                 
2018 / tables目录现在包含大多数标准的链接gg / tables中的文件和process.defaults的这些文件的副本,sestbl。,station.info(完整的MIT版本)和autcln.cmd。该sites.defaults文件已经在目录中,因此没有被sh_setup覆盖。

检查sites.defaults以注意它已设置为ftp远程存档(默认为CDDIS)10个IGS站的RINEX文件并且'xstnfo'设置为避免自动更新station.info
在处理过程中。请注意,sittbl。从gg / tables复制建立以对IGS核心站施加适度限制以支持GAMIT中的模糊解决方案;示例中包含的四个是
绰绰有余。

编辑process.defaults将“mailto”更改为您自己的电子邮件地址接收sh_gamit摘要文件。 (如果保留为null,则恢复为'whoami',所以可能不需要这种改变。)注意'aprf',使用初始化lfile。对于GAMIT指向itrf14_comb.apr,其中已被sh_setup复制到/ 2018 / tables目录。

如果您还没有从早期的发行版下载或复制过海洋潮汐装载的(大)网格文件,不想使用它,在sestb中更改'Tides applied'。从'31'到'23'。

使用。构造一个小的,特定于实验的station.info文件以下2018 / tables目录中的过程:
                                            

sh_upd_stnfo -l sd


将使用MIT station.info创建'station.info.new'sites.defaults中列出的网站。 (此步骤将需要一段时间MIT全局station.info文件太长了。)检查后,重命名它到'station.info'(覆盖不再有用的MIT station.info)。

 mv ./station.info.new ./station.info

(在您自己的处理中,如果您有不在麻省理工学院的电台station.info文件,您可以从RINEX标头添加它们的条目通过在/ tables中运行'sh_upd_stnfo -files ../rinex/*.18o'。 (见章节如果RINEX标题条目是什么,请参阅Intro手册的2.4非标准)。

                                     

GAMIT目前支持每天分别处理每个GNSS每个目录创建时添加一个单字母来指定系统,例如095G,095R,095E,095C。如果您希望测试完整的GNSS

功能,键入/ 2018级别

                 

 sh_gamit -expt eura -gnss G -s 2018 095 097 -pres ELEV -orbit igsf -copt x k p -dopts c ao>&! sh_gamit_2018G.log
 sh_gamit -expt eura -gnss R -s 2018 095 097 -jclock sp3 -pres ELEV -orbit codm -copt x k p -dopts c ao>&! sh_gamit_2018R.log
 sh_gamit -expt eura -gnss E -s 2018 095 097 -pres ELEV -orbit codm -copt x k p -dopts c ao>&! sh_gamit_2018E.log
 sh_gamit -expt eura -gnss C -s 2018 095 097 -pres ELEV -orbit codm -copt x k p -dopts c ao>&! sh_gamit_2018C.log


(使用bash使用'> sh_gamit.log&'进行重定向。)

                                     

对于仅限GPS,请省略-gnss,以便日目录名称仅为年日。

                         

每天完成执行时,应通过电子邮件发送摘要文件。检查这些文件的站数(6),Postfit RMS(3-15 mm,无0),postfit nrms(~0.2),模糊度分辨率和坐标调整(<30 cm)。对于这个大型网络,2018年的模糊度分辨率为80-90%伽利略,北斗为40-50%,但由于轨道不良,格兰纳斯接近零,条件也反映在/ figs目录中的天空图中。允许Glonass轨道调整(sestbl中的RELAX模式)将改善非常合适。查看GAMIT轨道模型拟合的结果输入轨道,请参阅/ igs中的rms文件。查看完整的编辑统计信息在日期目录中查看autcln.prefit.sum和autcln.post.sum;为了最小二乘调整,见qeuraa。[ddd]。


(可选)从日期目录中删除x,k和p文件以节省空间:

 sh_cleanup -d 2018 095 096 097 -dopts x k p


注意:

                                           

1.对于大型或复杂的数据集,实用程序sh_get_times可能会有所帮助确定要处理的天数和会话跨度。


2.在为自己的实验创建station.info时,检查是很重要的它从RINEX标头更新后,除非你确定这些标头是正确。在处理中,station.info条目总是覆盖任何东西

在RINEX标头或x文件中。另一种创建条目的方法对于调查模式站点,使用交互式程序'make_stnfo',然后使用sh_up_stnfo将此文件与从MIT station.info创建的文件合并连续网站的文件。调查模式文件的格式较短station.info格式,但在合并时将被转换连续文件,应该在其中列为引用(-ref)调用sh_upd_stnfo。


这个例子是为了使用海洋潮汐负荷('使用otl.grid = Y')sestbl。),它要求你先前下载到gg / tables中来自everest.mit.edu上的匿名ftp目录的OTL文件并已链接这个文件到otl.grid也在gg / tables中。 IERS / IGS标准模型是otl_FES2004.grid(730 Mb)。但是,您可以替换较小的(45 Mb)otl_CSR4.grid,或者您可以通过设置关闭海洋潮汐负荷'tides applied = 23'和'使用otl.grid = N'在sestbl中。 (注意链接到其他网格和列表文件(met.grid,met.list,map.grid等)

可以保持为空以运行示例和大多数处理。

          

4.如果您希望或需要在没有Internet访问的情况下运行该示例运行时,您可以将RINEX,导航和轨道文件预加载到 / rinex,/ brdc和/ igs目录。


                                                               

第2步:运行GLOBK以获得可重复性和跨度的组合H文件


在/ 2018年:


  sh_glred -cmd

获取从gg / tables复制到2018 / gsoln的globk.cmd和glorg.cmd文件,

然后

  sh_glred -s 2018 095 2018 097 -expt eura -gnss G R E C -opt R H G T>&! sh_glred.log

                                     

命令脚本将每天目录中的GAMIT ascii h文件转换为GLOBK二进制文件(.glx)并将它们放入/ glbf目录(H选项);为列出h文件的每一天创建一个gdl文件,使用globk.cmd和glorg.cmd每天运行GLOBK(G选项);并生成时间序列图

(T选项使用pos文件,程序tssum和sh_plot_pos)。


对于此测试,由'sh_glred -cmd'在/ gsoln中创建的globk / glred和glorg命令文件无需编辑即可运行;但是,一般情况下,您可能需要进行更改,例如,仅使用平移而不是转换和旋转来定义参考框架,以及更改列表

参考(稳定)站点。


可以在gsoln / plots_2018_095-2018_097中找到每天结合GNSS测量结果的每日图。您可以查看* .org文件以查看globkdid,注意稳定迭代(删除了哪些网站),'实验列表'中的chi2增量,以及

稳定('POS STATISTICS')。在您自己的运行处理许多天,通过grep'ing * .org与'USED'(来自chi2列表)和'POS STAT'获得快速摘要是很方便的。


接下来结合几天来获得3天跨度的单个H文件,用于多年可重复性(时间序列)和速度。

 sh_glred -s 2018 095 2018 097 -expt eura -netext G R E C -ncomb 3 -globk_cmd_prefix COMB -opt G>&! sh_glred_comb.log

                                                                           


第3步:重复2016年和2014年的步骤1和2:


  在/ 2016:

sh_setup -yr 2016

  在/ 2016 / tables,编辑mailto的process.defaults

  在/ 2016 / tables:

sh_upd_stnfo -l sd; mv station.info.new station.info

  在/ 2016:

    sh_gamit -expt eura -gnss G -s 2016 101 102 -pres ELEV -orbit igsf -copt x k p -dopts c ao>&! sh_gamit_2016G.log
   sh_gamit -expt eura -gnss R -s 2016 101 102 -jclock sp3 -pres ELEV -orbit codm -copt x k p -dopts c ao>&! sh_gamit_2016R.log
   sh_gamit -expt eura -gnss E -s 2016 101 102 -pres ELEV -orbit codm -copt x k p -dopts c ao>&! sh_gamit_2016E.log
   sh_gamit -expt eura -gnss C -s 2016 101 102 -pres ELEV -orbit codm -copt x k p -dopts c ao>&! sh_gamit_2016C.log
   sh_glred -cmd
   sh_glred -s 2016 101 2016 102 -expt eura -opt H G T>&! sh_glred_2016.log
   sh_glred -s 2016 101 2016 102 -expt eura -netext G R E C -ncomb 2 -globk_cmd_prefix COMB -opt G>&! sh_glred_comb.log

   

  在/ 2014:sh_setup -yr 2014

  在/ 2014 / tables:编辑mailto的process.defaults

  在/ 2014 / tables:sh_upd_stnfo -l sd; mv station.info.new station.info

  在/ 2014

   sh_gamit -expt eura -gnss G -s 2014 113 114 -pres ELEV -orbit igsf -copt x k p -dopts c ao>&! sh_gamit_2014G.log
   sh_gamit -expt e -gnss R -s 2014 113 114 -jclock sp3 -pres ELEV -orbit comf -copt x k p -dopts c ao>&! sh_gamit_2014R.log
   sh_glred -cmd
   sh_glred -s 2014 113 2014 114 -expt scal -opt H G T>&! sh_glred_2014.log
   sh_glred -s 2014 113 2014 114 -expt eura -netext G R -ncomb 2 -globk_cmd_prefix COMB -opt G>&! sh_glred_comb.log

                                              

我们在2016年和2014年仅指定了2天以节省运行时间,并且我们仅包括GPS,因为北斗和伽利略的卫星太少,所以2014年和格洛纳斯一样那个时候获得准确的解决方案。 


步骤4:运行GLOBK以获得3纪元(4年)的可重复性和速度


此步骤的关键用户特定控件(也包含在每年内运行的sh_glred中)是用于定义参考框架的站点列表,以及这些站点的先验坐标,这里我们使用相同的站点

并且在单年解决方案中协调文件(itrf14_comb.apr),但情况可能并非总是如此。对于多年可重复性,我们将使用在year / gsoln目录中创建的组合H文件(.GLX)。因此,在顶层创建一个解决方案目录

  mkdir vsoln


然后cd vsoln并为globk创建一个组合h-fiel的列表:

 ls ../????/gsoln/H*GLX>! eura.gdl


对于大型或复杂的数据集,此时运行glist将有助于检查错误并为您提供所有使用的站点及其跨度的列表。然后,程序glist2cmd可以帮助建立use_site列表。它也可能是理想的

汇总每年的数据,既可以节省多年组合的时间,又可以提供更具代表性的长期统计数据。执行此操作的过程在ftp区域或网站的文档目录中的文件sGPS_recipe.txt中给出。

   

 \ rm globk_replong.org globk_replong.log
 glred 6 globk_replong.prt globk_replong.log eura.gdl globk.cmd>&! glred.out

                                        

与sh_glred不同,这里可能需要rm命令,因为glred会将新的日志和org文件与之前的任何t文件连接起来

创建用于绘图的pos文件:

tssum。 mit.final.itrf14 -R globk_replong.org


创建多年图:

sh_plot_pos -f * .pos -r -t NONE -u -t1 2014-001 -t2 2018-180


获得4年的速度:


 \ rm globk_vel.org globk_vel.log
 globk 6 globk_vel.prt globk_vel.log eura.gdl globk.cmd VEL>&! globk.out

                                                             

最后的'VEL'标记告诉globk取消注释

以'VEL'开头的globk_long.cmd。


检查实验列表下的globk_vel.org文件中的chi2增量(<1.0)堆叠3个H文件和稳定统计数据'VEL STATISTICS'和'POS STATISTICS'系列。对于这个网络而言ITRF 2014 apr文件,9个稳定站适合参考框架速度在~1毫米/年水平和5毫米/年垂直水平大约一半有效值的不确定性。


绘制速度使用

 sh_plotvel -ps euratest -f globk_vel.org -R-60/50 / -5 / 55 -factor 0.5 -arrow_value 10 -page L

     

(在没有参数的情况下输入'sh_plotvel'可以查看完整的选项范围制作更具吸引力的地图。)


------------


从头开始,每年运行一次sh_clean_test和一次用于速度解决方案。 (不输入脚本的名称查看文档的参数。)




阅读次数: 3563

下一篇: Gamit 海潮文件更新
上一篇: Gamit 10.70 安装笔记,基于Ubuntu16.04 /12.04/18.10

尚无评论!

返回上一页面