Precise point position and real-time kinematic (PPP-RTK) positioning method and device

Abstract

A precise point position and real-time kinematic (PPP-RTK) positioning method, including: when direct emission signals broadcast by a multi-system navigation satellite and a low-earth-orbit constellation are detected, determining raw observation data (S11); receiving navigation satellite augmentation information broadcast by the low-earth-orbit constellation, and a low-earth-orbit satellite precise orbit and precise clock difference (S12); using the navigation satellite augmentation information, the low-earth-orbit satellite precise orbit and precise clock difference and the raw observation data for precise point positioning (S13); or when comprehensive ground-based augmentation error correction information is received, using the navigation satellite augmentation information, the low-earth-orbit satellite precise orbit and precise clock difference, the raw observation data and the comprehensive ground-based augmentation error correction information for precise point positioning of ground-based augmentation (S13′). The present application further relates to a precise point position and real-time kinematic (PPP-RTK) positioning device, a computer-readable storage medium and a processor.

Claims

1. A positioning method of Precise Point Positioning-Real Time Kinematic PPP-RTK, comprising: determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a Low Earth Orbit (LEO) constellation is detected via a receiver in a navigation apparatus, wherein the navigation apparatus is arranged at a user equipment; receiving by the receiver navigation satellite augmentation information broadcasted by the LEO constellation, as well as an LEO satellite precision orbit and a precision clock bias; performing Precise Point Positioning of the user equipment by the navigation apparatus using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, and the original observation data; or when ground based augmentation comprehensive error correction information is received, performing ground based augmented Precise Point Positioning of the user equipment by the navigation apparatus using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, the original observation data and ground based augmentation comprehensive error correction information.

2. The method as claimed in claim 1, wherein the determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a Low Earth Orbit (LEO) constellation is detected via a receiver in a navigation apparatus comprises: tracking the direct signal and determining time length for tracking when the direct signal broadcasted by the multi-system navigation satellite and the LEO constellation is detected; and measuring the direct signal to determine the original observation data when the time length for tracking is larger than a preset time length.

3. The method as claimed in claim 1, wherein the multi-system navigation satellite comprises at least one of the following: US Global Positioning System GPS, China Beidou, EU Galileo, Russian Global Navigation Satellite System GLONASS.

4. The method as claimed in claim 1, wherein the navigation satellite augmentation information comprises at least one of the following: a navigation satellite precision orbit and clock bias, a navigation satellite fractional phase bias correction, a LEO satellite fractional phase bias correction, and ionospheric model parameter information.

5. The method as claimed in claim 1, wherein the original observation data comprises at least one of the following: navigation satellite and LEO satellite pseudo-range observation data, navigation satellite and LEO satellite carrier phase observation data, and navigation satellite and LEO satellite Doppler observation data.

6. The method as claimed in claim 1, wherein a processing mode for the Precise Point Positioning comprises at least one of the following: a LEO satellite augmented ambiguity-float solution mode and a LEO satellite augmented ambiguity-fixed solution mode.

7. The method as claimed in claim 1, wherein the ground based augmentation comprehensive error correction information comprises at least one of the following: an undifferenced pseudo-range observation comprehensive error and an undifferenced carrier phase observation comprehensive error.

8. The method as claimed in claim 2, wherein the original observation data comprises at least one of the following: navigation satellite and LEO satellite pseudo-range observation data, navigation satellite and LEO satellite carrier phase observation data, and navigation satellite and LEO satellite Doppler observation data.

9. The method as claimed in claim 7, wherein the ground based augmentation comprehensive error correction information further comprises: ionosphere delays and troposphere delays in a direction of each satellite, and satellite-related channel delays, satellite clock biases, and satellite orbit errors.

10. An electronic device, comprising: at least one processor; and a memory in communication with the at least one processor; wherein the memory stores executable instructions that upon execution cause the at least one processor to perform operations, the operations comprising: determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a Low Earth Orbit (LEO) constellation is detected via a receiver in a navigation apparatus, wherein the navigation apparatus is arranged at a user equipment; receiving by the receiver navigation satellite augmentation information broadcasted by the LEO constellation, as well as an LEO satellite precision orbit and a precision clock bias; performing Precise Point Positioning of the user equipment by the navigation apparatus using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, and the original observation data; or when ground based augmentation comprehensive error correction information is received, performing ground based augmented Precise Point Positioning of the user equipment by the navigation apparatus using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, the original observation data and ground based augmentation comprehensive error correction information.

11. The electronic device as claimed in claim 10, wherein the determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a Low Earth Orbit (LEO) constellation is detected via a receiver in a navigation apparatus comprises: tracking the direct signal and determining time length for tracking when the direct signal broadcasted by the multi-system navigation satellite and the LEO constellation is detected; and measuring the direct signal to determine the original observation data when the time length for tracking is larger than a preset time length.

