Multi-antenna GNSS positioning method and system

Information

  • Patent Grant
  • 8686900
  • Patent Number
    8,686,900
  • Date Filed
    Thursday, January 8, 2009
    15 years ago
  • Date Issued
    Tuesday, April 1, 2014
    10 years ago
Abstract
A multi-antenna GNSS system and method provide earth-referenced GNSS heading and position solutions. The system and method compensate for partial blocking of the antennas by using a known attitude or orientation of the structure, which can be determined by an orientation device or with GNSS measurements. Multiple receiver units can optionally be provided and can share a common clock signal for processing multiple GNSS signals in unison. The system can optionally be installed on fixed or slow-moving structures, such as dams and marine vessels, and on mobile structures such as terrestrial vehicles and aircraft.
Description
BACKGROUND OF THE INVENTION

The invention relates generally to Global Navigation Satellite System (GNSS) receivers and more particularly to a method and an apparatus for computing multiple precise locations using differential carrier phases of a GNSS satellite signal by synchronizing the clocks between the master receiver and the slave receiver. It further describes a technique of connecting a plurality of antennas to the slave receiver, which can be switched on to measure each antenna's relative location to the master antenna for monitoring long-term deformation.


GNSS Background

GNSS includes the Global Positioning System (GPS), which was established by the United States government and employs a constellation of 24 or more satellites in well-defined orbits at an altitude of approximately 26,500 km. These satellites continually transmit microwave L-band radio signals in two frequency bands, centered at 1575.42 MHz and 1227.6 MHz, denoted as L1 and L2 respectively. These signals include timing patterns relative to the satellite's onboard precision clock (which is kept synchronized by a ground station) as well as a navigation message giving the precise orbital positions of the satellites. GPS receivers process the radio signals, compute ranges to the GPS satellites, and by triangulating these ranges, the GPS receiver determines its position and its internal clock error. Different levels of accuracies can be achieved depending on the techniques deployed. This invention specifically targets the sub-centimeter accuracies achievable on a remote and possibly mobile GPS receiver by processing carrier phase observations both from the remote receiver and from one or more fixed-position reference stations. This procedure is often referred to as Real Time Kinematic or RTK.


To gain a better understanding of the accuracy levels achievable by using the GPS system, it is necessary to understand the two types of signals available from the GPS satellites. The first type of signal includes both the Coarse Acquisition (C/A), which modulates the L1 radio signal and precision (P) code, which modulates both the L1 and L2 radio signals. These are pseudorandom digital codes that provide a known pattern that can be compared to the receiver's version of that pattern. By measuring the time-shift required to align the pseudorandom digital codes, the GPS receiver is able to compute an unambiguous pseudo-range to the satellite. Both the C/A and P codes have a relatively long “wavelength,” of about 300 meters (1 microsecond) and 30 meters (0.1 microsecond), respectively. Consequently, use of the C/A code and the P code yield position data only at a relatively coarse level of resolution.


The second type of signal utilized for position determination is the carrier signal. The term “carrier”, as used herein, refers to the dominant spectral component which remains in the radio signal after the spectral content caused by the modulated pseudorandom digital codes (C/A and P) is removed. The L1 and L2 carrier signals have wavelengths of about 19 and 24 centimeters, respectively. The GPS receiver is able to “track” these carrier signals, and in doing so, make measurements of the carrier phase to a small fraction of a complete wavelength, permitting range measurement to an accuracy of less than a centimeter.


In stand-alone GPS systems that determine a receiver's position coordinates without reference to a nearby reference receiver, the process of position determination is subject to errors from a number of sources. These include errors in the satellite's clock reference, the location of the orbiting satellite, ionospheric refraction errors (which delay GPS code signals but advance GPS carrier signals), and tropospheric induced delay errors. Prior to May 2, 2002, a large portion of the satellite's clock error, referred to as Selective Availability (SA) was purposefully induced by the U.S. Department of Defense to limit GPS accuracy to non-authorized users. SA would often cause positioning errors exceeding 40 meters, but even today, with SA off; errors caused by the ionosphere can be tens of meters. The above mentioned error sources (satellite clock and satellite position errors, ionosphere refraction, tropospheric delay and SA) are common-mode errors for two receivers that are nearby. That is, the errors caused by these sources are nearly the same for each receiver


Another error source, which is present in the carrier phase measurements, is the clock difference between the two receivers. This clock difference applies to all satellite measurements equally, and as such, can be eliminated by what is known as double differencing. This is where one of the satellites is used as a reference and the other satellite measurements are compared to it. This reduces the number of usable satellite measurements by one. As will be explained later, the more measurements available the better the final solution.


To overcome the common-mode errors of the stand-alone GPS system, many kinematic positioning applications make use of multiple GPS receivers. A reference receiver located at a reference site having known coordinates receives the satellite signals simultaneously with the receipt of signals by a remote receiver. Depending on the separation distance, the common-mode errors mentioned above will affect the satellite signals equally for the two receivers. By taking the difference between signals received both at the reference site and at the remote location, common-mode errors are effectively eliminated. This facilitates an accurate determination of the remote receiver's coordinates relative to the reference receiver's coordinates.


The technique of differencing signals is known in the art as differential GPS (DGPS) or differential GNSS (DGNSS). The combination of DGPS with precise measurements of carrier phase leads to position accuracies of less than one centimeter root-mean-squared (centimeter-level positioning). When DGPS/DGNSS positioning utilizing carrier phase is done in real-time while the remote receiver is potentially in motion, it is often referred to as Real-Time Kinematic (RTK) positioning.


One of the difficulties in performing RTK positioning using carrier signals is the existence of an inherent ambiguity that arises because each cycle of the carrier signal looks exactly alike. Therefore, the range measurement based upon carrier phase has an ambiguity equivalent to an integer number of carrier signal wavelengths. Various techniques are used to resolve the ambiguity, which usually involves some form of double-differencing of the carrier measurements. Once ambiguities are solved, however, the receiver continues to apply a constant ambiguity correction to a carrier measurement until loss of lock on that carrier signal or partial loss of lock that results in a carrier cycle slip.


Regardless of the technique deployed, the problem of solving integer ambiguities, in real-time, is always faster and more robust if there are more measurements upon which to discriminate the true integer ambiguities. Robust means that there is less chance of choosing an incorrect set of ambiguities. The degree to which the carrier measurements collectively agree to a common location of the GPS receiver is used as a discriminator in choosing the correct set of ambiguities. The more carrier phase measurements that are available, the more likely it is that the best measure of agreement will correspond to the true (relative to the reference GPS) position of the remote GPS receiver. One method, which effectively gives more measurements, is to use carrier phase measurements on both L1 and L2. The problem though is that it is relatively difficult to track L2 because it is modulated only by P code and the United States Department of Defense has limited access to P code modulation by encrypting the P code prior to transmission. Some receivers are capable of applying various cross-correlation techniques to track the P code on L2, but these are usually more expensive receivers than L1 only capable receivers.


Other approaches have been employed to gain additional measurements on GPS receivers utilizing additional satellites and other types of satellite systems such as the GLONASS system, pseudolites, or Low Earth Orbit (LEO) satellite signals in an attempt to enhance RTK. Nevertheless, it is often desired to perform RTK on low-cost L1 only receivers that do not have access to the GLONASS system, pseudolites, or LEO satellite signals.


SUMMARY OF THE INVENTION

Disclosed herein in an exemplary embodiment is a method for measuring relative position of fixed or slow-moving points in close proximity comprising: receiving a set of satellite signals with a first receiver corresponding to a first position; receiving a related set of satellite signals with a second receiver corresponding to a second position; and computing a position of the second position based on at least one of code phase and carrier phase differencing techniques. At least one of: a clock used in the first receiver and a clock used in the second receiver are synchronized to eliminate substantial clock variation between the first receiver and the second receiver; and the first receiver and the second receiver share a common clock.


Also disclosed herein in another exemplary embodiment is a system for measuring relative position of fixed or slow-moving points in close proximity comprising: a first receiver in operable communication with a first antenna configured to receive a first plurality of satellite signals at a first position; and a second receiver in operable communication with a second antenna configured to receive a second plurality of satellite signals at a second position; and at least one of the first receiver and the second receiver computing a position corresponding to a position of the second antenna based on at least one of code phase and carrier phase differencing techniques. At least one of: a clock used in the first receiver and a clock used in the second receiver are synchronized to eliminate clock variation between the first receiver and the second receiver, and the first receiver and the second receiver share a common clock.


Further, disclosed herein in yet another exemplary embodiment is a system for measuring relative position of fixed or slow-moving points in close proximity comprising: a means for receiving a set of satellite signals with a first receiver corresponding to a first position; a means for receiving a related set of satellite signals with a second receiver corresponding to a second position; and a means for computing a position of the second position based on at least one of code phase and carrier phase differencing techniques. At least one of: a clock used in the first receiver and a clock used in the second receiver are synchronized to eliminate clock variation between the first receiver and the second receiver, and the first receiver and the second receiver share a common clock.


Also disclosed herein in yet another exemplary embodiment is a storage medium encoded with a machine-readable computer program code, the code including instructions for causing a computer to implement the abovementioned method for measuring relative position of fixed or slow-moving points in close proximity.


Further disclosed herein in yet another exemplary embodiment is a computer data signal, the computer data signal comprising code configured to cause a processor to implement the abovementioned method for measuring relative position of fixed or slow-moving points in close proximity. Still further, applications are disclosed for the system and method using L1 receivers and RTK. Additional applications include marine vessels, terrestrial vehicles and aircraft with partially blocked GNSS antennas. Multipath effects are mitigated by averaging the GNSS signals received by multiple antennas or by an antenna(s) and multiple locations.





BRIEF DESCRIPTION OF THE DRAWINGS

Referring now to the drawings wherein like elements are numbered alike in the several figures:



FIG. 1 is a diagram of a system embodying an aspect of the invention and including combined master and slave receivers;



FIG. 1A is a diagram of a system embodying an alternative aspect of the invention and including separate master and slave receivers;



FIG. 2 is a vertical, cross-sectional view of an application of the invention, shown in connection with a dam for monitoring the locations of various points thereon;



FIG. 3 is top plan view of another application of the invention, shown in connection with a marine vessel;



