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

博文

GAMIT基线结果评价以及问题解决

已有 1542 次阅读 2018-9-4 23:23 |个人分类:GAMIT/GLOBK|系统分类:科研笔记

GAMIT批量处理(sh_gamit)完成以后,在每日解文件夹中会生成一个统计信息文件:

sh_gamit_<DDD>.summary

一、好的结果如下:

Input options -d 2002 30 31 32 33 -expt ncar -pres ELEV -yrext -netext a

Processing 2002 031 GPS week 1151 4 Raw 2 

/data51/tah/SENH02/glob02/suomi/2002_031a

Disk Usage:  12678.4  Free  76447.4 Mbyte. Used 15%


Number of stations used 4 Total xfiles 4

Postfit RMS rms, to and by satellite

RMS  IT Site   All  01  02  03  04  05  06  07  08  09 ...

RMS  20 ALL    4.8   4   5   6   5   5   4   5   4   5 ...

Best and Worst two sites:

RMS  20 TMGO   3.2   3   3   4   4   4   3   3   3   4 ...

RMS  20 SA09   4.6   4   4   5   4   5   4   4   4   5 ...

RMS  20 PLTC   5.4   4   5   5   6   5   4   5   5   6 ...

RMS  20 SA13   5.5   5   5   6   5   5   5   5   5   6 ...


Double difference statistics

 Prefit nrms:  0.31280E+03    Postfit nrms: 0.21324E+00            (Constrained free)

 Prefit nrms:  0.31185E+03    Postfit nrms: 0.21818E+00            (Constrained fixed)

 Prefit nrms:  0.31272E+03    Postfit nrms: 0.20470E+00            (Loose free)

 Prefit nrms:  0.31185E+03    Postfit nrms: 0.20756E+00            (Loose fixed)

 Number of double differences:  12447

 Numbers of WL and NL biases 120   Percent fixed 95% WL  85% NL

Any large adjustments to positions (> 0.3 m) ...


(1)测站数与x文件数匹配

(2)Site postfit RMS values 3–10 mm范围

(3)没有测站的RMS = 0 ( 意味着autcln没有保留数据 ) 

(4)solve结果中Postfit nrms   ~ 0.2 不论约束解还是松弛解

(5)模糊度解算 (70–85%一般, > 90%很好)


二、差的结果如下

Input options -expt 0001 -d 2018 243 -orbit IGSF -pres ELEV -noftp -dopt D ao c

 

Processing 2018 243 GPS week 2016 5 Using node: cc Started at: 18_09_04_23:05:54

 

Processing directory: /home/chaoshu/project/taizhou/0001/243

Disk Usage:  140303  Free  80515.2 Mbyte. Used 64% 

  

Number of stations used 6 Total xfiles 6

  

Postfit RMS rms, total and by satellite

RMS  IT Site   All   01   02   03   05   06   07   08   09   10   11   12   13   14   15   16   17   18   19   20   21   22   23   24   25   26   27   28   29   30   31   32 

RMS  16 ALL    7.2   49    0    0    0    0   46   56  103    0   42    0    0    0    0  129    0   46    0    0    0   94    0    0    0    0  100   96    0   83    0    0  

Best and Worst two sites:

RMS  16 PI39   5.1    4    0    0    0    0    4    5    8    0    3    0    0    0    0    0    0    3    0    0    0    0    0    0    0    0    7    8    0    6    0    0  

RMS  16 PI38   6.0    3    0    0    0    0    4    4    0    0    5    0    0    0    0   11    0    4    0    0    0    8    0    0    0    0    8    9    0    7    0    0  

RMS  16 JM05   8.2    0    0    0    0    0    5    5    9    0    4    0    0    0    0   19    0    4    0    0    0   10    0    0    0    0   12    0    0    7    0    0  

RMS  16 II39   9.3    7    0    0    0    0    6    9   12    0    5    0    0    0    0   12    0    7    0    0    0   11    0    0    0    0   10   11    0   14    0    0  

 

Double difference statistics

 Prefit nrms:  0.32817E+01    Postfit nrms: 0.92325E+00

 Prefit nrms:  0.32304E+01    Postfit nrms: 0.91131E+00

 Prefit nrms:  0.32817E+01    Postfit nrms: 0.92298E+00

 Prefit nrms:  0.32304E+01    Postfit nrms: 0.90988E+00

 Number of double differences:   1618

 

 Phase ambiguities (Total  WL-fixed   NL-fixed): 42 42 41  AUTCLN (Total Expected)  42 42

 Phase ambiguities WL fixed 100.0% NL fixed  97.6% 

 

List of sites updated in the tables lfile. based on adjustments to coordinates > .3 m in the Prefit Solution

II39 GEOC LAT   -0.3184

PI37 GEOC LAT   -0.3935


