DIY编程器网

 找回密码
 注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 640|回复: 0
打印 上一主题 下一主题

[GPS] 区域CORS网络增强PPP天顶对流层延迟内插建模

[复制链接]
跳转到指定楼层
楼主
发表于 2014-10-1 06:14:14 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
第38卷第6期
2013年6月
武汉大学学报· 信息科学版
Geomatics and Information Science of Wuhan University
Vol.38No.6
June 2013
收稿日期:2013-01-28。
项目来源:国家自然科学基金资助项目(41074024);中央高校基本科研业务费专项资金资助项目(2012214020207)。
文章编号:1671-8860(2013)06-0679-05 文献标志码:A
区域CORS网络增强PPP
天顶对流层延迟内插建模
张小红1 朱 锋1 李 盼1 翟 广1
(1 武汉大学测绘学院,武汉市珞喻路129号,430079)
摘 要:研究了利用区域CORS网基站上解算的高精度ZPD内插流动站ZPD的空间回归模型的适用性和内
插效果,证明了H1QM3模型与MLCM 模型等价。在此基础上,提出利用高差约束改善内插精度的空间回归
模型。实验表明,H1QX1模型和H1QM3模型稳定,在平缓区域能获得10mm以内的插值精度,在起伏地区
可获得20mm以内的插值精度;使用高差约束后,高程变化剧烈地区ZPD的内插精度可改善60%;内插参考
站数量为9或10时效果最佳。
关键词:CORS;ZPD;空间内插;增强;PPP
中图法分类号:P228.42
  精密单点定位(PPP)技术利用IGS发布的精
密轨道和钟差数据,使用双频载波相位和伪距组
合观测值,可提供高精度的静态定位和动态定位
服务[1]。然而,PPP定位受卫星星座几何结构、
观测条件、数据质量等影响,传统的PPP定位收
敛时间长,至少需要30min才能获得10cm的精
度[2-3]。为了加快PPP定位的收敛速度并提高其
定位精度,不少学者提出了利用IGS跟踪站数据
估计载波相位的小数周偏差以实现非差模糊度固
定解的方法[4]。在此基础上,Li等提出利用区域
CORS网络增强PPP的方法[3,5-6],并逐步发展成
目前的PPP-RTK技术[5,7]。区域CORS网络增
强PPP方法的核心是利用区域CORS网络参考
站数据实时计算参考站天顶对流层参数,采用合
适的数学模型内插出流动用户的天顶对流层延迟
改正数,增强流动站PPP定位,加快PPP模糊度
固定并缩短定位收敛时间[1,3]。误差改正内插模
型在网络RTK 中研究较多,目前主要有线性内
插模型(LIM)、低阶趋势面模型(LSM)、反距离
内插模型(IDW)、线性组合模型(LCM)以及最小
二乘配置模型(LSC)。Dai等通过研究和比较,认
为这些方法的内插效果并没有明显区别[8]。上述
模型基本上是针对网络RTK技术中双差误差改
正提出的。熊永良提出含高程因子的非差对流层
内插模型[9],Li等提出改进的LCM 方法(MLCM)
[3]用于非差对流层内插。
本文着重分析空间回归模型的适用性及其内
插效果,在此基础上提出高差约束的空间回归模
型。通过实验,给出空间回归模型中所需的最优
参考站数量配置方案。
1 ZPD空间回归模型
在局部区域,大气气象因素在高程方向的变
化较水平变化更为显著。图1给出了北美CORS
网共1 295个测站计算得到的单天ZPD均值与
测站高程的关系。从图中可以看出,测站ZPD与
高程存在明显的负相关,在建立ZPD回归模型
时,需要考虑测站高程因子。
图1 PPP-ZPD随高程变化趋势
Fig.1 Correlation of ZPD and Station Height
武汉大学学报· 信息科学版2013年6月
现有文献及笔者大量实验表明,表1中的6
类空间回归模型在进行流动站天顶对流层延迟误
差内插时效果较好,本文将详细比较这6类模型
的适用性和内插效果。
表1 6类ZPD空间回归模型
Tab.1 Six Spatial Regression Models for ZPD
模型表达式
H1QM3 ZPD=a0+a1B+a2L+a3H
H1QM4 ZPD=a0+a1B+a2L+a3BL+a4H
H1QM6
ZPD=a0+a1B+a2L+a3BL+a4B2+
 a5L2+a6H
