PurpleSky的个人博客分享 http://blog.sciencenet.cn/u/PurpleSky

博文

使用GLOBK批处理命令进行网平差

已有 8624 次阅读 2015-4-3 09:52 |个人分类:GAMIT/GLOBK|系统分类:科研笔记| GLOBK, 网平差

1. 说明

首先,你要先完成基线解算。

sh_glred 在工程级目录执行但在解目录gsoln中运行,该目录也可以由 process.defaults 文件中 glbpth 参数指定。

 

2. 参数配置

2.1 globk_comb.cmd glorg_comb.cmd

若这两个文件不存在,sh_glred 尽管可以生成他们,但推荐手动创建

运行命令“sh_glred -cmd”,将在gsoln文件夹下生成上述配置的模板文件以供修改。

需在其中配置先验坐标文件(.apr)、解算中用到测站、定义参考框架时用到的测站和约束等(若使用IGS站点作为固定点,则可以使用自带的itrf 文件作为先验坐标)。

2.2 sites.dafaults

每步处理中所使用的测站都由sites.dafaults 中标记指定:在 glred 组合中包含的测站由 glrepu 标记,在glorg 中用来定义参考框架的测站 glreps 标记,要生成时间序列图的测站由 glts 标记。默认使用所有测站并绘图。

注意:若在sites.defaults文件中使用glreps标记平差时的固定点,则在glorg_combNaNd中需将“stab_site clear”选项注释掉(前边加“x”)。

 

3. 命令执行

完成上述操作之后,你就可以进行批处理网平差了。在工程目录执行:

sh_glred -s<yr1 doy1 yr2doy2>  -expt <expt>  -net <networks>  -local  -nettext <char>  -yrext<year> -ncomb <num>  -stnfo<station.info>  -cmd  -opt <A F H L U G E K C R>

其中-s 指定开始时间(yr1 doy1)和结束时间(yr2 doy2),expt 4 个字符的工程名,network 表示包含SOPAC H-文件,SOPAC默认选项是all所有),或者仅组合某些解假如-local 被指定了,脚本会仅仅处理可以获得数据范围内的天;否则它会自动去 SOPAC 查找 H-文件,处理所有指定天数内的数据。要处理的天也可用 -d yr days或者采用-r days 详细指定。可选参数-ncomb用来指定每一个组合包含多少天默认为 1 天,但也可以用脚本产生每周或每月的本地或全球文件的平均值

下面是opt选项可选参数的详细释义(摘自GAMIT 10.35中文文档)

  • <R> 在开始之前把旧的 h*.gl?文件从 glfpth 目录下删除掉

  • <F> SOPAC 网站上下载全球 h 文件,根据<networks>参数项提供的名称搜索文件。

  • <H> 对所有现有的或者链接到 glfpth(通常是 procdir/glbf) ascii 文件运行 htoglb

  • <LA> 链接所有本地存档的 ascii H 文件(即所有的 h[net|expt]?.yyddd)到组合中去。

  • <LB> 链接所有本地存档的 ascii H 文件(即所有的 h[net|expt]?.yyddd)到组合中去。从 htnd 开始向下搜索,其中 htnd 是系统默认的。

  • <LC> 链接本地生成的组合二进制 H 文件(名为*.GLX)

  • <U> 对二进制 h 文件运行 hfupd(尚未试)

  • <G> 运行 glred 以进行组合或重新组合

  • <E> 运行 ensum sh_baseline 进行绘图

  • <C> 后压缩 ascii H 文件,删除所有的链接,分别复制 ascii H 文件和二进制 H 文件为 glbpth/ascii_yyyy     glbpth/bin_yyyy, 其中 glbpth 是由 process.defaults 指定的。  

  • <A> 完成所有选项


例:  % sh_glred -s $year1 \$sdoy \$year2 \$edoy -expt $expt  -local -opt H G E

$year1$sdoy为要进行网平差的开始年与年积日,$year2$edoy为结束年与年积日。

 

4. 解算结果

命令运行完成后,在gsoln/文件夹下即可找到平差结果文件*.prt和*.org。

固定点有“*”号标记。

PS:虽然这种方法可以快速得到网平差结果和站速度,但是在GAMIT/GLOBK软件的官方指南中对这种平差方法是不推荐的。推荐方法是只使用该方法合并GPS网得到站坐标,然后使用globk命令调用globk_vel.cmdglorg_vel.cmd得到站速度(因整理时间较晚,有些可能记得不对,请读者自行查阅 GLOBK_Ref)。



https://blog.sciencenet.cn/blog-1365943-879492.html

上一篇:使用GAMIT进行批处理基线解算
下一篇:使用 install_updates 升级 GAMIT/GLOBK
收藏 IP: 119.167.70.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-12-23 14:10

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部