FIG. 4 is a diagram of a real-time kinematic (RTK) system embodying another aspect of the present invention and using single frequency (L1) receivers;



FIGS. 5A-C show constrained, multiple-antenna GNSS positioning systems with satellite-blocking objects;



FIGS. 6 and 6A show another application of the invention in connection with an automobile;



FIG. 7 shows another application of the invention in connection with an aircraft.



FIG. 8 shows another application of the invention with an antenna ring;



FIG. 9 shows another application of the invention with a linked pair of antennas driven by a stepper motor and located in a radome; and



FIG. 10 shows another application of the invention with a continuously-rotating linked pair of antennas.





DETAILED DESCRIPTION OF PREFERRED ASPECTS OF THE INVENTION

As required, detailed embodiments of the present invention are disclosed herein; however, it is to be understood that the disclosed embodiments are merely exemplary of the invention, which may be embodied in various forms. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a basis for the claims and as a representative basis for teaching one skilled in the art to variously employ the present invention in virtually any appropriately detailed structure. Certain terminology will be used in the following description for convenience in reference only and will not be limiting. For example, up, down, front, back, right and left refer to the invention as oriented in the view being referred to. The words “inwardly” and “outwardly” refer to directions toward and away from, respectively, the geometric center of the embodiment being described and designated parts thereof. Global navigation satellite systems (GNSS) are broadly defined to include the global positioning system (GPS, U.S.), Galileo (proposed), GLONASS (Russia), Beidou (China), Compass (proposed), the Indian Regional Navigational Satellite System (IRNSS), QZSS (Japan, proposed) and other current and future positioning technology using signals from satellites, with or without augmentation from terrestrial sources. Inertial navigation systems (INS) include gyroscopic (gyro) sensors, accelerometers and similar technologies for providing output corresponding to the inertia of moving components in all axes, i.e. through six degrees of freedom (positive and negative directions along transverse X, longitudinal Y and vertical Z axes). Yaw, pitch and roll refer to moving component rotation about the Z, X and Y axes respectively. Said terminology will include the words specifically mentioned, derivatives thereof and words of similar meaning.


This aspect of the invention discloses the use of two receivers, which either share the same clock, or have a clock synchronization technique to eliminate the receiver clock errors. The reference receiver (herein called the master) is connected to a single antenna whereas the slave receiver, which is clock synchronized with the master, has a multitude of antennas connected to it, which are switched in and out to take a measurement at each antenna location.


The GPS slave (e.g., rover) receiver computes the location vector from a double or single difference of the GPS rover and reference carrier phases for a plurality of GNSS satellites. As the receivers are either co-located or have a link, the raw measurements from the slave antennas are sent to the master for computation (of course any receiver or even a separate computer could perform this computation). This eliminates the need for a radio link between the master and slave receivers as is required in prior art RTK.


According to a more specific aspect of the present invention, in order to solve the integer ambiguity problem, the master selects the slave antenna to be measured based on the GPS satellite almanac to provide the best geometry (or one of the best) and based on its time slot. The master also has the slave antenna's position stored to provide an immediate calculation of the carrier cycle ambiguity to each satellite. Position calculation then follows conventional RTK GPS practice of using single or double difference equations involving the total phase distance to each satellite to solve the relative location of the slave antenna with respect to the master antenna. As previously described, there is no clock difference between the two receivers (or the clock difference is known and nearly constant) so double differencing may not be required. There may however be a significant delay through the coaxial cable to each slave antenna. This also can be stored and the delay removed to the measurements. A temperature drift may be noticed which will gradually change the delay, but this too can be eliminated by the addition of a thermocouple to determine the ambient temperature around the cable and antennas. By doing this, all satellite measurements may be used in the solution.


Another advantage of eliminating double differencing is that ambiguity search routines will not have to form linear combinations to decorrelate the measurement data. When it is possible to use single differences, they are generally preferred over double difference equations. The double difference cross-correlations are more difficult to deal with mathematically, say in a measurement covariance matrix of a Kalman filter. Single difference equations result in a measurement covariance matrix having zero cross correlation. However, accuracy can be achieved with both approaches.


Referring now to FIG. 1, a simplified block diagram of a GNSS positioning system 10 embodying an aspect of the present invention is depicted. In an exemplary embodiment, a method and system use two receivers, which either share the same clock, or include a clock synchronization technique to eliminate a receiver clock. Further, the reference receiver (hereinafter also called the master or base receiver) 12 is connected to a master antenna 16 (Am), whereas the rover or slave receiver 14, which is clock synchronized with the master receiver 12, has a multitude of antennas 18 (A1-An) connected to it, which are switched in and out to take a measurement at each antenna location. As shown in FIG. 1, the master receiver 12 and the slave receiver 14 are combined in a single receiver unit 11 on a common printed circuit board (PCB), which also includes a central processing unit (CPU) 13, a temperature sensor 15, a clock 17 and an antenna switch control 23. Collocating the receivers 12, 14 facilitates communication between them. It will be appreciated that while an exemplary embodiment is described and illustrated with respect to measuring movement of a dam, dike or beam, the disclosed invention is readily applicable to other applications where fixed or slow moving phenomena are tracked. Such applications may include roadways, bridges, building motion, glacier and iceberg travels and the like. It is also applicable to conventional RTK applications that require relatively short distances between master and slave receivers and where it is desirable to take advantage of a common clock for added robustness and the elimination of a radio for cost and robustness. For example, one application is local surveying or measuring distance at a construction site, or leveling (such as required for foundation placement) at that site.


The master and slave receivers 12, 14 are configured to either share the same clock 17, or include a clock synchronization system (SYNC connection). This technique facilitates elimination of the receiver clock errors. The CPU 13 computes a location vector based on a double or single difference of the GNSS code and/or carrier phases for the master receiver 12, the slave receiver 14 and a plurality of GPS satellites 8. As the master and slave receivers 12 and 14 are either co-located or have a link, the raw measurements from the slave antennas are sent to the CPU 13 for computation (of course any receiver or even a separate computer could perform this computation). This eliminates the need for a radio link between the master and slave receivers 12, 14 as is required in existing RTK applications. Moreover, in another exemplary embodiment, satellite signals from multiple antennas with a known dimensional separation may be combined to achieve receiving an optimal set of satellite 8 signals for a given location. Such an approach will be beneficial for instances when insufficient data is available from a single antenna or a less desirable set of satellite 8 signals is all that is available. In this way, a location may still be computed despite poor satellite 8 geometry, obstructions, and the like.


Advantageously, in an exemplary embodiment, rather than increasing the number of measurements, a reduction in the number of unknowns is achieved by eliminating the clock errors between the reference receiver 12 and the rover 14 (or master and slave). This approach yields an even greater advantage than adding measurements, unless a substantial number of measurements could readily be added. In addition, an exemplary embodiment as disclosed herein significantly improves the ability to calculate the integer ambiguities to each satellite 8. It will be appreciated that because the slave antennas 18 are presumed to move far less than a fraction of a carrier cycle (e.g., 19 cm) between measurements, the positions of each slave antenna 18 location may be stored and then later retrieved as needed to facilitate the immediate calculation of the integer ambiguities.


In order to solve the integer ambiguity problem with current RTK applications, the master receiver 12 selects a particular slave antenna 18 to be measured based on the GPS satellite almanac to provide the best geometry (or one of the best) and based on its time slot. The master receiver 12 also has the slave antenna's position stored (as stated above) to provide an immediate calculation of the carrier cycle ambiguity to each satellite 8. Position calculation then follows RTK GNSS practice of using single or double difference equations involving the total phase distance to each satellite 8 to solve the relative location of slave antenna 18 with respect to the master antenna 16. One such methodology for GNSS positioning employing RTK is taught by Whitehead, U.S. Pat. No. 6,469,663 the contents of which are incorporated by reference herein in their entirety. As previously described, there is no clock difference between the two receivers 12 and 14 (or the clock difference is known and nearly constant) so double differencing may not be required. It will however, be readily appreciated that there may be a significant delay through a coaxial cable 20 to each slave antenna 18. This delay is dependent upon the selected position for each antenna relative to the master (e.g., the length of cable 20 to reach each antenna 18). Advantageously, the delay may readily be measured and stored and the delay mathematically removed to correct the measurements. Moreover, selected antennas 18 may exhibit a temperature drift that may result in a gradual change of the expected delay. However, advantageously, this too may be readily eliminated by the addition of a temperature sensor 15 connected to a thermocouple 22 to determine the ambient temperature around the cable 20 and the antennas 16 and 18. Advantageously, by employing the abovementioned correction and compensation schemes, all satellite 8 measurements may be used to formulate the solution.


Another advantage of eliminating double differencing is that ambiguity search routines will not have to form linear combinations to decorrelate the measurement data. When it is possible to use single differences, they are generally preferred over double difference equations. The double difference cross-correlations are more difficult to deal with mathematically, say in a measurement covariance matrix of a Kalman filter. Single difference equations result in a measurement covariance matrix with zero cross-correlation, which facilitates computation of the ambiguities. The accuracy of both approaches should be substantially similar. However, single differencing is an easier process.


Yet another exemplary embodiment as an enhancement to the abovementioned embodiments uses the capability to take advantage of the slow dynamics of antenna motion by averaging over periods of time, thereby reducing multipath contributions (which are time varying) and poor satellite 8 geometries. In fact, it will be appreciated that the master receiver 12 is constantly tracking the satellites 8 and may further be employed to select the best time of day, e.g., the best constellation (the GNSS satellites 8 orbit in a 12 hour cycle), to perform the measurements based on its knowledge of the slave antennas' 18 positions and the satellites currently visible. Additionally the master receiver 12 may select two separate times of day to provide two independent satellite position constellations for performing the measurements. This would reduce the amount of averaging time required, yet still provide the multipath and poor satellite geometry reduction benefits. Overall, such an approach may be employed to reduce power consumption requirements as the receiver 12 would not have to be averaging continuously for a twelve hour period. Power consumption reduction can be beneficial, especially at remote sites.


