感谢支持
我们一直在努力

使用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.cmd和glorg_vel.cmd得到站速度(因整理时间较晚,有些可能记得不对,请读者自行查阅 GLOBK_Ref)。

GAMIT相关阅读:http://www.linuxidc.com/GAMIT/

GAMIT10.4及最新版本的GAMIT下载在Linux公社5号FTP服务器,具体下载见 http://www.linuxidc.net/thread-1186-1-1.html

本文永久更新链接地址:http://www.linuxidc.com/Linux/2016-02/128279.htm

赞(0) 打赏
转载请注明出处:服务器评测 » 使用GLOBK批处理命令进行网平差
分享到: 更多 (0)

听说打赏我的人,都进福布斯排行榜啦!

支付宝扫一扫打赏

微信扫一扫打赏