(1)Postfit RMS rms出现了大量0和>10的情况;

(2)Postfit nrms > 0.5 (或者说>0.3)

(3)lfile.初始坐标误差太大,造成更新前后大于0.3m


说明

(1)Number of stations used 表示参与数据处理的测站数,Total xfiles表示创建x文件的总数。通常情况下,两者应该相等。如果参与数据处理的测站数小于创建x文件的总数,说明虽然某些x文件已经被创建, 但是由于它们观测时间太短,不满足process.defaults 中关于最短观测时间的设置,因此这部分站点不会参与数据处理。当然,如果是在sites.defaults中设置了xsite项,排除了某些站参与数据处理,也会造成创建的 x文件数大于参与数据处理的测站数。

(2)RMS那部分中,显示的是各个测站和卫星的均方根误差,紧随其后的是其中质量最好的两个测站和质量最差的两个测站相应的均方根误差。这部分中,第一列的数据单位是mm,后面的所有列的单位都是0. 1 mm。

一般最好的测站的RMS应该在3 mm〜5 mm; 最差的测站RMS应该在7 mm ~ 9 mm。如果在10 mm〜15 mm之间说明该测站有很大的噪声;如果RMS大于15 mm,说明该站的观测存在问题。这种问题可能是接收机的某种故障造成的,也有可能是因为测站附近多路径误差太大或者天气很恶劣。除此之外,如果测站的初始坐标有较大的误差或者观测时间太短以致引起收敛问题也会导致RMS值过大。

需要注意的是,如果最“好“的测站RMS值为0, 说明该站的数据已经移除。这种情况极有可能是因为测站的先验坐标质量太差,那么在atucln. Post.sum文件的最开始就可以看到伪距的RMS值很大。

(3)双差数据统计那四行表示的是不同情况下的nrms,依次是:约束浮动解;约朿同定解;松弛浮动解; 松弛同定解。其中左边的是先验值,右边的是最后确定值。

(4)最后一个部分是表示宽巷模糊度(WL)和窄巷模糊度(NL)固定的百分比。一般WL应该大于 90%,NL应该大于80%。



查看平均RMS和最好最坏的两个测站的RMS值,RMS值一般应该小于8mm,但是不能为0,最大不超过15mm。若最好的测站RMS值为0或者最差的RMS值大于15,则存在数据质量问题或者其他原因。


最好的两个站在3-5mm,最差的站点在7-9mm,若是在10-15mm说明是一个高的噪声,大于15mm代表一个很差的接收机,一个多路径影响很厉害,严重天气问题,或者autcln处理短间隔的问题应该查看天空图 与高度图,查看噪音来源,若两个最差站点噪音都大于15,则看autcln.post.rms查看其他的站点的rms

  grep RMS 027/sh*

在antcln.post.sum中看到伪距的rms范围在1-2m,大于5m代表这是坏的站点


A summary file should be emailed to you as each day completes execution.Check these files for number of stations (6), Postfit RMS (3-15 mm, none 0),postfit nrms (~0.2), ambiguity resolution and coordinate adjustments (< 30 cm).For this large network, the 2018 ambiguity resolution is 80-90% for GPS and Galileo, 40-50% for Beidou, but near zero for Glonass due to poor orbits),a condition also reflected in the sky plots in the /figs directory. Allowing the Glonass orbits to adjust (RELAX mode in the sestbl.) will improve the fit signficantly. To view the result of the fit of GAMIT's orbit model to
the input orbit, see the rms files in /igs. To view the complete editing statistics see autcln.prefit.sum and autcln.post.sum in the day directories; for the least-squares adjustment, see qeuraa.[ddd].



解决办法:

一旦你发现了解的结果有问题,你应该做出正确的行为并从合适的起点重新运行程序。如果问题所在是由于不好的数据或卫星引起,在用 autcln.cmd 中的edit_site_sv 命令剔除了所有不好的数据之后,你可以开始使用 autcln。如果问题是精度很差的先验坐标引起的,你应该在更新 L 文件或者 T 文件之后开始 model。如果你怀疑轨道不好,你可以排除一些卫星后重新 solve 来检验你的假设,然后添加 edit_site_sv 在 autcln.cmd 文件中,以便提供永久的记录和确保排除项。很少的情况下,不好的数据意味着 rinex 文件的问题或者表文件的问题。


edit_site_sv 

第一个变量是 4 个字母的台站代码,或者“all”;

第二个是卫星的 PRN 号,0 代表所有的卫星;

第三个和第四个是历元的跨度,在这个跨度之内数据需要被删除。



http://blog.sciencenet.cn/blog-858128-1132851.html

上一篇:北斗静态数据解算
下一篇:GNSS工程控制网解算——LGO基线+CosaGPS平差

0

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

数据加载中...

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

GMT+8, 2019-6-19 19:57

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部