H1QX1 ZPD=a0+a1H
H3QM3
ZPD=a0+a1H+a2L+a3BH+a4LH+
 a5H
H3QX1 ZPD=a0+a1BH+a2LH+a3H
  表1中,ai
为回归系数,B 为大地纬度,L 为
大地经度,H 为大地高。内插时,以流动站为中
心,建立局部坐标系,B、L、H 均为局部坐标系下
的值。显然,流动站的B、L、H 为0,因此,
ZPDu =a0 (1)
  事实上,MLCM 与H1QM3模型等价。
证明 MLCM 模型表达式为:
Σ n
i=1
αi =1,Σ

i=1
X =0,Σ

i=1
α2i
= min (2)
ZPDu =Σ

i=1
αiZPDi (3)
其中,X 表示B、L、H,则有矩阵形式:
1 1 … 1 1
B1 B2 … Bn-1 Bn
L1 L2 … Ln-1 Ln
H1 H2 … Hn-1 H



燅n
α1
α2

α



n燅









(4)
以向量替代,式(4)改写为:
S·α =e1,αTα = min,ZPDu =αTZPD (5)
对式(5)前两式使用拉格朗日极值法可求得α 并
代入第三式,结果为:
ZPDu =eT
1 (S·S ) T -1 S·ZPD (6)
H1QM3模型表达式:
ZPD1
ZPD2

ZPD



燅n

1 B1 L1 H1
1 B2 L2 H2
   
1 Bn Ln H



燅n
a0
a1
a2




燅3
(7)
ZPDu =a0 =eT·a (8)
使用最小二乘估计可得:
a= (S·S ) T -1 S·ZPD (9)
则:
ZPDu =eT (S·S ) T -1 S·ZPD (10)
  结果表达式一致,两者等价。从本质上讲,它
们都是位置的响应函数,但H1QM3模型比MLCM
具有更明确的数学意义。
一般情况下,对流层延迟在区域分布中是符
合连续平缓性、空间相关性和各向同性的,本文讨
论的区域都是基于这种假设。对于特殊气候条件
或特殊局部环境区域下的内插效果,需要进一步
分析验证。
2 高差平缓地区ZPD内插实验
选择北美CORS网北纬41°~44°、西经82°~
87°区域(300km×350km)作为实验区。该地区
地势平坦,除去部分数据不足的站点,剩下共86
个站点,选取5个站点作为流动站(圆点),其余为
参考站(三角形),如图2所示。
图2 高差平缓区域测站分布图
Fig.2 Distribution of CORS Stations in a Flat Region
用户站ZPD内插时,通过距离来选择参考站,
每个历元选择距离最近的10个有效参考站,平均
距离在40km左右,分别使用以上6种模型进行局
部拟合求得用户站ZPD内插值。将PPP计算得到
的用户站的ZPD作为真值[10,11],计算内插外符合精
度(RMS),同时给出了模型的内符合精度(IP)和交
叉验证值(CV)[12],见图3。其中内符合精度是指参
考站拟合残差的均方根。交叉验证是指略去一个参
考站,利用其余9个参考站内插该参考站的ZPD值,
得到内插值与真值的偏差。对所有参考站使用该法
共得到10个偏差,其均方根即为交叉验证值。
图3 6类模型内插结果比较
Fig.3 Comparison of Interpolation Results
between Six Models
680
 第38卷第6期张小红等:区域CORS网络增强PPP天顶对流层延迟内插建模