Referring once again to FIG. 1, the system 10 is shown configured with a plurality of slave antennas 18 (also denoted as A1, A2 . . . An) connected to the slave receiver 14. Each slave antenna 18 is switched (except the last one which is selected when all switches are connected through to it) with a switch box 24 (also denoted as S1, S2 . . . ). The switch(es) 24 are activated and the antennas 18 selected by an antenna switch controller 23, which can be incorporated on the receiver unit 11. The antenna switch controller 23 can send a tone or some other control signal 30 on the cable 20 to activate a particular desired switch 24 and thereby activate the slave antenna 18 connected thereto. It will be appreciated that in order to provide fault protection, the switch(es) 24 may be designed and configured so that in the event a switch 24 fails, the connection through to the next switch 24 is made. Advantageously, in this way, if one switch 24 should fail, it will still permit measurements on the remaining slave antennas 18. Smart reset circuitry can be employed to insure that the master receiver 12 and the slave receiver 14 will start up at the same instant and therefore the samples will be aligned as well. This approach substantially eliminates any receiver clock biases.


As mentioned previously, phase drift and delay can result from the coaxial cables 20, which may be removed and/or compensated by using a temperature sensor 15 connected to a thermocouple 22 to measure the temperature. A look-up table may be employed by the CPU 13 that has stored (alternatively a simple formula may be used to save memory) phase delay difference versus ambient temperature. An alternative embodiment could use equivalent coaxial cable 20 lengths to all antennas 16, 18 so that any temperature or other loss and drift effects would be matched and therefore cancelled in the single difference calculation.


Normally in order to solve for integer ambiguities from GNSS satellite 8 signals, double differencing is used to bring forth the integer nature of the ambiguities by removing other non-integer sources of error such as clock and atmospheric delays from the measurements. To illustrate, consider four equations describing pseudo-ranges resulting from measurements of carrier phase on receivers denoted m and n for the slave and master, respectively:

φmi=Rmi+τsvi+Ai+Bm+Nmi
φni=Rni+τsvi+Ai+Bn+Nni
φmk=Rmk+τsvk+Ak+Bm+Nmk
φnk=Rnk+τsvk+Ak+Bn+Nnk


Here φmi is the measured pseudorange from rover receiver m to satellite i, φni is the measured pseudorange from reference receiver n to satellite i, φmk is the measured pseudorange from rover receiver m to satellite k, and φnk is the measured pseudorange from reference receiver n to satellite k. Each pseudorange is actually a measure of the summation a number of different physical quantities all of which shall be expressed in units of carrier cycles at L1 (roughly 19 cm).


Specifically, in the first of these equations, the term Rmi is the true geometric range from receiver m to satellite i, τsvi is the clock error of satellite i, Ai is the atmospheric delays, which are associated with satellite i, Bm is the clock error of receiver m, and Nmi is the integer ambiguity in the range measurement from receiver m to satellite i. Similar notation applies to the remaining three equations. For simplicity, these equations do not show noise effects such as errors caused by receiver thermal noise or multipath noise.


Consider first applying the single difference. If the first two equations are differenced:

φmi−φni=Rmi−Rni+Bm−Bn+Nmi−Nni

Similarly, differencing the second two equations yields

φmk−φnk=Rmk−Rnk+Bm−Bn+Nmk−Nnk


The satellite common errors, such as satellite clock, τsvi and atmosphere, Ai (atmosphere is common if we assume relative close proximity of receivers m and n) are removed in the single difference. Furthermore, assume the receivers are clock-synchronized so that Bm=Bn, leaving:

φmi−φni=Rmi−Rni+Nmi−Nni
and
φmk−φnk=Rmk−Rnk+Nmk−Nnk


Since the ambiguities are all integers that can be lumped together into a single term, these equations may be written:

φmi−φni=Rmi<Rni+Nmni
φmk−φnk=Rmk−Rnk+Nmnk
where
Nmni=Nmi−Nni and Nmnk=Nmk−Nnk


This shows that single differencing the carrier phase pseudorange measurements removes common atmospheric errors from the equations while leaving simple combinations of the geometric ranges and integer ambiguities, and clock errors drop out due to the synchronization of the two receivers. For N satellites in common view of the master (reference) and slave (remote) receivers 12 and 14 respectively, there are N such single-difference equations that can be formed without causing mathematical redundancy. Whereas double differencing, to eliminate clock biases in receivers, which are not clock synchronous, results in only N−1 equations. This gives rise to N unknown integer ambiguities that must be solved in addition to the 3 unknown coordinates (X,Y,Z) of the GPS receiver. Note that each geometric range term, for example Rmi is a function only of the receiver's position and the transmitting satellite's position. Specifically:







R
m
i

=




(


Xrecv
m

-

Xsat
i


)

2

+


(


Yrecv
m

-

Ysat
i


)

2

+


(


Zrecv
m

-

Zsat
i


)

2







where Xrecvm, YrecvmZrecvm are the Cartesian coordinates of the receiver m at the time reception of the signal from satellite i, whose coordinates are Xsati, Ysati, Zsati at the time of signal transmission. In the problem at hand, only the rover's position is unknown. Once the ambiguities are determined, only the rover's 3-coordinates of position are unknown and these are easily solved using a mathematical approach such as Least-Squares.


Every time a new slave antenna 18 is selected, the integer ambiguities must be solved. This is a complex process and can be very time consuming if the position is unknown. However, in this instance, it will be appreciated that the movements to be measured are on the order of less than a quarter of a wavelength (5 cm) between measurements. This limitation permits a rapid calculation of the integer ambiguities since the master receiver 12 or the CPU 13 “knows” the satellite's position and the selected antenna's position well enough to directly calculate ambiguities. Such an approach will greatly reduce the time utilized to solve for the integer from up to 10 minutes to a second or less. Cycle slips, which result usually from motion which the receiver failed to track properly and therefore slipped from one ambiguity to another is also greatly reduced due to the very low dynamics of the selected antenna location. An added benefit of the low dynamics is the receiver can integrate the measurements over a long period of time and narrow the carrier tracking loop bandwidth to reduce noise.


As mentioned previously, it should be appreciated that another source of error in applying RTK positioning, especially when solving for integer ambiguities over long baselines, is non-common atmospheric propagation delays on the signals received by the slave (rover) 14 and master (reference) receivers 12. Since differencing cannot eliminate these non-common delays, the next best alternative is to estimate or model their effects. However, in an exemplary RTK embodiment, the slave antennas 18 and the master antenna 16 will most likely be within 5 kilometers of each other and at this distance the atmospheric effects are minimal and may readily be ignored. An optional orientation device 32, such as a compass or some other non-GNSS orientation device, can be affixed to a structure of interest to determine its attitude or orientation and to provide a corresponding signal to the CPU 13 for processing in connection with GNSS ranging data received via the receivers 12, 14. Alternatively, the orientation device 32 can comprise two or more antennas 16, 18 providing input to the master and slave receivers 12, 14 whereby attitude or orientation information can be derived by the processor 13.


A further advantage of this technique is that it permits a carrier phase based solution even when a large portion of the sky, and therefore the visible satellites, are obscured by a wall, dam (FIG. 2) or other structure. This is because, as described above, the receivers 12, 14 will still have one more measurement than previously due to the utilization of single differencing rather than double differencing techniques. In addition, the fixed or very slow moving nature of the problem permits long-term measurements.



FIG. 1A shows a GNSS positioning system 40 comprising an alternative aspect of the present invention with a master receiver unit 42 and a separate slave receiver unit 44, which can be connected by a clock-synchronizing connection (SYNC) 46 of the receivers 12, 14, a clock-sharing connection 48 and a link 50, which can comprise a cable or an RF connection between suitable transmitters and receivers. An optional orientation device 32 can be connected to either or both of the receiver units 42, 44, e.g., to the CPU 13 and/or an optional CPU 52 of the slave unit 42. Optionally, the slave unit 44 can include a clock 54, which can be synchronized with the master receiver unit clock 17. The slave receiver 14 is connected to a slave antenna array, which can comprise a single antenna or a multiple-antenna array as shown in FIG. 1.


Referring now to FIG. 2, in a GNSS dam-monitoring positioning system 60 comprising yet another exemplary embodiment, a technique is employed to utilize and take advantage of the master receiver's 12 knowledge of the satellites' locations in the sky, and a preprogrammed knowledge of the visibility of the sky for selected slave antennas 18. For example, FIG. 2 shows a configuration of satellites 8 and slave antennas 18 whereby the slave antenna A1 receives ranging data transmissions from all four satellites (SAT 1-4), but slave antenna A2 only receives transmissions from satellites SAT 2-4 and slave antenna A 3 only receives transmissions from satellite SAT 1. The master receiver 12 and/or the CPU 13 may then choose the best time, that is, the time with the most satellites visible to the selected slave antenna 18, to perform the measurement at that location. The receiver(s) can then dwell for some time (say one half hour) to integrate and reduce noise, then move on to another slave antenna 18. Moreover, it will be appreciated that the master receiver 12 and/or the CPU 13 may direct that the slave receiver 14 return to the same location after some duration, e.g. a few hours, when another optimal/desirable geometry is available, which is uncorrelated to the first. By taking measurements at two (or more) different times (and geometries), and averaging the two (or more) measurements, multipath and atmospheric induced errors, typically correlated over time, will be reduced. This method will allow monitoring of the face of a dam or berm, or even a valley wall.


Further assumptions may be made of the anticipated motion of the monitoring point at the selected slave antenna 18 to further reduce the number of measurements required. For example, the motion is of a darn is generally horizontally away from the pressure excerpted by the body of water behind it. By performing the calculation only in this direction, a single satellite may be enough to perform a measurement. This is obvious when looking at this equation:







R
m

i
=












(


Xrecv
m

-

Xsat
i


)

2

+


(


Yrecv
m

-

Ysat
i


)

2

+


(


Zrecv
m

-

Zsat
i


)

2







As explained previously the satellite position (Xsat, Ysat and Zsat) are known, and if the receiver assumes there is minimal motion in Y and Z, then there is only one unknown left. Of course, additional satellites are highly desired to reduce noise and errors and to help detect any false or erroneous readings from throwing the solution off.


Another area of concern for running a long length of coaxial cable 20 to the antennas 16, 18, other than phase delay, which was addressed earlier, is attenuation. In yet another exemplary embodiment, the slave antennas 18 may be configured as active antennas, e.g., antennas that include an internal Low Noise Amplifier (LNA). In a receiver design, the noise figure is often important, and comprises a combination of the noise temperature before the first LNA, the LNA noise figure and subsequent losses divided by the LNA gain. Subsequent amplifier gains will reduce following noise temperature (T) contributions by their gain as is shown in the equation below:

Tt=T(preLNA)+T(LNA)+T(lna2)/(CL×.Glna1)+T(lna3)/(CL×Glna1×Glna2)+T(lna4)/(CL×Glna1×Glna2×Glna3) etc.


where: CL refers to cable losses in linear terms, that is −10 dB is 0.1,


Glnan refers to gain of LNAn in linear terms so a gain of 20 dB is 100,


T(LNAn) refers to the noise temperature in Kelvin of stage n.


Noise Figure (F) is related to noise temperature by:

F(dB)=10×LOG((1+T)/Tamb)


Where Tamb refers to the reference temperature, typically 290 K (20 Celsius).


As an example, a typical low loss coaxial cable (RG6 type) has 20 dB (CL=0.01) of attenuation every 100 meters. The noise temperature of the antenna and LNA is 170 K (2 dB noise figure), the gain of the first LNA is 30 dB (or 1000). Subsequent LNA's have the same noise temperature and a gain of 12 dB (15.8). If each antenna is 50 meters apart the losses are −10 dB. After five stages the noise temperature of the system is:

T5=T1+T2/(CL1×.G1)+T3/(CL1×C12×G1×.G2)+T4/(CL1×CL2×CL3×GG2G3)+T5/(CL1×C12×C13×C14×GGGG4)
T5=190+190/100+190/158+190/250+190/395


T5=194 K


F5=2.22 dB


This is compared to the first stage, which would have a noise figure of 2 dB. A GPS receiver such as the master receiver 12, or slave receiver 14 can operate with a noise figure of up to 3.5 dB without suffering significant degradation. As can be seen, additional stages will have diminishing contributions. The total gain will be increasing by only 2 dB each step, so after 1 km, in this example, the maximum gain will be 68 dB, the gain of the first stage is 30 dB, the Automatic Gain Control of the receiver can remove this difference easily. Also after 20 stages (1 km) the total noise temperature in this example would be T(1 km)=194.7 K, an insignificant increase.


Further, in a marine vessel (e.g., barge) 71 positioning system 70 comprising another exemplary embodiment shown in FIG. 3, multiple antennas 18 (A1, A2) could be used to compute a solution of a single point on a rigid body to which they are attached, using known geometry and distances. Such an approach may be employed, for example, when not any one antenna 18 provides enough useful information (satellites 8) to compute a location solution due to obstructions, e.g., a superstructure 72, but the constellation of satellites 8 could provide sufficient positioning data. In other words, the superstructure 72 partially blocks the antennas 16 (A1) and 18 (A2) from views of satellites 8 whereby each antenna 16, 18 receives GNSS positioning signals from some, but not all, of the satellite 8 constellation. It will be appreciated that the antennas 16, 18 are positioned in a predetermined, known relation (i.e. spacing and orientation) for determining attitude comprising yaw and roll, although the primary concern would be yaw because the antennas 16, 18 would normally remain at a relatively constant level in a marine application. The antennas 16, 18 are connected to a receiver unit 74, which can be similar to the receiver units described above.


Advantageously, a position solution employing this approach would not necessarily have to utilize carrier-phase based differencing (it could be code phase). An application might include positioning on a marine vessel 71, such as a barge, where the location of a reference point is needed but there are cranes, towers and/or a superstructure 72 blocking the satellite view so that there is not one optimum GNSS location. However, by placing an antenna 18 on either side of the barge 71, enough satellites 8 could be tracked by the combined antenna 16, 18 arrangement that a solution of the location of some reference point on the barge 71 could still be obtained. Furthermore, on a barge 71, the orientation device 32, such as a compass, could also be used to give orientation, thus removing another unknown from the relative location of two receivers (e.g., 12, 14 in the receiver unit 74) rather than solving a relative location of one receiver with respect to the other by using the combined receivers 12, 14 to produce one non-relative location. The system shown in FIG. 3 can also include an optional base unit 76 for differential corrections, e.g., RTK applications. A position solution microprocessor (CPU) can be provided in the receiver unit 74 for calculating position solutions.



FIG. 4 shows a GNSS positioning system 80 comprising another alternative aspect of the present invention, with a base receiver unit 82 and a rover receiver unit 84, which can be configured similarly to the master and slave receiver units 42, 44 described above. The base and rover receiver units 82, 84 include base and rover GNSS receivers 90, 92, clocks 93, 95 and CPUs 94, 96 respectively. The base and rover receiver units 82, 84 can be connected by a cable link 85 for relatively close operation, such as surveying. Alternatively, the base and rover receiver units 82, 84 can be wirelessly connected via a base RF transmitter 86 and a rover RF receiver 88. An optional orientation device 32 can be connected to the rover CPU 92 for providing orientation and attitude information with respect to a rover vehicle or piece of equipment.


The receivers 90, 92 can comprise L1-only receivers, which tend to be less expensive than the more sophisticated dual frequency (L1/L2) receivers, which are commonly used for RTK fine positioning. The base receiver unit 82 is preferably equipped with a single base antenna 98 and the rover receiver unit 84 is preferably equipped with at least two rover antennas 100. Although the base receiver unit 82 could be equipped with multiple antennas (e.g., for reducing multipath errors) and the rover receiver unit 84 could be equipped with a single antenna, the normal preferred configuration includes multiple rover antennas 100 whereby the attitude (orientation) of the rover can be determined using GNSS ranging techniques. Attitude/orientation of the base is not generally needed for rover positioning.


The rover attitude information facilitates resolving integer ambiguities in GNSS positioning solutions by first solving for locations of the rover antennas 100 with respect to each other (an attitude solution). Next, using the known rover antenna relative locations, and nonrelative ambiguities, the system 80 solves for the global ambiguities using observations taken at each antenna 98, 100. The number of observations is thereby significantly increased over conventional RTK systems. The global ambiguities are the ambiguities of one rover antenna 100 that allow it be located in a global sense, i.e. relative to the base receiver unit 82.


The steps of the GNSS positioning method using the system 80 comprise:

  • 1. Transmit code and carrier phase data from a base station to a multiple antenna rover system as would be done in a conventional RTK system that uses only one rover antenna.
  • 2. At the rover side, determine the relative locations and relative ambiguities of the multiple antennas using an attitude solution that takes advantage of known constraints in geometry or clock as described in U.S. Pat. No. 7,388,539. The attitude solution is usually done much quicker than conventional RTK ambiguity resolution due to the use of geometry constraints and/or a common clock.
  • 3. Optionally store off the attitude solution (locations and ambiguities) in step 2 for later retrieval so that the data can be time-tag matched with the data from the base station. Also store off the current GPS observations (carrier phase) for the same purpose. This step is not necessary, but time tag matching of base and rover data improves results by avoiding extrapolation errors (not so critical with SA off, except for possibly on the WAAS satellites which have been observed to have high phase jitter).
  • 4. Form single or double difference equations and solve for the global ambiguities using knowledge of relative antenna locations/clocks and relative ambiguities.


An example is provided below.


Example Using the Two Antenna Rover System 80


At antenna one (1) of the rover, we can write the equation:

R1=[A]×1−N1,

    • Where R1 is a carrier phase observation vector (single or double difference) at antenna (1), A is a design matrix, X1 is the location vector of antenna 1 (may include clock if single differencing is used), and N1 is an ambiguity vector for antenna (1).
    • Similarly, at antenna two (2) we can write:

      R2=[A]×2−N2,
    • Where R2 is a carrier phase observation vector at antenna (1), A is a design matrix, X2 is the location vector of antenna 2, and N2 is an ambiguity vector for antenna (2).
    • In this example the design matrix A is taken to be the same in both antenna equations. However, this is true only if both antennas see the same satellites. A more general example would use separate A1 and A2 for the two equations.
    • Solving an attitude solution (see U.S. Pat. No. 7,388,539), we find the relative antenna displacement V, and the relative ambiguity M where:

      V=x2−x1
      and
      M=N2−N1
    • Thus, combining the above equations, we have:

      R1=[A]×1−N1
      R2=[A](x1+V)−(N1+M)
    • Re-arranging gives:

      R1=[A]×1−N1
      R2−[A]V+M=[A]×1−N1
    • And, combining into a single vector equation gives:

      R=[A}×1−N
      Where
      R=[R1,R2−[A]V+M]T and N=[N1,N1]T
    • Where ‘T’ denotes transpose
    • (end of example)


Compared to conventional RTK techniques, the method described above provides twice as many equations (for a rover with two antennas 100) for the same number of unknowns X1 and N1. N1 is referred to as the global ambiguity, which is easier to solve for with more equations. For example, see Whitehead U.S. Pat. No. 6,469,663 for Method and System for GPS and WAAS Carrier Phase Measurements for Relative Positioning, which is incorporated herein by reference.



FIGS. 5A, 5B and 5C show constrained, multiple-antenna systems 101, 102, 103 respectively illustrating operation with different combinations of visible satellites 8. System 101 shows master and slave antennas 16 (Am, As) on a constraining structure 104 (i.e., common rigid body such as a wall, a dam, a marine vessel, a terrestrial vehicle, an aircraft, etc.) whereon the antennas 16 are mounted in fixed relation constrained by distance, geometry and time (by providing a common clock or synchronize clocks). Knowing the orientation of the stricture 103 (e.g. with a compass or other orientation device 32 and known antenna 16 heights) the antennas' locations can be solved in terms of each other. In other words, both antennas 16 locations can be determined in absolute coordinates of an earth-referenced location system with four unique satellites A, B, C and D. By contrast, unconstrained systems would require four satellites to be seen by each antenna for a complete position solution. FIG. 5A shows antenna measurements divided as follows: Am (A, B) and As (C, D). Any other satellite/antenna combination would also provide an attitude/position solution provided four satellites collectively are seen by the antennas 16. With eight unknowns in the position equation for two unconstrained antennas, the constraining structure 104, the time (clock) constraint and the in the orientation device collectively remove four unknowns from the solution equation, enabling positioning operation as shown in FIG. 5A. For example, a satellite-blocking object 105 on the stricture 104 could comprise a crane on a marine vessel, which would otherwise seriously impede GNSS navigation.



