|||
处理过程中,为了自己的方便使用,需要编写自己的脚本(用红色标出)
求解基线估计也可以通过globk实现
1、在prt_opt中输入这两个命令
· BLEN -- Baselinelengths
· BRAT -- baselinerates when velocities estimated
2、然后编辑完globkcomd的cmd后,直接运行
glred6 glred_20150811.prt glred_20150811.log 20.gdl cmd
3、计算基线的wrms与nrms解
sh_globk_scatter -f sina.prt
4、提取基线长度以及生成mb文件
multibase val.sina.prt
5、生成GMT的基线重复性验证解通过前面提取得到:
sh_baseline -f mb* -cols 1 -y -30 30
当然在gamit中也可以得到基线重复性的解:
主要利用o文件,o*a文件,然后进行求解
可以选择制定的基线解:然后求解
sh_gamit_baseline -o o* -b PIN2_VNDP PIN2_YAM2"
另外:可以利用其它的工具进行求解基线重复性
gamitcookbook中有这样的脚本
#!/bin/csh
# generate repeatability statistics
# collect the O-files:
# $1 is name of run
set files = ../3*/o$1?.???
#grep Normalized $files
# collect the baselines in NEU
if (-e tmp.bsl) then
/bin/rm tmp.bsl
endif
# Get the bias free solution
grep -h _ $files | grep 'R N' | sort >! $1\.r.neu.bsl
# Get the bias free solution XYZ
grep -h _ $files | grep 'R X' | sort >! $1\.r.xyz.bsl
# Get the bias fixed solution
grep -h _ $files | grep 'X N' | sort >! $1\.x.neu.bsl
cat $1\.r.neu.bsl | wbslfilt >! $1\.r.neu.wbsl
cat $1\.x.neu.bsl | wbslfilt >! $1\.x.neu.wbsl
cat $1\.r.neu.wbsl | grep RMS | colrm 132 | sort -rn +14 >! $1\.r.neu.rms
cat $1\.x.neu.wbsl | grep RMS | colrm 132 | sort -rn +14 >! $1\.x.neu.rms
#Plot files for free solutions
cat $1.r.neu.rms | awk '{print$17/1000. , $7*1000.} ' >! $1.r.n
cat $1.r.neu.rms | awk '{print$17/1000. ,$11*1000.} ' >! $1.r.e
cat $1.r.neu.rms | awk '{print$17/1000. ,$15*1000.} ' >! $1.r.u
cat $1.r.neu.rms | awk '{print$17/1000. ,$19*1000.} ' >! $1.r.l
#Plot files for fixed solutions
cat $1.x.neu.rms | awk '{print$17/1000. , $7*1000.} ' >! $1.x.n
cat $1.x.neu.rms | awk '{print$17/1000. ,$11*1000.} ' >! $1.x.e
cat $1.x.neu.rms | awk '{print$17/1000. ,$15*1000.} ' >! $1.x.u
cat $1.x.neu.rms | awk '{print$17/1000. ,$19*1000.} ' >! $1.x.l
#plot on screen
cat rep.gnu | sed s/NAME/$1/ | sed s/F/r/ | sed s/TITLE/"Repeatability for $1 (free) `date`"/ >! r.gnu; gnuplot
r.gnu
cat rep.gnu | sed s/NAME/$1/ | sed s/F/x/ | sed s/TITLE/"Repeatability for $1 (fixed) `date`"/ >! x.gnu;
gnuplot x.gnu
#plot on paper
grep -v continue r.gnu >! r.lgnu
#lasergnu -b -f r.lgnu -p -Psparc
grep -v continue x.gnu >! x.lgnu
#lasergnu -b -f x.lgnu -p –Psparc
当然这段脚本是不能成立的: 所以我们将其中的改成可以实现的部分,最后得到的$1.r.n $1.r.e $1.r.u $1.r.l利用gmt psxy进行画图 可以得到free solution的基线长度与基线长度分量的RMS关系以及fixed solution的关系。
其中主要的实现功能是wbslfilt功能
以上都可以利用求解基线重复性。
总结
1、使用globk的程序脚本
2、使用gamit的baeline脚本
3、使用wbslfilt命令
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 21:48
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社