6类模型RMS均在10mm 以下,能满足区
域CORS网增强PPP的精度要求,其中只顾及高
程因子的H1QX1模型最佳,为6.2mm;H1QM3
模型次之,但同时顾及了平面因素。内符合精度
IP表征了内部拟合程度。采用不同的内插模型,
IP值与RMS值相关性不强;对相同数量测站的
拟合,模型参数越多,IP值越小。因此,IP值作为
内插评价指标并不可靠。而CV 值变化基本与
RMS一致,可以辅助评价内插效果。
取各站H1QX1模型内插序列,单日共2 880
个(30s采样),统计绝对偏差在10mm、20mm、
30mm以内的百分比,如表2所示。
表2 H1QX1模型绝对偏差统计/%
Tab.2 Satistics of Absolute Bias for H1QX1Model(%)
绝对偏差bayr hruf mico mipp mipv
<10mm 91.2 88.3 90.1 83.2 83.22
<20mm 100 98.7 98.2 99.2 98.2
<30mm 100 100 100 100 100
  所有内插站绝对偏差不超过30mm,大部分
在10 mm 以内,其中mipv 站稍差,但使用
H1QM3模型时各为92.5%、100%、100%,精度
提高显著。且此时的CV值同样小于H1QX1模
型,说明同一站可以备选几种模型,通过CV值来
自动选择最佳模型。bayr站的内插精度最高,从
图4中可以看出,ZPD 内插序列平滑,与PPPZPD
基本一致,其精度基本在10mm 以内,精度
最差处为18mm,在历元50~150这段时间的精
度稍差,超过10mm,大部分时间内内插精度优
于10mm。
图4 bayr站ZPD内插日变化
Fig.4 Daily Change of Interpolated ZPD
on bayr Station
3 高差剧烈地区ZPD内插实验
选择北美CORS站北纬32°~36°、西经114°
~122°区域(400km×700km)作为实验区,该地
区地势变化较为剧烈,最大高差为2 574m。除
去部分数据不足的站点,剩下共67个站点,选取
5个站点作为流动站(圆点),其余为参考站(三角
形),如图5所示。
图5 高差剧烈区域测站分布图
Fig.5 Distribution of CORS Stations in a
Undulated Region
内插仍为局部内插,有效参考站数量为10
个,平均站间距为65km,平均最大高差在1 500
m,高差变化剧烈。6种模型内插结果如图6。
图6 6类模型内插结果比较
Fig.6 Comparison of Interpolation Results
between Six Models
高差剧烈结果比高差平缓结果差,表明高差
因子对内插有显著影响。其中H1QM4精度最
好,为16.7mm;H1QM3与H1QX1次之,均为
17.9mm。
表3 H1QM4模型绝对偏差统计/%
Tab.3 Statistics of Absolute Bias for H1QM4Model(%)
绝对偏差bar1 cru1 monp p282 pin1
<10mm 57.5 35.5 4.5 94.0 60.0
<20mm 96.2 66.7 31.6 100 86.3
<30mm 100 88.6 84.3 100 99.2
  monp站精度较差,偏差几乎在10mm以上,
