US20200309961A1 - Improved gnss receiver using velocity integration - Google Patents

Improved gnss receiver using velocity integration Download PDF

Info

Publication number
US20200309961A1
US20200309961A1 US16/089,735 US201716089735A US2020309961A1 US 20200309961 A1 US20200309961 A1 US 20200309961A1 US 201716089735 A US201716089735 A US 201716089735A US 2020309961 A1 US2020309961 A1 US 2020309961A1
Authority
US
United States
Prior art keywords
location
gnss receiver
error budget
error
receiver
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.)
Abandoned
Application number
US16/089,735
Inventor
François-Xavier MARMET
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Centre National dEtudes Spatiales CNES
Original Assignee
Centre National dEtudes Spatiales CNES
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Centre National dEtudes Spatiales CNES filed Critical Centre National dEtudes Spatiales CNES
Assigned to CENTRE NATIONAL D'ETUDES SPATIALES reassignment CENTRE NATIONAL D'ETUDES SPATIALES ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: MARMET, François-Xavier
Publication of US20200309961A1 publication Critical patent/US20200309961A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/421Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/426Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/22Multipath-related issues
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/32Multimode operation in a single same satellite system, e.g. GPS L1/L2
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/396Determining accuracy or reliability of position or pseudorange measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/46Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being of a radio-wave signal type
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/47Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/52Determining velocity

