linux批处理访问网站,使用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  -expt   -net   -local  -nettext   -yrext -ncomb   -stnfo  -cmd  -opt

其中-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中文文档)

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

• 从 SOPAC 网站上下载全球 h 文件，根据参数项提供的名称搜索文件。

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

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

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

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

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

• 运行 glred 以进行组合或重新组合。

• 运行 ensum 和 sh_baseline 进行绘图。

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

• 完成所有选项

例：  % 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)。

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