20mm以内的偏差也只有31.6%,这是因为使用
的参考站高差较大。通过高差约束来选择参考
站,在距离最近的若干参考站中,选择高差小于
珨H
的参考站,约束后,monp最大高差由2 000m
降为850m,但站间距离由81km 拉长到340
km。
高差约束前后的内插结果比较如图7所示。
在高差约束后,monp站精度大幅提高,3种模型
681
武汉大学学报· 信息科学版2013年6月
的内插效果改善显著,RMS 分别由约束前的
24.5mm、30.1mm、31.3mm改善为17.9mm、
13.7mm、12.5mm,其中H1QX1模型最佳,偏
差百分比改善分别为61.0%、76.3%和52.3%。
因此,在高程变化剧烈的局部区域,通过高差约
束,可以显著改善ZPD的内插效果。当然,选取
参考站时,高差约束和参考站平面距离两个因素
需要综合考虑,因为高差约束越紧,在需要相同数
量的参考站进行内插时,参考站间的平面距离就
会越大,对流层的空间相关性也会变弱。
图7 高差约束与无约束结果比较
Fig.7 Comparison between Height Constraint and
Without Height Constraint
无论在高差平缓还是高差起伏区域,所选的
流动站是随机的,其参考站几何分布都不同,在内
插参数相同的条件下,精度并无多大差异。从理
论上讲,几何分布对内插方程的系数阵没有影响,
因此,参考站的分布对于内插精度没有很大影响。
4 参考站数量对内插结果的影响
参考站数量对内插结果也有影响。基于上面
的分析,整体上H1QX1模型和H1QX3模型内
插效果最优。因此,本文基于H1QX1模型来分
析参考站数量对内插结果的影响。分别计算了高
差平缓区域的bayr站与高差剧烈区域的monp
站的结果,统计其RMS、最大偏差指标(Max-
Bias)和内插误差的百分比,如图8和图9所示。
图8 bayr站不同数量参考站的内插效果比较
Fig.8 Comparison of Interpolation Performance
between Different Number of Reference Stations
at bayr Station
图9 monp站不同数量参考站的内插效果比较
Fig.9 Comparison of Interpolation Performance
between Different Number of Reference Stations
at monp Station
内插评判标准是:RMS与MaxBias值越小
越好,内插误差小的历元所占的百分比值则越大
越好。从图8可以看出,在平缓地区,参考站数大
于10时,精度不再发生变化,参考站数为9时,<
10mm 的百分比较高,而参考站数为10时的
MaxBias、<20mm以及<30mm 指标较优。总
体来说,参考站数为9或10时较好。图9中的高
差剧烈地区也有相同的结果,所有指标在参考站
数为9或10时达到最优,而两边为下降趋势。从
上述结果来看,无论是高差平缓还是剧烈地区,使
用9或10个参考站内插流动用户的天顶对流层
效果最好。对于平缓地区,继续增加参考站数量,
精度指标不再变化,却会增加计算量;而高差剧烈
地区,新加入的参考站距离拉长,空间相关性随之
降低,内插精度会受其影响。
5 结 语
本文使用空间回归模型内插ZPD取得了较
好的结果,并得出以下结论:
1)在局部区域,ZPD变化与高程因子具有强
相关性,平面距离越短,高差越小,参考站与流动
站相关性越强,内插效果越好。
2)H1QM3模型与MLCM 模型等价,且具
有更明确的数学意义,H1QM3模型求解简单。
3)无论是高差平缓还是高程变化剧烈地区,
H1QX1模型和H1QM3模型都有较高的内插精
度,且相对稳定。高差平缓地区为6.2mm和6.7
mm,高差剧烈地区为17.9mm和18.4mm,因此
建议局部区域使用这两种模型进行拟合内插。
4)若参考站高差过大,流动站内插效果不
佳,建议使用高差约束,减小参考站的高差,同时
综合考虑平面距离。高差约束后,H1QX1 和
H1QM3模型内插精度改善明显。
682
 第38卷第6期张小红等:区域CORS网络增强PPP天顶对流层延迟内插建模
5)参考站数量过多、过少对内插精度都有影
响,当数量为9或10时,内插效果最佳。
参 考 文 献
[1] Bisnath S,Gao Y.Current State of Precise Point
Positioning and Future Prospects and Limitations
[J].Observing Our Changing Earth,2008:615-
623
[2] Mervart L,Zdenek L,Rocken C,et al.Precise Point
Positioning with Ambiguity Resolution in Real-time
[C].GPS Solutions Inc,Boulder Co,USA,2008
[3] Li X,Zhang X,Ge M.Regional Reference Network
Augmented Precise Point Positioning for Instantaneous
 Ambiguity Resolution[J].Journal of Geodesy,
2011,85(3):151-158
[4] Ge M,Gendt G,Rothacher M,et al.Resolution of
GPS Carrier-phase Ambiguities in Precise Point Positioning(
PPP)with Daily Observations[J].Journal
 of Geodesy,2008,82(7):389-399
