US20210181357A1 - System and method for validating gnss ambiguities - Google Patents
System and method for validating gnss ambiguities Download PDFInfo
- Publication number
- US20210181357A1 US20210181357A1 US17/119,823 US202017119823A US2021181357A1 US 20210181357 A1 US20210181357 A1 US 20210181357A1 US 202017119823 A US202017119823 A US 202017119823A US 2021181357 A1 US2021181357 A1 US 2021181357A1
- Authority
- US
- United States
- Prior art keywords
- ambiguity
- ambiguities
- determining
- fixed
- satellite
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/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/23—Testing, monitoring, correcting or calibrating of receiver elements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Definitions
- This invention relates generally to the satellite positioning field, and more specifically to a new and useful system and method in the satellite positioning field.
- a receiver position is determined based on a set of satellite observations within a given time frame. These satellite observations are used to calculate the receiver position.
- satellite observations can, for example, include outliers and noisy data that can decrease the positioning accuracy. To achieve high levels of precision in the position, it can be desirable to validate the result.
- This invention provides such new and useful system and method for validating GNSS ambiguities.
- FIG. 1 is a schematic representation of the system.
- FIG. 2 is a schematic representation of the method.
- FIG. 3 is a schematic representation of an example of transforming the state vectors.
- FIG. 4A is a block diagram representation of an example of the method.
- FIG. 4B is a block diagram representation of an example of validating the state vectors.
- a method includes receiving satellite observations, determining state vectors, validating the state vectors, and determining the receiver position.
- the method can optionally include correcting the satellite observations.
- Validating the state vectors can include determining one or more transformations, transforming the ambiguities, and comparing the ambiguities.
- the system includes one or more receivers and a computing system.
- the system can optionally include one or more sensors, and/or any suitable components.
- the system is preferably communicably connected to one or more data sources (e.g., reference stations, satellites, satellite constellations, etc.).
- the system and method preferably function to estimate a position of a GNSS receiver and/or an external system.
- Embodiments of the system and/or method can be used, for example, in autonomous or semi-autonomous vehicle guidance (e.g., for unmanned aerial vehicles (UAVs), unmanned aerial systems (UAS), self-driving cars, agricultural equipment, robotics, rail transport/transit systems, autonomous trucking, last mile delivery, etc.), GPS/GNSS research, surveying systems, user devices, mobile applications, internet-of-things (TOT) devices, and/or may be used in any other suitable application.
- UAVs unmanned aerial vehicles
- UAS unmanned aerial systems
- UAS unmanned aerial systems
- self-driving cars agricultural equipment
- robotics robotics
- rail transport/transit systems autonomous trucking, last mile delivery, etc.
- GPS/GNSS research GPS/GNSS research
- surveying systems e.g., surveying systems, user devices, mobile applications, internet-of-things (TOT) devices
- system can be coupled to any suitable external system such as a vehicle (e.g., UAV, UAS, car, truck, etc.), robot, railcar, user device (e.g., cell phone), and/or any suitable system, and can provide positioning data, integrity data (e.g., protection level data), and/or other data to said system.
- a vehicle e.g., UAV, UAS, car, truck, etc.
- robot e.g., robot
- railcar e.g., cell phone
- user device e.g., cell phone
- integrity data e.g., protection level data
- the inventors have discovered that in some applications, to achieve high accuracy and/or high integrity positioning solutions based on GNSS techniques, in addition to validation steps performed when fixing the integer ambiguities, validation steps between integer ambiguity sets (e.g., fixed integer ambiguities) can be beneficial.
- the inter-ambiguity set validation can increase a confidence in the resulting ambiguities and can lead to more accurate and/or higher integrity positioning solutions.
- variants of the technology can be used independent of the data sources (e.g., satellites, satellite constellations, satellite frequencies, sensors, reference stations, etc.) and/or relationships between or within data collected from data sources.
- the technology can validate the state vectors (e.g., state vectors that are determined for different data sets and/or by different filters), regardless of what transformation(s) are applied to each state vector.
- variants of the technology can determine an optimal or near-optimal number of states that can be compared between state vectors, which can improve the probability of successful validation between the state vectors, can provide an optimal or near-optimal set of validated states, and/or otherwise improve the state vector determination.
- the system preferably functions to determine the position of a receiver and/or external system, determine a plurality of state vectors, and validate the state vectors.
- the system preferably uses a set of data collected by one or more data sources.
- Data sources can include: receivers, sensors (e.g., located onboard the receiver, the external system, a reference stations, etc.), databases, reference stations, satellites, and/or any other suitable data source. Examples of data that can be used include: satellite observations, sensor observations, and/or any other suitable data.
- the receiver e.g., mobile receiver
- the receiver preferably functions to receive a set of satellite observations (e.g., satellite signals) from one or more satellites.
- the receiver can determine the location (e.g., pseudorange) of the receiver (and/or external system) based on the satellite observations.
- the receiver is preferably in communication with the computing system.
- the receiver can be integrated with the computing system, and/or the receiver and computing system can be arranged in any suitable manner.
- the receiver is preferably a stand-alone device (e.g., a GNSS receiver, antenna).
- the receiver can be integrated into the external system (e.g., be a component of an automobile, aero vehicle, nautical vehicle, etc.), can be a user device (e.g., smart phone, laptop, cell phone, smart watch, etc.), and/or can be configured in any suitable manner.
- the external system e.g., be a component of an automobile, aero vehicle, nautical vehicle, etc.
- a user device e.g., smart phone, laptop, cell phone, smart watch, etc.
- the set of satellite observations can include orbital data (e.g., ephemeris), timestamp, code data, carrier phase data, pseudocode data, and/or any suitable data.
- the set of satellite observations preferably includes satellite observations corresponding to satellites from more than one satellite constellation (e.g., Global Positioning System (GPS), GLObal Navigation Satellite System (GLONASS), BeiDou positioning System (BDS), Galileo, etc.).
- GPS Global Positioning System
- GLONASS GLObal Navigation Satellite System
- BDS BeiDou positioning System
- Galileo Galileo
- the set of satellite observations can correspond to satellites from a single satellite constellation, can include data from an augmentation system (e.g., Satellite Based Augmentation System (SBAS) such as Wide Area Augmentation System (WAAS), European Geostationary Navigation Overlay Service (EGNOS), Multi-Functional Satellite Augmentation System (MSAS), Omnistar, StarFire, etc.; Ground Based Augmentation Systems (GBAS) such as Local Area Augmentation System (LAAS); etc.), and/or can include any suitable data.
- SBAS Satellite Based Augmentation System
- WAAS Wide Area Augmentation System
- GNOS European Geostationary Navigation Overlay Service
- MSAS Multi-Functional Satellite Augmentation System
- GBAS Ground Based Augmentation Systems
- LAAS Local Area Augmentation System
- Each satellite observation from the set of satellite observations preferably corresponds to a common time window (e.g., epoch).
- each satellite observation can be associated with a time stamp (e.g., time of transmission, time of receipt
- each receiver can be configured to receive satellite observations associated a satellite constellation, a carrier frequency (e.g., the L1, L2, L5, E1, E5a, E5b, Eab, E6, G1, G3, B1, B2, B3, LEX, etc. frequencies), and/or corresponding to any suitable data.
- a receiver can include a first antennae configured to receive satellite signals from GPS satellites and a second antennae configured to receive satellite signals from Galileo satellites.
- the computing system preferably functions to process the data (e.g., satellite observations, reference station observations, sensor data, etc.) received by the receiver.
- the computing system can: aggregate the data (e.g., combine the receiver satellite observations, reference station satellite observations, satellite corrections, and/or sensor data; reorganize the receiver satellite observations, reference station satellite observations, and/or sensor data such as based on the time stamp, time of transmission, time of receipt, etc.; etc.), filter the data (e.g., to calculate state vectors, ambiguities such as phase ambiguities, etc.
- the computing system can be local (e.g., to the external system, to the receiver, to the sensor, etc.), remote (e.g., cloud computing, server, networked, etc.), and/or otherwise distributed.
- the computing system is preferably communicably coupled to the receiver and/or the sensors, but can be communicably coupled to any suitable data sources.
- the computing system is preferably colocalized with (e.g., integrated into) the receiver (and/or external system), but the computing system can be remote from the receiver (and/or external system), and/or configured in any suitable manner.
- the sensors preferably function to measure sensor data (e.g., auxiliary data) associated with the external system (and/or the receiver).
- the sensor data can be used to assist (e.g., speed-up, correct, refine, etc.) the calculation (e.g., calculating the state vector, estimating the phase ambiguity).
- the sensor data can be used to estimate the receiver position at times when satellite observations are not received (e.g., in an urban canyon; due to an obstruction such as a billboard, weather, etc.; etc.).
- the sensor data can be used at any suitable time and in any suitable manner.
- the sensors are preferably in communication with the computing system.
- the sensors can be: on-board the external system, on-board a separate external system, integrated into the receiver, separate from the receiver, and/or otherwise associated with the receiver.
- the sensor data can include: inertial data (e.g., velocity, acceleration), odometry, pose (e.g., position, orientation), mapping data (e.g., images, point clouds), temperature, pressure, ambient light, and/or any other suitable data.
- the sensors can include one or more of: inertial moment unit (IMU), accelerometer, gyroscope, magnetometer, odometer (e.g., visual odometry, wheel odometry, etc.), and/or any suitable sensor.
- IMU inertial moment unit
- the method functions to determine the position of a receiver (e.g., an external system coupled to the receiver) based on one or more sets of satellite observations.
- the method preferably validates the position of the receiver and/or one or more intermediate (e.g., state vector) used to determine the position of the receiver.
- One or more instance of the method and/or steps of the method can be performed in series and/or in parallel (e.g., concurrently).
- Receiving one or more satellite observations S 100 preferably functions to measure (e.g., at the receiver, at reference stations, etc.) and/or access one or more set of satellite observations (e.g., carrier phase measurements. pseudo-range measurements, etc.) from one or more observed satellites.
- the satellite observations can be measured and/or received by a receiver, retrieved from a database (e.g., retrieve stored satellite observations; retrieve stored corrections; retrieve an almanac such as of weather conditions, tides, etc.; etc.), and/or be otherwise received.
- S 100 can include receiving Doppler measurement data, sensor data, and/or any suitable data.
- the satellite observations can include signals from one or more satellite constellations. Each set of satellite observations preferably corresponds to a plurality of satellite constellations. However, one or more set of satellite observations can correspond to a single satellite constellation.
- S 100 is preferably performed before correcting the satellite observations S 200 ; however, S 100 can be performed at the same time as S 200 .
- S 100 is preferably performed before determining the state vectors S 300 ; however, S 100 and S 300 can be performed at the same time and/or S 300 can be performed before S 100 (e.g., during different time windows).
- S 100 is preferably performed by a receiver. However, S 100 can additionally or alternatively be performed by one or more reference stations, by one or more sensors, and/or by any suitable component.
- the data (e.g., satellite observations, sensor data, etc.) received in S 100 is preferably measured during a time window (e.g., time period, epoch).
- a time window e.g., time period, epoch.
- Each of the sets of satellite observations can correspond to the same and/or to a different time window.
- the method can include repeating S 100 during a second time window.
- the second time window can be the same duration as the time window and/or a different duration.
- the second time window preferably immediately follows the time window; however, the second time window can be delayed relative to the time window by any amount.
- the (first and second) time windows are preferably separate and distinct (e.g., consecutive, nonconsecutive, etc.), but can alternatively overlap.
- the second time window can additionally or alternatively be started in response to a trigger event (e.g., convergence of one or more calculation; collective convergence of the calculations; receiver position; loss of satellite observations such as for a specific satellite, for an amount of time, etc.; etc.), at a predetermined time, and/or at any suitable time.
- a trigger event e.g., convergence of one or more calculation; collective convergence of the calculations; receiver position; loss of satellite observations such as for a specific satellite, for an amount of time, etc.; etc.
- the time window can be extended and/or any suitable action can occur in response to the trigger event.
- Correcting the satellite observations S 200 preferably functions to correct the satellite observations for local and/or global errors.
- S 200 is preferably performed by a computing system (e.g., a correction module of a computing system), but can be performed by any system or component.
- the corrections are preferably received from a computing system (e.g., a remote computing system, a cloud computing system, corrections service, etc.), but can additionally or alternatively be generated by the computing system (e.g., based on reference station data, based on almanac data, based on a model, etc.), received from a reference station, and/or be otherwise accessed.
- S 200 is preferably performed before determining the state vector S 300 . However, S 200 can be performed at the same time as and/or after S 300 .
- S 200 preferably corrects all of the satellite observations within each set of satellite observations. However, any suitable subset of satellite observations can be corrected. Corrections can be applied in response to (e.g., to correct): hardware bias, timing errors (e.g., satellite clock, receiver clock, reference station clock, etc.), atmospheric effects (e.g., ionosphere effects, tropospheric effects), relativistic effects, local effects (e.g., multipath errors), global effects, noise (e.g., outliers), and/or any suitable effects.
- reference station satellite observations can be used to correct the receiver satellite observations such as by calculating double differenced satellite observations.
- effects e.g., atmospheric effects
- the satellite observations can be corrected.
- correcting the satellite observations can be performed in any suitable manner.
- the corrections can be generated as disclosed in and/or received from a corrections service as disclosed in U.S. patent application Ser. No. 16/589,932 filed 1 Oct. 2019 entitled “SYSTEMS AND METHODS FOR DISTRIBUTED DENSE NETWORK PROCESSING OF SATELLITE POSITIONING DATA,” which is incorporated in its entirety by this reference.
- the corrections can be otherwise generated.
- Determining the state vector S 300 preferably functions to process the data (e.g., satellite observations, corrected satellite observations, sensor data) to determine phase ambiguities (e.g., fixed integer phase ambiguities) associated with each set of satellite observations.
- States of the state vector preferably include a carrier phase ambiguity (e.g., calculated as float ambiguities, calculated as an integer ambiguity, etc.) associated with one or more satellites (e.g., each satellite that provided one or more observation in the set of satellite observations) and phase ambiguity covariances (and/or correlations) associated with each one or more pairs of satellites (e.g., pairs of satellites that provided observations in the set of satellite observations).
- the state vector can additionally or alternatively include residuals (e.g., for the phase ambiguities, for the covariances), a validity of one or more states, a receiver position, an external system location, pseudorange, atmospheric error or effect, estimated error (e.g., estimated error in one or more states in the state vector), derivatives of one or more states (e.g., first, second, etc. derivative such as with respect to time, with respect to position, etc.), and/or any suitable state.
- S 300 is preferably performed by a computing system (e.g., a processing module of a computing system), but can be performed by any component.
- S 300 is preferably performed before validating the state vectors S 400 .
- S 300 can be performed at the same time as and/or after S 400 (e.g., when a subset of the state vector is validated and said subset is used to determine the remainder of the state vector).
- S 300 is preferably performed twice for each instance of the method (e.g., to generate a first state vector and a second vector associated with a first and second set of satellite observations, respectively).
- S 300 can be performed once for each instance of the method (e.g., wherein one instance of the method is compared or validated relative to another instance of the method), and/or can be performed any number of times.
- S 300 can include: determining a set of floating phase ambiguity hypotheses, determining a set of integer phase ambiguity hypotheses from the set of floating phase ambiguity hypotheses, performing hypothesis testing to validate the set of integer phase ambiguities hypotheses, and/or any suitable steps.
- S 300 can include steps and/or be performed as disclosed in U.S. patent application Ser. No. 16/685,927 filed 15 Nov. 2019 entitled “SYSTEM AND METHOD FOR SATELLITE POSITIONING” and/or U.S. patent application Ser. No. 16/817,196 filed 12 Mar. 2020 entitled “SYSTEMS AND METHODS FOR REAL TIME KINEMATIC SATELLITE POSITIONING,” each of which is incorporated in its entirety by this reference.
- Determining a set of floating phase ambiguity hypotheses preferably functions to determine a phase ambiguity (typically in a floating point representation, but other representations can be used) associated with each satellite represented in the set of satellite observations.
- the set of floating phase ambiguity hypotheses are preferably determined using a filter.
- the filter is preferably a Kalman filter.
- the filter can be an extended Kalman filter, an unscented Kalman filter, a Bierman-Thornton filter, a particle filter, Monte Carlo simulation, and/or any suitable mean-square error filter and/or sensor fusion algorithm.
- the inputs to the filter can include: satellite observations (e.g., receiver satellite observations, corrected satellite observations, reference station satellite observations, etc.), sensor data, and/or any suitable data.
- the output of the filter is preferably a state vector and/or one or more states of a state vector, but can include any suitable information.
- Determining a set of integer phase ambiguity hypotheses preferably functions to constrain (and/or otherwise determine from) the set of floating phase ambiguity hypotheses to a set of integer phase ambiguity hypotheses. Integer phase ambiguities can be beneficial for improving an accuracy of the determined receiver position as compared to the receiver position determined using floating ambiguities. Determining the set of integer phase ambiguity hypotheses can include reducing a correlation between ambiguities of the set of phase ambiguities (e.g., floating phase ambiguity hypotheses) and performing a search (e.g., a least squares search, an absolute value search, etc.) to identify possible integer phase ambiguities to include in the set of integer phase ambiguity hypotheses.
- a search e.g., a least squares search, an absolute value search, etc.
- Least-Squares Ambiguity Decorrelation Adjustment LAMBDA
- modified-LAMBDA MLAMBDA
- LLL reduction algorithm a whitening transformation, a coloring transformation, a decorrelation transformation, rounding
- the Ambiguity Function Method AFM
- FARA Fast Ambiguity Resolution Approach
- LSAST Least-Squares Ambiguity Search Techniques
- integer bootstrapping and/or any suitable decorrelation or reduction algorithm
- the set of integer phase ambiguity hypotheses preferably includes integer phase ambiguities that satisfy a hypothesis criterion, but can include any suitable integer phase ambiguities.
- Example of hypotheses criterion include a threshold sum of least squares, a threshold sum of absolute differences, a threshold likelihood of the correct solution, and/or any suitable criterion.
- the set of integer phase ambiguity hypotheses can be generated in any manner.
- Performing hypothesis testing preferably functions to determining the integer phase ambiguities from the set of integer phase ambiguity hypotheses that are most likely to be correct.
- the hypothesis test is preferably a Bayesian inference, but can additionally or alternatively be based on a statistical confidence, significance testing, and/or any suitable hypothesis testing. Examples of hypothesis tests can include difference test, ratio test, projector test, f-test GIA test, and/or any suitable hypothesis test.
- performing hypothesis testing can include determining a probability (e.g., a likelihood, a log likelihood, etc.) associated with subsets of integer phase ambiguities of the integer phase ambiguity hypotheses; computing a ratio of the probability between two subsets of integer phase ambiguities (e.g., within the same ambiguity set, across different ambiguity sets); when the ratio between a most-likely and a next most likely subset of integer phase ambiguities exceeds a threshold, storing the most-likely subset of integer phase ambiguities as the integer phase ambiguities and ceasing hypothesis testing.
- a probability e.g., a likelihood, a log likelihood, etc.
- the most-likely subset of integer phase ambiguities can be referred to as internally-validated integer phase ambiguities (e.g., integer phase ambiguities that are validated against the remaining integer phase ambiguities from the same set of integer phase ambiguity hypotheses).
- the internally-validated integer phase ambiguities can be contrasted with cross-validated integer phase ambiguities (e.g., integer phase ambiguities as validated in step S 400 ).
- S 300 can be performed more than once to generate at least two state vectors.
- Each state vector can be determined consecutively, iteratively, sequentially, delayed, and/or with any suitable timing relative to the other state vectors.
- the plurality of state vectors is preferably independent of one another, but can be correlated, dependent on each other, and/or otherwise be related.
- Each state vector can include a subset of common states (e.g., common integer phase ambiguities such as those associated with the same data source), a subset of distinct states (e.g., integer phase ambiguities associated with different data sources, uncommon states, etc.), and/or the states in the first and second state vectors can be otherwise related.
- Different state vectors can differ in: satellite composition, frequency, sampling epoch, otherwise differ, or not differ (e.g., share all or a portion of satellite observations).
- the states within each state vector can share: satellites, frequencies, sampling epochs, other parameters, or no parameters.
- two (or more) state vectors can be associated with (e.g., generated from) a two (or more) distinct sets of satellite observations (e.g., a first and a second set of satellite observations).
- the distinct sets of satellite observations preferably include observations from at least three common satellites, but can include observations from any number of common satellites.
- the first and second set of satellite observations can correspond to sets of satellite observations that are sent to distinct Kalman filters as disclosed in U.S.
- two (or more) state vectors can be generated using two (or more) different filters (and/or methods).
- a first state vector can be generated using phase hypotheses generated using a Kalman filter and a second state vector can be generated using phase hypotheses generated using a Bierman-Thornton filter.
- a first and second state vector can each be generated with different initial state guesses.
- the state vectors can be generated using any suitable data and/or any suitable method.
- the state vector(s) can be associated with one or more transformations.
- the transformation can improve the phase ambiguity determination (e.g., enable the phase ambiguity to be determined quickly, enable the phase ambiguity to be determined with greater accuracy, etc.), improve the receiver position estimate (e.g., improve the accuracy or integrity of the determined receiver position), and/or otherwise be used.
- the output of the filter e.g., the state vector
- the transformed phase ambiguity can be determined from the state vector (e.g., by transforming the phase ambiguities of the state vector), and/or the transformed phase ambiguity can be otherwise related to the state vector.
- transformations include: single differenced combinations, double differenced combinations, triple differenced combinations, wide lane combinations, middle lane combinations, narrow lane combinations, extra-wide lane combinations, ionosphere free combinations, geometry free combinations, noise reduced combinations, ionospheric delay reduced combinations, noise reducing combinations, linear combinations (e.g., between two or more carrier frequencies of satellite observations, between two or more data sources such as satellite observations and reference station observations, etc.), and/or any suitable combination of data source observations.
- linear combinations e.g., between two or more carrier frequencies of satellite observations, between two or more data sources such as satellite observations and reference station observations, etc.
- Each transformation is preferably an n ⁇ m matrix where n is less than or equal to m, where m corresponds to the number of ambiguities in the respective ambiguity set (e.g., ambiguity vector, state vector, etc.), and where n corresponds to the number of ambiguities in the transformed space.
- n can be greater than m.
- the elements of the transformation matrices are preferably integers, but can be rational, real, complex, and/or have any values.
- the transformations can include: transforming to wide-lane ambiguities, transforming to decrease the correlation between the ambiguities, transforming to a rotated state space, transforming to a scaled state space, transforming to a truncated state space, and/or transforming in any suitable manner.
- the transformed phase ambiguity and the underlying phase ambiguity are preferably both integers, but either can be rational numbers, floating values, real valued, and/or otherwise by related. However, the transformed phase ambiguity can be otherwise determined.
- each state vector can be associated with the same or a distinct transformation.
- Validating the state vectors S 400 preferably functions to cross validate one or more states from two or more state vectors (e.g., determine what states of two or more state vectors match). S 400 preferably functions to cross-validate a subset of integer phase ambiguities (e.g., internally-validated integer phase ambiguities) from the first state vector with a subset of integer phase ambiguities from the second state vector. However, any suitable state(s) can be cross-validated between the first and second state vector. S 400 is preferably performed by a computing system (e.g., a validation module of a computing system). However, S 400 can be performed by any suitable component. S 400 is preferably performed before determining the receiver position S 500 .
- a computing system e.g., a validation module of a computing system.
- S 400 can be performed by any suitable component. S 400 is preferably performed before determining the receiver position S 500 .
- S 400 can be performed at the same time as and/or after (e.g., in examples where the receiver position is used to validate the state vectors) S 500 .
- S 400 is preferably performed after S 300 , such that the states (e.g., the phase ambiguities) that S 400 is analyzing are fixed; however, S 400 can optionally analyze unfixed states (e.g., unfixed ambiguities).
- the state vectors are preferably validated based on a collective convergence criterion (e.g., a validation metric, validation score).
- the collective convergence criterion can include: one or more subsets of transformations associated with the state vectors matching, subsets of ambiguities matching across the state vectors, a residual threshold, one or more states (e.g., from each of the calculations) matching, a timeout, a covariance threshold, a receiver position threshold, and/or any suitable criterion.
- S 400 can include: combining the first and second state vectors 420 , determining an intersection transformation S 450 , comparing the state vectors in the intersection space S 480 , and/or any suitable steps.
- the first and second state vectors preferably at least partially share states from satellite observations from common: satellites, frequencies, and/or epochs, but can alternatively include entirely unrelated satellite observations.
- the first and second state vectors share states from observations from common satellites, frequencies, and epochs.
- the first and second state vectors share states from observations from common satellites and frequencies, but different epochs.
- the S 420 preferably functions to combine the first and second state vector to generate a combined state vector.
- the combined state vector preferably includes the non-transformed integer phase ambiguities (e.g., internally-validated integer phase ambiguity, fixed integer phase ambiguities, etc.) from the first and second state vector, but can additionally or alternatively include the transformed integer phase ambiguities, a subset of integer phase ambiguities, and/or any suitable states from the first and/or second state vectors.
- the combined state vector can be an augmented state vector (e.g., including all or a majority of integer phase ambiguities associated with the first and the second state vector), a common state vector (e.g., including the common integer phase ambiguities from the first and second state vector), an uncommon state vector (e.g., including uncommon ambiguities or unshared ambiguities), and/or any suitable state(s) can be included in the combined state vector.
- an augmented state vector e.g., including all or a majority of integer phase ambiguities associated with the first and the second state vector
- a common state vector e.g., including the common integer phase ambiguities from the first and second state vector
- an uncommon state vector e.g., including uncommon ambiguities or unshared ambiguities
- any suitable state(s) can be included in the combined state vector.
- S 420 can be performed using two fixed integer phase ambiguity sets.
- the fixed integer phase ambiguity sets are preferably not linearly independent (e.g., include one or more common fixed integer phase ambiguity, two fixed integer phase ambiguity sets, a 1 and a 2 can be expressed as
- each fixed integer phase ambiguity set can be associated with a transformation, T 1 and T 2 associated with the first and second integer phase ambiguity set, respectively.
- the transformations can be written as
- the transformed ambiguities can have any suitable form.
- the augmented integer ambiguity set (e.g., augmented integer ambiguity space, augmented state vector, etc.) preferably functions to ensure that the transformed states have the same dimensions.
- the augmented integer ambiguity set can be written as
- a ⁇ [ a 1 ′ a 2 ′ a c ] .
- the augmented integer ambiguity set can be associated with a first and second augmented transformation that transforms the augmented integer ambiguity set to the first and second integer ambiguity sets.
- the augmented transformations can be written as:
- the augmented transformations ⁇ circumflex over (T) ⁇ 1 and ⁇ circumflex over (T) ⁇ 2 are preferably re-ordered to align with the states of the augmented ambiguity vector.
- the augmented integer ambiguity set can include any suitable states and/or be arranged in any manner.
- Determining an intersection transformation S 450 functions to determine a transformation that transforms the state vectors to a shared intersection space.
- the intersection space preferably functions to enable the ambiguities to be compared in the same basis while retaining as many linearly independent rows or ambiguities as possible.
- the intersection space is preferably determined using the augmented transformations (augmented ambiguity transformations). However, the intersection space can be determined from the transformed ambiguities (e.g., z i ), the initial ambiguities (e.g., a i ), the initial transformations (e.g., T i ), the augmented integer phase ambiguities (e.g., a), and/or any suitable vector and/or transformation can be used.
- the row space of T i , row[ ⁇ circumflex over (T) ⁇ i ] is preferably the set of vectors that result from the possible combinations of rows of matrix T i .
- the row space of the intersection of T 1 and T 2 , row[ ⁇ circumflex over (T) ⁇ ⁇ ], can then be determined from (e.g., is the same as) the intersection of the row space of T 1 with the row space of T 2 , where the intersection of T 1 and T 2 can have the minimum number of rows that is possible while the row space condition is adhered to.
- the intersection transformation can be otherwise defined and/or determined.
- the intersection space can be determined, for example, by concatenating the transformations into a concatenated matrix (e.g., S), and searching for the null space of S.
- S a concatenated matrix
- S [ ⁇ circumflex over (T) ⁇ 1 T ⁇ circumflex over (T) ⁇ 2 T ], wherein
- the elements of S are preferably integers. However, the elements of S can be natural, rational, real, and/or complex valued.
- Each of the state vectors can be transformed to the intersection space, a subset of the state vectors (e.g., a 1 and a 2 , z 1 and z 2 , etc.) can be transformed to the intersection space, and/or any suitable states can be transformed to the intersection space.
- the second set of fixed carrier ambiguities using the intersection transformation ⁇ circumflex over (T) ⁇ ⁇ V 2 ⁇ circumflex over (T) ⁇ 2 .
- the ambiguity vectors can be transformed to any suitable space using any suitable transformation.
- Comparing the state vectors in the intersection space S 480 functions to determine whether the state vectors (or subsets thereof) match.
- the state vectors are preferably compared in the intersection space, but can additionally or alternatively be compared in the transformed space (e.g., comparing z 1 and z 2 ), in the untransformed space (e.g., comparing a 1 and a 2 ), in a different comparison space, and/or in any suitable space.
- S 480 preferably compares the integer phase ambiguities associated with the same data source between the two state vectors. However, any suitable states can be compared.
- each row of the first state vector can be compared to the row of the second state vector with the same row index to identify matches between the respective integer phase ambiguities.
- common ambiguities can be identified in the transformed space (e.g., by identifying the ambiguities determined from the same satellite observations, such as by using a unique satellite observation identifier) and compared.
- the state vectors (and/or subsets thereof) preferably match (e.g., are cross-validated) when a validation criterion is satisfied.
- the validation criterion e.g., collective convergence criterion
- the validation criterion preferably includes the number of differing ambiguities (e.g., between one or more integer phase ambiguity in the intersection space, transformed space, untransformed space, etc.) being less than a threshold, but can additionally or alternatively include a residual threshold, one or more states (e.g., from each of the state vectors) matching (e.g., having an identical value such as in the intersection space, the transformation space, the untransformed space, etc.), a difference between two ambiguities (e.g., a difference threshold), a timeout, a quorum of matching transformed ambiguities, and/or any suitable criterion.
- the ambiguities can be compared in any suitable manner.
- the validation criterion can be based on (e.g., depend on) the number of ambiguities (e.g., a threshold number of ambiguities), depend on the data source(s) (e.g., number of satellites, number of satellite constellations, identity of the satellite constellation(s), etc.), and/or depend on any suitable properties.
- the first and/or second state vectors can be validated (e.g., cross-validated; considered validated or cross-validated), a subset of the state vectors (e.g., the state(s) that meet the validation criterion, common states, uncommon states, etc.) can be validated (e.g., cross-validated), and/or any suitable state(s) can be validated, where validated ambiguities can subsequently be used to correct the carrier measurements and/or otherwise be used to determine the receive and/or external system location.
- the state vectors e.g., cross-validated; considered validated or cross-validated
- a subset of the state vectors e.g., the state(s) that meet the validation criterion, common states, uncommon states, etc.
- any suitable state(s) can be validated, where validated ambiguities can subsequently be used to correct the carrier measurements and/or otherwise be used to determine the receive and/or external system location.
- a new intersection space can be determined, a new set of state vectors can be compared, a flag can be issued that the state vectors cannot be cross-validated, transmitting an error message, calculating the receiver position (e.g., with non-validated ambiguities), starting a new time window, extending the time window, adding states (e.g., ambiguities) to the ambiguity vector(s), nullifying the ambiguity (and/or state) vectors, the external system can be operated at a lower accuracy and/or integrity level location determined, additional states can be determined (e.g., by performing S 300 ), the method can be reinitiated, and/or any suitable response can be performed.
- states e.g., ambiguities
- S 480 can include determining transformations using a different pair of ambiguity vectors, reporting an error, starting a new time window (e.g., restarting the method, returning to S 100 ), extending the current time window, starting a new filter (e.g., determining a new state vector such as in S 300 ), modifying a filter (e.g., change inputs to the filter), and/or can perform any suitable steps.
- matching preferably corresponds to a matching criterion (e.g., score, metric, quality).
- the matching criterion can be a threshold (e.g., a maximum difference between each value in V 1 ⁇ circumflex over (z) ⁇ 1 and V 2 ⁇ circumflex over (z) ⁇ 2 such as each row, each column, each element, etc.), a threshold distance, a residual (e.g., a target residual, a maximum residual, etc.), and/or any suitable value.
- a threshold e.g., a maximum difference between each value in V 1 ⁇ circumflex over (z) ⁇ 1 and V 2 ⁇ circumflex over (z) ⁇ 2 such as each row, each column, each element, etc.
- a threshold distance e.g., a target residual, a maximum residual, etc.
- a residual e.g., a target residual, a maximum residual, etc.
- matching can require the two state vectors (e.g., corresponding ambiguities, ambiguities in the intersection space, transformed ambiguities, etc.) to be the same.
- the cross-validated fixed carrier ambiguity set preferably includes at least 3 carrier phase ambiguities, but can contain fewer than 3 carrier phase ambiguities and/or any number of carrier phase ambiguities.
- the subset of carrier ambiguities that do not match can be nullified (e.g., cleared from the state vectors), can be cross-validated (e.g., considered valid at the same time the remaining states are considered valid), undergo additional validation testing, and/or otherwise be processed.
- any suitable response can be taken in response to a subset of the ambiguity vectors not matching.
- the fixed ambiguity vector e.g., validated ambiguity vector
- the fixed ambiguity vector can be determined in any suitable manner.
- Determining the receiver position S 500 preferably functions to calculate the position of the receiver (and/or external system) based on the fixed state vector (e.g., the fixed ambiguity vector).
- S 400 is preferably performed by the computing system; however, S 400 can additionally and/or alternatively be performed by the receiver, a reference station, and/or by any suitable component.
- the position of the receiver position is preferably determined from the cross-validated fixed carrier ambiguities; however, the receiver position can additionally or alternatively be determined using the first carrier phase ambiguity set (or a subset thereof), the second carrier phase ambiguity set (or a subset thereof), combinations of the above carrier phase ambiguities, and/or using any suitable states and/or data.
- the union of the state vectors (e.g., the validated states from each of the state vectors, the augmented carrier phase ambiguities, etc.) can be used to determine the receiver position.
- any power set of the state vectors e.g., validated in S 400 , fixed or otherwise determined in S 300
- any suitable states and/or state vectors
- the baseline vector b corresponding to the position of the receiver relative to a reference station, can be determined.
- the baseline vector can be determined based on the fixed integer ambiguity.
- the baseline vector can be determined in any suitable manner.
- the position of the receiver e.g., relative to the reference station(s), absolute position
- S 500 can determine (e.g., calculate) the receiver position in any suitable manner.
- S 500 can include determining an external system (e.g., vehicle) position and/or location based on the receiver position.
- the external system position can be determined based on a transformation from a receiver coordinate to an external system coordinate, based on a known relative position of the receiver to the external system, and/or otherwise be determined.
- the external system position can be a relative position (e.g., relative to the receiver, relative to a reference location, etc.) and/or an absolute position (e.g., coordinates).
- S 500 can include generating instructions for operation of the external system based on the external system location.
- operation instructions for a vehicle can include making lane adjustments to remain within a lane of traffic, turning instructions, and/or other instructions for traversing a path, obeying laws, and/or other purposes.
- the method can include receiving one or more sets of satellite observations, correcting each set of the satellite observations, determining a state vector associated with each set of satellite observations, fixing the carrier phase ambiguities for each state vector, and comparing the phase ambiguities associated with each of the state vectors.
- the phase ambiguities can be cross-validated.
- the cross-validated carrier phase ambiguities e.g., cross-validated fixed carrier phase ambiguities
- the method can include any suitable steps.
- comparing the state vectors can use two of the state vectors (e.g., from the set of N state vectors) to determine an intersection transformation to bring the set of state vectors to an intersection space.
- a transformation T 1 can be associated with the first state vector and a transformation T 2 can be associated with the second state vector can be determined.
- T 1 and T 2 are preferably different; however, T 1 and T 2 can be the same.
- the transformations can be used to convert the state vectors (e.g., ambiguity vectors of the state vectors) to a transformed spaced (e.g., transformed state vectors).
- the transformations can be used to determine an intersection space and associated intersection transformation that transforms the state vectors (or transformed state vectors) to the intersection space.
- the intersection space can correspond to a portion of (up to and including) all of the transformed space.
- the intersection space is preferably nontrivial (e.g., contains elements other than 0 and/or null values).
- a new transformation e.g., V 1 and V 2
- Both of the transformed state vectors can be transformed into the intersection space, for example to generate intersection space state vectors.
- the intersection space state vectors can be compared for example based on a validation criterion.
- intersection space transformation can be applied to the other state vectors.
- the state vectors can be nullified and/or a different pair of state vectors can be compared.
- comparing the state vectors can be performed in any suitable manner.
- the methods of the preferred embodiment and variations thereof can be embodied and/or implemented at least in part as a machine configured to receive a computer-readable medium storing computer-readable instructions.
- the instructions are preferably executed by computer-executable components integrated with a system for GNSS PVT generation.
- the computer-readable medium can be stored on any suitable computer-readable media such as RAMs, ROMs, flash memory, EEPROMs, optical devices (CD or DVD), hard drives, floppy drives, or any suitable device.
- the computer-executable component is preferably a general or application specific processor, but any suitable dedicated hardware or hardware/firmware combination device can alternatively or additionally execute the instructions.
- Embodiments of the system and/or method can include every combination and permutation of the various system components and the various method processes, wherein one or more instances of the method and/or processes described herein can be performed asynchronously (e.g., sequentially), concurrently (e.g., in parallel), or in any other suitable order by and/or using one or more instances of the systems, elements, and/or entities described herein.
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
A system and method for determining a mobile receiver position includes receiving a first and a second set of satellite observations; determining a first and second fixed ambiguity set associated with a first and second transformation respectively; determining cross-validated ambiguities between the first and second fixed ambiguity sets; and determining the mobile receiver position based on at least one of the first or second fixed ambiguity sets.
Description
- This application claims the benefit of U.S. Provisional Application No. 62/946,576, filed 11 Dec. 2019, which is incorporated in its entirety by this reference.
- This invention relates generally to the satellite positioning field, and more specifically to a new and useful system and method in the satellite positioning field.
- Traditionally, a receiver position is determined based on a set of satellite observations within a given time frame. These satellite observations are used to calculate the receiver position. However, satellite observations can, for example, include outliers and noisy data that can decrease the positioning accuracy. To achieve high levels of precision in the position, it can be desirable to validate the result. Thus, there is a need in the satellite positioning field to create a new and useful system and method. This invention provides such new and useful system and method for validating GNSS ambiguities.
-
FIG. 1 is a schematic representation of the system. -
FIG. 2 is a schematic representation of the method. -
FIG. 3 is a schematic representation of an example of transforming the state vectors. -
FIG. 4A is a block diagram representation of an example of the method. -
FIG. 4B is a block diagram representation of an example of validating the state vectors. - The following description of the preferred embodiments of the invention is not intended to limit the invention to these preferred embodiments, but rather to enable any person skilled in the art to make and use this invention.
- As shown in
FIG. 2 , a method includes receiving satellite observations, determining state vectors, validating the state vectors, and determining the receiver position. The method can optionally include correcting the satellite observations. Validating the state vectors can include determining one or more transformations, transforming the ambiguities, and comparing the ambiguities. - As shown in
FIG. 1 , the system includes one or more receivers and a computing system. The system can optionally include one or more sensors, and/or any suitable components. The system is preferably communicably connected to one or more data sources (e.g., reference stations, satellites, satellite constellations, etc.). - The system and method preferably function to estimate a position of a GNSS receiver and/or an external system. Embodiments of the system and/or method can be used, for example, in autonomous or semi-autonomous vehicle guidance (e.g., for unmanned aerial vehicles (UAVs), unmanned aerial systems (UAS), self-driving cars, agricultural equipment, robotics, rail transport/transit systems, autonomous trucking, last mile delivery, etc.), GPS/GNSS research, surveying systems, user devices, mobile applications, internet-of-things (TOT) devices, and/or may be used in any other suitable application. In specific examples, the system (and/or components) can be coupled to any suitable external system such as a vehicle (e.g., UAV, UAS, car, truck, etc.), robot, railcar, user device (e.g., cell phone), and/or any suitable system, and can provide positioning data, integrity data (e.g., protection level data), and/or other data to said system.
- Variations of the technology can confer several benefits and/or advantages.
- First, the inventors have discovered that in some applications, to achieve high accuracy and/or high integrity positioning solutions based on GNSS techniques, in addition to validation steps performed when fixing the integer ambiguities, validation steps between integer ambiguity sets (e.g., fixed integer ambiguities) can be beneficial. The inter-ambiguity set validation can increase a confidence in the resulting ambiguities and can lead to more accurate and/or higher integrity positioning solutions.
- Second, variants of the technology can be used independent of the data sources (e.g., satellites, satellite constellations, satellite frequencies, sensors, reference stations, etc.) and/or relationships between or within data collected from data sources. In specific examples, the technology can validate the state vectors (e.g., state vectors that are determined for different data sets and/or by different filters), regardless of what transformation(s) are applied to each state vector.
- Third, variants of the technology can determine an optimal or near-optimal number of states that can be compared between state vectors, which can improve the probability of successful validation between the state vectors, can provide an optimal or near-optimal set of validated states, and/or otherwise improve the state vector determination.
- However, variants of the technology can confer any other suitable benefits and/or advantages.
- The system preferably functions to determine the position of a receiver and/or external system, determine a plurality of state vectors, and validate the state vectors.
- The system preferably uses a set of data collected by one or more data sources. Data sources can include: receivers, sensors (e.g., located onboard the receiver, the external system, a reference stations, etc.), databases, reference stations, satellites, and/or any other suitable data source. Examples of data that can be used include: satellite observations, sensor observations, and/or any other suitable data.
- The receiver (e.g., mobile receiver) preferably functions to receive a set of satellite observations (e.g., satellite signals) from one or more satellites. In variants, the receiver can determine the location (e.g., pseudorange) of the receiver (and/or external system) based on the satellite observations. The receiver is preferably in communication with the computing system. However, the receiver can be integrated with the computing system, and/or the receiver and computing system can be arranged in any suitable manner. The receiver is preferably a stand-alone device (e.g., a GNSS receiver, antenna). However, the receiver can be integrated into the external system (e.g., be a component of an automobile, aero vehicle, nautical vehicle, etc.), can be a user device (e.g., smart phone, laptop, cell phone, smart watch, etc.), and/or can be configured in any suitable manner.
- The set of satellite observations can include orbital data (e.g., ephemeris), timestamp, code data, carrier phase data, pseudocode data, and/or any suitable data. The set of satellite observations preferably includes satellite observations corresponding to satellites from more than one satellite constellation (e.g., Global Positioning System (GPS), GLObal Navigation Satellite System (GLONASS), BeiDou positioning System (BDS), Galileo, etc.). However, the set of satellite observations can correspond to satellites from a single satellite constellation, can include data from an augmentation system (e.g., Satellite Based Augmentation System (SBAS) such as Wide Area Augmentation System (WAAS), European Geostationary Navigation Overlay Service (EGNOS), Multi-Functional Satellite Augmentation System (MSAS), Omnistar, StarFire, etc.; Ground Based Augmentation Systems (GBAS) such as Local Area Augmentation System (LAAS); etc.), and/or can include any suitable data. Each satellite observation from the set of satellite observations preferably corresponds to a common time window (e.g., epoch). However, each satellite observation can be associated with a time stamp (e.g., time of transmission, time of receipt, time of processing, etc.), and/or the satellite observations can have any suitable timing.
- In variants of the system including more than one receiver, each receiver can be configured to receive satellite observations associated a satellite constellation, a carrier frequency (e.g., the L1, L2, L5, E1, E5a, E5b, Eab, E6, G1, G3, B1, B2, B3, LEX, etc. frequencies), and/or corresponding to any suitable data. In an illustrative example, a receiver can include a first antennae configured to receive satellite signals from GPS satellites and a second antennae configured to receive satellite signals from Galileo satellites.
- The computing system preferably functions to process the data (e.g., satellite observations, reference station observations, sensor data, etc.) received by the receiver. The computing system can: aggregate the data (e.g., combine the receiver satellite observations, reference station satellite observations, satellite corrections, and/or sensor data; reorganize the receiver satellite observations, reference station satellite observations, and/or sensor data such as based on the time stamp, time of transmission, time of receipt, etc.; etc.), filter the data (e.g., to calculate state vectors, ambiguities such as phase ambiguities, etc. associated with the data), calculate the receiver position (e.g., based on the ambiguities), calculate the external system location, correct the data (e.g., correct the satellite observations for clock errors, hardware bias, atmospheric effects, etc.), and/or can process the data in any suitable manner. The computing system can be local (e.g., to the external system, to the receiver, to the sensor, etc.), remote (e.g., cloud computing, server, networked, etc.), and/or otherwise distributed.
- The computing system is preferably communicably coupled to the receiver and/or the sensors, but can be communicably coupled to any suitable data sources. The computing system is preferably colocalized with (e.g., integrated into) the receiver (and/or external system), but the computing system can be remote from the receiver (and/or external system), and/or configured in any suitable manner.
- In specific variants of the system including one or more sensors, the sensors preferably function to measure sensor data (e.g., auxiliary data) associated with the external system (and/or the receiver). The sensor data can be used to assist (e.g., speed-up, correct, refine, etc.) the calculation (e.g., calculating the state vector, estimating the phase ambiguity). In a specific example, the sensor data can be used to estimate the receiver position at times when satellite observations are not received (e.g., in an urban canyon; due to an obstruction such as a billboard, weather, etc.; etc.). However, the sensor data can be used at any suitable time and in any suitable manner. The sensors are preferably in communication with the computing system. The sensors can be: on-board the external system, on-board a separate external system, integrated into the receiver, separate from the receiver, and/or otherwise associated with the receiver. The sensor data can include: inertial data (e.g., velocity, acceleration), odometry, pose (e.g., position, orientation), mapping data (e.g., images, point clouds), temperature, pressure, ambient light, and/or any other suitable data. The sensors can include one or more of: inertial moment unit (IMU), accelerometer, gyroscope, magnetometer, odometer (e.g., visual odometry, wheel odometry, etc.), and/or any suitable sensor.
- The method functions to determine the position of a receiver (e.g., an external system coupled to the receiver) based on one or more sets of satellite observations. The method preferably validates the position of the receiver and/or one or more intermediate (e.g., state vector) used to determine the position of the receiver. One or more instance of the method and/or steps of the method can be performed in series and/or in parallel (e.g., concurrently).
- Receiving one or more satellite observations S100 preferably functions to measure (e.g., at the receiver, at reference stations, etc.) and/or access one or more set of satellite observations (e.g., carrier phase measurements. pseudo-range measurements, etc.) from one or more observed satellites. The satellite observations can be measured and/or received by a receiver, retrieved from a database (e.g., retrieve stored satellite observations; retrieve stored corrections; retrieve an almanac such as of weather conditions, tides, etc.; etc.), and/or be otherwise received. S100 can include receiving Doppler measurement data, sensor data, and/or any suitable data. The satellite observations can include signals from one or more satellite constellations. Each set of satellite observations preferably corresponds to a plurality of satellite constellations. However, one or more set of satellite observations can correspond to a single satellite constellation.
- S100 is preferably performed before correcting the satellite observations S200; however, S100 can be performed at the same time as S200. S100 is preferably performed before determining the state vectors S300; however, S100 and S300 can be performed at the same time and/or S300 can be performed before S100 (e.g., during different time windows). S100 is preferably performed by a receiver. However, S100 can additionally or alternatively be performed by one or more reference stations, by one or more sensors, and/or by any suitable component.
- The data (e.g., satellite observations, sensor data, etc.) received in S100 is preferably measured during a time window (e.g., time period, epoch). Each of the sets of satellite observations can correspond to the same and/or to a different time window. In variants, after the time window has expired, the method can include repeating S100 during a second time window. The second time window can be the same duration as the time window and/or a different duration. The second time window preferably immediately follows the time window; however, the second time window can be delayed relative to the time window by any amount. The (first and second) time windows are preferably separate and distinct (e.g., consecutive, nonconsecutive, etc.), but can alternatively overlap. The second time window can additionally or alternatively be started in response to a trigger event (e.g., convergence of one or more calculation; collective convergence of the calculations; receiver position; loss of satellite observations such as for a specific satellite, for an amount of time, etc.; etc.), at a predetermined time, and/or at any suitable time. However, the time window can be extended and/or any suitable action can occur in response to the trigger event.
- Correcting the satellite observations S200 preferably functions to correct the satellite observations for local and/or global errors. S200 is preferably performed by a computing system (e.g., a correction module of a computing system), but can be performed by any system or component. The corrections are preferably received from a computing system (e.g., a remote computing system, a cloud computing system, corrections service, etc.), but can additionally or alternatively be generated by the computing system (e.g., based on reference station data, based on almanac data, based on a model, etc.), received from a reference station, and/or be otherwise accessed. S200 is preferably performed before determining the state vector S300. However, S200 can be performed at the same time as and/or after S300. S200 preferably corrects all of the satellite observations within each set of satellite observations. However, any suitable subset of satellite observations can be corrected. Corrections can be applied in response to (e.g., to correct): hardware bias, timing errors (e.g., satellite clock, receiver clock, reference station clock, etc.), atmospheric effects (e.g., ionosphere effects, tropospheric effects), relativistic effects, local effects (e.g., multipath errors), global effects, noise (e.g., outliers), and/or any suitable effects. In a specific example, reference station satellite observations can be used to correct the receiver satellite observations such as by calculating double differenced satellite observations. In a second example, effects (e.g., atmospheric effects) on satellite observations can be modelled. Based on the model, the satellite observations can be corrected. However, correcting the satellite observations can be performed in any suitable manner. In an illustrative example, the corrections can be generated as disclosed in and/or received from a corrections service as disclosed in U.S. patent application Ser. No. 16/589,932 filed 1 Oct. 2019 entitled “SYSTEMS AND METHODS FOR DISTRIBUTED DENSE NETWORK PROCESSING OF SATELLITE POSITIONING DATA,” which is incorporated in its entirety by this reference. However, the corrections can be otherwise generated.
- Determining the state vector S300 preferably functions to process the data (e.g., satellite observations, corrected satellite observations, sensor data) to determine phase ambiguities (e.g., fixed integer phase ambiguities) associated with each set of satellite observations. States of the state vector preferably include a carrier phase ambiguity (e.g., calculated as float ambiguities, calculated as an integer ambiguity, etc.) associated with one or more satellites (e.g., each satellite that provided one or more observation in the set of satellite observations) and phase ambiguity covariances (and/or correlations) associated with each one or more pairs of satellites (e.g., pairs of satellites that provided observations in the set of satellite observations). However, the state vector can additionally or alternatively include residuals (e.g., for the phase ambiguities, for the covariances), a validity of one or more states, a receiver position, an external system location, pseudorange, atmospheric error or effect, estimated error (e.g., estimated error in one or more states in the state vector), derivatives of one or more states (e.g., first, second, etc. derivative such as with respect to time, with respect to position, etc.), and/or any suitable state. S300 is preferably performed by a computing system (e.g., a processing module of a computing system), but can be performed by any component. S300 is preferably performed before validating the state vectors S400. However, S300 can be performed at the same time as and/or after S400 (e.g., when a subset of the state vector is validated and said subset is used to determine the remainder of the state vector). S300 is preferably performed twice for each instance of the method (e.g., to generate a first state vector and a second vector associated with a first and second set of satellite observations, respectively). However, S300 can be performed once for each instance of the method (e.g., wherein one instance of the method is compared or validated relative to another instance of the method), and/or can be performed any number of times.
- S300 can include: determining a set of floating phase ambiguity hypotheses, determining a set of integer phase ambiguity hypotheses from the set of floating phase ambiguity hypotheses, performing hypothesis testing to validate the set of integer phase ambiguities hypotheses, and/or any suitable steps. In an illustrative example, S300 can include steps and/or be performed as disclosed in U.S. patent application Ser. No. 16/685,927 filed 15 Nov. 2019 entitled “SYSTEM AND METHOD FOR SATELLITE POSITIONING” and/or U.S. patent application Ser. No. 16/817,196 filed 12 Mar. 2020 entitled “SYSTEMS AND METHODS FOR REAL TIME KINEMATIC SATELLITE POSITIONING,” each of which is incorporated in its entirety by this reference.
- Determining a set of floating phase ambiguity hypotheses preferably functions to determine a phase ambiguity (typically in a floating point representation, but other representations can be used) associated with each satellite represented in the set of satellite observations. The set of floating phase ambiguity hypotheses are preferably determined using a filter. The filter is preferably a Kalman filter. However, the filter can be an extended Kalman filter, an unscented Kalman filter, a Bierman-Thornton filter, a particle filter, Monte Carlo simulation, and/or any suitable mean-square error filter and/or sensor fusion algorithm. The inputs to the filter can include: satellite observations (e.g., receiver satellite observations, corrected satellite observations, reference station satellite observations, etc.), sensor data, and/or any suitable data. The output of the filter is preferably a state vector and/or one or more states of a state vector, but can include any suitable information.
- Determining a set of integer phase ambiguity hypotheses preferably functions to constrain (and/or otherwise determine from) the set of floating phase ambiguity hypotheses to a set of integer phase ambiguity hypotheses. Integer phase ambiguities can be beneficial for improving an accuracy of the determined receiver position as compared to the receiver position determined using floating ambiguities. Determining the set of integer phase ambiguity hypotheses can include reducing a correlation between ambiguities of the set of phase ambiguities (e.g., floating phase ambiguity hypotheses) and performing a search (e.g., a least squares search, an absolute value search, etc.) to identify possible integer phase ambiguities to include in the set of integer phase ambiguity hypotheses. For example, Least-Squares Ambiguity Decorrelation Adjustment (LAMBDA) algorithm, modified-LAMBDA (MLAMBDA) algorithm, LLL reduction algorithm, a whitening transformation, a coloring transformation, a decorrelation transformation, rounding, the Ambiguity Function Method (AFM), Fast Ambiguity Resolution Approach (FARA), Least-Squares Ambiguity Search Techniques (LSAST), integer bootstrapping, and/or any suitable decorrelation or reduction algorithm can be used to fix and/or identify integer phase ambiguities. The set of integer phase ambiguity hypotheses preferably includes integer phase ambiguities that satisfy a hypothesis criterion, but can include any suitable integer phase ambiguities. Example of hypotheses criterion include a threshold sum of least squares, a threshold sum of absolute differences, a threshold likelihood of the correct solution, and/or any suitable criterion. However, the set of integer phase ambiguity hypotheses can be generated in any manner.
- Performing hypothesis testing preferably functions to determining the integer phase ambiguities from the set of integer phase ambiguity hypotheses that are most likely to be correct. The hypothesis test is preferably a Bayesian inference, but can additionally or alternatively be based on a statistical confidence, significance testing, and/or any suitable hypothesis testing. Examples of hypothesis tests can include difference test, ratio test, projector test, f-test GIA test, and/or any suitable hypothesis test. In an illustrative example, performing hypothesis testing can include determining a probability (e.g., a likelihood, a log likelihood, etc.) associated with subsets of integer phase ambiguities of the integer phase ambiguity hypotheses; computing a ratio of the probability between two subsets of integer phase ambiguities (e.g., within the same ambiguity set, across different ambiguity sets); when the ratio between a most-likely and a next most likely subset of integer phase ambiguities exceeds a threshold, storing the most-likely subset of integer phase ambiguities as the integer phase ambiguities and ceasing hypothesis testing. However, any hypothesis testing can be performed.
- In variants, the most-likely subset of integer phase ambiguities can be referred to as internally-validated integer phase ambiguities (e.g., integer phase ambiguities that are validated against the remaining integer phase ambiguities from the same set of integer phase ambiguity hypotheses). The internally-validated integer phase ambiguities can be contrasted with cross-validated integer phase ambiguities (e.g., integer phase ambiguities as validated in step S400).
- In embodiments of S300, S300 can be performed more than once to generate at least two state vectors. Each state vector can be determined consecutively, iteratively, sequentially, delayed, and/or with any suitable timing relative to the other state vectors. The plurality of state vectors is preferably independent of one another, but can be correlated, dependent on each other, and/or otherwise be related. Each state vector can include a subset of common states (e.g., common integer phase ambiguities such as those associated with the same data source), a subset of distinct states (e.g., integer phase ambiguities associated with different data sources, uncommon states, etc.), and/or the states in the first and second state vectors can be otherwise related. Different state vectors can differ in: satellite composition, frequency, sampling epoch, otherwise differ, or not differ (e.g., share all or a portion of satellite observations). The states within each state vector can share: satellites, frequencies, sampling epochs, other parameters, or no parameters.
- In a first variant, two (or more) state vectors can be associated with (e.g., generated from) a two (or more) distinct sets of satellite observations (e.g., a first and a second set of satellite observations). The distinct sets of satellite observations preferably include observations from at least three common satellites, but can include observations from any number of common satellites. In an illustrative example, the first and second set of satellite observations can correspond to sets of satellite observations that are sent to distinct Kalman filters as disclosed in U.S. patent application Ser. No. 16/685,927 filed 15 Nov. 2019 entitled “SYSTEM AND METHOD FOR SATELLITE POSITIONING” and incorporated in its entirety by this reference. In a second variant, two (or more) state vectors can be generated using two (or more) different filters (and/or methods). For example, a first state vector can be generated using phase hypotheses generated using a Kalman filter and a second state vector can be generated using phase hypotheses generated using a Bierman-Thornton filter. In a second example of the second variant, a first and second state vector can each be generated with different initial state guesses. However, the state vectors can be generated using any suitable data and/or any suitable method.
- In embodiments of S300, the state vector(s) can be associated with one or more transformations. The transformation can improve the phase ambiguity determination (e.g., enable the phase ambiguity to be determined quickly, enable the phase ambiguity to be determined with greater accuracy, etc.), improve the receiver position estimate (e.g., improve the accuracy or integrity of the determined receiver position), and/or otherwise be used. The output of the filter (e.g., the state vector) can include the transformed phase ambiguity, the transformed phase ambiguity can be determined from the state vector (e.g., by transforming the phase ambiguities of the state vector), and/or the transformed phase ambiguity can be otherwise related to the state vector. Examples of transformations include: single differenced combinations, double differenced combinations, triple differenced combinations, wide lane combinations, middle lane combinations, narrow lane combinations, extra-wide lane combinations, ionosphere free combinations, geometry free combinations, noise reduced combinations, ionospheric delay reduced combinations, noise reducing combinations, linear combinations (e.g., between two or more carrier frequencies of satellite observations, between two or more data sources such as satellite observations and reference station observations, etc.), and/or any suitable combination of data source observations.
- Each transformation is preferably an n×m matrix where n is less than or equal to m, where m corresponds to the number of ambiguities in the respective ambiguity set (e.g., ambiguity vector, state vector, etc.), and where n corresponds to the number of ambiguities in the transformed space. However, n can be greater than m. The elements of the transformation matrices are preferably integers, but can be rational, real, complex, and/or have any values. In a set of specific examples, the transformations can include: transforming to wide-lane ambiguities, transforming to decrease the correlation between the ambiguities, transforming to a rotated state space, transforming to a scaled state space, transforming to a truncated state space, and/or transforming in any suitable manner.
- In an illustrative example, the transformed phase ambiguity (e.g., transformed integer phase ambiguity) can be determined according to: {circumflex over (z)}i={circumflex over (T)}iâi where zi is the transformed phase ambiguity, Ti is the transformation relating the underlying phase ambiguity ai to the transformed phase ambiguity. The transformed phase ambiguity and the underlying phase ambiguity are preferably both integers, but either can be rational numbers, floating values, real valued, and/or otherwise by related. However, the transformed phase ambiguity can be otherwise determined. When more than one state vector is determined, each state vector can be associated with the same or a distinct transformation.
- Validating the state vectors S400 preferably functions to cross validate one or more states from two or more state vectors (e.g., determine what states of two or more state vectors match). S400 preferably functions to cross-validate a subset of integer phase ambiguities (e.g., internally-validated integer phase ambiguities) from the first state vector with a subset of integer phase ambiguities from the second state vector. However, any suitable state(s) can be cross-validated between the first and second state vector. S400 is preferably performed by a computing system (e.g., a validation module of a computing system). However, S400 can be performed by any suitable component. S400 is preferably performed before determining the receiver position S500. However, S400 can be performed at the same time as and/or after (e.g., in examples where the receiver position is used to validate the state vectors) S500. S400 is preferably performed after S300, such that the states (e.g., the phase ambiguities) that S400 is analyzing are fixed; however, S400 can optionally analyze unfixed states (e.g., unfixed ambiguities). The state vectors are preferably validated based on a collective convergence criterion (e.g., a validation metric, validation score). The collective convergence criterion can include: one or more subsets of transformations associated with the state vectors matching, subsets of ambiguities matching across the state vectors, a residual threshold, one or more states (e.g., from each of the calculations) matching, a timeout, a covariance threshold, a receiver position threshold, and/or any suitable criterion.
- S400 can include: combining the first and second state vectors 420, determining an intersection transformation S450, comparing the state vectors in the intersection space S480, and/or any suitable steps.
- The first and second state vectors preferably at least partially share states from satellite observations from common: satellites, frequencies, and/or epochs, but can alternatively include entirely unrelated satellite observations. In a first example, the first and second state vectors share states from observations from common satellites, frequencies, and epochs. In a second example, the first and second state vectors share states from observations from common satellites and frequencies, but different epochs.
- S420 preferably functions to combine the first and second state vector to generate a combined state vector. The combined state vector preferably includes the non-transformed integer phase ambiguities (e.g., internally-validated integer phase ambiguity, fixed integer phase ambiguities, etc.) from the first and second state vector, but can additionally or alternatively include the transformed integer phase ambiguities, a subset of integer phase ambiguities, and/or any suitable states from the first and/or second state vectors. The combined state vector can be an augmented state vector (e.g., including all or a majority of integer phase ambiguities associated with the first and the second state vector), a common state vector (e.g., including the common integer phase ambiguities from the first and second state vector), an uncommon state vector (e.g., including uncommon ambiguities or unshared ambiguities), and/or any suitable state(s) can be included in the combined state vector.
- In a specific example, as shown in
FIG. 3 , S420 can be performed using two fixed integer phase ambiguity sets. The fixed integer phase ambiguity sets are preferably not linearly independent (e.g., include one or more common fixed integer phase ambiguity, two fixed integer phase ambiguity sets, a1 and a2 can be expressed as -
- where ac are common ambiguities and a1′ and a2′ are unique ambiguities associated with a1 and a2 respectively, etc.). However, the fixed integer phase ambiguity sets can be linearly independent. The common ambiguities (e.g., shared ambiguities) preferably correspond to ambiguities associated with a common data source (e.g., satellite, reference station, satellite constellation, satellite frequency, sensor, etc.). The unique ambiguities (e.g., in the set of unique ambiguities associated with a1 and/or a2) preferably correspond to ambiguities associated with data sources that are represented in a1 or a2 but not both. However, the shared and/or unique ambiguities can correspond to any suitable ambiguities. Similarly, each fixed integer phase ambiguity set can be associated with a transformation, T1 and T2 associated with the first and second integer phase ambiguity set, respectively. The transformations can be written as
-
- In this specific example, the transformed ambiguities z1 and z2 can be determined from z1=T1a1 and z2=T2 a2. However, the transformed ambiguities can have any suitable form.
- In variants that include an augmented integer ambiguity set, the augmented integer ambiguity set (e.g., augmented integer ambiguity space, augmented state vector, etc.) preferably functions to ensure that the transformed states have the same dimensions. In an illustrative example, the augmented integer ambiguity set can be written as
-
- In this illustrative example, the augmented integer ambiguity set can be associated with a first and second augmented transformation that transforms the augmented integer ambiguity set to the first and second integer ambiguity sets. In this illustrative example, the augmented transformations can be written as:
-
- In this specific example, the augmented transformations {circumflex over (T)}1 and {circumflex over (T)}2 are preferably re-ordered to align with the states of the augmented ambiguity vector. However, the augmented integer ambiguity set can include any suitable states and/or be arranged in any manner.
- Determining an intersection transformation S450 functions to determine a transformation that transforms the state vectors to a shared intersection space. The intersection space preferably functions to enable the ambiguities to be compared in the same basis while retaining as many linearly independent rows or ambiguities as possible. The intersection space is preferably determined using the augmented transformations (augmented ambiguity transformations). However, the intersection space can be determined from the transformed ambiguities (e.g., zi), the initial ambiguities (e.g., ai), the initial transformations (e.g., Ti), the augmented integer phase ambiguities (e.g., a), and/or any suitable vector and/or transformation can be used.
- Determining the intersection transformation (e.g., shared space) preferably determines the row[{circumflex over (T)}∩]=row[{circumflex over (T)}1]∩row[{circumflex over (T)}2], where row[T∩] is the shared space (e.g., the space of transformed ambiguities that is shared between {circumflex over (T)}1 and {circumflex over (T)}2), and ∩ represents the intersection. The row space of Ti, row[{circumflex over (T)}i], is preferably the set of vectors that result from the possible combinations of rows of matrix Ti. The row space of the intersection of T1 and T2, row[{circumflex over (T)}∩], can then be determined from (e.g., is the same as) the intersection of the row space of T1 with the row space of T2, where the intersection of T1 and T2 can have the minimum number of rows that is possible while the row space condition is adhered to. However, the intersection transformation can be otherwise defined and/or determined.
- The intersection space can be determined, for example, by concatenating the transformations into a concatenated matrix (e.g., S), and searching for the null space of S. In a specific example, S=[{circumflex over (T)}1 T−{circumflex over (T)}2 T], wherein | can represent a concatenation operation. In this example, the elements of S are preferably integers. However, the elements of S can be natural, rational, real, and/or complex valued.
- In a specific example, the null space of S (e.g., null[S]) includes the set of vectors v=[V1 V2], where V1{circumflex over (T)}1−V2{circumflex over (T)}2=0. In this example, the null space of S can be partitioned into parts such as null [S]=[V1 V2]. In this example, the parts of the null space can be used to define the intersection space row [{circumflex over (T)}∩]=V1{circumflex over (T)}1=V2{circumflex over (T)}2.
- Each of the state vectors can be transformed to the intersection space, a subset of the state vectors (e.g., a1 and a2, z1 and z2, etc.) can be transformed to the intersection space, and/or any suitable states can be transformed to the intersection space. In a specific example, transforming the ambiguity vectors can include transforming the first set of fixed carrier ambiguities using the intersection transformation {circumflex over (T)}∩=V1{circumflex over (T)}1 so that the first set of carrier ambiguities are in the intersection space. Similarly, the second set of fixed carrier ambiguities using the intersection transformation {circumflex over (T)}∩=V2 {circumflex over (T)}2. However, the ambiguity vectors can be transformed to any suitable space using any suitable transformation.
- Comparing the state vectors in the intersection space S480 functions to determine whether the state vectors (or subsets thereof) match. The state vectors are preferably compared in the intersection space, but can additionally or alternatively be compared in the transformed space (e.g., comparing z1 and z2), in the untransformed space (e.g., comparing a1 and a2), in a different comparison space, and/or in any suitable space. S480 preferably compares the integer phase ambiguities associated with the same data source between the two state vectors. However, any suitable states can be compared. In an illustrative example, each row of the first state vector can be compared to the row of the second state vector with the same row index to identify matches between the respective integer phase ambiguities. In a second illustrative example, common ambiguities can be identified in the transformed space (e.g., by identifying the ambiguities determined from the same satellite observations, such as by using a unique satellite observation identifier) and compared.
- The state vectors (and/or subsets thereof) preferably match (e.g., are cross-validated) when a validation criterion is satisfied. The validation criterion (e.g., collective convergence criterion) preferably includes the number of differing ambiguities (e.g., between one or more integer phase ambiguity in the intersection space, transformed space, untransformed space, etc.) being less than a threshold, but can additionally or alternatively include a residual threshold, one or more states (e.g., from each of the state vectors) matching (e.g., having an identical value such as in the intersection space, the transformation space, the untransformed space, etc.), a difference between two ambiguities (e.g., a difference threshold), a timeout, a quorum of matching transformed ambiguities, and/or any suitable criterion. However, the ambiguities can be compared in any suitable manner. The validation criterion can be based on (e.g., depend on) the number of ambiguities (e.g., a threshold number of ambiguities), depend on the data source(s) (e.g., number of satellites, number of satellite constellations, identity of the satellite constellation(s), etc.), and/or depend on any suitable properties.
- When the validation criterion is satisfied, the first and/or second state vectors can be validated (e.g., cross-validated; considered validated or cross-validated), a subset of the state vectors (e.g., the state(s) that meet the validation criterion, common states, uncommon states, etc.) can be validated (e.g., cross-validated), and/or any suitable state(s) can be validated, where validated ambiguities can subsequently be used to correct the carrier measurements and/or otherwise be used to determine the receive and/or external system location. When the validation criterion is not satisfied, a new intersection space can be determined, a new set of state vectors can be compared, a flag can be issued that the state vectors cannot be cross-validated, transmitting an error message, calculating the receiver position (e.g., with non-validated ambiguities), starting a new time window, extending the time window, adding states (e.g., ambiguities) to the ambiguity vector(s), nullifying the ambiguity (and/or state) vectors, the external system can be operated at a lower accuracy and/or integrity level location determined, additional states can be determined (e.g., by performing S300), the method can be reinitiated, and/or any suitable response can be performed.
- In a specific example, comparing the state vectors includes comparing the first and second fixed integer ambiguity set by checking whether V1{circumflex over (z)}1==V2{circumflex over (z)}2 is valid for each phase ambiguity. In this example, when the equality is valid (e.g., when V1{circumflex over (z)}1 matches V2{circumflex over (z)}2), a new transformation ({circumflex over (T)}∩=V1{circumflex over (T)}1) and corresponding ambiguities (e.g., fixed ambiguities, {circumflex over (z)}∩ 1=V1{circumflex over (z)}1) can be computed. However, when the equality is not valid (e.g., when V1{circumflex over (z)}1 does not match V2{circumflex over (z)}2), S480 can include determining transformations using a different pair of ambiguity vectors, reporting an error, starting a new time window (e.g., restarting the method, returning to S100), extending the current time window, starting a new filter (e.g., determining a new state vector such as in S300), modifying a filter (e.g., change inputs to the filter), and/or can perform any suitable steps. In a variation of this example, matching preferably corresponds to a matching criterion (e.g., score, metric, quality). The matching criterion can be a threshold (e.g., a maximum difference between each value in V1{circumflex over (z)}1 and V2{circumflex over (z)}2 such as each row, each column, each element, etc.), a threshold distance, a residual (e.g., a target residual, a maximum residual, etc.), and/or any suitable value. In a second variation of this example, matching can require the two state vectors (e.g., corresponding ambiguities, ambiguities in the intersection space, transformed ambiguities, etc.) to be the same. However, any suitable matching criterion can be used.
- In specific variants, when a subset of the transformed ambiguity vectors match, that subset of ambiguities can be cross-validated (e.g., considered cross-validated) to generate a cross-validated fixed carrier ambiguity set. The cross-validated fixed carrier ambiguity set preferably includes at least 3 carrier phase ambiguities, but can contain fewer than 3 carrier phase ambiguities and/or any number of carrier phase ambiguities. In these variants, the subset of carrier ambiguities that do not match can be nullified (e.g., cleared from the state vectors), can be cross-validated (e.g., considered valid at the same time the remaining states are considered valid), undergo additional validation testing, and/or otherwise be processed. However, any suitable response can be taken in response to a subset of the ambiguity vectors not matching. However, the fixed ambiguity vector (e.g., validated ambiguity vector) can be determined in any suitable manner.
- Determining the receiver position S500 preferably functions to calculate the position of the receiver (and/or external system) based on the fixed state vector (e.g., the fixed ambiguity vector). S400 is preferably performed by the computing system; however, S400 can additionally and/or alternatively be performed by the receiver, a reference station, and/or by any suitable component. The position of the receiver position is preferably determined from the cross-validated fixed carrier ambiguities; however, the receiver position can additionally or alternatively be determined using the first carrier phase ambiguity set (or a subset thereof), the second carrier phase ambiguity set (or a subset thereof), combinations of the above carrier phase ambiguities, and/or using any suitable states and/or data. In a specific example, the union of the state vectors (e.g., the validated states from each of the state vectors, the augmented carrier phase ambiguities, etc.) can be used to determine the receiver position. However, any power set of the state vectors (e.g., validated in S400, fixed or otherwise determined in S300) and/or any suitable states (and/or state vectors) can be used.
- In specific variants including reference stations, the baseline vector b, corresponding to the position of the receiver relative to a reference station, can be determined. The baseline vector can be determined based on the fixed integer ambiguity. However, the baseline vector can be determined in any suitable manner. In these variants, the position of the receiver (e.g., relative to the reference station(s), absolute position) can be determined by applying b to the reference station location. However, S500 can determine (e.g., calculate) the receiver position in any suitable manner.
- S500 can include determining an external system (e.g., vehicle) position and/or location based on the receiver position. The external system position can be determined based on a transformation from a receiver coordinate to an external system coordinate, based on a known relative position of the receiver to the external system, and/or otherwise be determined. The external system position can be a relative position (e.g., relative to the receiver, relative to a reference location, etc.) and/or an absolute position (e.g., coordinates).
- S500 can include generating instructions for operation of the external system based on the external system location. For example, operation instructions for a vehicle can include making lane adjustments to remain within a lane of traffic, turning instructions, and/or other instructions for traversing a path, obeying laws, and/or other purposes.
- In a specific example, as shown in
FIG. 4A , the method can include receiving one or more sets of satellite observations, correcting each set of the satellite observations, determining a state vector associated with each set of satellite observations, fixing the carrier phase ambiguities for each state vector, and comparing the phase ambiguities associated with each of the state vectors. When a validation criterion is met, the phase ambiguities can be cross-validated. The cross-validated carrier phase ambiguities (e.g., cross-validated fixed carrier phase ambiguities) can be used to determine the receiver position. However, the method can include any suitable steps. - In a specific example, as shown in
FIG. 4B , comparing the state vectors can use two of the state vectors (e.g., from the set of N state vectors) to determine an intersection transformation to bring the set of state vectors to an intersection space. In this example, a transformation T1 can be associated with the first state vector and a transformation T2 can be associated with the second state vector can be determined. T1 and T2 are preferably different; however, T1 and T2 can be the same. The transformations can be used to convert the state vectors (e.g., ambiguity vectors of the state vectors) to a transformed spaced (e.g., transformed state vectors). The transformations (and/or transformed state vectors) can be used to determine an intersection space and associated intersection transformation that transforms the state vectors (or transformed state vectors) to the intersection space. The intersection space can correspond to a portion of (up to and including) all of the transformed space. The intersection space is preferably nontrivial (e.g., contains elements other than 0 and/or null values). After determining the intersection space, a new transformation (e.g., V1 and V2) can be determined for each transformed state vector. Both of the transformed state vectors can be transformed into the intersection space, for example to generate intersection space state vectors. The intersection space state vectors can be compared for example based on a validation criterion. When the intersection space state vectors match, the intersection space transformation can be applied to the other state vectors. In this specific example, when the intersection space state vectors do not match (and/or T∩ is a trivial transformation), the state vectors can be nullified and/or a different pair of state vectors can be compared. However, comparing the state vectors can be performed in any suitable manner. - The methods of the preferred embodiment and variations thereof can be embodied and/or implemented at least in part as a machine configured to receive a computer-readable medium storing computer-readable instructions. The instructions are preferably executed by computer-executable components integrated with a system for GNSS PVT generation. The computer-readable medium can be stored on any suitable computer-readable media such as RAMs, ROMs, flash memory, EEPROMs, optical devices (CD or DVD), hard drives, floppy drives, or any suitable device. The computer-executable component is preferably a general or application specific processor, but any suitable dedicated hardware or hardware/firmware combination device can alternatively or additionally execute the instructions.
- Embodiments of the system and/or method can include every combination and permutation of the various system components and the various method processes, wherein one or more instances of the method and/or processes described herein can be performed asynchronously (e.g., sequentially), concurrently (e.g., in parallel), or in any other suitable order by and/or using one or more instances of the systems, elements, and/or entities described herein.
- As a person skilled in the art will recognize from the previous detailed description and from the figures and claims, modifications and changes can be made to the preferred embodiments of the invention without departing from the scope of this invention defined in the following claims.
Claims (20)
1. A method for determining a mobile receiver position comprising:
receiving a first and a second set of satellite observations;
determining a first fixed ambiguity set associated with the first set of satellite observations and a second fixed ambiguity set associated with the second set of satellite observations, wherein the first fixed ambiguity set is associated with a first transformation and the second fixed ambiguity set is associated with a second transformation;
determining a first and a second intersection transformation from the first and second transformations, wherein the first and second intersection transformations transform the first and second fixed ambiguity sets to a shared intersection space;
determining cross-validated ambiguities from the first and second fixed ambiguity sets, wherein a cross-validated ambiguity is identified when a number of differing ambiguities from the first and second fixed ambiguity sets in the shared intersection space is less than a threshold; and
determining the mobile receiver position based on at least one of the first or second fixed ambiguity sets.
2. The method of claim 1 , further comprising:
generating an augmented ambiguity set from the first and second fixed ambiguity sets;
determining augmented first and second transformations based on the augmented ambiguity set and the first and second transformation; and
wherein the first and second intersection transformation are determined based on the augmented first and second transformations.
3. The method of claim 2 , wherein the first and second fixed ambiguity sets comprise common fixed ambiguities, wherein the augmented ambiguity set includes the common fixed ambiguities.
4. The method of claim 3 , wherein the augmented ambiguity set includes the uncommon ambiguities from the first and second fixed ambiguity sets.
5. The method of claim 1 , further comprising determining a union of ambiguities between the first and second fixed ambiguity sets, wherein the union of ambiguities is used to determine the mobile receiver position.
6. The method of claim 1 , further comprising determining a position of an external system based on the mobile receiver position.
7. The method of claim 6 , further comprising operating the external system based on to the external system position.
8. The method of claim 1 , wherein the first and second transformations are distinct.
9. The method of claim 1 , wherein determining each of the first and second fixed ambiguity set comprises:
determining a set of integer phase ambiguity hypotheses associated with the respective set of satellite observations; and
internally-validating the set of integer phase ambiguities.
10. The method of claim 9 , wherein internally-validating the set of integer phase ambiguities comprises:
performing hypothesis testing on the set of integer phase ambiguity hypotheses, wherein the hypothesis testing comprises determining a ratio between likelihoods of two integer phase ambiguity hypotheses; and
ceasing hypothesis testing when the ratio exceeds a threshold value.
11. The method of claim 9 , wherein determining the set of integer phase ambiguity hypotheses for each of the first and second ambiguity set comprises:
determining a set of floating integer phase ambiguities using a Kalman filter; and
determining the set of integer phase ambiguity hypotheses using at least one of a LAMBDA or MLAMBDA search algorithm.
12. The method of claim 1 , further comprising correcting the first and second set of satellite observations using a set of corrections received from a correction service, wherein the corrected first and second set of satellite observations are used to determine the first and second fixed ambiguity sets.
13. A system comprising:
a receiver comprising a processing module configured to:
receive a first and a second set of satellite observations;
determine a first fixed ambiguity set associated with the first set of satellite observations and a second fixed ambiguity set associated with the second set of satellite observations, wherein the first fixed ambiguity set is associated with a first transformation and the second fixed ambiguity set is associated with a second transformation;
generate an augmented ambiguity set from the first and second fixed ambiguity sets;
determine augmented first and second transformations based on the augmented ambiguity set and the first and second transformation;
determine a first and a second intersection transformation that transforms the first and second fixed ambiguity sets, respectively, to a shared intersection space;
determine cross-validated ambiguities from the first and second fixed ambiguity sets, wherein a cross-validated ambiguity is identified when a number of differing ambiguities from the first and second fixed ambiguity sets in the shared intersection space is less than a threshold; and
determine the receiver position based on at least one of the first or second fixed ambiguity sets.
14. The system of claim 13 , wherein the processing module is further configured to:
generate an augmented ambiguity set from the first and second fixed ambiguity sets;
determine augmented first and second transformations based on the augmented ambiguity set and the first and second transformation; and
wherein the first and second intersection transformation are determined based on the augmented first and second transformations.
15. The system of claim 14 , wherein the processing module is further configured to determine a union of ambiguities between the first and second fixed ambiguity sets, wherein the union of ambiguities is used to determine the receiver position.
16. The system of claim 13 , further comprising a second receiver, wherein the receiver is configured to receive satellite observations of the first and second set of satellite observations associated with a first satellite constellation and the second receiver is configured to receive satellite observations of the first and second set of satellite observations associated with a second satellite constellation.
17. The system of claim 13 , further comprising a vehicle wherein the receiver is mounted to the vehicle, wherein the processing module is further configured to determine a location of the vehicle based on the receiver position.
18. The system of claim 17 , wherein the processing module is further configured to generate instructions for operation of the vehicle based on the vehicle location.
19. The system of claim 13 , wherein determining each of the first and second fixed ambiguity set comprises:
determining a set of integer phase ambiguity hypotheses associated with the respective set of satellite observations;
performing hypothesis testing on the set of integer phase ambiguity hypotheses, wherein the hypothesis testing comprises determining a ratio between likelihoods of two integer phase ambiguity hypotheses; and
ceasing hypothesis testing when the ratio exceeds a threshold value.
20. The system of claim 19 , wherein determining a set of integer phase ambiguity hypotheses for each of the first and second ambiguity set comprises:
determining a set of floating integer phase ambiguities using a Kalman filter; and
determining the set of integer phase ambiguity hypotheses using at least one of a LAMBDA or MLAMBDA search algorithm.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/119,823 US11016199B1 (en) | 2019-12-11 | 2020-12-11 | System and method for validating GNSS ambiguities |
US17/235,333 US11681050B2 (en) | 2019-12-11 | 2021-04-20 | System and method for validating GNSS ambiguities |
US18/135,568 US20230296786A1 (en) | 2019-12-11 | 2023-04-17 | System and method for validating gnss ambiguities |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201962946576P | 2019-12-11 | 2019-12-11 | |
US17/119,823 US11016199B1 (en) | 2019-12-11 | 2020-12-11 | System and method for validating GNSS ambiguities |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US17/235,333 Continuation US11681050B2 (en) | 2019-12-11 | 2021-04-20 | System and method for validating GNSS ambiguities |
Publications (2)
Publication Number | Publication Date |
---|---|
US11016199B1 US11016199B1 (en) | 2021-05-25 |
US20210181357A1 true US20210181357A1 (en) | 2021-06-17 |
Family
ID=75981819
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US17/119,823 Active US11016199B1 (en) | 2019-12-11 | 2020-12-11 | System and method for validating GNSS ambiguities |
US17/235,333 Active US11681050B2 (en) | 2019-12-11 | 2021-04-20 | System and method for validating GNSS ambiguities |
US18/135,568 Pending US20230296786A1 (en) | 2019-12-11 | 2023-04-17 | System and method for validating gnss ambiguities |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US17/235,333 Active US11681050B2 (en) | 2019-12-11 | 2021-04-20 | System and method for validating GNSS ambiguities |
US18/135,568 Pending US20230296786A1 (en) | 2019-12-11 | 2023-04-17 | System and method for validating gnss ambiguities |
Country Status (2)
Country | Link |
---|---|
US (3) | US11016199B1 (en) |
WO (1) | WO2021119493A1 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4222609A1 (en) * | 2020-12-17 | 2023-08-09 | Swift Navigation, Inc. | System and method for fusing dead reckoning and gnss data streams |
CN113534210B (en) * | 2021-06-07 | 2022-05-31 | 湖南北斗微芯产业发展有限公司 | Ambiguity fixing method based on hybrid Kalman filtering |
EP4119988A1 (en) * | 2021-07-16 | 2023-01-18 | u-blox AG | Gnss positioning with fixing of carrier range ambiguities |
Family Cites Families (69)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4806940A (en) | 1986-04-30 | 1989-02-21 | Honeywell Inc. | Navigation mode selection apparatus |
US5451964A (en) | 1994-07-29 | 1995-09-19 | Del Norte Technology, Inc. | Method and system for resolving double difference GPS carrier phase integer ambiguity utilizing decentralized Kalman filters |
US5610614A (en) | 1995-09-13 | 1997-03-11 | Trimble Navigation Limited | Real-time kinematic initialization test system |
US5825326A (en) | 1996-07-09 | 1998-10-20 | Interstate Electronics Corporation | Real-time high-accuracy determination of integer ambiguities in a kinematic GPS receiver |
US5867411A (en) | 1996-12-19 | 1999-02-02 | The Aerospace Corporation | Kalman filter ionospheric delay estimator |
US5991691A (en) | 1997-02-20 | 1999-11-23 | Raytheon Aircraft Corporation | System and method for determining high accuracy relative position solutions between two moving platforms |
US5935196A (en) | 1997-06-11 | 1999-08-10 | Itt Manufacturing Enterprises | Technique for the use of GPS for high orbiting satellites |
US6127968A (en) | 1998-01-28 | 2000-10-03 | Trimble Navigation Limited | On-the-fly RTK positioning system with single frequency receiver |
US6205400B1 (en) | 1998-11-27 | 2001-03-20 | Ching-Fang Lin | Vehicle positioning and data integrating method and system thereof |
US6453237B1 (en) | 1999-04-23 | 2002-09-17 | Global Locate, Inc. | Method and apparatus for locating and providing services to mobile devices |
US6628231B2 (en) | 2000-02-17 | 2003-09-30 | Lockheed Martin Corp. | Location of radio frequency emitting targets |
US6408245B1 (en) | 2000-08-03 | 2002-06-18 | American Gnc Corporation | Filtering mechanization method of integrating global positioning system receiver with inertial measurement unit |
US6427122B1 (en) | 2000-12-23 | 2002-07-30 | American Gnc Corporation | Positioning and data integrating method and system thereof |
US6424914B1 (en) | 2000-12-26 | 2002-07-23 | American Gnc Corporation | Fully-coupled vehicle positioning method and system thereof |
US6735264B2 (en) | 2001-08-31 | 2004-05-11 | Rainmaker Technologies, Inc. | Compensation for non-linear distortion in a modem receiver |
JP4116792B2 (en) | 2001-12-19 | 2008-07-09 | 古野電気株式会社 | Carrier phase relative positioning device |
US20040006424A1 (en) | 2002-06-28 | 2004-01-08 | Joyce Glenn J. | Control system for tracking and targeting multiple autonomous objects |
US6856905B2 (en) | 2003-04-29 | 2005-02-15 | Garmin At, Inc. | Systems and methods for fault detection and exclusion in navigational systems |
US7148843B2 (en) | 2003-07-02 | 2006-12-12 | Thales North America, Inc. | Enhanced real time kinematics determination method and apparatus |
US7432853B2 (en) | 2003-10-28 | 2008-10-07 | Trimble Navigation Limited | Ambiguity estimation of GNSS signals for three or more carriers |
US20050114023A1 (en) | 2003-11-26 | 2005-05-26 | Williamson Walton R. | Fault-tolerant system, apparatus and method |
US20060074558A1 (en) | 2003-11-26 | 2006-04-06 | Williamson Walton R | Fault-tolerant system, apparatus and method |
US8013789B2 (en) | 2004-10-06 | 2011-09-06 | Ohio University | Systems and methods for acquisition and tracking of low CNR GPS signals |
FR2881008B1 (en) | 2005-01-20 | 2007-04-20 | Thales Sa | SATELLITE POSITIONING RECEIVER WITH IMPROVED INTEGRITY AND CONTINUITY |
US7292183B2 (en) | 2005-06-08 | 2007-11-06 | Trimble Navigation Limited | GPS reference system providing synthetic reference phases for controlling accuracy of high integrity positions |
US7312747B2 (en) * | 2005-10-03 | 2007-12-25 | Trimble Navigation Limited | Multiple-GNSS and FDMA high precision carrier-phase based positioning |
US7768449B2 (en) | 2006-01-10 | 2010-08-03 | Qualcomm Incorporated | Global navigation satellite system |
EP1972959B1 (en) | 2007-03-22 | 2012-08-15 | DLR Deutsches Zentrum für Luft- und Raumfahrt e.V. | Method for processing a set of signals of a global navigation satellite system with at least three carriers |
JP2010528320A (en) * | 2007-05-31 | 2010-08-19 | ナヴコム テクノロジー インコーポレイテッド | Reduction of distance-dependent error in real-time kinematic (RTK) positioning |
WO2009000314A1 (en) | 2007-06-22 | 2008-12-31 | Trimble Terrasat Gmbh | Position tracking device and method |
US9651667B2 (en) | 2007-06-22 | 2017-05-16 | Trimble Inc. | Combined cycle slip indicators for regionally augmented GNSS |
US8027413B2 (en) | 2007-07-11 | 2011-09-27 | The Aerospace Corporation | Ultratight coupling prefilter detection block |
WO2009058213A2 (en) | 2007-10-30 | 2009-05-07 | Trimble Navigation Limited | Generalized partial fixing |
US8694250B2 (en) | 2008-01-09 | 2014-04-08 | Trimble Navigation Limited | Processing multi-GNSS data from mixed-type receivers |
US8085190B2 (en) | 2008-03-31 | 2011-12-27 | Intel Corporation | Method and apparatus for faster global positioning system (GPS) location using a pre-computed spatial location for tracking GPS satellites |
DE112009002026B4 (en) * | 2008-08-19 | 2022-07-28 | Trimble Navigation Limited | Methods and apparatus for processing GNSS signals with selection of candidate sets |
JP5423036B2 (en) | 2009-02-18 | 2014-02-19 | セイコーエプソン株式会社 | Position calculation method and position calculation apparatus |
DE112010003681B4 (en) | 2009-09-19 | 2024-06-13 | Trimble Navigation Limited | GNSS signal processing for estimating MW biases |
US8760343B2 (en) | 2009-11-17 | 2014-06-24 | Topcon Positioning Systems, Inc. | Detection and correction of anomalous measurements and ambiguity resolution in a global navigation satellite system receiver |
US8659474B2 (en) | 2011-01-12 | 2014-02-25 | Navcom Technology, Inc. | Navigation system and method for resolving integer ambiguities using double difference ambiguity constraints |
CN103620444B (en) | 2011-03-22 | 2015-07-29 | 天宝导航有限公司 | Use the GNSS signal process being used for the ionosphere bridge joint of again restraining |
US8510041B1 (en) | 2011-05-02 | 2013-08-13 | Google Inc. | Automatic correction of trajectory data |
US20130050020A1 (en) | 2011-08-23 | 2013-02-28 | Raytheon Company | Method to handle single failure gps faults in high integrity relative positioning systems |
US9405012B2 (en) | 2012-04-12 | 2016-08-02 | Trimble Navigation Limited | Advanced global navigation satellite systems (GNSS) positioning using precise satellite information |
NL2009695C2 (en) | 2012-10-25 | 2014-05-06 | Fugro N V | Ppp-rtk method and system for gnss signal based position determination. |
US9612341B2 (en) | 2012-12-28 | 2017-04-04 | Trimble Inc. | GNSS receiver positioning system |
US10841741B2 (en) | 2015-07-07 | 2020-11-17 | Crowdcomfort, Inc. | Systems and methods for providing error correction and management in a mobile-based crowdsourcing platform |
US9784844B2 (en) | 2013-11-27 | 2017-10-10 | Honeywell International Inc. | Architectures for high integrity multi-constellation solution separation |
NL2013473B1 (en) | 2014-09-15 | 2016-09-28 | Fugro N V | Precise GNSS positioning system with improved ambiguity estimation. |
US9817129B2 (en) | 2014-10-06 | 2017-11-14 | Sierra Nevada Corporation | Monitor based ambiguity verification for enhanced guidance quality |
US9933528B2 (en) * | 2014-10-27 | 2018-04-03 | Swift Navigation, Inc. | Systems and methods for real time kinematic satellite positioning |
FR3030057B1 (en) * | 2014-12-12 | 2017-01-27 | Thales Sa | METHOD AND SYSTEM FOR VALIDATION OF SATELLITE GEOLOCATION |
US10114126B2 (en) | 2015-04-30 | 2018-10-30 | Raytheon Company | Sensor installation monitoring |
EP3109672B1 (en) * | 2015-06-24 | 2018-12-12 | Centre National d'Etudes Spatiales | Gnss receiver with a capability to resolve ambiguities using an uncombined formulation |
US20170192102A1 (en) | 2016-01-04 | 2017-07-06 | Qualcomm Incorporated | eLORAN POSITIONING VIA CROWDSOURCING |
US10802160B2 (en) | 2016-03-18 | 2020-10-13 | Deere & Company | Rapid determination of precise position by aiding data |
US10422885B2 (en) | 2016-03-18 | 2019-09-24 | Deere & Company | Rapid recovery of precise position after temporary signal loss |
US10191157B2 (en) | 2016-03-18 | 2019-01-29 | Deere & Company | Precise low-latency GNSS satellite clock estimation |
US20180091939A1 (en) | 2016-09-23 | 2018-03-29 | Qualcomm Incorporated | Geofenced access point measurement data collection |
EP3339908B1 (en) | 2016-12-23 | 2019-10-02 | u-blox AG | Distributed kalman filter architecture for carrier range ambiguity estimation |
EP3563181B1 (en) | 2016-12-30 | 2021-02-24 | u-blox AG | Gnss receiver protection levels |
US10338233B2 (en) * | 2017-04-12 | 2019-07-02 | Coherent Technical Services, Inc. | Assured validation of carrier-phase integer ambiguities for safety-of-life applications |
US10690775B2 (en) | 2017-06-29 | 2020-06-23 | Novatel Inc. | Crowdsourcing atmospheric correction data |
CN107390510B (en) * | 2017-08-22 | 2019-08-02 | 中国科学院国家授时中心 | Beidou navigation satellite precise time service unified approach based on carrier phase observable |
EP3462213A1 (en) * | 2017-09-28 | 2019-04-03 | Technische Universität München | Method for precise point positioning in a satellite navigation system |
FR3076354B1 (en) | 2017-12-28 | 2019-11-22 | Thales | METHOD FOR VERIFYING THE ENTIRE POSITION ESTIMATION OF A MOBILE CARRIER IN A SATELLITE POSITIONING MEASUREMENT SYSTEM |
US10996345B2 (en) * | 2018-06-11 | 2021-05-04 | Honeywell International Inc. | Signal fault detection for global navigation satellite system using multiple antennas |
CN109581453B (en) * | 2018-12-28 | 2023-06-02 | 北方信息控制研究院集团有限公司 | GNSS segmentation smoothing filtering method based on cycle slip detection and repair |
EP3809208A1 (en) * | 2019-10-18 | 2021-04-21 | IMEC vzw | Holographic imaging device and method |
-
2020
- 2020-12-11 WO PCT/US2020/064612 patent/WO2021119493A1/en active Application Filing
- 2020-12-11 US US17/119,823 patent/US11016199B1/en active Active
-
2021
- 2021-04-20 US US17/235,333 patent/US11681050B2/en active Active
-
2023
- 2023-04-17 US US18/135,568 patent/US20230296786A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
US20210239844A1 (en) | 2021-08-05 |
WO2021119493A1 (en) | 2021-06-17 |
US11681050B2 (en) | 2023-06-20 |
US11016199B1 (en) | 2021-05-25 |
US20230296786A1 (en) | 2023-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11543541B2 (en) | Systems and methods for high-integrity satellite positioning | |
US11681050B2 (en) | System and method for validating GNSS ambiguities | |
US20220187475A1 (en) | System and method for reconverging gnss position estimates | |
US12013472B2 (en) | System and method for fusing dead reckoning and GNSS data streams | |
US11480690B2 (en) | System and method for satellite positioning | |
US20240142637A1 (en) | System and method for gaussian process enhanced gnss corrections generation | |
US11105933B2 (en) | Method for determining a data profile for the satellite-based determination of a position of a vehicle | |
US20230184956A1 (en) | System and method for correcting satellite observations | |
US20230341563A1 (en) | System and method for computing positioning protection levels | |
US20240094415A1 (en) | System and method for detecting outliers in gnss observations | |
US12013468B2 (en) | System and method for determining GNSS corrections | |
US20220276388A1 (en) | Method for Determining a Model for Describing at least one Environment-Specific GNSS Profile | |
Stevanovic | Receiver Autonomous Integrity Monitoring Against Orbit Ephemeris Faults in Carrier Phase Differential GPS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
AS | Assignment |
Owner name: SWIFT NAVIGATION, INC., CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:GRGICH, PAUL;COLE, ANTHONY;CARCANAGUE, SEBASTIEN;REEL/FRAME:055338/0961 Effective date: 20210219 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |