GNSS signal processing methods and apparatus with geometric filter

• US 8,614,642 B2
• Filed: 02/16/2010
• Issued: 12/24/2013
• Est. Priority Date: 02/22/2009
• Status: Active Grant

• Current Assignee: Trimble Navigation Limited (Trimble Inc.)
• Original Assignee: Trimble Navigation Limited (Trimble Inc.)
• Inventors: Talbot, Nicholas Charles, Vollath, Ulrich
• Primary Examiner(s): LIU, HARRY K
• Application Number:
  US13/143,194
• Publication Number:
  US 20110267228A1
• Time in Patent Office: 1,407 Days
• Field of Search: 342/357.73
• US Class Current: 342/357.73
• CPC Class Codes: G01S 19/07 providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
• CPC Class Codes: G01S 19/32 Multimode operation in a single same satellite system, e.g. GPS L1/L2

Abstract

Methods and apparatus are provided for processing a set of GNSS signal data derived from signals of a first set of satellites having at least three carriers and signals of a second set of satellites having two carriers. A geometry filter uses a geometry filter combination to obtain an array of geometry-filter ambiguity estimates for the geometry filter combination and associated statistical information. Ionosphere filters use a two-frequency ionospheric combination to obtain an array of ionosphere-filter ambiguity estimates for the two-frequency ionospheric combinations and associated statistical information. Each two-frequency ionospheric combination comprises a geometry-free two-frequency ionospheric residual carrier-phase combination of observations of a first frequency and observations of a second frequency. Auxiliary ionosphere filters use an auxiliary ionospheric combination to obtain an array of auxiliary-ionosphere-filter ambiguity estimates for the auxiliary ionospheric combinations and associated statistical information. Each auxiliary ionospheric combination uses carrier-phase observations of a third frequency and carrier-phase observations of at least one of the first frequency and the second frequency. A combined array of ambiguity estimates is prepared for all carrier phase observations and associated statistical information by combining the arrays of the geometry filter and the ionosphere filters and the auxiliary ionosphere filters.

Claims

1. A method of processing a set of GNSS signal data derived from signals of a first set of satellites having at least three carriers and signals of a second set of satellites having two carriers, comprising:

a. applying to the set of GNSS signal data a geometry filter using a geometry filter combination to obtain an array of geometry-filter ambiguity estimates for the geometry filter combination and associated statistical information, wherein the geometry filter combination comprises one of;

a single-frequency carrier-phase and code combination, and a single-frequency carrier-phase and GRAPHIC combination;

b. applying to the set of GNSS signal data a bank of ionosphere filters each using a two-frequency ionospheric combination to obtain an array of ionosphere-filter ambiguity estimates for the two-frequency ionospheric combinations and associated statistical information, wherein each said two-frequency ionospheric combination comprises a geometry-free two-frequency ionospheric residual carrier-phase combination of observations of a first frequency and observations of a second frequency; and

c. preparing a combined array of ambiguity estimates for all carrier phase observations and associated statistical information by combining the arrays of the geometry filter and the ionosphere filters.

2. The method of claim 1, wherein the single-frequency carrier-phase and code combination is a combination of GPS L1 carrier-phase and GPS L1 code.

3. The method of claim 1, wherein the single-frequency of the single-frequency carrier-phase and GRAPHIC combination is the GPS L1 carrier frequency.

4. The method of claim 1, wherein the bank of ionosphere filters comprises one said ionosphere filter per satellite of the second set of satellites.

5. The method of claim 1, further comprising applying to the set of GNSS signal data at least one code filter using a respective geometry-free code-carrier combination to obtain an array of ambiguity estimates for the code-carrier combinations and associated statistical information.

6. The method of claim 5, wherein preparing a combined array comprises combining said arrays of said at least one code filter with the arrays of the geometry filter and the ionosphere filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information.

7. The method of claim 5, wherein each geometry-free code-carrier combination comprises one of:

(1) a combination of a first-frequency code observation with a first- and second-frequency carrier phase combination in which ionospheric bias of the first- and second-frequency carrier phase combination is matched to ionospheric bias of the first-frequency code observation; and

(2) a two-frequency narrow-lane code combination with a two-frequency wide-lane carrier-phase combination.

8. The method of claim 5, wherein said at least one code filter comprises one said code filter per satellite of the second set of satellites.

9. The method of claim 5, further comprising applying to the set of GNSS signal data derived from the first set of satellites at least one bank of additional filters, wherein each additional filter uses a geometry-free and ionosphere-free carrier-phase combination of at least three frequencies to obtain an array of ambiguity estimates for the geometry-free and ionosphere-free carrier-phase combination and associated statistical information, and wherein preparing a combined array comprises one of:

(1) combining said arrays of said additional filters with the arrays of the geometry filter and the ionosphere filters obtain said combined array of ambiguity estimates for all carries phase observations and associated statistical information; and

(2) combining said arrays of said additional filters with the arrays of said at least one code filter and the geometry filter and the ionosphere filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information.

10. The method of claim 9, wherein said at least one bank of additional filters comprises at least one said additional filter per satellite of the first set of satellites.

11. The method of claim 9, further comprising applying to the set of GNSS signal data at least one auxiliary code filter using a respective geometry-free auxiliary code-carrier combination to obtain an array of ambiguity estimates for the auxiliary code-carrier combination and associated statistical information, and wherein preparing a combined array comprises one of:

(1) combining said arrays of said additional filters with the arrays of the geometry filter and the ionosphere filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information, (2) combining the arrays of said at least one auxiliary code filter with the arrays of the geometry filter and the ionosphere filters and the at least;

one code filter to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information, and (3) combining the arrays of said at least one auxiliary code filter with the arrays of the geometry filter and the ionosphere filters and the at least one code filter and the additional filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information.

12. The method of claim 11, wherein each auxiliary code-carrier combination comprises one of:

(1) a combination of a third-frequency code observation with a first-and third-frequency carrier phase combination in which ionospheric bias of the first-and third-frequency carrier phase combination is matched to ionospheric bias of the third-frequency code observation, (2) a combination of a second- and third-frequency code combination with a second- and third-frequency carrier phase combination in which ionospheric bias of the second- and third-frequency carrier phase combination is matched to ionospheric bias of the second- and third-frequency code combination, and (3) a combination of a three-frequency narrow-lane code combination with a three-frequency wide-lane carrier phase combination in which ionospheric bias of the three-frequency wide-lane carrier-phase combination is matched to ionospheric bias of the three-frequency code combination.

13. The method of claim 11, wherein said at least one auxiliary code filter comprises one said auxiliary code filter per satellite of the first set of satellites.

14. The method of claim 1, wherein the satellites are satellites of the Global Positioning System (GPS), wherein the first set of satellites have GPS carriers L1, L2 and L5 and wherein the second set of satellites have GPS carriers L1 and L2.

15. A non-transitory computer-usable medium having a computer readable program embodied therein, said computer readable program code adapted to be executed to implement the method according to claim 1.

16. Apparatus for processing a set of GNSS signal data derived from signals of a first set of satellites having at least three carriers and signals of a second set of satellites having two carriers, comprising at least one processor with associated memory and a plurality of filters comprising:

a. a geometry filter using a geometry filter combination to produce from the set of GNSS signal data an array of geometry-filter ambiguity estimates for the geometry filter combination and associated statistical information, wherein the geometry filter combination comprises one of;

a single-frequency carrier-phase and code combination, and a single-frequency carrier-phase and GRAPHIC combination,

b. a bank of ionosphere filters each using a two-frequency ionospheric combination to obtain from the set of GNSS signal data an array of ionosphere-filter ambiguity estimates for the two-frequency ionospheric combinations and associated statistical information, wherein each said two-frequency ionospheric combination compresses a geometry-free two-frequency ionospheric residual carrier-phase combination of observations of a first frequency and observations of a second frequency; and

c. a combiner to prepare a combined array of ambiguity estimates for all carrier phase observations and associated statistical information by combining the arrays of the geometry filter and the ionosphere filters.

17. The apparatus of claim 16, wherein the singie-frequency carrier-phase and code combination is a combination of GPS L1 carrier-phase and GPS L1 code.

18. The apparatus of claim 16, wherein the single-frequency of the single-frequency carrier-phase and GRAPHIC combination is the GPS L1 carrier frequency.

19. The apparatus of claim 16, wherein the bank of ionosphere filters comprises one said ionosphere filter per satellite of the second set of satellites.

20. The apparatus of claim 16, further comprising means for applying to the set of GNSS signal data at least one code filter using a respective geometry-free code-carrier combination to obtain an array of ambiguity estimates for the code-carrier combinations and associated statistical information.

21. The apparatus of claim 20, wherein the combiner comprises means for combining said arrays of said at least one code filter with the arrays of the geometry filter and the ionosphere filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information.

22. The apparatus of claim 20, wherein each geometry-free code-carrier combination comprises one of:

(1) a combination of a first-frequency code observation with a first- and second-frequency carrier phase combination in which ionospheric bias of the first- and second-frequency carrier phase combination is matched to ionospheric bias of the first-frequency code observation; and

(2) a two-frequency narrow-lane code combination with a two-frequency wide-lane carrier-phase combination.

23. The apparatus of claim 20, wherein said at least one code filter comprises one said code filter per satellite of the second set of satellites.

24. The apparatus of claim 20, further comprising means for applying to the set of GNSS signal data derived from the first set of satellites at least one bank of additional filters, wherein each additional filter uses a geometry-free and ionosphere-free carrier-phase combination of at least three frequencies to obtain an array of ambiguity estimates for the geometry-free and ionosphere-free carrier-phase combination and associated statistical information, and wherein the combiner comprises one of:

(1) means for combining said arrays of said additional filters with the arrays of the geometry filter and the ionosphere filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information; and

(2) means for combining said arrays of said additional filters with the arrays of said at least one code filter and the geometry filter and the ionosphere filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information.

25. The apparatus of claim 24, wherein said at least one bank of additional filters comprises at least one said additional filter per satellite of the first set of satellites.

26. The apparatus of claim 24, further comprising means for applying to the set of GNSS signal data at least one auxiliary code filter using a respective geometry-free auxiliary code-carrier combination to obtain an array of ambiguity estimates for the auxiliary code-carrier combination and associated statistical information, and wherein the combiner comprises one of:

(1) means for combining said arrays of said additional filters with the arrays of the geometry filter and the ionosphere filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information, (2) means for combining the arrays of said at least one auxiliary code filter with the arrays of the geometry filter and the ionosphere filters and the at least one code filter to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information, and (3) means for combining the arrays of said at least one auxiliary code filter with the arrays of the geometry filter and the ionosphere filters and the at least one code filter and the additional filters to obtain said combined array of ambiguity estimates for all carrier phase observations and associated statistical information.

27. The apparatus of claim 26, wherein each auxiliary code-carrier combination comprises one of:

(1) a combination of a third-frequency code observation with a first-and third-frequency carrier phase combination in which ionospheric bias of the first-and third-frequency carrier phase combination is matched to ionospheric bias of the third-frequency code observation, (2) a combination of a second- and third-frequency code combination with a second- and third-frequency carrier phase combination in which ionospheric bias of the second- and third-frequency carrier phase combination is matched to ionospheric bias of the second- and third-frequency code combination, and (3) a combination of a three-frequency narrow-lane code combination with a three-frequency wide-lane carrier phase combination in which ionospheric bias of the three-frequency wide-lane carrier-phase combination is matched to ionospheric bias of the three-frequency code combination.

28. The apparatus of claim 26, comprising one said auxiliary code filter per satellite of the first set of satellites.

29. The apparatus of claim 16, wherein the satellites are satellites of the Global Positioning System (CPS), wherein the first set of satellites have GPS carriers L1, L2 and L5 and wherein the second set of satellites have GPS carriers L1, and L2.