[5] Geng J,Teferle F N,Meng X,et al.Towards PPPRTK:
Ambiguity Resolution in Real-time Precise
Point Positioning[J].Advances in Space Research,
2011,47(10):1 664-1 673
[6] Li H,Chen J,Wang J,et al.Network Based Realtime
 Precise Point Positioning[J].Advances in
Space Research,2010,46(9):1 218-1 224
[7] Bena G,Schmitz M,Bagge A.PPP-RTK:Precise
Point Positioning Using State-space Representation
in RTK Networks[C].ION GNSS-05,Long Beach,
California,2005
[8] Dai L.Augmentation of GPS with Glonass and
Pseudolite Signals for Carrier Phase Based Kinematic
 Positioning[M].University of New South Wales,
2002
[9] 熊永良,黄丁发,丁晓利,等.基于多个GPS基准站
的对流层延迟改正模型研究[J].工程勘察,2005
(5):55-57
[10]叶世榕,张双成,刘经南.精密单点定位方法估计对
流层延迟精度分析[J].武汉大学学报·信息科学
版,2008(8):788-791
[11]张小红,何锡扬,郭博峰,等.基于GPS非差观测值
估计大气可降水量[J].武汉大学学报·信息科学
版,2010(7):806-810
[12]陈正昌.多变量分析方法统计软件应用[M].北
京:中国税务出版社,2005:104-106
第一作者简介:张小红,教授,博士生导师。现主要从事GNSS精
密定位及GNSS地学应用等方面的教学和科研工作。
E-mail:xhzhang@sgg.whu.edu.cn
Zenith Troposphere Delay Interpolation Model for
Regional CORS Network Augmented PPP
ZHANG Xiaohong1 ZHU Feng1 LI Pan1 ZHAI Guang1
(1 School of Geodesy and Geomatics,Wuhan University,129Luoyu Road,Wuhan 430079,China)
Abstract:Taking apriori precise ZPD derived from CORS network to augment PPP can reduce
 the number of unknown parameters,accelerate its ambiguity resolution in PPP and
shorten the convergence time.We investigate the applicability and interpolation effect of spatial
 regression model in which precise ZPD derived from regional CORS network are used to
interpolate ZPD for user site.We also demonstrate that H1QM3model and MLCM model
are equivalent,and propose the height constraint strategy to improve interpolation precision.
It is shown that the H1QX1model and H1QM3model are superior to other models,interpolation
 precision is 10mm in the flat region and 20mm in the undulated region.With height
constraint,the precision of intepolated ZPD can be improved by 60%in the undulated region.
The optimal number of reference stations for intepolation is 9or 10.
Key words:CORS;ZPD;spatial interpolation;augmentation;PPP
About the first author:ZHANG Xiaohong,professor,Ph.D,Ph.D supervisor.He is concentrated on the research and education in GNSS precise
 positioning and its geoscience applications.
E-mail:xhzhang@sgg.whu.edu.cn
683
区域CORS网络增强PPP天顶对流层延迟内插建模.pdf (863.87 KB, 下载次数: 0)
分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友 微信微信
收藏收藏 分享分享 支持支持 反对反对
您需要登录后才可以回帖 登录 | 注册

本版积分规则

小黑屋|文字版|手机版|DIY编程器网 ( 桂ICP备14005565号-1 )

GMT+8, 2024-5-3 23:02 , 耗时 0.095677 秒, 21 个查询请求 , Gzip 开启.

各位嘉宾言论仅代表个人观点,非属DIY编程器网立场。

桂公网安备 45031202000115号

DIY编程器群(超员):41210778 DIY编程器

DIY编程器群1(满员):3044634 DIY编程器1

diy编程器群2:551025008 diy编程器群2

QQ:28000622;Email:libyoufer@sina.com

本站由桂林市临桂区技兴电子商务经营部独家赞助。旨在技术交流,请自觉遵守国家法律法规,一旦发现将做封号删号处理。

快速回复 返回顶部 返回列表