Definitions

  • the invention relates to a GNSS receiver. More specifically the GNSS receiver of the invention provides improved navigation accuracy and smoothness, notably in urban canyons, without a need to use complex multipath mitigation techniques and/or external aiding data.
  • navigation receivers becomes increasingly pervasive in everyday life. It is all the more common that cars' on-board electronics, smart phones, tablets include navigation receivers, and that applications running thereon capture as input information on position and trajectory of the user of the terminal.
  • Navigation receivers rely on L-Band RF signals transmitted by Medium Earth Orbiting satellites, which are generally included in constellations comprising tens of them to cover most of the surface of the earth, such as the GPSTM (US), GalileoTM (Europe), GlonassTM (Russia) and BeidouTM (China). These constellations are designated under the generic acronym of GNSS (Global Navigation Satellite System).
  • GNSS Global Navigation Satellite System
  • GNSS carrier signals are modulated by a pseudo-random code and a navigation message that allow calculation of a pseudo-range between the receiver and a specific satellite.
  • a pseudo-random code is modulated by a pseudo-random code and a navigation message that allow calculation of a pseudo-range between the receiver and a specific satellite.
  • PVT Position, Velocity and Time
  • PVT measurements are affected by a number of errors, some of which are intrinsic to the principle of measurement used (i.e, due to the deviation in trajectory of the RF signals through the atmosphere—ionosphere and troposphere—due to variations in the orbits of the satellites), intrinsic to the receiver and satellites imperfections (clock biases for instance), or intrinsic to some configurations of the satellites in view at a moment in time (i.e. elevation of the satellites over the horizon; low dispersion of visible satellites—high Dilution of Precision or DOP).
  • a number of corrections may be used to mitigate these errors, with the use of specific processing techniques which are only available to certain types of receivers.
  • bi-frequency receivers can mitigate ionospheric errors with a gain of precision from a few tens of meters to a few meters, and even better when combined with precise satellite orbits and clocks which then providing Precise Point Positioning (PPP)—a precision of a few tens of centimetres.
  • PPP Precise Point Positioning
  • Differential GPS and Real Time Kinematics solutions provide similar precision from integration of outside information (relative positioning vis-à-vis a number of fixed reference stations with known positions).
  • the degradation in UERE is due to the signal impairments of the specific satellite which is acquired or tracked by a tracking loop. Tracking of a satellite relies on a maximization of a correlation function between the acquired code signal and a number of local replicas generated by the receiver of the code signals which are specific to each satellite. The correlation functions will be corrupted by multipath and the satellite may not be correctly acquired or may be lost. Even if the signal tracking is still achievable, signal impairments will affect the correlation function's shape, thus degrading the pseudo-range estimation, and the UERE.
  • the present invention discloses a solution to overcome the previously cited drawbacks.
  • the invention discloses notably a GNSS receiver which is configured to use velocity measurements which are readily available whenever a condition is met whereby it is certain or likely that the velocity measurements will produce a position of the receiver with a better expected precision than the direct use of the position data.
  • this invention discloses a GNSS receiver comprising: first computing logic configured to calculate a first location and a second location of the GNSS receiver, the first location based on position data in a PVT vector, the second location based on a result of an integration of velocity data in the PVT vector; second computing logic configured to estimate a first error budget of the first location and a second error budget of the second location and to cause the first computing logic to output a calculated location of the GNSS receiver as a weighted average of the first location and the second location, wherein a first weight of the first location and a second weight of the second location in the calculated location are calculated by one or more of the first computing logic or the second computing logic, based on the first error budget and the second error budget.
  • a first weight of the first location in the calculated location equals a proportion of the first error budget in the total of the first error budget and the second error budget and a second weight of the second location in the calculated location equals a proportion of the second error budget in the total of the first error budget and the second error budget.
  • the first weight equals a hundred percent if the first error budget is lower than or equal to a combination of the second error budget and a first threshold value
  • the second weight equals a hundred percent if the second error budget is lower than a combination of the first error budget and a second threshold value.
  • the second location uses, based on a decision criterion, one or more of: i) a calculated location of the GNSS receiver at an epoch before a current epoch as an initial value of the integration of the velocity data in the PVT vector or ii) a location of the GNSS receiver at a previous epoch input from positioning aids as an initial value of the integration of the velocity data in the PVT vector.
  • the positioning aids comprise one or more of a map matching algorithm, an output of an inertial navigation system or a Wi-Fi localization system.
  • the estimated error budgets depend upon one or more of a configuration of said GNSS receiver, a geometry of GNSS satellites in view of said GNSS receiver or a topology in an environment of said GNSS receiver.
  • the estimated error budgets are estimated based on one or more of a Dilution Of Precision (DOP), a User Equivalent Range Error (UERE), a User Equivalent Range Rate Error (UERRE) or a probability estimate of a random walk of the integrated velocity along its path between a previous epoch and a current epoch.
  • DOP Dilution Of Precision
  • UERE User Equivalent Range Error
  • URRRE User Equivalent Range Rate Error
  • the first error budget equals the DOP multiplied by the UERE.
  • the second error budget equals a product of the DOP multiplied by the UERRE, times the duration during which velocity has been integrated.
  • the probability estimate of the random walk of the integrated velocity is modeled by an ellipsoid.
  • one or more of the UERE and the UERRE are a combination of one or more of a parameter value, a carrier to noise ratio dependent value, an elevation dependent value or a signal dependent value.
  • one or more of the first error budget or the second error budget are forecast estimated error budgets, said forecast being based on a combination of an estimated future trajectory of the GNSS receiver and a probability on said estimated future trajectory of one or more of a configuration of said GNSS receiver, a geometry of GNSS satellites in view of said GNSS receiver or a topology in an environment of said GNSS receiver.
  • one or more of the UERE and the UERRE are predefined and stored in a database accessible by the first computing logic.
  • the GNSS receiver of the invention further comprises a third computing logic configured to calculate a heading information from the PVT vector.
  • the invention also discloses a localization method using a GNSS receiver comprising: calculating one or more of a first location or a second location of the GNSS receiver, the first location based on position data in a PVT vector calculated by the GNSS receiver, the second location based on a result of an integration of velocity data in the PVT vector; estimating a first error budget of the first location and a second error budget of the second location and to cause the first computing logic to output a calculated location of the GNSS receiver as a weighted average of the first location and the second location; wherein a first weight of the first location and a second weight of the second location in the calculated location are calculated based on the first error budget and the second error budget.
  • the invention is easy to implement on any GNSS receiver. It is not even necessary to modify the receiver itself. Access to a number of data points which are normally readily accessible on standard receivers in well-known data formats (e.g. NMEA0183 or NMEA2000 protocols) will allow a GNSS receiver to practice the invention.
  • the invention brings position accuracy improvements, both in severe and benign environments. In doing so, it does not need external data, such as augmentation data, real time kinematics or information on the topology of the environment of the receiver.
  • FIG. 1 illustrates the multipath problem as viewed from a GNSS receiver of the prior art
  • FIGS. 2 a and 2 b illustrate typical errors impairing velocity and velocity-based position measurements for in a GNSS receiver of the invention in a number of its embodiments;
  • FIGS. 3 a and 3 b display examples of parameters used in computations of estimated error budgets of two position and velocity measurements of a standard GNSS receiver of the invention in a number of its embodiments;
  • FIG. 4 represents a display of a car navigation system of the prior art
  • FIG. 5 represents an architecture of a GNSS receiver modified to implement the invention in a number of its embodiments
  • FIG. 6 represents a flow chart of processing steps of computing logics of a receiver of the invention in a number of its embodiments
  • FIG. 7 illustrates a first example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments
  • FIG. 8 illustrates a second example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments.
  • FIG. 1 illustrates the multipath problem as viewed from a GNSS receiver of the prior art.
  • a receiver 110 is positioned between two buildings 120 . Certain signals 130 reach the receiver in a direct line. Other signals 140 reach the receiver after having been reflected from the buildings. A standard receiver is not able to distinguish between direct signals 130 and reflected signals 140 , thereby leading to a positioning error. Multipath errors at a given location of the receiver 110 will depend on the characteristics of the obstacles at the location (height, roughness of their surface, etc. . . . ) but also on the elevation of the satellites in view, and therefore on the time of the day, and on the meteorological conditions.
  • receivers which are configured or aided to cope easily with other types of errors, notably ionospheric errors, like bi-frequency receivers, are not efficient in a multipath environment, notably when a cold start is needed.
  • FIGS. 2 a and 2 b illustrate typical errors impairing velocity and velocity-based position measurements for in a GNSS receiver of the invention in a number of its embodiments.
  • FIG. 2 a displays the variations of GNSS velocity over time.
  • Lines 210 a illustrate the boundaries of the signal at 3 standard deviations (SD).
  • SD standard deviations
  • GNSS velocity being a direct measurement from the pseudo-ranges between the receiver and the satellites in view, it typically has a lower error budget than velocity output by mass-market inertial sensors, this velocity coming from acceleration integration.
  • hybridization solutions between a GNSS receiver and an inertial navigation system are complex to implement (especially as a retrofit solution), and for the time being are mostly confined to professional receivers, notably in the world of avionics.
  • FIG. 2 b displays a position calculated by integrating the velocity measurements of FIG. 2 a .
  • Line 210 b illustrates the same 3 SD boundary as lines 210 a. It can be seen that the position has the behavior of a random walk. Indeed, the comparison of the error budgets should be made after integration of the velocity to determine a position.
  • FIGS. 3 a and 3 b display examples of parameters used in computations of estimated error budgets of two position and velocity measurements of a standard GNSS receiver of the invention in a number of its embodiments
  • FIG. 3 a is taken from https://en.wikipedia.org/wiki/Error_analysis_for_the_Global_Positioning_System, to the author of which credit is given. It displays the factors which impact the position error budget, or UERE and the order of magnitude of these errors. It can be seen that the larger factor is the ionospheric effect with a possible error of 10 m ( ⁇ 5/+5 m), then the code tracking jitter ( ⁇ 3/+3 m), the satellite ephemerides, the satellite clocks, multipath distortion, tropospheric effect and carrier phase.
  • FIG. 3 b is taken from «Precise Velocity and Acceleration Determination Using a Standalone GPS Receiver in Real Time», J. ZHANG, PhD thesis (Table 6.2), to the author of which credit is given. It displays the factors which impact the velocity error budget, or User Equivalent Range Rate Error (UERRE) and the order of magnitude of these errors.
  • URRRE User Equivalent Range Rate Error
  • the error budget on velocity is two orders of magnitude lower than the UERE. This is essentially due to the fact that velocity is measured from the Doppler effect. Also, even for errors which are significant, their variation is much lower. For instance, the ionospheric effect drops from an error of up to 10 m to an error of an order of 1 cm/s. The multipath error becomes negligible.
  • velocity is integrated once to create a new calculated location which has, in general, a lower error budget.
  • FIG. 4 represents a display of a car navigation system of the prior art.
  • FIG. 4 is extracted from http://www.gpspower.net/garmin-tutorials/141052-secret-startup-commands-hidden-menus-garmin-devices.html, to the author of which credit is given.
  • NMEA National Marine Electronics Association
  • Two standards (NMEA 0183 and NMEA 2000) apply.
  • the receiver notably displays (and outputs at a port): Location, 410 , Velocity, 420 , C/N0, 430 , for each satellite in view which is designated by a number which is characteristic of its constellation, DOP, 450 , and Elevation 460 .
  • Other data may be available, like azimuth, raw pseudo-range measurements, estimated PVT residuals, among others.
  • FIG. 5 represents an architecture of a GNSS receiver modified to implement the invention in a number of its embodiments.
  • the receiver 510 is a core GNSS receiver in accordance with the present invention with one or more antennas, 511 , a number of radiofrequency processing channels, 512 and a digital processing capability, 513 .
  • the exact composition of the core receiver may vary depending on the use case (military or professional application, vehicle navigation, pedestrian dedicated device or smart phone). It may be capable of processing signals from multiple constellations (GPS, Glonass, Beidou, Galileo, etc. . . . ), and for each constellation on one, two, three or more frequencies.
  • the core receiver may be implemented in an FPGA or in an ASIC, or as a software-defined receiver, in multiple chips or in a single chipset.
  • the receiver 510 may output messages 531 , 532 compliant with the NMEA format, or compliant with other formats, including proprietary formats defined by manufacturers.
  • the manufacturers export their messages in their format or in NMEA format. Most of them publish the structure of the messages they output. As long as the output messages include the PVT calculated by the receiver, there will be enough information to practice the invention.
  • Output 531 and 532 have been split with contents mentioned hereunder, but this can be physically done through a single interface.
  • the processing logic 520 may be implemented in a hardware add-on module or in a software module loaded in a memory included in the receiver or added thereto.
  • the software may also be obtained as an application downloaded and/or installed on a smart phone. It may also be fully integrated with the core GNSS receiver, 510 , on a same board or in a same FPGA or a same ASIC.
  • Parameters are input, 521 , both from the receiver and possibly from tables stored in a memory contained in or adjacent to the processing logic 520 .
  • Parameters input from the core receiver through port 531 may include part or all of those listed in the description of FIG. 4 and other parameters which may be found useful to mix the first location and the second location, as further described below in relation with FIG. 6 .
  • the tables may include for instance, parameters to compute the error budgets, a first error budget for the GNSS calculated position (UERE) and a second error budget (UERRE) for the GNSS calculated velocity, or an error budget for a location calculated by integrating velocity.
  • These parameters may be preset values of UERE and UERRE. They may also comprise preset values of UERE and UERRE that depend on the satellites' elevations and/or their current C/NO and/or the GNSS signal type. They may also comprise a configuration of the receiver (type or model of the antenna, type or model of the RF processing stage, type or model of the digital processing stage). They may also comprise tables of multipath effects at locations on a trajectory of the receiver.
  • They may also comprise different values of thresholds which will be used to compute weights of a first location, calculated from P in PVT, and a second location, calculated from V in PVT, as provided by the invention.
  • a maximum value of UERE of 50 m may be set so that when the actual value is higher, the P measurements are completely discarded.
  • a maximum value of UERRE of 1 m/s For users requiring more precision, the thresholds will be lowered.
  • the first computing logic, 522 extracts P and V from the PVT data output through port or interface 532 from the core receiver.
  • P is a first location and a second location is computed by integrating V.
  • the distance calculated from the integration of velocity will be added to the location previously calculated by the first computing logic 522 , to compute the second location.
  • a second location may be input from an external source, 540 .
  • the external source may be an inertial module collocated with the GNSS receiver, for instance an arrangement comprising one or more of an accelerometer, a gyroscope and a magnetometer, as available in a smart phone.
  • the external source may also be the output of a map matching algorithm, as available in a car navigation system, or any other source of a position of a user or a device.
  • the second computing logic, 523 calculates the error budgets at the epoch from data which may combine predetermined information coming from the parameters input with information coming from the GNSS receiver, for instance a configuration of the receiver (type or model of the antenna, type or model of the RF processing stage, type or model of the digital processing stage), a geometry of GNSS satellites in view of said GNSS receiver, as provided by a NMEA message, or a topology in an environment of said GNSS receiver, depending on an evaluation of how the received signals are affected by multipath, based on a map of multipath effects available in the table.
  • the second computing logic may also use the C/N0 values of the channels processing the satellites in view, said values provided by NMEA messages.
  • the DOP may also use the DOP, also provided by a NMEA message. It may also use the probability of a random walk of the integrated velocity. In an embodiment, the random walk of the integrated velocity may be modeled by an ellipsoid.
  • the first error budget is computed as the DOP multiplied by the UERE.
  • the second error budget is calculated as a product of the DOP multiplied by the UERRE, times the duration during which velocity has been integrated.
  • One or more of the first computing logic, 522 , or of the second computing logic, 523 is configured to compute a first weight of the first error budget and a second weight of the second error budget, calculated for instance as a percentage of the total error budget.
  • the computing may be split between the first computing logic and the second computing logic to optimize the computation, without departing from the scope of the invention as claimed.
  • the weights may be calculated as a proportion of each error budget in the total of the first and the second error budgets.
  • the first weight equals a hundred percent if the first error budget is lower than or equal to a combination of the second error budget and a first threshold.
  • the first threshold may be static and drawn from the tables in the parameters values, 521 .
  • the first threshold may also be modified dynamically as a function of the conditions of reception at the GNSS receiver drawn from NMEA messages.
  • the second weight equals a hundred percent if the second error budget is lower than a combination of the first error budget and a second threshold.
  • the second threshold may be static and drawn from the tables in the parameters values, 521 .
  • the second threshold may also be modified dynamically as a function of the conditions of reception at the GNSS receiver drawn from NMEA messages.
  • the first computing logic is then caused to output a location of the GNSS receiver as a weighted average of the first location and of the second location.
  • the output location is the navigation solution.
  • only the first location or only the second location is taken into account to output the location of the GNSS receiver, based on relative values of the first and the second error budgets.
  • FIG. 6 represents a flow chart of processing steps of computing logics of a receiver of the invention in a number of its embodiments.
  • steps, 611 are performed in the first computing logic, 522 .
  • time reference t is obtained from the GNSS receiver at step 6111 .
  • time reference t will be an epoch after which a first P value with a confidence interval better than a threshold is obtained.
  • P is obtained from the PVT calculated by the digital processing stage of the GNSS receiver at step 6112 .
  • V is obtained from the PVT calculated by the digital processing stage of the GNSS receiver at step 6113 .
  • V is integrated once at step 6116 .
  • the improved position may thus be computed as:
  • Pinit being P or P′ at the starting time of improved position computation, or an averaged position obtained from several P or P′ values over time.
  • a position P′ is obtained at step 6115 from an external source, as explained above in relation to FIG. 5 .
  • a logic may be included in the first computing logic, to implement a decision criterion to select P or P′ at time t.
  • a decision criterion may take into account an accuracy of the GNSS position transmitted by an NMEA message.
  • a first location and a second location are computed.
  • the location which will be output as the navigation solution will be computed taking into account the first location, the second location and the weighting calculated either in the first computing logic (embodiment not represented on the figure) or in the second computing logic (as represented on the figure).
  • steps 612 are performed in the second computing logic, 523 .
  • parameters values are obtained from the parameters input 521 , either through the interface to the GNSS core receiver, 531 , or through an external input, 540 .
  • a first error budget and a second error budget are computed, as explained above in relation with FIG. 5 .
  • the weights of the error budgets are computed to be passed to the first computing logic.
  • the error budgets may be transferred to the first computing logic for it to compute the weights.
  • FIG. 7 illustrates a first example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments.
  • FIG. 7 The case illustrated by FIG. 7 is a trajectory of an automotive vehicle in a suburb of a city where the density of buildings is rather low.
  • Four constellations are acquired and tracked by the receiver (GPS, Glonass, Galileo, Beidou). The duration of the trajectory is about 50 minutes.
  • Simplistic assumptions are made on the error budgets: the UERE is set at a fixed value of 1 m and the UERRE is set at a fixed value 0.1 m/s. These assumptions are optimistic for the UERE and pessimistic for the UERRE given the local environment. In this case, no dynamic modification of the error budgets is added. In this example, a plain selection of P or integrated V is made based on the lower error budget:
  • the trajectory of the vehicle calculated by a known unaided GNSS receiver is represented by line 710 .
  • the trajectory of the vehicle calculated by a GNSS receiver to which an addition of the arrangement 520 configured according to the invention is represented by line 720 .
  • line 710 demonstrates a number of errors 730 , 740 , 750 in the trajectory, which does not follow the road.
  • This use case illustrates some of the advantages brought to the solutions of the prior art by a simplistic implementation of the invention, in some of its embodiments.
  • FIG. 8 illustrates a second example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments.
  • the case illustrated by FIG. 8 is a trajectory of an automotive vehicle in a centre of a city where the density of buildings is rather high.
  • the receiver used and the GNSS constellations acquired and tracked are the same as for the case of FIG. 7 .
  • the duration of the trajectory in this case is about 1 h 45 minutes.
  • Simple assumptions are made on the error budgets: the UERE is set at a fixed value of 2 m and the UERRE is set at a fixed value 0.2 m/s. These assumptions are optimistic for the UERE and pessimistic for the UERRE.
  • the other conditions of the case of FIG. 8 are the same as those of FIG. 7 .
  • the trajectory of the vehicle calculated by an unaided GNSS receiver of the prior art is represented by line 810 .
  • the trajectory of the vehicle calculated by a GNSS to which an addition of the arrangement 520 configured according to the invention is represented by line 820 .
  • line 810 demonstrates a number of errors 830 , 840 , 850 in the trajectory, which does not follow the road.
  • This use case further illustrates some of the advantages of the invention over the solutions of the prior art in a constrained environment by an implementation of the invention, in some of its embodiments.
  • the invention provides not only a more accurate navigation solution, but also a smoother trajectory, since it is not affected by errors 830 , 840 and 850 .
  • Standard receivers are, for example, the ones marked as «recreational», «handheld», «navigation» «vehicle tracking» or «board for OEM apps» in the Receiver Survey 2015 database by GPSWorld, available at http://gpsworld.com/wp-content/uploads/2015/01/GPSWorld_2015ReceiverSurvey.pdf)
  • the invention was described above in a context of receivers using GNSS signals. It may also be applied with receivers using signals transmitted by pseudolites, for examples those marketed by LocataTM.
  • vector V which is used to calculate the second location may also be used to determine a navigation heading which can be used either alone or in combination with outputs from other heading sensors, such as a magnetometer, or a combination of accelerometers and gyroscopes, or a map matching algorithm.
  • the navigation heading may be used by applications, for instance to give a direction of a Point of Interest relative to the navigation heading. It can also be used to calibrate other heading sensors. This can be applied to a receiver used by a pedestrian or located in a manned or unmanned car, aircraft, or boat.

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)
  • Navigation (AREA)

Abstract

The invention discloses an improved GNSS receiver which determines a location of the receiver by combining a first location determined either from the standard PVT and/or from a positioning aid like a map matching algorithm, inertial navigation system, WiFi localization system or other, and a second location determined by integrating the velocity from the standard PVT. The combination is based on a weighting of the error budgets of the first position and the second location. The improved receiver is preferably based on a standard receiver with an add-on software module which receives and processes data transmitted from the standard receiver by, for example, NMEA messages. The improved receiver allows a determination of a more precise and smoother trajectory in a simple way.

Description

    FIELD OF THE INVENTION
  • The invention relates to a GNSS receiver. More specifically the GNSS receiver of the invention provides improved navigation accuracy and smoothness, notably in urban canyons, without a need to use complex multipath mitigation techniques and/or external aiding data.
  • BACKGROUND
  • The use of navigation receivers becomes increasingly pervasive in everyday life. It is all the more common that cars' on-board electronics, smart phones, tablets include navigation receivers, and that applications running thereon capture as input information on position and trajectory of the user of the terminal.
  • Navigation receivers rely on L-Band RF signals transmitted by Medium Earth Orbiting satellites, which are generally included in constellations comprising tens of them to cover most of the surface of the earth, such as the GPS™ (US), Galileo™ (Europe), Glonass™ (Russia) and Beidou™ (China). These constellations are designated under the generic acronym of GNSS (Global Navigation Satellite System).
  • GNSS carrier signals are modulated by a pseudo-random code and a navigation message that allow calculation of a pseudo-range between the receiver and a specific satellite. With a minimum of four pseudo-ranges, it is possible to calculate Position, Velocity and Time (PVT) of the receiver. In receivers of the type used by consumers, the position information is the one which is directly used to compute the navigation solution.
  • PVT measurements are affected by a number of errors, some of which are intrinsic to the principle of measurement used (i.e, due to the deviation in trajectory of the RF signals through the atmosphere—ionosphere and troposphere—due to variations in the orbits of the satellites), intrinsic to the receiver and satellites imperfections (clock biases for instance), or intrinsic to some configurations of the satellites in view at a moment in time (i.e. elevation of the satellites over the horizon; low dispersion of visible satellites—high Dilution of Precision or DOP). A number of corrections may be used to mitigate these errors, with the use of specific processing techniques which are only available to certain types of receivers. For instance, bi-frequency receivers can mitigate ionospheric errors with a gain of precision from a few tens of meters to a few meters, and even better when combined with precise satellite orbits and clocks which then providing Precise Point Positioning (PPP)—a precision of a few tens of centimetres. Differential GPS and Real Time Kinematics solutions provide similar precision from integration of outside information (relative positioning vis-à-vis a number of fixed reference stations with known positions).
  • It is more difficult to mitigate in a consistent and efficient manner some errors which depend on the position of the receiver, notably when this position is surrounded by a number of objects which reflect the navigation RF signals and/or mask a number of the satellites which should be in line of sight (LOS) at a moment in time. In such conditions, the precision of the calculation of the PVT may be quite poor, all other causes of errors being equal, both at the time of acquiring a GNSS signal and at the time of tracking said signal.
  • In urban canyons (i.e. streets in between tall buildings), multipath will not only increase the error in the determination of the pseudo-range of a satellite (User Equivalent Range Error or UERE), but also the (Geometric) Dilution Of Precision, (GDOP or DOP), because the field of view of the antenna will be narrower thus limiting the increase in precision due to the use of additional satellites.
  • The degradation in UERE is due to the signal impairments of the specific satellite which is acquired or tracked by a tracking loop. Tracking of a satellite relies on a maximization of a correlation function between the acquired code signal and a number of local replicas generated by the receiver of the code signals which are specific to each satellite. The correlation functions will be corrupted by multipath and the satellite may not be correctly acquired or may be lost. Even if the signal tracking is still achievable, signal impairments will affect the correlation function's shape, thus degrading the pseudo-range estimation, and the UERE.
  • Most of the mitigation techniques therefore rely on specific processing applied at the correlation stage. For a general overview of these prior art techniques which can be applied, see for instance Bhuiyan et al, Advanced Multipath Mitigation Techniques for Satellite-Based Positioning Applications, International Journal of Navigation and Observation, Volume 2010, Article ID 412393. Among these techniques are: using narrow correlators (i.e. correlators spaced by much less than a chip—or code length); using numerous correlators per tracking channel instead of a standard number of three and performing a double delta calculation; Multipath Estimation Delay Lock Loop, which uses the output of the correlation functions to estimate a full navigation signal of a definite satellite comprising LOS and non LOS signals.
  • All these solutions are complex and costly to implement either at the hardware level or at the software level. In any case, they cannot be currently implemented in a standard consumer receiver of the type used in a smart phone or a car navigation system.
  • It is therefore important to find a solution to mitigate the effects of multipath or other kinds of impairments to a navigation solution (e.g. strong local ionospheric disturbances impairing signal from one or several satellites) which could be implemented in a standard receiver available to the general public.
  • The present invention discloses a solution to overcome the previously cited drawbacks.
  • SUMMARY OF THE INVENTION
  • To this effect, the invention discloses notably a GNSS receiver which is configured to use velocity measurements which are readily available whenever a condition is met whereby it is certain or likely that the velocity measurements will produce a position of the receiver with a better expected precision than the direct use of the position data.
  • To this effect, this invention discloses a GNSS receiver comprising: first computing logic configured to calculate a first location and a second location of the GNSS receiver, the first location based on position data in a PVT vector, the second location based on a result of an integration of velocity data in the PVT vector; second computing logic configured to estimate a first error budget of the first location and a second error budget of the second location and to cause the first computing logic to output a calculated location of the GNSS receiver as a weighted average of the first location and the second location, wherein a first weight of the first location and a second weight of the second location in the calculated location are calculated by one or more of the first computing logic or the second computing logic, based on the first error budget and the second error budget.
  • Advantageously, a first weight of the first location in the calculated location equals a proportion of the first error budget in the total of the first error budget and the second error budget and a second weight of the second location in the calculated location equals a proportion of the second error budget in the total of the first error budget and the second error budget.
  • Advantageously, the first weight equals a hundred percent if the first error budget is lower than or equal to a combination of the second error budget and a first threshold value,
  • Advantageously, the second weight equals a hundred percent if the second error budget is lower than a combination of the first error budget and a second threshold value.
  • Advantageously, the second location uses, based on a decision criterion, one or more of: i) a calculated location of the GNSS receiver at an epoch before a current epoch as an initial value of the integration of the velocity data in the PVT vector or ii) a location of the GNSS receiver at a previous epoch input from positioning aids as an initial value of the integration of the velocity data in the PVT vector.
  • Advantageously, the positioning aids comprise one or more of a map matching algorithm, an output of an inertial navigation system or a Wi-Fi localization system.
  • Advantageously, the estimated error budgets depend upon one or more of a configuration of said GNSS receiver, a geometry of GNSS satellites in view of said GNSS receiver or a topology in an environment of said GNSS receiver.
  • Advantageously, the estimated error budgets are estimated based on one or more of a Dilution Of Precision (DOP), a User Equivalent Range Error (UERE), a User Equivalent Range Rate Error (UERRE) or a probability estimate of a random walk of the integrated velocity along its path between a previous epoch and a current epoch.
  • Advantageously, the first error budget equals the DOP multiplied by the UERE.
  • Advantageously, the second error budget equals a product of the DOP multiplied by the UERRE, times the duration during which velocity has been integrated.
  • Advantageously, the probability estimate of the random walk of the integrated velocity is modeled by an ellipsoid.
  • Advantageously, one or more of the UERE and the UERRE are a combination of one or more of a parameter value, a carrier to noise ratio dependent value, an elevation dependent value or a signal dependent value.
  • Advantageously, one or more of the first error budget or the second error budget are forecast estimated error budgets, said forecast being based on a combination of an estimated future trajectory of the GNSS receiver and a probability on said estimated future trajectory of one or more of a configuration of said GNSS receiver, a geometry of GNSS satellites in view of said GNSS receiver or a topology in an environment of said GNSS receiver.
  • Advantageously, one or more of the UERE and the UERRE are predefined and stored in a database accessible by the first computing logic.
  • Advantageously, the GNSS receiver of the invention further comprises a third computing logic configured to calculate a heading information from the PVT vector.
  • The invention also discloses a localization method using a GNSS receiver comprising: calculating one or more of a first location or a second location of the GNSS receiver, the first location based on position data in a PVT vector calculated by the GNSS receiver, the second location based on a result of an integration of velocity data in the PVT vector; estimating a first error budget of the first location and a second error budget of the second location and to cause the first computing logic to output a calculated location of the GNSS receiver as a weighted average of the first location and the second location; wherein a first weight of the first location and a second weight of the second location in the calculated location are calculated based on the first error budget and the second error budget.
  • The invention is easy to implement on any GNSS receiver. It is not even necessary to modify the receiver itself. Access to a number of data points which are normally readily accessible on standard receivers in well-known data formats (e.g. NMEA0183 or NMEA2000 protocols) will allow a GNSS receiver to practice the invention. The invention brings position accuracy improvements, both in severe and benign environments. In doing so, it does not need external data, such as augmentation data, real time kinematics or information on the topology of the environment of the receiver.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • The invention and its advantages will be better understood upon reading the following detailed description of a particular embodiment, given purely by way of non-limiting example, this description being made with reference to the accompanying drawings in which:
  • FIG. 1 illustrates the multipath problem as viewed from a GNSS receiver of the prior art;
  • FIGS. 2a and 2b illustrate typical errors impairing velocity and velocity-based position measurements for in a GNSS receiver of the invention in a number of its embodiments;
  • FIGS. 3a and 3b display examples of parameters used in computations of estimated error budgets of two position and velocity measurements of a standard GNSS receiver of the invention in a number of its embodiments;
  • FIG. 4 represents a display of a car navigation system of the prior art;
  • FIG. 5 represents an architecture of a GNSS receiver modified to implement the invention in a number of its embodiments;
  • FIG. 6 represents a flow chart of processing steps of computing logics of a receiver of the invention in a number of its embodiments;
  • FIG. 7 illustrates a first example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments;
  • FIG. 8 illustrates a second example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments.
  • DETAILED DESCRIPTION
  • FIG. 1 illustrates the multipath problem as viewed from a GNSS receiver of the prior art.
  • A receiver 110 is positioned between two buildings 120. Certain signals 130 reach the receiver in a direct line. Other signals 140 reach the receiver after having been reflected from the buildings. A standard receiver is not able to distinguish between direct signals 130 and reflected signals 140, thereby leading to a positioning error. Multipath errors at a given location of the receiver 110 will depend on the characteristics of the obstacles at the location (height, roughness of their surface, etc. . . . ) but also on the elevation of the satellites in view, and therefore on the time of the day, and on the meteorological conditions.
  • It is quite difficult to correct multipath errors, notably because they are location and time dependent and therefore require significant memory and/or processing power. Also, receivers which are configured or aided to cope easily with other types of errors, notably ionospheric errors, like bi-frequency receivers, are not efficient in a multipath environment, notably when a cold start is needed.
  • It has thus always been an issue to find a solution to this problem, as with other types of impairments which impact the error budget of a navigation solution. And even more so in a simple way which can be applied with standard GNSS receivers.
  • FIGS. 2a and 2b illustrate typical errors impairing velocity and velocity-based position measurements for in a GNSS receiver of the invention in a number of its embodiments.
  • FIG. 2a displays the variations of GNSS velocity over time. Lines 210 a illustrate the boundaries of the signal at 3 standard deviations (SD). GNSS velocity being a direct measurement from the pseudo-ranges between the receiver and the satellites in view, it typically has a lower error budget than velocity output by mass-market inertial sensors, this velocity coming from acceleration integration. Also, hybridization solutions between a GNSS receiver and an inertial navigation system are complex to implement (especially as a retrofit solution), and for the time being are mostly confined to professional receivers, notably in the world of avionics.
  • FIG. 2b displays a position calculated by integrating the velocity measurements of FIG. 2a . Line 210 b illustrates the same 3 SD boundary as lines 210 a. It can be seen that the position has the behavior of a random walk. Indeed, the comparison of the error budgets should be made after integration of the velocity to determine a position.
  • FIGS. 3a and 3b display examples of parameters used in computations of estimated error budgets of two position and velocity measurements of a standard GNSS receiver of the invention in a number of its embodiments
  • FIG. 3a is taken from https://en.wikipedia.org/wiki/Error_analysis_for_the_Global_Positioning_System, to the author of which credit is given. It displays the factors which impact the position error budget, or UERE and the order of magnitude of these errors. It can be seen that the larger factor is the ionospheric effect with a possible error of 10 m (−5/+5 m), then the code tracking jitter (−3/+3 m), the satellite ephemerides, the satellite clocks, multipath distortion, tropospheric effect and carrier phase.
  • FIG. 3b is taken from «Precise Velocity and Acceleration Determination Using a Standalone GPS Receiver in Real Time», J. ZHANG, PhD thesis (Table 6.2), to the author of which credit is given. It displays the factors which impact the velocity error budget, or User Equivalent Range Rate Error (UERRE) and the order of magnitude of these errors. We can see that, in general, the error budget on velocity is two orders of magnitude lower than the UERE. This is essentially due to the fact that velocity is measured from the Doppler effect. Also, even for errors which are significant, their variation is much lower. For instance, the ionospheric effect drops from an error of up to 10 m to an error of an order of 1 cm/s. The multipath error becomes negligible.
  • According to the invention, then velocity is integrated once to create a new calculated location which has, in general, a lower error budget.
  • There are conditions under which the error on integrated velocity may become larger than the error on position. This may be the case for instance, when the product of the UERE×DOP becomes lower than the standard deviation of the integrated velocity, notably when the integration time is long. It is therefore quite advantageous to be able to then switch from velocity measurements to position measurements. It is also advantageous to use a precise position as an initial position for the integrated position.
  • FIG. 4 represents a display of a car navigation system of the prior art.
  • FIG. 4 is extracted from http://www.gpspower.net/garmin-tutorials/141052-secret-startup-commands-hidden-menus-garmin-devices.html, to the author of which credit is given.
  • The National Marine Electronics Association (NMEA) has defined standard sentences of a format for outputting from a GNSS receiver a large amount of information which can be used to implement the invention using a standard receiver. Two standards (NMEA 0183 and NMEA 2000) apply.
  • The receiver notably displays (and outputs at a port): Location, 410, Velocity, 420, C/N0, 430, for each satellite in view which is designated by a number which is characteristic of its constellation, DOP, 450, and Elevation 460. Other data may be available, like azimuth, raw pseudo-range measurements, estimated PVT residuals, among others.
  • FIG. 5 represents an architecture of a GNSS receiver modified to implement the invention in a number of its embodiments.
  • The receiver 510 is a core GNSS receiver in accordance with the present invention with one or more antennas, 511, a number of radiofrequency processing channels, 512 and a digital processing capability, 513. The exact composition of the core receiver may vary depending on the use case (military or professional application, vehicle navigation, pedestrian dedicated device or smart phone). It may be capable of processing signals from multiple constellations (GPS, Glonass, Beidou, Galileo, etc. . . . ), and for each constellation on one, two, three or more frequencies. The core receiver may be implemented in an FPGA or in an ASIC, or as a software-defined receiver, in multiple chips or in a single chipset.
  • As indicated above, the receiver 510 may output messages 531, 532 compliant with the NMEA format, or compliant with other formats, including proprietary formats defined by manufacturers. Generally, the manufacturers export their messages in their format or in NMEA format. Most of them publish the structure of the messages they output. As long as the output messages include the PVT calculated by the receiver, there will be enough information to practice the invention. Output 531 and 532 have been split with contents mentioned hereunder, but this can be physically done through a single interface.
  • These messages are fed to computing logic, 520 which will compute the navigation solution according to the invention, i.e. from a weighted average of a first location being the Position of the calculated PVT, and a second location which is the Velocity of the calculated PVT. The processing logic 520 may be implemented in a hardware add-on module or in a software module loaded in a memory included in the receiver or added thereto. The software may also be obtained as an application downloaded and/or installed on a smart phone. It may also be fully integrated with the core GNSS receiver, 510, on a same board or in a same FPGA or a same ASIC.
  • Parameters are input, 521, both from the receiver and possibly from tables stored in a memory contained in or adjacent to the processing logic 520.
  • Parameters input from the core receiver through port 531 may include part or all of those listed in the description of FIG. 4 and other parameters which may be found useful to mix the first location and the second location, as further described below in relation with FIG. 6.
  • The tables may include for instance, parameters to compute the error budgets, a first error budget for the GNSS calculated position (UERE) and a second error budget (UERRE) for the GNSS calculated velocity, or an error budget for a location calculated by integrating velocity. These parameters may be preset values of UERE and UERRE. They may also comprise preset values of UERE and UERRE that depend on the satellites' elevations and/or their current C/NO and/or the GNSS signal type. They may also comprise a configuration of the receiver (type or model of the antenna, type or model of the RF processing stage, type or model of the digital processing stage). They may also comprise tables of multipath effects at locations on a trajectory of the receiver. They may also comprise different values of thresholds which will be used to compute weights of a first location, calculated from P in PVT, and a second location, calculated from V in PVT, as provided by the invention. By way of example, for a pedestrian just looking for a coarse positioning, a maximum value of UERE of 50 m may be set so that when the actual value is higher, the P measurements are completely discarded. Likewise, a maximum value of UERRE of 1 m/s. For users requiring more precision, the thresholds will be lowered.
  • The first computing logic, 522, extracts P and V from the PVT data output through port or interface 532 from the core receiver. P is a first location and a second location is computed by integrating V. In an embodiment, the distance calculated from the integration of velocity will be added to the location previously calculated by the first computing logic 522, to compute the second location. But some variants are possible, for instance, based on a decision criterion, e.g. a threshold of the error budget or a determination that the receiver is deemed to be located in an area with a clear sky or not, a second location may be input from an external source, 540. The external source may be an inertial module collocated with the GNSS receiver, for instance an arrangement comprising one or more of an accelerometer, a gyroscope and a magnetometer, as available in a smart phone. The external source may also be the output of a map matching algorithm, as available in a car navigation system, or any other source of a position of a user or a device.
  • The second computing logic, 523, calculates the error budgets at the epoch from data which may combine predetermined information coming from the parameters input with information coming from the GNSS receiver, for instance a configuration of the receiver (type or model of the antenna, type or model of the RF processing stage, type or model of the digital processing stage), a geometry of GNSS satellites in view of said GNSS receiver, as provided by a NMEA message, or a topology in an environment of said GNSS receiver, depending on an evaluation of how the received signals are affected by multipath, based on a map of multipath effects available in the table. The second computing logic may also use the C/N0 values of the channels processing the satellites in view, said values provided by NMEA messages. It may also use the DOP, also provided by a NMEA message. It may also use the probability of a random walk of the integrated velocity. In an embodiment, the random walk of the integrated velocity may be modeled by an ellipsoid. In an embodiment, the first error budget is computed as the DOP multiplied by the UERE. In an embodiment, the second error budget is calculated as a product of the DOP multiplied by the UERRE, times the duration during which velocity has been integrated.
  • One or more of the first computing logic, 522, or of the second computing logic, 523, is configured to compute a first weight of the first error budget and a second weight of the second error budget, calculated for instance as a percentage of the total error budget. The computing may be split between the first computing logic and the second computing logic to optimize the computation, without departing from the scope of the invention as claimed. The weights may be calculated as a proportion of each error budget in the total of the first and the second error budgets.
  • In an embodiment, the first weight equals a hundred percent if the first error budget is lower than or equal to a combination of the second error budget and a first threshold. The first threshold may be static and drawn from the tables in the parameters values, 521. The first threshold may also be modified dynamically as a function of the conditions of reception at the GNSS receiver drawn from NMEA messages.
  • In an embodiment, the second weight equals a hundred percent if the second error budget is lower than a combination of the first error budget and a second threshold. The second threshold may be static and drawn from the tables in the parameters values, 521. The second threshold may also be modified dynamically as a function of the conditions of reception at the GNSS receiver drawn from NMEA messages.
  • According to the invention, the first computing logic is then caused to output a location of the GNSS receiver as a weighted average of the first location and of the second location. The output location is the navigation solution. In some embodiments, only the first location or only the second location is taken into account to output the location of the GNSS receiver, based on relative values of the first and the second error budgets.
  • It is to be noted that good results may be obtained without external data, even if this data may improve the navigation solution. This is due to the fact that parameters which are readily available from the GNSS receiver through the interface 531, 532, such as DOP, C/N0, elevation of the satellites, may be used to adapt preset UERE and UERRE to compute first and second error budgets adapted to the actual conditions or reception of the GNSS signals.
  • FIG. 6 represents a flow chart of processing steps of computing logics of a receiver of the invention in a number of its embodiments.
  • According to the invention, steps, 611, are performed in the first computing logic, 522.
  • A time reference t is obtained from the GNSS receiver at step 6111. On a cold start, time reference t will be an epoch after which a first P value with a confidence interval better than a threshold is obtained.
  • P is obtained from the PVT calculated by the digital processing stage of the GNSS receiver at step 6112.
  • At the same time reference, V is obtained from the PVT calculated by the digital processing stage of the GNSS receiver at step 6113. V is integrated once at step 6116.
  • The improved position may thus be computed as:
  • P improved ( t + 1 ) = P init + t t + 1 V ( τ ) × d τ
  • Pinit being P or P′ at the starting time of improved position computation, or an averaged position obtained from several P or P′ values over time.
  • In a variant of the invention, a position P′ is obtained at step 6115 from an external source, as explained above in relation to FIG. 5. In this variant, a logic may be included in the first computing logic, to implement a decision criterion to select P or P′ at time t. A decision criterion may take into account an accuracy of the GNSS position transmitted by an NMEA message.
  • Then, at step 6117 and 6118 respectively, a first location and a second location are computed.
  • At step 6119, the location which will be output as the navigation solution will be computed taking into account the first location, the second location and the weighting calculated either in the first computing logic (embodiment not represented on the figure) or in the second computing logic (as represented on the figure).
  • According to the invention, steps 612, are performed in the second computing logic, 523.
  • At step 6121, parameters values are obtained from the parameters input 521, either through the interface to the GNSS core receiver, 531, or through an external input, 540.
  • At steps 6122 and 6123, respectively, a first error budget and a second error budget are computed, as explained above in relation with FIG. 5.
  • At step 6124, the weights of the error budgets are computed to be passed to the first computing logic. Alternatively, the error budgets may be transferred to the first computing logic for it to compute the weights.
  • FIG. 7 illustrates a first example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments.
  • The case illustrated by FIG. 7 is a trajectory of an automotive vehicle in a suburb of a city where the density of buildings is rather low. Four constellations are acquired and tracked by the receiver (GPS, Glonass, Galileo, Beidou). The duration of the trajectory is about 50 minutes. Simplistic assumptions are made on the error budgets: the UERE is set at a fixed value of 1 m and the UERRE is set at a fixed value 0.1 m/s. These assumptions are optimistic for the UERE and pessimistic for the UERRE given the local environment. In this case, no dynamic modification of the error budgets is added. In this example, a plain selection of P or integrated V is made based on the lower error budget:
      • the first location computed from P at step 6117 is selected as the navigation solution if the first error budget is lower than the second error budget;
      • the second location computed at step 6118 (as the previous first location to which a distance calculated by integration of V is added) is selected, if the second error budget is lower than the first error budget.
  • In FIG. 7, the trajectory of the vehicle calculated by a known unaided GNSS receiver is represented by line 710. The trajectory of the vehicle calculated by a GNSS receiver to which an addition of the arrangement 520 configured according to the invention is represented by line 720.
  • It can be seen that line 710 demonstrates a number of errors 730, 740, 750 in the trajectory, which does not follow the road. This use case illustrates some of the advantages brought to the solutions of the prior art by a simplistic implementation of the invention, in some of its embodiments.
  • FIG. 8 illustrates a second example of a comparison of calculations of positions of a GNSS receiver of the prior art and of a GNSS receiver of the invention in a number of its embodiments.
  • The case illustrated by FIG. 8 is a trajectory of an automotive vehicle in a centre of a city where the density of buildings is rather high. The receiver used and the GNSS constellations acquired and tracked are the same as for the case of FIG. 7. The duration of the trajectory in this case is about 1 h 45 minutes. Simple assumptions are made on the error budgets: the UERE is set at a fixed value of 2 m and the UERRE is set at a fixed value 0.2 m/s. These assumptions are optimistic for the UERE and pessimistic for the UERRE. The other conditions of the case of FIG. 8 are the same as those of FIG. 7.
  • In FIG. 8, the trajectory of the vehicle calculated by an unaided GNSS receiver of the prior art is represented by line 810. The trajectory of the vehicle calculated by a GNSS to which an addition of the arrangement 520 configured according to the invention is represented by line 820.
  • It can be seen that line 810 demonstrates a number of errors 830, 840, 850 in the trajectory, which does not follow the road. This use case further illustrates some of the advantages of the invention over the solutions of the prior art in a constrained environment by an implementation of the invention, in some of its embodiments.
  • It is remarkable that the invention provides not only a more accurate navigation solution, but also a smoother trajectory, since it is not affected by errors 830, 840 and 850.
  • The improvements described and claimed in some embodiments of the invention will still improve the navigation solution in comparison to the cases illustrated above. For instance, the integration of dynamic parameters representative of the conditions of reception at the GNSS receiver, or the weighting of the two location calculations based on their error budgets, or the integration of external navigation aids. It is remarkable that these improvements may be obtained with a standard GNSS receiver to which the invention may be added in a simple way.
  • Standard receivers are, for example, the ones marked as «recreational», «handheld», «navigation» «vehicle tracking» or «board for OEM apps» in the Receiver Survey 2015 database by GPSWorld, available at http://gpsworld.com/wp-content/uploads/2015/01/GPSWorld_2015ReceiverSurvey.pdf)
  • The invention was described above in a context of receivers using GNSS signals. It may also be applied with receivers using signals transmitted by pseudolites, for examples those marketed by Locata™.
  • According to the invention, vector V which is used to calculate the second location may also be used to determine a navigation heading which can be used either alone or in combination with outputs from other heading sensors, such as a magnetometer, or a combination of accelerometers and gyroscopes, or a map matching algorithm. The navigation heading may be used by applications, for instance to give a direction of a Point of Interest relative to the navigation heading. It can also be used to calibrate other heading sensors. This can be applied to a receiver used by a pedestrian or located in a manned or unmanned car, aircraft, or boat.
  • The examples disclosed in this specification are only illustrative of some embodiments of the invention. They do not in any manner limit the scope of said invention which is defined by the appended claims.