12. The electronic device as claimed in claim 10, wherein the multi-system navigation satellite comprises at least one of the following: US Global Positioning System GPS, China Beidou, EU Galileo, Russian Global Navigation Satellite System GLONASS.

13. The electronic device as claimed in claim 10, wherein the navigation satellite augmentation information comprises at least one of the following: a navigation satellite precision orbit and clock bias, a navigation satellite fractional phase bias correction, a LEO satellite fractional phase bias correction, and ionospheric model parameter information.

14. The electronic device as claimed in claim 10, wherein the original observation data comprises at least one of the following: navigation satellite and LEO satellite pseudo-range observation data, navigation satellite and LEO satellite carrier phase observation data, and navigation satellite and LEO satellite Doppler observation data.

15. The electronic device as claimed in claim 10, wherein a processing mode for the Precise Point Positioning comprises at least one of the following: a LEO satellite augmented ambiguity-float solution mode and a LEO satellite augmented ambiguity-fixed solution mode.

16. The electronic device as claimed in claim 10, wherein the ground based augmentation comprehensive error correction information comprises at least one of the following: an undifferenced pseudo-range observation comprehensive error and an undifferenced carrier phase observation comprehensive error.

17. The electronic device as claimed in claim 11, wherein the original observation data comprises at least one of the following: navigation satellite and LEO satellite pseudo-range observation data, navigation satellite and LEO satellite carrier phase observation data, and navigation satellite and LEO satellite Doppler observation data.

18. The electronic device as claimed in claim 16, wherein the ground based augmentation comprehensive error correction information further compromises: ionosphere delays and troposphere delays in a direction of each satellite, and satellite-related channel delays, satellite clock biases, and satellite orbit errors.

19. A non-volatile computer readable storage medium storing at least one executable instruction, wherein the executable instruction is configured to cause a processor to perform operations of: determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a Low Earth Orbit (LEO) constellation is detected via a receiver in a navigation apparatus, wherein the navigation apparatus is arranged at a user equipment; receiving by the receiver navigation satellite augmentation information broadcasted by the LEO constellation, as well as an LEO satellite precision orbit and a precision clock bias; performing Precise Point Positioning of the user equipment by the navigation apparatus using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, and the original observation data; or when ground based augmentation comprehensive error correction information is received, performing ground based augmented Precise Point Positioning of the user equipment by the navigation apparatus using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, the original observation data and ground based augmentation comprehensive error correction information.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) By reading the detailed description of the preferred embodiments below, various other advantages and benefits will become clear to those ordinary skills in the art. The attached drawings are only for the purpose of illustrating the preferred embodiments, and are not considered to be restrictions of the present application. And throughout the drawings, the same reference symbols are used to denote the same parts. In the drawings:

(2) FIG. 1 shows a method flowchart of a PPP-RTK positioning method based on a LEO constellation navigation augmentation system of an embodiment of the present application;

(3) FIG. 2 shows a schematic diagram of realization principle of the PPP-RTK positioning method based on the LEO constellation navigation augmentation system of an embodiment of the present application.

(4) FIG. 3 shows a specific method flowchart of a PPP-RTK positioning based on a LEO constellation navigation augmentation system of another embodiment of the present application;

(5) FIG. 4 shows a structural block diagram of a device provided by an embodiment of the present application.

DESCRIPTION OF SYMBOLS FOR MAIN COMPONENTS

(6) 100 navigation constellation 101 Beidou navigation satellite 102 GPS navigation satellite 103 GLONASS navigation satellite 104 Other navigation system satellites 110 LEO augmentation constellation 111 LEO satellite 120 Navigation signals, including navigation satellite navigation signal and LEO satellite navigation signal. 130 Global regions 131 Navigation apparatus using a PPP-RTK positioning method based on LEO constellation navigation augmentation system 140 ground-based augmented region 141 ground-based augmentation monitoring station 142 ground-based augmentation information broadcasting device

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

(7) Hereinafter, exemplary embodiments of the present application will be described in more detail with reference to the accompanying drawings. Although the drawings show exemplary embodiments of the present disclosure, it should be understood that the present application can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present application to those skilled in the art.

(8) FIG. 1 shows a method flowchart of a PPP-RTK positioning method based on a LEO constellation navigation augmentation system of an embodiment of the present application. With reference to FIG. 1, the PPP-RTK positioning method based on a LEO constellation navigation augmentation system provided by an embodiment of the present application includes step S11, step S12 and step S13 (or a step S13′) as follows.

(9) In the step S11, original observation data is determined when a direct signal broadcasted by a multi-system navigation satellite and a LEO constellation is detected.

(10) In the step S12, navigation satellite augmentation information broadcasted by the LEO constellation, as well as a LEO satellite precision orbit and a precision clock bias are received.

