Bearing calculation

Information

  • Patent Grant
  • 10667145
  • Patent Number
    10,667,145
  • Date Filed
    Monday, June 20, 2016
    8 years ago
  • Date Issued
    Tuesday, May 26, 2020
    4 years ago
Abstract
There is provided an apparatus comprising a location receiver to receive a signal indicative of a location of the location receiver. A movement mechanism rotates the location receiver about an axis and provides an angle of the location receiver about the axis from a known point. Calculation circuitry produces an output value indicative of an absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.
Description
CROSS-REFERENCE TO OTHER APPLICATIONS

This application claims priority under 35 U.S.C. 119(b) to UK Application No. 1511200.6, filed 25 Jun. 2015 and entitled “STEERABLE ANTENNA SYSTEM”, and to UK Application No. 1514938.8, filed 21 Aug. 2015 and entitled “BEARING CALCULATION”. Both applications are incorporated by reference herein.


BACKGROUND

The present technique relates to the field of data processing and more particularly the calculation of bearings.


It is known to use global navigation system signals in order to determine the location of a device capable of receiving such signals. For example, such navigation systems signals include GPS (The Global Positioning System), GLONASS (The Russian Global Navigation Satellite System), the planned European Union Galileo positioning system, India's Indian Regional Navigation Satellite System and the Chinese BeiDou Navigation Satellite System as well as ground-based navigation system signals such as is used in A-GPS (assisted-GPS). These systems typically provide an absolute location of the receiving device. For example, a GPS receiver can translate received GPS signals into longitudes and latitudes.


Such navigation systems, although highly precise, are often slightly inaccurate. For example, a GPS signal may provide a location that is several metres away from the correct location.


Additionally, although such systems typically provide a location of the receiver, it would be desirable to be able to use these systems to produce an absolute bearing, e.g. a bearing relative to true north or true east (90 degrees clockwise of true north). Knowledge of an absolute bearing can be particularly useful in communications systems. For example, knowing the absolute bearing of a directional antenna may have make it possible to estimate which devices are likely to receive a transmission.


Viewed from a first example configuration, there is provided an apparatus comprising: a location receiver to receive a signal indicative of a location of the location receiver; a movement mechanism to rotate the location receiver about an axis and to provide an angle of the location receiver about the axis from a known point; and calculation circuitry to produce an output value indicative of an absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.


Viewed from a second example configuration, there is provided a method comprising the steps: receiving a signal indicative of a location of the location receiver; rotating the location receiver about an axis; providing an angle of the location receiver about the axis from a known point; and producing an output value indicative of an absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.


Viewed from a third example configuration, there is provided An apparatus comprising: means for receiving a signal indicative of a location of the location receiver; means for rotating the location receiver about an axis; means for providing an angle of the location receiver about the axis from a known point; and means for producing an output value indicative of an absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.





BRIEF DESCRIPTION OF THE DRAWINGS

The present technique will be described further, by way of example only, with reference to embodiments thereof as illustrated in the accompanying drawings, in which:



FIG. 1 illustrates an example GPS error from a stationary receiver;



FIG. 2 illustrates a power spectrum density estimate of the position error shown in FIG. 1;



FIG. 3 schematically illustrates an apparatus in accordance with one embodiment;



FIG. 4 shows calculation circuitry in accordance with one embodiment;



FIG. 5 illustrates characteristics of a suitable low-pass filter for use in one embodiment;



FIG. 6 shows a flow chart illustrating an embodiment of determining a bearing in accordance with one embodiment;



FIG. 7 illustrates a power spectrum density estimate of the position error obtained when using an apparatus in accordance with one embodiment;



FIG. 8 is a graph showing how the bearing estimation error varies depending on the number of samples taken;



FIG. 9 is an example illustrating the difference between bearing and tilt;



FIG. 10 schematically illustrates calculation circuitry in accordance with one embodiment;



FIGS. 11A, 11B, 11C and 11D are graphs showing how bearing and tilt estimation errors vary depending on the number of samples taken; and



FIG. 12 is a diagram illustrating an example of an apparatus comprising an antenna array in accordance with one embodiment;





Before discussing the embodiments with reference to the accompanying figures, the following description of embodiments and associated advantages is provided.


In accordance with one example configuration there is provided an apparatus comprising: a location receiver to receive a signal indicative of a location of the location receiver; a movement mechanism to rotate the location receiver about an axis and to provide an angle of the location receiver about the axis from a known point; and calculation circuitry to produce an output value indicative of an absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.


In accordance with the above, a location receiver is rotated (e.g. in a circle) about an axis. At two or more angles (measured in radians) of the location receiver about the axis from the known point (e.g. at two or more angular displacements from the known point), a reading is taken from the location receiver in order to arrive at an angle/reading pair. The calculation circuitry uses two or more of these pairs in order to determine an absolute bearing of the known point. By determining the absolute bearing in this manner, it is possible to compensate for inaccuracy in the location signal while determining the absolute bearing of the known point, which can be used to establish how the apparatus is oriented.


In some embodiments, the apparatus further comprises: a noise filter, wherein the signal indicative of the location of the location receiver comprises a noise component; and the noise filter is to filter out the noise component. The filter may be used to remove a noise component from the received signal that is used to determine a location of the receiver. For example, the filter may be a low-pass filter. The filter can be used to further increase the accuracy of the absolute bearing determined by the calculation circuitry. The noise filter may be implemented at a variety of different points within the apparatus.


Although the calculation circuitry may have a number of different components, in some embodiments, the calculation circuitry comprises: a correlator to correlate the angle of the location receiver about the axis from a known point and output a correlation result; and a phase angle extractor to produce the output value from the correlation result.


