科学网 › 标签 › GLOBK

标签: GLOBK

相关帖子	版块	作者	回复/查看	最后发表

没有相关内容

相关日志

GLOBK分布/批处理平差处理流程: tao164411096 2018-10-15 15:42; 解算过程是从陈超的腾讯课程 https://ke.qq.com/course/283024 .以及尹业彪网上上传的公开视频学习的。感谢两位前辈。最重要的还是自己看说明书实际操作慢慢理解学习的。哈哈哈哈哈哈哈哈哈哈哈写一个流程保存下来，以防以后使用时忘记。 sh_glred批处理：修改 sites.dafaults ：在 glred 组合中包含的测站由 glrepu 标记，在 glorg 中用来定义参考框架的测站由 glreps 标记，要生成时间序列图的测站由 glts 标记。默认使用所有测站并绘图。注意：若在sites.defaults文件中使用glreps标记平差时的固定点，则在glorg_combNaNd中需将“stab_site clear”选项注释掉（前边加“x”）。分布处理流程 ######################################################################### 创建H文件文件夹/hfiles/ ，将所有H文件以及tables里的svnav.dat文件复制到hfiles里，使用如下命令查找所有h文件并复制到hfiles里，也可手动复制。将h文件转换成globk识别的二进制文件 htoglb 》》》短距离用glx（松弛固定解），长距离大型网用glr（松弛浮点解）拷贝控制文件。 globk.cmd中指定已知点坐标apr文件，指定约束apr_site（XYZ）apr_neu（NEU），已知点约束0.1mm。 glorg.cmd中修改已知apr文件，指定固定站。制作h文件列表gdl文件 5. 修改两个控制文件后运行 glred/globk 6 globk.prt globk.log expt.gdl globk.cmd sh_plotcrd 调用sh_baseline 和sh_tshist 绘制时间序列 sh_plotk 从k文件从计算钟参数 GLOBK 正常结算过程输入的已知点apr文件确定坐标框架和历元。 +++++++++++++++++++++++++++++++++++++++++++++ + 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); 个人分类: GAMIT/GLOBK|5123 次阅读|0 个评论

使用GLOBK批处理命令进行网平差: PurpleSky 2015-4-3 09:52; 1. 说明首先，你要先完成基线解算。 sh_glred 在工程级目录执行，但在解目录（ gsoln ）中运行，该目录也可以由 process.defaults 文件中 glbpth 参数指定。 2. 参数配置 2.1 globk_comb.cmd 和 glorg_comb.cmd 若这两个文件不存在， sh_glred 尽管可以生成他们，但推荐手动创建。运行命令“ sh_glred -cmd ”，将在 gsoln 文件夹下生成上述配置的模板文件以供修改。需在其中配置先验坐标文件（ .apr）、解算中要用到的测站、定义参考框架时用到的测站和约束等（若使用 IGS 站点作为固定点，则可以使用自带的 itrf 文件作为先验坐标）。 2.2 sites.dafaults 每步处理中所使用的测站都由 sites.dafaults 中标记指定：在 glred 组合中包含的测站由 glrepu 标记，在 glorg 中用来定义参考框架的测站由 glreps 标记，要生成时间序列图的测站由 glts 标记。默认使用所有测站并绘图。注意：若在sites.defaults文件中使用glreps标记平差时的固定点，则在glorg_combNaNd中需将“stab_site clear”选项注释掉（前边加“x”）。 3. 命令执行完成上述操作之后，你就可以进行批处理网平差了。在工程目录执行： sh_glred -syr1 doy1 yr2doy2 -expt expt -net networks -local -nettext char -yrextyear -ncomb num -stnfostation.info -cmd -opt A F H L U G E K C R 其中 -s 指定开始时间（ yr1 doy1 ）和结束时间（ yr2 doy2 ）， expt 是 4 个字符的工程名， network 表示包含 SOPAC 的 H- 文件， SOPAC 默认选项是 all （所有），或者仅组合某些解，假如 -local 被指定了，脚本会仅仅处理可以获得数据范围内的天；否则它会自动去 SOPAC 查找 H- 文件，处理所有指定天数内的数据。要处理的天也可用 -d yr days 或者采用 -r days 详细指定。可选参数 -ncomb 用来指定每一个组合包含多少天，默认为 1 天，但也可以用脚本产生每周或每月的本地或全球文件的平均值。下面是 opt 选项可选参数的详细释义（摘自GAMIT 10.35中文文档） R 在开始之前把旧的 h*.gl? 文件从 glfpth 目录下删除掉。 F 从 SOPAC 网站上下载全球 h 文件，根据 networks 参数项提供的名称搜索文件。 H 对所有现有的或者链接到 glfpth( 通常是 procdir/glbf) 的 ascii 文件运行 htoglb 。 LA 链接所有本地存档的 ascii H 文件 ( 即所有的 h ?.yyddd) 到组合中去。 LB 链接所有本地存档的 ascii H 文件 ( 即所有的 h ?.yyddd) 到组合中去。从 htnd 开始向下搜索，其中 htnd 是系统默认的。 LC 链接本地生成的组合二进制 H 文件 ( 名为 *.GLX) 。 U 对二进制 h 文件运行 hfupd( 尚未试 ) 。 G 运行 glred 以进行组合或重新组合。 E 运行 ensum 和 sh_baseline 进行绘图。 C 后压缩 ascii H 文件，删除所有的链接，分别复制 ascii H 文件和二进制 H 文件为 glbpth/ascii_yyyy 和 glbpth/bin_yyyy, 其中 glbpth 是由 process.defaults 指定的。 A 完成所有选项例： % sh_glred -s \$year1 \$sdoy \$year2 \$edoy -expt $expt -local -opt H G E $year1 和 $sdoy 为要进行网平差的开始年与年积日， $year2 和 $edoy 为结束年与年积日。 4. 解算结果命令运行完成后，在gsoln/文件夹下即可找到平差结果文件*.prt和*.org。固定点有“ * ”号标记。 PS：虽然这种方法可以快速得到网平差结果和站速度，但是在 GAMIT/GLOBK 软件的官方指南中对这种平差方法是不推荐的。推荐方法是只使用该方法合并 GPS 网得到站坐标，然后使用 globk 命令调用 globk_vel.cmd 和 glorg_vel.cmd 得到站速度（因整理时间较晚，有些可能记得不对，请读者自行查阅 GLOBK_Ref）。; 个人分类: GAMIT/GLOBK|8044 次阅读|0 个评论