(11) In the step S13, Precise Point Positioning is performed by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, and the original observation data.

(12) In the step S13′, when the ground based augmentation comprehensive error correction information is received, ground based augmented Precise Point Positioning is performed by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, the original observation data and ground based augmentation comprehensive error correction information.

(13) In the PPP-RTK positioning method disclosed by the embodiment of the present application, the fast-moving characteristics of LEO satellites is used to broadcast navigation signals, and in the meantime multiple systems are used to increase the number of navigation satellite observations, so that geometric configurations of user observation space are comprehensively improved, and near real-time initialization can be realized globally.

(14) In an embodiment, in a region with a ground-based augmentation monitoring network, the station density of the ground augmentation monitoring network is effectively reduced by this method, and by correcting the user observation errors by receiving comprehensive error information of the ionosphere, troposphere and the like in the current region, real-time initialization can be realized with a unified PPP calculation and processing mode.

(15) FIG. 2 shows a schematic diagram of realization principle of the PPP-RTK positioning method of an embodiment of the present application. The processing procedure thereof is as shown in FIG. 3, which includes: step 201-step 210.

(16) In the step 201, navigation direct signals broadcasted by the multi-system navigation satellite and the LEO constellation are received and the direct signals are captured and tracked.

(17) In the step 202, the navigation direct signals are measured at each epoch, to generate pseudo-range, carrier phase and Doppler observation data.

(18) In the step 203, under the premise of stable signal tracking, telegram parameters of the LEO satellite direct signals are demodulated and navigation satellite augmentation information and LEO satellite precision orbit and clock bias are acquired, wherein the navigation satellite augmentation information includes navigation satellite precision orbit, clock bias, fractional phase bias, and global ionospheric model parameters.

(19) That is, determining original observation data when a direct signal broadcasted by the multi-system navigation satellite and the LEO constellation is detected includes: tracking the direct signal and determining time length for tracking when the direct signal broadcasted by the multi-system navigation satellite and the LEO constellation is detected; measuring the direct signal to determine the original observation data when the time length for tracking is larger than a preset time length.

(20) In the step 204, observation equations is established by using the original observation data, and by taking one kind of satellite navigation system as reference, normalization processing is performed on observation data from other satellite navigation systems and LEO satellites, to obtain observation equations with unified time reference.

(21) In the step 206, if in a non-ground-based-augmented region, observation is corrected by using the navigation satellite augmentation information and the LEO satellite precision orbit and the clock bias

(22) In the step 207, if in a ground-based augmented region, undifferenced comprehensive correction information broadcasted by a ground-based monitoring network is received through a communication link.

(23) In the step 208, error correction parameters of an approximate position of the use relative to each navigation satellite and LEO satellite are calculated according to the received undifferenced comprehensive correction information.

(24) In the step 209, observation data are corrected by using the navigation satellite augmentation information and the LEO satellite precision orbit and the clock bias, and the above calculated error correction parameters.

(25) In the step 210, positioning processing is performed by using a Precise Point Positioning mode, to obtain PPP-RTK positioning, timing and velocity measurement results and carrier phase ambiguity parameters and the like based on the LEO constellation navigation augmentation system.

(26) In the PPP-RTK positioning method provided by the present application, the fast-moving characteristics of LEO satellites is used to broadcast navigation signals, and in the meantime multiple systems are used to increase the number of navigation satellite observations, so that geometric configurations of user observation space are comprehensively improved, and near real-time precise positioning, velocity measurement and timing results can be obtained globally, real-time positioning in the order of centimeter, velocity measurement and timing results can be obtained in the ground-based augmented region, and seamless switching can be performed between the ground-based-augmented region and the rest non-ground-based-augmented region of the world.

(27) The technical solution of the present application will be described in detail below through a specific embodiment.

(28) The main procedure of a rapid PPP processing in the global region based on the LEO constellation navigation augmentation system is as follows.

(29) (1) Establishing Observation Equations by Using Original Observation Data

(30) The original observation data generated from the navigation signals received by a receiver include pseudo-range, carrier phase and Doppler observation data of multiple constellations and multiple frequency points, wherein the observation equations for the pseudo-range and the carrier phase can be denoted as follows:

