处理过程中，为了自己的方便使用，需要编写自己的脚本（用红色标出）

求解基线估计也可以通过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命令