CN104656108A - 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 - Google Patents
一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 Download PDFInfo
- Publication number
- CN104656108A CN104656108A CN201510076735.7A CN201510076735A CN104656108A CN 104656108 A CN104656108 A CN 104656108A CN 201510076735 A CN201510076735 A CN 201510076735A CN 104656108 A CN104656108 A CN 104656108A
- Authority
- CN
- China
- Prior art keywords
- elevation
- reference station
- delay
- zenith tropospheric
- zenith
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000005436 troposphere Substances 0.000 title claims abstract description 14
- 230000000694 effects Effects 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 239000005433 ionosphere Substances 0.000 claims abstract description 4
- 230000002596 correlated effect Effects 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 5
- 230000002706 hydrostatic effect Effects 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 description 3
- 230000002708 enhancing effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001061260 Emmelichthys struhsakeri Species 0.000 description 1
- 244000287680 Garcinia dulcis Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法。常规区域网络增强PPP(Precise Point Position,精密单点定位)建模方法,未考虑高程差异对于内插精度的影响并依赖于一定数量的参考站个数。本发明方法考虑到高程差异与天顶对流层延迟之间的关系,通过无电离层组合精密单点滤波计算得到参考站天顶对流层湿延迟值,并通过模型计算天顶对流层干延迟值,依据对流层延迟值与高程因子之间的关系将干(湿)延迟值划分为高程强相关部分和高程弱相关部分,在此基础上提出一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法。本发明方法考虑到高程差异对于对流层延迟值的影响,适用于稀疏参考站区域网,并且方法对于参考站数量依赖性低。
Description
技术领域
本发明涉及定位与监测领域,尤其涉及一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法,是区域参考站增强PPP(Precise Point Position,精密单点定位)高精度快速定位研究的重要部分。
背景技术
随着卫星定位***的完善与发展,定位与监测目标的精度和可靠性要求越来越高。在高精度定位和导航计算中,通常采用相对定位的方法,但当流动站间的距离较远时,差分定位的方法受到误差的影响,无法快速准确得到精密定位结果。因此对于大范围作业,为了提高定位,增加参考站的个数就无疑增加了工程成本。因此技术产生,并逐渐成为国内外学者的研究热点。
精密单点定位仅采用一台接收机进行数据采集,即节约成本又作业方便自由,且数据处理相对简单。目前,精密单点定位技术利用国际全球卫星定位导航服务组织(International GNSS Service,IGS)提供的精密星历和精密钟差文件,处理一台双频接收机的观测数据,采用非差模型进行精密单点定位。而精密单点定位主要的研究都是基于双频观测值的基础上,已经取得了很高的定位精度。对于PPP用户,在得到高精度卫星轨道、钟差的同时,由于大气延迟等环境因素的制约,仍需要至少30分钟才能达到10cm以内的收敛定位结果,大大限制了PPP技术的应用推广。
利用大气增强信息可显著提高PPP的收敛速度。利用坐标精确已知且观测连续的参考站,可以实时计算得到站点精确的大气延迟信息,将大气延迟信息进行区域内插建模,可以为区域内的PPP用户提供大气延迟误差先验值,PPP用户利用大气误差先验值即可大大缩短定位初始化时间,并提高定位结果精度。
对流层作为大气误差的主要影响因素,因为对流层的气象特征,天顶对流层延迟在高程方向上有强烈的垂直混合,而在水平方向上分布不均匀。测站天顶对流层与测站高程之间存在强烈的相关性,如图2,因此,在区域对流层延迟内插建模时,必须考虑测站高程的影响。
常规的天顶对流层延迟内插模型在网络RTK中研究较多,但基本上是针对网络RTK技术中双差误差改正提出的。在PPP领域,目前的区域增强的PPP天顶对流层误差建模方法大多基于区域连续运行参考***(Continuous Operational Reference System,CORS),对于参考站数量及空间距离有一定要求,对于偏远地区依赖稀疏参考站增强的区域适用性较差。因此需要建立一种适用于稀疏参考站区域的顾及高程差异的对流层延迟内插建模方法。
发明内容
发明目的:为了克服现有技术中存在的不足,提出了一种顾及高程差异的稀疏参考站网络的天顶对流层延迟建模新方法。
技术方案:一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法,包括如下步骤:
(1),基于区域参考站双频原始观测值,利用区域内各参考站点的精确坐标,通过无电离层组合精密单点滤波获得各参考站天顶对流层湿延迟值ZWDr,并通过对流层天顶延迟模型UNB3m获得各参考站天顶对流层干延迟值ZHDr,其中r=1,2,....,n为参考站序号,n为参考站总个数;
(2),根据UNM3m经验模型中对流层干延迟、湿延迟与高程因子的不同关系,将各参考站天顶对流层湿延迟值ZWDr分别分离得到高程强相关天顶对流层湿延迟值ZWDHr和高程弱相关天顶对流层湿延迟值ZWDDr;将各参考站天顶对流层干延迟值ZHDr分别分离得到高程强相关天顶对流层干延迟值ZHDHr和高程弱相关天顶对流层干延迟值ZHDDr;
(3),对于每个参考站,将其高程强相关天顶对流层湿延迟值ZWDHr和高程强相关天顶对流层干延迟值ZHDHr,得到参考站高程强相关天顶对流层延迟值ZTDHr;将其高程弱相关天顶对流层湿延迟值ZWDDr和高程弱相关天顶对流层干延迟值ZHDDr相加,得到参考站高程弱相关天顶对流层延迟值ZTDDr,如式子(1)所示:
ZTDDr=ZHDDr+ZWDDr
(1);
ZTDHr=ZTDHr+ZWDHr
(4),根据用户站与各参考站之间的位置关系,对高程强相关天顶对流层延迟值ZTDHr和高程弱相关天顶对流层延迟值ZTDDr分别建模,得到用户站高程强相关天顶对流层延迟值ZTDHu和高程弱相关天顶对流层延迟值ZTDDu;然后将用户站高程强相关天顶对流层延迟值ZTDHu和高程弱相关天顶对流层延迟值ZTDDu相加,得到用户站天顶对流层延迟值ZTDu。
作为本发明的优选方案,所述步骤(1)包括如下具体步骤:
步骤(11),建立基于无电离层模型的精密单点定位估计方程:
式(2)中,k表示参考站接收机j表示卫星,为卫星j与参考站k的距离,c为光速,为站点接收机k观测得到的卫星j的无电离层组合伪距观测值,为站点接收机k观测得到的卫星j的无电离层组合载波观测值,δtk为参考站接收机k的接收机钟差,δtj表示卫星j的钟差改正数,为站点接收机k与卫星j之间的对流层延迟,hk为接收机k的天线相位改正,为站点接收机k观测得到的卫星j的卫星高度角,为站点接收机k观测得到的卫星j的卫星轨道误差,为站点接收机k的无电离层组合伪距观测值的接收机硬件延迟,为卫星j的无电离层组合伪距观测值的卫星硬件延迟,为站点接收机k的无电离层组合载波观测值的接收机载波初始模糊度偏差,为卫星j的无电离层组合载波观测值的卫星载波初始模糊度偏差,为站点接收机k与卫星j之间无电离层组合伪距观测值的多路径效应影响、为站点接收机k与卫星j之间无电离层组合载波观测值的多路径效应影响,为站点接收机k与卫星j之间其他与频率无关的误差,包括地球自转,相对论效应,潮汐改正等,为站点接收机k与卫星j之间无电离层组合伪距观测值的伪距观测噪声,为站点接收机k与卫星j之间无电离层组合载波观测值的载波观测噪声,λif为无电离层组合波长,Nif为无电离层组合模糊度;
步骤(12),由于参考站坐标精确已知,所述步骤(11)中未知参数为接收机钟差、天顶对流层湿延迟、浮点无电离层模糊度,通过卡尔曼滤波方法对式(2)中的未知参数进行滤波估计后,获得各参考站天顶对流层湿延迟值ZWDr;
步骤(13),根据对流层天顶延迟模型UMB3m得到对流层天顶干延迟、湿延迟与高程H的关系ZHD(H)、ZWD(H):
式(3)、(4)中,P0为参考站大气压气象参数值,T0为参考站温度气象参数值,R为参考站大气相对湿度气象参数值,β为参考站温度变化率气象参数值,λ为参考站水汽压变化率气象参数值,e0为参考站水汽压气象参数值,H为参考站高程,gm为参考站重力加速度,g为标准重力加速度,λ′=λ+1,K1、K2′、K3均为UMB3m模型常系数;
当式(3)中的H为参考站的高程h时,通过式(3)可得到参考站天顶对流层干延迟值ZHDr=ZHD(H=h)。
作为本发明的优选方案,步骤(2)所述参考站对流层干、湿延迟值的分离方法具体为:
根据所述式(3)和式(4)中天顶对流层干、湿延迟值与高程H之间的关系:假设当高程H=0时,根据所述式(3)和式(4)得到天顶对流层干、湿延迟值与高程弱相关分量为:ZHD(H=0)、ZWD(H=0);当H=h时,根据所述式(3)和式(4)得到天顶对流层干、湿延迟值为:ZHD(H=h)、ZWD(H=h);天顶对流层干、湿延迟值与高程强相关分量分别通过天顶对流层干、湿延迟值与天顶对流层干、湿延迟值与高程弱相关分量相减得到为;对于天顶对流层干延迟而言,因为干延迟值的计算通过经验模型得到,故可直接进行加减获得天顶对流层延迟值高程弱相关分量ZHDDr与天顶对流层延迟值高程强相关分量ZHDHr,具体为:
对于某参考站的天顶对流层干延迟值ZHDr,其高程弱相关天顶对流层干延迟值ZHDDr与高程强相关天顶对流层干延迟值ZHDHr分别为:
ZHDDr=ZHD(H=0)
(5)
ZHDHr=ZHD(H=h)-ZHDDr=ZHDr-ZHDDr
对于某参考站的天顶对流层湿延迟值ZWDr,其高程弱相关天顶对流层湿延迟值ZWDDr与高程强相关天顶对流层湿延迟值ZWDHr分别为:
(6)。
ZWDHr=ZWDr-ZWDDr
作为本发明的优选方案,所述步骤(4)包括如下具体步骤:
步骤(41),对高程弱相关天顶对流层延迟值ZTDDr(r=1,2,....,n)内插建模,具体为:
令ZTDDr=a1+a2xr+a3yr,其中a1、a2、a3均为内插系数,内插时以用户站为中心,建立局部坐标系,(xr,yr)为局部坐标系下参考站平面坐标,用户站的(xu,yu)为(0,0);因此,用户站高程弱相关天顶对流层延迟值ZTDDu为:
ZTDDu=a1, (7)
内插系数ai(1=1.2.3)满足式(13):
式(8)中, 其中(xi,yi)i=1.2.3...n为局部坐标系下参考站平面坐标,ZTDDr(r=1,2,....,n)为各参考站高程弱相关天顶对流层延迟值;
步骤(42),对高程强相关天顶对流层延迟值ZTDHr(r=1,2,....,n)分别内插建模,具体为:
令ZTDHu=b1ZTDH1+b2ZTDH2…+brZTDHr…+bnZTDHn,其中,bi(i=1.2.3...n)为内插系数,满足hi(i=1.2.3...n)为局部坐标系下各参考站高程,n为参考站数量;内插时以用户站为中心,建立局部坐标系;
内插系数bi(i=1.2.3...n)解算满下式:
式(14)中,
步骤(43),根据步骤(41)得到的用户站高程弱相关天顶对流层延迟值ZTDDu和步骤(42)得到的用户站高程强相关天顶对流层延迟值ZTDHu相加得到用户站的天顶对流层延迟值ZTDu:
ZTDu=ZTDHu+ZTDDu (10)。
有益效果:本发明提出的顾及高程差异的稀疏参考站天顶对流层延迟内插建模方法,利用UNB3m经验模型中高程因子和测站天顶干、湿延迟之间的关系,对测站天顶对流层延迟进行分离,得到高程强相关分量和高程弱相关分量。分别对两个分量进行区域内插建模,得到用户站的天顶对流层延迟。内插方法简便,计算量小,且适用性广泛,利用稀疏参考站区域少量参考站即可进行区域内插建模,为用户站定位提供增强的大气先验信息,加快用户站定位的收敛时间,有利于PPP级数的应用推广。
附图说明
图1是本发明方法流程图;
图2是测站天顶对流层与测站高程之间相关性;
图3是实施例测站站点网分布图;
图4是利用3个参考站建模时模型内插值与真值偏差值变化图;
图5是利用4个参考站建模时模型内插值与真值偏差值变化图;
图6是利用5个参考站建模时模型内插值与真值偏差值变化图;
图7是利用6个参考站建模时模型内插值与真值偏差值变化图;
图8是利用7个参考站建模时模型内插值与真值偏差值变化图;
图9是不同参考站数量下建模结果绝对偏差与中误差RMS统计图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
如图1所示,一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法,包括如下步骤:
(1),基于区域参考站双频原始观测值,利用区域内各参考站点的精确坐标,通过无电离层组合精密单点滤波获得各参考站天顶对流层湿延迟值ZWDr,并通过对流层天顶延迟模型UNB3m获得各参考站天顶对流层干延迟值ZHDr,其中r=1,2,....,n为参考站序号,n为参考站总个数;具体步骤为:
步骤(11),建立基于无电离层模型的精密单点定位估计方程:
式(1)中,k表示参考站接收机j表示卫星,为卫星j与参考站k的距离,c为光速,为站点接收机k观测得到的卫星j的无电离层组合伪距观测值,为站点接收机k观测得到的卫星j的无电离层组合载波观测值,δtk为参考站接收机k的接收机钟差,δtj表示卫星j的钟差改正数,为站点接收机k与卫星j之间的对流层延迟,hk为接收机k的天线相位改正,为站点接收机k观测得到的卫星j的卫星高度角,为站点接收机k观测得到的卫星j的卫星轨道误差,为站点接收机k的无电离层组合伪距观测值的接收机硬件延迟,为卫星j的无电离层组合伪距观测值的卫星硬件延迟,为站点接收机k的无电离层组合载波观测值的接收机载波初始模糊度偏差,为卫星j的无电离层组合载波观测值的卫星载波初始模糊度偏差,为站点接收机k与卫星j之间无电离层组合伪距观测值的多路径效应影响、为站点接收机k与卫星j之间无电离层组合载波观测值的多路径效应影响,为站点接收机k与卫星j之间其他与频率无关的误差,包括地球自转,相对论效应,潮汐改正等,为站点接收机k与卫星j之间无电离层组合伪距观测值的伪距观测噪声,为站点接收机k与卫星j之间无电离层组合载波观测值的载波观测噪声,λif为无电离层组合波长,Nif为无电离层组合模糊度;
步骤(12),由于参考站坐标精确已知,步骤(11)中未知参数为接收机钟差、天顶对流层湿延迟、浮点无电离层模糊度,通过卡尔曼滤波方法对式(1)中的未知参数进行滤波估计后,获得各参考站天顶对流层湿延迟值ZWDr;
步骤(13),以参考站的大地纬度和测量时间为依据,按照对流层天顶延迟模型UMB3m中对流层天顶方向干、湿延迟的气象参数格网值进行内插求得UNB3m的气象参数值。其中,气象参数值包括大气压P0(mbar)、温度T0(K)、大气相对湿度RH、温度变化率β(K/m)和水汽压变化率λ(mbar/m),内插公式如式(2)所示:
式(2)中,φ为参考站处的大地纬度,φi、φi+1为与φ相差最近的格网大地纬度,m=(φ-φi)/(φi+1-φi),t为年积日,ξ为各气象参数的内插值,ξavg(·)为气象参数年度平均值,ξamp(·)为气象参数年度振幅;
步骤(14),根据解算得到的大气相对湿度RH,得到对流层天顶延迟模型UMB3m中的水汽压e0(mbar)气象参数值:
其中,fw=1.00062+3.14×10-6P0+5.6×10-7(T0-273.15)2 (4)
步骤(15),根据上述气象参数值,得到对流层天顶延迟模型UMB3m的干延迟、湿延迟与高程H的关系ZHD(H)、ZWD(H):
式(6)、(7)中,K1=77.60k·mbar-1,R=287.054J·kg-1·K-1,H为高程(单位m),gm=9.784(1-2.66×10-3cos(2φ)-2.8×10-7H)(单位mn/s2),g=9.80665m/s2,K2′=16.6k·mbar-1,K3=377600k·mbar-1,λ′=λ+1;
当式(6)中的H为参考站的高程h时,通过式(6)可得到参考站天顶对流层干延迟值ZHDr=ZHD(H=h)。
式(7)的湿延迟与高程H的关系不能用来直接计算对流层湿延迟值,这是因为其精度不够,对流层湿延迟量是一个时变化量,不能通过经验模型精确得到,必须通过步骤(12)的滤波实时计算。而对流程干延迟是一个相对时间变化较慢的量,所以可以通过经验模型计算得到,并且模型可以保证结果需要的精度。
(2),根据UNM3m经验模型中对流层干延迟、湿延迟与高程因子的不同关系,将各参考站天顶对流层湿延迟值ZWDr分别分离得到高程强相关天顶对流层湿延迟值ZWDHr和高程弱相关天顶对流层湿延迟值ZWDDr;将各参考站天顶对流层干延迟值ZHDr分别分离得到高程强相关天顶对流层干延迟值ZHDHr和高程弱相关天顶对流层干延迟值ZHDDr;具体步骤如下:
根据式(6)和式(7)中天顶对流层干、湿延迟值与高程H之间的关系:不难发现,参考站高程和大地纬度对于参考站重力加速度gm的影响很小,因此在计算时可以近似取值gm=9.784,在此基础上,当参考站大地纬度确定后,式(6)和式(7)中的自变量仅有高程H,其余参数全部确定,并且考虑到高程的基准面的高程值为0,因此假设当高程H=0时,根据式(6)和式(7)得到对流层天顶延迟模型UMB3m下的天顶对流层干、湿延迟值与高程弱相关分量为:ZHD(H=0)、ZWD(H=0);当H=h时,根据式(6)和式(7)可得到天顶对流层干、湿延迟值为:ZHD(H=h)、ZWD(H=h);天顶对流层干、湿延迟值与高程强相关分量分别通过天顶对流层干、湿延迟值与天顶对流层干、湿延迟值与高程弱相关分量相减得到,具体为:
对于某参考站的天顶对流层干延迟值ZHDr,因为干延迟值的计算是通过经验模型UMB3m得到,故可直接进行加减获得天顶对流层延迟值高程弱相关分量ZHDDr与天顶对流层延迟值高程强相关分量ZHDHr为:
ZHDDr=ZHD(H=0)
(9)
ZHDHr=ZHD(H=h)-ZHDDr=ZHDr-ZHDDr
对于某参考站的天顶对流层湿延迟ZWDDr,因为湿延迟的计算通过滤波参数估计得到,故无法直接进行加减,需利用经验模型进行比例分配得到天顶对流层延迟值高程弱相关分量ZWDDr,然后进行加减得到天顶对流层延迟值高程强相关分量ZWDHr:
(10)
ZWDHr=ZWDr-ZWDDr
需要说明的是,对于对流层干延迟来说,它的计算获得自始至终都是通过经验模型UNM3m得到的,首先通过经验模型得到总的天顶对流层干延迟,在分离对流层干延迟的时候,只能通过经验模型计算得到与高程弱相关部分ZHDDr,然后通过减法得到干延迟强相关量ZHDHr。步骤(12)通过滤波得到的精确的实时变化的对流程湿延迟值,但没有直接的方法对湿延迟值进行分离,只能参照经验模型中湿延迟强相关量与弱相关量之间比例关系,将步骤(12)得到的湿延迟值进行一个按比例分配。
(3),对于每个参考站,将其高程强相关天顶对流层湿延迟值ZWDHr和高程强相关天顶对流层干延迟值ZHDHr,得到参考站高程强相关天顶对流层延迟值ZTDHr;将其高程弱相关天顶对流层湿延迟值ZWDDr和高程弱相关天顶对流层干延迟值ZHDDr相加,得到参考站高程弱相关天顶对流层延迟值ZTDDr,如式子(11)所示:
ZTDDr=ZHDDr+ZWDDr
(11)。
ZTDHr=ZTDHr+ZWDHr
(4),根据用户站与各参考站之间的位置关系,对高程强相关天顶对流层延迟值ZTDHr和高程弱相关天顶对流层延迟值ZTDDr分别建模,得到用户站高程强相关天顶对流层延迟值ZTDHu和高程弱相关天顶对流层延迟值ZTDDu;然后将用户站高程强相关天顶对流层延迟值ZTDHu和高程弱相关天顶对流层延迟值ZTDDu相加,得到用户站天顶对流层延迟值ZTDu;具体步骤如下:
步骤(41),对高程弱相关天顶对流层延迟值ZTDDr(r=1,2,....,n)内插建模,具体为:
令ZTDDr=a1+a2xr+a3yr,其中a1、a2、a3均为内插系数,内插时以用户站为中心,建立局部坐标系,(xr,yr)为局部坐标系下参考站平面坐标,用户站的(xu,yu)为(0,0);因此,用户站高程弱相关天顶对流层延迟值ZTDDu为:
ZTDDu=a1, (12)
内插系数ai(1=1.2.3)满足式(13):
式(13)中, 其中(xi,yi)i=1.2.3...n为局部坐标系下参考站平面坐标,ZTDDr(r=1,2,....,n)为各参考站高程弱相关天顶对流层延迟值;
步骤(42),对高程强相关天顶对流层延迟值ZTDHr(r=1,2,....,n)分别内插建模,具体为:
令ZTDHu=b1ZTDH1+b2ZTDH2…+brZTDHr…+bnZTDHn,其中,bi(i=1.2.3...n)为内插系数,满足hi(i=1.2.3...n)为局部坐标系下各参考站高程,n为参考站数量;内插时以用户站为中心,建立局部坐标系;
内插系数bi(i=1.2.3...n)解算满下式:
式(14)中,
步骤(43),根据步骤(41)得到的用户站高程弱相关天顶对流层延迟值ZTDDu和步骤(42)得到的用户站高程强相关天顶对流层延迟值ZTDHu相加得到用户站的天顶对流层延迟值ZTDu:
ZTDu=ZTDHu+ZTDDu (15)
用户站利用内插建模得到的天顶对流层延迟值,可以用于用户站的精密单点定位大气先验信息,提高自身定位收敛速度。
本实例采用美国CORS的参考站组网进行试验。选择美国CORS网北纬34°~36°、西经118°~121°区域作为实验区域,选取八个站点进行实验,中间站作为用户站(三角形),其余站作为参考站(圆圈),如图3所示。参考站高程最大值为1526.119m,最小值为56.588m,用户站高程为1709.067m,用户站与参考站之间平均距离为79.25km。实验数据为2013年8月8日24小时的观测数据,采样率为15秒。
选取站点PPP计算时,模糊度收敛后的2000个历元的ZTD值采用本发明方法进行实时区域建模,用户站PPP计算得到的ZTD值作为真值。分别选取不同数量的参考站进行区域内插建模,将内插建模结果与真值进行比较,结果如图4-图8。并针对2000个历元数据统计不同内插建模情况的绝对偏差及中误差RMS值,如图9所示。
从图4-图9也不难发现本发明方法在高程差异较大区域的稀疏参考站情况下,本发明方法建模结果随着参考站数量的增加,会提高模型的内插精度,特别是绝对偏差的精度会相对提高较大。并且本发明方法利用三个参考数量即可进行内插建模,三个参考站的内插建模精度就已经达到了对流层延迟建模的需要,绝对偏差小于0.03m就达到93.85%,基本保证小于0.03m,并且内插中误差也保证在0.02m以内。本发明方法考虑到对流层建模时高程因素对其的影响,采用较为简便的方法进行内插建模。建模结果受参考站数量影响较弱,且内插精度较高,可作为大气先验信息提供给区域内的PPP用户使用。对于参考站稀疏且高程变化剧烈地区的天顶对流层延迟建模有明显参考应用价值。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法,其特征在于,包括如下步骤:
(1),基于区域参考站双频原始观测值,利用区域内各参考站点的精确坐标,通过无电离层组合精密单点滤波获得各参考站天顶对流层湿延迟值ZWDr,并通过对流层天顶延迟模型UNB3m获得各参考站天顶对流层干延迟值ZHDr,其中r=1,2,....,n为参考站序号,n为参考站总个数;
(2),根据UNM3m经验模型中对流层干延迟、湿延迟与高程因子的不同关系,将各参考站天顶对流层湿延迟值ZWDr分别分离得到高程强相关天顶对流层湿延迟值ZWDHr和高程弱相关天顶对流层湿延迟值ZWDDr;将各参考站天顶对流层干延迟值ZHDr分别分离得到高程强相关天顶对流层干延迟值ZHDHr和高程弱相关天顶对流层干延迟值ZHDDr;
(3),对于每个参考站,将其高程强相关天顶对流层湿延迟值ZWDHr和高程强相关天顶对流层干延迟值ZHDHr,得到参考站高程强相关天顶对流层延迟值ZTDHr;将其高程弱相关天顶对流层湿延迟值ZWDDr和高程弱相关天顶对流层干延迟值ZHDDr相加,得到参考站高程弱相关天顶对流层延迟值ZTDDr,如式子(1)所示:
(4),根据用户站与各参考站之间的位置关系,对高程强相关天顶对流层延迟值ZTDHr和高程弱相关天顶对流层延迟值ZTDDr分别建模,得到用户站高程强相关天顶对流层延迟值ZTDHu和高程弱相关天顶对流层延迟值ZTDDu;然后将用户站高程强相关天顶对流层延迟值ZTDHu和高程弱相关天顶对流层延迟值ZTDDu相加,得到用户站天顶对流层延迟值ZTDu。
2.根据权利要求1所述的一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法,其特征在于,所述步骤(1)包括如下具体步骤:
步骤(11),建立基于无电离层模型的精密单点定位估计方程:
式(1)中,k表示参考站接收机j表示卫星,为卫星j与参考站k的距离,c为光速,为站点接收机k观测得到的卫星j的无电离层组合伪距观测值,为站点接收机k观测得到的卫星j的无电离层组合载波观测值,δtk为参考站接收机k的接收机钟差,δtj表示卫星j的钟差改正数,为站点接收机k与卫星j之间的对流层延迟,hk为接收机k的天线相位改正,为站点接收机k观测得到的卫星j的卫星高度角,为站点接收机k观测得到的卫星j的卫星轨道误差,为站点接收机k的无电离层组合伪距观测值的接收机硬件延迟,为卫星j的无电离层组合伪距观测值的卫星硬件延迟, 为站点接收机k的无电离层组合载波观测值的接收机载波初始模糊度偏差,为卫星j的无电离层组合载波观测值的卫星载波初始模糊度偏差,为站点接收机k与卫星j之间无电离层组合伪距观测值的多路径效应影响、为站点接收机k与卫星j之间无电离层组合载波观测值的多路径效应影响,为站点接收机k与卫星j之间其他与频率无关的误差,包括地球自转,相对论效应,潮汐改正等,为站点接收机k与卫星j之间无电离层组合伪距观测值的伪距观测噪声,为站点接收机k与卫星j之间无电离层组合载波观测值的载波观测噪声,λif为无电离层组合波长,Nif为无电离层组合模糊度;
步骤(12),由于参考站坐标精确已知,所述步骤(11)中未知参数为接收机钟差、天顶对流层湿延迟、浮点无电离层模糊度,通过卡尔曼滤波方法对式(2)中的未知参数进行滤波估计后,获得各参考站天顶对流层湿延迟值ZWDr;
步骤(13),根据对流层天顶延迟模型UMB3m得到对流层天顶干延迟、湿延迟与高程H的关系ZHD(H)、ZWD(H):
式(3)、(4)中,P0为参考站大气压气象参数值,T0为参考站温度气象参数值,R为参考站大气相对湿度气象参数值,β为参考站温度变化率气象参数值,λ为参考站水汽压变化率气象参数值,e0为参考站水汽压气象参数值,H为参考站高程,gm为参考站重力加速度,g为标准重力加速度,λ′=λ+1,K1、K2′、K3均为UMB3m模型常系数;
当式(3)中的H为参考站的高程h时,通过式(3)可得到参考站天顶对流层干延迟值ZHDr=ZHD(H=h)。
3.根据权利要求2所述的一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法:步骤(2)所述参考站对流层干、湿延迟值的分离方法具体为:
根据所述式(3)和式(4)中天顶对流层干、湿延迟值与高程H之间的关系:假设当高程H=0时,根据所述式(3)和式(4)得到天顶对流层干、湿延迟值与高程弱相关分量为:ZHD(H=0)、ZWD(H=0);当H=h时,根据所述式(3)和式(4)得到天顶对流层干、湿延迟值为:ZHD(H=h)、ZWD(H=h);天顶对流层干、湿延迟值与高程强相关分量分别通过天顶对流层干、湿延迟值与天顶对流层干、湿延迟值与高程弱相关分量相减得到为;对于天顶对流层干延迟而言,因为干延迟值的计算通过经验模型得到,故可直接进行加减获得天顶对流层延迟值高程弱相关分量ZHDDr与天顶对流层延迟值高程强相关分量ZHDHr,具体为:
对于某参考站的天顶对流层干延迟值ZHDr,其高程弱相关天顶对流层干延迟值ZHDDr与高程强相关天顶对流层干延迟值ZHDHr分别为:
对于某参考站的天顶对流层湿延迟值ZWDr,其高程弱相关天顶对流层湿延迟值ZWDDr与高程强相关天顶对流层湿延迟值ZWDHr分别为:
ZWDHr=ZWDr-ZWDDr。
4.根据权利要求3所述的一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法,其特征在于:所述步骤(4)包括如下具体步骤:
步骤(41),对高程弱相关天顶对流层延迟值ZTDDr(r=1,2,....,n)内插建模,具体为:
令ZTDDr=a1+a2xr+a3yr,其中a1、a2、a3均为内插系数,内插时以用户站为中心,建立局部坐标系,(xr,yr)为局部坐标系下参考站平面坐标,用户站的(xu,yu)为(0,0);因此,用户站高程弱相关天顶对流层延迟值ZTDDu为:
ZTDDu=a1, (7)
内插系数ai(1=1.2.3)满足式(8):
式(8)中,其中(xi,yi)i=1.2.3...n为局部坐标系下参考站平面坐标,ZTDDr(r=1,2,....,n)为各参考站高程弱相关天顶对流层延迟值;
步骤(42),对高程强相关天顶对流层延迟值ZTDHr(r=1,2,....,n)分别内插建模,具体为:
令ZTDHu=b1ZTDH1+b2ZTDH2…+brZTDHr…+bnZTDHn,其中,bi(i=1.2.3...n)为内插 系数,满足为局部坐标系下各参考站高程,n为参考站数量;内插时以用户站为中心,建立局部坐标系;
内插系数bi(i=1.2.3...n)解算满下式:
式(9)中,
步骤(43),根据步骤(41)得到的用户站高程弱相关天顶对流层延迟值ZTDDu和步骤(42)得到的用户站高程强相关天顶对流层延迟值ZTDHu相加得到用户站的天顶对流层延迟值ZTDu:
ZTDu=ZTDHu+ZTDDu (10)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510076735.7A CN104656108B (zh) | 2015-02-12 | 2015-02-12 | 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510076735.7A CN104656108B (zh) | 2015-02-12 | 2015-02-12 | 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104656108A true CN104656108A (zh) | 2015-05-27 |
CN104656108B CN104656108B (zh) | 2017-05-10 |
Family
ID=53247460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510076735.7A Expired - Fee Related CN104656108B (zh) | 2015-02-12 | 2015-02-12 | 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104656108B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105785407A (zh) * | 2016-02-23 | 2016-07-20 | 东南大学 | 一种适用于中国地区的无气象参数对流层延迟改正方法 |
CN106814373A (zh) * | 2017-01-22 | 2017-06-09 | 武汉大学 | 大气加权平均温度估算及对流层延迟积分方法 |
CN107037470A (zh) * | 2017-05-04 | 2017-08-11 | 南京康帕斯导航科技有限公司 | 基于逆向rtd和行人航位推算融合的米级定位方法及*** |
CN109782275A (zh) * | 2019-03-14 | 2019-05-21 | 中国电建集团成都勘测设计研究院有限公司 | Gnss变形监测的参考点校验***及方法 |
CN110488332A (zh) * | 2018-05-14 | 2019-11-22 | 广州市中海达测绘仪器有限公司 | 基于网络rtk技术的定位信息处理方法和装置 |
CN110618438A (zh) * | 2019-09-09 | 2019-12-27 | 广州市中海达测绘仪器有限公司 | 大气误差计算方法、装置、计算机设备和存储介质 |
CN111896977A (zh) * | 2019-05-06 | 2020-11-06 | 千寻位置网络有限公司 | 对流层湿延迟精度计算方法及其***和定位方法及其*** |
CN112987058A (zh) * | 2021-05-20 | 2021-06-18 | 长江空间信息技术工程有限公司(武汉) | 一种利用地表气象站增强短距离大高差rtk定位的方法 |
CN113093241A (zh) * | 2021-03-12 | 2021-07-09 | 东南大学 | 一种考虑高度角的单测站对流层斜延迟计算方法 |
CN113267793A (zh) * | 2021-05-26 | 2021-08-17 | 中国电子科技集团公司第五十四研究所 | 一种基于外部增强信息的gbas对流层参数生成方法 |
CN113325448A (zh) * | 2021-04-21 | 2021-08-31 | 中铁第一勘察设计院集团有限公司 | 顾及对流层延迟改造的大高差cors网解算方法 |
CN114527497A (zh) * | 2022-02-08 | 2022-05-24 | 国汽大有时空科技(安庆)有限公司 | 一种基于ppp-rtk的定位增强信息传输方法 |
CN115144878A (zh) * | 2022-07-13 | 2022-10-04 | 武汉大学 | 基于ppp的短距离大高差nrtk对流层延迟改正方法 |
CN116626730A (zh) * | 2023-07-24 | 2023-08-22 | 山东科技大学 | 一种顾及nwp的海上区域cors增强ppp方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027624A1 (en) * | 2003-04-17 | 2007-02-01 | Secretary Of State For Defence | Correction of troposhere induced errors in global positioning systems |
CN101943749A (zh) * | 2010-09-10 | 2011-01-12 | 东南大学 | 基于星型结构的虚拟参考站网络rtk定位方法 |
CN102298151A (zh) * | 2011-07-20 | 2011-12-28 | 东南大学 | 一种gnss网络差分定位***中的误差改正方法 |
-
2015
- 2015-02-12 CN CN201510076735.7A patent/CN104656108B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070027624A1 (en) * | 2003-04-17 | 2007-02-01 | Secretary Of State For Defence | Correction of troposhere induced errors in global positioning systems |
CN101943749A (zh) * | 2010-09-10 | 2011-01-12 | 东南大学 | 基于星型结构的虚拟参考站网络rtk定位方法 |
CN102298151A (zh) * | 2011-07-20 | 2011-12-28 | 东南大学 | 一种gnss网络差分定位***中的误差改正方法 |
Non-Patent Citations (1)
Title |
---|
杨徉等: "一种综合的对流层延迟模型算法", 《东南大学学报(自然科学版)》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105785407B (zh) * | 2016-02-23 | 2017-12-22 | 东南大学 | 一种适用于中国地区的无气象参数对流层延迟改正方法 |
CN105785407A (zh) * | 2016-02-23 | 2016-07-20 | 东南大学 | 一种适用于中国地区的无气象参数对流层延迟改正方法 |
CN106814373A (zh) * | 2017-01-22 | 2017-06-09 | 武汉大学 | 大气加权平均温度估算及对流层延迟积分方法 |
CN106814373B (zh) * | 2017-01-22 | 2019-09-10 | 武汉大学 | 大气加权平均温度估算及对流层延迟积分方法 |
CN107037470A (zh) * | 2017-05-04 | 2017-08-11 | 南京康帕斯导航科技有限公司 | 基于逆向rtd和行人航位推算融合的米级定位方法及*** |
CN110488332B (zh) * | 2018-05-14 | 2021-09-10 | 广州市中海达测绘仪器有限公司 | 基于网络rtk技术的定位信息处理方法和装置 |
CN110488332A (zh) * | 2018-05-14 | 2019-11-22 | 广州市中海达测绘仪器有限公司 | 基于网络rtk技术的定位信息处理方法和装置 |
CN109782275A (zh) * | 2019-03-14 | 2019-05-21 | 中国电建集团成都勘测设计研究院有限公司 | Gnss变形监测的参考点校验***及方法 |
CN109782275B (zh) * | 2019-03-14 | 2023-04-28 | 中国电建集团成都勘测设计研究院有限公司 | Gnss变形监测的参考点校验***及方法 |
CN111896977A (zh) * | 2019-05-06 | 2020-11-06 | 千寻位置网络有限公司 | 对流层湿延迟精度计算方法及其***和定位方法及其*** |
CN111896977B (zh) * | 2019-05-06 | 2024-03-29 | 千寻位置网络有限公司 | 对流层湿延迟精度计算方法及其***和定位方法及其*** |
CN110618438A (zh) * | 2019-09-09 | 2019-12-27 | 广州市中海达测绘仪器有限公司 | 大气误差计算方法、装置、计算机设备和存储介质 |
CN113093241A (zh) * | 2021-03-12 | 2021-07-09 | 东南大学 | 一种考虑高度角的单测站对流层斜延迟计算方法 |
CN113093241B (zh) * | 2021-03-12 | 2024-02-02 | 东南大学 | 一种考虑高度角的单测站对流层斜延迟计算方法 |
CN113325448A (zh) * | 2021-04-21 | 2021-08-31 | 中铁第一勘察设计院集团有限公司 | 顾及对流层延迟改造的大高差cors网解算方法 |
CN112987058A (zh) * | 2021-05-20 | 2021-06-18 | 长江空间信息技术工程有限公司(武汉) | 一种利用地表气象站增强短距离大高差rtk定位的方法 |
CN113267793A (zh) * | 2021-05-26 | 2021-08-17 | 中国电子科技集团公司第五十四研究所 | 一种基于外部增强信息的gbas对流层参数生成方法 |
CN114527497A (zh) * | 2022-02-08 | 2022-05-24 | 国汽大有时空科技(安庆)有限公司 | 一种基于ppp-rtk的定位增强信息传输方法 |
CN115144878A (zh) * | 2022-07-13 | 2022-10-04 | 武汉大学 | 基于ppp的短距离大高差nrtk对流层延迟改正方法 |
CN116626730A (zh) * | 2023-07-24 | 2023-08-22 | 山东科技大学 | 一种顾及nwp的海上区域cors增强ppp方法 |
CN116626730B (zh) * | 2023-07-24 | 2023-10-10 | 山东科技大学 | 一种顾及nwp的海上区域cors增强ppp方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104656108B (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104656108A (zh) | 一种顾及高程差异的稀疏参考站网络天顶对流层延迟建模方法 | |
CN108415049B (zh) | 提高网络rtk双差宽巷模糊度固定正确率的方法 | |
EP3130943B1 (en) | Navigation satellite system positioning involving the generation of tropospheric correction information | |
CN104597471B (zh) | 面向时钟同步多天线gnss接收机的定向测姿方法 | |
CN102331583B (zh) | 利用模糊度固定的gnss大气估计 | |
Li et al. | Determination of the differential code bias for current BDS satellites | |
Lu et al. | Improving BeiDou real-time precise point positioning with numerical weather models | |
Pan et al. | Integrating multi-GNSS to improve the performance of precise point positioning | |
CN107942346B (zh) | 一种高精度gnss电离层tec观测值提取方法 | |
CN104483691B (zh) | 一种gnss组合精密单点定位方法 | |
CN110007326B (zh) | 一种用于星基增强***的双频测距误差参数生成方法 | |
CN111596315B (zh) | 一种用于实时监测双频多星座星基增强***性能的方法 | |
CN104316943B (zh) | 一种伪距离和多普勒组合差分定位***及方法 | |
CN105510945A (zh) | 一种应用于卫导着陆外场检测的ppp定位方法 | |
CN111796309A (zh) | 导航卫星单频数据同步确定大气水汽和总电子含量的方法 | |
Godah et al. | COMPARISON OF VERTICAL DEFORMATIONS OF THE EARTH’S SURFACE OBTAINED USING GRACE-BASED GGMS AND GNSS DATA–A CASE STUDY OF SOUTH-EASTERN POLAND. | |
Xu et al. | Influence of meteorological data and horizontal gradient of tropospheric model on precise point positioning | |
CN108459334A (zh) | 一种顾及***间偏差的gps/bds双***单钟差定位方法 | |
Wielgosz et al. | Research on GNSS positioning and applications in Poland in 2015–2018 | |
Li et al. | Inter-system biases solution strategies in multi-GNSS kinematic precise point positioning | |
CN114994724A (zh) | 一种gnss伪距差分定位性能评估方法及*** | |
Abdelazeem | AFRC-trop: new real-time zenith tropospheric delay model over Africa | |
CN103389502B (zh) | 基于多个地面基站高精度确定载体加速度的方法 | |
Abdallah et al. | Accuracy Assessment Study of GNSS Precise Point Positioning for Kinematic Positioning | |
Gao et al. | Real-time sensing atmospheric water vapor using precise GPS orbit and clock products |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170510 |