(31) $\begin{matrix} P_{i, r}^{G, S} = ρ_{r}^{G, S} + {cdt}_{r} - {cdt}^{G, S} + {dTrop}_{r}^{G, S} + \frac{f_{1}^{2}}{f_{i}^{2}} {dIon}_{r}^{G, S} + b_{i, r} - b_{i}^{G, S} + {.Math.}_{i, r}^{G, S} & (1) \\ L_{i, r}^{G, S} = ρ_{r}^{G, S} + {cdt}_{r} - {cdt}^{G, S} + {dTrop}_{r}^{G, S} - \frac{f_{1}^{2}}{f_{i}^{2}} {dIon}_{r}^{G, S} + λ_{i} N_{i, r}^{G, S} + h_{i, r} + {δϕ}_{i, r} - h_{i}^{G, S} - {δϕ}_{i}^{G, S} + δ_{i, r}^{G, S} & (2) \end{matrix}$
in which,
G: denotes the satellite navigation system and the LEO augmentation system;
i: denotes signal frequency identification, i=1, 2, 3;
r, S: respectively denotes receiver and satellite identification;
P.sub.i,r.sup.G,S, L.sub.i,r.sup.G,S: respectively denotes pseudo-range and carrier phase observation values in meters;
ρ.sub.r.sup.G,S: is the geometric distance from the satellite to the receiver;
dt.sub.r, dt.sup.G,S: respectively are clock biases of the receiver and the satellite;
dTrop.sub.r.sup.G,S: is the tropospheric delay;
dIon.sub.r.sup.G,S: is the ionospheric delay on frequency 1;
b.sub.i,r, b.sub.i.sup.G,S: respectively are pseudo-range hardware channel delays of the receiver and the satellite;
N.sub.i,r.sup.G,S: is the integer ambiguity;
h.sub.i,r, δϕ.sub.i,r: respectively are the carrier phase channel delay and the initial phase deviation of the receiver;
h.sub.i.sup.G,S, δϕ.sub.i.sup.G,S: respectively are the carrier phase channel delay and the initial phase deviation of the satellite;
ε.sub.i,r.sup.G,S, δ.sub.i,r.sup.G,S: respectively are the unmodeled residual errors in the pseudo-range and carrier phase observations.

(32) Because the carrier phase channel delay and the initial phase deviation cannot be separated, these two are usually combined and referred to as uncalibrated hardware delay, and at the sides of the receiver and the satellite, are respectively denoted as:
B.sub.i,r=h.sub.i,r+δϕ.sub.i,r (3)
B.sub.i.sup.G,S=h.sub.i.sup.G,S+δϕ.sub.i.sup.G,S (4)

(33) And the observation equation for the carrier phase can be denoted as:

(34) $\begin{matrix} L_{i, r}^{G, S} = ρ_{r}^{G, S} + {cdt}_{r} - {cdt}^{G, S} + {dTrop}_{r}^{G, S} - \frac{f_{1}^{2}}{f_{i}^{2}} {dIon}_{r}^{G, S} + λ_{i} N_{i, r}^{G, S} + B_{i, r} - B_{i}^{G, S} + δ_{i, r}^{G, S} & (5) \end{matrix}$

(35) (2) Constructing Observation Value of Ionosphere-Free Combination

(36) Ionosphere-free combination observations are constructed by using dual frequency observation data, to eliminate the influence of the ionospheric delay of the first order and reduce unknown parameters, and the specific combined model is:

(37) $\begin{matrix} P_{3, r}^{G, S} = \frac{f_{1}^{2} P_{1, r}^{G, S}}{f_{1}^{2} - f_{2}^{2}} - \frac{f_{2}^{2} P_{2, r}^{G, S}}{f_{1}^{2} - f_{2}^{2}} = ρ_{r}^{G, S} + {cdt}_{r} - {cdt}^{G, S} + {dTrop}_{r}^{G, S} + b_{3, r} - b_{3}^{G, S} + {.Math.}_{3, r}^{G, S} & (6) \\ L_{3, r}^{G, S} = \frac{f_{1}^{2} L_{1, r}^{G, S}}{f_{1}^{2} - f_{2}^{2}} - \frac{f_{2}^{2} L_{2, r}^{G, S}}{f_{1}^{2} - f_{2}^{2}} = ρ_{r}^{G, S} + {cdt}_{r} - {cdt}^{G, S} + {dTrop}_{r}^{G, S} + λ_{3} N_{3, r}^{G, S} + B_{3, r} - B_{3}^{G, S} + δ_{3, r}^{G, S} & (7) \end{matrix}$
wherein

(38) $b_{3, r} = \frac{f_{1}^{2} b_{1, r}}{f_{1}^{2} - f_{2}^{2}} - \frac{f_{2}^{2} b_{2, r}}{f_{1}^{2} - f_{2}^{2}}$ $b_{3}^{G, S} = \frac{f_{1}^{2} b_{1}^{G, S}}{f_{1}^{2} - f_{2}^{2}} - \frac{f_{2}^{2} b_{2}^{G, S}}{f_{1}^{2} - f_{2}^{2}}$ $B_{3, r} = \frac{f_{1}^{2} B_{1, r}}{f_{1}^{2} - f_{2}^{2}} - \frac{f_{2}^{2} B_{2, r}}{f_{1}^{2} - f_{2}^{2}}$ $B_{3}^{G, S} = \frac{f_{1}^{2} B_{1}^{G, S}}{f_{1}^{2} - f_{2}^{2}} - \frac{f_{2}^{2} B_{2}^{G, S}}{f_{1}^{2} - f_{2}^{2}}$ $λ_{3} = \frac{c}{f_{1} + f_{2}}$ $N_{3, r}^{G, S} = N_{1, r}^{G, S} + \frac{f_{2}}{f_{1} - f_{2}} (N_{1, r}^{G, S} - N_{2, r}^{G, S})$

