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

博文

[转载]gamit基线解算说明

已有 8357 次阅读 2017-8-11 14:58 |个人分类:GAMIT/GLOBK|系统分类:科研笔记|文章来源:转载

sites.defaults准备介绍

sites.defaults文件是站处理控制文件。格式如下:
all_sites expt xstinfo  
brus_gps expt ftprnx ftpraw localrx xstinfo xsite glrepu glreps glts
其中
ftprnx ftpraw localrx xstinfo xsite为gamit选项,glrepu glreps glts为GLOBK选项,格式具体含义如下:
expt为4个字符的项目名称(工程名称)
ftprnx表示该站点的
rinex数据从FTP站点下载
ftpraw表示该站点的
raw数据从FTP站点下载
localrx表示站点的
rinex数据位本地数据
xstinfo表示该站点不自动更新
station.info表中的信息
xsite该站点的所有天或指定日期的数据不处理
glrepu表示站点用于计算重复解(默认全部)  
glreps用于定义站点的在
GLRED处理中的参考框架(默认为IGS站点列表)
glts在GLRED处理中以时间序列绘制站点图(默认全部)
 
一般通过自己编制的sites.defaults文件可以设置所有站点不自动更新station.info表中的信息。示例如下:
all_sites test xstinfo(test为工程名)
LHAS_GPS test localrx glrepu glreps  

SHAO_GPS test localrx glrepu


sittbl.编辑

打开sittble文件删除无关站点,增加工程站点设置各站点的约束属性。
对于IGS框架点或你需要约束的点设置为0.050 0.050 0.050,其他点100. 100. 100.
一般来说约束点为IGS站点且不少于3个。

要想点位置精度高则设置IGS站的约束小一些,即认为IGS站已经很准,强约束,若目的是求基线后面平差则可以将IGS站的约束设置大些。


process.defaults文件可修改部分的解释
process.default文件中需要修改的有采样间隔,最大历元数,开始时间等,示例如下:
set sint = '30' 采样间隔依据需要和rinex的情况可修改为15 30 60

set nepc = '2880' 最大历元数依据观测时长计算调整。
set stime = '0 0' 开始时间(小时分)依据rinex文件修改,一般要求时间与rinex时间差不大于0.1h否则报错。
set use_rxc = "N" 设置为N不依据rinex更新lfile文件。
set aprf = itrf08.apr 采用08框架数据(默认itrf2000) 如不采用GLOBK平差可不设置。
set minxf = '300' Xfile文件大小的最小值,小于该值的将被抛弃。依据历元数设置。历元数少时应减小该值。
set mailto = '***@163.com' 设置为有效邮箱以便接受gamit结算结果。

sestbl.文件依据实际需要进行修改。示例如下:
Choice of Experiment = BASELINE(计算高精度基线)
对于Choice of Experiment,选择BASELINE时将固定轨道并在GAMIT处理中和输出h-文件时忽略轨道参数;选择RELAX时将采用松弛解,合并全球IGS h-文件时需要。要想点位置精度高用RELAX;若目的是求基线后面平差则用BASELINE。在此实例中采用默认的BASELINE。
Choice of Observable = LC_AUTCLN 观测值类型选择

对于L1L2 _INDENT和L1 _ONLY 和LC_AUTCLN三种解,由于短基线两端观测环境相关性极强,
短基线宜采用差分模式消除观测误差的L1L2 _INDENT和L1 _ONLY 解,而不宜采用模型改正的LC_AUTCLN解。对于15km以下的短基线应选择L1,L2_INDEPENDENT 或L1_ONLY

Type of analysis = 0-iter 迭代多少次[计算结果作为初始值再次计算一般设置1-iter]
Zenith delay = N 天顶延迟对于短基线选N
Interval zen = 2  ;2 hrs = 13 knots/day (default is 1 ZD per day)天顶延迟修改为1h,默认2h


若软件限制了最大测站数据,最大卫星数目,最大对流层延迟参数,最大解算历元数目。更改于Makefile.config文件。然后重新编。


判断解算质量主要有两个方面:
一是数据利用率(即是否有足够多的可用于进行合理估计的数据参与计算);
二是数据的噪声水平是否与选用的模型匹配。(这一项主要有两个指标,基线分量误差的大小和nrms(标准化的均方根误差),对于基线分量误差的判断主要依据预先设定的站点坐标和轨道参数的先验约束,如果大于先验可以通过查看Q文件或者autcln.sum文件,是否有大量被autcln丢弃的数据。对于nrms一般应在0.25左右。如果大于0.5就意味着周跳没有消除,或出现了模型错误(固定站点坐标错误或者未模型化的卫星)等。解算结果文件autcln.sum.post中的rms残差为显示控制网站点的相对质量。打开sh_gamit_doy.summary,otesta.doy, qtesta.doy 文件查看解算总结和基线解算情况,解算成功与否判断:Postfit nrms: *****E+00 是否在0.25左右,最差不能大于0.5,基线分量的误差大小是否小于预设值。至此基线解算全部结束。

格式问题
o文件名需要小写,并且内容格式符合UINX 文件名改小写:
sh_casefold -dir d -files *
格式转UNIX: sh_dos2unix *


基线长度对基准点的选择
基线处理时应注意如下3点:  
1)对于短基线(小于10km), 其基准站坐标精度应优于1m,只需将基准站与国家A、B级网点联测0.5h即可得到。
2)对于中长基线(小于100km),其基准站坐标精度应优于0.1m,应取国家A级网点作为基准点,并与其同步联测24h以上。
3)对于长基线(大于100km),其基准站坐标精度应优于0.05m,应取全球站作为基准站,并与其同步联测48h以上。


各测站的近似坐标。GAMIT要求是球面坐标,可以先在测站的O文件里提取近似坐标XYZ,生成lfile.apr文件,然后在/tables目录下输入命令:gapr_to_l lfile.apr lfile. “” yr doy


不同测站坐标约束对基线(BJFS_KUNM)精度的影响:测站约束越严格,精度越高。
不同天顶延迟参数对基线精度的影响:一般情况下,天顶延迟参数取为13,即每2h产生一个天顶延迟估值是比较合适的。天顶延迟参数变化对长基线的精度影响不是非常明显。
不同截止高度角对基线精度的影响:不过高度角太低时路径延迟对水汽变化非常敏感,因此在一般的解算情况下选择10°的截止高度角比较适宜。




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

上一篇:坐标系转换投影方法图
下一篇:[转载]globk通用讲解
收藏 IP: 171.113.236.*| 热度|

0

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

数据加载中...
扫一扫,分享此博文

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

GMT+8, 2024-12-24 00:24

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部