CN113805212A - 一种自适应的gnss载波相位差分滑坡监测方法 - Google Patents

一种自适应的gnss载波相位差分滑坡监测方法 Download PDF

Info

Publication number
CN113805212A
CN113805212A CN202111128231.7A CN202111128231A CN113805212A CN 113805212 A CN113805212 A CN 113805212A CN 202111128231 A CN202111128231 A CN 202111128231A CN 113805212 A CN113805212 A CN 113805212A
Authority
CN
China
Prior art keywords
difference
carrier phase
epoch
satellite
value
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
Application number
CN202111128231.7A
Other languages
English (en)
Other versions
CN113805212B (zh
Inventor
孙希延
梁维彬
鞠涛
杜洋
付文涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanning Guidian Electronic Technology Research Institute Co ltd
Guilin University of Electronic Technology
Original Assignee
Nanning Guidian Electronic Technology Research Institute Co ltd
Guilin University of Electronic Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanning Guidian Electronic Technology Research Institute Co ltd, Guilin University of Electronic Technology filed Critical Nanning Guidian Electronic Technology Research Institute Co ltd
Priority to CN202111128231.7A priority Critical patent/CN113805212B/zh
Publication of CN113805212A publication Critical patent/CN113805212A/zh
Application granted granted Critical
Publication of CN113805212B publication Critical patent/CN113805212B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种自适应的GNSS载波相位差分滑坡监测方法，通过接收监测站和参考站GNSS观测数据，得到观测到的卫星的伪距、载波相位、载噪比，将观测数据中的载波相位观测量作双差组成载波相位双差观测方程组，通过抗差卡尔曼滤波进行解算，得到监测站位置的浮点解，并通过剔除观测质量较差卫星来对载波相位的整周模糊度参数进行筛选，采用MLAMBDA算法搜索整周模糊度，最终得到监测站位置的固定解，另外通过统计所有卫星的连续锁定历元数大于设定值的比例，并与设定比例值比较，自动调整滤波策略，将连续解算与单历元解算的优势结合，保证恶劣观测环境下的固定率及定位精度。

Description

一种自适应的GNSS载波相位差分滑坡监测方法
技术领域
本发明涉及卫星导航定位技术领域，尤其涉及一种自适应的GNSS载波相位差分滑坡监测方法。
背景技术
随着全球卫星导航系统(GNSS，GlobalNavigation Satellite System)的发展，高精度差分定位技术已广泛应用于山体滑坡等高精度监测领域。传统高精度差分定位技术通常构建载波相位双差观测方程通过卡尔曼滤波技术求的高精度定位结果，但当观测环境较差时使用传统卡尔曼滤波技术与全部模糊度固定的方法，通常会出现周跳频繁、固定困难、精度达不到要求等情况。。
发明内容
本发明的目的在于提供一种自适应的GNSS载波相位差分滑坡监测方法，旨在解决传统高精度差分定位技术在恶劣观测环境下的固定率及定位精度较差的技术问题。
为实现上述目的，本发明采用一种自适应的GNSS载波相位差分滑坡监测方法，包括下列步骤：
步骤1、在监测区域放置参考站接收机和监测站接收机，接收GNSS观测数据，获得观测到的卫星数据；
步骤2、根据所述GNSS观测数据计算监测站接收机的位置，并计算每颗卫星仰角；
步骤3、选取参考站与监测站间的共视卫星，构建载波相位双差观测方程；
步骤4、求解所述载波相位双差观测方程，计算得到监测站位置及站间单差整周模糊度的浮点状态矩阵和浮点状态协方差矩阵；
步骤5、筛选站间单差整周模糊度并进行单差转双差处理，搜索整周模糊度的最优估计值；
步骤6、对所述最优估计值判定确认，输出整周模糊度固定解或浮点解；
步骤7、将双差模糊度恢复成单差模糊度，进入下一历元解算。
其中，观测到的卫星数据包括卫星的伪距、载波相位和载噪比。
其中，在根据所述GNSS观测数据计算监测站接收机的位置，并计算每颗卫星仰角的过程中，根据接收到的GNSS观测数据，使用单点定位计算监测站接收机的位置，为抗差卡尔曼滤波的时间更新过程提供位置初始值，并计算每颗卫星仰角。
其中，在步骤4中，具体为：使用抗差卡尔曼滤波算法进行时间和状态的更新，包括下列步骤：
统计当前历元连续锁定历元数大于设定值的卫星占卫星总数的比例，
当比例大于或等于设定值时，采用连续解算模式；
当比例小于设定值时，采用单历元解算模式；
根据时间更新过程的状态预测值及观测值，得到新的浮点状态矩阵和浮点状态协方差矩阵。
统计当前历元连续锁定历元数大于设定值的卫星占卫星总数比例的计算为抗差卡尔曼的时间更新，根据根据时间更新过程的状态预测值及观测值可以进行抗差卡尔曼的状态更新。
其中，所述连续解算模式的更新过程包括下列步骤：
若当前历元为首历元，则根据监测站单点定位位置及载波相位单差初始值来初始化浮点状态矩阵和浮点状态协方差矩阵；
若当前历元为非首历元且上一历元的解为固定解，则根据上一历元的固定解状态矩阵和状态协方差矩阵来更新当前历元的浮点状态矩阵和浮点状态协方差矩阵；
若当前历元为非首历元且上一历元的解为浮点解，则根据上一历元的浮点解状态矩阵和状态协方差矩阵来更新当前历元的浮点状态矩阵及浮点状态协方差矩阵。
其中，采用所述单历元解算模式的更新时，根据当前历元监测站单点定位位置及载波相位单差初始值来更新浮点状态矩阵和浮点状态协方差矩阵，并且增加滤波迭代次数。
其中，在筛选站间单差整周模糊度并进行单差转双差处理，搜索整周模糊度的最优估计值的过程中，包括下列步骤：
将各卫星剔除标志置0；
根据预设的卫星仰角门限值、载噪比门限值以及卫星连续锁定历元数门限值，将低于门限值的卫星剔除标志置1；
将卫星剔除标志为0的卫星在状态矩阵中的单差载波相位和状态协方差矩阵中的对应部分提取出来；
对提取出来的状态矩阵和状态协方差矩阵进行单差转双差处理；
采用MLAMBDA算法搜索整周模糊度的最优估计值。
其中，在步骤6中，采用Ratio算法对步骤5获得的整周模糊度最优估计值进行确认，判定成功则输出整周模糊度固定解，判定失败则输出浮点解。
步骤1至步骤6为单次历元解算过程，步骤7作为连续解算的准备步骤，根据计算要求自适应转回步骤1进入循环解算操作。
本发明的一种自适应的GNSS载波相位差分滑坡监测方法，通过接收监测站和参考站GNSS观测数据，得到观测到的卫星的伪距、载波相位、载噪比，将观测数据中的载波相位观测量作双差组成载波相位双差观测方程组，通过抗差卡尔曼滤波进行解算，得到监测站位置的浮点解，并通过剔除观测质量较差卫星来对载波相位的整周模糊度参数进行筛选，采用MLAMBDA算法搜索整周模糊度，最终得到监测站位置的固定解，另外通过统计所有卫星的连续锁定历元数大于设定值的比例，并与设定比例值比较，自动调整滤波策略，将连续解算与单历元解算的优势结合，保证恶劣观测环境下的固定率及定位精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明的一种自适应的GNSS载波相位差分滑坡监测方法的流程示意图。
图2是本发明的一种自适应的GNSS载波相位差分滑坡监测方法的具体执行流程图。
图3是本发明的抗差卡尔曼滤波过程流程图。
图4是本发明的模糊度筛选流程图。
图5是本发明的载波相位双差模型示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
请参阅图1，本发明提出了一种自适应的GNSS载波相位差分滑坡监测方法，包括下列步骤：
S1、在监测区域放置参考站接收机和监测站接收机，接收GNSS观测数据，获得观测到的卫星数据；
S2、根据所述GNSS观测数据计算监测站接收机的位置，并计算每颗卫星仰角；
S3、选取参考站与监测站间的共视卫星，构建载波相位双差观测方程；
S4、求解所述载波相位双差观测方程，计算得到监测站位置及站间单差整周模糊度的浮点状态矩阵和浮点状态协方差矩阵；
S5、筛选站间单差整周模糊度并进行单差转双差处理，搜索整周模糊度的最优估计值；
S6、对所述最优估计值判定确认，输出整周模糊度固定解或浮点解；
S7、将双差模糊度恢复成单差模糊度，进入下一历元解算。
详细执行流程请参阅图2，所述自适应的GNSS载波相位差分滑坡监测方法基于复杂环境下的高精度要求，采用抗差卡尔曼滤波与部分模糊度固定策略，并根据所有卫星的周跳情况将滤波策略调整，使得单历元解算对周跳不敏感和连续解算的高精度优势充分结合，保证恶劣观测环境下的固定率及精度。
以下本发明基于实施步骤和具体参数进一步说明：
步骤S1：通过放置在监测区域的参考站接收机和监测站接收机，同时接收GNSS观测数据，得到观测到的卫星的伪距
Figure BDA0003279491930000041
与
Figure BDA0003279491930000042
载波相位
Figure BDA0003279491930000043
与
Figure BDA0003279491930000044
载噪比
Figure BDA0003279491930000045
与
Figure BDA0003279491930000046
步骤S2：根据接收到的GNSS数据，使用单点定位计算监测站接收机的位置(xs,r,ys,r,zs,r)，为抗差卡尔曼滤波的时间更新过程提供位置初始值，并计算每颗卫星仰角ELi；
具体地，计算监测站接收机位置的单点定位数学表达式为：
Figure BDA0003279491930000051
Figure BDA0003279491930000052
其中，
Figure BDA0003279491930000053
表示监测站接收机观测到的卫星i的伪距值；
Figure BDA0003279491930000054
表示监测站与卫星i的几何距离；δtmoni表示监测站接收机钟差，为待求量；δti表示卫星i的钟差，可通过星历参数求得；I、T分别表示电离层误差和对流层误差，可通过建模得到；
Figure BDA0003279491930000055
表示伪距的随机性噪声误差；(xs,r,ys,r,zs,r)表示待求的监测站接收机位置；
Figure BDA0003279491930000056
可通过卫星i的星历求得，表示卫星i的位置。
根据数学表达式可知，由于有4个未知参数，包括3个接收机位置参数和1个接收机钟差参数，则至少观测到4颗卫星组成方程组后，即可得到单点定位结果，解算方法可采用最小二乘法。
具体地，计算各颗卫星仰角ELi的数学表达式为：
Figure BDA0003279491930000057
Figure BDA0003279491930000058
Figure BDA0003279491930000059
Figure BDA00032794919300000510
其中，(Δx,Δy,Δz)为观测向量，S为变换矩阵，ELi为卫星i的仰角，Lat和Lon分别为监测站接收机坐标对应的纬度和经度，可由接收机坐标(xs,r,ys,r,zs,r)经坐标变换得到。
步骤S3：选取参考站与监测站间的共视卫星i，构建载波相位双差观测方程；
具体地，对于共视卫星i，监测站接收机和参考站接收机的载波相位观测值分别可表示为：
Figure BDA0003279491930000061
Figure BDA0003279491930000062
其中，
Figure BDA0003279491930000063
和
Figure BDA0003279491930000064
分别表示监测站接收机和参考站接收机的载波相位观测值；λ表示波长；f表示频率；
Figure BDA0003279491930000065
和
Figure BDA0003279491930000066
分别表示监测站接收机和参考站接收机与共视卫星i的几何距离，计算方法同式(2)；δtmoni和δtref分别表示监测站接收机钟差和参考站接收机钟差；δti表示共视卫星i的卫星钟差；Imoni、Iref和Tmoni、Tref表示监测站和参考站的电离层及对流层误差；
Figure BDA0003279491930000067
和
Figure BDA0003279491930000068
分别表示监测站接收机和参考站接收机对卫星i观测时载波相位的整周模糊度；
Figure BDA0003279491930000069
和
Figure BDA00032794919300000610
表示载波相位观测的随机性噪声误差。
假设共视卫星i为当前历元仰角最高的卫星，将其作为参考为卫星。对于共视卫星i，在载波相位观测值
Figure BDA00032794919300000611
和
Figure BDA00032794919300000612
之间作站间单差，得到站间单差值
Figure BDA00032794919300000613
再对于共视卫星j，在载波相位观测值
Figure BDA00032794919300000614
和
Figure BDA00032794919300000615
同样作站间单差得到
Figure BDA00032794919300000616
可消除卫星钟差，短基线情况下可消除电离层误差及对流层误差。
Figure BDA00032794919300000617
Figure BDA00032794919300000618
在两颗卫星i和j的站间单差值之间作差，得到双差载波相位观测方程：
Figure BDA00032794919300000619
经过双差后，可以进一步消除原观测方程中的接收机钟差以及部分随机误差。
步骤S4，使用抗差卡尔曼滤波算法求解载波相位双差观测方程，计算得到监测站位置及站间单差整周模糊度的浮点状态矩阵和浮点状态协方差矩阵；
具体地，如图3所示，首先，统计当前历元连续锁定历元数大于设定值的卫星占卫星总数的比例Rlock，将比例Rlock与设定值Rlockmin进行比较，根据比较结果，抗差卡尔曼滤波的时间更新过程分为两种情况。
情况1、比例Rlock小于设定值Rlockmin，采用单历元解算模式，增加滤波迭代次数，根据当前历元监测站单点定位位置及载波相位单差初始值来更新浮点状态矩阵及浮点状态协方差矩阵。
情况2、比例Rlock大于或等于设定值Rlockmin，采用连续解算模式，抗差卡尔曼滤波时间更新过程可分为3种情况。
情况1、若当前历元为首历元，则根据监测站单点定位位置及载波相位单差初始值来初始化浮点状态矩阵及浮点状态协方差矩阵。
情况2、若当前历元为非首历元且上一历元的解为固定解，则根据上一历元的固定解状态矩阵及状态协方差矩阵来更新当前历元的浮点状态矩阵及浮点状态协方差矩阵。
情况3、若当前历元为非首历元且上一历元的解为浮点解，则根据上一历元的浮点解状态矩阵及状态协方差矩阵来更新当前历元的浮点状态矩阵及浮点状态协方差矩阵。
抗差卡尔曼滤波的状态更新过程：根据时间更新过程的状态预测值及观测值，得到新的浮点估计值及其协方差矩阵，其过程为：
Figure BDA0003279491930000071
Figure BDA0003279491930000072
Figure BDA0003279491930000073
Figure BDA0003279491930000074
Figure BDA0003279491930000075
其中，Xk为系统的状态向量，Kk为滤波增益矩阵，Zk为系统的观测向量，Hk为观测矩阵，wk,i为权因子，vi为观测向量Z对应的残差向量，
Figure BDA0003279491930000076
为vi的均方差，σ0为方差因子，k0和ki可根据需要分别取1.5～2.5、3.0～5.0，
Figure BDA0003279491930000077
为等价权阵，
Figure BDA0003279491930000078
为等价协方差阵可由等价权阵
Figure BDA0003279491930000079
求逆得到。
步骤S5：对站间单差整周模糊度进行筛选，选择出符合预设条件的站间单差整周模糊度进行单差转双差处理，采用MLAMBDA算法搜索整周模糊度的最优估计值；
具体地，如图4所示，为了避免质量较差卫星对整周模糊度的搜索效率的影响，对步骤S4中得到的单差整周模糊度进行筛选。首先将每颗卫星的剔除标志置0，然后将每颗卫星的仰角与预设值进行比较，将小于预设值的卫星剔除标志设置为1，将每颗卫星的载噪比与预设值进行比较，将小于预设值的卫星剔除标志设置为1，将每颗卫星的连续锁定历元数与预设值进行比较，将小于预设值的卫星剔除标志设置为1，然后将卫星剔除标志为0的卫星在状态矩阵中的单差载波相位及协方差阵中的对应部分提取出来，对筛选后的状态矩阵及其协方差矩阵进行单差转双差处理，最后采用MLAMBDA算法搜索整周模糊度的最优估计值。
步骤S6：采用Ratio算法对步骤S5得到的整周模糊度最优估计值进行确认，判定成功则输出整周模糊度固定解，判定失败则输出浮点解。对模糊度进行确认的公式为：
Figure BDA0003279491930000081
其中，为N为经过抗差卡尔曼滤波之后的模糊度浮点值，QN为对应的协方差矩阵，R0为预设的门限值，
Figure BDA0003279491930000082
为步骤5中模糊度的次优值，
Figure BDA0003279491930000083
为步骤5中模糊度的最优值。
步骤S7：得到浮点解或者固定解之后，将双差模糊度恢复成单差模糊度，进入下一历元解算，双差模糊度恢复成单差模糊度的数学表达式为：
Figure BDA0003279491930000084
以上所揭露的仅为本发明一种较佳实施例而已，当然不能以此来限定本发明之权利范围，本领域普通技术人员可以理解实现上述实施例的全部或部分流程，并依本发明权利要求所作的等同变化，仍属于发明所涵盖的范围。

Claims (8)

1.一种自适应的GNSS载波相位差分滑坡监测方法，其特征在于，包括下列步骤：
步骤1、在监测区域放置参考站接收机和监测站接收机，接收GNSS观测数据，获得观测到的卫星数据；
步骤2、根据所述GNSS观测数据计算监测站接收机的位置，并计算每颗卫星仰角；
步骤3、选取参考站与监测站间的共视卫星，构建载波相位双差观测方程；
步骤4、求解所述载波相位双差观测方程，计算得到监测站位置及站间单差整周模糊度的浮点状态矩阵和浮点状态协方差矩阵；
步骤5、筛选站间单差整周模糊度并进行单差转双差处理，搜索整周模糊度的最优估计值；
步骤6、对所述最优估计值判定确认，输出整周模糊度固定解或浮点解；
步骤7、将双差模糊度恢复成单差模糊度，进入下一历元解算。
2.如权利要求1所述的自适应的GNSS载波相位差分滑坡监测方法，其特征在于，观测到的卫星数据包括卫星的伪距、载波相位和载噪比。
3.如权利要求1所述的自适应的GNSS载波相位差分滑坡监测方法，其特征在于，在根据所述GNSS观测数据计算监测站接收机的位置，并计算每颗卫星仰角的过程中，根据接收到的GNSS观测数据，使用单点定位计算监测站接收机的位置，为抗差卡尔曼滤波的时间更新过程提供位置初始值，并计算每颗卫星仰角。
4.如权利要求1所述的自适应的GNSS载波相位差分滑坡监测方法，其特征在于，在步骤4中，具体为：使用抗差卡尔曼滤波算法进行时间和状态的更新，包括下列步骤：
统计当前历元连续锁定历元数大于设定值的卫星占卫星总数的比例，
当比例大于或等于设定值时，采用连续解算模式；
当比例小于设定值时，采用单历元解算模式；
根据时间更新过程的状态预测值及观测值，得到新的浮点状态矩阵和浮点状态协方差矩阵。
5.如权利要求4所述的自适应的GNSS载波相位差分滑坡监测方法，其特征在于，所述连续解算模式的更新过程包括下列步骤：
若当前历元为首历元，则根据监测站单点定位位置及载波相位单差初始值来初始化浮点状态矩阵和浮点状态协方差矩阵；
若当前历元为非首历元且上一历元的解为固定解，则根据上一历元的固定解状态矩阵和状态协方差矩阵来更新当前历元的浮点状态矩阵和浮点状态协方差矩阵；
若当前历元为非首历元且上一历元的解为浮点解，则根据上一历元的浮点解状态矩阵和状态协方差矩阵来更新当前历元的浮点状态矩阵及浮点状态协方差矩阵。
6.如权利要求4所述的自适应的GNSS载波相位差分滑坡监测方法，其特征在于，采用所述单历元解算模式的更新时，根据当前历元监测站单点定位位置及载波相位单差初始值来更新浮点状态矩阵和浮点状态协方差矩阵，并且增加滤波迭代次数。
7.如权利要求1所述的自适应的GNSS载波相位差分滑坡监测方法，其特征在于，在筛选站间单差整周模糊度并进行单差转双差处理，搜索整周模糊度的最优估计值的过程中，包括下列步骤：
将各卫星剔除标志置0；
根据预设的卫星仰角门限值、载噪比门限值以及卫星连续锁定历元数门限值，将低于门限值的卫星剔除标志置1；
将卫星剔除标志为0的卫星在状态矩阵中的单差载波相位和状态协方差矩阵中的对应部分提取出来；
对提取出来的状态矩阵和状态协方差矩阵进行单差转双差处理；
采用MLAMBDA算法搜索整周模糊度的最优估计值。
8.如权利要求1所述的自适应的GNSS载波相位差分滑坡监测方法，其特征在于，在步骤6中，采用Ratio算法对步骤5获得的整周模糊度最优估计值进行确认，判定成功则输出整周模糊度固定解，判定失败则输出浮点解。
CN202111128231.7A 2021-09-26 2021-09-26 一种自适应的gnss载波相位差分滑坡监测方法 Active CN113805212B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111128231.7A CN113805212B (zh) 2021-09-26 2021-09-26 一种自适应的gnss载波相位差分滑坡监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111128231.7A CN113805212B (zh) 2021-09-26 2021-09-26 一种自适应的gnss载波相位差分滑坡监测方法

Publications (2)

Publication Number Publication Date
CN113805212A true CN113805212A (zh) 2021-12-17
CN113805212B CN113805212B (zh) 2023-11-07

Family

ID=78896678

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111128231.7A Active CN113805212B (zh) 2021-09-26 2021-09-26 一种自适应的gnss载波相位差分滑坡监测方法

Country Status (1)

Country Link
CN (1) CN113805212B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114859389A (zh) * 2022-04-18 2022-08-05 华力智芯(成都)集成电路有限公司 一种gnss多系统抗差自适应融合rtk解算方法
CN115480278A (zh) * 2022-08-01 2022-12-16 北方雷科(安徽)科技有限公司 一种适应于复杂多径环境下的双天线定向抗差算法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108317949A (zh) * 2018-02-07 2018-07-24 桂林电子科技大学 一种rtk高精度差分定位形变监测系统及方法
CN108871179A (zh) * 2018-05-07 2018-11-23 桂林电子科技大学 基于载波相位差分动静态融合的形变监测定位方法及装置
CN110132121A (zh) * 2019-05-10 2019-08-16 国网浙江省电力有限公司信息通信分公司 北斗三号双频非组合rtk定位的输电杆塔形变监测方法
WO2019228439A1 (zh) * 2018-06-01 2019-12-05 浙江亚特电器有限公司 基于gnss-rtk的定位方法
CN111025357A (zh) * 2019-12-31 2020-04-17 中国科学院国家授时中心 一种高精度gnss虚拟snr加权定位方法
CN111610308A (zh) * 2020-04-27 2020-09-01 杭州电子科技大学 一种基于rtk技术的双层滑坡监测系统及方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108317949A (zh) * 2018-02-07 2018-07-24 桂林电子科技大学 一种rtk高精度差分定位形变监测系统及方法
CN108871179A (zh) * 2018-05-07 2018-11-23 桂林电子科技大学 基于载波相位差分动静态融合的形变监测定位方法及装置
WO2019228439A1 (zh) * 2018-06-01 2019-12-05 浙江亚特电器有限公司 基于gnss-rtk的定位方法
CN110132121A (zh) * 2019-05-10 2019-08-16 国网浙江省电力有限公司信息通信分公司 北斗三号双频非组合rtk定位的输电杆塔形变监测方法
CN111025357A (zh) * 2019-12-31 2020-04-17 中国科学院国家授时中心 一种高精度gnss虚拟snr加权定位方法
CN111610308A (zh) * 2020-04-27 2020-09-01 杭州电子科技大学 一种基于rtk技术的双层滑坡监测系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王守华 等: "基于GPS 载波相位的边坡形变监测方法", 科学技术与工程, vol. 17, no. 25 *
雷波;李青;李雄;: "RTK技术在滑坡监测中的应用及系统设计", 自动化仪表, no. 06 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114859389A (zh) * 2022-04-18 2022-08-05 华力智芯(成都)集成电路有限公司 一种gnss多系统抗差自适应融合rtk解算方法
CN115480278A (zh) * 2022-08-01 2022-12-16 北方雷科(安徽)科技有限公司 一种适应于复杂多径环境下的双天线定向抗差算法

Also Published As

Publication number Publication date
CN113805212B (zh) 2023-11-07

Similar Documents

Publication Publication Date Title
CN107422354B (zh) 一种模糊度固定的ppp/sins紧组合定位定姿方法
CN109738917B (zh) 一种北斗变形监测中的多路径误差削弱方法及装置
EP2156214B1 (en) Partial search carrier-phase integer ambiguity resolution
WO2017160619A1 (en) Navigation satellite orbit and low latency clock determination with wide-lane and narrow-lane bias corrections
WO2012096773A1 (en) Navigation system and method for resolving integer ambiguities using double difference ambiguity constraints
CN111578935A (zh) 一种利用惯导位置增量辅助gnss模糊度固定的方法
CN113805212B (zh) 一种自适应的gnss载波相位差分滑坡监测方法
WO2017160618A1 (en) Navigation satellite wide-lane bias determination and over-range adjustment system and method
CN113466903B (zh) 一种顾及观测值系统误差的部分模糊度固定算法
WO2017160621A1 (en) Navigation satellite orbit and clock determination with low latency clock corrections
CN112835082A (zh) 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质
CN111077550A (zh) 一种应用于智能终端rtd定位的粗差探测方法及系统
CN112285745B (zh) 基于北斗三号卫星导航系统的三频模糊度固定方法及系统
CN105929430B (zh) 一种gnss零基线参考站间模糊度快速固定方法
CN109212562A (zh) 一种生成载波伪距观测量的方法
CN116299623B (zh) 一种城市复杂场景下的ppp和ins紧组合方法与系统
CN115932923A (zh) 一种基于v2v的共享gnss车辆增强协作定位方法
CN113671551B (zh) Rtk定位解算方法
CN113970772B (zh) 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法
CN112987048B (zh) 一种适用于Andriod智能终端的高精度定位方法、智能终端以及存储介质
CN110208841B (zh) 一种改进的面向非重叠频率的gnss紧组合方法
CN114114357A (zh) 一种多星座动态数据高性能质量控制方法
CN110941002A (zh) 一种自适应抗差的序贯最小二乘精密单点定位方法
Cellmer On-the-fly ambiguity resolution using an estimator of the modified ambiguity covariance matrix for the GNSS positioning model based on phase data
US20240118435A1 (en) Post-processing of global navigation satellite system (gnss) data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant