First, the relative maximal amplitudes
(g
Dx ,
g
Dy , g
Dz ,
g
Tx , g
Ty ,
g
Tz ) induced by the three
magnetic fields (
x,
y,
z) and measured
from the two search coils (D, direction coil; T, torsion coil) are
computed by dividing the maximal voltages
(v
Dx ,
v
Dy , v
Dz ,
v
Tx , v
Ty ,
v
Tz ) by the values induced by
the
z-field (v
Dz ,
v
Tz ):
\[\mathbf{g}_{\mathbf{D}}{=}\left(\begin{array}{l}g_{\mathrm{D}_{x}}\\g_{\mathrm{D}_{y}}\\g_{\mathrm{D}_{z}}\end{array}\right)\ {=}\left(\begin{array}{l}v_{\mathrm{D}_{x}}\\v_{\mathrm{D}_{y}}\\v_{\mathrm{D}_{z}}\end{array}\right)\ {\div}{\vert}v_{\mathrm{D}_{z}}{\vert}\ \mathbf{g}_{\mathbf{T}}{=}\left(\begin{array}{l}g_{\mathrm{T}_{x}}\\g_{\mathrm{T}_{y}}\\g_{\mathrm{T}_{z}}\end{array}\right)\ {=}\left(\begin{array}{l}v_{\mathrm{T}_{x}}\\v_{\mathrm{T}_{y}}\\v_{\mathrm{T}_{z}}\end{array}\right)\ {\div}{\vert}v_{\mathrm{T}_{z}}{\vert}\]
Each sample of the signal from the two coils during the eye
movement recording consisted of six voltages: three values from the
direction coil (s
Dx ,
s
Dy ,
s
Dz ) and three values from the
torsion coil (s
Tx ,
s
Ty ,
s
Tz ). These voltages are
divided by the field gains (
g D ) and
(
g T ):
\[\mathbf{d}{=}\left(\begin{array}{l}d_{x}\\d_{y}\\d_{z}\end{array}\right)\ {=}\left(\begin{array}{l}s_{\mathrm{D}_{x}}\\s_{\mathrm{D}_{y}}\\s_{\mathrm{D}_{z}}\end{array}\right)\ {\div}\left(\begin{array}{l}g_{\mathrm{D}_{x}}\\g_{\mathrm{D}_{y}}\\g_{\mathrm{D}_{z}}\end{array}\right)\ \mathbf{t}{=}\left(\begin{array}{l}t_{x}\\t_{y}\\t_{z}\end{array}\right)\ {=}\left(\begin{array}{l}s_{\mathrm{T}_{x}}\\s_{\mathrm{T}_{y}}\\s_{\mathrm{T}_{z}}\end{array}\right)\ {\div}\left(\begin{array}{l}g_{\mathrm{T}_{x}}\\g_{\mathrm{T}_{y}}\\g_{\mathrm{T}_{z}}\end{array}\right)\]
Because the maximal amplitudes of the induced voltages during
calibration are measured according to the right-hand rule (see Appendix I), this rule is automatically preserved by this division, independent
of the signs of the raw signals.
Each sample now consists of two vectors that are nonorthogonal, because
the effective surfaces of the search coils in the annulus usually do
not form a perfect 90° angle. Thus, the next step is to orthogonalize
the sensitivity vector of the torsion coil (
t) to the
sensitivity vector of the direction coil (
d). First, we
normalize the vector of the direction coil:
\[\mathbf{f}^{\mathbf{1}}{=}\mathbf{d}{\div}{\parallel}\mathbf{d}{\parallel}\]
Then
t is projected onto
f 1 by the scalar product:
\[\mathbf{t{^\prime}}{=}(\mathbf{t}{\cdot}\mathbf{f}^{\mathbf{1}})\mathbf{f}^{\mathbf{1}}\]
The orthogonalized vector of the torsion coil
(
t″) is calculated by the vectorial subtraction:
\[\mathbf{t{^{\prime\prime}}}{=}\mathbf{t}-\mathbf{t{^\prime}}\]
Thereafter,
t″ is normalized:
\[\mathbf{f}^{\mathbf{2}}{=}\mathbf{t{^{\prime\prime}}}{\div}{\parallel}\mathbf{t{^{\prime\prime}}}{\parallel}\]
A third sensitivity vector
(
f 3 ) orthogonal to the other two
orthogonal sensitivity vectors is computed by the cross product:
\[\mathbf{f}^{\mathbf{3}}{=}\mathbf{f}^{\mathbf{1}}{\times}\mathbf{f}^{\mathbf{2}}\]
The three orthogonal normalized vectors form the eye-fixed base
vectors (
f 1 ,
f 2 ,
f 3 ). Hence, the three-dimensional
eye position at each moment in time is represented by the matrix
f that consists of the three normalized orthogonal vectors:
\[f{=}{[}\mathbf{f}^{\mathbf{1}}\mathbf{\ f}^{\mathbf{2}}\mathbf{\ f}^{\mathbf{3}}{]}\]
Next, the rotation matrix (
R) describing the rotation
of the three orthogonal vectors from the reference position
(
ref f) to the
instantaneous position (
f) is computed:
\[\mathit{R}{=}f{\cdot}_{\mathrm{ref}}f^{\mathrm{T}}\]
For each rotation matrix (
R), the corresponding
rotation vector (
r) is extracted:
\[\mathbf{r}{=}\left(\begin{array}{l}r_{x}\\r_{y}\\r_{z}\end{array}\right)\ {=}\left(\begin{array}{l}R_{32}-R_{23}\\R_{13}-R_{31}\\R_{21}-R_{12}\end{array}\right)\ {\div}(1{+}R_{11}{+}R_{22}{+}R_{33})\]