FIG. 5B shows a similar system 102 with an unknown clock bias, i.e. independent clocks associated with separate receivers. In this case a total of five GNSS measurements are needed from four unique satellites, one of which (C) is commonly seen by both antennas 16. FIG. 5B shows antenna measurements divided as follows: Am (A, B, C) and As (C, D). Any other satellite/antenna combination would also provide an attitude/position solution provided five satellites collectively are seen by the antennas 16, four unique and one common.



FIG. 5C shows a similar system 103 with common or synchronized clocks (clock constraint) and GNSS-determined attitude using additional measurements at the antennas 16 in place of an orientation device 32. In this case a total of six GNSS measurements are needed from four unique satellites, two of which (C1, C2) are commonly seen by both antennas 16. FIG. 5C shows antenna measurements divided as follows: Am (A, C1, C2) and As (B, C1, C2). Any other satellite/antenna combination would also provide an attitude/position solution provided six satellites collectively are seen by the antennas 16, four unique and two common. In this configuration four measurements provide a location for the master antenna Am and the other two measurements provide a bearing to the slave antenna As. In other words, such constraints remove two unknowns from the position/attitude solution equation enabling a complete solution from six measurements, as compared to the eight measurements required for an unconstrained system. Because of the distance/geometry/clock constraints, the six measurements can be received at the antennas 16 in various combinations.


It will be appreciated that with the systems 101, 102, 103, the common satellites could have been seen in the past and inertial navigation system (INS) devices, such as gyroscopes and/or accelerometers, could be used for guidance during GNSS signal loss. Such INS devices would be previously aligned and calibrated using GNSS-based attitude and effectively provide short-term replacement measurement input. Thus, the common measurements need not be received simultaneously. Upon GNSS signal reacquisition, the INS devices are corrected, aligned and calibrated with GNSS-defined attitude solutions. Moreover, multipath mitigation is provided by both (multiple) antennas 16 receiving measurements from the same satellites. Measurements are generally strengthened by providing more measurements, provided the antennas 16 are separated by more than a wavelength (i.e. greater than a carrier phase wavelength of 0.19 meters) whereby multipath will be uncorrelated and multipath errors will not accumulate linearly. Therefore, the measurement to multipath ratio will improve by approximately the square root of the number of duplicate measurements.



FIGS. 6 and 6A show a GNSS positioning system 110 comprising another alternative aspect of the present invention installed on a vehicle, e.g., an automobile 108, with a first receiver 112 connected to a first antenna 116 (A1) and a second receiver 114 connected to a second antenna 118 (A2). As shown, the antennas 116, 118 are placed inside the windshield (A1) and inside the rear window (A2). The antennas 116, 118 would not interfere with the vehicle's aerodynamics, whereby they would be suitable for racecars and similar vehicles. The antennas are positioned in a predetermined, known relation (i.e. constrained spacing and orientation) for determining attitude comprising yaw and pitch, although the primary concern would normally be yaw for terrestrial vehicle navigation purposes. The use of multiple antennas in a predetermined spacing is particularly advantageous because one or more of the antennas may be partially blocked by the vehicle or other objects, such as buildings and other structures. Moreover, sloping antenna mounts on windshields and rear windows can reduce satellite visibility, which can be compensated for by providing multiple antennas.



FIG. 6A shows the dimensional relationship between the antennas A1 and A2. D1 is the longitudinal (Y axis) spacing between A1 and A2, Such predetermined spacings and relative orientations (i.e., constraints) can be stored in the CPU 113 for use in connection with GNSS-position calculations. FIG. 6A also shows a configuration where antennas A1 and A2 are partially blocked from the satellite 8 constellation whereby the CPU calculates a position solution using partial GNSS ranging data and predetermined antenna positions and orientations (attitudes).



FIG. 7 shows a GNSS positioning system 130 comprising another alternative aspect of the invention with a master antenna 136 (Am) located on top of the fuselage and slave antennas 138 (As1, As2) located on opposite sides of the aircraft fuselage 132. The top-mounted master antenna 136 (Am) will normally see all available satellites, although the slave antennas 138 (As1, As2) may be partially blocked, for example, with the aircraft is in a banking attitude. The antennas Am, As1, As2 are positioned in a predetermined, known relation (i.e. constrained distance and geometry) for determining attitude comprising yaw, pitch and roll. With the aircraft in a banking (rolling) attitude as shown, the antennas see different combinations of satellites 8 because the fuselage is positioned to block the view of certain satellites 8, particularly on the downward side of the fuselage 132. However, collectively the antennas can provide a complete positioning solution, including the attitude/orientation of the aircraft, using the procedures discussed above. The antennas 136, 138 can be connected to receiver units similar to those described above. The positioning system 110 can be configured similarly to the positioning systems described above with a CPU 113, a clock 117 and an antenna switch control comprising a receiver unit 111 for switching between the slave antennas 118 (As1, As2). The receiver unit 111 can also include an optional orientation device 32.



FIG. 8 shows a GNSS system 152 comprising yet another alternative aspect or embodiment of the present invention and including a rover receiver unit 154 connected to an antenna array 156, which is mounted on top of a vehicle 158, such as a tractor with a cab 160. The receiver unit 154 includes a rover receiver 162 connected to a clock 164 and a microprocessor 166. An antenna switch control 168 is connected to the receiver 162, the microprocessor 166 and the antenna array 156 and functions to switch among individual antennas 170 comprising the antenna array 156.


The spacing between each pair of antennas 170 is preferably significant, i.e. greater than a carrier phase wavelength of 0.19 meters whereby multipath is uncorrelated, but the mean position is correlated. Therefore, the multipath error factor will increase by the square root of the number of antennas but the position signals and noise ratio will improve by the number of antennas. The position-to-multipath ratio thus improves by the square root of the number of antennas over configurations with fewer antennas. For example, increasing from two antennas to four theoretically results in a 2:2 improvement in accuracy from multipath.


An advantage of this configuration is that by using more antennas, less expensive individual antennas can be used with smaller spacing, and provide equal or better performance compared to fewer, more expensive antennas. For example, eight individual antennas 170 comprise the antenna array 156 in FIG. 8. In a relatively slow-moving vehicle, such as an agricultural tractor 158, the rover receiver unit 154 could switch among the antennas 170 via the switch control 168. Because multipath is a slow-changing error, the effective multi-antenna improvements in accuracy from multipath described above can be achieved. Moreover, a single rover receiver 162 can be utilized in place of multiple receivers, thereby eliminating any group temperature variation delays, which are generally more significant in receivers than in antennas.



FIG. 9 shows another application of the invention comprising a GNSS guidance system 172 with a two-antenna array 174 comprising a pair of antennas 176 connected by a rigid link 178 and driven by a stepper motor 180. A rover receiver unit 182 includes a stepper motor control 184. In operation, the system 172 functions similarly to the GNSS guidance system 152 described above and shown in FIG. 8, with the incremental positioning of the antenna array 174 replacing the antenna ring array 156 whereby readings can be obtained at predetermined angular increments with predetermined orientations with respect to a vehicle or other object of interest. The rotational speed of the stepper motor 180 can be adjustable and matched to the travel velocity of a vehicle whereby ranging signals can be received at suitable intervals as necessary to guide with the desired accuracy. A suitable radome 185 is placed over the antenna array 174.



FIG. 10 shows yet another application of the invention comprising a GNSS guidance system 186 with an antenna array 188, which is similar to the antenna array 174 described above. A motor 190 continuously rotates the antenna array 188 and is connected to a rotary encoder 192 producing output signals corresponding to the angular orientation of the array 188, and a gyroscope 194 producing output signals corresponding to its rate of turn. The operation of the motor 190, which can be adjusted corresponding to the velocity of the vehicle, is controlled by an antenna motor control 191 in a rover receiver unit 193. The output signals are processed by the processor 166, which can take into account and correct for yaw, pitch and roll for the system 186 and an associated vehicle. Ranging signals received by the antenna array 188 are correlated with various angular orientations and processed as described above. A suitable radome 185 is placed over the antenna array 174.


It will be appreciated that the satellite systems as discussed herein may include but not be limited to Wide Area Augmentation System (WAAS), Global Navigation Satellite System (GNSS) including GPS, GLONASS and other satellite ranging technologies. The term WAAS is used herein as a generic reference to all GNSS augmentation systems which, to date, include three programs: WAAS (Wide Area Augmentation System) in the USA, EGNOS (European Geostationary Navigation Overlay System) in Europe and MSAS (Multifunctional Transport Satellite Space-based Augmentation System) in Japan. Each of these three systems, which are all compatible, consists of a ground network for observing the GPS constellation, and one or more geostationary satellites.


It will also be appreciated that while a particular series of steps or procedures is described as part of the abovementioned process, no order of steps should necessarily be inferred from the order of presentation. For example, the process includes receiving one or more sets of satellite signals. It should be evident that the order of receiving the satellite signals is variable and could be reversed without impacting the methodology disclosed herein or the scope of the claims.


It should further be appreciated that while an exemplary partitioning functionality has been provided, it should be apparent to one skilled in the art, that the partitioning could be different. For example, the control of the master receiver 12 and slave receiver 14, could be integrated in any, or another unit. The processes may for ease of implementation, be integrated into a single unit. Such configuration variances should be considered equivalent and within the scope of the disclosure and Claims herein.


The disclosed invention may be embodied in the form of computer-implemented processes and apparatuses for practicing those processes. The present invention can also be embodied in the form of computer program code containing instructions embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other computer-readable storage medium, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus for practicing the invention. The present invention can also be embodied in the form of computer program code, for example, whether stored in a storage medium, loaded into and/or executed by a computer, or as data signal transmitted whether a modulated carrier wave or not, over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus for practicing the invention. When implemented on a general-purpose microprocessor, the computer program code segments configure the microprocessor to create specific logic circuits.


While the description has been made with reference to exemplary embodiments, it will be understood by those of ordinary skill in the pertinent art that various changes may be made and equivalents may be substituted for the elements thereof without departing from the scope of the disclosure. In addition, numerous modifications may be made to adapt the teachings of the disclosure to a particular object or situation without departing from the essential scope thereof. Therefore, it is intended that the Claims not be limited to the particular embodiments disclosed as the currently preferred best modes contemplated for carrying out the teachings herein, but that the Claims shall cover all embodiments falling within the true scope and spirit of the disclosure.

Claims
  • 1. A GNSS-based method of generating a navigation solution of a rover vehicle comprising the steps of: providing a base receiver unit at a fixed location with a base GNSS receiver having a first clock, a base RF transmitter, and a base processor;providing a base GNSS antenna connected to said base receiver;providing a rover receiver unit with a rover GNSS receiver having a second clock, a rover RF receiver, and a rover processor;providing a rover antenna array comprising multiple rover GNSS antennas, each connected to said rover receiver;providing a stepper motor mounted on said vehicle and connected to said rover antenna array;providing an encoder connected to said rover antenna array;providing a gyroscope connected to said rover antenna array;said rover processor receiving outputs from said encoder and said gyroscope;providing said rover antenna array with multiple incremental orientations each corresponding to a portion of a full revolution;programming said rover processor to rotate said stepper motor through said multiple incremental orientations relative to said vehicle;adjusting the rotational speed of the stepper motor based on the rover vehicle velocity to receive ranging signals at intervals timed to the multiple, incremental orientations of the antenna array;synchronizing the first and second clocks;receiving GNSS signals at each of the base GNSS antenna and the rover GNSS antennas and respectively measuring carrier phase observations therefrom;transmitting base receiver unit carrier phase observations using the base RF transmitter;receiving said base receiver unit carrier phase observations using the rover RF receiver;determining an attitude solution of said rover antenna array by solving for locations of the rover antennas with respect to each other on the basis of the base receiver unit carrier phase observations and the rover receiver unit carrier phase observations;using said rover antenna array attitude solution to resolve global, carrier phase GNSS integer ambiguities at each antenna;storing with said rover processor the attitudes and relative positions of said rover antenna array relative to said vehicle at said incremental orientations;time-tagging and storing the relative antenna locations, antenna array attitude solution and carrier phase observations; andsaid rover processor determining a position solution for said vehicle based on GNSS signals received by said rover antenna array, the time-tagged relative antenna locations, the antenna array attitude solution, received GNSS carrier phase observations, and said encoder and gyroscope outputs.
  • 2. The method of claim 1, which includes the additional steps of: storing in the rover processor the determined position solutions for said vehicle at a plurality of spaced points in time.
  • 3. The method of claim 2, which includes additional steps of: comparing the position and attitude solutions at first and second spaced points in time to provide an output representing movement of the vehicle based on rover vehicle antenna and antenna array position and attitude differences at said spaced points in time.
CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation-in-part of and claims the benefit of: U.S. patent application Ser. No. 10/804,758, filed Mar. 19, 2004, now U.S. Pat. No. 7,400,956 and Ser. No. 10/828,745, filed Apr. 21, 2004, now abandoned; and claims the benefit of U.S. Provisional patent application Ser. No. 60/456,146, filed Mar. 20, 2003 and No. 60/464,756, filed Apr. 23, 2003. This application is also related to U.S. patent application Ser. No. 12/554,741, filed Sep. 4, 2009, and Ser. No. 12/683,994, filed Jan. 7, 2010, now U.S. Pat. No. 8,138,970. The contents of all of the aforementioned applications are incorporated by reference herein in their entireties.