(39) Because the pseudo-range hardware delay at the side of the receiver will be absorbed by the receiver clock bias, let
cd{circumflex over (t)}.sub.r=cdt.sub.r+b.sub.3,r
and the above equation becomes
p.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r−cdt.sup.G,S+dTrop.sub.r.sup.G,S+ε.sub.3,r.sup.G,S (8)
L.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r−cdt.sup.G,S+dTrop.sub.r.sup.G,S+λ.sub.3N.sub.3,r.sup.G,S−b.sub.3,r+B.sub.3,r−B.sub.3.sup.G,S+δ.sub.3,r.sup.G,S (9)

(40) When the observation data of the multiple systems are processed in combination, the channel delay of the pseudo-range at the receiver side is absorbed by the clock bias parameter of the receiver and the channel delay is dependent on the signals, which therefore leads to that different systems correspond to different receiver clock biases cd{circumflex over (t)}.sub.r, and then the observation equations corresponding to the LEO satellite and other satellite navigation systems can be rewritten as:
P.sub.3,r.sup.L,S=ρ.sub.r.sup.L,S+cd{circumflex over (t)}.sub.r.sup.L−cdt.sup.L,S+dTrop.sub.r.sup.L,S−d.sub.3.sup.L,S+ε.sub.3,r.sup.L,S (10)
L.sub.3,r.sup.L,S=ρ.sub.r.sup.L,S+cd{circumflex over (t)}.sub.r.sup.L−cdt.sup.L,S+dTrop.sub.r.sup.L,S+λ.sub.3N.sub.3,r.sup.L,S−b.sub.3,r+B.sub.3,r−B.sub.3.sup.L,S+δ.sub.3,r.sup.L,S (11)
P.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r.sup.G−cdt.sup.G,S+dTrop.sub.r.sup.G,S−b.sub.3.sup.G,S+ε.sub.3,r.sup.G,S (12)
P.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r.sup.G−cdt.sup.G,S+dTrop.sub.r.sup.G,S−b.sub.3.sup.G,S+ε.sub.3,r.sup.G,S (12)
L.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r.sup.G−cdt.sup.G,S+dTrop.sub.r.sup.G,S+λ.sub.3N.sub.3,r.sup.G,S−b.sub.3,r+B.sub.3,r−B.sub.3.sup.G,S+δ.sub.3,r.sup.G,S (13)
in which cd{circumflex over (t)}.sub.r.sup.L is the receiver clock bias corresponding to the LEO satellite, and cd{circumflex over (t)}.sub.r.sup.G is the receiver clock bias corresponding to the respective navigation systems. The Frequency Division Multiple Access technology is used in the GLONASS, the receiver pseudo-range channel delays generated by the satellite signals of different frequencies are different, which results in not being completely absorbed by the receiver clock biases. However, if smaller weights are assigned to the GLONASS pseudo-range observation values, these channel delay differences can be approximately considered to be included in the residual errors. Therefore, these variables are no longer reflected in the observation model.

(41) (3) Performing Error Correction by Using Navigation Satellite Augmentation Information Broadcasted by the LEO Satellites and the Model

(42) In the augmentation information broadcasted by the LEO satellites, a unified spatial coordinate reference is used in all of the precise satellite orbit products, and a unified time reference is used in the precise satellite clock bias products. Therefore, there is no problem of inconsistent coordinate reference or time reference in the observation model. At the same time, the precise satellite clock bias products are generated by using the observation value of ionosphere-free combinations, which include the pseudo-range channel delay at the satellite side.

(43) In addition, tropospheric delay can usually be divided into two parts of dry component and wet component. The dry component can be corrected by the model, and the wet component is estimated as the parameters to be estimated. In order to reduce the number of parameters to be estimated, a mapping function can be used to project the oblique delay to the zenith direction, and only one zenith wet delay is estimated.

(44) Let
λ.sub.3{circumflex over (N)}.sub.3,r.sup.G,S=λ.sub.3N.sub.3,r.sup.G,S−(b.sub.3,r−b.sub.3.sup.L,S)+(B.sub.3,r−B.sub.3.sup.G,S) (14)

(45) By using the provided navigation satellite augmentation information, as well as the models such as relativistic effect, earth rotation, and antenna phase center, to correct the observation equations, to eliminate part of unknown parameters, and to ignore residual satellite orbits and clock bias errors, the corresponding observation equations of the LEO satellite and other satellite navigation systems become:
P.sub.3,r.sup.L,S=ρ.sub.r.sup.L,S+cd{circumflex over (t)}.sub.r.sup.L+M.sub.r.sup.L,SZ.sub.r+ε.sub.3,r.sup.L,S (15)
L.sub.3,r.sup.L,S=ρ.sub.r.sup.L,S+cd{circumflex over (t)}.sub.r.sup.L+M.sub.r.sup.L,SZ.sub.r+λ.sub.3{circumflex over (N)}.sub.3,r.sup.L,S+δ.sub.3,r.sup.L,S (16)
P.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r.sup.G+M.sub.r.sup.G,SZ.sub.r+ε.sub.3,r.sup.G,S (17)
L.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r.sup.G+M.sub.r.sup.G,SZ.sub.r+λ.sub.3{circumflex over (N)}.sub.3,r.sup.G,S+δ.sub.3,r.sup.G,S (18)
wherein M.sub.r.sup.L,S, M.sub.r.sup.G,S respectively are the mapping functions corresponding to the LEO satellite and other satellite navigation systems, and Z.sub.r is the tropospheric zenith wet delay.

(46) (4) Linearizing the Observation Equations

(47) By performing Taylor expansion at the approximate position of the receiver, and discarding the second-order terms, the linearized observation equations are obtained as follows:
P.sub.3,r.sup.L,S=R.sub.r,0.sup.L,S−lΔx.sub.r−mΔy.sub.r−nΔz.sub.r+cd{circumflex over (t)}.sub.r.sup.L+M.sub.r.sup.L,SZ.sub.r+ε.sub.3,r.sup.L,S (19)
L.sub.3,r.sup.L,S=R.sub.r,0.sup.L,S−lΔx.sub.r−mΔy.sub.r−nΔz.sub.r+cd{circumflex over (t)}.sub.r.sup.L+M.sub.r.sup.L,SZ.sub.r+λ.sub.3{circumflex over (N)}.sub.3,r.sup.L,S+δ.sub.3,r.sup.L,S (20)
P.sub.3,r.sup.G,S=R.sub.r,0.sup.L,S−lΔx.sub.r−mΔy.sub.r−nΔz.sub.r+cd{circumflex over (t)}.sub.r.sup.G+M.sub.r.sup.G,SZ.sub.r+ε.sub.3,r.sup.G,S (21)
L.sub.3,r.sup.G,S=R.sub.r,0.sup.L,S−lΔx.sub.r−mΔy.sub.r−nΔy.sub.r+cd{circumflex over (t)}.sub.r.sup.G+M.sub.r.sup.G,SZ.sub.r+λ.sub.3{circumflex over (N)}.sub.3,r.sup.G,S+δ.sub.3,r.sup.G,S (22)
in which

(48) $l = \frac{(x^{s} - x_{r, 0})}{\sqrt{{(x^{s} - x_{r, 0})}^{2} + {(y^{s} - y_{r, 0})}^{2} + {(z^{s} - z_{r, 0})}^{2}}}$ $m = \frac{(y^{s} - y_{r, 0})}{\sqrt{{(x^{s} - x_{r, 0})}^{2} + {(y^{s} - y_{r, 0})}^{2} + {(z^{s} - z_{r, 0})}^{2}}}$ $n = \frac{(z^{s} - z_{r, 0})}{\sqrt{{(x^{s} - x_{r, 0})}^{2} + {(y^{s} - y_{r, 0})}^{2} + {(z^{s} - z_{r, 0})}^{2}}}$
(x.sup.s,y.sup.s,z.sup.s) are the coordinates of the LEO satellite and navigation satellite precision orbit, and (x.sub.r,0,y.sub.r,0, z.sub.r,0) are the approximate position of the receiver. And then the observation equations can be simplified as:
V=AΔX+L (23)
in which V is the observation residual error, A is the coefficient matrix, ΔX is the unknown vector including receiver coordinate correction, receiver clock bias, tropospheric zenith wet delay, carrier phase ambiguity, L is the calculation vector.

(49) (5) Performing Parameter Estimation and Ambiguity Fixation Processing

(50) Kalman filtering is used for comprehensive PPP processing. In the filtering, it is necessary to provide a suitable observation value random model and a state vector dynamic model. A random model describes the statistical characteristics of the observation values, and is usually represented by the variance-covariance matrix of the observation values. It can be known from the observation equations that the observation values of ionosphere-free combination are a linear combination of the original observation values. Assuming that the observation values at different frequencies are uncorrelated, the initial variance of the observation values of ionosphere-free combination can be calculated by the law of error propagation. The specific variance can be defined as a function of the initial variance and the satellite altitude angle. Assuming that the observation values of different satellites and different systems are uncorrelated, and that the observation values of different types, that is, the pseudo-range and phase observation values, are uncorrelated, the variance-covariance matrix of the observation values can be obtained.

