JN Add DOIs to your references at manuscript stage!
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 98: 3095-3108, 2007. First published September 26, 2007; doi:10.1152/jn.00764.2007
0022-3077/07 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
98/5/3095    most recent
00764.2007v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kunin, M.
Right arrow Articles by Raphan, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kunin, M.
Right arrow Articles by Raphan, T.

Rotation Axes of the Head During Positioning, Head Shaking, and Locomotion

Mikhail Kunin1, Yasuhiro Osaki2,3, Bernard Cohen2 and Theodore Raphan1

1Institute for Neural and Intelligent Systems, Department of Computer and Information Science, City University of New York (Brooklyn College), Brooklyn, New York; 2Department of Neurology, Mount Sinai School of Medicine, New York, New York; and 3Department of Otorhinolaryngology, Osaka University School of Medicine, Osaka, Japan

Submitted 6 June 2007; accepted in final form 19 September 2007


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Static head orientations obey Donders’ law and are postulated to be rotations constrained by a Fick gimbal. Head oscillations can be voluntary or generated during natural locomotion. Whether the rotation axes of the voluntary oscillations or during locomotion are constrained by the same gimbal is unknown and is the subject of this study. Head orientation was monitored with an Optotrak (Northern Digital). Human subjects viewed visual targets wearing pin-hole goggles to achieve static head positions with the eyes centered in the orbit. Incremental rotation axes were determined for pitch and yaw by computing the velocity vectors during head oscillation and during locomotion at 1.5 m/s on a treadmill. Static head orientation could be described by a generalization of the Fick gimbal by having the axis of the second rotation rotate by a fraction, k, of the angle of the first rotation without a third rotation. We have designated this as a k-gimbal system. Incremental rotation axes for both pitch and yaw oscillations were functions of the pitch but not the yaw head positions. The pivot point for head oscillations was close to the midpoint of the interaural line. During locomotion, however, the pivot point was considerably lower. These findings are well explained by an implementation of the k-gimbal model, which has a rotation axis superimposed on a Fick-gimbal system. This could be realized physiologically by the head interface with the dens and occipital condyles during head oscillation with a contribution of the lower spine to pitch during locomotion.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Changes in head position are critical to direct gaze during natural motion. A better understanding of the coordinate frames and axes that govern head rotation both globally and incrementally would give insight into how fixation is accomplished. Understanding of head rotation axes under static and dynamic situations, however, is limited. One reason for this is that the head-neck system is a complex biomechanical linkage with multiple degrees of freedom (Ishii et al. 2004aGo,bGo), coupled with a complex musculature for moving the head on the trunk (Peterson et al. 2001Go). Despite this complexity, some simplifications have proven useful in examining the coordinate frames and axes that govern head movement.

One such simplification is to model the head as a rigid sphere rotating about its center of mass (Ceylan et al. 2000Go; Glenn and Vilis 1992Go; Stark et al. 1980Go, 1988Go), but constraints reduce the degrees of freedom of head movement (Ceylan et al. 2000Go; Crawford et al. 1999Go; Radau et al. 1994Go). A similar reduction in degrees of freedom govern eye positioning mechanisms. Although the eyes are capable of moving about roll, pitch, and yaw axes of the head in response to vestibular stimulation, the rotation axes during visual fixation lie in a plane, obeying Listing's law (Helmholz 1867Go; Tweed and Vilis 1987Go). Head orientations in three dimensions do not obey Listing's law (Glenn and Vilis 1992Go), but the axes are constrained to a curved surface by the less restrictive Donders law (Glenn and Vilis 1992Go; Haslwanter et al. 1991Go; Straumann et al. 1991Go).

Each surface shape, be it planar or twisted, has its own particular sensory and motor consequences (Crawford 1997Go; Crawford and Guitton 1997Go; DeSouza et al. 1997Go). By examining the relationship between the roll of the head in space as a function of pitch and yaw, it was inferred that the curved surface generated by this relationship corresponded to a Fick gimbal, where the roll was zero. That is, the yaw axis of rotation is fixed relative to the trunk, the pitch axis moves with the head, and there is no roll (Radau et al. 1994Go). Studies in the monkey indicate that the functional relationship between head roll and the corresponding pitch and yaw components give rise to a more planar surface, thereby obeying Listing's law (Crawford et al. 1999Go). In humans using pinhole goggles, however, the difference between the surfaces with and without using pinhole goggles was shown to be small (Ceylan et al. 2000Go). There has been considerable speculation as to why the head positioning mechanism might follow a Fick-like system. These range from optimizing perception of lines in the horizon (Glenn and Vilis 1992Go; Hore et al. 1992Go), to improving binocular alignment of the eyes with the horizon (Crawford and Vilis 1995Go), to reducing the workload on the neck muscles (Glenn and Vilis 1992Go; Radau et al. 1994Go). A satisfactory model has not been developed to explain the underlying organization of the system, however.

The physiological constraints that maintain Fick-like gimbal behavior depend on anatomical constraints. The C1inferior articular facets, which support the head on the spinal column, are circular, flattened and slightly convex and are directed downward and medialward, articulating with C2 (Gray 1918Go). The head pitches through movement of the occipital condyles on the superior facets of C1, while it rotates about the vertical bony component of the dens, the portion of C2. which protrudes through the odontoid process of C1 (Ishii et al. 2004aGo,bGo). We have recently hypothesized that the dens and condyles are critical for determining the incremental axes about which the head rotates (Kunin 2004Go). However, the cervical spine flexes, rotates, and extends the head relative to the trunk, which would move the yaw and pitch axes away from the occipital condyles and the dens (Dutia 1991Go; Graf et al. 1992Go; Ishii et al. 2004aGo; Vidal et al. 1988Go). This complex structure has suggested that the axes of gaze shifts to targets at different locations occupy different spatial locations (Medendorp et al. 1998Go). However, no model-based understanding of mathematical relationship between the rotation axes for the wide range of head movements in three dimensions has been determined.

The purpose of this study was to determine how the axes of rotation vary during voluntary and natural movements, and how these axes might be related to the anatomical structures. As noted in the preceding text, the axes of eye rotation are constrained for static positions of fixation and saccades so that they lie in Listing's plane, but they are unconstrained during head movements that activate the vestibuloocular reflex (VOR). We questioned whether there was a similar lack of constraint for oscillatory head movements that would activate the vestibular system. Our hypothesis was that the rotations of the head about roll, pitch, and yaw do not intersect at a common pivot point and that the axes would lie along lines in space that are unique to the direction about which the head rotates. We also hypothesized that these lines would vary as a function of head orientation on the body. To study this, we developed a new methodology to determine incremental rotation axes and their localizations in three dimensions using displacement matrices. The axes associated with incremental rotations were then compared with the axes of head orientation during natural locomotion with and without pinhole goggles to fixate points in the visual field binocularly. A model was then developed that shows how the constraints of the spinal column on which the head rotates explain the behavior of the axes during head positioning and locomotion.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Subjects

Three general conditions were studied in this report: head orientation, head oscillation in different orientations and head oscillation during locomotion. Six subjects were studied during both head orientation and head shaking (3 female and 3 male). Three of these subjects (3 males) as well as one additional subject (female) were studied during locomotion. Their ages ranged between 32 and 57. Subjects viewed visual targets with pinhole goggles to assume static head positions. They then made small voluntary head oscillations at each of these head orientations. Finally, their head movements were recorded during locomotion on a linear treadmill. The experiments were approved by the Institutional Review Boards of the Mount Sinai School of Medicine and the City University of New York and were performed in accordance with the ethical standards established in the 1964 Declaration of Helsinki. Subjects gave their informed consent prior to the inclusion in the study. Privacy and anonymity of the data associated with the subjects were maintained.

Data acquisition

Head and trunk movements were recorded using the OPTOTRAK 3020 video motion analysis system (Northern Digital, Ontario, Canada). Briefly, the OPTOTRAK system sends a pulse to each of the markers (light-emitting diodes, LEDs) and uses a bank of three cameras (the sensor) to triangulate the location of the LED in three dimensions. Each marker is pulsed in sequence. The Optotrak system then tracks the three-dimensional position of each LED and stores this information in a computer file. Twenty-one markers on the head, trunk, and legs were utilized in this study. This allowed measurement of the head and trunk in space in three dimensions at a sampling rate of 100 Hz both when stationary and while walking on a treadmill. The OPTOTRAK system was positioned 4 m from the subject; in this configuration, the accuracy of translational measurement was 0.3 mm (manufacturer's specification) with accuracy in pitch and yaw of 0.1° (See Hirasaki et al. 1999Go for a complete description of the OPTOTRAK characteristics.)

Marker placement and determination of head coordinate frame

Coordinate frames and rigid bodies were formed by combining groups of markers. The spatial coordinate frame to which movements of the head were referenced was obtained from a set of four markers that were embedded 24 cm apart in a square plastic plate (Fig. 1A). During head positioning and voluntary head oscillation, the plate was attached to the screen that held the visual targets. The sides of the plate were aligned with the spatial horizontal and vertical axes.


Figure 1
View larger version (39K):
[in this window]
[in a new window]

 
FIG. 1. A: spatial coordinate frame. The zS axis is aligned with gravity, the xS axis approximately coincides with the nasooccipital axis of the subject in primary position, and the yS axis completes right-handed triple (the origin was placed close to the central target). B: landmark-based head coordinate frame. The yH axis is aligned with the inter-aural axis of the head the landmarks of which are located at each external auditory meatus (EAM) (e, f). The xH axis is perpendicular to the line connecting the points e, f and passes through the nose landmark (d). The z axis completes right-handed triple and is through the top of the head. The origin of the head coordinate frame is a point on the line connecting points e, f closest to the landmark d. C: visual target array. Targets are at the points of intersection of the gaze line with the spatial y-z plane, when the gaze line is rotated according to a Fick gimbal and spaced 10° in pitch and yaw. D: head calibrator consisted of standard ear plugs with a plastic front that held a nose clip. The ear plugs had light-emitting diodes (LEDs) on each side (e, f) and 1 LED on the front of the plastic piece holding the nose clip (d). The ear plugs were inserted into each auditory meatus, and the nose clip was placed on the bridge of the nose. When worn in this manner, a plane was formed by d, e, f, which contained Reid's lines and was the basis for the head coordinate frame.

 
The spatial coordinate frame was obtained by choosing one of the LED locations as the origin (Fig. 1A, LED 1). A line connecting LEDs 1 and 2, aligned with the spatial vertical and directed upward, was the zS axis (Fig. 1A). A line connecting LEDs 1 and 4, aligned with the spatial horizontal in the direction of walking and orthogonal to the zS axis, was the yS direction (Fig. 1A). A line normal to the LED plate (dashed arrow) that formed a right-handed coordinate frame was the xS axis (Fig. 1A).

The rigid body for determining translation and rotation of the head consisted of 12 LEDs (markers) that were embedded in a lightweight headband (Fig. 1B). The head coordinate frame was determined using a head calibration device that was composed of a plastic frame that had three LEDs worn by the subject only at calibration time (Fig. 1D). The right and left LEDs were placed at the right and left external auditory meatii (EAM), and the central LED lay over the bridge of the nose. The plane formed by the LEDs was adjusted so that it intersected the inferior rim of the orbits. Thus this plane contained the Reid's baselines on either side of the head (Fig. 1B, e, f, d). The y axis (yH) was defined as the interaural axis with its positive direction from the right to the left ear (Fig. 1B, right). The origin of the head coordinate frame was chosen at the point where the LED over the nose projected onto the interaural axis. This was approximately at the midpoint of the line segment joining the points e and f (Fig. 1B, right). The positive x axis (xH) was along the line from the origin to the landmark on the nose. The positive z axis (zH) was chosen to complete the right-handed coordinate triple and is designated by the symbol {odot}.

Coordinate and rotation conventions in three dimensions

Various conventions and coordinate frames have been adapted to represent three-dimensional rotations (Goldstein 1980Go). In this study, head orientations in space are described by the Euler angles {phi}, {theta}, and {psi} in accordance with the Fick convention. This convention assigns the first angle {phi} to rotation about the spatial vertical axis, which coincides with the axis out of the top of the head. We refer to this rotation as yaw. The second angle {theta}, describes rotation about an interaural axis after the head has rotated in yaw. This rotation is pitch. The final rotation is roll by angle {psi} about the naso-occipital axis (Fig. 1C). Thus {phi}, {theta}, and {psi}, the Euler-based yaw, pitch, and roll angles, form a Fick gimbal system.

Head orientation could also be described by roll, pitch, and yaw about the spatial coordinates, xS, yS, and zS, respectively. This is the axis-angle or rotation-vector-based roll, pitch, and yaw components and has been referred to as Listing's convention. Each of these conventions leads to a unique description of head orientation in three dimensions. Whereas there is some ambiguity in the usage of roll, pitch, and yaw for these various conventions, they all follow the convention that roll is a rotation about some x axis, pitch about some y axis, and yaw about some z axis. The context will clarify the convention for the usage. The question that we address here is whether voluntary head rotations and oscillations as well as head oscillations during locomotion follow some consistent convention and whether it can be theoretically related to the structure of the spinal column and its interface to the head.

Calibration

Calibration was done by asking the subjects to stand erect on a platform while being slowly rotated for 30 s over 360°. Subjects held onto a bar connected to the rotator to maintain their stability. The coordinates of the LEDs on the head rigid body (Fig. 1B) and those associated with the landmarks were determined dynamically using a rigid-body building program (6DARCHITECT, Northern Digital). From this, the coordinates of the headband LEDs were determined with respect to a coordinate frame defined by the landmark LEDs as described in the preceding text. Following the calibration procedure and the determination of the head coordinate frame, the Reid's plane calibrator was removed, but all movements of the head could be determined relative to space or this coordinate frame.

Experimental protocol

Subjects were first asked to orient their heads by looking at an array of 49 targets placed 60 cm from the subject. The targets were arranged in accordance with a Fick gimbal projected onto a plane (Fig. 1C). The targets were 10° apart and ranged from –30° to +30° in pitch and yaw. This arrangement was based on the assumption that the head orientation system is close to a Fick gimbal. In later experiments, 25 targets were spaced 15° apart to reduce the experimental time. Subjects viewed the targets through binocular pinhole goggles. The pinholes were aligned so that the head was in the zero position when looking at the central target. This forced the subjects to align their naso-occipital axes (x axis) along their line of sight so that the head was pointed at a given target (Ceylan et al. 2000Go).

Three paradigms were executed for each of the head orientations. Subjects first fixed their heads by looking at one of the 25 targets for 3 s. These head orientations were represented by a rotation matrix or rotation vector relative to a primary orientation and were referred to as global head positions. Subjects then oscillated their heads for 3 s in yaw in each of these orientations and, on the same or a different occasion, in pitch. These head oscillations were represented by incremental rotations about the given head orientation. Head-oscillation axes were also determined during treadmill locomotion at walking velocities of 0.9, 1.2, 1.5, 1.8, 2.1, and 2.4 m/s while fixating a target at 25 cm. These involuntary head oscillations were also represented by incremental rotations during locomotion. The incremental rotations were the angular velocity vectors at those particular positions. The close viewing distance during locomotion was chosen so that the angular head velocity was maximal, which made computation of rotation axes more robust. In initial experiments, neither the goggles nor the targets were used, and subjects were directed to look into different directions. At each head orientation, they were given verbal feedback on the pitch and yaw coordinates of the head until the head orientation was correct. The data from the initial and subsequent experiments were not different and were combined in the analysis. The equivalence between the initial and subsequent experiments using the pinhole goggles verified that the goggles had produced the appropriate head positions and shortened the experimental time.

Data processing

ANALYSIS OF DATA FOR AXIS OF HEAD ROTATION. Orientations of the head were represented by rotation vectors (Haustein 1989Go) composed of 300 samples obtained from 3 s of data collection at each position. Using this convention, the relationship between the x (roll) and y-z (pitch-yaw) product was utilized to develop a methodology for determining the gimbal structure associated with the head (see Conceptual and mathematical basis).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
General properties of head orientation and oscillations

Figure 2 shows the yaw, pitch, and roll of the head as Euler angles and the character of the voluntary oscillations that were superimposed during the experimental protocol. The process by which the head orientations were achieved is considered in Head orientation in three dimensions. Oscillations of the head in yaw (Fig. 2A, top), were approximately sinusoidal and of constant amplitude in association with changes in head position in yaw (top) and pitch (middle). Roll oscillations increased when the head was positioned in extreme positions in pitch (bottom). Although the roll was on an expanded scale (note ordinate), the finding that the axis of yaw oscillations of the head developed a roll component as the head was pitched, indicates that the yaw axis of the oscillations did not move exactly with the head yaw axis although the deviation was small.


Figure 2
View larger version (14K):
[in this window]
[in a new window]

 
FIG. 2. A: yaw, pitch, and roll components of head orientation when the head is oriented in pitch and yaw positions and oscillated in yaw at these positions. B: yaw, pitch, and roll components of head orientation, when the head is oriented in pitch and yaw positions and oscillated in pitch at these positions. Data were obtained from a representative subject (MK).

 
Pitch oscillations, on the other hand, were sinusoidal and of constant amplitude (Fig. 2B, middle) regardless of the yaw or pitch head position (top and bottom), and there was little or no oscillation in yaw when the head was sinusoidally pitched. Thus the axes of rotation in pitch were close to the interaural axis of the head. To determine these axes more precisely, an algorithm was developed that computed the instantaneous axis of the head oscillation in three-dimensional space (see Incremental Head Rotation: Mathematical Basis).

Head orientation in three dimensions

CONCEPTUAL AND MATHEMATICAL BASIS. In general, three parameters are required to describe orientation of a rigid body (Goldstein 1980Go). Based on the idea that head orientation follows Donders’ law, two parameters are sufficient to describe the surface of head orientation (Ceylan et al. 2000Go). In a Fick-gimbal system the roll angle of which is zero, the parameterization of the head rotation in space can be given by a rotation about the spatial vertical, followed by a rotation about the rotated pitch axis. The roll component of the axis-angle representation of the head is determined by its pitch and yaw components, forming a two-dimensional surface. Ceylan et al. (2000)Go utilized a gimbal score to reflect this constraint using a ratio of the quaternion components

Formula 1(1)
where s is the gimbal score and q0, q1, q2, q3 are the components of the quaternion representing the rotation (Tweed and Vilis 1987Go). The gimbal score can be equivalently represented in terms of the components of a rotation vector describing the rotation (Haustein 1989Go) as

Formula 2(2)
where the rotation vector, r, is directed along the axis of the rotation and has the magnitude of the tangent of half the angle of the rotation (Haustein 1989Go). This representation allows the value of the gimbal score to be estimated given several orientations of a rigid body by using linear regression on the plots of rx versus the product ryrz. We demonstrated this model with a mechanical Fick-gimbal system where the roll was zero, i.e., a truncated Fick gimbal, on which an artificial head was mounted (Fig. 3A). The artificial head was positioned in the same orientations as during experiments and the relationship between rx versus the product ryrz was plotted (Fig. 3, B and C). For both extension (B) and flexion (C), the regression line had a slope equal to –1 and little variance (Fig. 3, B and C, circles). A similar plot for a Helmoltz gimbal, which is governed by a pitch about a spatial horizontal axis, followed by a rotation about the rotated yaw axis, would give a slope of +1 (Fig. 3, B and C, diamonds). A Listing's law mechanism, which is a rotation from primary position about axes that lie in a plane would give a slope of 0 (Fig. 3, B and C, squares) because roll is zero.


Figure 3
View larger version (18K):
[in this window]
[in a new window]

 
FIG. 3. A: head orientations according to a truncated Fick-gimbal model, i.e., where roll is equal to zero. B and C: for Helmholtz and Listing's models—simulations of roll (rx) vs. pitch-yaw product (ryrz) components of the rotation vector over a range of simulated head orientations. Plots of roll (rx) vs. pitch-yaw product (ryrz) components of the rotation vector were also generated for a mechanical model of a Fick gimbal. The plots were generated for extension (B) and flexion (C) separately.

 
k-GIMBAL MODEL. Although the gimbal score is a measure of how close the system is to a 2 df rotating gimbal system, there is no clear relationship of how the gimbal score reflects the actual rotatory kinematic system. To explore how close the head movement control is to a 2 df system, we generalized the notion of gimbal score and utilized a metric based on a coefficient, k (k-gimbal model), which reflects the fraction of the angle by which the second axis rotates with respect to the angle of rotation about the first axis. We then determined how close the data were to a k-gimbal model prediction.

For the k-gimbal model, there are two choices for the first rotation axis, either yaw or pitch as described in the preceding text. Based on the anatomy of the joints between C1 and C2, and C1 and the occiput, we chose the representation, in which the first angle of rotation is about the pitch axis given by {theta}. The yaw axis, about which the second rotation takes place, is oriented by a fraction, k, of the pitch angle of rotation. Analytically, the first rotation about the spatial y axis can be given as

Formula 3(3)
The magnitude of the second rotation is {varphi}, and it occurs about an axis that has been rotated by a fraction k of the first rotation angle, {theta}. The unit vector along this axis, âz' is given by

Formula 4(4)
where âz' is a unit vector along the z axis and Ry (k{theta}) is given by

Formula 5(5)
The second rotation about this axis, Rz', can therefore be described by a similarity transformation

Formula 6(6)
where Rz({varphi}) is the rotation about the spatial z axis by an angle {varphi}. Composing the transformation in (Eqs. 5 and 6), the total transformation, R ({theta}, {varphi}) can be given as

Formula 7(7)
It should be noted that if k = 0, then R({theta}, {varphi}) represents an orientation according to the first two rotations of a Fick gimbal. If k = 1, then R({theta}, {varphi}) represents a Helmholtz gimbal where the pitch about a spatial horizontal axis is the first rotation. For k = 0.5, R({theta}, {varphi}) represents a gimbal that satisfies Listing's law. Thus the k-gimbal model is capable of representing a wide range of models.

Equation 7 also gives the following relationship between the roll component of the rotation vector, rx and the pitch and yaw components of the rotation vector for each head orientation as

Formula 8(8)

Formula 9(9)
Incorporating a possible offset in roll (rx0) and combining (Eqs. 8 and 9), the relationship between the rotation axis component, rx as a function of rz and ry is

Formula 10(10)
Equation 10 can then be used to fit data on head orientation to determine the k value and the organization of the k-gimbal system.


    EXPERIMENTAL RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Fick angles, {theta}, {phi}, {psi} for data obtained from all subjects for head orientations at up-left (Fig. 4A, {phi} = 27.5 ± 4.0°; {theta} = –31.8 ± 5.7°; {psi}data = –3.9 ± 5.7°), up-right (Fig. 4B, {phi} = –29.5.5 ± 3.8°; {theta} = –32.3 ± 4.5°; {psi}data = 3.1 ± 5.1°), center (Fig. 4C, {phi} = –1.1 ± 2.9°; {theta} = –4.0 ± 4.9°; {psi}data = –1.5 ± 2.2°), down-left (Fig. 4D, {phi} = 28.2 ± 3.8°; {theta} = 25.9 ± 4.5°; {psi}data = 1.1 ± 5.0°), and down-right (Fig. 4E, {phi} = –31.3 ± 3.0°; {theta} = 26.8 ± 4.2°; {psi}data = –3.6 ± 4.6°) had substantial roll components at the tertiary positions. This did not match the predictions of the truncated Fick-gimbal model at any head orientation (Fig. 4, AE, {psi}Fick = 0). The head orientations also did not obey Listing's law (Fig. 4, AE, {psi}List = –8.0°; {psi}List = 8.7°; {psi}List = 0.0°; {psi}List = 6.6°; {psi}List = –7.6°).


Figure 4
View larger version (14K):
[in this window]
[in a new window]

 
FIG. 4. Comparison of amount of roll observed in the data ({psi}data) with what would be obtained if the head followed the Fick gimbal ({psi}Fick) or Listing's Law ({psi}List) when the head points up-left (A), up-right (B), center (C), down-left (D), and down-right (E). The angles {phi} and {theta} at each of the head orientations are the Fick yaw and pitch angles and their SDs over all trials and subjects. Note that the data ({psi}data) do not fit either the Fick or the Listing's law models for all head positions.

 
We next considered whether Eq. 10, i.e., the k-gimbal model, was capable of fitting these experimental data. We found k by obtaining the best fit to the data from the orientations of the head for each subject using the simplex method for nonlinear optimization (Matlab 7.0, Mathworks). Orientations corresponding to flexion and extension were fitted separately, and different k values were obtained. As a basis of comparison, we also estimated the gimbal score, s, by fitting a line through the x component of the rotation vector versus the product of y and z components as described above. Using Eqs. 8 and 9, the relationship between gimbal score and parameter k was approximately linear, given by

Formula 11(11)
We also estimated the value of k by first finding the gimbal score and then using Eq. 11. This approach gave the same results but had significantly better computational efficiency.

When data for different yaw ({phi}) components of the head rotation were considered, the amount of roll was different for complementary head orientations during flexion and extension and so the two were analyzed separately. Figure 5 shows a plot of rx versus the product of ry*rz, which can be compared with the idealized plots of Fig. 3. The different colors correspond to data taken in different yaw orientations of the head. The circles represent the actual data and the lines are the linear regressions associated with the data for each yaw head position. For head extension (Fig. 5, A–C), the averaged data over three trials for three typical subjects had a declining roll component, rx, as a function of the pitch-yaw product, ryrz, for head yaw ranging from –30 to 30°. The linear regression lines representing the gimbal scores for each yaw head position had similar negative slopes for extension (Fig. 5, AC) and flexion (Fig. 5, DF). The coefficient k for three of the subjects with five or more repeated trials were (mk-flexion: 0.26 ± 0.09; mk-extension: 0.17 ± 0.08; tr-flexion: 0.16 ± 0.08; tr-extension: 0.21 ± 0.06; yo-flexion: 0.27 ± 0.04; yo-extension: 0.17 ± 0.08). The coefficient k for extension over all subjects was 0.17 ± 0.08. Similar behavior was observed during flexion, with an average k over all subjects having a value of 0.22 ± 0.10. The k values for flexion and extension were significantly different (matched t-test, P > 0.05), although the difference was small.


Figure 5
View larger version (18K):
[in this window]
[in a new window]

 
FIG. 5. Linear regressions of roll (rx) vs. pitch-yaw product (ryrz) for subjects CC (A and D), MK (B and E), and YO (C and F) computed in extension (A–C) and flexion (D–F) separately. Due to coordinate misalignment in pitch, the regression lines were obtained separately for each yaw orientation.

 
Thus the k-gimbal model suggests that voluntary head orientations are achieved not by Fick, Listing, or Helmholtz gimbal systems (Ceylan et al. 2000Go) but rather by a gimbal system that is tied to fractional changes in axis orientation as the head is moved into different orientations. The value of k = 0.2, shows that the head rotates as a gimbal where the first rotation axis is about pitch, and the second axis is 20% of the angle of pitch rotation. This behavior is probably not only related to the articulation of the head on the spinal column but to varying involvement of lower portions of the cervical spinal column during pitch and yaw movements in different head orientations (see Implementation of k-gimbal model).

Incremental head rotations in three dimensions

MATHEMATICAL BASIS. To determine the incremental axes of head rotation during oscillation, we developed a methodology for filtering movements of the head displacements using twists and exponential coordinates (Murray et al. 1994Go). Let D–n, D–n+1, ... , D0, ... , Dn1, Dn be a sequence of displacements given in a spatial coordinate frame, sampled from a small portion of the rigid body trajectory. We will assume that the sampled trajectory is "smooth" so that incremental displacements between, Dk and Dk+1, could be approximated by a single incremental displacement matrix e{Xi}, obtained from the twist {Xi} for k ranging from –n to n–1. This allowed us to define a matrix B over the window –n, ... , 0, ... , n, relative to the spatial frame of reference, such that each Dk could be obtained as this incremental displacement from B

Formula 12(12)
The incremental displacement, {Delta}k = DkDk–1, could then be approximated as

Formula 13(13)
and

Formula 14(14)
independent of B. To find the matrix, {Xi}, that gives the "best" approximation over the range of displacement matrices, D–n, D–n+1, ... , D0, ... , Dn–1, Dn, we utilized a procedure that is similar to minimizing the mean squared error for one-dimensional data. The value of the displacement matrix at some sample, Dk, can be multiplied by the inverse of the predicted displacement, (ek{Xi}B)–1. If the data were predicted perfectly by the approximation, then the matrix multiplication should be the identity transformation, I. Therefore an appropriate error at each sample displacement is the norm of the difference between the matrix product Dk(ek{Xi}B)–1 and the identity matrix I. Thus the total error, E, can be defined as

Formula 15(15)
where we choose || . ||, as the quadratic norm for matrices (Blum 1972Go; Press et al. 1992Go). The matrices {Xi} and B that minimize E can then be determined.

A straightforward minimization procedure for finding B and {Xi} would involve utilization of nonlinear optimization techniques (Press et al. 1992Go). The minimization procedure was simplified, however, by minimizing the residual norm (Blum 1972Go) for the incremental displacement, {Delta}k, which has the effect of eliminating the matrix B from the minimization procedure and is given as follows

Formula 16(16)
Because each {Delta}k is a displacement in its own right, it is generated by some twist, which is a function of {Delta}k, denoted by {Xi}({Delta}k). Thus instead of minimizing Eq. 16, we can minimize the function T, given by

Formula 17(17)
The optimal estimate of {Xi} in the least squares sense can be given as follows (Kunin 2004Go)

Formula 18(18)
The original sum of the differences over the window given by

Formula 19(19)
can serve as a measure of the goodness of fit. The smaller the value of E, the better the fit. Obtained in this way, the optimal twist estimate, {Xi}, yields both the angular velocity estimate {omega}, and pivot point estimate, q0, which can be given as

Formula 20(20)
where v is the differential translation part of the twist. This point in three-dimensional space q0, is a point on the rotational axis corresponding to the incremental displacement, e{Xi}.

VOLUNTARY HEAD OSCILLATION. For yaw head movements, the angular velocity, {omega}, in head coordinates, oscillated approximately sinusoidally and was maintained close to the yaw axis, with a small roll and still smaller pitch components (Fig. 6A). The point of intersection of the axes with the yaw (x-y) plane was close to the origin of the head coordinate frame when the head was directed at the center target (Fig. 6B). The zoomed inset shows the compactness of the data values. There were small variations of the point of intersection of the axes as the head was moved from center (Fig. 6B, x = –5.4 mm, y = –0.1 mm) to up-right (Fig. 6C, x = –8.1 mm, y = 4.1 mm), up-left (Fig. 6D, x = –11.7 mm, y = –0.6 mm), and down-right (Fig. 6E, x = –18.5 mm, y = 0.7 mm). The concentration of points close to a single axis is shown by the dense color in each of these four positions. A principal component analysis of the distribution of the intersection points indicated that the ellipse of the variance had a maximum singular value of 21.62 mm (Fig. 6C, up-right) and an eccentricity of 2 (variance of data points along x being twice that along y) with the major axis close to the naso-occipital axis. Over all subjects and trials, however, the location of the intersection point was not correlated with head orientation in any direction in 12 of the 15 experiments. This indicates, there was no consistent correlation between intersection point and head orientation, and that the pivot point of the rotation was within 40 mm of the interaural axis in the mid-sagittal plane.


Figure 6
View larger version (21K):
[in this window]
[in a new window]

 
FIG. 6. A: roll, pitch, and yaw components of angular velocity during voluntary head oscillation in yaw. B: distribution of instantaneous rotation axes in the yaw plane while looking at the center target in the representative subject. The line between the green diamonds and the line perpendicular to it from the red diamond are the y and x axes for the head coordinate frame, respectively. The ellipse represents the variance of the point distribution defined by its principal components (red and blue axes). The expanded view shows the tightness of the distribution of the rotation axes. C–E: distributions of the instantaneous rotation axes in the yaw (x-y) plane during voluntary head oscillations in yaw, while looking up-right (C), up-left (D), and down-right (E). The axes were similarly distributed close to the origin. Data were obtained from a representative subject (YO).

 
The average vector along the axis of head oscillation over all head orientations and all subjects was (0.27, 0.0, 0.96) in head coordinates, which was within 16° of the z axis of the head (z = 0.96) (Fig. 7A). There was a substantial x component (x = 0.27), indicating that the z-axis oscillation was tilted forward relative to the stereotaxic vertical. When the head was upright, the yaw axis of the head was along z1 and the velocity axis was along {omega}1 (Fig. 7A). When the head was pitched forward, the yaw axis of the head was along z2, and the velocity axis was along {omega}2. The general finding was that the angle between the velocity vector along the head oscillation axis, {omega}, and the z axis of the head decreased with head pitch forward and increased with head pitch back (Fig. 7, BD). When the angle at 0° of head tilt was equalized across all subjects, there was a linear relationship between head pitch and change in oscillation axis with small SD ({sigma} = 3.7°; Fig. 7E). Of note, when the head was tilted forward or back, the axis of yaw oscillation was tilted less than the pitch angle of the head. This difference was small, and did not conform to predications of the k-gimbal model, which governed static head orientations. Thus there must be transient effects during dynamic head oscillations that tend to bring the axis of rotation closer to the head yaw axis than would be predicted from the k-gimbal model.


Figure 7
View larger version (18K):
[in this window]
[in a new window]

 
FIG. 7. A: pictorial representation of the angle between the incremental axis of rotation ({omega}) and the yaw axis of the head (z). When the head is pitched forward by an angle {alpha}1 relative to the upright (Red z re zS), the angle is {phi}1. When the head is pitched back by an angle {alpha}2 relative to the upright (Black z re zS), the angle is f2. B–E: for each of 3 subjects and all subjects combined, there was a linear decline in angle between the incremental axis of rotation and the head yaw. In all instances, the incremental axis of rotation was close to the head yaw axis.

 
Pitch oscillation of the head had almost no roll or yaw components (Fig. 8A). The points of intersection of the pitch axes with the pitch plane had an eccentric distribution that was more elongated than the distribution of the intersection of the yaw oscillation axis in the yaw plane (Fig. 6, BE). The average intersection points during pitch oscillation also deviated to a greater extent from the origin of the head coordinate frame (Fig. 8B, x = –21.6, z = –15.8) to up-right (Fig. 8C, x = –10.5, z = 1.0), bottom-right (Fig. 8D, x = –32.6, –27.1), and up-left (Fig. 8E, x = –36.4, z = –30.6) than during yaw oscillation. The average point of intersection of the instantaneous pitch axis with the pitch plane of the head was {approx}22 mm behind the interaural axis and 16 mm below the origin of the head coordinate frame. The average x coordinates of the point of intersection of the pitch oscillation axis with the pitch plane was negatively correlated with the angle of pitch of the head in individual subjects (Fig. 9, AC). The data for all subjects was normalized by setting pivot x (Fig. 9, AC) for each subject's regression line to zero. When this was done, the x coordinate of the pivot point declined linearly with small SD ({sigma} = 7.6 mm; Fig. 9D). That is, as the head pitched down, the intersection was further back in head coordinates along the naso-occipital axis (Fig. 9). The average vector along the axis of head oscillation over all head orientations and all subjects was (0, 1.0, 0) in head coordinates, supporting the conclusion that the pitch axis in head coordinates was unaffected by any head position or dynamics.


Figure 8
View larger version (17K):
[in this window]
[in a new window]

 
FIG. 8. A: roll, pitch and yaw components of angular head velocity during voluntary head oscillation in pitch. B–E: distribution of instantaneous rotation axes in the pitch (x-z) plane while looking at the center (B), up-right (C), up-left (D), and down-right (E). Data were obtained from a representative subject (MK).

 

Figure 9
View larger version (15K):
[in this window]
[in a new window]

 
FIG. 9. A–C: change in x component of the mean value of the pivot point of head rotation as a function of pitch angle of the head for 3 subjects. D: cumulative change in the X component of the mean of the pivot points of head rotation as a function of pitch angle over all subjects. The X component of the pivot moves toward the back of the head as the head pitches forward.

 
HEAD OSCILLATION DURING LOCOMOTION. Treadmill walking at 1.5 m/s, while wearing the pinhole goggles and fixating a target at 25 cm, generated roll and yaw oscillations of the head at the stride frequency and pitch oscillations related to the step frequency (Hirasaki et al. 1999Go; Moore et al. 1999Go; Raphan et al. 2001Go). The head also translated vertically and horizontally (Hirasaki et al. 1999Go). There was greater variation in the oscillation amplitude and frequency, probably due to the presence of the pinhole goggles. The latter were used, however, to equalize the locomotor conditions to those during head oscillation. Walking without the goggles was also done to provide further comparison with natural locomotion.

The intersection point of the pitch axes while walking with the goggles and fixating at 25 cm (Fig. 10B, x = –20 mm, z = –147 mm) was considerably lower than during voluntary head oscillation with goggles (Fig. 9). The point of intersection was altered somewhat when walking slower at 0.9 m/s (Fig. 10C, x = –42 mm, z = –123 mm), but was still much lower along the z axis than during voluntary head oscillation. When walking at 1.5 m/s and viewing a target at 25 cm without goggles, the intersection point was consistent with walking with goggles and at slower walking velocities (Fig. 10D, x = –36 mm, z = –130 mm). Fixating a target without goggles at 1 m while walking generated an intersection point that was similar to the other conditions (Fig. 10E, x = –28 mm, z = –123 mm). The average intersection point over all subjects was (x = –30 mm, z = –149 mm). The average vector along the axis of head oscillation for all subjects was (0.116, 0.990, –0.073) in head coordinates, showing that the dominant direction of the axis about which the head rotates was aligned with the pitch axis of the head. There was no correlation among the walking velocity, viewing conditions, direction of the average axis of head rotation, and the location of the point about which the head rotated. Thus during locomotion the pitch axis of the head was aligned with the interaural axis and was lower on the spinal column than when pitching the head voluntarily.


Figure 10
View larger version (24K):
[in this window]
[in a new window]

 
FIG. 10. A: roll, pitch, and yaw components of angular velocity during locomotion on a linear treadmill at 1.5 m/s while wearing goggles and fixating a target at 25 cm. B–E: distribution of the intersection of the instantaneous rotation axes of the head with the pitch plane of the head during locomotion at different velocities and fixating targets at different distances. B: walking at 1.5 m/s while wearing goggles and fixating a space-fixed point at 25 cm from the head with. C: walking at 0.9 m/s while wearing goggles and fixating at 25 cm. D: walking at 1.5 m/s without wearing goggles and fixating a target at 25 cm. E: walking at 1.5 m/s without goggles and fixating a target at 1 m. In all instances, the distribution of the axes were below that for voluntary pitch head oscillation and more widely and evenly distributed. Data were obtained from a representative subject (TR).

 
Implementation of the k-gimbal model for head rotation

One realization for the k-gimbal model that is responsible for rotating the head would be to have a pitch rotation axis superimposed on a Fick-gimbal system as shown in Fig. 11A. This is consistent with how the head interfaces with the dens, occipital condyles, and the lower spinal column. Let {theta}1 be the angle of pitch rotation of the lower joint (Fig. 11B), {theta}2 the pitch angle of the head rotation about the upper joint, and {phi} the yaw angle of the head rotation. The equations governing the rotations can be given as follows

Formula 20

Formula 21(21)

Formula 22(22)
where Eq. 22 is equivalent to Eq. 7.


Figure 11
View larger version (12K):
[in this window]
[in a new window]

 
FIG. 11. A: model of the k-gimbal implementation. The upper pitch and yaw rotation implements a Fick gimbal. The addition of the pitch rotation below the Fick gimbal implements the k-gimbal system. B: angular orientations, {theta}2, {phi}, and {theta}1 associated with the k-gimbal orientation. See text for full description.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This study demonstrates that the rotation axes for positioning the head from a primary position follow a simple "k-gimbal model." The k-gimbal can be implemented by a sequence of two pitch joints with an interposed yaw joint (Fig. 11). The upper axes for this structure implement a two-dimensional Fick gimbal the movements of which are in pitch and yaw. The k-gimbal realization arises due to the linkage of the upper pitch and yaw joints to the chain of cervical vertebrae that implement the lower pitch axis. We also developed a mathematical methodology to determine incremental head rotation axes in three dimensions using displacement matrices (Kunin 2004Go; Kunin et al. 2006aGo,bGo). Using this methodology, which filters the translations and rotation in three dimensions, it was possible to obtain instantaneous axes of rotation as a function of time. From this, we were able to show that voluntary head oscillations about static head positions in pitch and yaw, and oscillations during locomotion had velocities that were not aligned with the axes of the k-gimbal. Rather, consistent with findings of others (Medendorp et al. 1998Go), the pitch and yaw incremental rotation axes were more closely aligned with the pitch and yaw axes of the head, regardless of head orientation. Moreover, the pivot points about which the head rotated were different for the different head rotations. Thus the model of head rotation axes developed from this study links the skeletal geometry to the constraints of Donders’ law during voluntary head positioning and also explains the consequences of head perturbations that occur during natural movement.

The finding that the points of intersection of the incremental pitch rotation axes during voluntary head oscillations were close to the origin of the head coordinate frame is consistent with previous suggestions that head movements are essentially aligned with the axes determined by vertebra and joints from C2 and above (Dutia 1991Go; Richmond and Vidal 1988Go; Vidal et al. 1988Go). The finding that incremental rotations of the head during locomotion were significantly lower than C2 supports the idea that the lower vertebrae play a more significant role in moving the head during locomotion than during voluntary head rotations and may be a result of the fore-aft movement of the body during locomotion. Engaging many levels of vertebral motion during locomotion would increase the radius of curvature of head rotation in space and could give the head movement controller more flexibility in maintaining heading (forward gaze) and aligning the head yaw axis with the spatial vertical. Maintenance of invariant yaw and roll in space would provide a stable platform for the visual system to process information during locomotion.

The different coordinate schema for global and voluntary incremental rotation axes may be related to the functions of the different types of head-positioning maneuvers. During global positioning, the head must be oriented relative to space so that the eyes can be positioned for generating gaze without requiring ocular compensatory mechanisms. Incremental rotations, on the other hand, are more apt to occur during perturbations of the head, which activate the vestibular apparatus to induce compensatory eye movements. The gain of the aVOR is maximal for yaw and pitch, being close to 1.0, and is only ~0.6 for roll (Henn 1988Go; Tweed et al. 1994Go). By maintaining the axes of voluntary incremental head rotations close to the pitch and yaw axes of the head, optimal ocular compensation could be maintained during incremental head movement, regardless of where one pointed the head during the movement without necessitating active roll compensation. This implies that incremental rotations that occur during voluntary head oscillations and during locomotion are governed by other factors, such as dynamic muscle action, that alter the axes of head rotation in relation to the k-gimbal axes, which are more spatially oriented.

Although it has been widely postulated that rotations of the head are in accordance with the skeletal geometry of the vertebra (Dutia 1989Go, 1991Go; Richmond and Vidal 1988Go; Vidal et al. 1986Go, 1988Go), more recent investigations have suggested a neuro-muscular basis for head-positioning geometry (Ceylan et al. 2000Go; Radau et al. 1994Go). Based on differences of orientation observed when subjects were asked to direct the head toward targets with and without monocular pin-hole goggles, it was suggested that there is a transformation from Fick-like to Listing-like behavior (Ceylan et al. 2000Go). Our data showed no difference in the distribution of head pointing with and without goggles, however (Fig. 5), and the k-gimbal value of k {approx} 0.22 that we obtained was equivalent to the gimbal score of s = 0.6 found by Ceylan et al. (2000)Go for comparable sized head movements. Thus the differences between wearing and not wearing goggles in the Ceylan et al. (2000)Go study may be attributed to the different-sized head movements for the two conditions as well as to an increase in the torsional thickness of the Donders’ surface associated with larger head movements (Ceylan et al. 2000Go).

An appealing aspect of the generalized k-gimbal is that it can be directly related to anatomical structures of the head and neck. In the model, global head positioning is essentially governed by a Fick gimbal, implemented anatomically by the joints of head positioning control at and above the level of C2 (Fig. 11). These include the occipital condyles that link C1 to the skull, producing pitch. Yaw rotation of the head is produced by rotation of C1 about the dens, the protruberent process of C2 (Ishii et al. 2004aGo,bGo). The lower vertebrae implement an additional pitch axis by flexing forward and extending back to generalize the Fick gimbal to the k-gimbal. There were also significant differences in extension and flexion in defining the k-gimbal axes about which the head rotated because the gimbal score for flexion was consistently smaller than for extension (Fig. 5). Thus our data are in agreement with the idea that the lower vertebrae contribute more to head orientation when the head is flexed than when it is extended (de Waele et al. 1988Go; Graf et al. 1992Go; Vidal et al. 1988Go). Taken together, these results generate the holonomic constraints (Eq. 22) that implement Donders’ law. This supports the idea that head fixation is governed by a k-gimbal that is directly related to the anatomy of the cervical spinal column and is not dependent on central commands, which modulate the velocity of "Donders’ operators" (Ceylan et al. 2000Go).

Because the head and eyes generally move together in generating gaze, why then is head orientation governed by a twisted Donders’ surface, whereas eye movements are governed by Listing's law, a planar realization of Donders’ law? The answer may lie in the clear differences in the pivot points about which the eye and head rotate as well as in how the central commands are organized in relation to the muscles that rotate the head and eyes (de Waele et al. 1988Go; Graf et al. 1992Go). There is considerable evidence that saccades are driven by two-dimensional control signals (Van Opstal et al. 1991Go) that activate the muscles and muscle pulleys that rotate the eyes about a central pivot in three dimensions (Raphan 1997Go, 1998Go; Schnabolk and Raphan 1994Go; Thurtell et al. 2000Go). Thus Listing's law is simply a representation of a two-dimensional control of a three-dimensional rotating plant (Raphan 1997Go, 1998Go). In contrast, head movements are constrained by a joint structure, represented by the occipital condyles and the dens, and by flexion and extension of the cervical spine that maintain the k-gimbal architecture. There are also different axes and pivots, which govern global head orientation than voluntary head oscillation and perturbations during locomotion. Thus if central commands to the head muscles are aligned in accordance with the spinal column structure that has many different pivot points and axes of rotation, head position would be naturally governed by a Donders’ surface not by a Listing's plane. The natural spinal structures may also explain why Donders’ law for head movement is obeyed even when muscle stiffness caused by spasmodic torticollis compromises head orientation (Medendorp et al. 1999Go) and eye-head coordination (Maurer et al. 2001Go).

There is a considerable literature on the physiology governing head movements (Freedman and Sparks 1997Go; Freedman et al. 1996Go; Harris 1980Go; Isa and Sasaki 2002Go; Roucoux and Guitton 1980Go; Wilson and Melvill-Jones 1979Go). However, the neural control mechanisms that govern head movements, especially in three dimensions, are not well known, and only recently have been studied in the primate (Crawford et al. 2003Go; Isa and Sasaki 2002Go; Klier et al. 2001Go, 2002aGo,bGo). The midbrain interstitial nucleus of Cajal (INC) functions as a neural integrator for holding head orientation (Klier et al. 2001Go), just as the neural velocity-position integrator encodes eye orientation during eye movements (Robinson 1971Go; 1968Go). Driving these integrators, the superior colliculus encodes gaze commands in retinal coordinates (Klier et al. 2001Go). If the superior colliculus encodes gaze shifts as a two-dimensional rotational command, then the vectorial command must be separated downstream in the brain stem to hold the head and eyes in a particular orientation. The knowledge that the head movements are coded according to the k-gimbal model, as determined in this study, would therefore be important for defining the neural signals that determine head trajectory and for understanding how the head velocity-position integrator activates the neck muscles in maintaining the final head orientation.


    GRANTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This work was supported by National Institute of Deafness and Other Communications Disorders Grants DC-05222 and DC-05204.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank L. Tundis of the Brooklyn College Machine Shop for fabricating the calibration device, pinhole goggles, and the target array that were utilized in this study.


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: T. Raphan, Dept. of Computer and Information Science, Brooklyn College of the City University of New York, 2900 Bedford Ave., Brooklyn, New York 11210 (E-mail: raphan{at}nsi.brooklyn.cuny.edu)


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 EXPERIMENTAL RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Blum EK. Numerical Analysis. Reading, MA: Addison-Wesley, 1972.

Ceylan M, Henriques DY, Tweed DB, Crawford JD. Task-dependent constraints in motor control: pinhole goggles make the head move like an eye. J Neurosci 20: 2719–2730, 2000.[Abstract/Free Full Text]

Crawford JD. Geometric Transformations in the Visual-Motor Interface for Saccades. In: Three-Dimensional Kinematics of Eye, Head and Limb Movements, edited by Fetter M, Haslwanter Th, Misslisch H, and Tweed D. Australia: Harwood Academic Publishers, 1997, p. 85–98.

Crawford JD, Ceylan MZ, Klier EM, Guitton D. Three-dimensional eye-head coordination during gaze saccades in the primate. J Neurophysiol 81: 1760–1782, 1999.[Abstract/Free Full Text]

Crawford JD, Guitton D. Visual-motor transformations required for accurate and kinematically correct saccades. J Neurophysiol 78: 1447–1467, 1997.[Abstract/Free Full Text]

Crawford JD, Martinez-Trujillo JC, Klier EM. Neural control of three-dimensional eye and head movements. Curr Opin Neurobiol 13: 655–662, 2003.[CrossRef][Web of Science][Medline]

Crawford JD, Vilis T. How does the brain deal with the problem of rotational movement. J Mot Behav 27: 89–99, 1995.[Web of Science]

de Waele C, Graf W, Berthoz A, Vidal PP. Vestibular control of skeletal geometry. In: Posture and Gait: Development, Adaptation, and Modulation, edited by Amblard B, Berthoz A, Clarac F. New York: Wiley, 1988, p. 423–431.

DeSouza JFX, Nicolle DA, Vilis T. Task-dependent changes in the shape and thickness of Listing's plane. Vision Res 37: 2271–2282, 1997.[CrossRef][Web of Science][Medline]

Dutia MB. Mechanisms of head stabilization. News Physiol Sci 4: 101–103, 1989.[Abstract/Free Full Text]

Dutia MB. The muscles and joints of the neck: their specialisation and role in head movement. Prog Neurobiol 37: 165–178, 1991.[CrossRef][Web of Science][Medline]

Freedman EG, Sparks DL. Activity of cells in the deeper layers of the superior colliculus of the rhesus monkey: evidence for a gaze displacement command. J Neurophysiol 78: 1669–1690, 1997.[Abstract/Free Full Text]

Freedman EG, Stanford TR, Sparks DL. Combined eye-head gaze shifts produced by electrical stimulation of the superior colliculus. J Neurophysiol 76: 927–952, 1996.[Abstract/Free Full Text]

Glenn B, Vilis T. Violations of Listing's law after large eye and head gaze shifts. J Neurophysiol 68: 309–318, 1992.[Abstract/Free Full Text]

Goldstein H. Classical Mechanics. Reading, MA: Addison-Wesley, 1980.

Graf W, de Waele C, Vidal PP. Skeletal geometry in verterbrates and its relation to the vestibular end organs. In: The Head-Neck Sensory Motor System, edited by Berthoz A, Graf W, Vidal PP. New York: Wiley, 1992, p. 129–134.

Gray H. Anatomy of the Human Body. Philadelphia: Lea and Febiger, 1918.

Harris LR. The superior colliculus and movements of the head and eye in the cats. J Physiol 300: 367–391, 1980.[Abstract/Free Full Text]

Haslwanter T, Straumann D, Hepp K, Hess BJM, Henn V. Smooth pursuit eye movements obey Listing's law in the monkey. Exp Brain Res 87: 470–472, 1991.[Web of Science][Medline]

Haustein W. Consideration on Listing's law and the primary position by means of a matrix description of eye position control. Biol Cybern 60: 411–420, 1989.[Web of Science][Medline]

Helmholz Hv. Handbuch der Physiologichen Optik. Leipzig, Germany: Leopold Voss, 1867.

Henn V. Representation of three-dimensional space in the vestibular, oculomotor and visual systems: Introduction. In: Representation of Three-Dimensional Space in the Vestibular, Oculomotor, and Visual Systems, edited by Cohen B, Henn V. New York: New York Academy of Sciences, 1988, p. 1–9.

Hirasaki E, Moore ST, Raphan T, Cohen B. Effects of walking velocity on vertical head and body movements during locomotion. Exp Brain Res 127: 117–130, 1999.[CrossRef][Web of Science][Medline]

Hore J, Watts S, Vilis T. Constraints on arm position when pointing in three dimensions: Donder's law and the Fick gimbal strategy. J Neurophysiol 68: 374–383, 1992.[Abstract/Free Full Text]

Isa T, Sasaki S. Brain stem control of head movements during orienting; organization of the premotor circuits. Prog Neurobiol 66: 205–241, 2002.[CrossRef][Web of Science][Medline]

Ishii T, Mukai Y, Hosono N, Sakaura H, Fujii R, Nakajima Y, Tamura S, Sugamoto K, Yoshikawa H. Kinematics of the subaxial cervical spine in rotation in vivo three-dimensional analysis. Spine 29: 2826–2831, 2004a.[CrossRef][Web of Science][Medline]

Ishii T, Mukai Y, Hosono N, Sakaura H, Nakajima Y, Sato Y, Sugamoto K, Yoshikawa H. Kinematics of the upper cervical spine in rotation: in vivo three-dimensional analysis. Spine 29: E139–144, 2004b.[CrossRef][Web of Science][Medline]

Klier EM, Wang H, Constantin AG, Crawford JD. Midbrain control of three-dimensional head orientation. Science 295: 1314–1316, 2002a.[Abstract/Free Full Text]

Klier EM, Wang H, Crawford JD. Neural mechanisms of three-dimensional eye and head movements. Ann NY Acad Sci 956: 512–514, 2002b.[Web of Science][Medline]

Klier EM, Wang H, Crawford JD. The superior colliculus encodes gaze commands in retinal coordinates. Nat Neurosci 4: 627–632, 2001.[CrossRef][Web of Science][Medline]

Kunin M. Computational approach for modeling body and head movements during locomotion. Doctoral Dissertation, Dept. of Computer Science. Graduate Center of City University of New York, 2004, p. 174.

Kunin M, Osaki Y, Cohen B, Raphan T. Head rotation in three dimensions during voluntary gaze shifts, oscillations in tertiary positions, and during locomotion. Soc for Neurosc Abstr 349.318, 2006a.

Kunin M, Osaki Y, Cohen B, Raphan T. Relationship between global orientations and incremental rotations of the head in three dimensions. In: ARO, 29th Midwinter Meeting. Denver: 2006b, p. 69.

Maurer C, Mergner T, Luking CH, Becker W. Adaptive changes of saccadic eye-head coordination resulting from altered head posture in torticollis spasmodicus. Brain 124: 413–426, 2001.[Abstract/Free Full Text]

Medendorp WP, Melis BJM, Gielen CCAM, van Gisbergen JAM. Off-centric rotation axes in natural head movements: implications for vestibular reafference and kinematic redundancy. J Neurophysiol 79: 2025–2039, 1998.[Abstract/Free Full Text]

Medendorp WP, van Gisbergen JAM, Horstink MWIM, Gielen CCAM. Donders’ law in torticollis. J Neurophysiol 82: 2833–2838, 1999.[Abstract/Free Full Text]

Moore ST, Hirasaki E, Cohen B, Raphan T. Effect of viewing distance on the generation of vertical eye movements during locomotion. Exp Brain Res 129: 347–361, 1999.[CrossRef][Web of Science][Medline]

Murray RM, Li Z, Sastry SS. A Mathematical Introduction to Robotic Manipulation. Boca Raton: CRC, 1994.

Peterson BW, Choi H, Hain T, Keshner E, Peng GCY. Dynamic and kinematic strategies for head movement control. Ann NY Acad Sci 942: 381–393, 2001.[Web of Science][Medline]

Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes in C: The Art of Scientific Computing. New York: Cambrige University Press, 1992.

Radau P, Tweed D, Vilis T. Three-dimensional eye head and chest orientations after large gaze shifts and the underlying neural strategies. J Neurophysiol 72: 2840–2852, 1994.[Abstract/Free Full Text]

Raphan T. Modeling control of eye orientation in three dimensions. In: Three-Dimensional Kinematics of Eye, Head, and Limb Movements, edited by Fetter M, Haslwanter T, Misslisch H, Tweed D. Amsterdam: Harwood Academic Publishers, 1997, p. 359–374.

Raphan T. Modeling control of eye orientation in three dimensions. I. Role of muscle pulleys in determining saccadic trajectory. J Neurophysiol 79: 2653–2667, 1998.[Abstract/Free Full Text]

Raphan T, Imai T, Moore ST, Cohen B. Vestibular compensation and orientation during locomotion. Ann NY Acad Sci 942: 128–138, 2001.[Web of Science][Medline]

Richmond FJR, Vidal PP. The motor system: joints and muscles of the neck. In: Control of Head Movement, edited by Peterson BW, Richmond FJR. New York: Oxford, 1988, p. 1–21.

Robinson DA. The oculomotor control system. Proc IEEE 56: 1032–1049, 1968.[CrossRef]

Robinson DA. Models of oculomotor neural organization. In: The Control of Eye Movements, edited by Bach-y-Rita P, Collins CC. New York: Academic, 1971, p. 519–538.

Roucoux A, Guitton D. Stimulation of the superior colliculus in the alert cat. II. Eye and head movements evoked when the head is unrestrained. Exp Brain Res 39: 75–85, 1980.[Web of Science][Medline]

Schnabolk C, Raphan T. Modeling three dimensional velocity-to-position transformation in oculomotor control. J Neurophysiol 71: 623–638, 1994.[Abstract/Free Full Text]

Stark L, Zangemeister WH, Edwards J, Grinberg J, Jones A, Lehman S, Lubock P, Narayan V, Nystrom M. Head rotation trajectories compared with eye saccades by main sequence relationships. Invest Ophthalmol Vis Sci 19: 986–988, 1980.[Abstract/Free Full Text]

Stark L, Zangmeister WH, Hannaford B. Head movement models, optimal control, and clinical application. In: Control of Head Movement, edited by Peterson BW and Richmond F. New York: Oxford, 1988, p. 245–260.

Straumann D, Haslwanter T, Hepp-Reymond MC, Hepp K. Listing's law for eye, head and arm movements and their synergistic control. Exp Brain Res 86: 209–215, 1991.[Web of Science][Medline]

Thurtell MJ, Kunin M, Raphan T. Role of muscle pulleys in producing eye position-dependence in the angular vestibuloocular reflex: a model-based study. J Neurophysiol 84: 639–650, 2000.[Abstract/Free Full Text]

Tweed D, Sievering D, Misslisch H, Fetter M, Zee D, Koenig E. Rotational kinematics of the human vestibuloocular reflex. I. Gain matrices. J Neurophysiol 72: 2467–2479, 1994.[Abstract/Free Full Text]

Tweed D, Vilis T. Implications of rotational kinematics for the oculomotor system in three dimensions. J Neurophysiol 58: 832–849, 1987.[Abstract/Free Full Text]

Van Opstal AJ, Hepp K, Hess B, Straumann D, Henn V. Two-, rather than three-dimensional representation of saccades in monkey superior colliculus. Science 252: 1313–1315, 1991.[Abstract/Free Full Text]

Vidal PP, de Waele C, Graf W, Berthoz A. Skeletal geometry underlying head movements. Ann NY Acad Sci 545: 228–238, 1988.[CrossRef][Medline]

Vidal PP, Graf W, Berthoz A. The orientation of the cervical vertebral column in unrestrained, awake animals. I. Resting posture. Exp Brain Res 61: 549–559, 1986.[Web of Science][Medline]

Wilson VJ, Melvill-Jones G. Mammalian Vestibular Physiology. New York: Plenum, 1979.




This article has been cited by other articles:


Home page
BMJHome page
D. Patton and A. McIntosh
Head and neck injury risks in heavy metal: head bangers stuck between rock and a hard bass
BMJ, December 17, 2008; 337(dec17_2): a2825 - a2825.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
Y. Xiang, S. B. Yakushin, M. Kunin, T. Raphan, and B. Cohen
Head Stabilization by Vestibulocollic Reflexes During Quadrupedal Locomotion in Monkey
J Neurophysiol, August 1, 2008; 100(2): 763 - 780.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
98/5/3095    most recent
00764.2007v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kunin, M.
Right arrow Articles by Raphan, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kunin, M.
Right arrow Articles by Raphan, T.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2007 by the The American Physiological Society.