Claims (16)

1. A GNSS receiver comprising:
first computing logic configured to calculate a first location and a second location of the GNSS receiver, the first location based on position data in a PVT vector, the second location based on a result of an integration of velocity data in the PVT vector;
second computing logic configured to estimate a first error budget of the first location and a second error budget of the second location and to cause the first computing logic to output a calculated location of the GNSS receiver as a weighted average of the first location and the second location, wherein a first weight of the first location and a second weight of the second location in the calculated location are calculated by one or more of the first computing logic or the second computing logic, based on the first error budget and the second error budget.
2. The GNSS receiver of claim 1, wherein a first weight of the first location in the calculated location equals a proportion of the first error budget in the total of the first error budget and the second error budget and a second weight of the second location in the calculated location equals a proportion of the second error budget in the total of the first error budget and the second error budget.
3. The GNSS receiver of claim 1, wherein the first weight equals a hundred percent if the first error budget is lower than or equal to a combination of the second error budget and a first threshold value.
4. The GNSS receiver of claim 1, wherein the second weight equals a hundred percent if the second error budget is lower than a combination of the first error budget and a second threshold value.
5. The GNSS receiver of claim 1, wherein the second location uses, based on a decision criterion, one or more of: i) a calculated location of the GNSS receiver at an epoch before a current epoch as an initial value of the integration of the velocity data in the PVT vector or ii) a location of the GNSS receiver at a previous epoch input from positioning aids as an initial value of the integration of the velocity data in the PVT vector.
6. The GNSS receiver of claim 5, wherein the positioning aids comprise one or more of a map matching algorithm, an output of an inertial navigation system or a Wi-Fi localization system.
7. The GNSS receiver of claim 6, wherein the estimated error budgets depend upon one or more of a configuration of said GNSS receiver, a geometry of GNSS satellites in view of said GNSS receiver or a topology in an environment of said GNSS receiver.
8. The GNSS receiver of claim 6, wherein the estimated error budgets are estimated based on one or more of a Dilution Of Precision (DOP), a User Equivalent Range Error (UERE), a User Equivalent Range Rate Error (UERRE) or a probability estimate of a random walk of the integrated velocity along its path between a previous epoch and a current epoch.
9. The GNSS receiver of claim 8, wherein the first error budget equals the DOP multiplied by the UERE.
10. The GNSS receiver of claim 9, wherein the second error budget equals a product of the DOP multiplied by the UERRE, times the duration during which velocity has been integrated.
11. The GNSS receiver of claim 10, wherein the probability estimate of the random walk of the integrated velocity is modeled by an ellipsoid.
12. The GNSS receiver of claim 8, wherein one or more of the UERE and the UERRE are a combination of one or more of a parameter value, a carrier to noise ratio dependent value, an elevation dependent value or a signal dependent value.
13. The GNSS receiver of claim 8, wherein one or more of the first error budget or the second error budget are forecast estimated error budgets, said forecast being based on a combination of an estimated future trajectory of the GNSS receiver and a probability on said estimated future trajectory of one or more of a configuration of said GNSS receiver, a geometry of GNSS satellites in view of said GNSS receiver or a topology in an environment of said GNSS receiver.
14. The GNSS receiver of claim 8, wherein one or more of the UERE and the UERRE are predefined and stored in a database accessible by the first computing logic.
15. The GNSS receiver of claim 1, further comprising a third computing logic configured to calculate a heading information from the PVT vector.
16. A localization method using a GNSS receiver comprising:
calculating a first location and a second location of the GNSS receiver, the first location based on position data in a PVT vector calculated by the GNSS receiver, the second location based on a result of an integration of velocity data in the PVT vector;
estimating a first error budget of the first location and a second error budget of the second location and to cause the first computing logic to output a calculated location of the GNSS receiver as a weighted average of the first location and the second location;
wherein a first weight of the first location and a second weight of the second location in the calculated location, are calculated based on the first error budget and the second error budget.
US16/089,735 2016-04-01 2017-03-29 Improved gnss receiver using velocity integration Abandoned US20200309961A1 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP16305387.9A EP3226034A1 (en) 2016-04-01 2016-04-01 Improved gnss receiver using velocity integration
EP16305387.9 2016-04-01
PCT/EP2017/057377 WO2017167791A1 (en) 2016-04-01 2017-03-29 Improved gnss receiver using velocity integration

