EP3189349A1 - Method of collaborative determination of positioning errors of a satellite-based navigation system - Google Patents
Method of collaborative determination of positioning errors of a satellite-based navigation systemInfo
- Publication number
- EP3189349A1 EP3189349A1 EP15763847.9A EP15763847A EP3189349A1 EP 3189349 A1 EP3189349 A1 EP 3189349A1 EP 15763847 A EP15763847 A EP 15763847A EP 3189349 A1 EP3189349 A1 EP 3189349A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- receivers
- receiver
- data
- error
- contributing
- 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 description 48
- 238000004364 calculation method Methods 0.000 claims description 28
- 238000012545 processing Methods 0.000 claims description 23
- 230000002123 temporal effect Effects 0.000 claims description 18
- 238000004891 communication Methods 0.000 claims description 9
- 230000004927 fusion Effects 0.000 claims description 8
- 230000002706 hydrostatic effect Effects 0.000 claims description 6
- 230000001413 cellular effect Effects 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 2
- 230000004807 localization Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 abstract description 10
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 24
- 230000009471 action Effects 0.000 description 15
- 239000005433 ionosphere Substances 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 239000005436 troposphere Substances 0.000 description 5
- 230000003416 augmentation Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 239000000969 carrier Substances 0.000 description 4
- 230000000875 corresponding effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000037361 pathway Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000029305 taxis Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
-
- 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/40—Correcting position, velocity or attitude
Definitions
- the present invention applies to systems for aiding satellite-based navigation. More precisely, the aim of the invention is to produce information making it possible to determine the positioning errors resulting from the use of a satellite- based navigation signal.
- the first constellation of satellites transmitting positioning signals was put in place for military applications by the American State (Global Positioning System or GPS) at the start of the 1980s. Since then, GPS signals have been used by professional civil applications (management of fleets of lorries, aids to aerial navigation, geodesic surveys, etc.), and henceforth for general-public applications (automobile navigation with onboard terminals and pedestrian navigation with terminals of smartphone type). Other constellations were put in place by the Russian State (GLONASS) and the Chinese State (Ba ' idou). A constellation of European satellites (Galileo) is undergoing deployment. Generally, these navigation systems are designated by the acronym GNSS (Global Navigation Satellite Systems).
- GLONASS Russian State
- Ba ' idou Chinese State
- GNSS Global Navigation Satellite Systems
- the basic principle of aiding satellite-based positioning and navigation is the calculation by a receiver furnished with electronic processing circuits of position, velocity and time (PVT) data on the basis of electromagnetic signals of centimetric wavelength transmitted by satellites in orbit.
- PVT position, velocity and time
- the calculation of the PVT data by a receiver on the basis of the signals of the satellites is affected by numerous errors of various types: impact of the crossing by the electromagnetic signals of the various layers of the atmosphere (troposphere, ionosphere), errors due to the reflections of the signals on objects in the neighbourhood of the receiver (multipaths), clock errors, errors of the electronic processing circuits, etc.
- P(Y) code of GPS reserved carriers
- Specific means of multi-sensor processing and fusion are furthermore generally envisaged in order to guarantee precision and integrity of measurements intended for critical uses. But these solutions are restricted and expensive.
- the errors due to the crossing of the ionosphere by the signals transmitted by the radio navigation satellites weigh heavily in the global toll of the positioning errors (4 standard deviations, according to the publication by A Angrisano et al., "Ionospheric models comparison for single-frequency GNSS positioning", ENC 201 1 , 12/201 1 , http://panq.uniparthenope.it/node/64).
- Several types of techniques for correcting these ionospheric errors may be employed in the state of the art.
- a dual-frequency receiver can thus use a linear combination of the pseudo- distances calculated on the basis of the signals of each of the frequencies.
- the ionospheric error being strongly correlated with the frequencies, it will be able to be eliminated by said combination.
- dual-frequency receivers are not yet widespread among the general public. Furthermore, the convergence time to a stabilized measurement is relatively long (possibly reaching as much as half an hour).
- a method making it possible to guarantee at one and the same time fast convergence, precision and integrity is the acquisition of specialized signals containing corrections calculated on the basis of the differences between the known positions of fixed stations belonging to a network and the positions calculated on the basis of the navigation signals of a GNSS constellation.
- augmentation systems Small Based Augmentation Systems or SBAS
- SBAS Stelite Based Augmentation Systems
- the operational systems comprise notably EGNOS in Europe (European Geostationary Navigation Overlay Service), and WAAS in the United States (Wide Area Augmentation System).
- RTK Real Time Kinematics
- a limitation which is common to the differential approaches of SBAS and RTK type is that of requiring the use of reference stations whose position is known very precisely and of calculation algorithms which process the entirety or the major part of the errors with equivalent precision. This limitation restricts the access of users of general-public terminals furnished with standard GNSS signal acquisition capabilities to greater positioning precision.
- the invention affords a solution to this problem by generating a map of the ionospheric errors of sufficient precision on the basis of terminals for receiving GNSS signals whose position may not be known precisely, and optionally whose onboard processing capability may be limited.
- the invention discloses a method of determination by a computer server of errors in calculating a position calculated by user receivers of positioning signals received from at least one constellation of satellites in a determined geographical zone, said method comprising: determining a list of contributing receivers Ri of positioning signals situated in said geographical zone Zj, a position of said contributing receivers not being known a priori; extracting, for each of the contributing receivers Ri, code Ci and phase ⁇ ; data of at least one channel of the receiver; a step of calculating an atmospheric error li applicable in the zone to the contributing receivers Ri , code Ci and phase ⁇ ; data of at least one channel of the contributing receiver; and calculating an atmospheric error l(Zj) applicable in the zone Zj, on the basis of the errors li applicable to the contributing receivers Ri.
- Advantageously calculating the atmospheric error li consists in calculating only the ionospheric error for at least one contributing receiver Ri.
- calculating the atmospheric error li consists in calculating at least the tropospheric error for at least one contributing receiver Ri.
- the tropospheric error for the at least one contributing receiver Ri is calculated as the sum of a first delay termed the wet delay (SWD) and of a second delay termed the hydrostatic delay SHD.
- calculating the atmospheric error I (Zj) receives as input the ionospheric errors of the contributing receivers Ri and a tropospheric error applicable in the geographical zone Zj originating from an external source of provisioning information.
- determining at least one part of the list of contributing receivers Ri (2221 , 2222) situated in said geographical zone uses data of localization of said receivers in relation to at least one base station of a cellular radiocommunications network.
- the method of the invention furthermore comprises classifyinng pairs of code Ci and phase ⁇ ; data of channels of a receiver Ri, said classifying step using a weighted combination of criteria comprising criteria representative of one or more of a signal-to-noise ratio of the channel or an elevation of the satellite axis of the channel.
- calculating an ionospheric error li applicable to the receiver Ri comprises subtracting said code Ci and phase ⁇ ; data.
- the method of the invention further comprises calculating a quality index IQ(li) of the ionospheric error li applicable to the receivers Ri.
- the quality index IQ(li) is calculated by a weighted combination of criteria comprising criteria representative of one or more of a signal- to-noise ratio of a channel of the receiver Ri, an elevation of a satellite axis acquired by the receiver Ri, a number of satellite axes acquired.
- the quality index IQ(li) further depends on parameters characterizing one or more of a hardware or a software configuration of the receiver Ri.
- the errors li of the receivers Ri whose quality index IQ(li) is less than a threshold value Threshl have a zero weighting in calculating an ionospheric error l(Zj) applicable in the zone Zj.
- the ionospheric error I (Zj) applicable in the zone Zj is a linear combination of the errors li weighted by the quality indices IQ(li).
- the method of the invention further comprises filtering values I (Zj) by a model of ionospheric errors.
- the method of the invention further comprises determining a confidence indicator of the ionospheric error I (Zj) .
- the confidence indicator of the ionospheric error I (Zj) is dependent on a number Nj of contributing receivers Ri whose quality index IQ(li) is greater than Threshl .
- the ionospheric error I (Zj) is replaced, when its confidence indicator is less than a predetermined threshold value Thresh2, by a value l(Z) calculated by one or more of a spatial or a temporal interpolation within a set of values of ionospheric errors calculated for neighbouring geographical zones or periods.
- the method of the invention further comprises making values of ionospheric errors I (Zj) available to the user receivers.
- the method of the invention further comprises acquiring by the user receivers of a position determined by a server (230).
- the method of the invention further comprises calculating a position corrected by fusion of the PVT data of the user receiver and of the values of ionospheric errors l(Zj).
- the invention also discloses a method for correcting a position calculation performed by a user receiver of positioning signals received from at least one constellation of satellites in a determined geographical zone, said method comprising at least: a acquiring from a server an atmospheric error l(Z) calculated by a procedure of fusion by said server of atmospheric errors determined from a plurality of contributing receivers situated in the geographical zone, a position of said contributing receivers not being known a priori; fusioning the atmospheric error acquired with the PVT point calculated by the user receiver.
- the invention also discloses a collaborative system for aiding the positioning of user receivers of positioning signals received from at least one constellation of satellites in at least one geographical zone, said system comprising a plurality of hardware and computer code elements configured to: determine a list of contributing receivers Ri of positioning signals situated in said geographical zone Zj, a position of said contributing receivers not being known a priori; extract, for each of the contributing receivers Ri, code Ci and phase ⁇ ; data of at least one channel of the receiver; calculate an atmospheric error li applicable to the contributing receivers Ri from said code Ci and phase ⁇ ; data; calculate an atmospheric error l(Zj) applicable in the zone Zj, from the errors li of the contributing receivers Ri; make available to the user receivers at least one result of the calculation of atmospheric error in said at least one geographical zone.
- a processor of at least one contributing receiver Ri is configured to perform the calculation of atmospheric error li applicable to said contributing receiver Ri from said code Ci and phase ⁇ ; data, the result of said calculation being transmitted to a server (230) through a communication link.
- code Ci and phase ⁇ ; data of at least one channel of at least one contributing receiver Ri are transmitted to a server (230) through a communication link, a processor of said server being configured to perform the calculation of atmospheric error li applicable to said contributing receiver Ri from said code Ci and phase ⁇ ; data.
- the invention also discloses a server of positioning data, said server comprising: circuit logic configured to determine a list of contributing receivers Ri situated in a geographical zone; processing logic configured to acquire with a determined frequency temporal sequences of data transmitted by the contributing receivers Ri of said list, said temporal sequences comprising one or more of: i) tcode and phase of at least one channel of each of the contributing receivers Ri, or, ii) a datum calculated from said code and said phase; data logic configured to calculate a temporal sequence of atmospheric errors applicable in said geographical zone; a link to make available to user receivers said temporal sequence of atmospheric errors.
- the invention also discloses a receiver of positioning signals, comprising: a memory for storing temporal sequences of data comprising: i) one or more of code and phase of at least one channel of the receiver, or, ii) a datum calculated from the code and of the phase, and; iii) at least one quality index of said data, said index being chosen from a group comprising a signal-to-noise ratio of the channel, an elevation of a satellite axis corresponding to said channel, and a number of satellite axes acquired; and a communication link configured to teletransmit said temporal sequences stored in a memory of the receiver to a pre-identified remote server with a predetermined frequency.
- the invention also discloses a receiver of positioning signals, comprising: a communication link configured to tele-receive with a predetermined frequency temporal sequences of data of the server of positioning data according to the invention; and circuit logic configured to calculate a corrected position from a PVT point calculated by the receiver and said temporal sequences of data.
- Another advantage of the invention is to allow faster convergence of the calculation of precise position (quasi-instantaneous convergence can be envisaged) by virtue of the information which would be shared about the ionosphere and optionally the troposphere.
- the invention in certain embodiments, it is possible to determine at one and the same time, separately or in combination, the ionospheric errors and the tropospheric errors (each of these errors constituting the atmospheric errors).
- Another advantage of the invention is that of providing a calculation of position corrected of the ionospheric errors, making it possible to obtain, under most conditions of use, a global precision of positioning of the same order of magnitude as that provided by SBAS systems, when the number of terminals in a given geographical zone is sufficient (without however the same integrity guarantees).
- Another advantage of the invention is that of authorizing the recombination of the corrections of ionospheric errors obtained of the invention with positioning data of variable precision obtained moreover so as to obtain a global precision level which will be able to depend, on the one hand, on the configuration of the receiver possibly allowing it to acquire and process the data of this or that precise positioning service, and on the other hand on the reception conditions just where said receiver is located.
- the position data used to initialize the position calculation in a GNSS receiver can for example be data calculated by another means provided either by a radiocommunication, cellular or WiFi-terminal-based network, or by a cartographic means.
- Another advantage of the invention, coupled with these initialization means, is then that of allowing faster acquisition of the GNSS signals by a receiver entering a given geographical zone or exiting a masking zone.
- the invention is sufficiently versatile for its implementation to cope with several types of architectures making it possible to obtain a contribution of correction data according to different modalities of receivers having variable onboard calculation capabilities.
- FIG. 1 represents a basic view of a system of RTK type according to the prior art
- FIG. 2 represents an architecture diagram for implementing the invention, in several of its embodiments
- FIG. 3 represents a general flowchart of the processing actions for determining the ionospheric errors, according to several embodiments of the invention
- FIG. 4 represents a general flowchart of the processing actions for a use of the data of a server of positioning data, according to several embodiments of the invention.
- Figure 1 represents a basic view of an RTK system according to the prior art.
- a user 1 00 receives positioning signals for at least one constellation of satellites, 1 14, 1 19, 121 , 1 23.
- the positioning signals for said satellites are also received by a reference station 1 30.
- the user 1 00 can be mobile.
- the reference station 1 30 must be fixed, of precisely known position and situated in a radius of about 1 0 to 20 km from the user (see http://www.navipedia.net/index.php/RTK Fundamentals).
- An algorithm of RTK type can be implemented by using the code, the phase or a combination of the two, if appropriate with variable weightings.
- the precision is better when the phase is used, since its measurement is less noisy than that of the code, at the price of greater complexity, notably in order to remove the cycle ambiguity.
- p is the geometric expression of the pseudo-distance calculated as a function of the cartesian coordinates of the satellite (x Sat , y Sat , z Sat ) and of the receiver (3 ⁇ 4 , y ft , z ft ) ;
- / is the delay of the signal due to crossing the ionosphere
- T r is the delay of the signal due to crossing the troposphere
- c is the velocity of light in vacuo
- b s is the bias of the clock of the receiver with respect to the GNSS reference time
- b Sat is the bias of the clock of the satellite with respect to the GNSS reference time
- N is the ambiguity of the phase of the carrier
- /l is the nominal wavelength of the carrier
- ⁇ ⁇ summarizes the measurement noise, notably the multipaths.
- the precision of the position measurement of the receiver of the user 100 will be of the same order of magnitude as the precision of the various terms of equation (3), or of the equivalent equations for the code and code/phase combination variants, with a precision equivalent to the precision of the positioning which must be provided.
- To obtain metric precision it is thus necessary to ascertain the position of the reference station with identical precision, thereby in practice excluding the possibility of it being mobile. It is also necessary to calculate the various error terms with the same precision, this being very difficult for the ionospheric error, having regard to the multiple factors which affect it.
- Figure 2 represents an architecture diagram for implementing the invention, in several of its embodiments.
- the invention can be implemented with the satellites of the constellation of the prior art, 1 14, 1 19, 121 , 123.
- Geographical zones 210, 220 are defined within which the measurements of receivers 21 1 1 , 21 12, 21 13, 221 1 , 2212, 2213, notably, which are active in said geographical zones, contribute to a global calculation of ionospheric errors.
- the geographical zones 210, 220 have a mesh which may be of the order of a few kilometres (for example between 1 and 10 km).
- a server 230 managed by a provider of positioning data 240 can receive data from the receivers and transmit same thereto.
- the receivers can be of all types, mono-frequency or multifrequency, able to pick up the signals of a single GNSS constellation or of several, furnished or not with antennas suitable for attenuating the effects of the multipaths, and aided or not by an inertial hybridization or a mapping. They are also coupled to a module allowing them to transmit and receive data by RF pathway.
- simple smartphones comprising a standard GPS microchip, can be integrated into the measures.
- the receivers 21 1 1 , 21 12, 221 1 , 2212 may be solely contributors, that is to say may provide the server 230, and/or other receivers, with raw or summarized data making it possible to calculate the map of the ionospheric error in the zones 210, 220, without being users thereof. This may in particular be the case for vehicles of commercial or administrative fleets (public transport, taxis, ambulances, postal vehicles, lorries, etc.) equipped with means of precise positioning and receiving precision data.
- the receivers 21 13, 2213 may be solely users, that is to say may receive positioning data without contributing their own position data.
- the receivers may also be contributors and users simultaneously or successively.
- the receivers must be able to provide as a minimum the code and the phase of the raw navigation signals. If the position of the receiver must be determined on the basis of the GNSS signal, the minimum number of satellite axes must be four. If the position of the receiver is determined by another means, as indicated hereinbelow, this is not necessary. A selection of the best signals will nonetheless be carried out. It is also possible to generate a summary signal, according to the modalities described hereinbelow as commentary to Figure 3. The raw signals or the summary signal are thereafter transmitted by RF pathway to the server 230.
- the zones 210 and 220 differ from one another in that the second comprises base stations 2221 , 2222, 2223 of a cellular radiocommunications network which are equipped with transmit relay antennas. These antennas may belong to different operators, or optionally be rented by operators that do not have their own network of base stations.
- the operator of a radiocommunications network is able to ascertain at any instant the base station to which an active terminal of one of its subscribers is attached. The precision of the position of the subscriber therefore varies as a function of the density of the cells, therefore from a hundred or so metres to a few kilometres.
- the provider of positioning data 240 is therefore able to ascertain the receivers that are active in the geographical zone 220, thereby making it possible to initialize the method of the invention, without the precise position of the receiver being known, as indicated hereinbelow. In this case, the transmission or the calculation of a GNSS position is therefore not necessary. It must on the other hand be possible to identify receivers that are located in the geographical zone 210, not covered by cells of a radiocommunications network of an operator with whom the positioning operator has negotiated an agreement to share data.
- Figure 3 represents a general flowchart of the processing actions for determining the ionospheric errors according to several embodiments of the invention.
- the contributors to the cooperative production of the positioning data according to the invention by the provider 240 are furnished with communicating positioning receivers, 21 1 1 , etc., designated in Figure 3 by Ri, having the characteristics indicated hereinabove.
- the contributors are registered by the provider 240, which also registers the identifiers of their GNSS receiver and of their radiocommunication receiver, as well as, optionally, the configuration characteristics of said receivers.
- the contributors must give, at the moment of their registration, their consent to the permanent collection of their location data for the needs of the service. In return for this consent, the data provider 240 will make the necessary legal and contractual undertakings in compliance with the provisions of the countries in which he operates regarding the protection of personal data. Optionally, the data provider will be able to collect the geo-location coordinates of the cellular network base stations onto which the receiver Ri is successively locked.
- processing actions for carrying out the invention are, on the one hand, implemented for each of the receivers Ri of the contributors located in a zone Zj (processing actions 310), and on the other hand are implemented by the entirety of the zones Zj served by the provider of positioning data (processing action 320).
- the positioning data produced by the processing actions 320, and optionally by external processing actions, are thereafter made available to the users and contributors/users of the service (processing action 330).
- each receiver Ri can be carried out according to variants of the invention, either in the receiver itself, or in the server 230. Likewise certain processing actions applicable in a concomitant manner to all the receivers Ri of a zone Zi can also be carried out either on each of the receivers, or on the server.
- a reasonable frequency of the processing actions for updating the data of the server 230, and therefore for extracting and/or calculating the data required at the level of the receivers and/or of the server, is of the order of a few minutes (between 1 and 5 min, for example, for the information relating to the ionosphere, and probably 1 min for the information about the troposphere).
- a position of the receiver Ri is determined; this position may be that provided to the receiver by the GNSS constellation or a position provided by the contributor's radiocommunications operator, as indicated above; in the first case, it must be transmitted to the server 230 so that the receiver Ri is assigned to the zone Zj in the time interval considered;
- the code and phase of one or more carriers on a certain number of satellite axes are extracted from the GNSS receiver Ri; the code and phase of one or more carriers on the satellite axis giving the best quality (for example, the best signal-to-noise ratio and the most significant elevation) are preserved so as to be processed in the course of the following steps;
- step 3120 can be processed in the receiver or transmitted raw to the server 230 so as to be processed therein; the processing 3130 consists in calculating the ionospheric error at the position Pi by performing the sum of the code and of the phase of the positioning signal for the satellite axis retained for the carrier or carriers, as explained further on;
- the quality data corresponding to the best satellite axis, as well as summary data, such as, for example, the number of satellite axes locked on, are extracted so as to be processed, either in the receiver, or on the server in the course of a step 3150 of calculating a synthetic quality index; if it is produced by the receiver, the synthetic quality index IQ(li) of the receiver Ri is transmitted to the server 230; by default, as a variant, the raw quality data themselves are transmitted to the server.
- e c is the measurement noise relating to the code.
- the integer ambiguity ⁇ can be neglected (it is of the order of 10 to 50 cm) in the case of a mono-frequency receiver.
- the residual value of the ambiguity ⁇ will be processed by the positioning filter of the receiver.
- the integer ambiguity is resolved by the standard processing of the receiver.
- the estimator for the measurement noise in the code and the phase is Gaussian noise of zero mean, thereby making it possible, according to the invention, to delete said noise from the calculations without loss of precision in the calculation of an ionospheric correction applicable in the zone Zj, on condition that it is possible to use a sufficient number of measurements of receivers in the zone Zj.
- a step 3210 by eliminating from the subsequent calculations all the receivers whose quality index IQ(li) is less than an acceptance threshold; this threshold is defined on the basis of the criteria indicated above, or of others specific to the zone Zj or to the measurement period; likewise, the acceptability level may be different as a function of the characteristics of the zone and of the period; by way of example, a signal-to- noise ratio of less than 30 dB or an elevation of the satellite axis of less than 15° will make it possible to discard the corresponding measurements; another example could be an indicator as to the multipath environment quality detected in the signal;
- step 3220 it may, in the course of a step 3220, be decided that no ionospheric error measurement point will be provided for the zone Zj because of an insufficiency in the number of acceptable measurements in said zone; in this case, it will be possible despite everything to provide a point l(Z) calculated, in the course of a step 3230, on the basis of an interpolation of the measurements of the neighbouring zones or epochs, of interpolation techniques discussed further on;
- the ionospheric error l(Zj) is calculated, in the course of a step 3240, by taking the mean of the li output by equation (5) of the receivers Ri that are active in the zone Zj, whose measurements satisfy the quality indicator threshold.
- the measurements aggregated during step 3240 can be weighted by their quality indicator.
- a second weighting can be performed as a function of an indicator calculated on the server, for each receiver featuring in the database of the provider, on the basis of the configuration parameters of said receiver. It is for example possible to use the characteristics of the antennas, of the microchip, notably the code noise indicated by the provider of the receiver.
- the atmospheric error is normally the sum of the ionospheric and tropospheric errors. It is possible, however, to decide to neglect one of the two, the atmospheric error taken into account then being reduced to the other.
- the tropospheric error is customarily decomposed into two terms: the wet delay or SWD, for Slant Wet Delay which represents the contribution of water vapour and corresponds to about 10% of the total tropospheric delay and the hydrostatic delay (SHD, for Slant Hydrostatic Delay) which represents the remaining contribution and can be calculated on the basis of the equation of hydrostatic equilibrium.
- SWD wet delay
- SHD Slant Hydrostatic Delay
- the drawback of a sparse network is that the tropospheric error will be precise only in proximity to the stations. It is therefore important to have a dense network of stations in order to make a precise evaluation thereof. Organisations such as the IGN in France which operates this kind of dense network carry out a determination of the tropospheric error. Access to this error customarily requires a commercial agreement. [0081]
- the description of the invention given for the ionospheric error is adapted in the following manner to determine at least one of the contributions (wet delay, hydrostatic delay) of the tropospheric error.
- the step of determining the error based on the precise position of the receiver, comes the step of transmitting this error which is done without transmitting the precise position. Indeed, since the troposphere varies only very little on aDiac scale, a position of an order of magnitude of a kilometre is generally sufficient.
- the tropospheric component is determined by an external means such as that of the IGN network in France. The atmospheric error is then the sum of this tropospheric error with the ionospheric error determined according to the invention.
- the tropospheric error alone is determined according to the invention.
- the tropospheric and ionospheric error is determined, simultaneously or successively by the same receivers.
- a combination of the previous variants is implemented so as to exploit the environment and peculiarities and capability of the various receivers.
- PPP precise positioning
- Figure 4 represents a general flowchart of the processing actions for a use of the data of a server of positioning data, according to several embodiments of the invention.
- the data made available on the server 230 can be transmitted to the receivers of the users in multicast mode, or withdrawn on the server, either by an automatic procedure, or by a request procedure of said receivers addressed to said server.
- the distribution of the data of the server can be performed either by message or in http or https mode. A combination of the various modes of distribution is possible, according to the invention.
- Step 410 consists of this broadcasting in multicast mode or the request addressed to the server.
- the receiver acquires the satellite axis or axes in sight.
- step 430 the receiver position determined by the server according to a variant of the invention, on the basis of data of the mobile radiocommunications operator to which the user of the receiver is subscribed.
- This step makes it possible, combined with the ionospheric error data l(Z) applicable to the receiver, to shorten the delay in acquiring the satellite axes and the convergence delay of the following step of calculating the PVT on the basis of the satellite axes.
- the uncorrected PVT positioning data are calculated by the receiver.
- the PVT data are corrected by fusion with the data l(Z) transmitted by the server.
- the precision of the position measurement provided by the GNSS system is improved by fusion with the data relating to the ionospheric errors transmitted by the server 230.
- the time for acquisition and the time for convergence of the calculations are also improved.
Landscapes
- Engineering & Computer Science (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
Description
Claims
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1458336A FR3025610B1 (en) | 2014-09-05 | 2014-09-05 | METHOD FOR COLLABORATIVE DETERMINATION OF POSITIONING ERRORS OF A SATELLITE NAVIGATION SYSTEM |
PCT/EP2015/070047 WO2016034622A1 (en) | 2014-09-05 | 2015-09-02 | Method of collaborative determination of positioning errors of a satellite-based navigation system |
Publications (2)
Publication Number | Publication Date |
---|---|
EP3189349A1 true EP3189349A1 (en) | 2017-07-12 |
EP3189349B1 EP3189349B1 (en) | 2024-03-20 |
Family
ID=53008569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP15763847.9A Active EP3189349B1 (en) | 2014-09-05 | 2015-09-02 | Method of collaborative determination of positioning errors of a satellite-based navigation system |
Country Status (6)
Country | Link |
---|---|
US (1) | US10641902B2 (en) |
EP (1) | EP3189349B1 (en) |
KR (1) | KR20170039718A (en) |
CN (1) | CN107076856A (en) |
FR (1) | FR3025610B1 (en) |
WO (1) | WO2016034622A1 (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10012738B2 (en) * | 2016-03-24 | 2018-07-03 | Hitachi Zosen Corporation | Positioning method and positioning apparatus using satellite positioning system |
EP3226034A1 (en) | 2016-04-01 | 2017-10-04 | Centre National d'Etudes Spatiales | Improved gnss receiver using velocity integration |
EP3316001B1 (en) | 2016-10-25 | 2024-06-12 | Centre National d'Etudes Spatiales | Collaborative improvement of a vehicle's positioning |
US10304342B2 (en) * | 2016-11-08 | 2019-05-28 | Ge Aviation Systems Llc | Ground-based data acquisition system |
DE102017206275A1 (en) * | 2017-04-12 | 2018-10-18 | Robert Bosch Gmbh | A method of operating a correction service system and correction service system |
FR3085082B1 (en) * | 2018-08-17 | 2020-07-17 | Continental Automotive France | ESTIMATION OF THE GEOGRAPHICAL POSITION OF A ROAD VEHICLE FOR PARTICIPATORY PRODUCTION OF ROAD DATABASES |
IL269263B (en) | 2019-09-10 | 2021-09-30 | Veeride Geo Ltd | Cellular-based navigation method |
CN111045042B (en) * | 2019-12-20 | 2022-03-04 | 西安空间无线电技术研究所 | PPP-RTK enhancement method and system based on 'cloud-end' framework |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5867411A (en) * | 1996-12-19 | 1999-02-02 | The Aerospace Corporation | Kalman filter ionospheric delay estimator |
US6397147B1 (en) * | 2000-06-06 | 2002-05-28 | Csi Wireless Inc. | Relative GPS positioning using a single GPS receiver with internally generated differential correction terms |
TWI269052B (en) * | 2002-02-13 | 2006-12-21 | Matsushita Electric Works Ltd | Ionospheric error prediction and correction in satellite positioning systems |
US6701253B2 (en) * | 2002-02-19 | 2004-03-02 | Eride, Inc. | Total correction strategy |
RU2348052C2 (en) * | 2003-04-17 | 2009-02-27 | Секретэри Оф Стейт Фор Дефенс Дстл | Correction of errors caused by troposphere in global systems of location detection |
WO2006073490A2 (en) * | 2004-07-23 | 2006-07-13 | California Institute Of Technology | Generating high precision ionospheric ground-truth measurements |
US7656352B2 (en) * | 2005-09-20 | 2010-02-02 | Novariant, Inc. | Troposphere corrections for ground based positioning systems |
FR2929015B1 (en) * | 2008-03-21 | 2010-04-23 | Thales Sa | NETWORK AND METHOD FOR CALCULATING IONOSPHERIC CORRECTIONS |
US20120119944A1 (en) * | 2010-05-30 | 2012-05-17 | Trimble Navigation Limited | Gnss atmospheric estimation with ambiguity fixing |
US8803735B2 (en) * | 2010-11-19 | 2014-08-12 | Agjunction Llc | Portable base station network for local differential GNSS corrections |
DE102011051100A1 (en) * | 2011-06-16 | 2012-12-20 | Universität Kassel | Method for providing correction signals for terminal of satellite supported position determining system, involves determining position of mobile reference device for unit of position determining system in form of position data |
CN102426374A (en) * | 2011-09-02 | 2012-04-25 | 东南大学 | GPS (Global Positioning System) mobile base station rapid positioning and resolving method |
US20130127658A1 (en) * | 2011-11-22 | 2013-05-23 | Radio Systems Corporation | Method and Apparatus to Determine Actionable Position and Speed in GNSS Applications |
CN103018758A (en) * | 2012-12-03 | 2013-04-03 | 东南大学 | Method for moving differential base station based on global positioning system (GPS)/inertial navigation system (INS)/assisted global positioning system (AGPS) |
DE112013007301B4 (en) * | 2013-08-07 | 2023-08-31 | Topcon Positioning Systems, lnc. | Attenuation of scintillations in signals from global navigation satellite systems caused by ionospheric anomalies |
CN103592653B (en) * | 2013-10-17 | 2016-08-17 | 中国科学院光电研究院 | Ionosphere delay modification method for local area single-frequency Satellite navigation users |
CN103645482A (en) * | 2013-12-25 | 2014-03-19 | 上海华测导航技术有限公司 | Arbitrary point based set achievement method of GPS (Global Position System) real-time dynamic measurement of base station |
EP3106897A1 (en) * | 2015-06-19 | 2016-12-21 | Centre National d'Etudes Spatiales | Gnss receiver with an on-board capability to implement an optimal error correction mode |
EP3226034A1 (en) * | 2016-04-01 | 2017-10-04 | Centre National d'Etudes Spatiales | Improved gnss receiver using velocity integration |
WO2018034580A1 (en) * | 2016-08-15 | 2018-02-22 | Limited Liability Company "Topcon Positioning Systems" | Adaptive gnss positioning method capable of adjusting to ionospheric delay |
US10690775B2 (en) * | 2017-06-29 | 2020-06-23 | Novatel Inc. | Crowdsourcing atmospheric correction data |
DK201870263A1 (en) * | 2018-05-02 | 2019-12-03 | Intelligent Marking Aps | Method for marking a ground surface using a robot unit and a local base station, the system therefore and use thereof |
-
2014
- 2014-09-05 FR FR1458336A patent/FR3025610B1/en active Active
-
2015
- 2015-09-02 US US15/508,296 patent/US10641902B2/en active Active
- 2015-09-02 WO PCT/EP2015/070047 patent/WO2016034622A1/en active Application Filing
- 2015-09-02 CN CN201580047618.6A patent/CN107076856A/en active Pending
- 2015-09-02 KR KR1020177006063A patent/KR20170039718A/en active Search and Examination
- 2015-09-02 EP EP15763847.9A patent/EP3189349B1/en active Active
Non-Patent Citations (2)
Title |
---|
None * |
See also references of WO2016034622A1 * |
Also Published As
Publication number | Publication date |
---|---|
FR3025610B1 (en) | 2022-01-14 |
KR20170039718A (en) | 2017-04-11 |
US10641902B2 (en) | 2020-05-05 |
EP3189349B1 (en) | 2024-03-20 |
US20170307761A1 (en) | 2017-10-26 |
CN107076856A (en) | 2017-08-18 |
WO2016034622A1 (en) | 2016-03-10 |
FR3025610A1 (en) | 2016-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10641902B2 (en) | Method of collaborative determination of positioning errors of a satellite-based navigation system | |
US11709274B2 (en) | Determining correct location in the presence of GNSS spoofing | |
Orabi et al. | Opportunistic navigation with Doppler measurements from Iridium Next and Orbcomm LEO satellites | |
EP3115807B1 (en) | Gnss receiver with an on-board capability to implement an optimal error correction mode | |
US6243648B1 (en) | Fast acquisition position reporting system | |
US7800530B2 (en) | Method and system for providing assistance data for A-GPS location of handsets in wireless networks | |
CN111045034A (en) | GNSS multi-system real-time precise time transfer method and system based on broadcast ephemeris | |
CA2483318A1 (en) | Position determination for a wireless terminal in a hybrid position determination system | |
US11460583B2 (en) | Method and apparatus for providing correction data for satellite navigation | |
JP2014523519A (en) | A differential correction system reinforcement means for characterizing local errors using a target receiver capable of using non-GPS auxiliary PN & T signals | |
US20220291395A1 (en) | Polarization configurable gnss smartphone antenna | |
WO2014209967A2 (en) | Positioning using dtv broadcast signaling | |
KR20190070696A (en) | PPP-RTK Service Method for Port Navigation using State Space Representation in DGNSS Medium Frequency Wave | |
US6970785B2 (en) | Device for a mobile terminal for determining position by filtering integrity data from an augmentation device | |
Zhang et al. | Evaluation of qzss centimeter level augmentation system (clas): Open-sky to urban environments and geodetic to low-cost receivers | |
US11686850B2 (en) | Ionosphere grid history and compression for GNSS positioning | |
Chen et al. | Signal‐in‐Space and Positioning Performance of BDS Open Augmentation Service | |
Fakhreddine et al. | Data fusion for hybrid and autonomous time-of-flight positioning | |
EP4099061A1 (en) | Method for generating and providing a precise positioning solution of a mobile receiver in a gnss system by a central computation unit and a software product and its dissemination | |
Truong-Ngoc et al. | Robust GPS/Galileo/GLONASS Data Fusion Using Extended Kalman Filter | |
CN117368951A (en) | Carrier phase ambiguity rapid fixing system and method | |
CN116338751A (en) | Multi-satellite system and ground network lead-through fusion high-precision positioning module | |
CN117250649A (en) | Distance and Doppler combined quick positioning method based on single low-orbit satellite | |
CN116449403A (en) | Method and device for realizing precise single-point positioning and GNSS receiver | |
CN114076965A (en) | Method, apparatus and system for providing satellite positioning correction data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20170404 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1240652 Country of ref document: HK |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
17Q | First examination report despatched |
Effective date: 20200417 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
P01 | Opt-out of the competence of the unified patent court (upc) registered |
Effective date: 20230526 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
INTG | Intention to grant announced |
Effective date: 20231025 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE PATENT HAS BEEN GRANTED |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP Ref country code: HK Ref legal event code: WD Ref document number: 1240652 Country of ref document: HK |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 602015087989 Country of ref document: DE |