Gentle''的个人博客分享 http://blog.sciencenet.cn/u/tao164411096 九分耕耘,一份收获。

博文

GLOBK分布/批处理平差处理流程

已有 1186 次阅读 2018-10-15 15:42 |个人分类:GAMIT/GLOBK|系统分类:科研笔记| GLOBK

 解算过程是从陈超的腾讯课程https://ke.qq.com/course/283024.以及尹业彪网上上传的公开视频学习的。   感谢两位前辈。

 最重要的还是自己看说明书实际操作慢慢理解学习的。哈哈哈哈哈哈哈哈哈哈哈

 写一个流程保存下来,以防以后使用时忘记。

sh_glred批处理:

修改sites.dafaults :在 glred 组合中包含的测站由 glrepu 标记,在glorg 中用来定义参考框架的测站由 glreps 标记,要生成时间序列图的测站由 glts 标记。默认使用所有测站并绘图。

注意:若在sites.defaults文件中使用glreps标记平差时的固定点,则在glorg_combNaNd中需将“stab_site clear”选项注释掉(前边加“x”)。

分布处理流程

  #########################################################################

  1. 创建H文件文件夹/hfiles/ ,将所有H文件以及tables里的svnav.dat文件复制到hfiles里,使用如下命令查找所有h文件并复制到hfiles里,也可手动复制。1.png



  2. 将h文件转换成globk识别的二进制文件  htoglb  》》》短距离用glx(松弛固定解),长距离大型网用glr(松弛浮点解)

    2.png

  3. 拷贝控制文件。3.png

    globk.cmd中指定已知点坐标apr文件,指定约束apr_site(XYZ)apr_neu(NEU),

    已知点约束0.1mm。 glorg.cmd中修改已知apr文件,指定固定站。

  4. 制作h文件列表gdl文件

4.png

  5. 修改两个控制文件后运行

 glred/globk  6  globk.prt globk.log expt.gdl globk.cmd


sh_plotcrd    调用sh_baseline 和sh_tshist 绘制时间序列

sh_plotk       从k文件从计算钟参数

5.png


GLOBK 正常结算过程

输入的已知点apr文件确定坐标框架和历元。

g1.png

g2.png






 +++++++++++++++++++++++++++++++++++++++++++++

 + GLOBAL KALMAN FILTER             Ver 5.30 +

 +++++++++++++++++++++++++++++++++++++++++++++



STATUS :181017:1322:48.0 GLOBK/main: Start: Processing list file (Name expt.gdl)

 Global    1

Cleaning and creating sort_file expt.srt

 READING APR_FILE            1  ./xyz.apr

**WARNING** 1083_GPS not in apr_files

**WARNING** WUHN_GPS not in apr_files

STATUS :181017:1322:48.0 GLOBK/main: Running forward Kalman filter

 Maximum derivatives found   5 with   3 coordinate and   0 EQ log terms

 Allocating         0.15 Mbytes of memory for run

 Initial memory mapped


 Updating SV ephemeris epoch by   113.6016 years

 Global    1 using     0.2 Mb. Running time     0.00 Scaling by      1.000 1.00000000

STATUS :181017:1322:48.0 GLOBK/glfor: File    1 using     0.15 Mb. Scaled by      1.000 1.00000000 (Name ../glbf/h1308081200_expt.glx)

 Retaining single entry multiday PMU

 For mexpta.220 Glbf ../glbf/h1308081200_expt.glx             Chi**2 NP    12 is      0.054

STATUS :181017:1322:48.0 GLOBK/glfor: Chi**2 Increment for    12 dof   0.054 (Name ../glbf/h1308081200_expt.glx)

 Orient_adj. (mas) 2013  8  8   0.01    9.95    0.02    9.95      -0.0      9.95


STATUS :181017:1322:48.0 GLOBK/main: Outputting filter run to file (Name globk.prt)

STATUS :181017:1322:48.0 GLOBK/main: Running glorg to output file (Name globk.org)

STATUS :181017:1322:48.0 GLOBK/glorg: Frame Realization

 Iteration   1 of Position stabilization expt.gdl


 Pivot element too small for row   6

 Value is    0.539499D-15, tolerancc is    0.100000D-11

 Iteration   2 of Position stabilization expt.gdl

 Iteration   3 of Position stabilization expt.gdl

 Iteration   4 of Position stabilization expt.gdl

STATUS :181017:1322:48.0 GLOBK/glorg: Applying equates

STATUS :181017:1322:48.0 GLOBK/glorg: Outputting solution

STATUS :181017:1322:48.0 GLOBK/glorg: Normal finish

 DELETING COM, SOL, SRT and MAKE_SVS FILES (del_scra)




http://blog.sciencenet.cn/blog-3391834-1140957.html

上一篇:画区域底图
下一篇:sh_gamit 的一些命令参数

0

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

数据加载中...

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2019-8-22 11:20

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部