CN107421434B - 多基站Multi-GNSS长基线近实时形变监测方法 - Google Patents
多基站Multi-GNSS长基线近实时形变监测方法 Download PDFInfo
- Publication number
- CN107421434B CN107421434B CN201710672999.8A CN201710672999A CN107421434B CN 107421434 B CN107421434 B CN 107421434B CN 201710672999 A CN201710672999 A CN 201710672999A CN 107421434 B CN107421434 B CN 107421434B
- Authority
- CN
- China
- Prior art keywords
- ambiguity
- station
- processing
- gnss
- reference station
- 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.)
- Active
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000012545 processing Methods 0.000 claims abstract description 89
- 238000012216 screening Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 13
- 239000005436 troposphere Substances 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000007499 fusion processing Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000012958 reprocessing Methods 0.000 abstract 1
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供了一种基于多基站的Multi‑GNSS长基线近实时形变监测方法,为了提高长基线解算精度,本发明利用超快速Multi‑GNSS卫星轨道产品减少轨道误差对基线解算的影响;利用多次后验残差编辑迭代,尽量提出数据异常对基线解算的影响;利用滑动窗口批处理模式,可以实现近实时高精度基线解算;利用多基准站联合处理模式,可以提高基线解算的精度与稳定性;采用创新的模糊度处理策略,进一步提高基线解算精度。为了提高数据处理效率,本发明利用法方程叠加技术,充分利用以前法方程信息,规避解算重复性;利用共用基准站处理,规避相同基准站的重复处理;利用非差处理模式,可以充分利用分布式计算优势,提高数据处理效率。
Description
技术领域
本发明涉及形变监测技术领域,具体涉及一种基于多基站的Multi-GNSS长基线近实时形变监测方法。
背景技术
GNSS(Global navigation satellite system,全球卫星定位***)形变监测是关乎人民生命财产安全的重要技术,在危房监测、滑坡监测、大坝监测等工况有大量应用。传统形变监测只是为了解决具体业务开展的,随着国民经济的发展,大范围形变监测业务需求日趋明显,如全国范围的危房监测、长距离铁路沿线滑坡监测、长距离电力铁塔沉降监测等。
GNSS形变监测需要布设大量基准站,支持全国范围的毫米形变监测服务,传统近实时监测服务基线长度在15公里左右。为满足全国范围的危房监测、长距离铁路沿线滑坡监测、长距离电力铁塔沉降监测等大区域跨度的近实时形变监测业务需求,目前传统的中短基线长度的算法不能满足大范围的业务需求,并且传统处理算法对于上万监测点的工况,处理效率不高。此外,传统形变监测是为具体需求开展的业务,用户量较少。
发明内容
本发明提出了利用长基线进行近实时毫米级形变监测,为了提高近实时长基线处理精度,从后验残差编辑迭代、滑动窗口处理、多基准站联合处理、模糊度固定策略处理方面进行了改进;为了提高海量用户高效率处理,除了分布式处理以外,还从非差模式并行处理、共用基准站处理、法方程叠加方面进行了优化。
本发明解决了以下技术问题:
1)可以有效减少国内基准站部署数量,充分利用现有基准站资源。
2)并提出了一种高效率数据处理策略,包括短弧法方程合并与非差分布式处理模式可以有效减少计算资源,提供全国范围内近实时形变监测服务,并降低服务成本。
本发明采用的技术方案是:
一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,包括以下步骤:
根据监测点坐标,结合基准站分析,利用基准站选择策略,给出参与解算的基准站列表;
根据参与解算的监测点数据进行子任务划分,执行云平台分布式处理;
云平台分布式处理利用超快速Multi-GNSS卫星轨道产品减少轨道误差对Multi-GNSS近实时长基线解算的影响。
进一步地,所述基准站选择策略具体包括以下步骤:
步骤1,根据监测点坐标,筛选距离100km以内的基准站;
步骤2,根据基准站坐标稳定信息,筛选绝对坐标点位重复率小于水平3mm、垂直6mm的基准站;
步骤3,根据基准站数据质量信息,筛选满足准入标准的基准站;
步骤4,根据方位角,将候选的基准站划分到六个区,在每个区选择距离监测点最近的基准站;
步骤5,给出参与解算的基准站列表。
进一步地,候选的基准站按照以下方法划分到六个区:
1区,0°~60°;
2区,60°~120°;
3区,120°~180°;
4区,180°~240°;
5区,240°~300°;
6区,300°~360°。
进一步地,所述Multi-GNSS近实时长基线解算具体包括以下步骤:
利用超快速Multi-GNSS卫星轨道产品,对基准站和监测点进行非差PPP处理,并进行后验残差编辑迭代,剔除异常数据,生成当前弧段的法方程信息,当前弧段的法方程之前有m-1个法方程;
利用滑动窗口批处理,将当前弧段的法方程与之前m-1个法方程合并成长弧法方程;
采用模糊度固定策略,多基准站采用独立基线对,基线对上卫星对也采用独立卫星对,对基线对上的卫星对的宽巷/窄巷模糊度固定概率进行筛选;
采用宽巷与窄巷组合模式进行模糊度固定,将模糊度固定信息加入到长弧法方程中更新长弧法方程,重复进行宽巷/窄巷模糊度筛选与模糊度固定处理。
进一步地,对基准站和监测点进行非差PPP处理需要估计的参数包括包括:卫星钟差、接收机钟差、坐标、天顶对流层延迟、卫星***间延迟、卫星频率间延迟。
进一步地,所述模糊度固定策略采用双差模糊度约束加入长弧法方程,实现模糊度固定与多基准站联合处理。
进一步地,所述双差模糊度约束等效于对四个非差模糊度加上约束条件,定义双差模糊度的虚拟观测值VB:
其中,为固定的无电离层组合的双差模糊度,PB为双差模糊度虚拟观测值的权值,D是非差模糊度的投影矩阵,X是待估参数向量;
其中,和代表两颗卫星,代表两个接收机,为测站与卫星之间的非差模糊度,为测站与卫星之间的非差模糊度,为测站与卫星之间的非差模糊度,为测站与卫星之间的非差模糊度;
基于公式(11)和公式(12)的约束,最小二乘解的形式改写为公式(13),实现从非差到双差基线处理的转换:
其中,A是系数矩阵,P是权值阵,L是观测值减去模型计算值,NX=W,N是法方程,Xnew为加了双插模糊度约束的X,Nnew为加了双插模糊度约束的法N,Wnew为加了双插模糊度约束的W。
进一步地,如果施加多条基线的双差模糊度约束,即可实现多基准站联合融合处理。
进一步地,在模糊度固定处理中模糊度需要达到95%。
进一步地,所述云平台分布式处理具体包括以下步骤:
设置n个监测点参与非差PPP处理子任务划分,设置q个基准站参与非差PPP处理子任务划分,设置p台云服务器参与非差PPP处理子任务划分;
通过非差PPP处理子任务划分模块将监测点与基准站的处理任务划分到不同的云处理服务器上进行并行处理;
对于监测点与基准站的短弧法方程回收并合并,划分到云服务器上进行并行处理;
进行模糊度固定处理。
本发明的有益效果如下:
1)提出了近实时长基线解算算法,解决了长基线近实时处理精度不高的问题,核心创新点:利用超快速Multi-GNSS(Mulit-Global navigation satellite system,多模全球卫星定位***)卫星轨道产品减少轨道误差对基线解算的影响、利用多次残差迭代可以提高基线处理精度与稳定性、利用滑动窗口批处理模式可以实现近实时高精度基线解算、利用多基准站联合处理模式可以提高基线解算的精度与稳定性、采用创新的模糊度处理策略可以提高长基线模糊度固定率进一步提高精度。
2)基于非差处理模式,通过虚拟观测值方式施加双差模糊度约束,可以灵活实现多基准站处理技术,有利于更细粒度进行分布式子任务划分,可以充分利用分布式计算优势,子任务划分灵活,提高数据处理效率;提出了短弧法方程合并成长弧法方程的方法,可以避免滑动窗口批处理中的冗余处理,进一步提高解算效率。
附图说明
图1是本发明多基站Multi-GNSS近实时长基线形变监测总体流程图。
图2是本发明基准站选择流程图。
图3是本发明基准站选择示意图。
图4是本发明Multi-GNSS近实时长基线解算流程图。
图5是本发明分布式解算原理示意图。
具体实施方式
为了提高长基线解算的稳定性与精度,本发明采用多基准站联合处理模式。由于单基准站数据中断、卫星观测几何周期性变化、基准站稳定性不同、基准站数据质量不同、基准站与监测点几何距离不同、基准站相对于监测点的几何分布不同,选择不合适的基准站对监测精度影响比较大。本发明提出了一种监测点周边多基准站选择策略,并根据一定阈值自动选择参与解算基准站,初步保证长基线解算的稳定性与精度。
为了提高长基线解算精度,本发明利用超快速Multi-GNSS卫星轨道产品减少轨道误差对基线解算的影响;利用多次后验残差编辑迭代,尽量提出数据异常对基线解算的影响;利用滑动窗口批处理模式,可以实现近实时高精度基线解算;利用多基准站联合处理模式,可以提高基线解算的精度与稳定性;采用创新的模糊度处理策略,进一步提高基线解算精度。
为了提高数据处理效率,本发明利用法方程叠加技术,充分利用以前法方程信息,规避解算重复性;利用共用基准站处理,规避相同基准站的重复处理;利用非差处理模式,可以充分利用分布式计算优势,提高数据处理效率。
下文中,结合附图和实施例对本发明作进一步阐述。
图1是本发明多基站Multi-GNSS近实时长基线形变监测总体流程图,包括以下步骤:根据监测点坐标,结合基准站分析,利用基准站选择策略,给出参与解算的基准站列表;根据参与解算的监测点数据进行子任务划分,执行云平台分布式处理,提高处理效率。
图2是本发明基准站选择流程图,包括以下步骤:第一步根据监测点坐标,筛选距离100km以内的基准站;第二步根据基准站坐标稳定信息,筛选绝对坐标点位重复率小于水平3mm、垂直6mm的基准站;第三步根据基准站数据质量信息,筛选满足准入标准的基准站;第四步根据方位角,将候选的基准站划分到六个区,如图3所示(1区,0°~60°;2区,60°~120°;3区,120°~180°;4区,180°~240°;5区,240°~300°;6区,300°~360°),在每个区选择距离监测点最近的基准站;最后给出参与解算的基准站列表。
本发明优选实施例:70公里基线,监测精度为水平2mm、垂直4mm,解算更新频率为1小时。
图4是本发明Multi-GNSS近实时长基线解算流程图。首先利用超快速Multi-GNSS卫星轨道产品,对基准站和监测点进行非差PPP处理,估计参数包括:卫星钟差、接收机钟差、坐标、天顶对流层延迟、卫星***间延迟(ISB)、卫星频率间延迟(IFB),并进行后验残差编辑迭代,剔除异常数据,生成当前弧段的法方程信息(公式、矩阵详细描述、各个重要参数的具体描述):
其中,为GPS无电离层载波相位观测值,为GLONASS(俄罗斯格洛纳斯卫星导航***)无电离层载波相位观测值,为BDS(中国北斗卫星导航***)无电离层载波相位观测值,G代表GPS,R代表GLONASS,C代表BDS,i代表接收机,j代表卫星编号,ρ是卫星到接收机的几何距离,是卫星位置[xj yj zj]T与接收机坐标[xi yi zi]T的函数,M是对流层模型映射函数,T是天顶对流层延迟,c是光速,λ是波长,是接收机GPS时下钟差,dtj是卫星钟差,是卫星j到接收机i之间的信号传播时间,ISBR-G与ISBC-G是***间偏差,lj是GLONASS卫星的频率标号,IFBi是接收机载波相位频间偏差,bi与bj分别为接收机与卫星的相位小数偏差,B是模糊度参数。
本发明在做非差PPP处理时,采用超快速Multi-GNSS卫星轨道产品对[xj yj zj]T进行固定;模糊度固定采用双差约束模式,将bi与bj统一到非差模糊度参数中,以后提到模糊度包含了硬件相位小数偏差项。非差PPP处理需要估计的参数包括,接收机坐标[xi yizi]T、接收机钟差卫星钟差dtk、***间偏差ISBR-G与ISBC-G、频间偏差IFBi、天顶对流层延迟T、模糊度参数B。
公式(2)是误差方程,A是系数矩阵,X是待估参数向量,V是残差向量,是先验值的虚拟观测值,L为观测值减去模型计算值,P是权值阵,X0与是参数向量的先验值和取值阵,相应的法方程为,NX=W, N是法方程。为了提高解的精度与稳定性,需求对后验残差进行迭代编辑,剔除异常数据,基本原理如下,第k个历元的验后残差相对于第k-1历元残差的变化量可表示为,
Δvk=v(k)-v(k-1) (3)
连续观测弧段内各历元的残差平均值和均方差可表示为,
剔除最大的残差变化量max(Δv)后,对各历元的残差变化再次求均值与均方差,是观测值个数,
如果满足以下条件,
δ>k0δnew
max(Δv)>J0 (6)
其中k0、J0是设定的阈值,则历元k可能发生周跳和粗差,继续检查第k+1个历元。如果第k+1个历元也超限,而且残差变化量也超限,则认为k历元为粗差。否则,认为在第k个历元发生了周跳。把第k个历元作为新残差弧段的开始历元,重新进行判断。
经过多次非差PPP(Precise Point Positioning,精密单点定位)处理与后验残差编辑,可以得到干净的法方程N,为了增强解的稳定性,需要进行长弧段解算,为了提高解算效率避免冗余处理,利用滑动窗口批处理,利用之前m-1个法方程与当前法方程合并,相邻两个方法的边界约束条件如公式(7)所示,其中Vk,k+1是边界约束条件的虚拟观测值,Xk与Xk+1是两个弧段的估计参数改正值,Pk,k+1是虚拟观测值的权重值,Lk,k+1是观测值减去模型计算值。相邻法方程合并实际上是带附加虚拟观测值约束的平差问题,采用最小方差估计,性能函数最小,满足方程(8)。
Vk,k+1=Φ(tk+1,tk)Xk-Xk+1,Pk,k+1-Lk,k+1 (7)
其中,tk是历元k的时间,tk+1是历元k+1的时间,Φ(tk+1,tk)是状态转移矩阵,Vk与Pk是k历元残差向量与权重值,Vk+1与Pk+1是k+1历元残差向量与权重值,min代表让这个公式表述的值最小。
根据可得,
其中,Nk是k历元的法方程,Nk+1是k+1历元的法方程。考虑k历元法方程、k+1历元法方程以及历元间约束信息的,合并后的法方程为公式(10)所示,利用该合并方法可以充分利用长弧解算的稳定性并可以避免冗余非差PPP处理提高处理效率。
为了避免单基准站故障,减弱卫星几何对结果的影响,增强结果的稳定性,采用多基准站模式,采用双差模糊度约束加入的法方程中实现模糊度固定与多基准站联合处理效果。双差模糊度问题等效于对四个非差模糊度加上适当的约束条件,定义双差模糊度的虚拟观测值VB,
其中,为固定的无电离层组合的双差模糊度,PB为虚拟观测值的权值,D是非差模糊度的投影矩阵,
其中,和代表两颗卫星,代表两个接收机,为测站与卫星之间的非差模糊度。考虑以上约束,最小二乘解的形式可改写为公式(13),实现从非差到双差基线处理的转换,如果施加多条基线的双差模糊度约束,即可实现多基准站的融合处理模式,增强处理结果的可靠性。
其中,Xnew、Nnew、Wnew分别为加了双插模糊度约束的X、N、W。
与其它双差模糊度固定策略相比,利用本发明的算法,固定模糊度可以从法方程中消去,计算效率可以得到大幅度提升。
采用创新的模糊度固定策略,多基准站采用独立基线对,基线对上卫星对也采用独立卫星对,采用类似选择独立基线的策略进行独立卫星对的选取,对基线对上的卫星对的宽巷(WL)/窄巷(NL)模糊度固定概率进行筛选。
采用宽巷与窄巷组合模式进行模糊度固定,首先固定宽巷模糊度Bw,B1是L1频率上的模糊度,B2是L2频率上的模糊度,宽巷双差模糊度由四个非差模糊度组成,如公式(15)所示,宽巷双差模糊度方差值由误差传播公式(16)得出。根据公式(17)可以得出双差宽巷模糊度的固定概率Probw,其中B是模糊度参数,σ是估计方差值,是B的最近整数值。
Bw=B1-B2 (14)
利用公式(19)与(20),可以得到窄巷双差模糊的估计值以及估计方差其中与是无电离层组合双差模糊度值与估计方差值,可以通过公式(17)得到窄巷双差模糊度固定概率Probn,
Probw>Prob0,Probn>Prob0 (21)
对通过筛选公式(21)的卫星对进行模糊度固定,其中Prob0是模糊度固定概率阈值,将模糊度固定信息加入到法方程中更新法方程,重复进行模糊度筛选与固定处理直到95%的模糊度被固定。
图5是本发明分布式解算原理示意图。如图5所示,设置参与某次处理的监测点Mon1,Mon2,...,Monn共计n个,参与处理的基准站Ref1,Ref2,...,Refq共计q个,参与处理的云服务器Ser1,Ser2,...,SerP共p台,由于本发明采用非差PPP处理模式,可以通过非差PPP处理子任务划分模块将监测点与基准站的处理任务划分到不同的云处理服务器上进行并行处理,对于监测点与基准站的短弧法方程合并有可以进行划分到云处理服务器上进行并行处理,在模糊度固定阶段可以根据公式(11)~(13)施加双差模糊度约束实现单个监测点的多基准站联合坐标解算,通过以上策略,可以充分利用分布式计算优势,更适用于海量数据处理的场景。
本发明核心处理代码利用C/C++,自动化控制利用Python语言。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
Claims (9)
1.一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,包括以下步骤:
根据监测点坐标,结合基准站分析,利用基准站选择策略,给出参与解算的基准站列表;
根据参与解算的监测点数据进行子任务划分,执行云平台分布式处理;
云平台分布式处理利用超快速Multi-GNSS卫星轨道产品减少轨道误差对Multi-GNSS近实时长基线解算的影响;
所述Multi-GNSS近实时长基线解算具体包括以下步骤:
利用超快速Multi-GNSS卫星轨道产品,对基准站和监测点进行非差PPP处理,并进行后验残差编辑迭代,剔除异常数据,生成当前弧段的法方程信息,当前弧段的法方程之前有m-1个法方程;
利用滑动窗口批处理,将当前弧段的法方程与之前m-1个法方程合并成长弧法方程;
采用模糊度固定策略,多基准站采用独立基线对,基线对上卫星对也采用独立卫星对,对基线对上的卫星对的宽巷/窄巷模糊度固定概率进行筛选;
采用宽巷与窄巷组合模式进行模糊度固定,将模糊度固定信息加入到长弧法方程中更新长弧法方程,重复进行宽巷/窄巷模糊度筛选与模糊度固定处理。
2.如权利要求1所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,所述基准站选择策略具体包括以下步骤:
步骤1,根据监测点坐标,筛选距离100km以内的基准站;
步骤2,根据基准站坐标稳定信息,筛选绝对坐标点位重复率小于水平3mm、垂直6mm的基准站;
步骤3,根据基准站数据质量信息,筛选满足准入标准的基准站;
步骤4,根据方位角,将候选的基准站划分到六个区,在每个区选择距离监测点最近的基准站;
步骤5,给出参与解算的基准站列表。
3.如权利要求2所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,候选的基准站按照以下方法划分到六个区:
1区,0°~60°;
2区,60°~120°;
3区,120°~180°;
4区,180°~240°;
5区,240°~300°;
6区,300°~360°。
4.如权利要求1所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,对基准站和监测点进行非差PPP处理需要估计的参数包括包括:卫星钟差、接收机钟差、坐标、天顶对流层延迟、卫星***间延迟、卫星频率间延迟。
5.如权利要求1所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,所述模糊度固定策略采用双差模糊度约束加入长弧法方程,实现模糊度固定与多基准站联合处理。
6.如权利要求5所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,所述双差模糊度约束等效于对四个非差模糊度加上约束条件,定义双差模糊度的虚拟观测值VB:
其中,为固定的无电离层组合的双差模糊度,PB为双差模糊度虚拟观测值的权值,D是非差模糊度的投影矩阵,X是待估参数向量;
其中,和代表两颗卫星,代表两个接收机,为测站与卫星之间的非差模糊度,为测站与卫星之间的非差模糊度,为测站与卫星之间的非差模糊度,为测站与卫星之间的非差模糊度;
基于公式(11)和公式(12)的约束,最小二乘解的形式改写为公式(13),实现从非差到双差基线处理的转换:
其中,A是系数矩阵,P是权值阵,L是观测值减去模型计算值,NX=W,N是法方程,Xnew为加了双插模糊度约束的X,Nnew为加了双插模糊度约束的法N,Wnew为加了双插模糊度约束的W。
7.如权利要求6所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,如果施加多条基线的双差模糊度约束,即可实现多基准站联合融合处理。
8.如权利要求7所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,在模糊度固定处理中模糊度需要达到95%。
9.如权利要求7所述的一种基于多基站的Multi-GNSS长基线近实时形变监测方法,其特征在于,所述云平台分布式处理具体包括以下步骤:
设置n个监测点参与非差PPP处理子任务划分,设置q个基准站参与非差PPP处理子任务划分,设置p台云服务器参与非差PPP处理子任务划分;
通过非差PPP处理子任务划分模块将监测点与基准站的处理任务划分到不同的云处理服务器上进行并行处理;
对于监测点与基准站的短弧法方程回收并合并,划分到云服务器上进行并行处理;
进行模糊度固定处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710672999.8A CN107421434B (zh) | 2017-08-08 | 2017-08-08 | 多基站Multi-GNSS长基线近实时形变监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710672999.8A CN107421434B (zh) | 2017-08-08 | 2017-08-08 | 多基站Multi-GNSS长基线近实时形变监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107421434A CN107421434A (zh) | 2017-12-01 |
CN107421434B true CN107421434B (zh) | 2019-09-24 |
Family
ID=60437565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710672999.8A Active CN107421434B (zh) | 2017-08-08 | 2017-08-08 | 多基站Multi-GNSS长基线近实时形变监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107421434B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108196282A (zh) * | 2017-12-06 | 2018-06-22 | 千寻位置网络有限公司 | Gnss定位的执行方法及***、定位装置、服务终端、存储器 |
CN109540095A (zh) * | 2018-12-29 | 2019-03-29 | 北方信息控制研究院集团有限公司 | 基于卫星导航及最小二乘的路基沉降监测方法 |
CN112444187B (zh) * | 2019-08-28 | 2022-09-06 | 千寻位置网络有限公司 | 形变监测方法及装置 |
CN111220065B (zh) * | 2019-12-11 | 2020-10-09 | 湖南联智科技股份有限公司 | 一种多模式联合的北斗变形监测数据处理***及方法 |
CN111045042B (zh) * | 2019-12-20 | 2022-03-04 | 西安空间无线电技术研究所 | 一种基于“云-端”架构的ppp-rtk增强方法与*** |
CN111258999B (zh) * | 2020-01-16 | 2021-07-16 | 山东大学 | 一种基于slr的bds卫星轨道近实时检核服务*** |
CN111610308A (zh) * | 2020-04-27 | 2020-09-01 | 杭州电子科技大学 | 一种基于rtk技术的双层滑坡监测***及方法 |
CN113670252B (zh) * | 2020-05-13 | 2024-01-26 | 千寻位置网络有限公司 | 基于多源信息融合的自适应弧长形变监测方法及其装置 |
CN112013756B (zh) * | 2020-08-27 | 2022-05-17 | 桂林电子科技大学 | 一种双层基线边坡形变监测方法 |
CN112381309B (zh) * | 2020-11-23 | 2022-04-12 | 珠江水利委员会珠江水利科学研究院 | 水库大坝安全监测预警方法、装置、***及存储介质 |
CN112444832A (zh) * | 2020-12-03 | 2021-03-05 | 深圳大学 | 一种一机多天线接收机的周跳修复方法 |
CN117232380B (zh) * | 2023-11-13 | 2024-01-30 | 武汉大学 | 一种带状工程变形监测方法、装置、设备及可读存储介质 |
CN117630988A (zh) * | 2023-11-14 | 2024-03-01 | 中交一公局厦门工程有限公司 | 动态调整优化的北斗测量数据处理云平台 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014065664A1 (en) * | 2012-10-25 | 2014-05-01 | Fugro N.V. | Ppp-rtk method and system for gnss signal based position determination |
CN105158780A (zh) * | 2015-07-24 | 2015-12-16 | 北京跟踪与通信技术研究所 | 一种基于多种导航卫星可互换的导航定位方法 |
CN106371119A (zh) * | 2016-08-23 | 2017-02-01 | 武汉羲和科技有限公司 | 基于云定位的gnss全球高精度增强定位方法 |
CN106871776A (zh) * | 2017-02-14 | 2017-06-20 | 千寻位置网络有限公司 | 一种基于gnss的实时变形监测*** |
-
2017
- 2017-08-08 CN CN201710672999.8A patent/CN107421434B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014065664A1 (en) * | 2012-10-25 | 2014-05-01 | Fugro N.V. | Ppp-rtk method and system for gnss signal based position determination |
CN105158780A (zh) * | 2015-07-24 | 2015-12-16 | 北京跟踪与通信技术研究所 | 一种基于多种导航卫星可互换的导航定位方法 |
CN106371119A (zh) * | 2016-08-23 | 2017-02-01 | 武汉羲和科技有限公司 | 基于云定位的gnss全球高精度增强定位方法 |
CN106871776A (zh) * | 2017-02-14 | 2017-06-20 | 千寻位置网络有限公司 | 一种基于gnss的实时变形监测*** |
Non-Patent Citations (1)
Title |
---|
大规模GNSS数据的分布式处理与实现;陈正生 等;《武汉大学学报(信息科学版)》;20150331;第40卷(第3期);384-388 * |
Also Published As
Publication number | Publication date |
---|---|
CN107421434A (zh) | 2017-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107421434B (zh) | 多基站Multi-GNSS长基线近实时形变监测方法 | |
CN110231037B (zh) | 一种附加钟差模型约束的gnss机动卫星定轨方法 | |
KR102547522B1 (ko) | 쾌속 정밀 측위 방법 및 시스템 | |
CN108415049B (zh) | 提高网络rtk双差宽巷模糊度固定正确率的方法 | |
WO2019228439A1 (zh) | 基于gnss-rtk的定位方法 | |
CN106468774B (zh) | 一种应用于星基增强***的星历星钟改正参数及空间信号完好性参数方法 | |
CN110007320B (zh) | 网络rtk解算方法 | |
Fey et al. | The second extension of the international celestial reference frame: ICRF-Ext. 1 | |
CN105699999B (zh) | 一种固定北斗地基增强***基准站窄巷模糊度的方法 | |
CN108549095B (zh) | 一种区域cors网非差并行增强方法及*** | |
CN107064981B (zh) | 基于gnss的差分定位方法及***、服务终端 | |
CN112817023B (zh) | 基于星基增强服务的无依托式低轨导航增强***与方法 | |
CN107942346B (zh) | 一种高精度gnss电离层tec观测值提取方法 | |
Li et al. | Review of PPP–RTK: Achievements, challenges, and opportunities | |
CN110109158A (zh) | 基于gps、glonass和bds多***的事后超快速rtk定位算法 | |
CN103698790B (zh) | 北斗与gps双***宽巷载波相位混频星间差分组合方法 | |
CN111381264A (zh) | 网络rtk中长基线模糊度固定方法和平台 | |
CN115933356B (zh) | 一种虚拟原子钟的高精度时间同步***和方法 | |
Xia et al. | Observation of BDS-2 IGSO/MEOs yaw-attitude behavior during eclipse seasons | |
CN111308515A (zh) | 基于精密轨道和地面台站数据的北斗卫星机动及异常探测方法 | |
Dai et al. | Real-time precise orbit and clock estimation of multi-GNSS satellites with undifferenced ambiguity resolution | |
CN117665871A (zh) | 一种载波距支持下的gnss非差固定解精密钟差确定方法 | |
Li et al. | Robust Kalman filter for RTK positioning under signal-degraded scenarios | |
Yu et al. | Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products | |
CN116299586B (zh) | 基于广播星历的精密单点定位方法、接收机、设备和介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211026 Address after: 200438 9 / F, 10 / F, 11 / F, 12 / F, 38 Lane 1688, Guoquan North Road, Yangpu District, Shanghai Patentee after: QIANXUN SPATIAL INTELLIGENCE Inc. Patentee after: QIANXUN SI NETWORK (ZHEJIANG) Co.,Ltd. Address before: Room j165, 1st floor, building 64, 1436 Jungong Road, Yangpu District, Shanghai, 200433 Patentee before: QIANXUN SPATIAL INTELLIGENCE Inc. |