|||
解完基线后的平差,我一般使用globk分布处理。
写了个批处理的globk脚本,根据输入的参数glred/globk,可以分别得到时间序列和速度并画图。
使用前需新建glred和globk文件夹,并将提前准备好的两个cmd和apr文件放到glred/globk文件夹内。
将脚本放到工程目录执行即可。
#!/bin/bash ############# ##提前准备好两个cmd,apr文件放到globk/glred文件里,之后执行脚本 if [ $# -eq 0 ] then echo '##########################################' echo '# #' echo '# Usage: bash sh_globk globk/glred #' echo '# #' echo '##########################################' exit fi ########### run=$1 if [ ${run} == 'glred' -o ${run} == 'globk' ];then if [ -e ./${run}/globk_简.cmd ];then mkdir hfiles cp ./tables/svnav.dat hfiles find ./ -name "hexpt*" | xargs -i cp {} ./hfiles/ htoglb ./glbf/ ./tables/svs_myexp.svs ./hfiles/h* cd ./${run} ls ../glbf/h*.glx >> expt.gdl ${run} 6 ${run}.prt ${run}.log expt.gdl globk_简.cmd #### glred for TimeSeries ################# if [ ${run} == 'glred' ];then echo -e "\033[42;37m "++++++++ TimeSeries +++++++"\033[0m" echo -e "\033[42;37m "++++++++ TimeSeries +++++++"\033[0m" mkdir plot cp ./${run}.org plot cd plot mkdir pos_xyz tssum ./pos_xyz cors -R glred.org xyzsave 0 ./pos_xyz glred.org sh_plot_pos -f ${run}.org gmt psconvert -Tt -A0.5c *.ps rm *.ps ##### globk for Vel ################### elif [ ${run} == 'globk' ];then echo -e "\033[42;37m "++++++++ globk +++++++"\033[0m" echo -e "\033[42;37m "++++++++ globk +++++++"\033[0m" grep 'Unc.' globk.org >> result.xyz mkdir vel cp globk.org vel cd vel sh_org2vel -file globk.org sh_plotvel -file globk.vel fi ###################################### else echo -e "\033[42;37m "+++++ do not exist ${run} or CMD files +++++"\033[0m" fi else echo . echo -e "\033[42;37m "+++++ please input glred/globk +++++"\033[0m" echo . fi
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-25 16:07
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社