|||
这里要用的linux的程序,包括cat、grep、awk、sort、sed
在windows下,可以先安装cygwin,然后安装上面需要的程序.下面的命令是cygwin下执行的的。
1、合并多天gamit的o文件
cat *o????a.* >step1.txt
2、提取基线长信息
grep 'X X' step1.txt|gawk '{print $1,$16,$18}' >step2.txt
3、根据基线名排序
sort step2.txt >step3.txt
4、把空格替换为逗号,方便计算重复率程序处理
sed 's/ /,/'g step3.txt > step4.txt
最终结果如下:
A033_GMSD,1176008.1036,0.0028
A033_GMSD,1176008.1040,0.0031
A033_GMSD,1176008.1041,0.0030
A033_GMSD,1176008.1044,0.0028
A033_GMSD,1176008.1046,0.0029
A033_GMSD,1176008.1046,0.0030
A033_GMSD,1176008.1062,0.0032
A033_GMSD,1176008.1079,0.0034
A033_KUNM,1695123.7931,0.0040
A033_KUNM,1695123.7944,0.0041
A033_KUNM,1695123.7947,0.0039
A033_KUNM,1695123.7953,0.0039
A033_KUNM,1695123.7953,0.0041
A033_KUNM,1695123.7961,0.0040
第一列为基线名(相同的基线列在一起),
第二列为基线长(单位m),
第三列为基线长误差(单位m)
经基线重复率计算后得
A033_GMSD, 1176008.10479296, 0.00123866390486845,1.05327837437525E-09
A033_KUNM, 1695123.79590439, 0.00204319966805489,1.20533949968226E-09
A033_SHAO, 480733.898299702, 0.00053035414908655,1.10321770726456E-09
第一列为基线名,
第二列为加权基线长(单位m),
第三列为绝对基线重复性(单位m),
第四列为相对基线重复性(单位m)
GLOBK 自带的 sh_globk_scatter 得到相应的基线向量的重复性
sh_globk_scatter -f globk_rep.prt
multibase val.globk_rep.prt -s baseline –d
baseline是基线站点名单文件,详细参考globk手册
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-16 16:44
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社