|
|
||||||||
Department of Neurology with Center for Sensorimotor Research, Klinikum Grosshadern, Ludwig-Maximilians University, Munich, Germany
Submitted 22 June 2006; accepted in final form 30 January 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
If, however, the rotation involves more than a single axise.g., first a rotation around the subject's yaw axis and subsequently a rotation around the pitch axisa simple integration of angular velocity no longer suffices because of the kinematics of three-dimensional (3D) rotations. Instead, to accurately compute self-orientation in space from angular velocity signals, the instantaneous orientation has to be taken into account because rotations do not commute. In other words, the final orientation depends on the order of the applied rotations.
It has been shown that the brain accounts for the noncommutativity of rotations when the vestibuloocular reflex stabilizes gaze during whole body rotations (Tweed et al. 1999
) or during the remapping of visual targets during eye or head movements (Crawford et al. 2004
; Medendorp et al. 2002
; Smith and Crawford 2001
). It is not known at present whether this also holds for the updating of self-orientation. A straightforward way to test this hypothesis would be to submit subjects to the same passive whole body rotations, but in different order, thus yielding different final orientations. This approach was used before for the vestibuloocular reflex (Tweed et al. 1999
).
Here we used an alternative approach that is equally well suited to test whether noncommutativity is taken into account by the updating process: if different rotations starting at the same initial orientation and yielding the same final orientation are applied, any commutative operation should result in differences in perceived final self-orientation. This approach can also be regarded as a test of whether estimation of self-orientation is path independent. The APPENDIX (Noncommutativity and path independence) shows that our approach can indeed be used as an equivalent test of whether the path-integration process1 was aware of the noncommutativity of rotations.
On earth, rotations around more than one axis necessarily involve changes with respect to gravity. The straightforward approach of using the same rotations in different order would thus lead to different final orientations with respect to gravity, making it possible to infer the criterion of inequality of final orientations exclusively from the final graviceptive cues. The approach used here avoids this possibility: using different sequences of rotation, subjects were rotated from the same initial orientation (upright; see Fig. 1, left) to various supine orientations, which were identical with respect to gravity (Fig. 1, right). Therefore gravity cues in the final position could not be used to distinguish orientations. However, final subject orientations differed with respect to the angle
within the earth-horizontal plane (e.g., different sequences of rotation from upright facing north to supine with the head pointing either west or east). Subjects were then instructed to indicate their remembered initial orientation in the earth-horizontal plane. We hypothesized that if the required path integration process takes into account the noncommutativity of rotations, the same final orientation should correctly be perceived as being equal despite different paths taken to get there, but different final orientations should be distinguished despite equal graviceptive cues.
|
| METHODS |
|---|
|
|
|---|
Ten healthy subjects (two females, eight males, ages 2541 yr) without a history of vestibular or ocular motor problems gave their informed consent to participate in the present experiment, which conformed to the standards set by the Declaration of Helsinki. Subjects sat upright in a two-axis rotating chair (SEGA), initially always facing the same space-fixed direction in the horizontal plane (Fig. 1, left, defined as
= 0°). During the experiment, subjects were rotated in complete darkness to six supine orientations (Fig. 1, right). The supine orientations differed by the final orientation
in the horizontal plane, i.e., by the angle between the initial facing direction and the final longitudinal body axis. The final angles
with respect to the initial orientation were ±80, ±90, and ±100°. The two-axis rotating chair consists of a gimbal system with the inner subject-fixed axis parallel to the subject's nasooccipital roll axis and the outer earth-fixed axis is earth-horizontal (Fig. 1). Subjects, secured by safety belts and a head restraint, sat within a cabin that was closed throughout the experiment to shield subjects from airflow and any residual external light. Additionally, the laboratory was darkened throughout the experiment. To mask external sounds, subjects wore headphones, which delivered white noise.
To test for path independence, three different rotations were performed to bring subjects from the initial orientation into the final orientations (Fig. 2) : 1) pitch-roll: the subject was first pitched back by 90° into the supine position and then rolled around the nasooccipital axis by the angle
(±80, ±90, and ±100°) to reach the final orientation; 2) roll-yaw: the subject was first rolled around the angle
(±80, ±90, and ±100°) and subsequently rotated around the earth-horizontal axis by +90° or 90° to reach the supine position (for
= ±90°, this is a pure yaw rotation); 3) combined: the rotation was performed by rotating both axes of the chair at approximately the same time; this yielded an oblique time-varying axis of rotation resulting from the different peak velocities of the two axes (inner axis peak velocity 100°/s, outer axis 60°/s) and final orientations equal to those of pitch-roll and roll-yaw. All rotations were done with trapezoidal velocity profiles (see Fig. 2). Rotations were performed in randomized order.
|
Before each rotation, subjects were asked to remember their orientation in the earth-horizontal plane (see Fig. 1, left), i.e., their heading or straight-ahead direction (e.g., looking north while being in the initial position). After each rotation (duration <5 s; see Fig. 2) the supine subjects had to indicate their self-orientation in space by rotating an indicator, a computer-controlled luminous line polarized by an arrowhead and displayed on a dark TV screen, around the nasooccipital axis by remote control. The TV screen was mounted inside the chair's cabin in front of the subject, i.e., the screen was always in the same position with respect to the subject. The edges of the screen were masked by a blackboard with a circular hole, and the luminance of the screen was dimmed so that the interior of the cabin was not visible. The subjects were instructed to rotate the line within 30 s so that it pointed along the remembered initial orientation, as if the indicator were a compass needle (i.e., the indicator should have the same orientation in space as the subjects' straight-ahead before the rotation, e.g., when the initial heading was north, the correct response would have been to rotate the indicator so that it points north). After 30 s, the chair was quickly rotated back to the initial position and the next trial commenced after the subjects were reminded to memorize their orientation.
Chair position and indicator setting were recorded (125 Hz) for further analysis. Statistical analysis was done on the final indicator orientation pooled for positive and negative final orientations using a repeated-measures ANOVA with two within-subjects factors: orientation |
| (80, 90, 100°) and path (pitch-roll, roll-yaw, combined) using Statistica 6.1 (Statsoft, Tulsa, OK).
Model simulations
To show the difference in expected responses for commutative versus noncommutative processing, we simulated both cases. As implementation of the noncommutative model, any mathematically correct transformation from angular velocity to angular orientation is possible. In the simulations shown, we used a rotation vector model. A rotation vector (underline indicates a vector) describes angular orientation using only three components (Haustein 1989
)
![]() | (1) |
is the angle of rotation around an axis n. The differential equation to compute the angular orientation r from angular velocity
is given by (Hepp 1994
![]() | (2) |
denotes the inner or dot product and x the vector cross-product. Equation 2 was used for the noncommutative simulations.
Because a commutative model is mathematically incorrect for arbitrary rotations, the model chosen was designed to fulfill two criteria: 1) The commutative model should be correct for rotations around a single axis and 2) the angular orientation resulting from the commutative model should be a rotation vector (Eq. 1) for easy comparison. This results in
![]() | (3) |
| RESULTS |
|---|
|
|
|---|
Figure 3 depicts the model simulations for the noncommutative model (filled symbols) and the commutative model (open symbols). As expected, the noncommutative model shows veridical responses, i.e., simulated perceived orientation coincides with physical orientation (dashed line). Note that, because of the path independence of rotations, final orientations for the three paradigms overlap. In contrast, the simulations of the commutative model, which simply integrates the velocity stimulus, show a clear dependency on paradigm: whereas simulated responses for roll-yaw and combined rotations are comparable, as expected from their similar velocity profiles (cf. Fig. 2, middle and bottom), predicted responses for pitch-roll differ from these by about 80°.
|
Figure 4 shows the final indicator orientation for all ten subjects, rotational paths (circles: pitch-roll; squares: roll-yaw; diamonds: combined), and final orientations. Seven subjects adjusted the indicator close to veridical responses (diagonal lines in Fig. 4), thus clearly showing a response that reflected their final orientation. The three remaining subjects differed in that they either responded as if they had been pitched back to the supine position only (S5, S7), i.e., they adjusted the line around 180°, or showed a mixture of responses (S8). S8 reported to be confused by the experiment and considered herself unable to solve the task. Another subject (S1) adjusted the indicator twice to the 180° position.
|
|
| DISCUSSION |
|---|
|
|
|---|
Three of ten subjects (two females, one male) were not able to correctly indicate their self-orientation. Two subjects responded as if they had been tilted back by only 90°, but did so independently of condition (that is, even for roll-yaw, which involved no pitch rotation). One subject was confused; her responses varied between correct indication, mirror indication (side error), and pitch-back response. Because none of these subjects suffered from vestibular deficits, we hypothesize that their inability arises from a deficiency in the central processing of vestibular self-velocity information, such as in the ability to update self-orientation by 3D path integration.
The orientation estimates of the remaining subjects showed errors of about 20° on average. This raises the question of whether path-independent orientation estimates are compatible with such large estimation errors. An underestimation of angular velocity in one of the components (such as in roll) would lead to large differences between conditions (Fig. 5C, open symbols), i.e., to path-dependent responses even though the processing is still noncommutative (see the APPENDIX, Modified noncommutative models). However, in the present case, a solution that accounts for both path independence and estimations errors can be found as shown in Fig. 5C (filled symbols) by assuming differences in processing roll versus pitch and yaw information (see the APPENDIX, Modified noncommutative models). Thus the experimentally determined path independence and estimation errors are compatible with a noncommutative model. The weighting of the derivative of roll orientation used here is similar to that proposed earlier to explain the eye-position dependency of the torsional vestibuloocular reflex (Tweed 1997
). For visuospatial updating, it was recently shown that the processing of roll and yaw rotations differs indeed with respect to processing of gravity as an additional cue (Klier et al. 2006
).
Because 3D rotations on earth necessarily involve changes in orientation with respect to gravity, there is an alternative to the exclusive use of angular velocity information. Subjects could have used gravitational cues that are sensed by the otoliths, truncal graviceptors (Mittelstaedt 1996
), or somatosensation to determine their orientation with respect to gravity. However, note that the gravity cues available at the end of the rotation are not sufficient to determine the orientation in the earth-horizontal plane: subjects were in the supine position in all occasions and thus their final orientation with respect to gravity was always the same. Thus using gravity cues for this task also requires an ongoing updating process, i.e., keeping track of the change of the direction of gravity with respect to the head. Such a solution, if implemented correctly, allows accurate estimates of self-orientation and also involves noncommutative processing of both graviceptive and angular velocity cues (see the APPENDIX, Projection model; predictions of this "projection model" are equivalent to the noncommutative model shown in Fig. 5B). Using such processing, updating of the orientation in the earth-horizontal plane using only the changing gravity cues would be possible for the roll-yaw condition, but not for pitch-roll: because the second rotation in the pitch-roll condition was a turn around the earth-vertical axis (coinciding with the subject's nasooccipital axis), it should have gone unnoticed. In contrast, roll-yaw consisted of two rotations around earth-horizontal axes, thus making it possible to infer the final orientation solely from keeping track of changing graviceptive cues.
The importance of gravitational cues for orienting movements toward remembered targets was emphasized by various recent studies (e.g., Klier et al. 2005
, 2006
; Prieur et al. 2005; Van Pelt et al. 2005
). For example, subjects could make accurate saccades to remembered target locations after roll rotations around an earth-horizontal, but not around an earth-vertical axis (Klier et al. 2005
). However, such a dissociation was not observed in the present experiments. Another study (Van Pelt et al. 2005
) showed that, after intervening whole body roll rotations, errors in saccades to previously seen visual targets were not related to the intervening rotation in roll, but to body orientation with respect to gravity. This finding evidently points toward the use of gravitational cues in updating visual space, suggesting that the same may be true for estimation of self-orientation.
We conclude that whatever sensory input is used, the majority of subjects seem to be able to correctly integrate self-velocity or graviceptive cues into a consistent estimate of self-orientation in the horizontal plane, i.e., an estimate that is independent of the path taken to get to that orientation. As we have shown, such a path-independent estimate is possible only if the computation of self-orientation is noncommutative.
| APPENDIX |
|---|
|
|
|---|
Rotations are noncommutative, except in the one-dimensional case. This means, if two subsequent rotations are performed in reverse order, the final orientation of the rotated object will not be the same (e.g., Tweed et al. 1999
).
Consequently, temporal integration of a time-varying angular velocity vector alone does not yield angular position, except in the one-dimensional case. This is because temporal integration can be expressed as the cumulative summation of the integrated variable over infinitely small time steps. Because summation is commutative, so is integration.
Thus for temporal integration to yield angular orientation, it is not angular velocity, but the derivative of angular orientation that needs to be integrated. This derivative is computed from angular velocity and orientation (e.g., Haslwanter 1995
).
In the present experiments, two sequences of different rotations (pitch-roll and roll-yaw) were performed consecutively, yielding the same final orientation. During a particular rotation (e.g., roll), the orientation of the rotation axis with respect to the subject did not change, and thus each rotation was one-dimensional. Thus for each rotation alone, one could argue that a commutative operation would suffice to estimate the final orientation with respect to the starting orientation before this particular rotation. To answer the question of whether noncommutativity is known to the brain, the concatenation of the two rotations is crucial.
The final orientation r (underline indicates a vector) of the chair following the rotations ri was
![]() |
denotes the concatenation of two rotations, e.g., the product of two rotation matrices or two quaternions (e.g., Haslwanter 1995We will show now that, if the brain mistakenly assumed that rotations are commutative, the final orientation estimates would not be equal. That is, if the estimates of orientation for the same rotations, but in different order, were equal, then the estimates for the two sequences of different rotations shown earlier would not be the same.
In the following, the sign "^" indicates such estimates from commutative computations. The estimate of orientation after the sequence r1 followed by r2 can thus be written as
![]() |
![]() |
to denote any possible commutative concatenation of rotations.
Because commutativity is assumed, the order of rotations does not matter, and thus
![]() |
1,2 and
2,3 can be computed from rotation r2 followed by a second rotation (r1 or r3, respectively). This, however, shows that the estimate
1,2 could be equal to
2,3 only if rotations r1 and r3 were equal or zero, or if the brain were mistakenly unable to distinguish between both. The first was not the case in the present experiment and the second can be safely excluded because it would mean that the brain cannot distinguish a 90° pitch-back rotation from a 90° yaw rotation. Thus assuming commutativity, the two sequences would yield different orientation estimates, as had to be shown. Note that the proof makes no assumptions about how the integration of angular velocity or concatenation of successive rotations is actually performed by the brain, i.e., the experiment shows that the brain must necessarily take into account the noncommutativity of rotations to arrive at the same final estimate of orientation. This is also the case if the computations would be performed in an earth-fixed rather than a subject-centered coordinate system. To show this, one can simply translate the subject-centered rotations into earth-fixed rotations (pitch-roll translates to pitch-yaw, and roll-yaw translates to roll-pitch) and adopt the same line of reasoning as described earlier.
Modified noncommutative models
For simulations shown in Fig. 5C, the model in Eq. 2 (METHODS) was modified. For the underestimation of angular velocity (Fig. 3C, open symbols), the roll component of angular velocity was multiplied by 0.75.
For the simulations showing estimation errors similar to actual data together with path independence (filled symbols), the model was modified as follows
![]() |
![]() |
Projection model
Here, a mathematically correct model is derived for computing self-orientation in the earth-horizontal plane from two sensory inputs: angular velocity and gravity. Assuming that gravity g is known from graviceptive sensory input or a central gravity estimator (e.g., Glasauer and Merfeld 1997
; Merfeld and Zupan 2002
), the z-axis of the space-fixed coordinate system z can be expressed in subject-centered coordinates by
![]() |
The space-fixed coordinate axis x (within the earth-horizontal plane) expressed in subject-centered coordinates can then be computed from angular velocity
as
![]() |
into two components parallel and perpendicular to the z-axis, the x-axis can be expressed as
![]() | (A1) |
onto the z-axis
![]() | (A2) |
denotes the inner or dot product. In contrast, the cross-product 
x x, which points along the z-axis, can be expressed solely by using the gravity-derived z-axis and its temporal derivative
![]() | (A3) |
=
x z. Inserting Eqs. A2 and A3 into Eq. A1 yields
![]() | (A4) |
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
1 Path integration usually refers to the temporal integration of self-velocity to yield self-displacement. However, in accordance with the original definition (Mittelstaedt and Mittelstaedt 1980
), we suggest that any process is termed path integration, if it yields self-displacement from idiothetic cues such as self-velocity (temporal integration) or instantaneous displacement (spatial integration). ![]()
Address for reprint requests and other correspondence: S. Glasauer, Center for Sensorimotor Research, Department of Neurology, Marchioninistr. 23, 81377 Munich, Germany (E-mail: sglasauer{at}nefo.med.uni-muenchen.de)
| REFERENCES |
|---|
|
|
|---|
Glasauer S, Mayr A, Brandt T. Estimation of heading in three dimensions : is there a 3D velocity to position integrator? Abstracts for the XIX Barany Soc. Meeting, Sydney, Australia, 1996. J Vestib Res 6/4S: 68, 1996.
Glasauer S, Merfeld DM. Modeling three-dimensional vestibular responses during complex motion stimulation. In: Three-Dimensional Kinematic Principles of Eye-, Head-, and Limb Movements in Health and Disease, edited by Fetter M, Haslwanter T, Misslisch H, Tweed D. Amsterdam: Harwood, 1997, p. 387398.
Guedry FE. Psychophysics of vestibular sensation. In: Handbook of Sensory Physiology: Vestibular System, edited by Kornhuber HH. New York: Springer-Verlag, 1974, vol. VI/2, p. 1154.
Haslwanter T. Mathematics of three-dimensional eye rotations. Vision Res 35: 17271739, 1995.[CrossRef][ISI][Medline]
Haustein W. Considerations on Listing's Law and the primary position by means of a matrix description of eye position control. Biol Cybern 60: 411420, 1989.[ISI][Medline]
Hepp K. Oculomotor control: Listing's law and all that. Curr Opin Neurobiol 4: 862868, 1994.[CrossRef][Medline]
Israel I, Sievering D, Koenig E. Self-rotation estimate about the vertical axis. Acta Otolaryngol 115: 38, 1995.[Medline]
Jürgens R, Boss T, Becker W. Estimation of self-turning in the dark: comparison between active and passive rotation. Exp Brain Res 128: 491504, 1999.[CrossRef][ISI][Medline]
Klatzky RL, Loomis JM, Beall AC, Chance SS, Golledge RG. Spatial updating of self-position and orientation during real, imagined, and virtual locomotion. Psychol Sci 9: 293298, 1998.[CrossRef][ISI]
Klier EM, Angelaki DE, Hess BJ. Roles of gravitational cues and efference copy signals in the rotational updating of memory saccades. J Neurophysiol 94: 468478, 2005.
Klier EM, Hess BJ, Angelaki DE. Differences in the accuracy of human visuospatial memory after yaw and roll rotations. J Neurophysiol 95: 26922697, 2006.
Mach E. Grundlinien der Lehre von den Bewegungsempfindungen. Leipzig, Germany: Engelmann, 1875.
Marlinsky VV. Vestibular and vestibulo-proprioceptive perception of motion in the horizontal plane in blindfolded manII. Estimations of rotations about the earth-vertical axis. Neuroscience 90: 395401, 1999.[CrossRef][ISI][Medline]
Medendorp WP, Smith MA, Tweed DB, Crawford JD. Rotational remapping in human spatial memory during eye and head motion. J Neurosci 22: RC196, 2002.
Merfeld DM, Zupan LH. Neural processing of gravitoinertial cues in humans. III. Modeling tilt and translation responses. J Neurophysiol 87: 819833, 2002.
Mergner T, Siebold C, Schweigart G, Becker W. Human perception of horizontal trunk and head rotation in space during vestibular and neck stimulation. Exp Brain Res 85: 389404, 1991.[ISI][Medline]
Mittelstaedt H. Somatic graviception. Biol Psychol 42: 5374, 1996.[CrossRef][ISI][Medline]
Mittelstaedt H, Mittelstaedt ML. Homing by path integration in a mammal. Naturwissenschaften 67: 566, 1980.[CrossRef][ISI]
Mittelstaedt H, Mittelstaedt ML. The effect of centrifugal force on the perception of rotation about a vertical axis. Naturwissenschaften 84: 366369, 1997.[CrossRef][ISI][Medline]
Mittelstaedt ML, Glasauer S. Idiothetic navigation in gerbils and humans. Zool Jb Physiol 95: 427435, 1991.
Smith MA, Crawford JD. Implications of ocular kinematics for the internal updating of visual space. J Neurophysiol 86: 21122117, 2001.
Tweed D. Velocity-to-position transformation in the VOR and the saccadic system. In: Three-Dimensional Kinematic Principles of Eye-, Head-, and Limb Movements in Health and Disease, edited by Fetter M, Haslwanter T, Misslisch H, Tweed D. Amsterdam: Harwood, 1997, p. 375386.
Tweed DB, Haslwanter TP, Happe V, Fetter M. Non-commutativity in the brain. Nature 399: 261263, 1999.[CrossRef][Medline]
Van Pelt S, Van Gisbergen JA, Medendorp WP. Visuospatial memory computations during whole-body rotations in roll. J Neurophysiol 94: 14321442, 2005.
This article has been cited by other articles:
![]() |
E. M. Klier, B. J. M. Hess, and D. E. Angelaki Human Visuospatial Updating After Passive Translations in Three-Dimensional Space J Neurophysiol, April 1, 2008; 99(4): 1799 - 1809. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Tchelidze and B. J. M. Hess Noncommutative Control in the Rotational Vestibuloocular Reflex J Neurophysiol, January 1, 2008; 99(1): 96 - 111. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. M. Klier, D. E. Angelaki, and B. J. M. Hess Human Visuospatial Updating After Noncommutative Rotations J Neurophysiol, July 1, 2007; 98(1): 537 - 541. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |