US20080189037A1 - Method and apparatus for decompression of compressed long-term-orbit data - Google Patents
Method and apparatus for decompression of compressed long-term-orbit data Download PDFInfo
- Publication number
- US20080189037A1 US20080189037A1 US11/927,236 US92723607A US2008189037A1 US 20080189037 A1 US20080189037 A1 US 20080189037A1 US 92723607 A US92723607 A US 92723607A US 2008189037 A1 US2008189037 A1 US 2008189037A1
- Authority
- US
- United States
- Prior art keywords
- term
- long
- orbit data
- data
- orbit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/10—Artificial satellites; Systems of such satellites; Interplanetary vehicles
- B64G1/1014—Navigation satellites
-
- 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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/05—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing aiding data
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/258—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to the satellite constellation, e.g. almanac, ephemeris data, lists of satellites in view
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/27—Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
Definitions
- the present invention generally relates to generating satellite orbit information for earth orbiting satellites. More specifically, the invention relates to a method and apparatus for compressing long-term-orbit information (also known as “extended ephemeris” information) prior to distributing the information through a network or communications link.
- long-term-orbit information also known as “extended ephemeris” information
- a positioning receiver for a Global Navigation Satellite System such as the Global Positioning System (GPS) uses measurements from several satellites to compute a position.
- the process of acquiring the GNSS radio signal is enhanced in speed and sensitivity if the GNSS receiver has prior access to a model of the satellite orbit and clock.
- This model is broadcast by the GNSS satellites and is known as an ephemeris or ephemeris information. Each satellite broadcasts its own ephemeris once every 30 seconds. Once the GNSS radio signal has been acquired, the process of computing position requires the use of the ephemeris information.
- the broadcast ephemeris information is encoded in a 900 bit message within the GNSS satellite signal. It is transmitted at a rate of 50 bits per second, taking 18 seconds in all for a complete ephemeris transmission.
- the broadcast ephemeris information is typically valid for 2 to 4 hours into the future (from the time of broadcast). Before the end of the period of validity the GNSS receiver must obtain a fresh broadcast ephemeris to continue operating correctly and produce an accurate position. It is always slow (no faster than 18 seconds), frequently difficult, and sometimes impossible (in environments with very low signal strengths), for a GNSS receiver to download an ephemeris from a satellite.
- Previously available systems use a technique that collects ephemeris information at a GNSS reference station, and transmits the ephemeris to the remote GNSS receiver via a wireless transmission.
- This technique of providing the ephemeris, or equivalent data, to a GNSS receiver has become known as “Assisted-GNSS”. Since the source of ephemeris in Assisted-GNSS is the satellite signal, the ephemeris information remains valid for only a few hours. As such, the remote GNSS receiver must periodically connect to a source of ephemeris information whether that information is received directly from the satellite or from a wireless transmission. Without such a periodic update, the remote GNSS receiver will not accurately determine position.
- the present invention is a method and apparatus for generating satellite orbit data that is valid for extend periods of time into the future, i.e., long-term-orbits or LTO (also referred to as extended ephemeris), compressing the LTO before transmission to a GNSS receiver, and decompressing the LTO at the GNSS receiver.
- LTO may contain future satellite trajectory information and/or satellite clock information.
- the LTO is derived by receiving at one or more satellite tracking stations the signals from at least one satellite and determining satellite tracking information (STI) from the received signals.
- STI contains present satellite orbit trajectory data and satellite clock information (e.g., ephemeris).
- PRN pseudo-random noise
- LTO is often formatted as sequential blocks of data, each block valid for a period of time. Some of the data fields change little or not at all from block to block, and other data fields can be predicted from other data. Thus some of the data in the LTO may be described as redundant information. To compress the LTO, redundancy is reduced amongst the individual LTO associated with each satellite. In this manner, a data set containing a set of LTO for a plurality of satellites is substantially smaller than the uncompressed set of LTO for the same plurality of satellites.
- FIG. 1 depicts a system for creating and distributing long term satellite tracking data (LTO) to remote GNSS receivers;
- LTO long term satellite tracking data
- FIG. 2 depicts a method for forming the LTO from the satellite measurements made at satellite tracking stations
- FIG. 3 depicts a timeline of LTO data that conforms to the broadcast ephemeris format models as described in ICD-GPS-200C yet spans many hours;
- FIG. 4 depicts a flow diagram of a method that uses a least squares estimation technique to update parameters in an orbit trajectory model
- FIG. 5 depicts the error in the orbit model derived from the LTO, and compares the error to the error in the broadcast ephemeris;
- FIG. 6 depicts an example of a data table that could be used in an LTO database
- FIG. 7 depicts a flow diagram depicting a method of compressing a set of LTO in accordance with one embodiment of the invention
- FIG. 8 depicts a graphical illustration of an LTO compression data set in accordance with one embodiment of the invention.
- FIG. 9 depicts a flow diagram of a method of decompressing a compressed LTO data set in accordance with one embodiment of the invention.
- FIG. 1 depicts a block diagram of a system 100 for creating and distributing long term satellite tracking data (LTO).
- the global navigation satellite system may include the global positioning system (GPS), GLONASS, GALILEO, or other satellite systems that may use LTO to enhance the performance of the receiver.
- GPS global positioning system
- GLONASS global positioning system
- GALILEO global positioning system
- the following disclosure uses GPS as an illustrative system within which the invention operates. From the following disclosure, those skilled in the art will be able to practice the invention in conjunction with other satellite systems.
- a network of GNSS tracking stations 102 is used to collect measurement data from the GNSS satellites 104 .
- the network could comprise several tracking stations that collect satellite tracking information (STI) from all the satellites in the constellation, or a few tracking stations, or a single tracking station that only collects STI for a particular region of the world.
- An LTO collection and computation server 106 collects and processes the measurement data (this measurement data is referred to herein as satellite tracking information (STI) that includes at least one of: code phase measurements, carrier phase measurements, Doppler measurements, or ephemeris data).
- STI satellite tracking information
- measurement data is obtained from both the L 1 and L 2 frequencies on which the GPS satellites transmit.
- Alternative embodiments may use only one of these frequencies, and/or other frequencies used by other satellite systems or by future versions of the GNSS system.
- the server 106 produces: 1) accurate satellite tracking data (STD) (e.g., a trajectory of each satellite and/or a clock offset measurement) during the data collection period, 2) a prediction of the future STD of each satellite, and 3) models that match the future STD of each satellite to form an LTO for each satellite.
- the server 106 comprises a central processing unit (CPU) 118 , support circuits 122 , and memory 120 .
- the CPU 118 may be any one of the many CPUs available on the market to perform general computing.
- the CPU may be a specific purpose processor such as an application specific integrated circuit (ASIC) that is designed to process satellite tracking information.
- ASIC application specific integrated circuit
- the support circuits 122 are well known circuits such as clock circuits, cache, power supplies and the like.
- the memory 120 may be read only memory, random access memory, disk drive storage, removable storage or any combination thereof.
- the memory 120 stores executable software, e.g., LTO software 124 , that, when executed by the CPU 118 , causes the system 100 to operate in accordance with the present invention.
- the set of satellite trajectory and clock data produced by the LTO software 124 constitutes the STD information, and is stored in an STD database 108 .
- a distribution server 110 accesses the database 108 to gather the most recent set of data, formats the data using the trajectory conversion software 111 according to the relevant interface standard, and distributes the formatted data to GNSS devices 112 that require satellite orbit information.
- the distribution process may be by some form of wireless communications system 114 , or over the Internet 116 , or a combination of both, or by some other means of communication.
- the GNSS devices 112 Once the GNSS devices 112 have received the orbit data, they may operate continually for many days without needing to download fresh broadcast ephemeris from the satellites or from any other source.
- the orbit data distributed to the GNSS devices may be in the same format as the broadcast ephemeris or may be some other model format that is defined by the GPS device.
- this orbit data is generally referred to as a satellite tracking model (STM).
- STM satellite tracking model
- the loading of the STM into the GNSS receiver can be accomplished in many ways.
- PDA personal digital assistant
- direct connection to a network, or a wireless technology, such as Bluetooth or a cellular network are a few examples of how the ephemeris data can be transferred to the receiver.
- the transmission is generally accomplished by broadcasting the LTO (or a model representing all or a portion of the LTO) without knowledge of the specific location of the GNSS receiver.
- the distribution server does not require the GNSS receiver to send any information through the network to the distribution server.
- the data transmitted by the GNSS satellites can be used to determine the range, range-rate and clock offsets to the GNSS satellites from a set of tracking stations.
- This set of observations generated by the tracking stations 102 is used in the orbit determination process, and in the estimation of the satellite clock characteristics.
- the set of monitoring stations 102 could be a single station, a public network such as the Continuously Operating Reference System (CORS), or a privately owned and/or operated network.
- CORS Continuously Operating Reference System
- FIG. 2 illustrates the preferred embodiment of a process for computing LTO.
- the process begins at step 202 with the collection of satellite measurements from the network of tracking stations. Measurements such as code phase, (CP), carrier phase (CPH), and Doppler may be used for GNSS satellite tracking information.
- the measurements are used to compute the satellite trajectories and clock offsets over the periods during which the data was collected. This step is performed using standard GPS processing techniques and software packages well known in the art. Examples of this type of software are GIPSY from the Jet Propulsion Laboratory (JPL), GEODYN from NASA Goddard Space Flight Center (GSFC), and the commercial product, MicroCosm, from Van Martin Systems.
- JPL Jet Propulsion Laboratory
- GSFC NASA Goddard Space Flight Center
- the satellite trajectories and clock offsets from step 204 are propagated into the future with the same software package, using standard orbit models, such as gravity, drag, solar radiation pressure, tides, third body effects, precession, nutation, and other conservative and non-conservative forces effecting the satellite trajectory.
- standard orbit models such as gravity, drag, solar radiation pressure, tides, third body effects, precession, nutation, and other conservative and non-conservative forces effecting the satellite trajectory.
- These are normally the same force models that are used in the estimation of the satellite orbits during the data fit interval.
- a subset of these models, such as those for drag and solar radiation pressure are adjusted during the orbit estimation process described in step 204 to best fit the trajectory.
- This combination of known and estimated force models and parameters is used in the propagation 206 to provide the propagated orbit for time outside the data fit interval.
- the clock offsets for GPS satellites are typically very small, and change linearly over time. These clock offsets are propagated into the future using standard models, such as a second order model
- the propagated satellite trajectories and/or clock offsets are stored as STD in a database.
- the trajectory conversion software converts the LTO data into a model and format expected by the GNSS device to which the model is to be provided.
- the prescribed model or information is output.
- the preferred embodiment of the model is the GNSS ephemeris model as described in ICD-GPS-200 and an ephemeris model is generated from the LTO for each 4 hour period as illustrated in the timeline 300 of FIG. 3 , i.e., a different model 301 , 302 and so on is generated for each six hour period.
- the plurality of models 301 , 302 and so on cumulatively span the length of the available LTO.
- the satellite trajectories and clock offsets may be estimated using the data broadcast by the satellites and the standard equations given in ICD-GPS-200c.
- the orbit model is a mathematical representation of the satellite trajectory that describes the trajectory as a function of a small number of variables and eliminates the need to provide satellite position vectors explicitly as a table of time vs. satellite positions.
- An example of an ephemeris model is the classic six element Keplerian orbital model. Although this model lacks long term accuracy, it is a functional ephemeris model for providing satellite trajectory information as a function of a small number of variables.
- the model used to describe the trajectory is GNSS standard ephemeris, specified in ICD-GPS-200c, following the same conventions and units. This is the preferred method to provide maximum compatibility with existing GNSS receivers.
- other orbit models could also be used to represent the satellite trajectory. Orbit models can be selected to provide increased accuracy, longer duration fits, more compact representation of the trajectory, or other optimizations required in an application.
- This invention is different from the current art in that the orbit model provided to the GNSS device is not the ephemeris data broadcast by the GNSS satellites.
- Current art downloads the ephemeris broadcast from the GNSS satellites and retransmits that data to GNSS devices.
- the broadcast ephemeris data is not required at any stage and is not used in the preferred implementation.
- the broadcast ephemeris data provided by the GNSS satellites cover a specific time period (typically 4 hours) and the end of that time the information becomes unusable. For example, if a device receives a broadcast ephemeris that will expire in 5 minutes, the device would need the new broadcast ephemeris before operating outside that 5 minute interval.
- the STD may be formatted for the time period required by the device. This time period may be for the current time forward or may be for some time interval in the future.
- a device may request orbit information in the standard GPS ephemeris format for the current time. In this case, the ephemeris provided to the device would be valid for the next 6 hours. The device could request orbit information for the next 12 hours in the standard GPS format which could be supplied as two six hour ephemeris orbit models.
- different orbit models and formats that support different accuracies and standards can be generated from the LTO.
- Fitting the LTO to the desired orbit model can be accomplished in a number of mathematical methods.
- the preferred embodiment is a least-squares fit of the orbit model parameters to the trajectory data.
- Other methods such as Kalman filters or other estimators can also be used to obtain the orbit model parameters that best fit the trajectory data.
- FIG. 4 depicts a flow diagram of a method of generating an orbit model using a least squares estimation technique.
- LTO is a table representation of time, position, and clock offset for each satellite, as shown in FIG. 6 .
- the time, position, and clock offset can be in any time/coordinate system.
- the time/coordinate system is GPS time and Earth-Centered-Earth-Fixed (ECEF) position in the World Geodetic Survey 1984 (WGS-84) reference frame.
- ECEF Earth-Centered-Earth-Fixed
- the STD for the desired time interval is extracted from the STD database.
- the orbit model parameters are initialized to the orbit model values obtained by a similar process for the previous interval. This guarantees that the initial orbit model parameters are a good fit at least for the beginning of the desired time interval. The rest of the process 400 will ensure that the parameters are adjusted so that they become a good fit for the entire time interval.
- the orbit model is used to predict what the trajectory would be, the predicted data is denoted the “Model Trajectory Data” (MTD). If the model were perfect, the MTD would coincide exactly with the STD.
- the MTD and STD are compared to see how closely the orbit model fits the orbit data. In the preferred embodiment, the comparison step 406 is performed by summing the squares of the differences between each trajectory point in the STD and the corresponding point in the MTD, and comparing the resulting sum to a threshold. If the fit is “good”, the model parameters are deemed “good” and the process stops at step 410 . If the fit is not good then the model parameters are adjusted at step 408 .
- Steps 404 , 406 and 408 are repeated until the model parameters are found that fit the STD well.
- the preferred embodiment is to use orbit models with parameters that have been defined in well-known standards.
- the ephemeris parameters defined in the GPS interface control document, ICD-GPS-200c are used.
- the ICD-GPS-200c definition includes a bit that specifies a 4-hour fit or a 6-hour fit.
- the satellite data is broadcast in 4-hour fits and, by the time this data is obtained by the observer of the satellite, the data is often near the end of its fit interval.
- sequential 6 hour windows of STD are used to create 6-hour ephemeris models, using the technique described in FIG. 4 and the accompanying text.
- FIG. 5 shows an example of Satellite Tracking Data (STD) that was generated for a time interval of greater than six hours. Then, using the technique described by FIG. 4 and accompanying text, parameters of an ICD-GPS-200c ephemeris model were adjusted to give a best fit to 6 hours of the STD. The orbit modeled by this 6-hour ephemeris was then compared to the true trajectory and, for comparison; the true trajectory was also compared to the orbit modeled by the broadcast ephemeris. The results are shown in graph 500 of FIG. 5 .
- the vertical axis 502 represents position error in meters and the horizontal axis 504 represents time in hours.
- the graph 500 illustrates how the broadcast ephemeris 506 loses validity while the ephemeris 508 created by this invention maintains its validity with approximately one meter of error.
- the clock offset of GPS satellites is easily modeled by three parameters.
- the measured clock offset is modeled by the three parameters defined in ICD-GPS-200c. These parameters represent clock offset, drift, and drift rate. The parameters are adjusted in a similar way to the method 400 described above to give a model that best fits the measured data over the time interval.
- Alternative embodiments may use longer fit intervals, such as 8, 14, 26, 50, 74, 98, 122, or 146 hours for each ephemeris model. These fit intervals are envisaged in ICD-GPS-200c, but are seldom, if ever, available from the broadcast ephemeris. Under the current invention, models with these fit intervals may be generated even when the broadcast ephemeris is limited to a 4-hour fit interval.
- Alternative embodiments of the STD data may include observed satellite velocity, acceleration, clock drift, or clock drift rate and these terms may be used in the process of fitting a model in ways which are well known in the art.
- an orbit model uses the spare data bits in the current ephemeris format of a conventional GPS signal to provide additional model parameters that would improve the data fit over long time intervals. For example, subframe 1 has 87 spare bits that are available for additional parameters. This technique allows for more parameters to describe the orbital motion of the satellites without compromising the standard data format. This new ephemeris model is based on the current ephemeris model with additional correction terms used to augment the model to support the longer fit intervals with greater accuracy.
- Yet another embodiment of an orbit model is to develop a new set of orbital parameters that describe the satellite orbit which are different, in part or in their entirety, from the GPS ephemeris model parameters. With the goal of making the fit interval longer, different parameters may provide a better description of the satellite orbit. This new set of parameters could be defined such that they would fit into the existing data structures, however, their implementation and algorithms for use would be different.
- an orbit model would be to develop a new set of orbital parameters that would not fit into the existing GPS ephemeris model format. This new set of parameters would be developed to better address the trade-off between the number of parameters required, the fit interval, and the orbit accuracy resulting from the model.
- An example of this type of ephemeris parameter set is Brouwer's theory that could be used as is or modified to account for GPS specific terms. Brouwer's theory as described in Brouwer, D. “Solution of the Problem of Artificial Satellite Theory without Drag”, Astron J. 64: 378-397, November 1959 is limited to satellites in nearly circular orbits such as GPS satellites.
- Another embodiment is to use a subset of the standard ephemeris parameters defined in ICD-GPS-200c. This approach is particularly useful when bandwidth and/or packet size is limited in the communication link that will be used to convey the orbit model to the Remote GPS Receiver.
- the fifteen orbit parameters described above, and in ICD-GPS-200c may be reduced to a subset of 9 parameters, by setting all harmonic terms in the model to zero:
- Process 400 is then executed using this subset of parameters. This reduces the amount of data that must be sent to the Remote GPS Receiver.
- the receiver can then reconstruct a standard ephemeris model by setting the “missing” harmonic terms to zero.
- Some of these approaches reduce the ability of the model to fit the data over a period of time (e.g., six hours). In this case, the fit interval may be reduced (e.g. to four hours) to compensate.
- Each satellite in the plurality of GNSS satellites has an associated LTO. Together all the LTOs for a plurality of satellites form an LTO set.
- the present invention compresses the LTO set by reducing information redundancy from the set.
- the LTO set becomes a compressed set comprising an ephemeris structure (denoted lto — 0) plus a sequence of “deltas” to the ephemeris structure.
- client software in the GNSS receiver extracts the ephemeris structure of the lto — 0 and adds a corresponding delta to lto — 0 to form each LTO in the set.
- FIG. 7 depicts a flow diagram of one embodiment of a compression method 700 .
- the method 700 begins at step 702 and proceeds to step 704 , where the method 700 being executed by the server ( 106 of FIG. 1 ) accesses a set of LTO. Method 700 applies to each satellite separately.
- the method 700 designates one LTO as a reference LTO (lto — 0). In the one embodiment, lto — 0 is the first LTO in the set of LTO for this particular satellite.
- lto — 0 could be any LTO in the set of LTO, or it could be another structure with the same fields, such as the current satellite ephemeris, or almanac, In another alternative embodiment lto — 0 could derived by using information from the set of LTO, (such as forming lto — 0 from the average LTO, the median LTO, or some other derived LTO). In yet another alternative embodiment lto — 0 could be arbitrarily assigned.
- variable “i” is set to 1 and, at step 710 , the method selects the next LTO (lto_i).
- the difference ( ⁇ lto_i) on a field-by-field basis between the reference LTO (lto — 0) and the current selected LTO (lto_i) is computed at step 712 .
- fields that do not change from LTO to LTO have a value of zero to indicate no change. In this manner, the LTO of each satellite is compressed.
- the method 700 queries whether the method 700 has completed compressing all the satellite LTOs. If the query is negatively answered, the method 700 proceeds to step 720 , where “i” is increased by 1 to enable step 710 to select the next LTO for compression.
- step 722 the method 700 arranges lto — 0 and the ⁇ lto_i values into a data set.
- the method 700 also creates an overhead record—a dynamic range map—that is added to the compressed data set.
- the dynamic range map contains a dynamic range value (the number of bits) for each field in an LTO. For example, if a field contains no change, then the field in that ⁇ lto_i is assigned a dynamic range of zero and no information needs to be included in that particular field in the ⁇ lto_i.
- the method 700 ends at step 724 . It is understood that the data set could be implemented in different ways, for example: as a file, or a series of communications packets, or some other means of transferring the information.
- FIG. 8 depicts a graphical illustration of a compressed data set 800 comprising an uncompressed reference LTO (lto — 0) 802 , and a series of delta LTOs ( ⁇ lto_i) 806 .
- LTO uncompressed reference LTO
- ⁇ lto_i delta LTOs
- at least one of a dynamic range map 804 and a bitmask 808 may also form a portion of the C-LTO data set 800 .
- TOC (16 bits) changes with every LTO, it is identified in the dynamic range map as having a range of 16 bits and all 16 bits of TOC information appear in the ⁇ lto_i. However the af2 field does not change and has a range of zero and therefore does not appear at all in the ⁇ lto_i.
- Such compression to reduce redundancy is performed on all the fields.
- the dynamic range map 804 is provided in advance or agreed upon in a specification.
- the dynamic range map 804 there is no need to include the dynamic range map 804 as part of the LTO data being distributed because it is known.
- the number of bits of each delta LTO field is fixed.
- This technique has the advantage of simplicity and ease of encoding in a specification but it is less efficient because it has to use more bits per delta field to cover all possible variations of each parameter in the LTO data.
- a further refinement of this technique uses a bitmask 808 with one bit for each field indicating whether the compressed format has been used for each particular field. This technique provides a fallback mechanism for the case where one or more parameters exceed the specified dynamic range.
- Each of the parameters of the LTO has an associated scale factor.
- the LTO comprises standard ephemeris parameters
- one such parameter is i 0 (orbit inclination angle) which has units of semi-circles and a scale factor of 2 ⁇ 31 .
- the scale factors could be changed so that fewer bits are used to carry the information.
- certain terms may be eliminated from the set of delta-LTOs, and reconstructed.
- the four terms: M 0 , i 0 , ⁇ 0 and af0 do not have to be sent more than once per satellite per LTO, since, for each subsequent 6 hour segment, these terms can be reconstructed from the previously distributed values.
- the method propagates forward n times using the above equations, and computes the appropriate values on the right hand side of the equations at each step.
- FIG. 9 depicts a flow diagram of a method 900 that is performed by the GNSS receiver to unpack the compressed data set and recreate LTO.
- the method 900 starts at step 902 and proceeds to step 904 , wherein the compressed data set is received as described above with respect to a decompressed LTO.
- the reference LTO (lto — 0) is extracted from the data set and, at step 908 , the variable “i” is set to 1.
- the method 900 selects a ⁇ lto_i and, at step 912 , adds lto — 0 to ⁇ lto_i on a field by field basis.
- the dynamic range map is used to identify where the data is located within ⁇ lto_i and/or to identify bit sizes of each parameter with ⁇ lto_i.
- the value of lto_i is saved in memory.
- the method 900 queries whether all the ⁇ lto_i have been processed. If the query is negatively answered, the method 900 proceeds to step 922 to increase i by 1 and return to step 910 to select the next ⁇ lto_i for decompression. If, at step 916 , the decompression method 900 has processed all ⁇ lto_i, then the method 900 continues to step 918 . At step 918 , the method 900 arranges the stored LTO values (lto — 0 and all lto_i) into an LTO set that is then used by the GNSS receiver in the manner described above with respect to uncompressed LTO.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Power Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Description
- This application is a continuation of a co-pending U.S. application Ser. No. 11/724,627, filed Mar. 15, 2007, which is a continuation-in-part application of U.S. patent application Ser. No. 11/333,787 filed Jan. 17, 2006, which is a continuation-in-part application of co-pending U.S. patent application Ser. No. 09/993,335, filed Nov. 6, 2001, which is a continuation-in-part of U.S. patent application Ser. No. 09/884,874, filed Jun. 19, 2001, now U.S. Pat. No. 6,560,534, which is a continuation-in-part of U.S. patent application Ser. No. 09/875,809, filed Jun. 6, 2001, now U.S. Pat. No. 6,542,820. This application contains subject matter that is related to co-pending U.S. patent application Ser. No. 09/715,860, filed Nov. 17, 2000, now U.S. Pat. No. 6,417,801. Each of the aforementioned related patents and/or patent applications is herein incorporated by reference.
- 1. Field of the Invention
- The present invention generally relates to generating satellite orbit information for earth orbiting satellites. More specifically, the invention relates to a method and apparatus for compressing long-term-orbit information (also known as “extended ephemeris” information) prior to distributing the information through a network or communications link.
- 2. Description of the Related Art
- A positioning receiver for a Global Navigation Satellite System (GNSS), such as the Global Positioning System (GPS), uses measurements from several satellites to compute a position. The process of acquiring the GNSS radio signal is enhanced in speed and sensitivity if the GNSS receiver has prior access to a model of the satellite orbit and clock. This model is broadcast by the GNSS satellites and is known as an ephemeris or ephemeris information. Each satellite broadcasts its own ephemeris once every 30 seconds. Once the GNSS radio signal has been acquired, the process of computing position requires the use of the ephemeris information.
- The broadcast ephemeris information is encoded in a 900 bit message within the GNSS satellite signal. It is transmitted at a rate of 50 bits per second, taking 18 seconds in all for a complete ephemeris transmission. The broadcast ephemeris information is typically valid for 2 to 4 hours into the future (from the time of broadcast). Before the end of the period of validity the GNSS receiver must obtain a fresh broadcast ephemeris to continue operating correctly and produce an accurate position. It is always slow (no faster than 18 seconds), frequently difficult, and sometimes impossible (in environments with very low signal strengths), for a GNSS receiver to download an ephemeris from a satellite. For these reasons it has long been known that it is advantageous to send the ephemeris to a GNSS receiver by some other means in lieu of awaiting the transmission from the satellite. Previously available systems use a technique that collects ephemeris information at a GNSS reference station, and transmits the ephemeris to the remote GNSS receiver via a wireless transmission. This technique of providing the ephemeris, or equivalent data, to a GNSS receiver has become known as “Assisted-GNSS”. Since the source of ephemeris in Assisted-GNSS is the satellite signal, the ephemeris information remains valid for only a few hours. As such, the remote GNSS receiver must periodically connect to a source of ephemeris information whether that information is received directly from the satellite or from a wireless transmission. Without such a periodic update, the remote GNSS receiver will not accurately determine position.
- Until recently, the deficiency of the current art was that there was no source of satellite trajectory and clock information that is valid for longer than a few hours into the future, and it can be expensive to send the ephemeris information repeatedly to the many remote devices that may need it. Moreover, mobile devices may be out of contact from the source of the Assisted-GNSS information when their current ephemeris becomes invalid. Consequently, long-term-orbit models (sometimes referred to as extended ephemeris) are used to enable an assisted-GNSS receiver to operate for a long period of time before new ephemeris is required. However, long-term-orbit models that extend over a long period of time can contain many bits that form large files. These large files utilize extensive bandwidth when being transmitted through a network to a GNSS receiver.
- Therefore, there is a need in the art for a method and apparatus for providing satellite trajectory and clock information that is valid for an extended period into the future, e.g., many days into the future, and sending that information to a GNSS receiver in a compressed form.
- The present invention is a method and apparatus for generating satellite orbit data that is valid for extend periods of time into the future, i.e., long-term-orbits or LTO (also referred to as extended ephemeris), compressing the LTO before transmission to a GNSS receiver, and decompressing the LTO at the GNSS receiver. The LTO may contain future satellite trajectory information and/or satellite clock information. The LTO is derived by receiving at one or more satellite tracking stations the signals from at least one satellite and determining satellite tracking information (STI) from the received signals. STI contains present satellite orbit trajectory data and satellite clock information (e.g., ephemeris). The LTO is created for each satellite, as identified by its unique pseudo-random noise (PRN) code. LTO is often formatted as sequential blocks of data, each block valid for a period of time. Some of the data fields change little or not at all from block to block, and other data fields can be predicted from other data. Thus some of the data in the LTO may be described as redundant information. To compress the LTO, redundancy is reduced amongst the individual LTO associated with each satellite. In this manner, a data set containing a set of LTO for a plurality of satellites is substantially smaller than the uncompressed set of LTO for the same plurality of satellites.
- So that the manner in which the above recited features of the present invention are attained and can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to the embodiments thereof which are illustrated in the appended drawings.
- It is to be noted, however, that the appended drawings illustrate only typical embodiments of this invention and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective embodiments.
-
FIG. 1 depicts a system for creating and distributing long term satellite tracking data (LTO) to remote GNSS receivers; -
FIG. 2 depicts a method for forming the LTO from the satellite measurements made at satellite tracking stations; -
FIG. 3 depicts a timeline of LTO data that conforms to the broadcast ephemeris format models as described in ICD-GPS-200C yet spans many hours; -
FIG. 4 depicts a flow diagram of a method that uses a least squares estimation technique to update parameters in an orbit trajectory model; -
FIG. 5 depicts the error in the orbit model derived from the LTO, and compares the error to the error in the broadcast ephemeris; -
FIG. 6 depicts an example of a data table that could be used in an LTO database; -
FIG. 7 depicts a flow diagram depicting a method of compressing a set of LTO in accordance with one embodiment of the invention; -
FIG. 8 depicts a graphical illustration of an LTO compression data set in accordance with one embodiment of the invention; and -
FIG. 9 depicts a flow diagram of a method of decompressing a compressed LTO data set in accordance with one embodiment of the invention. -
FIG. 1 depicts a block diagram of asystem 100 for creating and distributing long term satellite tracking data (LTO). The global navigation satellite system (GNSS) may include the global positioning system (GPS), GLONASS, GALILEO, or other satellite systems that may use LTO to enhance the performance of the receiver. The following disclosure uses GPS as an illustrative system within which the invention operates. From the following disclosure, those skilled in the art will be able to practice the invention in conjunction with other satellite systems. - A network of
GNSS tracking stations 102 is used to collect measurement data from theGNSS satellites 104. Such a network is described in detail in U.S. patent application Ser. No. 09/615,105, filed Jul. 13, 2000. The network could comprise several tracking stations that collect satellite tracking information (STI) from all the satellites in the constellation, or a few tracking stations, or a single tracking station that only collects STI for a particular region of the world. An LTO collection andcomputation server 106 collects and processes the measurement data (this measurement data is referred to herein as satellite tracking information (STI) that includes at least one of: code phase measurements, carrier phase measurements, Doppler measurements, or ephemeris data). In the preferred embodiment, measurement data is obtained from both the L1 and L2 frequencies on which the GPS satellites transmit. Alternative embodiments may use only one of these frequencies, and/or other frequencies used by other satellite systems or by future versions of the GNSS system. Theserver 106 produces: 1) accurate satellite tracking data (STD) (e.g., a trajectory of each satellite and/or a clock offset measurement) during the data collection period, 2) a prediction of the future STD of each satellite, and 3) models that match the future STD of each satellite to form an LTO for each satellite. Theserver 106 comprises a central processing unit (CPU) 118,support circuits 122, andmemory 120. TheCPU 118 may be any one of the many CPUs available on the market to perform general computing. Alternatively, the CPU may be a specific purpose processor such as an application specific integrated circuit (ASIC) that is designed to process satellite tracking information. Thesupport circuits 122 are well known circuits such as clock circuits, cache, power supplies and the like. Thememory 120 may be read only memory, random access memory, disk drive storage, removable storage or any combination thereof. Thememory 120 stores executable software, e.g.,LTO software 124, that, when executed by theCPU 118, causes thesystem 100 to operate in accordance with the present invention. - The set of satellite trajectory and clock data produced by the
LTO software 124 constitutes the STD information, and is stored in anSTD database 108. Adistribution server 110 accesses thedatabase 108 to gather the most recent set of data, formats the data using thetrajectory conversion software 111 according to the relevant interface standard, and distributes the formatted data toGNSS devices 112 that require satellite orbit information. The distribution process may be by some form ofwireless communications system 114, or over theInternet 116, or a combination of both, or by some other means of communication. Once theGNSS devices 112 have received the orbit data, they may operate continually for many days without needing to download fresh broadcast ephemeris from the satellites or from any other source. The orbit data distributed to the GNSS devices may be in the same format as the broadcast ephemeris or may be some other model format that is defined by the GPS device. Herein this orbit data is generally referred to as a satellite tracking model (STM). The loading of the STM into the GNSS receiver can be accomplished in many ways. Using the cradle for a personal digital assistant (PDA), direct connection to a network, or a wireless technology, such as Bluetooth or a cellular network, are a few examples of how the ephemeris data can be transferred to the receiver. The transmission is generally accomplished by broadcasting the LTO (or a model representing all or a portion of the LTO) without knowledge of the specific location of the GNSS receiver. As such, the distribution server does not require the GNSS receiver to send any information through the network to the distribution server. - Since GNSS is a ranging system in and of itself, the data transmitted by the GNSS satellites can be used to determine the range, range-rate and clock offsets to the GNSS satellites from a set of tracking stations. This set of observations generated by the tracking
stations 102 is used in the orbit determination process, and in the estimation of the satellite clock characteristics. The set ofmonitoring stations 102 could be a single station, a public network such as the Continuously Operating Reference System (CORS), or a privately owned and/or operated network. -
FIG. 2 illustrates the preferred embodiment of a process for computing LTO. The process begins atstep 202 with the collection of satellite measurements from the network of tracking stations. Measurements such as code phase, (CP), carrier phase (CPH), and Doppler may be used for GNSS satellite tracking information. Atstep 204, the measurements are used to compute the satellite trajectories and clock offsets over the periods during which the data was collected. This step is performed using standard GPS processing techniques and software packages well known in the art. Examples of this type of software are GIPSY from the Jet Propulsion Laboratory (JPL), GEODYN from NASA Goddard Space Flight Center (GSFC), and the commercial product, MicroCosm, from Van Martin Systems. - At
step 206, the satellite trajectories and clock offsets fromstep 204 are propagated into the future with the same software package, using standard orbit models, such as gravity, drag, solar radiation pressure, tides, third body effects, precession, nutation, and other conservative and non-conservative forces effecting the satellite trajectory. These are normally the same force models that are used in the estimation of the satellite orbits during the data fit interval. A subset of these models, such as those for drag and solar radiation pressure, are adjusted during the orbit estimation process described instep 204 to best fit the trajectory. This combination of known and estimated force models and parameters is used in thepropagation 206 to provide the propagated orbit for time outside the data fit interval. The clock offsets for GPS satellites are typically very small, and change linearly over time. These clock offsets are propagated into the future using standard models, such as a second order model containing clock offset, drift, and drift rate. - At
step 208, the propagated satellite trajectories and/or clock offsets are stored as STD in a database. Atstep 210, the trajectory conversion software converts the LTO data into a model and format expected by the GNSS device to which the model is to be provided. Atstep 212, the prescribed model or information is output. For use with existing GNSS receivers, the preferred embodiment of the model is the GNSS ephemeris model as described in ICD-GPS-200 and an ephemeris model is generated from the LTO for each 4 hour period as illustrated in thetimeline 300 ofFIG. 3 , i.e., adifferent model models - In an alternate embodiment, at
step 204, the satellite trajectories and clock offsets may be estimated using the data broadcast by the satellites and the standard equations given in ICD-GPS-200c. - The orbit model is a mathematical representation of the satellite trajectory that describes the trajectory as a function of a small number of variables and eliminates the need to provide satellite position vectors explicitly as a table of time vs. satellite positions. An example of an ephemeris model is the classic six element Keplerian orbital model. Although this model lacks long term accuracy, it is a functional ephemeris model for providing satellite trajectory information as a function of a small number of variables. In the preferred embodiment, the model used to describe the trajectory is GNSS standard ephemeris, specified in ICD-GPS-200c, following the same conventions and units. This is the preferred method to provide maximum compatibility with existing GNSS receivers. However, other orbit models could also be used to represent the satellite trajectory. Orbit models can be selected to provide increased accuracy, longer duration fits, more compact representation of the trajectory, or other optimizations required in an application.
- This invention is different from the current art in that the orbit model provided to the GNSS device is not the ephemeris data broadcast by the GNSS satellites. Current art downloads the ephemeris broadcast from the GNSS satellites and retransmits that data to GNSS devices. In this invention, the broadcast ephemeris data is not required at any stage and is not used in the preferred implementation.
- The broadcast ephemeris data provided by the GNSS satellites, specifically GPS satellites, cover a specific time period (typically 4 hours) and the end of that time the information becomes unusable. For example, if a device receives a broadcast ephemeris that will expire in 5 minutes, the device would need the new broadcast ephemeris before operating outside that 5 minute interval. With this invention, the STD may be formatted for the time period required by the device. This time period may be for the current time forward or may be for some time interval in the future. For example, a device may request orbit information in the standard GPS ephemeris format for the current time. In this case, the ephemeris provided to the device would be valid for the next 6 hours. The device could request orbit information for the next 12 hours in the standard GPS format which could be supplied as two six hour ephemeris orbit models. In addition, different orbit models and formats that support different accuracies and standards can be generated from the LTO.
- Fitting the LTO to the desired orbit model can be accomplished in a number of mathematical methods. The preferred embodiment is a least-squares fit of the orbit model parameters to the trajectory data. Other methods, such as Kalman filters or other estimators can also be used to obtain the orbit model parameters that best fit the trajectory data. These techniques of fitting data to orbit models are well known to people skilled in the art of orbit determination and orbit modeling.
- The least squares technique provides an optimal fit of the trajectory data to the orbit model parameters.
FIG. 4 depicts a flow diagram of a method of generating an orbit model using a least squares estimation technique. One embodiment of LTO is a table representation of time, position, and clock offset for each satellite, as shown inFIG. 6 . The time, position, and clock offset can be in any time/coordinate system. For the purpose of simplicity and illustration, the time/coordinate system is GPS time and Earth-Centered-Earth-Fixed (ECEF) position in the World Geodetic Survey 1984 (WGS-84) reference frame. - At
step 402, the STD for the desired time interval is extracted from the STD database. The orbit model parameters are initialized to the orbit model values obtained by a similar process for the previous interval. This guarantees that the initial orbit model parameters are a good fit at least for the beginning of the desired time interval. The rest of theprocess 400 will ensure that the parameters are adjusted so that they become a good fit for the entire time interval. - In the preferred embodiment there are 15 orbital parameters to be adjusted:
-
- Square root of semi-major axis (meterŝ½)
- Eccentricity (dimensionless)
- Amplitude of sine harmonic correction term to the orbit radius (meters)
- Amplitude of cosine harmonic correction term to the orbit radius (meters)
- Mean motion difference from computed value (radians/sec)
- Mean anomaly at reference time (radians)
- Amplitude of cosine harmonic correction term to the argument of latitude (radians)
- Amplitude of sine harmonic correction term to the argument of latitude (radians)
- Amplitude of cosine harmonic correction term to the angle of inclination (radians)
- Amplitude of sine harmonic correction term to the angle of inclination (radians)
- Longitude of ascending node of orbit plane at weekly epoch (radians)
- Inclination angle at reference time (radians)
- Rate of inclination angle (radians/sec)
- Argument of perigee (radians)
- Rate of right ascension (radians/sec)
- Although it will be readily apparent that more terms may be used, for better fits, or, fewer terms may be used for a more compact model.
- At
step 404, the orbit model is used to predict what the trajectory would be, the predicted data is denoted the “Model Trajectory Data” (MTD). If the model were perfect, the MTD would coincide exactly with the STD. Atstep 406, the MTD and STD are compared to see how closely the orbit model fits the orbit data. In the preferred embodiment, thecomparison step 406 is performed by summing the squares of the differences between each trajectory point in the STD and the corresponding point in the MTD, and comparing the resulting sum to a threshold. If the fit is “good”, the model parameters are deemed “good” and the process stops atstep 410. If the fit is not good then the model parameters are adjusted atstep 408. There are many techniques well known in the art for adjusting model parameters to fit data. For example, inFIG. 5 , the six-hour ephemeris model was adjusted to fit six hours of STD using a subspace trust region method based on the interior-reflective Newton method described in Coleman, T. F., and Y. Li, “On the convergence of reflective Newton methods for large scale nonlinear minimization subject to bounds”, Mathematical Programming, Vol. 67,Number 2, pp. 189-224, 1994, and Coleman, T. F., and Y. Li, “An interior, trust region approach for nonlinear minimization subject to bounds”, SIAM Journal on Optimization, Vol. 6, pp. 418-445, 1996. There are standard computer packages, e.g., MATLAB Optimization Toolbox, which may be used to implement these methods. -
Steps - When fitting an orbit model to trajectory data, there are many choices of which orbit model to choose. The preferred embodiment is to use orbit models with parameters that have been defined in well-known standards. In one embodiment, the ephemeris parameters defined in the GPS interface control document, ICD-GPS-200c, are used. The ICD-GPS-200c definition includes a bit that specifies a 4-hour fit or a 6-hour fit. Typically, the satellite data is broadcast in 4-hour fits and, by the time this data is obtained by the observer of the satellite, the data is often near the end of its fit interval. In one embodiment of the current invention, sequential 6 hour windows of STD are used to create 6-hour ephemeris models, using the technique described in
FIG. 4 and the accompanying text. This produces a set of ephemeris models as illustrated inFIG. 3 . Although these particular 6-hour models are not available without this invention, the models nonetheless are defined using standard parameters (i.e. ICD-GPS-200c) and will be understood by any device that was designed to be compatible with said standard. -
FIG. 5 shows an example of Satellite Tracking Data (STD) that was generated for a time interval of greater than six hours. Then, using the technique described byFIG. 4 and accompanying text, parameters of an ICD-GPS-200c ephemeris model were adjusted to give a best fit to 6 hours of the STD. The orbit modeled by this 6-hour ephemeris was then compared to the true trajectory and, for comparison; the true trajectory was also compared to the orbit modeled by the broadcast ephemeris. The results are shown ingraph 500 ofFIG. 5 . Thevertical axis 502 represents position error in meters and thehorizontal axis 504 represents time in hours. Thegraph 500 illustrates how thebroadcast ephemeris 506 loses validity while theephemeris 508 created by this invention maintains its validity with approximately one meter of error. - The clock offset of GPS satellites is easily modeled by three parameters. In the preferred embodiment, the measured clock offset is modeled by the three parameters defined in ICD-GPS-200c. These parameters represent clock offset, drift, and drift rate. The parameters are adjusted in a similar way to the
method 400 described above to give a model that best fits the measured data over the time interval. - Alternative embodiments may use longer fit intervals, such as 8, 14, 26, 50, 74, 98, 122, or 146 hours for each ephemeris model. These fit intervals are envisaged in ICD-GPS-200c, but are seldom, if ever, available from the broadcast ephemeris. Under the current invention, models with these fit intervals may be generated even when the broadcast ephemeris is limited to a 4-hour fit interval.
- Alternative embodiments of the STD data may include observed satellite velocity, acceleration, clock drift, or clock drift rate and these terms may be used in the process of fitting a model in ways which are well known in the art.
- Another embodiment of an orbit model uses the spare data bits in the current ephemeris format of a conventional GPS signal to provide additional model parameters that would improve the data fit over long time intervals. For example,
subframe 1 has 87 spare bits that are available for additional parameters. This technique allows for more parameters to describe the orbital motion of the satellites without compromising the standard data format. This new ephemeris model is based on the current ephemeris model with additional correction terms used to augment the model to support the longer fit intervals with greater accuracy. - Yet another embodiment of an orbit model is to develop a new set of orbital parameters that describe the satellite orbit which are different, in part or in their entirety, from the GPS ephemeris model parameters. With the goal of making the fit interval longer, different parameters may provide a better description of the satellite orbit. This new set of parameters could be defined such that they would fit into the existing data structures, however, their implementation and algorithms for use would be different.
- Still a further embodiment of an orbit model would be to develop a new set of orbital parameters that would not fit into the existing GPS ephemeris model format. This new set of parameters would be developed to better address the trade-off between the number of parameters required, the fit interval, and the orbit accuracy resulting from the model. An example of this type of ephemeris parameter set is Brouwer's theory that could be used as is or modified to account for GPS specific terms. Brouwer's theory as described in Brouwer, D. “Solution of the Problem of Artificial Satellite Theory without Drag”, Astron J. 64: 378-397, November 1959 is limited to satellites in nearly circular orbits such as GPS satellites.
- Another embodiment is to use a subset of the standard ephemeris parameters defined in ICD-GPS-200c. This approach is particularly useful when bandwidth and/or packet size is limited in the communication link that will be used to convey the orbit model to the Remote GPS Receiver. In one such embodiment, the fifteen orbit parameters described above, and in ICD-GPS-200c, may be reduced to a subset of 9 parameters, by setting all harmonic terms in the model to zero:
-
- Square root of semi-major axis (meterŝ½)
- Eccentricity (dimensionless)
- Mean motion difference from computed value (radians/sec)
- Mean anomaly at reference time (radians)
- Longitude of ascending node of orbit plane at weekly epoch (radians)
- Inclination angle at reference time (radians)
- Rate of inclination angle (radians/sec)
- Argument of perigee (radians)
- Rate of right ascension (radians/sec)
-
Process 400 is then executed using this subset of parameters. This reduces the amount of data that must be sent to the Remote GPS Receiver. The receiver can then reconstruct a standard ephemeris model by setting the “missing” harmonic terms to zero. There are a large number of alternative embodiments to reduce the size of the data, while still providing a model that fits the STD, including: -
- Removing parameters from the model, and replacing them with a constant, such as zero—as done above—or some other predetermined value, which is either stored in the Remote GPS Receiver, or occasionally sent to the receiver.
- The resolution of the parameters may be restricted in the
process 400; this too reduces the amount of data that must be sent to the mobile GPS receiver. - Parameters, which are similar among two or more satellites, may be represented as a master value plus a delta, where the delta requires fewer bits to encode; an example of this is the parameter Eccentricity, which changes very little among different GPS satellites.
- Some of these approaches reduce the ability of the model to fit the data over a period of time (e.g., six hours). In this case, the fit interval may be reduced (e.g. to four hours) to compensate.
- Each satellite in the plurality of GNSS satellites has an associated LTO. Together all the LTOs for a plurality of satellites form an LTO set. The present invention compresses the LTO set by reducing information redundancy from the set. The LTO set becomes a compressed set comprising an ephemeris structure (denoted lto—0) plus a sequence of “deltas” to the ephemeris structure. To decompress the compressed set, client software in the GNSS receiver extracts the ephemeris structure of the
lto —0 and adds a corresponding delta tolto —0 to form each LTO in the set. -
FIG. 7 depicts a flow diagram of one embodiment of acompression method 700. Themethod 700 begins atstep 702 and proceeds to step 704, where themethod 700 being executed by the server (106 ofFIG. 1 ) accesses a set of LTO.Method 700 applies to each satellite separately. Atstep 706, themethod 700 designates one LTO as a reference LTO (lto—0). In the one embodiment,lto —0 is the first LTO in the set of LTO for this particular satellite. In alternative embodiments,lto —0 could be any LTO in the set of LTO, or it could be another structure with the same fields, such as the current satellite ephemeris, or almanac, In anotheralternative embodiment lto —0 could derived by using information from the set of LTO, (such as forming lto—0 from the average LTO, the median LTO, or some other derived LTO). In yet anotheralternative embodiment lto —0 could be arbitrarily assigned. Atstep 708, variable “i” is set to 1 and, atstep 710, the method selects the next LTO (lto_i). The difference (Δlto_i) on a field-by-field basis between the reference LTO (lto—0) and the current selected LTO (lto_i) is computed atstep 712. As such, fields that do not change from LTO to LTO have a value of zero to indicate no change. In this manner, the LTO of each satellite is compressed. - At
step 714, the Δlto_i is saved. Atstep 716, themethod 700 queries whether themethod 700 has completed compressing all the satellite LTOs. If the query is negatively answered, themethod 700 proceeds to step 720, where “i” is increased by 1 to enablestep 710 to select the next LTO for compression. - If the
method 700 has completed the compression process for all LTOs, then the method continues withstep 722. Atstep 722, themethod 700 arrangeslto —0 and the Δlto_i values into a data set. Themethod 700 also creates an overhead record—a dynamic range map—that is added to the compressed data set. The dynamic range map contains a dynamic range value (the number of bits) for each field in an LTO. For example, if a field contains no change, then the field in that Δlto_i is assigned a dynamic range of zero and no information needs to be included in that particular field in the Δlto_i. Once the compressed data set is complete, themethod 700 ends atstep 724. It is understood that the data set could be implemented in different ways, for example: as a file, or a series of communications packets, or some other means of transferring the information. -
FIG. 8 depicts a graphical illustration of acompressed data set 800 comprising an uncompressed reference LTO (lto—0) 802, and a series of delta LTOs (Δlto_i) 806. In various alternative embodiments, at least one of adynamic range map 804 and abitmask 808 may also form a portion of the C-LTO data set 800. Note that since TOC (16 bits) changes with every LTO, it is identified in the dynamic range map as having a range of 16 bits and all 16 bits of TOC information appear in the Δlto_i. However the af2 field does not change and has a range of zero and therefore does not appear at all in the Δlto_i. Such compression to reduce redundancy is performed on all the fields. Once complete, the compressed data set can be distributed to the GNSS receiver. - In an alternative embodiment of that shown in
FIG. 8 , thedynamic range map 804 is provided in advance or agreed upon in a specification. In this case, there is no need to include thedynamic range map 804 as part of the LTO data being distributed because it is known. As such, the number of bits of each delta LTO field is fixed. This technique has the advantage of simplicity and ease of encoding in a specification but it is less efficient because it has to use more bits per delta field to cover all possible variations of each parameter in the LTO data. Using this technique raises the risk that an instance occurs where one or more of the fields of the LTO deltas fall outside the fixed dynamic range that has been specified or agreed upon. A further refinement of this technique uses abitmask 808 with one bit for each field indicating whether the compressed format has been used for each particular field. This technique provides a fallback mechanism for the case where one or more parameters exceed the specified dynamic range. - Each of the parameters of the LTO has an associated scale factor. For example, if the LTO comprises standard ephemeris parameters, one such parameter is i0 (orbit inclination angle) which has units of semi-circles and a scale factor of 2−31. In an alternative embodiment, the scale factors could be changed so that fewer bits are used to carry the information.
- In an alternative embodiment, certain terms may be eliminated from the set of delta-LTOs, and reconstructed. For example, the four terms: M0, i0, Ω0 and af0 do not have to be sent more than once per satellite per LTO, since, for each subsequent 6 hour segment, these terms can be reconstructed from the previously distributed values.
-
- M0=Mean anomaly, the angular distance along the orbital plane at which the satellite is found at toe.
- i0=Orbital inclination, the angle of the orbital plane to the equator
- Ω0=Right ascension at weekly epoch, the longitude of the intersection of the orbital plane with the equator, at the start of the week.
- af0=clock offset
- Each of these terms can be computed from previous values as follows. Given a set of orbit parameters defined at some reference time toe=toc, the values of Ω0, af0, M0 and i0 can be computed for some other reference time (toe+dt) as follows:
-
- Ω0=Ω0+Ω_dot*dt;
- af0=af0+af1*dt+½af2*dt̂2
- M0=Mo+dt*√{square root over ((μ/(α3)))}
- i0=i0+iDot*dt
where: - μ=3.986005e14; %WGS-84 Universal gravitational parameter (m3/s2);
- We=7.2921151467e-5; %WGS-84 Earth Rotation Rate (rad/sec);
- a=semi-major axis of the orbit;
- af1=clock rate, af2=clock acceleration; and
- the new values are defined with respect to the other time t=toe+dt
- To recover Ω0, af0, M0 and i0 at
n 6 hour epochs after toe, the method propagates forward n times using the above equations, and computes the appropriate values on the right hand side of the equations at each step. -
FIG. 9 depicts a flow diagram of amethod 900 that is performed by the GNSS receiver to unpack the compressed data set and recreate LTO. Themethod 900 starts atstep 902 and proceeds to step 904, wherein the compressed data set is received as described above with respect to a decompressed LTO. Atstep 906, the reference LTO (lto—0) is extracted from the data set and, atstep 908, the variable “i” is set to 1. Atstep 910, themethod 900 selects a Δlto_i and, atstep 912, addslto —0 to Δlto_i on a field by field basis. The dynamic range map is used to identify where the data is located within Δlto_i and/or to identify bit sizes of each parameter with Δlto_i. Atstep 914, the value of lto_i is saved in memory. - At
step 916, themethod 900 queries whether all the Δlto_i have been processed. If the query is negatively answered, themethod 900 proceeds to step 922 to increase i by 1 and return to step 910 to select the next Δlto_i for decompression. If, atstep 916, thedecompression method 900 has processed all Δlto_i, then themethod 900 continues to step 918. Atstep 918, themethod 900 arranges the stored LTO values (lto —0 and all lto_i) into an LTO set that is then used by the GNSS receiver in the manner described above with respect to uncompressed LTO. - While the foregoing is directed to the preferred embodiment of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
Claims (32)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/927,236 US20080189037A1 (en) | 2001-06-06 | 2007-10-29 | Method and apparatus for decompression of compressed long-term-orbit data |
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/875,809 US6542820B2 (en) | 2001-06-06 | 2001-06-06 | Method and apparatus for generating and distributing satellite tracking information |
US09/884,874 US6560534B2 (en) | 2001-06-06 | 2001-06-19 | Method and apparatus for distributing satellite tracking information |
US09/993,335 US7053824B2 (en) | 2001-11-06 | 2001-11-06 | Method and apparatus for receiving a global positioning system signal using a cellular acquisition signal |
US11/333,787 US7443340B2 (en) | 2001-06-06 | 2006-01-17 | Method and apparatus for generating and distributing satellite tracking information |
US11/724,627 US8090536B2 (en) | 2001-06-06 | 2007-03-15 | Method and apparatus for compression of long term orbit data |
US11/927,236 US20080189037A1 (en) | 2001-06-06 | 2007-10-29 | Method and apparatus for decompression of compressed long-term-orbit data |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/724,627 Continuation US8090536B2 (en) | 2001-06-06 | 2007-03-15 | Method and apparatus for compression of long term orbit data |
Publications (1)
Publication Number | Publication Date |
---|---|
US20080189037A1 true US20080189037A1 (en) | 2008-08-07 |
Family
ID=46206133
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/724,627 Expired - Fee Related US8090536B2 (en) | 2001-06-06 | 2007-03-15 | Method and apparatus for compression of long term orbit data |
US11/927,236 Abandoned US20080189037A1 (en) | 2001-06-06 | 2007-10-29 | Method and apparatus for decompression of compressed long-term-orbit data |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/724,627 Expired - Fee Related US8090536B2 (en) | 2001-06-06 | 2007-03-15 | Method and apparatus for compression of long term orbit data |
Country Status (1)
Country | Link |
---|---|
US (2) | US8090536B2 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100052980A1 (en) * | 2008-09-02 | 2010-03-04 | Seiko Epson Corporation | Data compressing method, data providing method, and data compressing device |
US20110012779A1 (en) * | 2009-07-16 | 2011-01-20 | Charles Abraham | Method and system for a virtual wide area gnss reference network |
US20120119945A1 (en) * | 2008-07-25 | 2012-05-17 | Seiko Epson Corporation | Satellite orbital data compression method, satellite orbital data providing method, satellite orbital data expansion method, server, and positioning apparatus |
USRE48176E1 (en) | 2001-07-25 | 2020-08-25 | Avago Technologies International Sales Pte. Limited | Method and apparatus for generating and distributing satellite tracking information in a compact format |
US11428822B2 (en) | 2016-12-01 | 2022-08-30 | Google Llc | Methods and systems for location determination |
US12061275B2 (en) | 2019-07-08 | 2024-08-13 | Google Llc | Enhancing sensitivity to reflected GNSS signals |
US12140686B2 (en) | 2019-03-01 | 2024-11-12 | Google Llc | Determining velocity using a reflected positioning signal |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7839324B2 (en) * | 2007-02-12 | 2010-11-23 | Sirf Technology, Inc. | Efficient ephemeris coding |
US9007261B2 (en) | 2010-05-06 | 2015-04-14 | Mediatek Inc. | Method and apparatus for fast TTFF |
US20140266872A1 (en) * | 2013-03-15 | 2014-09-18 | Hackproof Technologies Inc. | Space Needles |
US20150247928A1 (en) * | 2014-02-28 | 2015-09-03 | Texas Instruments Incorporated | Cooperative location sensor apparatus and system for low complexity geolocation |
JP6476609B2 (en) * | 2014-06-30 | 2019-03-06 | カシオ計算機株式会社 | POSITIONING SYSTEM, POSITIONING DEVICE, STORAGE DEVICE, AND POSITIONING METHOD |
US10747747B2 (en) * | 2014-12-11 | 2020-08-18 | International Business Machines Corporation | Interpreting invalid data as valid data |
CA3045501A1 (en) * | 2016-12-22 | 2018-06-28 | Myriota Pty Ltd | System and method for generating extended satellite ephemeris data |
CN111046049B (en) * | 2019-12-20 | 2021-10-26 | 西南交通大学 | Truck GPS track data compression method |
Citations (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4894662A (en) * | 1982-03-01 | 1990-01-16 | Western Atlas International, Inc. | Method and system for determining position on a moving platform, such as a ship, using signals from GPS satellites |
US5225842A (en) * | 1991-05-09 | 1993-07-06 | Navsys Corporation | Vehicle tracking system employing global positioning system (gps) satellites |
US5365450A (en) * | 1992-12-17 | 1994-11-15 | Stanford Telecommunications, Inc. | Hybrid GPS/data line unit for rapid, precise, and robust position determination |
US5420592A (en) * | 1993-04-05 | 1995-05-30 | Radix Technologies, Inc. | Separated GPS sensor and processing system for remote GPS sensing and centralized ground station processing for remote mobile position and velocity determinations |
US5459473A (en) * | 1990-03-09 | 1995-10-17 | Sigtec Navigation Pty Ltd. | GPS receiver |
US5736960A (en) * | 1995-09-19 | 1998-04-07 | Northrop Grumman Corporation | Atomic clock augmented global positioning system receivers and global positioning system incorporating same |
US5828336A (en) * | 1996-03-29 | 1998-10-27 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Robust real-time wide-area differential GPS navigation |
US5841396A (en) * | 1996-03-08 | 1998-11-24 | Snaptrack, Inc. | GPS receiver utilizing a communication link |
US5874914A (en) * | 1995-10-09 | 1999-02-23 | Snaptrack, Inc. | GPS receiver utilizing a communication link |
US5877723A (en) * | 1997-03-05 | 1999-03-02 | Caterpillar Inc. | System and method for determining an operating point |
US5893044A (en) * | 1997-01-21 | 1999-04-06 | Motorola Inc. | Real time clock apparatus for fast acquisition or GPS signals |
US5945944A (en) * | 1996-03-08 | 1999-08-31 | Snaptrack, Inc. | Method and apparatus for determining time for GPS receivers |
US5957982A (en) * | 1998-08-18 | 1999-09-28 | Trw Inc. | Method and system for space navigation |
US6028537A (en) * | 1996-06-14 | 2000-02-22 | Prince Corporation | Vehicle communication and remote control system |
US6052081A (en) * | 1997-02-03 | 2000-04-18 | Snaptrack, Inc. | Method and apparatus for satellite positioning system based time measurement |
US6084544A (en) * | 1997-12-18 | 2000-07-04 | Ericsson Inc. | Method for determining the location of a GPS receiver using an estimated reference time |
US6133874A (en) * | 1996-03-08 | 2000-10-17 | Snaptrack, Inc. | Method and apparatus for acquiring satellite positioning system signals |
US6191731B1 (en) * | 1999-08-25 | 2001-02-20 | Trimble Navigation Limited | GPS receiver having a fast time to first fix |
US6215442B1 (en) * | 1997-02-03 | 2001-04-10 | Snaptrack, Inc. | Method and apparatus for determining time in a satellite positioning system |
US6285316B1 (en) * | 2000-06-02 | 2001-09-04 | Cellguide Ltd. | Locating a mobile unit using signals from both mobile beacons and stationary beacons |
US6313787B1 (en) * | 1999-11-12 | 2001-11-06 | Motorola, Inc. | Method and apparatus for assisted GPS protocol |
US20020005802A1 (en) * | 2000-05-08 | 2002-01-17 | Bryant Roderick C. | Satellite-based positioning system receiver for weak signal operation |
US6433739B1 (en) * | 1998-03-17 | 2002-08-13 | Qualcomm, Incorporated | Method and apparatus for synchronizing base stations using remote synchronizing stations |
US6438702B1 (en) * | 1999-12-21 | 2002-08-20 | Telcordia Technologies, Inc. | Method for providing a precise network time service |
US20020190898A1 (en) * | 2001-06-14 | 2002-12-19 | Global Locate Inc. | Method and apparatus for creating and distributing satellite orbit and clock data |
US20030023379A1 (en) * | 2001-07-25 | 2003-01-30 | Diggelen Frank Van | Method and apparatus for generating and distributing satellite tracking information in a compact format |
US6538600B1 (en) * | 1998-10-16 | 2003-03-25 | Lucent Technologies Inc. | Wireless assisted GPS using a reference location |
US6542820B2 (en) * | 2001-06-06 | 2003-04-01 | Global Locate, Inc. | Method and apparatus for generating and distributing satellite tracking information |
US6560534B2 (en) * | 2001-06-06 | 2003-05-06 | Global Locate, Inc. | Method and apparatus for distributing satellite tracking information |
US6573799B1 (en) * | 2001-12-12 | 2003-06-03 | Nokia Corporation | Clock system and corresponding method for providing a clock time accounting for systematic error |
US6618670B1 (en) * | 2000-09-15 | 2003-09-09 | Sirf Technology, Inc. | Resolving time ambiguity in GPS using over-determined navigation solution |
US20070103364A1 (en) * | 2004-09-14 | 2007-05-10 | Garin Lionel J | Determining position without use of broadcast ephemeris information |
US7679550B2 (en) * | 2006-02-21 | 2010-03-16 | Garrison James L | System and method for model-base compression of GPS ephemeris |
Family Cites Families (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3133404B2 (en) | 1991-08-20 | 2001-02-05 | ジーイー横河メディカルシステム株式会社 | Ultrasound diagnostic equipment |
US5430657A (en) * | 1992-10-20 | 1995-07-04 | Caterpillar Inc. | Method and apparatus for predicting the position of a satellite in a satellite based navigation system |
US5726893A (en) | 1992-12-17 | 1998-03-10 | Stanford Telecommunications, Inc. | Cellular telephone with voice-in-data modem |
US5506781A (en) | 1994-06-03 | 1996-04-09 | Itt Corporation | RF link control of satellite clocks |
US5537410A (en) * | 1994-09-15 | 1996-07-16 | Oki Telecom | Subsequent frame variable data rate indication method |
US5673256A (en) | 1995-07-25 | 1997-09-30 | Motorola, Inc. | Apparatus and method for sending data messages at an optimum time |
US5884214A (en) | 1996-09-06 | 1999-03-16 | Snaptrack, Inc. | GPS receiver and method for processing GPS signals |
US6272316B1 (en) | 1995-11-17 | 2001-08-07 | Globalstar L.P. | Mobile satellite user information request system and methods |
US5963167A (en) | 1996-03-13 | 1999-10-05 | California Institute Of Technology | Analyzing system for global positioning system and general satellite tracking |
US5862495A (en) * | 1996-09-18 | 1999-01-19 | Lockheed Martin Corp. | Real time position correction to ground generated spacecraft ephemeris |
US6411899B2 (en) | 1996-10-24 | 2002-06-25 | Trimble Navigation Ltd. | Position based personal digital assistant |
US5938721A (en) | 1996-10-24 | 1999-08-17 | Trimble Navigation Limited | Position based personal digital assistant |
US5983109A (en) | 1997-02-03 | 1999-11-09 | Northern Telecom Limited | Method and apparatus for using advanced positioning systems in cellular communications networks |
US6118977A (en) | 1997-09-11 | 2000-09-12 | Lucent Technologies, Inc. | Telecommunications-assisted satellite positioning system |
US6107960A (en) * | 1998-01-20 | 2000-08-22 | Snaptrack, Inc. | Reducing cross-interference in a combined GPS receiver and communication system |
AU2866899A (en) * | 1998-02-06 | 1999-08-23 | Government Of The United States Of America, As Represented By The Secretary Of The Navy, The | Orbit/covariance estimation and analysis (ocean) determination for satellites |
US20030189136A1 (en) * | 1998-05-20 | 2003-10-09 | Toshihide Maeda | Communication system, communication receiving device and communication terminal in the system |
US6058338A (en) | 1999-02-12 | 2000-05-02 | Qualcomm Incorporated | Method and apparatus for efficient GPS assistance in a communication system |
US6411892B1 (en) | 2000-07-13 | 2002-06-25 | Global Locate, Inc. | Method and apparatus for locating mobile receivers using a wide area reference network for propagating ephemeris |
US6453237B1 (en) | 1999-04-23 | 2002-09-17 | Global Locate, Inc. | Method and apparatus for locating and providing services to mobile devices |
US6211819B1 (en) * | 1999-08-27 | 2001-04-03 | Motorola, Inc. | Mobile station location determination in a radio communication system |
US6725012B1 (en) * | 2000-02-01 | 2004-04-20 | The Aerospace Corporation | Method for deploying an orbiting sparse array antenna |
US6429811B1 (en) * | 2000-02-15 | 2002-08-06 | Motorola, Inc. | Method and apparatus for compressing GPS satellite broadcast message information |
US6775802B2 (en) | 2000-08-08 | 2004-08-10 | Qualcomm Incorporated | Method, apparatus, and system for signal prediction |
US6437734B1 (en) | 2000-10-11 | 2002-08-20 | Seiko Epson Corporation | Satellite navigation receiver and method |
US6799116B2 (en) | 2000-12-15 | 2004-09-28 | Trimble Navigation Limited | GPS correction methods, apparatus and signals |
WO2002075349A2 (en) | 2001-03-15 | 2002-09-26 | Qualcomm Incorporated | Time acquisition in a wireless position determination system |
US6856282B2 (en) | 2002-02-08 | 2005-02-15 | Qualcomm Incorporated | Directly acquiring precision code GPS signals |
US7623871B2 (en) | 2002-04-24 | 2009-11-24 | Qualcomm Incorporated | Position determination for a wireless terminal in a hybrid position determination system |
US20030212821A1 (en) | 2002-05-13 | 2003-11-13 | Kiyon, Inc. | System and method for routing packets in a wired or wireless network |
US7349390B2 (en) | 2002-05-28 | 2008-03-25 | Ntt Docomo, Inc. | Packet transmission method and communication system |
-
2007
- 2007-03-15 US US11/724,627 patent/US8090536B2/en not_active Expired - Fee Related
- 2007-10-29 US US11/927,236 patent/US20080189037A1/en not_active Abandoned
Patent Citations (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4894662A (en) * | 1982-03-01 | 1990-01-16 | Western Atlas International, Inc. | Method and system for determining position on a moving platform, such as a ship, using signals from GPS satellites |
US5459473A (en) * | 1990-03-09 | 1995-10-17 | Sigtec Navigation Pty Ltd. | GPS receiver |
US5225842A (en) * | 1991-05-09 | 1993-07-06 | Navsys Corporation | Vehicle tracking system employing global positioning system (gps) satellites |
US5365450A (en) * | 1992-12-17 | 1994-11-15 | Stanford Telecommunications, Inc. | Hybrid GPS/data line unit for rapid, precise, and robust position determination |
US5420592A (en) * | 1993-04-05 | 1995-05-30 | Radix Technologies, Inc. | Separated GPS sensor and processing system for remote GPS sensing and centralized ground station processing for remote mobile position and velocity determinations |
US5736960A (en) * | 1995-09-19 | 1998-04-07 | Northrop Grumman Corporation | Atomic clock augmented global positioning system receivers and global positioning system incorporating same |
US5874914A (en) * | 1995-10-09 | 1999-02-23 | Snaptrack, Inc. | GPS receiver utilizing a communication link |
US6064336A (en) * | 1995-10-09 | 2000-05-16 | Snaptrack, Inc. | GPS receiver utilizing a communication link |
US5841396A (en) * | 1996-03-08 | 1998-11-24 | Snaptrack, Inc. | GPS receiver utilizing a communication link |
US6133874A (en) * | 1996-03-08 | 2000-10-17 | Snaptrack, Inc. | Method and apparatus for acquiring satellite positioning system signals |
US5945944A (en) * | 1996-03-08 | 1999-08-31 | Snaptrack, Inc. | Method and apparatus for determining time for GPS receivers |
US6150980A (en) * | 1996-03-08 | 2000-11-21 | Snaptrack, Inc. | Method and apparatus for determining time for GPS receivers |
US5828336A (en) * | 1996-03-29 | 1998-10-27 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Robust real-time wide-area differential GPS navigation |
US6028537A (en) * | 1996-06-14 | 2000-02-22 | Prince Corporation | Vehicle communication and remote control system |
US5893044A (en) * | 1997-01-21 | 1999-04-06 | Motorola Inc. | Real time clock apparatus for fast acquisition or GPS signals |
US6215442B1 (en) * | 1997-02-03 | 2001-04-10 | Snaptrack, Inc. | Method and apparatus for determining time in a satellite positioning system |
US6052081A (en) * | 1997-02-03 | 2000-04-18 | Snaptrack, Inc. | Method and apparatus for satellite positioning system based time measurement |
US5877723A (en) * | 1997-03-05 | 1999-03-02 | Caterpillar Inc. | System and method for determining an operating point |
US6084544A (en) * | 1997-12-18 | 2000-07-04 | Ericsson Inc. | Method for determining the location of a GPS receiver using an estimated reference time |
US6433739B1 (en) * | 1998-03-17 | 2002-08-13 | Qualcomm, Incorporated | Method and apparatus for synchronizing base stations using remote synchronizing stations |
US5957982A (en) * | 1998-08-18 | 1999-09-28 | Trw Inc. | Method and system for space navigation |
US6538600B1 (en) * | 1998-10-16 | 2003-03-25 | Lucent Technologies Inc. | Wireless assisted GPS using a reference location |
US6191731B1 (en) * | 1999-08-25 | 2001-02-20 | Trimble Navigation Limited | GPS receiver having a fast time to first fix |
US6313787B1 (en) * | 1999-11-12 | 2001-11-06 | Motorola, Inc. | Method and apparatus for assisted GPS protocol |
US6438702B1 (en) * | 1999-12-21 | 2002-08-20 | Telcordia Technologies, Inc. | Method for providing a precise network time service |
US20020005802A1 (en) * | 2000-05-08 | 2002-01-17 | Bryant Roderick C. | Satellite-based positioning system receiver for weak signal operation |
US6285316B1 (en) * | 2000-06-02 | 2001-09-04 | Cellguide Ltd. | Locating a mobile unit using signals from both mobile beacons and stationary beacons |
US6618670B1 (en) * | 2000-09-15 | 2003-09-09 | Sirf Technology, Inc. | Resolving time ambiguity in GPS using over-determined navigation solution |
US6560534B2 (en) * | 2001-06-06 | 2003-05-06 | Global Locate, Inc. | Method and apparatus for distributing satellite tracking information |
US6542820B2 (en) * | 2001-06-06 | 2003-04-01 | Global Locate, Inc. | Method and apparatus for generating and distributing satellite tracking information |
US20020190898A1 (en) * | 2001-06-14 | 2002-12-19 | Global Locate Inc. | Method and apparatus for creating and distributing satellite orbit and clock data |
US20030023379A1 (en) * | 2001-07-25 | 2003-01-30 | Diggelen Frank Van | Method and apparatus for generating and distributing satellite tracking information in a compact format |
US6829535B2 (en) * | 2001-07-25 | 2004-12-07 | Global Locate, Inc. | Method and apparatus for generating satellite tracking information in a compact format |
US6573799B1 (en) * | 2001-12-12 | 2003-06-03 | Nokia Corporation | Clock system and corresponding method for providing a clock time accounting for systematic error |
US20070103364A1 (en) * | 2004-09-14 | 2007-05-10 | Garin Lionel J | Determining position without use of broadcast ephemeris information |
US7679550B2 (en) * | 2006-02-21 | 2010-03-16 | Garrison James L | System and method for model-base compression of GPS ephemeris |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
USRE48176E1 (en) | 2001-07-25 | 2020-08-25 | Avago Technologies International Sales Pte. Limited | Method and apparatus for generating and distributing satellite tracking information in a compact format |
US20120119945A1 (en) * | 2008-07-25 | 2012-05-17 | Seiko Epson Corporation | Satellite orbital data compression method, satellite orbital data providing method, satellite orbital data expansion method, server, and positioning apparatus |
US20100052980A1 (en) * | 2008-09-02 | 2010-03-04 | Seiko Epson Corporation | Data compressing method, data providing method, and data compressing device |
US20110012779A1 (en) * | 2009-07-16 | 2011-01-20 | Charles Abraham | Method and system for a virtual wide area gnss reference network |
US8085194B2 (en) * | 2009-07-16 | 2011-12-27 | Broadcom Corporation | Method and system for a virtual wide area GNSS reference network |
US8339314B2 (en) | 2009-07-16 | 2012-12-25 | Broadcom Corporation | Method and system for a virtual wide area GNSS reference network |
US8638258B2 (en) | 2009-07-16 | 2014-01-28 | Broadcom Corporation | Method and system for a virtual wide area GNSS reference network |
US11428822B2 (en) | 2016-12-01 | 2022-08-30 | Google Llc | Methods and systems for location determination |
US11808863B2 (en) | 2016-12-01 | 2023-11-07 | Google Llc | Methods and systems for location determination |
US12140686B2 (en) | 2019-03-01 | 2024-11-12 | Google Llc | Determining velocity using a reflected positioning signal |
US12061275B2 (en) | 2019-07-08 | 2024-08-13 | Google Llc | Enhancing sensitivity to reflected GNSS signals |
Also Published As
Publication number | Publication date |
---|---|
US8090536B2 (en) | 2012-01-03 |
US20070282910A1 (en) | 2007-12-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8090536B2 (en) | Method and apparatus for compression of long term orbit data | |
US6542820B2 (en) | Method and apparatus for generating and distributing satellite tracking information | |
US9274225B2 (en) | Method and apparatus for generating and distributing satellite tracking information | |
US6560534B2 (en) | Method and apparatus for distributing satellite tracking information | |
USRE48176E1 (en) | Method and apparatus for generating and distributing satellite tracking information in a compact format | |
US20020190898A1 (en) | Method and apparatus for creating and distributing satellite orbit and clock data | |
EP2350681B1 (en) | Providing ephemeris data and clock corrections to a satellite navigation system receiver | |
US6285315B1 (en) | Positioning systems | |
US8362950B2 (en) | Method and device for predicting GNSS satellite trajectory extension data in mobile apparatus | |
US8134500B2 (en) | Method and apparatus for enhanced autonomous GPS | |
US7839324B2 (en) | Efficient ephemeris coding | |
US20110032147A1 (en) | Providing ephemeris data and clock corrections to a satellite navigation system receiver |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: GLOBAL LOCATE, INC., CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:VAN DIGGELEN, FRANK;RIBEN, MATT;DE SALES, JAVIER;REEL/FRAME:021114/0632;SIGNING DATES FROM 20080410 TO 20080418 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |
|
AS | Assignment |
Owner name: BROADCOM CORPORATION, CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:GLOBAL LOCATE, INC.;REEL/FRAME:036617/0654 Effective date: 20150908 |
|
AS | Assignment |
Owner name: BANK OF AMERICA, N.A., AS COLLATERAL AGENT, NORTH CAROLINA Free format text: PATENT SECURITY AGREEMENT;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:037806/0001 Effective date: 20160201 Owner name: BANK OF AMERICA, N.A., AS COLLATERAL AGENT, NORTH Free format text: PATENT SECURITY AGREEMENT;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:037806/0001 Effective date: 20160201 |
|
AS | Assignment |
Owner name: AVAGO TECHNOLOGIES GENERAL IP (SINGAPORE) PTE. LTD., SINGAPORE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:041706/0001 Effective date: 20170120 Owner name: AVAGO TECHNOLOGIES GENERAL IP (SINGAPORE) PTE. LTD Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:041706/0001 Effective date: 20170120 |
|
AS | Assignment |
Owner name: BROADCOM CORPORATION, CALIFORNIA Free format text: TERMINATION AND RELEASE OF SECURITY INTEREST IN PATENTS;ASSIGNOR:BANK OF AMERICA, N.A., AS COLLATERAL AGENT;REEL/FRAME:041712/0001 Effective date: 20170119 |