US Referenced Citations (420)
Number Name Date Kind
3585537 Rennick et al. Jun 1971 A
3596228 Reed, Jr. et al. Jul 1971 A
3727710 Sanders et al. Apr 1973 A
3815272 Marleau Jun 1974 A
3899028 Morris et al. Aug 1975 A
3987456 Gelin Oct 1976 A
4132272 Holloway et al. Jan 1979 A
4170776 MacDoran et al. Oct 1979 A
4180133 Collogan et al. Dec 1979 A
4398162 Nagai Aug 1983 A
4453614 Allen et al. Jun 1984 A
4529990 Brunner Jul 1985 A
4599620 Evans Jul 1986 A
4637474 Leonard Jan 1987 A
4667203 Counselman, III May 1987 A
4689556 Cedrone Aug 1987 A
4694264 Owens et al. Sep 1987 A
4710775 Coe Dec 1987 A
4714435 Stipanuk et al. Dec 1987 A
4739448 Rowe et al. Apr 1988 A
4751512 Longaker Jun 1988 A
4769700 Pryor Sep 1988 A
4785463 Janc et al. Nov 1988 A
4802545 Nystuen et al. Feb 1989 A
4812991 Hatch Mar 1989 A
4814036 Hatch Mar 1989 A
4858132 Holmquist Aug 1989 A
4864320 Munson et al. Sep 1989 A
4894662 Counselman Jan 1990 A
4916577 Dawkins Apr 1990 A
4918607 Wible Apr 1990 A
4963889 Hatch Oct 1990 A
5031704 Fleischer et al. Jul 1991 A
5100229 Lundberg et al. Mar 1992 A
5134407 Lorenz et al. Jul 1992 A
5148179 Allison Sep 1992 A
5152347 Miller Oct 1992 A
5155490 Spradley et al. Oct 1992 A
5155493 Thursby et al. Oct 1992 A
5156219 Schmidt et al. Oct 1992 A
5165109 Han et al. Nov 1992 A
5173715 Rodal et al. Dec 1992 A
5177489 Hatch Jan 1993 A
5185610 Ward et al. Feb 1993 A
5191351 Hofer et al. Mar 1993 A
5202829 Geier Apr 1993 A
5207239 Schwitalia May 1993 A
5239669 Mason et al. Aug 1993 A
5255756 Follmer et al. Oct 1993 A
5268695 Dentinger et al. Dec 1993 A
5293170 Lorenz et al. Mar 1994 A
5294970 Dornbusch et al. Mar 1994 A
5296861 Knight Mar 1994 A
5311149 Wagner et al. May 1994 A
5323322 Mueller et al. Jun 1994 A
5334987 Teach Aug 1994 A
5343209 Sennott et al. Aug 1994 A
5345245 Ishikawa et al. Sep 1994 A
5347286 Babitch Sep 1994 A
5359332 Allison et al. Oct 1994 A
5361212 Class et al. Nov 1994 A
5365447 Dennis Nov 1994 A
5369589 Steiner Nov 1994 A
5375059 Kyrtsos et al. Dec 1994 A
5390124 Kyrtsos Feb 1995 A
5390125 Sennott et al. Feb 1995 A
5390207 Fenton et al. Feb 1995 A
5416712 Geier et al. May 1995 A
5442363 Remondi Aug 1995 A
5444453 Lalezari Aug 1995 A
5451964 Babu Sep 1995 A
5467282 Dennis Nov 1995 A
5471217 Hatch et al. Nov 1995 A
5476147 Fixemer Dec 1995 A
5477228 Tiwari et al. Dec 1995 A
5477458 Loomis Dec 1995 A
5490073 Kyrtsos Feb 1996 A
5491636 Robertson Feb 1996 A
5495257 Loomis Feb 1996 A
5504482 Schreder Apr 1996 A
5511623 Frasier Apr 1996 A
5519620 Talbot et al. May 1996 A
5521610 Rodal May 1996 A
5523761 Gildea Jun 1996 A
5534875 Diefes et al. Jul 1996 A
5543804 Buchler et al. Aug 1996 A
5546093 Gudat et al. Aug 1996 A
5548293 Cohen Aug 1996 A
5561432 Knight Oct 1996 A
5563786 Torii Oct 1996 A
5568152 Janky et al. Oct 1996 A
5568162 Samsel et al. Oct 1996 A
5583513 Cohen Dec 1996 A
5589835 Gildea et al. Dec 1996 A
5592382 Colley Jan 1997 A
5596328 Stangeland et al. Jan 1997 A
5600670 Turney Feb 1997 A
5604506 Rodal Feb 1997 A
5608393 Hartman Mar 1997 A
5610522 Locatelli et al. Mar 1997 A
5610616 Vallot et al. Mar 1997 A
5610845 Slabinski Mar 1997 A
5612883 Shaffer et al. Mar 1997 A
5615116 Gudat et al. Mar 1997 A
5617100 Akiyoshi et al. Apr 1997 A
5617317 Ignagni Apr 1997 A
5621646 Enge et al. Apr 1997 A
5638077 Martin Jun 1997 A
5644139 Allen Jul 1997 A
5664632 Frasier Sep 1997 A
5673491 Brenna et al. Oct 1997 A
5680140 Loomis Oct 1997 A
5684696 Rao et al. Nov 1997 A
5706015 Chen et al. Jan 1998 A
5717593 Gvili Feb 1998 A
5725230 Walkup Mar 1998 A
5731786 Abraham et al. Mar 1998 A
5739785 Allison Apr 1998 A
5757316 Buchler May 1998 A
5765123 Nimura et al. Jun 1998 A
5777578 Chang et al. Jul 1998 A
5810095 Orbach et al. Sep 1998 A
5828336 Yunck et al. Oct 1998 A
5838562 Gudat et al. Nov 1998 A
5854987 Sekine et al. Dec 1998 A
5862501 Talbot et al. Jan 1999 A
5864315 Welles et al. Jan 1999 A
5864318 Cozenza et al. Jan 1999 A
5875408 Bendett et al. Feb 1999 A
5877725 Kalafus Mar 1999 A
5890091 Talbot et al. Mar 1999 A
5899957 Loomis May 1999 A
5906645 Kagawa et al. May 1999 A
5912798 Chu Jun 1999 A
5914685 Kozlov et al. Jun 1999 A
5917448 Mickelson Jun 1999 A
5918558 Susag Jul 1999 A
5919242 Greatline et al. Jul 1999 A
5923270 Sampo et al. Jul 1999 A
5926079 Heine et al. Jul 1999 A
5927603 McNabb Jul 1999 A
5928309 Korver et al. Jul 1999 A
5929721 Munn et al. Jul 1999 A
5933110 Tang Aug 1999 A
5935183 Sahm et al. Aug 1999 A
5936573 Smith Aug 1999 A
5940026 Popeck Aug 1999 A
5941317 Mansur Aug 1999 A
5943008 Van Dusseldorp Aug 1999 A
5944770 Enge et al. Aug 1999 A
5945917 Harry Aug 1999 A
5949371 Nichols Sep 1999 A
5955973 Anderson Sep 1999 A
5956250 Gudat et al. Sep 1999 A
5969670 Kalafus et al. Oct 1999 A
5987383 Keller Nov 1999 A
6014101 Loomis Jan 2000 A
6014608 Seo Jan 2000 A
6018313 Engelmayer et al. Jan 2000 A
6023239 Kovach Feb 2000 A
6023240 Sutton Feb 2000 A
6052647 Parkinson et al. Apr 2000 A
6055477 McBurney et al. Apr 2000 A
6057800 Yang et al. May 2000 A
6061390 Meehan et al. May 2000 A
6061632 Dreier May 2000 A
6062317 Gharsalli May 2000 A
6069583 Silvestrin et al. May 2000 A
6076612 Carr et al. Jun 2000 A
6081171 Ella Jun 2000 A
6100842 Dreier et al. Aug 2000 A
6104978 Harrison et al. Aug 2000 A
6122595 Varley et al. Sep 2000 A
6128574 Diekhans Oct 2000 A
6144335 Rogers Nov 2000 A
6191730 Nelson, Jr. Feb 2001 B1
6191733 Dizchavez Feb 2001 B1
6198430 Hwang et al. Mar 2001 B1
6198992 Winslow Mar 2001 B1
6199000 Keller et al. Mar 2001 B1
6205401 Pickhard et al. Mar 2001 B1
6211821 Ford Apr 2001 B1
6215828 Signell et al. Apr 2001 B1
6229479 Kozlov et al. May 2001 B1
6230097 Dance et al. May 2001 B1
6233511 Berger et al. May 2001 B1
6236916 Staub et al. May 2001 B1
6236924 Motz May 2001 B1
6253160 Hanseder Jun 2001 B1
6256583 Sutton Jul 2001 B1
6259398 Riley Jul 2001 B1
6266595 Greatline et al. Jul 2001 B1
6285320 Olster et al. Sep 2001 B1
6292132 Wilson Sep 2001 B1
6307505 Green Oct 2001 B1
6313788 Wilson Nov 2001 B1
6314348 Winslow Nov 2001 B1
6325684 Knight Dec 2001 B1
6336066 Pellenc et al. Jan 2002 B1
6345231 Quincke Feb 2002 B2
6356602 Rodal et al. Mar 2002 B1
6377889 Soest Apr 2002 B1
6380888 Kucik Apr 2002 B1
6389345 Phelps May 2002 B2
6392589 Rogers et al. May 2002 B1
6397147 Whitehead May 2002 B1
6415229 Diekhans Jul 2002 B1
6418031 Archambeault Jul 2002 B1
6421003 Riley et al. Jul 2002 B1
6424915 Fukuda et al. Jul 2002 B1
6431576 Viaud et al. Aug 2002 B1
6434462 Bevly et al. Aug 2002 B1
6441779 Bennett et al. Aug 2002 B1
6445983 Dickson et al. Sep 2002 B1
6445990 Manring Sep 2002 B1
6449558 Small Sep 2002 B1
6463091 Zhodzicshsky et al. Oct 2002 B1
6463374 Keller et al. Oct 2002 B1
6466871 Reisman et al. Oct 2002 B1
6469663 Whitehead et al. Oct 2002 B1
6484097 Fuchs et al. Nov 2002 B2
6501422 Nichols Dec 2002 B1
6515619 McKay, Jr. Feb 2003 B1
6516271 Upadhyaya et al. Feb 2003 B2
6539303 McClure et al. Mar 2003 B2
6542077 Joao Apr 2003 B2
6549835 Deguchi Apr 2003 B2
6553299 Keller et al. Apr 2003 B1
6553300 Ma et al. Apr 2003 B2
6553311 Aheam et al. Apr 2003 B2
6570534 Cohen et al. May 2003 B2
6577952 Geier et al. Jun 2003 B2
6587761 Kumar Jul 2003 B2
6606542 Hauwiller et al. Aug 2003 B2
6611228 Toda et al. Aug 2003 B2
6611754 Klein Aug 2003 B2
6611755 Coffee et al. Aug 2003 B1
6622091 Perlmutter et al. Sep 2003 B2
6631394 Ronkka et al. Oct 2003 B1
6631916 Miller Oct 2003 B1
6643576 O'Connor et al. Nov 2003 B1
6646603 Dooley et al. Nov 2003 B2
6657875 Zeng et al. Dec 2003 B1
6671587 Hrovat et al. Dec 2003 B2
6686878 Lange Feb 2004 B1
6688403 Bernhardt et al. Feb 2004 B2
6703973 Nichols Mar 2004 B1
6711501 McClure et al. Mar 2004 B2
6721638 Zeitler Apr 2004 B2
6732024 Rekow et al. May 2004 B2
6744404 Whitehead et al. Jun 2004 B1
6754584 Pinto et al. Jun 2004 B2
6774843 Takahashi Aug 2004 B2
6792380 Toda Sep 2004 B2
6819269 Flick Nov 2004 B2
6822314 Beasom Nov 2004 B2
6865465 McClure Mar 2005 B2
6865484 Miyasaka et al. Mar 2005 B2
6879283 Bird et al. Apr 2005 B1
6882312 Vorobiev et al. Apr 2005 B1
6900992 Kelly et al. May 2005 B2
6922635 Rorabaugh Jul 2005 B2
6931233 Tso et al. Aug 2005 B1
6961018 Heppe et al. Nov 2005 B2
6967538 Woo Nov 2005 B2
6990399 Hrazdera et al. Jan 2006 B2
7006032 King et al. Feb 2006 B2
7026982 Toda et al. Apr 2006 B2
7027918 Zimmerman et al. Apr 2006 B2
7031725 Rorabaugh Apr 2006 B2
7089099 Shostak et al. Aug 2006 B2
7142956 Heiniger et al. Nov 2006 B2
7155335 Rennels Dec 2006 B2
7162348 McClure et al. Jan 2007 B2
7191061 McKay et al. Mar 2007 B2
7221314 Brabec et al. May 2007 B2
7231290 Steichen et al. Jun 2007 B2
7248211 Hatch et al. Jul 2007 B2
7268727 Montgomery et al. Sep 2007 B2
7271766 Zimmerman et al. Sep 2007 B2
7277784 Weiss Oct 2007 B2
7292186 Miller et al. Nov 2007 B2
7324915 Altmann Jan 2008 B2
7358896 Gradincic et al. Apr 2008 B2
7373231 McClure et al. May 2008 B2
7388539 Whitehead et al. Jun 2008 B2
7395769 Jensen Jul 2008 B2
7428259 Wang et al. Sep 2008 B2
7437230 McClure et al. Oct 2008 B2
7451030 Eglington et al. Nov 2008 B2
7479900 Horstemeyer Jan 2009 B2
7489270 Lawrence et al. Feb 2009 B2
7505848 Flann et al. Mar 2009 B2
7522099 Zhodzishsky et al. Apr 2009 B2
7522100 Yang et al. Apr 2009 B2
7571029 Dai et al. Aug 2009 B2
7689354 Heiniger et al. Mar 2010 B2
20030014171 Ma et al. Jan 2003 A1
20030093210 Kondo et al. May 2003 A1
20030187560 Keller et al. Oct 2003 A1
20030208319 Ell et al. Nov 2003 A1
20040039514 Steichen et al. Feb 2004 A1
20040212533 Whitehead Oct 2004 A1
20050080559 Ishibashi et al. Apr 2005 A1
20050225955 Grebenkemper et al. Oct 2005 A1
20050242991 Montgomery et al. Nov 2005 A1
20050265494 Goodlings Dec 2005 A1
20060031664 Wilson et al. Feb 2006 A1
20060167600 Nelson et al. Jul 2006 A1
20060206246 Walker Sep 2006 A1
20060215739 Williamson et al. Sep 2006 A1
20060244656 Lawrence et al. Nov 2006 A1
20070078570 Dai et al. Apr 2007 A1
20070088447 Stothert et al. Apr 2007 A1
20070121708 Simpson May 2007 A1
20070205940 Yang et al. Sep 2007 A1
20070285308 Bauregger et al. Dec 2007 A1
20080129586 Martin Jun 2008 A1
20080204312 Euler Aug 2008 A1
20090171583 DiEsposti Jul 2009 A1
20090174597 DiLellio et al. Jul 2009 A1
20090174622 Kanou Jul 2009 A1
20090177395 Stelpstra Jul 2009 A1
20090177399 Park et al. Jul 2009 A1
20090259397 Stanton Oct 2009 A1
20090259707 Martin et al. Oct 2009 A1
20090262014 DiEsposti Oct 2009 A1
20090262018 Vasilyev et al. Oct 2009 A1
20090262974 Lithopoulos Oct 2009 A1
20090265054 Basnayake Oct 2009 A1
20090265101 Jow Oct 2009 A1
20090265104 Shroff Oct 2009 A1
20090273372 Brenner Nov 2009 A1
20090273513 Huang Nov 2009 A1
20090274079 Bhatia et al. Nov 2009 A1
20090274113 Katz Nov 2009 A1
20090276155 Jeerage et al. Nov 2009 A1
20090295633 Pinto et al. Dec 2009 A1
20090295634 Yu et al. Dec 2009 A1
20090299550 Baker Dec 2009 A1
20090322597 Medina Herrero et al. Dec 2009 A1
20090322598 Fly et al. Dec 2009 A1
20090322600 Whitehead et al. Dec 2009 A1
20090322601 Ladd et al. Dec 2009 A1
20090322606 Gronemeyer Dec 2009 A1
20090326809 Colley et al. Dec 2009 A1
20100013703 Tekawy et al. Jan 2010 A1
20100026569 Amidi Feb 2010 A1
20100030470 Wang et al. Feb 2010 A1
20100039316 Gronemeyer et al. Feb 2010 A1
20100039318 Kmiecik Feb 2010 A1
20100039320 Boyer et al. Feb 2010 A1
20100039321 Abraham Feb 2010 A1
20100060518 Bar-Sever et al. Mar 2010 A1
20100063649 Wu et al. Mar 2010 A1
20100084147 Aral Apr 2010 A1
20100085249 Ferguson et al. Apr 2010 A1
20100085253 Ferguson et al. Apr 2010 A1
20100103033 Roh Apr 2010 A1
20100103034 Tobe et al. Apr 2010 A1
20100103038 Yeh et al. Apr 2010 A1
20100103040 Broadbent Apr 2010 A1
20100106414 Whitehead Apr 2010 A1
20100106445 Kondoh Apr 2010 A1
20100109944 Whitehead et al. May 2010 A1
20100109945 Roh May 2010 A1
20100109947 Rintanen May 2010 A1
20100109948 Razoumov et al. May 2010 A1
20100109950 Roh May 2010 A1
20100111372 Zheng et al. May 2010 A1
20100114483 Heo et al. May 2010 A1
20100117894 Velde et al. May 2010 A1
20100117899 Papadimitratos et al. May 2010 A1
20100117900 van Diggelen et al. May 2010 A1
20100121577 Zhang et al. May 2010 A1
20100124210 Lo May 2010 A1
20100124212 Lo May 2010 A1
20100134354 Lennen Jun 2010 A1
20100149025 Meyers et al. Jun 2010 A1
20100149030 Verma et al. Jun 2010 A1
20100149033 Abraham Jun 2010 A1
20100149034 Chen Jun 2010 A1
20100149037 Cho Jun 2010 A1
20100150284 Fielder et al. Jun 2010 A1
20100152949 Nunan et al. Jun 2010 A1
20100156709 Zhang et al. Jun 2010 A1
20100156712 Pisz et al. Jun 2010 A1
20100156718 Chen Jun 2010 A1
20100159943 Salmon Jun 2010 A1
20100161179 McClure et al. Jun 2010 A1
20100161211 Chang Jun 2010 A1
20100161568 Xiao Jun 2010 A1
20100171660 Shyr et al. Jul 2010 A1
20100171757 Melamed Jul 2010 A1
20100185364 McClure Jul 2010 A1
20100185366 Heiniger et al. Jul 2010 A1
20100185389 Woodard Jul 2010 A1
20100188285 Collins Jul 2010 A1
20100188286 Bickerstaff et al. Jul 2010 A1
20100189163 Burgi et al. Jul 2010 A1
20100207811 Lackey Aug 2010 A1
20100210206 Young Aug 2010 A1
20100211248 Craig et al. Aug 2010 A1
20100211315 Toda Aug 2010 A1
20100211316 DaSilva Aug 2010 A1
20100220004 Malkos et al. Sep 2010 A1
20100220008 Conover et al. Sep 2010 A1
20100222076 Poon et al. Sep 2010 A1
20100225537 Abraham Sep 2010 A1
20100228408 Ford Sep 2010 A1
20100228480 Lithgow et al. Sep 2010 A1
20100231443 Whitehead Sep 2010 A1
20100231446 Marshall et al. Sep 2010 A1
20100232351 Chansarkar et al. Sep 2010 A1
20100235093 Chang Sep 2010 A1
20100238976 Young Sep 2010 A1
20100241347 King et al. Sep 2010 A1
20100241353 Park Sep 2010 A1
20100241441 Page et al. Sep 2010 A1
20100241864 Kelley et al. Sep 2010 A1
Foreign Referenced Citations (12)
Number Date Country
07244150 Sep 1995 JP
WO9836288 Aug 1998 WO
WO0024239 May 2000 WO
WO03019430 Mar 2003 WO
WO2005119386 Dec 2005 WO
WO2009066183 May 2009 WO
WO-2009082745 Jul 2009 WO
WO-2010005945 Jul 2009 WO
WO2009126587 Oct 2009 WO
WO2009148638 Dec 2009 WO
WO-2010104782 Sep 2010 WO
WO-2011014431 Feb 2011 WO
Non-Patent Literature Citations (39)
Entry
Lu, Gang. “Development of a GPS Multi-Antenna System for Attitude Determination,” Dept. of Geomatics Engineering, pp. 1-196, Jan. 1995.
Parkinson, Bradford W., et al., “Global Positioning System: Theory and Applications, vol. II”, Bradford W. Parkinson and James J. Spiker, Jr., eds., Global Postioning System: Theory and Applicaitons, vol. II, 1995, AIAA, Reston, VA, USA, pp. 3-50, (1995),3-50.
“Orthman Manufacturing Co., www.orthman.com/htm;guidance.htm”, 2004, regarding the “Tracer Quick-Hitch”.
Lin, Dai et al., “Real-time Attitude Determination fro Microsatellite by Lamda Method Combined with Kalman Filtering”, A Collection fof the 22nd AIAA International Communications Satellite Systems Conference and Exhibit Technical Paers vol. 1, Monetrey, California American Institute of Aeronautics and Astronautics, Inc., (May 2004),136-143.
Xu, Jiangning et al., “An EHW Architecture for Real-Time GPS Attitude Determination Based on Parallel Genetic Algorithm”, The Computer SocietyProceedings of the 2002 NASA/DOD Conference on Evolvable Hardware (EH'02), (2002).
Han, Shaowel et al., “Single-Epoch Ambiguity Resolution for Real-Time GPS Attitude Determination with the Aid of One-Dimensional Optical Fiber Gyro”, GPS Solutions, vol. 3, No. 1, pp. 5-12 (1999) John Wiley & Sons, Inc.
Park, Chansik et al., “Integer Ambiguity Resolution for GPS Based Attitude Determination System”, SICE , Jul. 29-31, 1998, Chiba, 1115-1120.
Last, J. D., et al., “Effect of skywave interference on coverage of radiobeacon DGPS stations”, IEEE Proc.-Radar, Sonar Navig.. vol. 144, No. 3, Jun. 1997, pp. 163-168.
“International Search Report and Written Opinion”, PCT/US2004/015678, filed May 17, 2004, Jun. 21, 2005.
“ISO”, 11783 Part 7 Draft Amendment 1 Annex Paragraphs B.6 and B.7.ISO 11783-7 2004 DAM1, ISO: Mar. 8, 2004.
Kaplan, E D., “Understanding GPS: Principles and Applications”, Artech House, MA, 1998.
Irsigler, M et al., “PPL Tracking Performance in the Presence of Oscillator Phase Noise”, GPS Solutions, vol. 5, No. 4, pp. 45-57 (2002).
Ward, Phillip W., “Performance Comparisons Between FLL, PLL and a Novel FLL-Assisted-PLL Carrier Tracking Loop Under RF Interference Conditions”, 11th Int. Tech Meeting of the Satellite Division of the U.S. Inst. of Navigation, Nashville, TN, Sep. 15-18, 1998, 783-795.
Bevly, David M., “Comparison of INS v. Carrier-Phase DGPS for Attitude Determination in the Control of Off-Road Vehicles”, ION 55th Annual Meeting; Jun. 28-30, 1999; Cambridge, Massachusetts; pp. 497-504.
Keicher, R. et al., “Automatic Guidance for Agricultural Vehicles in Europe”, Computers and Electronics in Agriculture, vol. 25, (Jan. 2000), 169-194.
“International Search Report”, PCT/US09/49776, (Aug. 11, 2009).
“International Search Report”, PCT/AU/2008/000002, (Feb. 28, 2008).
“International Search Report”, PCT/US09/33693, (Mar. 30, 2009).
“International Search Report”, PCT/US09/039686, (May 26, 2009).
“International Search Report”, PCT/US09/34376, (Nov. 2, 2009).
“ISR Notification & Written Opinion”, PCT/US10/26509, (Apr. 20, 2010),1-7.
“Notification Concerning Transmittal of International Report on Patentability (PCT)”, PCT/US2009/049776, (Jan. 20, 2011).
“Notification of Transmittal of InternatinalPrelim. Report of Patentability”, International Application No. PCT/US09/039686, (Oct. 21, 2010).
“International Search Report and Written Opinion”, PCT/US2010/043094, (Sep. 17, 2010).
“Notification of Publication of International Application”, WO 2011/014431, (Feb. 3, 2011).
“International Search Report and Written Opinion”, PCT/US08/81727, (Dec. 23, 2008).
Takac, Frank et al., “SmartRTK: A Novel Method of Processing Standardised RTCM Network RTK Information for High Precision Positioning”, Proceedings of ENC GNSS 2008, Toulouse, France,(Apr. 22, 2008).
“International Search Report”, PCT/US09/33567, (Feb. 9, 2009).
“International Search Report and Written Opinion”, PCT/IB2008/003796, (Jul. 15, 2009).
“International Search Report / Written Opinion”, PCT/US09/63594, (Jan. 11, 2010).
“International Search Report”, PCT/US09/60668, (Dec. 9, 2009).
“International Search Report”, PCT/US09/067693, (Jan. 26, 2010).
“International Search Report and Written Opinion”, PCT/US10/21334, (Mar. 12, 2010).
Rho, Hyundho et al., “Dual-Frequency GPS Precise Point Positioning with WADGPS Corrections”, [retrieved on May 18, 2010]. Retrieved from the Internet: ,URL:http://gauss.gge.unb.ca/papers.pdf/iongnss2005.rho.wadgps.pdf, (Jul. 12, 2006).
“Eurocontrol, Pegasus Technical Notes on SBAS”, report [online], Dec. 7, 2004 [retrieved on May 18, 2010]. Retrieved from the Internet: <URL: http://www.icao.int/icao/en/ro/nacc/meetings/2004/gnss/documentation/Pegasus/tn.pdf>, (Dec. 7, 2004), p. 89 paras [0001]-[0004].
“ARINC Engineering Services, Interface Specification IS-GPS-200, Revision D”, Online [retrieved on May 18, 2010]. Retrieved from the Internet;<URL: http://www.navcen.uscg.gov/gps/geninfo/IS-GPS-200D.pdf>, (Dec. 7, 2004), p. 168 para [0001].
Schaer, et al., “Determination and Use of GPS Differential Code Bias Values”, Presentation [online]. Revtrieved May 18, 2010. Retrieved from the Internet: <http://nng.esoc.esa.de/ws2006/REPR2.pdf>, (May 8, 2006).
“International Search Report”, PCT/US10/26509., (Apr. 20, 2010).
“International Search Report and Written Opinion”, International Searching Authortiy, PCT/US08/88070, Feb. 9, 009.
Related Publications (1)
Number Date Country
20090121932 A1 May 2009 US
Provisional Applications (2)
Number Date Country
60456146 Mar 2003 US
60464756 Apr 2003 US
Continuation in Parts (2)
Number Date Country
Parent 10804758 Mar 2004 US
Child 12350431 US
Parent 10828745 Apr 2004 US
Child 10804758 US