GAMIT使用手册Word文档格式.docx
- 文档编号:22435391
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:14
- 大小:3.60MB
GAMIT使用手册Word文档格式.docx
《GAMIT使用手册Word文档格式.docx》由会员分享,可在线阅读,更多相关《GAMIT使用手册Word文档格式.docx(14页珍藏版)》请在冰豆网上搜索。
GAMIT每个时段观测数据要求的周期最长为1个UTC天,即从UTC的0点到24点(北京时间8:
00~24:
00),原则上不要跨天作业。
GAMIT软件的组成结构见图,它由不同功能模块组成,主要包括数据准备、生成参考轨道、计算残差和偏导数、周跳检测与修复、最小二乘平差等模块,这些模块即可以单独运行,也可以用批处理命令联在一起运行,最大限度地减少人为操作,提高运算效率。
软件的执行程序放在/com、/kf/bin和/gamit/bin三个目录下。
图GAMIT/GLOBK软件组成结构图
二批处理实例
确认GAMIT/GLOBK软件安装成功之后,为了能正确处理数据和绘图,还需要做两件事情。
第一,更新表文件。
在ftp服务器上下载最新的tables表文件(,然后复制到安装目录下的tables中,遇到相同文件选择覆盖即可;
第二,安装并配置好绘图工具GMT(TheGenericMappingTools)。
下面的实例中,处理的数据为2010年第56到60共5天的3个IGS站(BJFS、和WUHN、KUNM)和2个待求站(chdu、pixi)的数据。
一、处理前的准备
1、在主文件夹内新建test项目文件夹,项目内新建brdc、igs、rinex三个文件夹,分别存放当天的广播星历、精密星历和观测值文件。
注意:
如果文件采用.z压缩格式,则应当使用gunzip命令进行解压;
如果解压出的观测值文件仍采用.d的压缩格式,则应使用crx2rnx命令将其转化为.o的标准RINEX格式(sh_crx2rnx-f*.*d可用于.d到.o文件的批量转换)。
2、终端进入test项目文件夹内链接tables,运行:
sh_setup-yr2010
3、生成文件。
将test/tables下的文件拷贝到rinex文件夹下,打开并编辑,仅保留以#或*开头的前几行,保存并关闭。
打开终端并进入rinex文件夹,运行
sh_upd_stnfo-files*.10o
运行成功以后文件便存放了此项目内的开始、结束时间,站名和接收机、天线类型等信息。
4、建立lfile.文件。
lfile.是测站的先验坐标文件。
打开终端并进入rinex文件夹,用批处理的方式生成lfile.文件可以分为三步:
a)提取观测值.o文件的先验XYZ坐标
grepPOSITION*.10o>
b)将.rnx文件转化为.apr文件
rx2apr2010056
c)由.apr文件生成lfile.文件
gapr_to_llfile.2010056
将和lfile.两个文件拷贝到test/tables文件夹,覆盖原文件。
5、sestbl.的配置。
test/tables下的sestbl.文件是测段分析策略文件,该文件内部有详细的说明。
一般来说采用默认配置即可,通常需要修改的三个地方是:
ChoiceofExperiment
选择处理方式
ChoiceofObservable
选择观测值类型和模糊度解算
Use
选择是否使用潮汐文件
对于ChoiceofExperiment,选择BASELINE时将固定轨道并在GAMIT处理中和输出h-文件时忽略轨道参数;
选择RELAX时将采用松弛解,合并全球IGSh-文件时需要。
要想点位置精度高用RELAX;
若目的是求基线后面平差则用BASELINE。
在此实例中采用默认的BASELINE。
对于ChoiceofObservable,选择LC_AUTCLN为采用宽巷模糊度值并用伪距在autcln中解算;
对于小于几公里的基线,用L1和L2独立载波相位观测值(L1,L2_INDEPENDENT)或者仅用LI(L1_ONLY),相比用无电离层组合(LC_HELP)可以减少噪声水平。
对于Use,由于这里已在ftp上更新tables,有了最新的文件,所以这里选择Y。
6、sittbl.的配置。
test/tables下的sittbl.文件对各个测站的先验坐标(或钟差、大气模型等)进行约束。
对高精度的已知坐标采取强约束,而对待求点采用松弛约束。
如IGS站的坐标分量约束在较小的1~75px,对未知点的约束可以到5~10m。
7、和的配置。
文件用来控制需要参与解算的测站。
在文档末尾可以根据提示编辑,来给定那些不参与解算的测站或是测站的某些天。
文件用来控制处理过程中的很多细节,比如samplinginterval,numberofepochs,starttimeforprocessing,defaultglobk.aprfile等等,根据需求和提示进行编辑。
在这里均使用它们的默认值。
二、利用GAMIT解算基线
用终端进入test的项目文件夹,输入批处理命令进行解算:
sh_gamit-expttest-s2010056060-noftp-doptDaocx>
&
参数说明:
-expt:
指定四个字符的项目名称
-d:
指定需要处理的指定日期,例如-d20105660,指的是处理2010年第56和60天。
-s:
指定需要处理的时间序列,例如-s20105660,指的是处理2010年第56到60天。
-orbit:
卫星轨道类型。
-yrext:
给日目录前添加年前缀,例如2010_006。
-noftp:
处理过程中不连接ftp下载数据。
-copt:
数据处理完成后待压缩的文件类型,例如-coptoqmkx。
-dopt:
数据处理完成后待删除的文件类型,例如-doptDaocx。
结果文件将存放在名称为年积日的文件夹内,此例中,可供参考的结果文件为:
解算总结
解算记录
解算记录的简略版,一般关注此文件
协防差矩阵、参数平差值
在上述文件中,可根据描述来判断解算结果是否符合相应的需求。
其中,基线解算结果(O文件,即这里的)中的postfit_nrms项优于左右时最佳;
如果大于,则表示此解存在问题。
三、利用GLOBK进行平差处理
用终端进入test的项目文件夹内,运行:
sh_glred-expttest-s20100502010065-yrext-optHGE>
H:
运行htoglb,把文件转换为二进制文件
G:
运行glred,合并文件
E:
绘图
运行成功后,进入gslon文件夹,以psbase开头的文件即为各个站三维坐标时间序列的图形表示。
用终端进入gsoln目录下,依次运行:
ls../glbf/h*glx>
glred6
生成的文件中就包含了解算点的三维坐标和相关参数。
至此,利用GAMIT-GLOBK解算GPS基线并进行平差的实例叙述完毕。
三分布处理实例
1.在桌面上新建一个test的文件夹,在test中新建名称为brdc(广播星历)、igs(存放精密星历)、rinex(下载的观测值文件和自己要解算的观测值文件)。
2.用doy在test的终端里测试,比如输入:
doy20130201查看要解算的是GPS第几周,第几天,以便下载数据用。
3.在终端输入gftp,下载相关的文件,一定要选对自己解算需要的数据。
一般下载*.文件(比较小,好下载)。
下载完所需要的数据后,要解压:
语句:
gunzip*.Z
4.新建一个文件夹,比如:
名字为032,把要解算的测站o文件和n文件,以及精密星历文件放进来。
把032和test进行如下链接。
5.打开rinex文件,输入下面的语句,会产生一个文件:
查看文件中各点的坐标是否正确:
不能出现坐标值为0。
6.输入rx2apr2013032,会多出来一个名为:
的文件,
接着输入:
gapr_llfile.“’2013032生成一个用纬度和精度来表示点的坐标的文件(lfile.),并把它复制到tables文件夹。
7.在tables文件夹中新建一个文件sitelist,在里面输入参考站点的名称,每行的第一列是空格,从第二行开始输。
8.在tabie文件夹启动一个终端,输入:
sh_upd_stnfo-lsitelist。
9.会生成一个文件,把它给成保存后,再输入:
sh_upd_stnfo-files../rinex/*13o
10.在032文件中查看的链接是否正确,然后做一个l文件的链接,打开终端输入:
ln-s../tables/lfile../到此准备工作完成。
11.在032文件打开终端,输入:
makexp然后按照以下输入。
12.按照提示输入:
历元间隔这里是30s,开始时间,如果从0000开始一天正好是2880个历元。
13.
精密星历
接下来按照提示输入:
sh_sp3fit-figsoigsft利用sp3精密星历生成g文件(上面已经下载好了),并且同时生产t-和u-文件
14.
sh_check_sess–sess32–typegfile-file
检查g-文件
15.makej利用导航文件生成j-文件
16.sh_check_sess–sess32–typejfile–file
检查j-文件
17.
makex生产x-文件
18.fixdrv生成批处理文件
19.csh执行批处理脚本,得到q-文件,即解算结果。
四精度评估
(1)GPS常用评价站坐标精度的指标是多时段基线重复性和多时段坐标重复性。
基线重复性计算模型:
()
上式中
为观测时段,
为基线s的重复性统计值,
为第
时段基线s的中误差,
时段基线解算结果,
为总的时段数,
为基线结果的加权平均值。
有
()
基线结果的评价还以基线长度与误差的关系来衡量,有两种形式:
这里的a、b、c、d四个参数可以根据基线长度和它的重复性统计值
拟合出来。
也有以基线的三个分量分别检验它的重复性的,这三个分量是南北分量
、东西分量
和高程分量
。
计算公式与基线重复性的公式()相同,只是用基线分量代替基线长度而已。
坐标分量重复性计算模型也与()相同,只是用坐标分量和坐标分量中误差代替基线长度和基线长度中误差计算而已。
(2)通过天空图评定精度
如上图:
同一个地点在各个的时间段有高残差预示着存在多路径效应。
同一个地点在特定的时间段有高残差预示着存在水波折射。
(3)通过相位与高度角评定精度
如上图,错误:
整体趋势线在有的高角度离开蓝色线,表明错误引用了天线相位模型,或卫星信号受到了电磁干扰。
如上图,正常情况下:
红色整体趋势线平顺,在蓝色线上下均匀波动。
(4)查看qexpta文件:
正常Adjust不等于0,若等于0,则该测点未参与计算。
(5)通过基线标准差、相关系数评定精度
正常NE<
1cm,U<
2cm,如果NEU>
3cm,则数据质量有问题。
(6)通过评定精度
标准均方根nrms解算结果为左右;
若是nrms大于,表明处理结果是有问题的(例如:
周跳没有修复;
测站的起算点坐标有问题等);
若是nrms小于,则认为成功解算,否则需要检查原因,重新处理。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- GAMIT 使用手册
![提示](https://static.bdocx.com/images/bang_tan.gif)