(51) For the dynamic model of the state vector, the static receiver coordinates can be represented as constants, and the dynamic receiver coordinates and the receiver clock biases can be represented as random walks or first-order Gaussian-Markov processes, and the tropospheric zenith wet delay can be represented as random walk processes, the carrier phase ambiguity parameter can be represented as a constant, and then the state equation is obtained:
X.sub.k=Φ(t.sub.k,t.sub.k-1)X.sub.k-1+w.sub.k-1 (24)
wherein, X is the parameter such as receiver coordinate correction to be estimated, receiver clock bias and the like, Φ is the state transition matrix, and w.sub.k-1 is the state transition noise. By integrating the observation equations and the state equation, the standard Kalman filtering process can be used for parameter estimation. Here, since the satellite fractional phase bias correction is not performed, only the carrier phase ambiguity-float solution result is obtained. If the satellite fractional phase bias contained in the LEO satellite augmentation information is further used to correct the observation equations, the integer characteristic of the ambiguity can be restored, the ambiguity fixation can be realized, and the carrier phase ambiguity-fixed solution result can be obtained, which further shortens the time for initialization and improves the accuracy of positioning, velocity measurement and timing.

(52) Due to the addition of the observation data of the navigation direct signals of the LEO constellation, observation geometry for the user is greatly improved by the fast-moving characteristics of the low-orbit satellites, thereby greatly reducing the time for PPP initialization.

(53) The main procedure of the PPP-RTK processing augmented based on LEO constellation navigation in the ground-based augmented region is:

(54) In the ground-based augmented region, all reference stations will be divided into several triangular subnets using the Delaunay method, and the comprehensive error correction information of each visible satellite will be constructed for each subnet according to the network RTK method based on undifferenced correction numbers, respectively, in which the ionosphere and troposphere in the direction of each satellite, and satellite-related channel delays, satellite clock biases, and satellite orbit errors are included, which are represented as:

(55) $\begin{matrix} {omc}_{P_{i, r}}^{G, S} = - {cdt}^{G, S} - {dOrb}_{r}^{G, S} + {dTrop}_{r}^{G, S} + \frac{f_{1}^{2}}{f_{i}^{2}} {dIon}_{r}^{G, S} - b_{i}^{G, S} & (25) \\ {omc}_{L_{i, r}}^{G, S} = - {cdt}^{G, S} - {dOrb}_{r}^{G, S} + {dTrop}_{r}^{G, S} - \frac{f_{1}^{2}}{f_{i}^{2}} {dIon}_{r}^{G, S} - B_{3}^{G, S} & (26) \end{matrix}$
omc.sub.P.sub.i,r.sup.G,S, omc.sub.L.sub.i,r.sup.G,S respectively denote the pseudo-range and carrier phase comprehensive error correction information.

(56) Plane fitting modeling is performed by the receiver on the comprehensive error correction information of at least 3 ground-based augmentation stations at the periphery according to the approximate position, and the pseudo-range and carrier phase observation values of the user are refined by using the local error correction information obtained by interpolation. After correction, the corresponding observation equations of the LEO satellites and other satellite navigation systems can be written as:
P.sub.3,r.sup.L,S=ρ.sub.r.sup.L,S+cd{circumflex over (t)}.sub.r.sup.L+ε.sub.3,r.sup.L,S (27)
L.sub.3,r.sup.L,S=ρ.sub.r.sup.L,S+cd{circumflex over (t)}.sub.r.sup.L+λ.sub.3N.sub.3,r.sup.L,S−b.sub.3,r+B.sub.3,r+δ.sub.3,r.sup.L,S (28)
P.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r.sup.G+ε.sub.3,r.sup.G,S (29)
L.sub.3,r.sup.G,S=ρ.sub.r.sup.G,S+cd{circumflex over (t)}.sub.r.sup.G+λ.sub.3N.sub.3,r.sup.G,S−b.sub.3,r+B.sub.3,r+δ.sub.3,r.sup.G,S (30)

(57) The unknown variables in the equations at this time include receiver position coordinates, receiver clock biases, receiver channel delays and carrier phase ambiguity parameters. The receiver clock biases and the channel delays can be further eliminated by using the single difference between satellites. And then receiver position estimation and carrier phase ambiguity fixation can be performed by using linearization processing strategy and parameter estimation method introduced above.

(58) Due to the addition of the observation data of the navigation direct signals of the LEO constellation, observation geometry for the user is greatly improved by the fast-moving characteristics of the LEO satellites, so that under the same initialization time and positioning accuracy requirements, the station density of the ground augmented monitoring network can be greatly reduced, thereby reducing the construction cost of the monitoring network.