Publications (1)

Publication Number Publication Date
US20200309961A1 true US20200309961A1 (en) 2020-10-01

Family

ID=56087213

Family Applications (2)

Application Number Title Priority Date Filing Date
US16/089,735 Abandoned US20200309961A1 (en) 2016-04-01 2017-03-29 Improved gnss receiver using velocity integration
US16/089,814 Active 2037-07-23 US10942281B2 (en) 2016-04-01 2017-03-31 GNSS receiver using a combination of velocity integration and precise point positioning

Family Applications After (1)

Application Number Title Priority Date Filing Date
US16/089,814 Active 2037-07-23 US10942281B2 (en) 2016-04-01 2017-03-31 GNSS receiver using a combination of velocity integration and precise point positioning

Country Status (4)

Country Link
US (2) US20200309961A1 (en)
EP (1) EP3226034A1 (en)
CN (2) CN109313272B (en)
WO (2) WO2017167791A1 (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3025610B1 (en) * 2014-09-05 2022-01-14 Centre Nat Etd Spatiales METHOD FOR COLLABORATIVE DETERMINATION OF POSITIONING ERRORS OF A SATELLITE NAVIGATION SYSTEM
EP3226034A1 (en) 2016-04-01 2017-10-04 Centre National d'Etudes Spatiales Improved gnss receiver using velocity integration
CN108919321B (en) * 2018-05-18 2019-05-10 长安大学 A kind of GNSS positioning Detection of Gross Errors method based on trial and error method
EP3833997A1 (en) * 2018-08-07 2021-06-16 Taber Innovations Group LLC Personnel location and monitoring system
CN109787678A (en) * 2019-01-30 2019-05-21 西安寰宇赛特信息科技有限公司 A kind of communication positioning integrated terminal leading to satellite based on Beidou and day
CN111045046B (en) * 2019-12-03 2021-10-19 东南大学 Short-term ionosphere forecasting method and device based on NARX
CN111812697B (en) * 2020-05-25 2024-05-31 浙江航天润博测控技术有限公司 Combined positioning gesture measurement data processing method based on multi-mode precise positioning
CN112630812B (en) * 2020-11-30 2024-04-09 航天恒星科技有限公司 Multisource navigation positioning method
CN115542354A (en) * 2021-06-30 2022-12-30 博通集成电路(上海)股份有限公司 Apparatus and method for calculating receiver position-velocity-time result
US11451955B2 (en) 2021-09-01 2022-09-20 Autonomous Roadway Intelligence, Llc V2X and vehicle localization by local map exchange in 5G or 6G
CN116088015B (en) * 2023-04-11 2023-06-06 中国民用航空飞行学院 Method for improving track precision of small-sized training machine after accident without air-ground equipment reconstruction
CN116545517B (en) * 2023-07-07 2023-09-26 青岛国数信息科技有限公司 Stable and effective Beidou communication-in-motion antenna dual-mode tracking control method

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5375059A (en) 1990-02-05 1994-12-20 Caterpillar Inc. Vehicle position determination system and method
US5877723A (en) * 1997-03-05 1999-03-02 Caterpillar Inc. System and method for determining an operating point
FR2779234A1 (en) 1998-05-28 1999-12-03 Alsthom Cge Alcatel POSITIONING BY CALCULATION OF PSEUDO-SPEEDS IN A SATELLITE NAVIGATION SYSTEM
US6525688B2 (en) * 2000-12-04 2003-02-25 Enuvis, Inc. Location-determination method and apparatus
JP2005300167A (en) * 2004-04-06 2005-10-27 Alpine Electronics Inc Satellite positioning system, and navigation system
JP2006317241A (en) * 2005-05-12 2006-11-24 Furuno Electric Co Ltd Positioning device and positioning method
US8064923B2 (en) * 2007-04-28 2011-11-22 Broadcom Corporation Wireless communications device and integrated circuits with global positioning and method for use therewith
JP4453728B2 (en) * 2007-08-07 2010-04-21 株式会社デンソー Position correction device
US9366763B2 (en) * 2009-02-04 2016-06-14 Qualcomm Incorporated Method and apparatus for position determination with hybrid SPS orbit data
US8368590B2 (en) * 2008-08-19 2013-02-05 Trimble Navigation Limited GNSS signal processing methods and apparatus with candidate set selection
US8531333B2 (en) * 2009-12-10 2013-09-10 Maxlinear, Inc. Intermittent tracking for GNSS
JP5521531B2 (en) * 2009-12-18 2014-06-18 セイコーエプソン株式会社 Position calculation method and position calculation system
US9829582B2 (en) * 2011-09-19 2017-11-28 Raytheon Company Method and apparatus for differential global positioning system (DGPS)-based real time attitude determination (RTAD)
IL216475A (en) * 2011-11-20 2016-09-29 Intel Corp Navigation system and method with location-aware accuracy and/or power adjustments
US8972119B2 (en) * 2013-03-15 2015-03-03 Novatel Inc. System and method for heavy equipment navigation and working edge positioning
FR3025610B1 (en) 2014-09-05 2022-01-14 Centre Nat Etd Spatiales METHOD FOR COLLABORATIVE DETERMINATION OF POSITIONING ERRORS OF A SATELLITE NAVIGATION SYSTEM
EP3226034A1 (en) 2016-04-01 2017-10-04 Centre National d'Etudes Spatiales Improved gnss receiver using velocity integration

Also Published As

Publication number Publication date
US20190094385A1 (en) 2019-03-28
EP3226034A1 (en) 2017-10-04
CN109313272A (en) 2019-02-05
CN109196382A (en) 2019-01-11
WO2017167791A1 (en) 2017-10-05
US10942281B2 (en) 2021-03-09
CN109196382B (en) 2023-12-29
WO2017167943A1 (en) 2017-10-05
CN109313272B (en) 2023-10-27

Similar Documents

Publication Publication Date Title
US20200309961A1 (en) Improved gnss receiver using velocity integration
US7724184B2 (en) System and method for detecting false navigation signals
US8447517B2 (en) Tightly-coupled GNSS/IMU integration filter having speed scale-factor and heading bias calibration
US8521179B2 (en) Mobile unit's position measurement apparatus and mobile unit's position measurement method
US8570214B2 (en) Position calculating method and position calculating device
US20080129586A1 (en) Satellite-Based Positioning Receiver with Improved Integrity and Continuity
US10012737B2 (en) Method for estimating the level of error in satellite geolocation measurements and for monitoring the reliability of said estimations and associated device
KR102188880B1 (en) Terminal, base station and location positioning method
US20180011200A1 (en) Satellite signal exclusion based on doppler information
US9817126B2 (en) Methods for identifying whether or not a satellite has a line of sight
US10830898B2 (en) Method and apparatus applicable to positioning in NLOS environment
US20110206090A1 (en) Method for Optimizing an Acquisition of a Spread-Spectrum Signal from a Satellite by a Mobile Receiver
Musumeci et al. Design of a very high sensitivity acquisition system for a space GNSS receiver
KR102428135B1 (en) Method for estimating multipath error of pseudo-range measurement and positioning method using the same
US11747487B2 (en) GNSS receiver clock frequency drift detection
Preston GPS Multipath Detection and Mitigation Timing Bias Techniques
JP4518096B2 (en) Mobile positioning device
US11493638B2 (en) Carrier phase based multipath (MP) detection
Wang et al. Research on the incomplete constellation GPS positioning algorithm by a combined altitude and clock bias model
US20200386896A1 (en) Single-epoch pseudo-range positioning under varying ionosphere delays
US20240027630A1 (en) Real-time ppe base station measurement uncertainty modeling for protection level computation
Kubo et al. Performance Evaluation of PPP/PPP-RTK in Semi-Urban Areas
GB2612114A (en) Positioning system and method
Sanromà Sánchez Positioning in urban environments
Castanié Safe navigation for vehicles

Legal Events

Date Code Title Description
AS Assignment

Owner name: CENTRE NATIONAL D'ETUDES SPATIALES, FRANCE

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:MARMET, FRANCOIS-XAVIER;REEL/FRAME:047578/0511

Effective date: 20181114

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: FINAL REJECTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE AFTER FINAL ACTION FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: ADVISORY ACTION MAILED

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION