Method for positioning using GPS in a restrictive coverage environment

• US 20060012516A1
• Filed: 07/13/2005
• Published: 01/19/2006
• Est. Priority Date: 07/15/2004
• Status: Active Grant

• Current Assignee: Novatel Incorporated (Hexagon AB)
• Original Assignee: Novatel Incorporated (Hexagon AB)
• Inventors: Ford, Thomas J., Kennedy, Sandy
• Granted Patent: US 7,528,770 B2
• Time in Patent Office: Days
• US Class Current: 342/357.590
• CPC Class Codes: G01S 19/44 Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
• CPC Class Codes: G01S 5/0036 of measured values, i.e. measurement on mobile and position calculation on base station

Abstract

A system consisting of a base GPS receiver with a clear view of the sky, one or more remote GPS receivers with restricted views of the sky and a processing center that batch processes range information provided by the GPS receivers determines the positions of the remote GPS receivers to within tight tolerances. The base GPS receiver and the remote GPS receivers produce range information based on the satellite signals that the respective receivers can track at a given time, and provide the range information to the processing center. The range information includes both code and carrier measurements for each of the signals that are being tracked by the respective GPS receivers. The center collects the range information over an extended period of time, for example, hours, days or weeks, and then batch processes the collected information in multiple passes through the data, to calculate the precise latitude, longitude and height of the receiver. The center batch processes the data using a floating ambiguity filter that, after a first pass through the data, is initialized with the position calculated in an earlier pass. The center also calculates the quality of the collected range information by fixing the position and position covariance of the floating ambiguity filter and using double differences with the base GPS receiver measurements, to ensure that the information used in subsequent position calculations is sufficiently reliable.

Claims

1. A GPS system including:

a base GPS receiver with an antenna that has a clear view of the sky;

one or more remote GPS receivers;

a data processing subsystem that collects range data from the base and remote GPS receivers over a plurality of hours and calculates the positions of the respective GPS receivers, the data processing center batch processing the data received from the respective remote receivers in multiple passes through the data using floating ambiguity filters that after the first pass are initialized with estimated positions that have been previously calculated.

2. The GPS system of claim 1 wherein the data processing center, for a second pass through the data,

further initializes the floating ambiguity filter with a previously calculated position covariance, and

calculates an estimated position using carrier phase measurements that are included in the range data.

3. The GPS system of claim 2 wherein the data processing center

updates the position and position covariance of the floating ambiguity filter at code epochs in which the remote GPS receiver has two or more GPS satellites in view, and

re-initializes the position and position covariance when a cycle slip occurs or when fewer than two satellites are in view for at least one code epoch.

4. The GPS system of claim 2 wherein the data processing center further, in one or more subsequent passes through the data,

holds the position and position covariance fixed and determines the validity of the range data for respective code epochs, and

uses only data that are determined to be valid in the calculation of the estimated position and position covariance.

5. The GPS system of claim 4 wherein the data processing center uses as the fixed position covariance a diagonal matrix with elements that are small fractions of a square meter.

6. The GPS system of claim 4 wherein the data processing center assigns weights to the data.

7. The GPS system of claim 4 wherein the data processing center determines that the range data are valid based on sums of squares of residuals of carrier phase measurements that are double differenced with corresponding measurements made by the base GPS receiver.

8. The GPS system of claim 7 wherein the data processing center estimates double difference carrier cycle ambiguities and resets the carrier ambiguities if a carrier cycle slip is detected.

9. A method of calculating GPS position in a restrictive coverage environment, the method including the steps of:

receiving range data associated with a clear view of the sky from one or more base GPS receivers;

receiving range data from one or more remote receivers that have restricted views of the sky;

determining the position of a given remote receiver by batch processing the range data received over a plurality of hours in multiple passes through the data using a floating ambiguity filter that after a first pass through the data is initialized with a previously calculated estimated position.

10. The method of claim 9 wherein the step of batch processing includes, for a second pass through the data, further initializing the floating ambiguity filter with a previously calculated position covariance and calculating the estimated position using carrier phase measurements that are included in the range data.

11. The method of claim 10 wherein the step of batch processing further includes

updating the position and position covariance of the floating ambiguity filter at code epochs in which the remote GPS receiver has two or more GPS satellites in view, and

re-initializing the position and position covariance when a cycle slip occurs or whenever fewer than two satellites are in view for at least one code epoch.

12. The method of claim 10 wherein the step of batch processing further includes in one or more subsequent passes through the data

holding the position and position covariance fixed and determining the validity of the range data for respective code epochs, and

using only data that are determined to be valid in the calculation of the estimated position and position covariance.

13. The method of claim 12 wherein the step of batch processing further includes using as the fixed position covariance a diagonal matrix with elements that are a small fraction of a square meter.

14. The method of claim 12 wherein the step of batch processing further includes assigning weights to the data.

15. The method claim 12 wherein the step of batch processing further includes determining that the range data is valid based on sums of squares of residuals of the carrier phase measurements that are double differenced with corresponding measurements made by the one or more base GPS receivers.

16. The method of claim 15 wherein the step of batch processing further includes estimating double difference carrier cycle ambiguities and resetting the carrier ambiguities if a carrier cycle slip is detected.

17. The method of claim 15 wherein the step of batch processing further includes determining that the range data is valid if

a. the time since the last cycle slip exceeds a first predetermined threshold,

b. the sum of squares of the residuals of the carrier measurement falls below a normalized threshold over an interval of continuously available measurements, and

c. the growth rate of the sum of squares of the residuals does not exceed a second predetermined threshold over an interval of continuously available measurements.

18. The method of claim 9 wherein the step of receiving range data associated with a clear view of the sky includes receiving range data from a plurality of GPS receivers that together have a clear view of the sky.

19. The method of claim 9 wherein the step of receiving range data associated with a clear view of the sky includes receiving range data from a base GPS receiver that has a clear view of the sky.