GAMIT-GLOBK处理实例——从零开始，到解算结束: 热度 2 POPOPOLAR 2014-8-28 09:55; 确认 GAMIT/GLOBK 软件安装成功之后，为了能正确处理数据和绘图，还需要做两件事情。第一，更新表文件。在 ftp 服务器上下载最新的 tables 表文件（ ftp://garner.ucsd.edu/pub/gamit/tables/ ），然后复制到安装目录下的 tables 中，遇到相同文件选择覆盖即可；第二，安装并配置好绘图工具 GMT （ The Generic Mapping Tools ）。友情链接： Linux/Ubuntu 环境下轻松安装 GAMIT/GLOBK—— http://blog.sciencenet.cn/blog-1217335-814640.html Linux/Ubuntu 环境下轻松安装配置 GMT—— http://blog.sciencenet.cn/blog-1217335-758694.html 下面的实例中，处理的数据为 2010 年第 56 到 60 共 5 天的 3 个 IGS 跟踪站（ cas1 、 dav1 和 mac1 ）和 2 个待求站（ grw1 、 zhn1 ）的数据。一、处理前的准备 1 、在主文件夹内新建 nanj 项目文件夹，项目内新建 brdc 、 igs 、 rinex 三个文件夹，分别存放当天的广播星历、精密星历和观测值文件。注意：如果文件采用 .z 压缩格式，则应当使用 gunzip 命令进行解压；如果解压出的观测值文件仍采用 .d 的压缩格式，则应使用 crx2rnx 命令将其转化为 .o 的标准 RINEX 格式（ sh_crx2rnx -f *.*d 可用于 .d 到 .o 文件的批量转换）。 2 、终端进入 nanj 项目文件夹内链接 tables ，运行： sh_setup -yr 2010 3 、生成 station.info 文件。将 nanj/tables 下的 station.info 文件拷贝到 rinex 文件夹下，打开并编辑，仅保留以 # 或 * 开头的前几行，保存并关闭。打开终端并进入 rinex 文件夹，运行 sh_upd_stnfo -files *.10o 运行成功以后 station.info 文件便存放了此项目内的开始、结束时间，站名和接收机、天线类型等信息。 4 、建立 lfile. 文件。 lfile. 是测站的先验坐标文件。打开终端并进入 rinex 文件夹，用批处理的方式生成 lfile. 文件可以分为三步： a ）提取观测值 .o 文件的先验 XYZ 坐标 grep POSITION *.10o lfile.rnx b ）将 .rnx 文件转化为 .apr 文件 rx2apr lfile.rnx 2010 056 c ）由 .apr 文件生成 lfile. 文件 gapr_to_l lfile.rnx.apr lfile. 2010 056 将 station.info 和 lfile. 两个文件拷贝到 nanj/tables 文件夹，覆盖原文件。 5 、 sestbl. 的配置。 nanj/tables 下的 sestbl. 文件是测段分析策略文件，该文件内部有详细的说明。一般来说采用默认配置即可，通常需要修改的三个地方是： Choiceof Experiment 选择处理方式 Choice of Observable 选择观测值类型和模糊度解算 Use otl.grid 选择是否使用潮汐文件对于 Choice of Experiment ，选择 BASELINE 时将固定轨道并在 GAMIT 处理中和输出 h- 文件时忽略轨道参数；选择 RELAX 时将采用松弛解，合并全球 IGS h- 文件时需要。要想点位置精度高用 RELAX ；若目的是求基线后面平差则用 BASELINE 。在此实例中采用默认的 BASELINE 。对于 Choice of Observable ，选择 LC_AUTCLN 为采用宽巷模糊度值并用伪距在 autcln 中解算；对于小于几公里的基线，用 L1 和 L2 独立载波相位观测值（ L1,L2_INDEPENDENT ）或者仅用 LI （ L1_ONLY ），相比用无电离层组合（ LC_HELP ）可以减少噪声水平。对于 Use otl.grid ，由于这里已在 ftp 上更新 tables ，有了最新的 otl.grid 文件，所以这里选择 Y 。 6 、 sittbl. 的配置。 nanj/tables 下的 sittbl. 文件对各个测站的先验坐标（或钟差、大气模型等）进行约束。对高精度的已知坐标采取强约束，而对待求点采用松弛约束。如 IGS 站的坐标分量约束在较小的 1 ～ 75px ，对未知点的约束可以到 5 ～ 10m 。 7 、 sites.defaults 和 process.defaults 的配置。 sites.defaults 文件用来控制需要参与解算的测站。在文档末尾可以根据提示编辑，来给定那些不参与解算的测站或是测站的某些天。 process.defaults 文件用来控制处理过程中的很多细节，比如 sampling interval, number of epochs, start time for processing ， default globk .apr file 等等，根据需求和提示进行编辑。在这里均使用它们的默认值。二、利用 GAMIT 解算基线用终端进入 nanj 的项目文件夹，输入批处理命令进行解算： sh_gamit -expt nanj -s 2010 056 060 -orbitIGSF -yrext -noftp -dopt D ao c x sh_gamit.log 参数说明： -expt ：指定四个字符的项目名称 -d ：指定需要处理的指定日期，例如 -d 2010 56 60 ，指的是处理 2010 年第 56 和 60 天。 -s ：指定需要处理的时间序列，例如 -s 2010 56 60 ，指的是处理 2010 年第 56 到 60 天。 -orbit ：卫星轨道类型。 -yrext ：给日目录前添加年前缀，例如 2010_006 。 -noftp ：处理过程中不连接 ftp 下载数据。 -copt ：数据处理完成后待压缩的文件类型，例如 -copt o q m k x 。 -dopt ：数据处理完成后待删除的文件类型，例如 -dopt D ao c x 。结果文件将存放在名称为年积日的文件夹内，此例中，可供参考的结果文件为： sh_gamit_2010_006.summary 解算总结 qnanja.006 解算记录 onanja.006 解算记录的简略版，一般关注此文件 hnanja.10006 协防差矩阵、参数平差值在上述文件中，可根据描述来判断解算结果是否符合相应的需求。其中，基线解算结果（ O 文件，即这里的 onanja.006 ）中的 postfit_nrms 项优于 0.3 左右时最佳；如果大于 1.0, 则表示此解存在问题。三、利用 GLOBK 进行平差处理用终端进入 nanj 的项目文件夹内，运行： sh_glred -expt nanj -s 2010 050 2010 065 -yrext -opt H G E sh_glred.log 参数说明： H ：运行 htoglb ，把文件转换为二进制文件 G ：运行 glred ，合并文件 E ：绘图运行成功后，进入 gslon 文件夹，以 psbase 开头的文件即为各个站三维坐标时间序列的图形表示。用终端进入 gsoln 目录下，依次运行： ls../glbf/h*glx nanj.gdl glred6 globk_comb.prt globk_comb.log nanj.gdl globk_comb.cmd 生成的 globk_comb.org 文件中就包含了解算点的三维坐标和相关参数。至此，利用 GAMIT-GLOBK 解算 GPS 基线并进行平差的实例叙述完毕。本文出于帮助初学者利用 GAMIT-GLOBK 软件解算 GPS 基线并平差的目的撰写，但受笔者水平限制，难免存在疏漏和错误，欢迎大家批评指正。; 个人分类: GAMIT|17512 次阅读|2 个评论

