chaoshu的个人博客分享 http://blog.sciencenet.cn/u/chaoshu

博文

[转载]利用gamit结果做基线重复率----预处理命令

已有 3553 次阅读 2015-7-6 12:56 |个人分类:GAMIT/GLOBK|系统分类:科研笔记|关键词:学者|文章来源:转载

这里要用的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手册




https://m.sciencenet.cn/blog-858128-903245.html

上一篇:GAMIT学习中遇到的问题
下一篇:NRMS and WRMS generated by gamit/globk

1 陆泽橼

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-18 15:55

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部