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

博文

rnx2rtkp解基线和gmt绘图

已有 2332 次阅读 2020-5-22 12:35 |个人分类:GNSS|系统分类:科研笔记

rnx2rtkp解算批处理脚本:get_fix_pos.sh


#众数求解
awk -F ' ' '{if($6==1) print substr($3,0,length($3))}' | #查找出需要比对的数据
sort | #数据排序
uniq -c | #统计重复出现个数
awk 'BEGIN {max=0}{if($1>max) {max=$1;content=$2}} END {print content}'` #重复出现最大的数据


gmt6绘图bat脚本

echo off
for /f "tokens=1,2,5" %%i in (pos_comb_LNGT_CK_01_LNGT_JC_03) do echo %%i %%j %%k >> E.txt
for /f "tokens=1,3,6" %%i in (pos_comb_LNGT_CK_01_LNGT_JC_03) do echo %%i %%j %%k >> N.txt  
for /f "tokens=1,4,7" %%i in (pos_comb_LNGT_CK_01_LNGT_JC_03) do echo %%i %%j %%k >> U.txt 
echo on
gmt begin pos_fix png,pdf
 gmt set FONT_TAG 15p FONT_LABEL 15p
 gmt subplot begin 3x1 -Fs10c/5c -A -M0.2c -SRl -SCb -BWSrt 
 gmt subplot set 0 
 gmt basemap -R2079/2095/1714.30/1714.35 -JX10c/5c -Byaf+l"east" 
 gmt plot E.txt -Sc0.1c -Ey -W0.5p -Gred
  gmt subplot set 1 
  gmt basemap -R2079/2095/-193.93/-193.88 -JX10c/5c -Byaf+l"north"
  gmt plot N.txt -Sc0.1c -Ey -W0.5p -Ggreen
  gmt subplot set 2
  gmt basemap -R2079/2095/-20.30/-20.24 -JX10c/5c -Byaf+l"up" -Bxa5f1+l"GPS Week"
  gmt plot U.txt -Sc0.1c -Ey -W0.5p -Gblue
  gmt subplot end
gmt end
REM delete file
del E_Temp.txt
del N_Temp.txt
del U_Temp.txt



pos_fix.png



https://blog.sciencenet.cn/blog-858128-1234396.html

上一篇:awk处理数据笔记
下一篇:基于VS2017编译调试rnx2rtkp
收藏 IP: 119.79.228.*| 热度|

0

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

数据加载中...

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

GMT+8, 2024-11-19 16:48

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部