Method of calibrating a compass heading

Abstract: A method for compensating for pitch and roll misalignment in 3-axis magnetometer-type compasses wherein the magnetometer during factory calibration is subjected to magnetic fields corresponding to at least two different latitudes, and correction coefficients for the measured x and y (or horizontal) magnetic field components are derived from signals derived from the z axis (or vertical) magnetic field sensor. The correction coefficients are permanently stored in the compass to provide a heading corrected for internal pitch and roll misalignment of the magnetometer field measuring elements relative to the true vertical that would normally cause errors in azimuth with changes in latitude. No further correction or misalignment is necessary in the field for compass use.

Description: The present invention relates to magnetometer-type compasses, in particular, methods and apparatus for correcting magnetometer compass errors arising from non-vertical disposition of the field sensing components therein relative to "down", or true gravitational vertical. In the type of magnetometer-type compasses which incorporate three orthogonally disposed providing output signals corresponding to the two horizontal components, Hx and Hy, and a signal corresponding to the vertical magnetic field component, Hz, the compass heading  is calculated according to the arc tan (Hy/Hx) and the Hz component is reserved for other positioned calculations. Such compasses as disclosed in US-A-4 109 199 typically dispose the horizontal field (Hx and Hy) sensors on a common planar structure and movable about both horizontal axes, and include a pendant weight disposed orthogonal to the plane having the Hx and Hy sensors so that the Hz component is directed radially inward towards the earth's center, and that the planar structure is disposed tangentially to the earth's surface. Such compasses are typically electronically and mechanically complex and are contained within a sealed, fluid-filled housing making accurate adjustment of the assembled elements difficult. When the horizontal field measuring components are not perfectly orthogonal to the Hz component measuring element, the heading as calculated by the equation (1) will in certain conditions, such as latitude extremes, include significant errors therein unless accurately recalibrated. Compass calibrations have typically involved field calculations wherein the platform is moved to a variety of different locations at the earth's surface and the calculated heading corrected according to errors between the known locations and the heading produced by the compass. Such calibration procedures are clearly inconvenient and frequently inaccurate. SUMMARY OF THE INVENTION The method of correcting the compass heading according to the present invention comprises the application of at least two different magnetic fields to the compass in the final steps of manufacture. Signals corresponding to the horizontal and vertical field components, Hx and Hy, and Hz, are stored for at least the two different applied magnetic fields, wherein differences between the measured horizontal field components (relative to the vertical components) for the two applied magnetic fields are processed to provide two correction coefficients which are subsequently permanently stored in the compass for inclusion in the calculation of the heading angle, The process is executed entirely within the laboratory or within the manufacturing environment and the calculated error correction coefficients are stored locally on each compass manufactured. The coefficients once calculated and stored on the compass need not be revised in the field after installation and reflect residual imbalances in pitch and roll. Furthermore, should the compass be subjected to physical shock causing the alignment of the horizontal field measuring components to change relative to the vertical field measuring component, the compass may be returned for factory recalibration without disassembling the housing or the mechanical structure. BRIEF DESCRIPTION OF THE DRAWING These and further features of the present invention will be better understood by reading the following detailed description, taken together with the drawing which shows a pictorial arrangement of the compass according to the present invention including a block diagram of the heading calculation and calibration elements. In a compass including a magnetometer assembly 20, as illustrated in the figure, which includes three orthogonal magnetic field sensing flux gates 22, 24 and 26, disposed along the x, y and z coordinates, respectively, and attached to a planar mounting surface 28 from which a pendant weight 30 is suspended, typically along the z coordinate. The planar surface 28 is gimballed to be freely moveable about the x and y coordinates. The weight 30 maintains the z axis radially directed toward the center of the earth, and the planar surface 28, generally tangential to the surface of the earth. The Hy, Hx and Hz signals as measured by the flux gates 22, 24 and 26 are received and stored by a signal sampling element 42 as part of the computational portion 40 of the compass, so that the heading angle ϑ is calculated in the element 44. In the system as described, the heading angle , if calculated according to equation (1), would be uncorrected for errors in alignment between the gravitational vertical and the horizontal sensing elements 22 and 24 on the planar surface 28. According to the present invention, the misalignment of the vertical sensing element 46 produces errors R and P as illustrated in the figure, which are reflected in the calculation of the heading angle , as illustrated by equations 2 and 3: Hx = H x , cosP + H z , sin P where H z , is the actual vertical component of the earth's field and H x , is the actual x horizontal component of the earth's field, and Hy = H y' cosR + H z' sin R where H y , is the actual y horizontal component of the earth's field. If R and P are small angles, i.e., less than 15 °, then Hz = Hz' and Hx = Hx' + Hz * P and Hy = Hy' + Hz * R Correction coefficients R and P are determined and correspond to the misalignment of the planar surface 28 as reflected in the spacings illustrated by R and P in the figure. The correction coefficients R and P are calculated by sequentially providing substantially uniform first and second magnetic fields through the area of the field sensing elements 22, 24 and 26, typically provided by Helmholz coils 50 as part of the laboratory or manufacturing assembly process. The Helmholz coils 50 are energized to provide the first and second magnetic fields to correspond to two different latitudes, wherein two different Hz' fields are provided, the signal from the horizontal component Hz is measured for the two different magnetic fields, according to equations 6 and 7, below, wherein the difference is provided by equation (8): Hx 1 = Hx' 1 + Hz' 1 * P Hx 2 = Hx' 2 + Hz' 2 * P (7 Hx 1 - Hx 2 = (Hz' 1 - Hz' 2 ) * P The difference in the Hy signal is similarly measured and the difference calculated by equation (9): Hy 1 - Hy 2 = (Hz' 1 - Hz' 2 ) * R The P and R correction coefficients are calculated as follows, reflecting the imbalance or misalignment of the magnetic field sensing elements in the compass: (Hx 1 - Hx 2 )/(Hz 1 - Hz 2 ) = dHx/dHz = P (Hy 1 - Hy 2 )/(Hz 1 - Hz 2 ) = dHy/dHz = R The resulting coefficients P and R are stored in the compass, wherein the two horizontal components, as illustrated by equations (12) and (13) are calculated and include the correction coefficients P and R.. Hxm (Hx measured) = H H cos = Hx - Hz * P Hym = H H sin = Hy - Hz * R Therefore, the correction heading of  COR is calculated according to equation (14):  COR = ATN (Hym/Hxm) = ATN [(Hy - Hz * R)/(Hx - Hz * P) or more exactly:  COR * = ATN (Hym/Hxm) = ATN[(Hy - Hz sin R)/(Hx - Hz sin P)] Also, H vector output signals can be provided according to equations (16) and (17): Hy (COR) = Hy - Hz sin R Hx (COR) = Hx - Hz sin P As illustrated in the figure, the stored values of Hx, Hy and Hz are received by the element 46 which calculates the R and P correction coefficients according to the equations described above. The resulting correction coefficients R and P are stored in a storage register 48 which comprises an EPROM, a ROM or equivalent device. The correction coefficient calculation element may be part of the calculation element 44 and resident with the magnetometer element 20. In alternate embodiments, the calculation of the correction coefficients R and P by element 46 is provided by a device external to the compass and connected only when the compass is calibrated during manufacture. As an example of the calibration procedure according to the present invention, the Helmholz coils are energized to provide two separate Hz corresponding magnetic fields of -500 mg (milligauss), +500 mg. Thus, for the calculation of P in equation (10), Hx 1 may equal 100, Hx 2 may equal 110, and Hz 1 - Hz 2 = 1000. Thus, P = dHx/dHz = .01 which corresponds to an imbalance or misalignment of .57°. Modifications and substitutions made by one of ordinary skill in the art are within the scope of the present invention as defined by the appended claims.

Claims: A method of calibrating a compass heading  for a compass comprising a magnetometer having means for sensing horizontal components of a local magnetic field as independent output signals H x and H y and means for sensing the vertical component of the local magnetic field as independent output signal H z , the method comprising the steps of: applying a first local magnetic field representative of the earth's magnetic field at a first latitude, the first local magnetic field having horizontal components and a vertical component of a first magnitude; storing the resulting H x , H y and H z first magnetometer output signals; applying a second local magnetic field representative of the earth's magnetic field at a second latitude, the second local magnetic field having a vertical component different from that of the first local magnetic field; storing the resulting H x , H y and H z second magnetometer output signals; calculating the difference between the stored first and second magnetometer output signals as dH x , dH y and dH z , calculating a first error coefficient P according to P = dH x /dH z ; calculating a second error coefficient R according to R = dH y /dH z ; computing a corrected heading  for any applied earth's magnetic fields as: = arctan [H y - H z sin R)/(H x - H z sin P)] The method of claim 1, wherein, in the step of computing a corrected heading , the heading  is instead computed as: = arctan [H y - H z · R)/(H x - H z · P)].

Inventors: FOWLER JOHN T, FOWLER, JOHN T.
Applicants: DIGICOURSE INC, DIGICOURSE, INC.
Application year: 1989
Main classification: G01C17
IPC codes: G01C17, G01C21