(59) The LEO constellation augmented fast velocity measurement, calculation process of timing for global region and ground-based augmented region is similar to the positioning process, and will not be repeated here.

(60) In the PPP-RTK positioning method provided by the embodiment of the present application, navigation direct signals are broadcasted through LEO satellites, geometric configurations of user observation space are comprehensively improved by using the fast-moving characteristics of the LEO satellites, and the time for initialization of the PPP of the user can be shortened to quasi real-time. In the ground-based augmented region, the time for initialization is further shortened to real-time by correcting the user observation errors with further receiving the comprehensive error information of the current region, and adopting the calculation processing mode unified with the PPP. Through this method, the station density of the ground augmented monitoring network can be effectively reduced and the seamless switching of multi-system PPP and RTK services augmented by LEO constellation can be realized through a unified model.

(61) FIG. 4 is a structural block diagram of an electronic device provided by the present embodiment. The device provided by the present embodiment may include a processor 401 and a memory 403, and may also include a communications interface 402 and a bus 404. Wherein the processor 401, the communication interface 402, and the memory 403 can communicate with each other through the bus 404. The communication interface 402 can be used for information transmission. The processor 801 may call logic instructions in the memory 403 to execute the positioning method of Precise Point Positioning-Real Time Kinematic PPP-RTK of the forgoing embodiment, including: determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a LEO constellation is detected; receiving navigation satellite augmentation information broadcasted by the LEO constellation, as well as a LEO satellite precision orbit and a precision clock bias; performing Precise Point Positioning by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, and the original observation data; or when the ground based augmentation comprehensive error correction information is received, performing ground based augmented Precise Point Positioning by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, the original observation data and ground based augmentation comprehensive error correction information.

(62) The embodiment of the present application also provides a non-volatile storage medium, which includes a stored program, wherein the method provided by the foregoing embodiment is executed when the program is running, and the method includes:

(63) determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a LEO constellation is detected; receiving navigation satellite augmentation information broadcasted by the LEO constellation, as well as a LEO satellite precision orbit and a precision clock bias; performing Precise Point Positioning by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, and the original observation data; or when the ground based augmentation comprehensive error correction information is received, performing ground based augmented Precise Point Positioning by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, the original observation data and ground based augmentation comprehensive error correction information.

(64) The present application also provides a processor, which is configured to run a program, wherein the program is configured to execute the method provided by the foregoing embodiment on the processor, and the method includes:

(65) determining original observation data when a direct signal broadcasted by a multi-system navigation satellite and a LEO constellation is detected; receiving navigation satellite augmentation information broadcasted by the LEO constellation, as well as a LEO satellite precision orbit and a precision clock bias; performing Precise Point Positioning by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, and the original observation data; or when the ground based augmentation comprehensive error correction information is received, performing ground based augmented Precise Point Positioning by using the navigation satellite augmentation information, the LEO satellite precision orbit and the precision clock bias, the original observation data and ground based augmentation comprehensive error correction information.

(66) Through the above description of the embodiments, those skilled in the art can clearly understand that multiple embodiments can be implemented by means of software plus a necessary general hardware platform, and of course, they can also be implemented by hardware. Based on such understanding, the above technical solution essentially or the part of the above technical solution that contributes to the related technology can be embodied in the form of a software product. The computer software product can be stored in a computer-readable storage medium, such as a Read-Only Memory (ROM)/Random Access Memory (RAM), magnetic disks, optical disks, etc., including several instructions to make a computer device (which may be a personal computer, server, or network device, etc.) execute multiple embodiments or methods described in some parts of the embodiments.

(67) In addition, those skilled in the art can understand that, although some embodiments herein include certain features included in other embodiments but not other features, the combination of features of different embodiments means that they are within the scope of the present application and form different embodiments. For example, in the following claims, any one of the claimed embodiments can be used in any combination.

Precise point position and real-time kinematic (PPP-RTK) positioning method and device

Assignee

Inventors

Cpc classification

Classification Explorer

G01S19/073

PHYSICS

Classification Explorer

G01S19/425

PHYSICS

Classification Explorer

G01S19/115

PHYSICS

Classification Explorer

G01S19/44

PHYSICS

Classification Explorer

G01S19/23

PHYSICS

Classification Explorer

G01S19/46

PHYSICS

Classification Explorer

G01S19/52

PHYSICS

Classification Explorer

G01S19/37

PHYSICS

Classification Explorer

G01S19/071

PHYSICS

International classification

Classification Explorer

G01S19/44

PHYSICS

Classification Explorer

G01S19/23

PHYSICS

Classification Explorer

G01S19/37

PHYSICS

Abstract

Claims

Description