globk GLORG/apply_cond_full: Failure to converge问题的解决: 热度 1 zhoufcumt 2011-10-25 20:18; 利用GAMIT/GLOBK软件运行sh_glred模块时，在工程目录下的gsoln目录出现一个GLORG.warning的警告文件，里面有一些类似WARNING:111025:1945:36.0 GLORG/apply_cond_full: Failure to converge这样的warning，应该是影响到kalman滤波的收敛，尽管是warning，但个人感觉会对结果精度产生影响。经过反复实验，发现问题出在glorg_comb.cmd命令文件下，当该文件下仅有 stab_site clear 时，将出现上述警告文件，但加上参考站时，如 stab_site clear stab_site albh algo brmu dav1 drao fair fort gold graz guam kosg mac1 # Set parameters to estimate in stabilization pos_org xrot yrot zrot xtran ytran ztran scale 上述警告文件将不出现，问题解决了！注（个人理解）：这个涉及到坐标框架的七参数转换，如果你不想进行坐标转换的话可以把上面的三项注释掉，这样运行也不会出现警告文件。涉及到坐标框架七参数转换的，必须要固定至少三个以上均匀分布的控制点！; 个人分类: Gamit|6398 次阅读|2 个评论

更多...

帐号		自动登录	找回密码
密码			注册

关闭 安全验证

标签: GLOBK

相关帖子

相关日志

关闭安全验证