|||
1、globk_comb.cmd
* This group must be first
eq_file ../tables/eq_renames
make_svs ../tables/sat1.apr
com_file comb.com
srt_file comb.srt
sol_file comb.sol
* Earth-rotation values (not needed if pmu free in final combination)
in_pmu ../tables/pmu.bull_b
* File(s) for a priori coordinates
apr_file ../tables/itrf00.apr
* Sites to include in the solution (default is all)
use_site clear
use_site mate ankr_gps kit3 nico sofi zwen kant
use_site –zwen
* (1) Max chi**2, (2) Max prefit diff, (3) Max rotation; defaults are 100 10000 10000
max_chi 30 50 2000.0
* Apply the pole tide whenever not applied in GAMIT
app_ptid ALL
* Invoke glorg for stabilization
org_cmd glorg_comb.cmd
org_opt PSUM GDLF CMDS
* Set minimal globk print options since using glorg output
prt_opt NOPR
* Allow the network to be loose since using glorg for stabilization
apr_neu all 10 10 10 0 0 0
* Satellites are loose if using global data -- ignored if GAMIT BASELINE mode
x apr_svs all 100 100 100 10 10 10 1R
* but tight if not combining with global data
x apr_svs all .05 .05 .05 .005 .005 .005 .01 .01 FR
* EOP loose if estimating rotation in glorg
apr_wob 100 100 10 10
apr_ut1 100 10
* Write out an h-file if needed for future combinations
x out_glb H------.GLX
上面第一组命令定义了glred/globk中使用的文件,这些命令必须放在控制文件的前面,这样程序才能很好的运行。对于时间序列计算指令,除了 eq_file所有的都可以认为是临时文件。下面将要介绍的eq_file使得我们可以排除测站或重命名测站以解释由于地震或仪器变化造成的地形突变。
紧接着的两个记录具体指定了用于地球旋转和坐标的先验文件。如果在此忽略的话,将会采用h文件中的,这样在很多情况也是可行的。
use_site命令允许我们选择解算中包含的测站。这里我们展示了多余的命令以举例说明它们的用法:默认的是all;use_site clear命令将删除所有测站,在它的后面应该跟着你希望解算中包含的测站清单;在列出的测站名称前加负号,将会在清单中删除它。对于从全球处理的h文件 中选择测站,我们可以执行use_pos以在一定的地理区域内使用或排除测站(查看globk的帮助文件或者globk参考手册的3.1节)
当数据是“坏”的时候 max_chi命令允许自动删除h文件,在globk中检查高卡方值探测,坐标相对于先验值调整较大,或者测站网较大的旋转。
包含glorg命令文件的名字会告诉globk运行完globk后执行glorg for stabilization(框架定义)。org_opt选项规定了glorg打印总结中包含的内容(org文件)。PSUM表示位置平差的总 结;GDLF表示列出输入的h文件的名字;CMDS表示显示globk 和glorg命令文件。如果你正在执行glorg,然后globk打印将不起作用,且通过prt_opt NOPR压缩。
apr_形式的命令告诉globk估计参数和指定使用的先验值。如果你通过在glorg中估计旋转和转换固定框架(通常是稳健网的情况下),然后坐标 (apr_neu)和EOP(apr_ut1 apr_wob)约束应该是松弛的,如上所示。如果只在glorg中估计转换,在小区域中可选的是少于4个站,然后应该严密约束EOP。为了应用有限的约 束,要为每个站设置合适的apr_neu,且不能执行glorg。
out_glb命令将会在执行时告诉globk将解算输出以h文件的形式输出。这对于聚集按周或按年的平均值和将自己的处理与全球或局部h文件联合解算很有用。
2、glorg_comb.cmd
* Stabilization using IGS stations
apr_file ../tables/itrf00.apr
* Regional stabilization using coordinates from a prior solution
x apr_file ../tables/globk_vel_990306.apr
* Define the stabilization frame
stab_site clear
stab_site mate kit3 nico sofi zwen
* Set parameters to estimate in stabilization
pos_org xtran ytran ztran xrot yrot zrot
* Default height ratio is 10; increase if heights questionable
cnd_hgtv 1000 1000 2.0 10.0
* Default iterations is 4, relative weight ratio 50%, horizontal sigma cuoff 4.0
x stab_ite 4 0.8 3.
定义框架通过在glorg中具体指定一系列先验坐标已知且满足框架定义要求的测站来实现(通常水平式1-2mm,垂直2-5mm),想要估计的参数是为了 使这些测站的坐标调整量最小。如果在glorg控制文件中包含apr文件,它的坐标将会取代gamit或globk中使用的坐标,所以对于三个程序没有必 要使用相同的先验坐标。对于glorg,着重关心的是稳定测站的坐标。
pos_org命令指定了估计的参数,只进行平移(忽略后三个记录) or 转换和旋转.stab_site的使用规则同globk控制文件中的use_site一样。
最后两个记录控制固定过程中的(in the stabilization)准观测值(测站坐标)的权。对于cnd_hgtv,前两个值具体指定相对于水平坐标的downweighting of height(两个值分别对应位置和速度)。默认的是10,假设高度的不确定性事水平坐标的三倍;但是如果你想定义一个水平框架或知道高度估计不可靠,你 可以提高这个比率。最后两个值限制高度不确定性,允许包含一个测站和删除不稳定的测站????。对于stab_ite,第一个值是 stabilization过程中迭代的次数,第二个值是迭代过程中重新加权的数量,第三个是在删除测站前允许的相对于不确定性的残差比率。你可以在定义 框架的过程中改变最后一个值以调整包含或排除的测站。这些命令的详细描述可看glorg的帮助文件和globk参考手册的3.2节。
在固定过程中允许平移和旋转是absorb空间相关噪声最普遍的方法,但是如果你只有少量的测站或异常的几何体时要小心。例如,如果你只有三个站,包含旋 转,你从九个观测值(每个站三个坐标)中估计六个参数,将会导致只有三度的冗余度。而且即使有很多测站,如果其中一个离其他的站很远的话,长基线会使 absorb所有的旋转自由度。如果测站可靠,不关心测站残差的话这并不是一个致命的缺点,但是你需要知道this chararistic of the stabilization.
如果你需要提供坐标或速度的step changes,或在某个时期排除测站,可以在globk中使用地震文件(eq_file)。这个文件有两种记录。第一个重命名测站(为globk创建新的逻辑上的测站)基于它们接近地震的程度:
# November 3, 2002, UT=22:12 Lat = 63.52 Lon -147.53 Depth=5 km M 7.9
eq_def DN 63.52 -147.53 700 5. 2002 11 3 22 12
eq_cosei DN 1.0 1.0 1.0 1.8 1.8 0.7
eq_post DN 100 0.1 0.1 0.1 1.8 1.8 0.7 Postseismic: 0.1 mm**2/day and spatial about same
eq_rename DN
eq_file文件中包含以上各行,开始和结束时间晚于2002.11.3的22:12的距离Denali地震700公里内的测站将会被重命名。 eq_cosei和eq_post命令描述了co-seismic和post-seismic变化应用的限制(查看globk帮助文件和globk参考手 册的3.1节)。第二种类型重命名测站的命令按照具体的说明解释数据质量和仪器的变化:
rename mate_gps mate_1ps 1999 6 18 0 0
rename kit3_gps kit3_xcl 1997 5 14 0 0 1997 5 18 0 0
rename madr_gps madr_xhi emed 1996 9 26 0 0 1996 10 27 0 0
第一行说明在Matera使用新接收机导致的位置变化;第二行是删除Kitab四天中的坏数据;第三行说明在h文件中名字包含字节emed的使用了不正确 的天线高度。可以使用glorg中的equate特征强制madr_gps和madr_xhi使用相同的水平调整(查看第四章)。也可以使用rename 命令改变测站位置而不用引入一个新的测站,应用已知的当地测量的偏差或地震模型中假定的偏移:
rename tigr_gni tigr_gni 2001 2 28 0 0 2100 1 1 0 0 0.0032 0.0033 0.0000 NEU
该记录使用了测站坐标北和东方向基于先验模型3mm的偏差。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-22 19:49
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社