第一范文网 - 专业文章范例文档资料分享平台

GAMIT-GLOBK处理实例 - 从零开始,到解算结束 文档

来源:用户分享 时间:2025/6/1 15:17:24 本文由loading 分享 下载这篇文档手机版
说明:文章内容仅供预览,部分内容可能不全,需要完整文档或者需要复制内容,请下载word后使用。下载word有问题请添加微信号:xxxxxxx或QQ:xxxxxx 处理(尽可能给您提供完整文档),感谢您的支持与谅解。

确认GAMIT/GLOBK软件安装成功之后,为了能正确处理数据和绘图,还需要做两件事情。第一,更新表文件。在ftp服务器上下载最新的tables表文件(ftp://garner.ucsd.edu/pub/gamit/tables/),然后复制到安装目录下的tables中,遇到相同文件选择覆盖即可;第二,安装并配置好绘图工具GMT(The Generic Mapping Tools)。 友情链接:

Linux/Ubuntu环境下轻松安装GAMIT/GLOBK—— Linux/Ubuntu环境下轻松安装配置GMT——

下面的实例中,处理的数据为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坐标 grepPOSITION *.10o lfile.rnx b)将.rnx文件转化为.apr文件 rx2aprlfile.rnx 2010 056 c)由.apr文件生成lfile.文件

gapr_to_llfile.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_combNaNd 生成的globk_comb.org文件中就包含了解算点的三维坐标和相关参数。 至此,利用GAMIT-GLOBK解算GPS基线并进行平差的实例叙述完毕。

本文出于帮助初学者利用GAMIT-GLOBK软件解算GPS基线并平差的目的撰写,但受笔者水平限制,难免存在疏漏和错误,欢迎大家批评指正。

GAMIT-GLOBK处理实例 - 从零开始,到解算结束 文档.doc 将本文的Word文档下载到电脑,方便复制、编辑、收藏和打印
本文链接:https://www.diyifanwen.net/c2zjqu77rp97u3cn9al3a_1.html(转载请注明文章来源)
热门推荐
Copyright © 2012-2023 第一范文网 版权所有 免责声明 | 联系我们
声明 :本网站尊重并保护知识产权,根据《信息网络传播权保护条例》,如果我们转载的作品侵犯了您的权利,请在一个月内通知我们,我们会及时删除。
客服QQ:xxxxxx 邮箱:xxxxxx@qq.com
渝ICP备2023013149号
Top