Navigation satellite system positioning involving the generation of correction information

Patent metadata

Field	Value
Publication number	US-10281587-B2
Application number	US-201615002247-A
Country	US
Kind code	B2
Filing date	Jan 20, 2016
Priority date	Dec 16, 2014
Publication date	May 7, 2019
Grant date	May 7, 2019

Abstract

The invention relates to generating correction information to be used to correct observations coming from a navigation satellite system (NSS) receiver in a region of interest. For each of a plurality of reference stations in said region, raw observations obtained by the reference station observing NSS multiple-frequency signals from a plurality of satellites over multiple epochs are received. Then, precise satellite information on the orbit position, clock offset, and biases of each satellite is obtained. For each reference station, ambiguities in the carrier phase of the received raw observations are estimated, using the precise satellite information and the position coordinates of the reference station. Geometric-free phase linear combination values are then computed based on the received raw observations together with the estimated ambiguities. The correction information is generated based on the computed geometric-free phase linear combination values. The correction information comprises, for each NSS satellite, a regional ionospheric delay function and its coefficients, those representing, per NSS satellite, the ionospheric delay in said region. The correction information is sent, for example to facilitate position determination, to the NSS receiver or to a server in charge of processing observations from the NSS receiver.

First claim

The invention claimed is:

1. Method, carried out by a computer or set of computers, for generating correction information, wherein the correction information comprises information for correcting observations useful for estimating at least one of phase ambiguities and a position of at least one global or regional navigation satellite system receiver, hereinafter abbreviated as NSS receiver, the NSS receiver being in a region of interest under a region of the ionosphere, the method comprising:
for each of a plurality of reference stations in the region of interest, receiving raw observations obtained by the reference station observing NSS multiple-frequency signals from a plurality of NSS satellites over multiple epochs;
obtaining information, hereinafter referred to as “precise satellite information”, on:
the orbit position of each one of the plurality of NSS satellites,
a clock offset of each one of the plurality of NSS satellites, and
a set of biases associated with each one of the plurality of NSS satellites, or information derived from the precise satellite information;

for each of the plurality of reference stations, estimating ambiguities in the carrier phase of the received raw observations for the reference station, using:
the precise satellite information or the information derived from the precise satellite information, and
the known position coordinates of the reference station;

computing geometric-free phase linear combination values based on the received raw observations together with the estimated ambiguities;
generating the correction information, based on the computed geometric-free phase linear combination values, wherein the correction information comprises, for each of the plurality of NSS satellites;
at least one mathematical function, each of which being hereinafter referred to as “regional ionospheric delay function”, and
coefficients of the regional ionospheric delay function, hereinafter referred to as “regional ionospheric delay coefficients”,

the at least one regional ionospheric delay function and the regional ionospheric delay coefficients representing, per NSS satellite, the ionospheric delay in the region of interest; and
sending the generated correction information.

2. Method of claim 1 , wherein the at least one regional ionospheric delay function and the regional ionospheric delay coefficients represent, per NSS satellite, the slant ionospheric delay in the region of interest.

3. Method of claim 1 , wherein the at least one regional ionospheric delay function and the regional ionospheric delay coefficients represent, per NSS satellite, the vertical ionospheric delay in the region of interest.

4. Method of claim 1 , wherein the step of generating the correction information based on the computed geometric-free phase linear combination values comprises:
computing, based on the computed geometric-free phase linear combination values from which the satellite ionospheric-phase bias has been removed using the precise satellite information, satellite-to-satellite single-differenced combination values to cancel out the effects of the reference stations' receiver phase biases;
estimating an ionosphere single layer model common to the plurality of NSS satellites, said ionosphere single layer model being hereinafter referred to as “regional ionosphere single layer model”, based on the computed satellite-to-satellite single-differenced combination values, wherein the regional ionosphere single layer model represents the vertical ionospheric delay in the region of interest;
estimating receiver ionospheric-phase biases for each reference station based on the differences between the computed geometric-free phase linear combination values from which the satellite ionospheric-phase bias has been removed using the precise satellite information, and the estimated regional ionosphere single layer model;
computing an ionospheric delay for each satellite-to-reference-station pair using the estimated receiver ionospheric-phase biases; and
for each of the plurality of NSS satellites, determining a regional ionospheric delay function and the regional ionospheric delay coefficients based on the computed ionospheric delay, wherein the regional ionospheric delay function and the regional ionospheric delay coefficients represent the ionospheric delay for the satellite, so as to generate the correction information comprising the at least one regional ionospheric delay function and the regional ionospheric delay coefficients representing the ionospheric delay in the region of interest.

5. Method of claim 4 , wherein the plurality of NSS satellites belong to at least two navigation satellite systems, and the step of estimating receiver ionospheric-phase biases comprises estimating receiver ionospheric-phase biases for each reference station and navigation satellite system.

6. Method of claim 4 , wherein the step of generating the correction information based on the computed geometric-free phase linear combination values further comprises:
for each of the plurality of NSS satellites, determining a mathematical function, hereinafter referred to as “regional accuracy indicating function”, and coefficients of the regional accuracy indicating function, hereinafter referred to as “regional accuracy indicating coefficients”, based on the residuals of the regional ionospheric delay coefficients determination, wherein the regional accuracy indicating function and the regional accuracy indicating coefficients represent the vertical accuracy of the ionospheric delay.

7. Method of claim 4 , wherein the step of determining the at least one regional ionosphere delay function and the regional ionospheric delay coefficients comprises:
improving the at least one regional ionosphere delay function and the regional ionospheric delay coefficients using the regional ionosphere single layer model associated with points in areas within the region of interest which are not covered by computed geometric-free phase linear combination values.

8. Method of claim 4 , wherein the step of estimating the regional ionosphere single layer model comprises:
improving the regional ionosphere single layer model using the global correction information associated with points provided in a grid within the region of interest to provide an absolute level of the ionospheric delay to the regional ionosphere single layer model.

9. Method of claim 1 , wherein the at least one regional ionospheric delay function is in the form of a spherical harmonics expansion.

10. Method of claim 1 , wherein the at least one regional ionospheric delay function is obtained by an Adjusted Spherical Cap Harmonics Analysis technique.

11. Method of claim 1 , wherein the aforementioned correction information is hereinafter referred to as “regional correction information”, the method further comprising receiving additional correction information, hereinafter referred to as “global correction information”, wherein the global correction information comprises at least one mathematical function and associated coefficients, the at least one mathematical function and the associated coefficients representing the ionospheric delay under the entire ionosphere.

12. Method of claim 1 , wherein the region of the ionosphere is a spherical cap.

13. Method of claim 1 , wherein the size of the region of the ionosphere comprises less than half of the ionosphere.

14. Computer program comprising instructions configured, when executed on a computer, to carry out a method according to claim 1 .

15. Method carried out by a global or regional navigation satellite system receiver, hereinaf

Assignees

• Trimble Inc

Inventors

• Drescher Ralf
• Chen Xiaoming

Classifications

• G01S19/073

  involving a network of fixed stations · CPC title

• G01S19/072 Primary

  Ionosphere corrections · CPC title

• G01S19/40 Primary

  Correcting position, velocity or attitude · CPC title

• G01S19/44

  Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method · CPC title

• G01S19/13

  Receivers · CPC title