In some embodiments, the correlator comprises: a multiplier to multiply the location of the location receiver by e−jφ, where j is the unit imaginary number and φ is the angle of the location receiver about the axis from the known point with which the location is associated, and to output a plurality of multiplication results; and an integrate and dump filter to operate on the plurality of multiplication results. The correlator in these embodiments considers an angle/reading pair (the reading taken at a particular angle from the known point about the axis). The particular angle is denoted by φ. A particular reading is multiplied by the conjugate, i.e. e−jφ, where j is the unit imaginary number (often referred to as i, or √(−1)). This produces a single multiplication result. By providing multiple angle/reading pairs, a plurality of multiplication results are produced. These are provided to an integrate and dump filter. The integrate and dump filter may, for example, average the provided values and accordingly, eliminate or reduce any noise component inherent in the location value obtained from the location receiver. The result from the integrate and dump filter is the output value and this is indicative of an absolute bearing of the known point about the axis.


The absolute bearing may be represented in a number of different ways. For example, in some embodiments, the absolute bearing of the known point about the axis is a bearing of the known point about the axis relative to true east. The term “true east” is used here to refer to 90 degrees clockwise from true north. For example, true east can be thought of as being along the line of zero latitude. In other embodiments, the absolute bearing may be a bearing relative to true north.


The received signal may be indicative of a location of the location receiver in a number of different ways. For example, in some embodiments, the signal is indicative of a longitude and latitude of the location receiver.


In some embodiments, the signal is a navigation system signal. For example, the signal may be a satellite based global navigation system signal that allows a receiver to estimate or determine its location based on the time taken for signals to travel from one or more satellites to the receiver. These timings may then be used to narrow down the possible locations of the receiver, for example, into an average or expected longitude and latitude. An example of such a satellite based global navigation system signal is a GPS signal. Other systems may include GLONASS (The Russian Global Navigation Satellite System), the planned European Union Galileo positioning system, India's Indian Regional Navigation Satellite System and the Chinese BeiDou Navigation Satellite System.


In some embodiments, the apparatus receives a further signal indicative of a location of the location receiver; and the further signal is a ground based global navigation system signal. Such systems may be used in order to further refine or aid the determination of the location of the location receiver. For example, assisted-GPS (A-GPS) is a ground-based system that can be used to improve the startup performance of a GPS receiver in order for a location receiver to determine its location more quickly.


In some embodiments, the apparatus further comprises: an antenna assembly to perform at least one of receiving or transmitting a further signal, wherein the antenna assembly defines a reference plane; the movement mechanism is to rotate the antenna assembly about the axis; the location receiver is directly or indirectly mounted to the antenna assembly at a distance from the axis equal to the radius such that when the movement mechanism rotates the antenna assembly about the axis, the location receiver rotates about the axis. In such embodiments, the location receiver is directly or indirectly mounted on the antenna assembly such that as the antenna assembly rotates, the location receiver rotates as well. Accordingly, once the absolute bearing of the known point has been determined, since the rotation of the location receiver/antenna assembly around the known point at the axis is known, it is possible to determine the absolute bearing of a reference plane of the antenna assembly. This can be useful where the antenna assembly prominently transmits or receives in a particular direction, in order to help orient the antenna assembly in order to improve reception and/or transmission.


In some embodiments, the apparatus further comprises: an antenna bearing calculation unit to determine an indication of an absolute bearing of the reference plane based on an absolute bearing of the known point about the axis, an angle of the location receiver about the axis from a known point and a position at which the location receiver is directly or indirectly mounted to the antenna assembly. The antenna bearing calculation unit may be used to determine the absolute bearing of the reference plane at a particular time.


There are a number of ways of representing an indication of the absolute bearing of the reference plane. However, in some embodiments, the indication of the absolute bearing of the reference plane is a vector perpendicular to the reference plane.


In some embodiments, the axis is perpendicular to the reference plane; and the axis passes through a mid-point of the antenna assembly. In this manner, the antenna assembly may rotate on the spot. For example, the orientation of the rotation assembly changes without changing the position of the assembly. Conversely, the location receiver continues to rotate at a radius, thereby altering its position as it rotates.


In some embodiments, the axis is upright. For example, the axis may be substantially perpendicular to the plane of the ground.


In some embodiments, the calculation circuitry is to produce the output value and a further output value indicative of a tilt of the location receiver, additionally based on a radius between the location receiver and the axis. The tilt of the location receiver may be defined as, for example, the angle between the horizontal plane (e.g. the plane tangent to the Earth's surface at the location of the apparatus) and the plane formed by the rotation of the location receiver about the axis. Calculating both the tilt and the absolute bearing can be achieved by additionally considering the radius between the location receiver and the axis (in addition to the plurality of given angles of the location receiver about the axis from the known point and the plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles).


There are a number of ways in which the tilt can be determined. In some embodiments, the calculation circuitry comprises cost minimising circuitry to calculate the output value and further output value by performing a cost minimising function. For example, the cost minimising function can be based on a least squares estimate of the output value and the further output value. For example, given a cost function C(absolute bearing, tilt), the cost minimising function may use a least squares estimate in order to seek the absolute bearing and tilt values that best fit the received location data from the location receiver.


In some other embodiments, the calculation circuitry comprises: a correlator to correlate the angle of the location receiver about the axis from a known point and output a correlation result; a phase angle extractor to produce the output value from the correlation result; and a tilt estimator to produce the further output value based on a magnitude of the correlation result.


Particular embodiments will now be described with reference to the figures.



FIG. 1 illustrates the error when a received GPS signal is translated into a longitude and latitude. A GPS receiver was placed in a stationary location. The graph shown in FIG. 1 is a trace that indicates, over time, how the error in calculated position changes. The trace shows that, in this case, the error in the y axis deviated from about 5 metres to −5 metres and the error in the x axis deviated from about −7 metres to just over 15 metres. Hence, it can be seen that the position determined from a GPS signal may be inaccurate. Although this figure relates to the error encountered when using a received GPS signal, it will be appreciated that similar errors may be encountered with other navigation system signals such as satellite based global navigation system signals. Ordinarily, GPS coordinates are given as a longitude ({tilde over (α)}k) and latitude ({tilde over (β)}k) for a reading number k. Small angle approximation can be used to convert longitude and latitude into a Cartesian coordinate system measured in metres where the centre of the Cartesian coordinate system is the average of the longitude and latitude measurements. Specifically, assuming the Earth has a constant radius of 6371000 metres, the xk and yk coordinates of a reading number k (out of K total readings) are:







x
k

=


(



β
~

k

-


K

-
1








k


=
1

K




β
~


k






)



π6371000
180



cos


(


π

180

K








k


=
1

K




α
~


k





)










y
k

=


(



α
~


k



-


K

-
1








k


=
1

K




α
~


k






)



π6371000
180







FIG. 2 illustrates a power spectrum density graph of the position error shown in FIG. 1. The profile of the graph shown in FIG. 2 is not flat, thereby indicating that the position error cannot be modelled as Additive White Gaussian Noise. A different process must therefore be used in order to overcome the noise component associated with the location estimation error.



FIG. 3 schematically illustrates an apparatus 100 in accordance with one embodiment. The apparatus comprises a location receiver 110. The location receiver is able to receive a signal that indicates the location of the location receiver 110. In the embodiment shown in FIG. 3, the signal is a GPS signal. However, the signal may be a navigation system signal such as a satellite based global navigation system signal. In some embodiments, the location receiver 110 is able to receive a further signal indicative of the location of the location receiver, from a ground based global navigation system. This can be used to further refine the determined location of the location receiver or may be used to assist in the determination of the location of the location receiver using the signal.


The location receiver 110 may be made to rotate about an axis 130 using a movement mechanism. The location receiver rotates at a radius r from the axis 130. In this embodiment, the location receiver 110 is directly attached to the movement mechanism 120 itself, which rotates about the axis 130. The location receiver 110 is therefore capable of clockwise or anticlockwise movement about axis 130. The location receiver 110 makes an angle with a known point 140 and the movement mechanism 120 is able to determine this angle. In FIG. 3, for example, the location receiver 110 makes a π/2 angle with the known point 140 at the axis 130. As the location receiver 110 rotates, it receives a plurality of signals in order to provide a plurality of locations. The calculation circuitry 150 is able to determine the absolute bearing of the known point 140 based on the plurality of locations and the plurality of angles made between the location receiver 110, the axis 130, and the known point 140 when each location was determined. For example, if the location receiver determines a first location L1 at π/2 radians and a second location L2 at π radians then these pairings are used by the calculation circuitry 150 in order to determine the absolute bearing of the known point 140. The absolute bearing, in this embodiment, is relative to true east (i.e. 90 degrees clockwise of true north at a latitude of 0). In other embodiments, other points of reference may be used in order to determine the absolute bearing.


Although the above apparatus describes the use of a GPS signal, the location receiver 110 may also use other satellite based systems. For example, the location receiver 110 may operating using a GLONASS (The Russian Global Navigation Satellite System) signal. The apparatus may alternatively use a ground-based navigation signal, or may use such a signal in addition to a satellite based system in order to assist the operation of or refine the estimate produced by, a satellite based system.



FIG. 4 schematically illustrates an example of the calculation circuitry 150 in accordance with one embodiment. The calculation circuitry comprises a correlator 240 and a phase angle extractor 250. The correlator 240 comprises a multiplier 210. The multiplier 210 receives location reading number k, Pk based on the signal received from the location receiver 110. Pk is a complex number, having a real component and an imaginary component with the real component corresponding with an estimate of the location along the x axis and the imaginary component corresponding with an estimate of the location along the y axis. Pk may be thought of as having the following relationship:

Pk=rej(θ+φk)+nk


Where r is the radius that the location receiver 110 makes with the axis 130, j is the unit imaginary number (i.e. √(−1)), θ is the absolute bearing of the known point 140 about the axis 130, φk is the bearing of the location receiver 110 from the known point 140 about the axis 130 at reading k, and nk is the noise component at reading k. This equation represents a conversion of a location reading into polar form. It also recognises that the location of the location receiver about the axis 130 is a function of absolute bearing of the known point 140 and the angle that the location receiver makes with the axis 130 and the known point 140.


The multiplier 210 multiplies Pk by the conjugate e−j(φk). This has the effect of cancelling out ej(φk):

Pk×e−j(φk)=rej(θ)+nke−j(φk)


Accordingly, the result of the multiplication is a function of θ, i.e. the absolute bearing of the known point 140 about the axis 130 and a modified noise component nke−j(φk).


In this embodiment, the multiplication result passes through a low pass filter 220. This filter is strictly optional and in other embodiments, the multiplication result may pass directly to the integrate and dump filter 230. In this embodiment, the low-pass filter performs a filtering process such that the noise component is reduced.


In any event, the multiplication result or filtered multiplication result passes to the integrate and dump filter 230. This determines a total sum (T) of all readings in order to provide an average value. The averaging performed by the integrate and dump filter 230 has the effect of reducing or eliminating the noise component, thereby producing the value rej(θ). From this, an angle estimation {circumflex over (θ)} is produced using a phase angle extractor 250.


In this embodiment, the angle estimate corresponds with the absolute bearing of the known point 140 relative to true east. In other words, an angle estimate of 0 means that the bearing of the known point about the axis 130 points true east (where true east is π/2 radians clockwise of true north, along a line of 0 latitude).



FIG. 5 illustrates an example configuration of suitable characteristics for the optional low-pass filter 220. It is clear from FIG. 5 that frequencies in the signal with a low frequency are allowed to pass whereas higher frequency components of the signal are dampened.



FIG. 6 shows a flow chart, which illustrates the procedure for collecting a plurality of location values as the location receiver 110 rotates about the axis 130. A total of N×(2π/A)×M readings are taken, in which N is the number of rotations performed by the location receiver 110 about the axis 130, A is the scanning angular resolution (e.g. the number of radians between each angular step), and M is the number of measurements per angular step.


The process starts at step 305 and at step 310, the variables N, A, and M are set to the desired values as specified above. At step 315, the quality of the GPS measurements is evaluated and, if necessary, the values of N, A, and M are changed. For example, if the GPS signal is particularly poor then the values of N, A, and M are altered so that more readings are taken. At step 320, the location receiver is rotated counter clockwise until a mechanical end stop is reached. In other words, the location receiver is rotated counter clockwise as far as it will go. Once this end point is reached, the orientation of the location receiver 110 at this position is assigned an angle. In this embodiment, the assigned angle is −188 degrees because at the position of the mechanical stop, the location receiver 110 will be −188 degrees from the known point about the axis 130. Note that the mechanical stop may be in another position, provided the position is known in advance. At step 325, the antenna is then rotated by 8 degrees until it is at an angle of −180 degrees from the known point. At step 330, the variable n (the current rotation number) is initialised to 0. At step 335, it is determined whether or not n (the current rotation number) is less than N (the number of rotations to be performed). If not, then at step 340, the readings can be used to calculate the bearing offset, and the process completes at step 345. Alternative, if n is not less than N then at step 350, n is incremented by 1. The value a (the current angle of rotation) is initialised to −180 degrees at step 355. At step 360, it is determined whether or not a (the current angle of rotation) is greater than 180 degrees. This represents the situation in which the location receiver has completed an entire rotation from −180 degrees to 180 degrees. If this is the case, then the process returns to step 335. Otherwise, the process proceeds to step 365 where m (the number of readings taken at this position) is initialised to 0. At step 370, m is incremented by one. At step 375, the longitude and latitude are determined from the GPS signal received by the location receiver 110. At step 380, it is determined whether or not m (the number of readings taken at this position) is less than M (the number of readings to take at each position). If so, then the process returns to step 370 where the value of m is incremented and another reading is taken. If not, then the process proceeds to step 385 where the value of a (the current angle of rotation) is incremented by the scanning resolution (A). At step 390, the location receiver 110 is rotated by A degrees and flow proceeds to step 360.


In effect, the flowchart shown in FIG. 6 represents a 3-nested loop. The outer-most loop (steps 330, 335, and 350) iterates through the number of rotations of the location receiver 110. Within each of those iterations, a second loop (steps 355, 360, and 385) iterates through each position at which a reading is to be taken. At each position, an inner-most loop (steps 365, 370, and 380) iterates through the number of readings to be taken at each position. The process therefore collects a set of readings (the number of which is N×(2π/A)×M), which is passed through the multiplier 210 and integrate and dump filter 230 as illustrated previously with reference to FIG. 4.



FIG. 7 illustrates a power spectrum density graph for the measurements taken using the process described with reference to FIG. 6 when an angular resolution (A) of 45 degrees is used. As can be seen, a corresponding spectral peak appears at (0.25 π radians/sample), which corresponds with A. This suggests that the use of the correlator 240 that uses a multiplier 210 followed by an integrate and dump filter 230 will yield an accurate estimate of the value of θ (i.e. the bearing that the location receiver 110 makes with the known point 140 about the axis 130. Note that if the angular resolution was π radians then the spectral peak would appear at 1×π radians per sample.



FIG. 8 illustrates a bearing estimation error for 96 readings (255), 200 readings (260), 400 readings (265), 800 readings (270), 1600 readings (275), and 3200 (280) readings. As shown, 90% of the time the bearing estimate is better than ±31.6, ±14.7, ±7.8, ±4.4, ±2.6 degrees, for 200, 400, 800, 1600, 3200 number of readings, respectively.


Using a similar methodology to that described above, it is possible to determine both the absolute bearing of a known point 140 and also the tilt of the location receiver relative to the plane tangent to the Earth's surface at the location of the location receiver 110. The different angles involved can be seen with reference to FIG. 9. In the figure, H represents the horizontal plane, which is tangent to the Earth's surface at the location of the location receiver 110. In this example, A is a plane made by the rotation of the location receiver 110 around the axis 130. The rotation of the location receiver 110 around the axis 130 is shown by a dotted line. The angle of tilt θ1 is the angle between these two planes. Meanwhile, the absolute bearing of the known point 140 about the axis 130 is given by θ0. In this example, the absolute bearing is given with reference to true east. The solid circle illustrates the projection of the path taken by the location receiver 110 on to the plane H. This forms an ellipse.


There are a number of ways of determining both θ0 and θ1. In one embodiment, cost minimising circuitry may be used to perform a cost minimising function based on θ0 and θ1. For example, consider the following quadratic cost function:







C


(


θ
0

,

θ
1


)


=




k
=
0


k
-
1








P
k

-



x
k



(

θ
1

)




e

j






θ
0







2






Where k is the number of readings, Pk is location reading number k (as previously discussed), j is the unit imaginary number (√(−1)), φk is the angle made by the location receiver 110 about the axis 130 at reading k and xk1)=r[cos(φk)+j sin(φk)cos(θ1)] denotes the trajectory of the ellipse. The aim is to minimise C(θ01) with respect to θ0 and θ1. In other words, this method seeks to compute the absolute bearing and tilt estimates that best fit the received location data.


Minimizing C(θ01) with respect to θ0 yields:









θ
^

0



(

θ
1

)


=











k
=
0


k
-
1





P
k




x
k
*



(

θ
1

)









Where ∠z denotes the angle of the complex number z and * refers to the complex conjugate. In other words, xk*(θ1)=r[cos(φk)−j sin(φk)cos(θ1)]. In the above equation, it can be seen that {circumflex over (θ)}0 is a function of θ1. Substituting {circumflex over (θ)}01) into C(θ01) yields:







C


(




θ
^

0



(

θ
1

)


,

θ
1


)


=




k
=
0


k
-
1








P
k

-



x
k



(

θ
1

)




e

j











k
=
0


k
-
1





P
k




x
k
*



(

θ
1

)










2






The least squared tilt estimate {circumflex over (θ)}1 is evaluated by firstly finding the minimum of C({circumflex over (θ)}01),θ1) and then the least square estimate of the absolute bearing angle {circumflex over (θ)}0 is derived via back substitution. Accordingly, the value of θ1 for which C({circumflex over (θ)}01),θ1) attains its minimum (arg min) is computed by:








θ
^

1

=

arg







min

θ
1







k
=
0


k
-
1








P
k

-



x
k



(

θ
1

)




e

j











k
=
0


k
-
1





P
k




x
k
*



(

θ
1

)










2








And the least squared estimate of θ0 is given by:








θ
^

0

=












k
=
0


k
-
1





P
k




x
k
*



(


θ
^

1

)









Another way to determine both the absolute bearing and tilt is schematically illustrated by the calculation circuitry 400 shown in FIG. 10. In this circuitry, after the accumulate the dump filter 230, the result z is passed to bearing estimation circuitry 250 to determine the absolute bearing as previously discussed, as well as tilt estimation circuitry 420, which is used to determine the tilt. The tilt estimation circuitry determines the tilt based on a function of the magnitude of z. In particular, as previously discussed:






