Method of calibrating a compass heading

• US 4,843,865 A
• Filed: 02/29/1988
• Issued: 07/04/1989
• Est. Priority Date: 02/29/1988
• Status: Expired due to Term

• Current Assignee: Input/Output Incorporated
• Original Assignee: Digicourse, Inc. (ION Geophysical Corporation)
• Inventors: Fowler, John T.
• Primary Examiner(s): Noland, Tom
• Application Number:
  US07/161,915
• Time in Patent Office: 491 Days
• Field of Search: 73/1 E, 73/356
• US Class Current: 73/1.76
• CPC Class Codes: G01C 17/38 Testing, calibrating, or compensating of compasses

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.

Claims

1. A method of calibrating compass heading θ

for a compass comprising a magnetometer having means for sensing horizontal components as independent output signals Hx and Hy and means for sensing the vertical component providing a corresponding Hz output signal, the method comprising the steps of;

applying a first magnetic field representative of the earth'"'"'s magnetic field at a first latitude;

storing the Hx, Hy and Hz signals corresponding to said first magnetic field;

applying a second magnetic field representative of the earth'"'"'s magnetic field at a second latitude;

storing the Hx, Hy and Hz signals corresponding to said second magnetic field;

providing the difference between the different stored instances of Hx, Hy and Hz as dHx, dHy and dHz signals, respectively;

calculating a first error coefficient signal P according to dHx/dHz;

calculating a second error coefficient signal R according to dHy/dHz; and

adding the error coefficients P and R to the calculation of the cosine and sine of the heading angle, wherein

the corrected heading angle is computed as;

##EQU1## wherein Hxm is Hx measured,

Hym is Hy measured,

Hx'"'"' is the actual x horizontal component of the earth'"'"'s field, and

Hz'"'"' is the actual vertical component of the earth'"'"'s field.

2. The method of claim 1, wherein the R and P terms comprise sin R and cos P, respectively, to provide the exact form of computed heading angle:

##EQU2##

3. The method of claim 1, wherein the calculation is provided externally to the compass.

4. A method of calibrating compass heading θ

for a compass comprising a magnetometer having means for sensing horizontal components as independent output signals Hx and Hy and means for sensing the vertical component providing a corresponding Hz output signal, the method comprising the steps of;

applying a first magnetic field representative of the earth'"'"'s magnetic field at a first latitude;

storing the Hx, Hy and Hz signals corresponding to said first magnetic field;

applying a second magnetic field representative of the earth'"'"'s magnetic field at a second latitude;

storing the Hx, Hy and Hz signals corresponding to said second magnetic field;

providing the difference between the difference stored instances of Hx, Hy and Hz as dHx, dHy and dHz signals, respectively;

calculating a first error coefficient signal P according to dHx/dHz;

calculating a second error coefficient signal R according to dHy/dHz; and

adding the error coefficients P and R to the calculation of the cosine and sine of the heading angle producing the following H vector output signals;

space="preserve" listing-type="equation">Hx(corrected)=Hx'"'"' cos θ

+Hz'"'"' sin P

space="preserve" listing-type="equation">Hy(corrected)=Hy'"'"' cos θ

+Hz'"'"' sin R,

wherein

Hx'"'"' is the actual x horizontal component of the earth'"'"'s field,

Hy'"'"' is the actual y horizontal component of the earth'"'"'s field, and

Hz'"'"' is the actual z vertical component of the earth'"'"'s field.