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