z
=




k
=
0


k
-
1





P
k



e


-
j







φ
k









Assuming that φk is sampled uniformly in a circle then by performing trigonometric and algebraic manipulation, it can be shown that:






z
=



k
2




re

j






θ
0





(

1
+

cos


(

θ
1

)



)



+
n





Where n again represents a noise component. The bearing and tilt estimates are therefore given by:








θ
^

0

=







z









θ
^

1

=

{






cos

-
1




(



2
rk




z



-
1

)


,





if







z



<

rk
2







0
,



otherwise








This method may be less computationally expensive than using the previously described cost minimising circuitry.



FIGS. 11A-11D illustrate the quality of the angle estimates as a function of the number of samples (k) using each of the two methods for determining absolute bearing and tilt. In particular, the figures show the pointing error, in degrees, as a cumulative density function for 96 samples (450), 200 samples (455), 400 samples (460), 800 samples (465), 1600 samples (470), 3200 samples (475), and 6400 samples (480). FIG. 11A illustrates the absolute bearing error when using the least squares (cost minimisation) method. FIG. 11B illustrates the absolute bearing error when using the correlation method (e.g. when using the bearing estimation circuitry and tilt estimation circuitry illustrated in the embodiment of FIG. 10). FIG. 11C illustrates the tilt error when using the least squares (cost minimisation) method and FIG. 11D illustrates the tilt error when using the correlation method (e.g. when using the bearing estimation circuitry and tilt estimation circuitry illustrated in the embodiment of FIG. 10). The figures show that getting absolute bearing and tilt estimates using the computationally simpler correlation method still yields satisfactory results.



FIG. 12 illustrates an embodiment in which the apparatus 500 additionally comprises an antenna array 510 defining a reference plane 520. The axis 130 is a longitudinal axis of the antenna array 510 and the rotation mechanism 120 rotates the antenna array 510 about this axis. The location receiver 110 is directly attached to the antenna array 520 at a position that is offset from the axis 130 such that rotation of the antenna array 510 about its axis 130 causes the location receiver 110 to rotate about the axis 130 at a radius. Calculation circuitry 150 and an antenna bearing calculation unit 530 are included within the main assembly. The calculation unit 150 determines the absolute bearing and/or tilt of the apparatus 500 as previously discussed. The antenna bearing calculation unit 530 can be used to determine the absolute bearing of the reference plane 520 based on an absolute bearing of the known point 140 about the axis 130, and an angle of the location receiver 110 about the axis 130 from the known point 140 together with knowledge of a position at which the location receiver 110 is directly or indirectly mounted to the antenna assembly 510. For example, if it is known that the known point 140 is 45 degrees clockwise of true east and if it is known that the location receiver 110 is 45 degrees clockwise of the known point 140 then it is known that the location receiver is 90 degrees clockwise of true east. If it is also known that the location receiver is mounted as illustrated in FIG. 12, then it can be known that the direction perpendicular to the reference plane 520 (illustrated by the dashed line in FIG. 12) is 90 degrees clockwise of the angle of the location receiver 110, in other words, that the reference plane faces 180 degrees from true east (i.e. true west). If the location receiver 110 were mounted on the other side of the antenna array 510 at 110′ then it would be known that the reference plane 520 faces an angle 90 degrees anti-clockwise of the angle of the location receiver 110. Accordingly, the true bearing of the antenna array 510 can be determined. It will be appreciated that a similar technique may also be used to determine the tilt of the antenna array 510.


In the present application, the words “configured to . . . ” are used to mean that an element of an apparatus has a configuration able to carry out the defined operation. In this context, a “configuration” means an arrangement or manner of interconnection of hardware or software. For example, the apparatus may have dedicated hardware which provides the defined operation, or a processor or other processing device may be programmed to perform the function. “Configured to” does not imply that the apparatus element needs to be changed in any way in order to provide the defined operation.


Although illustrative embodiments of the invention have been described in detail herein with reference to the accompanying drawings, it is to be understood that the invention is not limited to those precise embodiments, and that various changes, additions and modifications can be effected therein by one skilled in the art without departing from the scope and spirit of the invention as defined by the appended claims. For example, various combinations of the features of the dependent claims could be made with the features of the independent claims without departing from the scope of the present invention.

Claims
  • 1. An apparatus comprising: a location receiver to receive a signal indicative of a location of the location receiver;a movement mechanism to rotate the location receiver about an axis and to provide an angle of the location receiver about the axis from a known point of the apparatus, wherein an absolute bearing of the known point about the axis is initially unknown; andcalculation circuitry to calculate an output value indicative of the absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.
  • 2. An apparatus according to claim 1, the apparatus further comprising: a noise filter, whereinthe signal indicative of the location of the location receiver comprises a noise component; andthe noise filter is to filter out the noise component.
  • 3. An apparatus according to claim 2, wherein the noise filter is a low-pass filter.
  • 4. An apparatus according to claim 1, wherein the calculation circuitry comprises: a correlator to correlate the angle of the location receiver about the axis from the known point and output a correlation result; anda phase angle extractor to produce the output value from the correlation result.
  • 5. An apparatus according to claim 4, wherein the correlator comprises: a multiplier to multiply the location of the location receiver by e−jφ, where j is the unit imaginary number and φ is the angle of the location receiver about the axis from the known point with which the location is associated, and to output a plurality of multiplication results; andan integrate and dump filter to operate on the plurality of multiplication results.
  • 6. An apparatus according to claim 1, wherein the absolute bearing of the known point about the axis is a bearing of the known point about the axis relative to true east.
  • 7. An apparatus according to claim 6, wherein the signal is a navigation system signal.
  • 8. An apparatus according to claim 7, wherein the signal is a satellite based global navigation system signal.
  • 9. An apparatus according to claim 8, wherein the signal is a GPS signal.
  • 10. An apparatus according to claim 8, wherein the location receiver is to receive a further signal indicative of a location of the location receiver; andthe further signal is a ground based global navigation system signal.
  • 11. An apparatus according to claim 1, wherein the signal is indicative of a longitude and latitude of the location receiver.
  • 12. An apparatus according to claim 1, the apparatus further comprising: an antenna assembly to perform at least one of receiving or transmitting a further signal, whereinthe antenna assembly defines a reference plane;the movement mechanism is to rotate the antenna assembly about the axis;the location receiver is directly or indirectly mounted to the antenna assembly at a distance from the axis equal to the radius such that when the movement mechanism rotates the antenna assembly about the axis, the location receiver rotates about the axis.
  • 13. An apparatus according to claim 12, the apparatus further comprising: an antenna bearing calculation unit to determine an indication of an absolute bearing of the reference plane based on an absolute bearing of the known point about the axis, an angle of the location receiver about the axis from a known point and a position at which the location receiver is directly or indirectly mounted to the antenna assembly.
  • 14. An apparatus according to claim 13, wherein the indication of the absolute bearing of the reference plane is a vector perpendicular to the reference plane.
  • 15. An apparatus according to claim 12, wherein the axis is perpendicular to the reference plane; andthe axis passes through a mid-point of the antenna assembly.
  • 16. An apparatus according to claim 1, wherein the axis is upright.
  • 17. An apparatus according to claim 1, wherein the calculation circuitry is to produce the output value and a further output value indicative of a tilt of the location receiver, additionally based on a radius between the location receiver and the axis.
  • 18. An apparatus according to claim 17, wherein the calculation circuitry comprises cost minimising circuitry to calculate the output value and further output value by performing a cost minimising function.
  • 19. An apparatus according to claim 18, wherein the cost minimising function is based on a least squares estimate of the output value and the further output value.
  • 20. An apparatus according to claim 17, wherein the calculation circuitry comprises: a correlator to correlate the angle of the location receiver about the axis from a known point and output a correlation result;a phase angle extractor to produce the output value from the correlation result; anda tilt estimator to produce the further output value based on a magnitude of the correlation result.
  • 21. A method of operating an apparatus, the method comprising the steps: receiving a signal indicative of a location of a location receiver;rotating the location receiver about an axis;providing an angle of the location receiver about the axis from a known point of the apparatus, wherein an absolute bearing of the known point about the axis is initially unknown; andcalculating an output value indicative of the absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.
  • 22. An apparatus comprising: means for receiving a signal indicative of a location of the location receiver;means for rotating the location receiver about an axis;means for providing an angle of the location receiver about the axis from a known point of the apparatus, wherein an absolute bearing of the known point about the axis is initially unknown; andmeans for calculating an output value indicative of an absolute bearing of the known point about the axis, based on a plurality of given angles of the location receiver about the axis from the known point, and a plurality of associated locations of the location receiver indicated by the signal received when the location receiver is at each of the given angles.
Priority Claims (2)
Number Date Country Kind
1511200.6 Jun 2015 GB national
1514938.8 Aug 2015 GB national
US Referenced Citations (107)
Number Name Date Kind
3789415 Vickland Jan 1974 A
3844341 Bimshas, Jr. et al. Oct 1974 A
4436190 Wentzell Mar 1984 A
4599620 Evans Jul 1986 A
4633256 Chadwick Dec 1986 A
4959653 Ganz Sep 1990 A
5049891 Ettinger et al. Sep 1991 A
5125008 Trawick et al. Jun 1992 A
5357259 Nosal Oct 1994 A
6124832 Jeon et al. Sep 2000 A
6404385 Croq et al. Jun 2002 B1
6448930 Judd Sep 2002 B1
6486832 Abramov et al. Nov 2002 B1
6522898 Kohno et al. Feb 2003 B1
6621454 Reudink et al. Sep 2003 B1
6784830 Lawrence Aug 2004 B1
6934511 Lovinggood et al. Aug 2005 B1
6963747 Elliott Nov 2005 B1
7062294 Rogard et al. Jun 2006 B1
7403748 Keskitalo et al. Jul 2008 B1
7515916 Alexander Apr 2009 B1
7593693 Kasapi et al. Sep 2009 B1
7664534 Johnson Feb 2010 B1
7697626 Wang et al. Apr 2010 B2
8340580 Epstein Dec 2012 B1
8509724 D'Amico et al. Aug 2013 B2
8577416 Nandagopalan et al. Nov 2013 B2
8630267 Jin Jan 2014 B1
9173064 Spain, Jr. Oct 2015 B1
9179360 Vivanco Nov 2015 B1
9692124 Caimi et al. Jun 2017 B2
9698891 Larsson Jul 2017 B2
9706419 Bozier et al. Jul 2017 B2
9924385 Lysejko et al. Mar 2018 B2
9973943 Lysejko et al. May 2018 B2
20020042274 Ades Apr 2002 A1
20020142779 Goto et al. Oct 2002 A1
20020159409 Wolfe et al. Oct 2002 A1
20030195017 Chen et al. Oct 2003 A1
20030228857 Maeki Dec 2003 A1
20040043745 Najarian Mar 2004 A1
20040077354 Jason et al. Apr 2004 A1
20040106436 Ochi et al. Jun 2004 A1
20040204850 MacNeille Oct 2004 A1
20040233103 Toshev Nov 2004 A1
20040242274 Corbett et al. Dec 2004 A1
20050048921 Chung Mar 2005 A1
20050063340 Hoffmann et al. Mar 2005 A1
20050157684 Ylitalo et al. Jul 2005 A1
20050192037 Nanda et al. Sep 2005 A1
20050285784 Chiang et al. Dec 2005 A1
20060072518 Pan et al. Apr 2006 A1
20060292991 Abramov et al. Dec 2006 A1
20080005121 Lam et al. Jan 2008 A1
20080123589 Lee et al. May 2008 A1
20080287068 Etemad Nov 2008 A1
20090005121 Wong et al. Jan 2009 A1
20090032223 Zimmerman et al. Feb 2009 A1
20090046638 Rappaport et al. Feb 2009 A1
20090067333 Ergen et al. Mar 2009 A1
20090086864 Komninakis et al. Apr 2009 A1
20090103492 Altshuller et al. Apr 2009 A1
20090116444 Wang et al. May 2009 A1
20090207077 Hwang et al. Aug 2009 A1
20090252088 Rao et al. Oct 2009 A1
20090274076 Muharemovic et al. Nov 2009 A1
20090310554 Sun et al. Dec 2009 A1
20100071049 Bahr et al. Mar 2010 A1
20100130150 D'Amico et al. May 2010 A1
20100216477 Ryan Aug 2010 A1
20100240380 Yim et al. Sep 2010 A1
20110003554 Sekiya Jan 2011 A1
20110163905 Denis et al. Jul 2011 A1
20110235569 Huang et al. Sep 2011 A1
20110244808 Shiotsuki et al. Oct 2011 A1
20110312269 Judd et al. Dec 2011 A1
20120002598 Seo et al. Jan 2012 A1
20120015659 Kalyani et al. Jan 2012 A1
20120046026 Chande et al. Feb 2012 A1
20120119951 Vollath May 2012 A1
20120252453 Nagaraja et al. Oct 2012 A1
20120299765 Huang et al. Nov 2012 A1
20120329511 Keisu Dec 2012 A1
20130041549 Reeve Feb 2013 A1
20130203401 Ryan et al. Aug 2013 A1
20130215844 Seol et al. Aug 2013 A1
20140256376 Weissman et al. Sep 2014 A1
20140313080 Smith et al. Oct 2014 A1
20150065164 Hoseinitabatabaei Mar 2015 A1
20150078191 Jongren et al. Mar 2015 A1
20160037550 Barabell et al. Feb 2016 A1
20160255667 Schwartz Sep 2016 A1
20160262045 Yang et al. Sep 2016 A1
20160277087 Jo et al. Sep 2016 A1
20160377695 Lysejko et al. Dec 2016 A1
20160380353 Lysejko et al. Dec 2016 A1
20160380354 Bozier et al. Dec 2016 A1
20160380355 Lysejko et al. Dec 2016 A1
20160380363 Logothetis Dec 2016 A1
20160381570 Lysejko et al. Dec 2016 A1
20160381574 Dudzinski et al. Dec 2016 A1
20160381585 Dudzinski et al. Dec 2016 A1
20160381590 Lysejko et al. Dec 2016 A1
20160381591 Lysejko et al. Dec 2016 A1
20160381698 Grinshpun et al. Dec 2016 A1
20170111181 Zhou Apr 2017 A1
20170280314 Yang et al. Sep 2017 A1
Foreign Referenced Citations (32)
Number Date Country
1284794 Feb 2001 CN
104635203 May 2015 CN
0654915 May 1995 EP
1903348 Mar 2008 EP
2113145 Nov 2009 EP
2207267 Jul 2010 EP
2391157 Nov 2011 EP
2448144 May 2012 EP
2506625 Oct 2012 EP
2538712 Dec 2012 EP
2318914 May 1998 GB
2484377 Apr 2012 GB
H05188128 Jul 1993 JP
H06188802 Jul 1994 JP
101346062 Dec 2013 KR
9426001 Nov 1994 WO
199965105 Dec 1999 WO
0152447 Jul 2001 WO
0231908 Apr 2002 WO
03096560 Nov 2003 WO
2004095764 Nov 2004 WO
2004114546 Dec 2004 WO
2005064967 Jul 2005 WO
2007010274 Jan 2007 WO
2007069809 Jun 2007 WO
2008064696 Jun 2008 WO
2008111882 Sep 2008 WO
2008151057 Dec 2008 WO
2010077790 Jul 2010 WO
2011044947 Apr 2011 WO
2011060058 May 2011 WO
2014106539 Jul 2014 WO
Non-Patent Literature Citations (45)
Entry
U.S. Notice of Allowance from U.S. Appl. No. 15/186,134 dated Mar. 14, 2017, 8 pages.
U.S. Office Action in U.S. Appl. No. 15/182,209 dated Apr. 13, 2017, 17 pages.
Doi et al., “Low-Cost Antenna Array Via Antenna Switching for High Resolution 2-D DOA Estimation,” SIPS 2013 Proc. IEEE, Oct. 16, 2013, pp. 83-88.
Jung et al., “Attitude Sensing Using a GPS Antenna on a Turntable: Experimental Tests,” NAVIGATION, J. of the Institute of Navigation, Fairfax, VA, US, vol. 51, No. 3, Dec. 1, 2004, pp. 221-230.
PCT Search Report from PCT/GB2016/051234 (corresponding to U.S. Appl. No. 15/187,570), dated Aug. 5, 2016, 13 pgs.
PCT Search Report from PCT/GB2016/051360 (corresponding to U.S. Appl. No. 15/187,900), dated Aug. 19, 2016, 12 pgs.
PCT Search Report from PCT/GB2016/051615 (corresponding to U.S. Appl. No. 15/187,515), dated Aug. 12, 2016, 13 pgs.
PCT Search Report from PCT/GB2016/051617 (corresponding to U.S. Appl. No. 15/187,616), dated Aug. 23, 2016, 11 pgs.
PCT Search Report from PCT/GB2016/051618 (corresponding to U.S. Appl. No. 15/187,602), dated Aug. 12, 2016, 12 pgs.
PCT Search Report from PCT/GB2016/051759 (corresponding to U.S. Appl. No. 15/187,680), dated Sep. 14, 2016, 15 pgs.
PCT Written Opinion from PCT/GB2016/051195, dated May 17, 2017, 5 pgs.
PCT Written Opinion from PCT/GB2016/051205, dated May 16, 2017, 5 pgs.
PCT Written Opinion from PCT/GB2016/051234, dated May 17, 2017, 10 pgs.
PCT Written Opinion from PCT/GB2016/051285, dated May 10, 2017, 5 pgs.
PCT Written Opinion from PCT/GB2016/051360, dated Jun. 1, 2017, 6 pgs.
PCT Written Opinion from PCT/GB2016/051428, dated Jun. 2, 2017, 12 pgs.
PCT Written Opinion from PCT/GB2016/051615, dated May 17, 2017, 5 pgs.
PCT Written Opinion from PCT/GB2016/051759, dated May 18, 2017, 6 pgs.
Office Action in related Case U.S. Appl. No. 15/186,134 dated Nov. 22, 2016, 18 pages.
PCT Search Report from PCT/GB2016/051428 (corresponding to U.S. Appl. No. 15/187,188), dated Dec. 13, 2016, 19 pgs.
Office Action in related Case U.S. Appl. No. 15/187,574 dated May 15, 2018, 13 pages.
U.S. Appl. No. 15/187,570—Office Action dated Feb. 27, 2018, 16 pages.
U.S. Appl. No. 15/187,602—Office Action dated Feb. 22, 2018, 19 pages.
U.S. Notice of Allowance from U.S. Appl. No. 15/187,515 dated May 4, 2018, 8 pages.
U.S. Notice of Allowance from U.S. Appl. No. 15/187,616 dated Jan. 9, 2018, 12 pages.
U.S. Notice of Allowance from U.S. Appl. No. 15/187,680 dated May 21, 2018, 8 pages.
U.S. Office Action in U.S. Appl. No. 15/187,574 dated Dec. 15, 2017, 14 pages.
UK Search Report from GB 1514938.8, dated Apr. 18, 2016, 4 pgs.
UK Search Report from GB 1516901.4, dated Mar. 18, 2016, 5 pgs.
UK Search Report from GB 1518654.7, dated Mar. 24, 2016, 3 pgs.
UK Search Report from GB 1519216.4, dated Apr. 15, 2016, 5 pgs.
UK Search Report from GB 1519220.6, dated Apr. 4, 2016, 4 pgs.
UK Search Report from GB 1519228.9, dated Apr. 29, 2016, 4 pgs.
UK Search Report from GB 1519270.1, dated Apr. 25, 2016, 5 pgs.
UK Search Report from GB 1519272.7, dated Jun. 10, 2016, 3 pgs.
UK Search Report from GB 1519273.5, dated Apr. 27, 2016, 3 pgs.
PCT Search Report from PCT/GB2016/051195, dated Jul. 25, 2016, 13 pgs.
PCT Search Report from PCT/GB2016/051205, dated Jul. 25, 2016, 13 pgs.
Office Action in related Case U.S. Appl. No. 15/187,680 dated Nov. 3, 2017, 8 pages.
PCT Search Report from PCT/GB2016/051285, dated Jul. 13, 2016, 12 pgs.
UK Search Report from GB 1519237.0, dated Jun. 10 Apr. 2016, 3 pgs.
Office Action in related Case U.S. Appl. No. 15/187,515 dated Dec. 7, 2017, 9 pages.
U.S. Notice of Allowance from U.S. Appl. No. 15/182,209 dated Nov. 22, 2017, 13 pages.
U.S. Notice of Allowance from U.S. Appl. No. 15/187,574 dated Dec. 15, 2017, 14 pages.
U.S. Office Action from U.S. Appl. No. 15/187,570 dated Jun. 7, 2019, 17 pages.
Related Publications (1)
Number Date Country
20160377695 A1 Dec 2016 US