Charts on SO(3)

Last updated

In mathematics, the special orthogonal group in three dimensions, otherwise known as the rotation group SO(3), is a naturally occurring example of a manifold. The various charts on SO(3) set up rival coordinate systems: in this case there cannot be said to be a preferred set of parameters describing a rotation. There are three degrees of freedom, so that the dimension of SO(3) is three. In numerous applications one or other coordinate system is used, and the question arises how to convert from a given system to another.

Contents

The space of rotations

In geometry the rotation group is the group of all rotations about the origin of three-dimensional Euclidean space R3 under the operation of composition. [1] By definition, a rotation about the origin is a linear transformation that preserves length of vectors (it is an isometry) and preserves orientation (i.e. handedness) of space. A length-preserving transformation which reverses orientation is called an improper rotation. Every improper rotation of three-dimensional Euclidean space is a rotation followed by a reflection in a plane through the origin.

Composing two rotations results in another rotation; every rotation has a unique inverse rotation; and the identity map satisfies the definition of a rotation. Owing to the above properties, the set of all rotations is a group under composition. Moreover, the rotation group has a natural manifold structure for which the group operations are smooth; so it is in fact a Lie group. The rotation group is often denoted SO(3) for reasons explained below.

The space of rotations is isomorphic with the set of rotation operators and the set of orthonormal matrices with determinant +1. It is also closely related (double covered) with the set of quaternions with their internal product, as well as to the set of rotation vectors (though here the relation is harder to describe, see below for details), with a different internal composition operation given by the product of their equivalent matrices.

Rotation vectors notation arise from the Euler's rotation theorem which states that any rotation in three dimensions can be described by a rotation by some angle about some axis. Considering this, we can then specify the axis of one of these rotations by two angles, and we can use the radius of the vector to specify the angle of rotation. These vectors represent a ball in 3D with an unusual topology.

This 3D solid sphere is equivalent to the surface of a 4D disc, which is also a 3D variety. For doing this equivalence, we will have to define how will we represent a rotation with this 4D-embedded surface.

The hypersphere of rotations

Two rotations by different angles and different axes in the space of rotations. The length of the vector is related to the magnitude of the rotation. Space of rotations.png
Two rotations by different angles and different axes in the space of rotations. The length of the vector is related to the magnitude of the rotation.

Visualizing the hypersphere

It is interesting to consider the space as the three-dimensional sphere S3, the boundary of a disk in 4-dimensional Euclidean space. For doing this, we will have to define how we represent a rotation with this 4D-embedded surface.

The way in which the radius can be used to specify the angle of rotation is not straightforward. It can be related to circles of latitude in a sphere with a defined north pole and is explained as follows:

Beginning at the north pole of a sphere in three-dimensional space, we specify the point at the north pole to represent the identity rotation. In the case of the identity rotation, no axis of rotation is defined, and the angle of rotation (zero) is irrelevant. A rotation with its axis contained in the xy-plane and a very small rotation angle can be specified by a slice through the sphere parallel to the xy-plane and very near the north pole. The circle defined by this slice will be very small, corresponding to the small angle of the rotation. As the rotation angles become larger, the slice moves southward, and the circles become larger until the equator of the sphere is reached, which will correspond to a rotation angle of 180 degrees. Continuing southward, the radii of the circles now become smaller (corresponding to the absolute value of the angle of the rotation considered as a negative number). Finally, as the south pole is reached, the circles shrink once more to the identity rotation, which is also specified as the point at the south pole. Notice that a number of characteristics of such rotations and their representations can be seen by this visualization.

The space of rotations is continuous, each rotation has a neighborhood of rotations which are nearly the same, and this neighborhood becomes flat as the neighborhood shrinks.

Aliases

The hypersphere of rotations for the rotations that have a "horizontal" axis (in the xy-plane). Hypersphere of rotations.png
The hypersphere of rotations for the rotations that have a "horizontal" axis (in the xy-plane).

Also, each rotation is actually represented by two antipodal points on the sphere, which are at opposite ends of a line through the center of the sphere. This reflects the fact that each rotation can be represented as a rotation about some axis, or, equivalently, as a negative rotation about an axis pointing in the opposite direction (a so-called double cover). The "latitude" of a circle representing a particular rotation angle will be half of the angle represented by that rotation, since as the point is moved from the north to south pole, the latitude ranges from zero to 180 degrees, while the angle of rotation ranges from 0 to 360 degrees. (the "longitude" of a point then represents a particular axis of rotation.) Note however that this set of rotations is not closed under composition.

Two successive rotations with axes in the xy-plane will not necessarily give a rotation whose axis lies in the xy-plane, and thus cannot be represented as a point on the sphere. This will not be the case with a general rotation in 3-space, which do form a closed set under composition.

This visualization can be extended to a general rotation in 3-dimensional space. The identity rotation is a point, and a small angle of rotation about some axis can be represented as a point on a sphere with a small radius. As the angle of rotation grows, the sphere grows, until the angle of rotation reaches 180 degrees, at which point the sphere begins to shrink, becoming a point as the angle approaches 360 degrees (or zero degrees from the negative direction). This set of expanding and contracting spheres represents a hypersphere in four-dimensional space (a 3-sphere).

Just as in the simpler example above, each rotation represented as a point on the hypersphere is matched by its antipodal point on that hypersphere. The "latitude" on the hypersphere will be half of the corresponding angle of rotation, and the neighborhood of any point will become "flatter" (i.e. be represented by a 3D Euclidean space of points) as the neighborhood shrinks.

This behavior is matched by the set of unit quaternions: A general quaternion represents a point in a four-dimensional space, but constraining it to have unit magnitude yields a three-dimensional space equivalent to the surface of a hypersphere. The magnitude of the unit quaternion will be unity, corresponding to a hypersphere of unit radius.

The vector part of a unit quaternion represents the radius of the 2-sphere corresponding to the axis of rotation, and its magnitude is the sine of half the angle of rotation. Each rotation is represented by two unit quaternions of opposite sign, and, as in the space of rotations in three dimensions, the quaternion product of two unit quaternions will yield a unit quaternion. Also, the space of unit quaternions is "flat" in any infinitesimal neighborhood of a given unit quaternion.

Parametrizations

We can parameterize the space of rotations in several ways, but degenerations will always appear. For example, if we use three angles (Euler angles), such parameterization is degenerate at some points on the hypersphere, leading to the problem of gimbal lock. We can avoid this by using four Euclidean coordinates w,x,y,z, with w2 + x2 + y2 + z2 = 1. The point (w,x,y,z) represents a rotation around the axis directed by the vector (x,y,z) by an angle

This problem is similar to parameterize the bidimensional surface of a sphere with two coordinates, such as latitude and longitude. Latitude and longitude are ill-behaved (degenerate) at the north and south poles, though the poles are not intrinsically different from any other points on the sphere. At the poles (latitudes +90° and −90°), the longitude becomes meaningless. It can be shown that no two-parameter coordinate system can avoid such degeneracy.

The possible parametrizations candidates include:

Problems of the parametrizations

There are problems in using these as more than local charts, to do with their multiple-valued nature, and singularities. That is, one must be careful above all to work only with diffeomorphisms in the definition of chart. Problems of this sort are inevitable, since SO(3) is diffeomorphic to real projective space P3(R), which is a quotient of S3 by identifying antipodal points, and charts try to model a manifold using R3.

This explains why, for example, the Euler angles appear to give a variable in the 3-torus, and the unit quaternions in a 3-sphere. The uniqueness of the representation by Euler angles breaks down at some points (cf. gimbal lock), while the quaternion representation is always a double cover, with q and −q giving the same rotation.

If we use a skew-symmetric matrix, every 3 × 3 skew-symmetric matrix is determined by 3 parameters, and so at first glance, the parameter space is R3. Exponentiating such a matrix results in an orthogonal 3 × 3 matrix of determinant 1 – in other words, a rotation matrix, but this is a many-to-one map. Note that it is not a covering map – while it is a local homeomorphism near the origin, it is not a covering map at rotations by 180 degrees. It is possible to restrict these matrices to a ball around the origin in R3 so that rotations do not exceed 180 degrees, and this will be one-to-one, except for rotations by 180 degrees, which correspond to the boundary S2, and these identify antipodal points – this is the cut locus. The 3-ball with this identification of the boundary is P3(R). A similar situation holds for applying a Cayley transform to the skew-symmetric matrix.

Axis angle gives parameters in S2 × S1; if we replace the unit vector by the actual axis of rotation, so that n and −n give the same axis line, the set of axis becomes P2(R), the real projective plane. But since rotations around n and −n are parameterized by opposite values of θ, the result is an S1 bundle over P2(R), which turns out to be P3(R).

Fractional linear transformations use four complex parameters, a, b, c, and d, with the condition that adbc is non-zero. Since multiplying all four parameters by the same complex number does not change the parameter, we can insist that adbc=1. This suggests writing (a,b,c,d) as a 2 × 2 complex matrix of determinant 1, that is, as an element of the special linear group SL(2,C). But not all such matrices produce rotations: conformal maps on S2 are also included. To only get rotations we insist that d is the complex conjugate of a, and c is the negative of the complex conjugate of b. Then we have two complex numbers, a and b, subject to |a|2+|b|2=1. If we write a+bj, this is a quaternion of unit length.

Ultimately, since R3 is not P3(R), there will be a problem with each of these approaches. In some cases, we need to remember that certain parameter values result in the same rotation, and to remove this issue, boundaries must be set up, but then a path through this region in R3 must then suddenly jump to a different region when it crosses a boundary. Gimbal lock is a problem when the derivative of the map is not full rank, which occurs with Euler angles and Tait–Bryan angles, but not for the other choices. The quaternion representation has none of these problems (being a two-to-one mapping everywhere), but it has 4 parameters with a condition (unit length), which sometimes makes it harder to see the three degrees of freedom available.

Applications

One area in which these considerations, in some form, become inevitable, is the kinematics of a rigid body. One can take as definition the idea of a curve in the Euclidean group E(3) of three-dimensional Euclidean space, starting at the identity (initial position). The translation subgroup T of E(3) is a normal subgroup, with quotient SO(3) if we look at the subgroup E+(3) of direct isometries only (which is reasonable in kinematics). The translational part can be decoupled from the rotational part in standard Newtonian kinematics by considering the motion of the center of mass, and rotations of the rigid body about the center of mass. Therefore, any rigid body movement leads directly to SO(3), when we factor out the translational part.

These identifications illustrate that SO(3) is connected but not simply connected. As to the latter, in the ball with antipodal surface points identified, consider the path running from the "north pole" straight through the center down to the south pole. This is a closed loop, since the north pole and the south pole are identified. This loop cannot be shrunk to a point, since no matter how you deform the loop, the start and end point have to remain antipodal, or else the loop will "break open". In terms of rotations, this loop represents a continuous sequence of rotations about the z-axis starting and ending at the identity rotation (i.e. a series of rotations through an angle φ where φ runs from 0 to 2π).

Surprisingly, if you run through the path twice, i.e., from north pole down to south pole and back to the north pole so that φ runs from 0 to 4π, you get a closed loop which can be shrunk to a single point: first move the paths continuously to the ball's surface, still connecting north pole to south pole twice. The second half of the path can then be mirrored over to the antipodal side without changing the path at all. Now we have an ordinary closed loop on the surface of the ball, connecting the north pole to itself along a great circle. This circle can be shrunk to the north pole without problems. The Balinese plate trick and similar tricks demonstrate this practically.

The same argument can be performed in general, and it shows that the fundamental group of SO(3) is cyclic group of order 2. In physics applications, the non-triviality of the fundamental group allows for the existence of objects known as spinors, and is an important tool in the development of the spin–statistics theorem.

The universal cover of SO(3) is a Lie group called Spin(3). The group Spin(3) is isomorphic to the special unitary group SU(2); it is also diffeomorphic to the unit 3-sphere S3 and can be understood as the group of unit quaternions (i.e. those with absolute value 1). The connection between quaternions and rotations, commonly exploited in computer graphics, is explained in quaternions and spatial rotations. The map from S3 onto SO(3) that identifies antipodal points of S3 is a surjective homomorphism of Lie groups, with kernel {±1}. Topologically, this map is a two-to-one covering map.

See also

Related Research Articles

Sphere Geometrical object that is the surface of a ball

A sphere is a geometrical object that is a three-dimensional analogue to a two-dimensional circle. A sphere is the set of points that are all at the same distance r from a given point in three-dimensional space. That given point is the centre of the sphere, and r is the sphere's radius. The earliest known mentions of spheres appear in the work of the ancient Greek mathematicians.

3-sphere Mathematical object

In mathematics, a 3-sphere is a higher-dimensional analogue of a sphere. It may be embedded in 4-dimensional Euclidean space as the set of points equidistant from a fixed central point. Analogous to how the boundary of a ball in three dimensions is an ordinary sphere, the boundary of a ball in four dimensions is a 3-sphere. A 3-sphere is an example of a 3-manifold and an n-sphere.

Quaternion Noncommutative extension of the real numbers

In mathematics, the quaternion number system extends the complex numbers. Quaternions were first described by the Irish mathematician William Rowan Hamilton in 1843 and applied to mechanics in three-dimensional space. Hamilton defined a quaternion as the quotient of two directed lines in a three-dimensional space, or, equivalently, as the quotient of two vectors. Multiplication of quaternions is noncommutative.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition.

Unit quaternions, known as versors, provide a convenient mathematical notation for representing spatial orientations and rotations of elements in three dimensional space. Specifically, they encode information about an axis-angle rotation about an arbitrary axis. Rotation and orientation quaternions have applications in computer graphics, computer vision, robotics, navigation, molecular dynamics, flight dynamics, orbital mechanics of satellites, and crystallographic texture analysis.

Elliptic geometry Non-Euclidean geometry

Elliptic geometry is an example of a geometry in which Euclid's parallel postulate does not hold. Instead, as in spherical geometry, there are no parallel lines since any two lines must intersect. However, unlike in spherical geometry, two lines are usually assumed to intersect at a single point. Because of this, the elliptic geometry described in this article is sometimes referred to as single elliptic geometry whereas spherical geometry is sometimes referred to as double elliptic geometry.

Rotation (mathematics) Motion of a certain space that preserves at least one point

Rotation in mathematics is a concept originating in geometry. Any rotation is a motion of a certain space that preserves at least one point. It can describe, for example, the motion of a rigid body around a fixed point. Rotation can have sign : a clockwise rotation is a negative magnitude so a counterclockwise turn has a positive magnitude. A rotation is different from other types of motions: translations, which have no fixed points, and (hyperplane) reflections, each of them having an entire (n − 1)-dimensional flat of fixed points in a n-dimensional space.

Real projective plane Compact non-orientable two-dimensional manifold

In mathematics, the real projective plane is an example of a compact non-orientable two-dimensional manifold; in other words, a one-sided surface. It cannot be embedded in standard three-dimensional space without intersecting itself. It has basic applications to geometry, since the common construction of the real projective plane is as the space of lines in passing through the origin.

Euler angles Description of the orientation of a rigid body

The Euler angles are three angles introduced by Leonhard Euler to describe the orientation of a rigid body with respect to a fixed coordinate system.

Hopf fibration Fiber bundle of the 3-sphere over the 2-sphere, with 1-spheres as fibers

In the mathematical field of differential topology, the Hopf fibration describes a 3-sphere in terms of circles and an ordinary sphere. Discovered by Heinz Hopf in 1931, it is an influential early example of a fiber bundle. Technically, Hopf found a many-to-one continuous function from the 3-sphere onto the 2-sphere such that each distinct point of the 2-sphere is mapped from a distinct great circle of the 3-sphere. Thus the 3-sphere is composed of fibers, where each fiber is a circle — one for each point of the 2-sphere.

In linear algebra, a rotation matrix is a transformation matrix that is used to perform a rotation in Euclidean space. For example, using the convention below, the matrix

Eulers rotation theorem Movement with a fixed point is rotation

In geometry, Euler's rotation theorem states that, in three-dimensional space, any displacement of a rigid body such that a point on the rigid body remains fixed, is equivalent to a single rotation about some axis that runs through the fixed point. It also means that the composition of two rotations is also a rotation. Therefore the set of rotations has a group structure, known as a rotation group.

In mathematics, the group of rotations about a fixed point in four-dimensional Euclidean space is denoted SO(4). The name comes from the fact that it is the special orthogonal group of order 4.

Orientation (geometry) Notion of pointing in a direction

In geometry, the orientation, angular position, attitude, or direction of an object such as a line, plane or rigid body is part of the description of how it is placed in the space it occupies. More specifically, it refers to the imaginary rotation that is needed to move the object from a reference placement to its current placement. A rotation may not be enough to reach the current placement. It may be necessary to add an imaginary translation, called the object's location. The location and orientation together fully describe how the object is placed in space. The above-mentioned imaginary rotation and translation may be thought to occur in any order, as the orientation of an object does not change when it translates, and its location does not change when it rotates.

In mathematics, a versor is a quaternion of norm one. The word is derived from Latin versare = "to turn" with the suffix -or forming a noun from the verb. It was introduced by William Rowan Hamilton in the context of his quaternion theory.

Screw axis Geometric axis of rotation and translation

A screw axis is a line that is simultaneously the axis of rotation and the line along which translation of a body occurs. Chasles' theorem shows that each Euclidean displacement in three-dimensional space has a screw axis, and the displacement can be decomposed into a rotation about and a slide along this screw axis.

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

Axis–angle representation

In mathematics, the axis–angle representation of a rotation parameterizes a rotation in a three-dimensional Euclidean space by two quantities: a unit vector e indicating the direction of an axis of rotation, and an angle θ describing the magnitude of the rotation about the axis. Only two numbers, not three, are needed to define the direction of a unit vector e rooted at the origin because the magnitude of e is constrained. For example, the elevation and azimuth angles of e suffice to locate it in any particular Cartesian coordinate frame.

Six-dimensional space is any space that has six dimensions, six degrees of freedom, and that needs six pieces of data, or coordinates, to specify a location in this space. There are an infinite number of these, but those of most interest are simpler ones that model some aspect of the environment. Of particular interest is six-dimensional Euclidean space, in which 6-polytopes and the 5-sphere are constructed. Six-dimensional elliptical space and hyperbolic spaces are also studied, with constant positive and negative curvature.

In mathematics and mechanics, the Euler–Rodrigues formula describes the rotation of a vector in three dimensions. It is based on Rodrigues' rotation formula, but uses a different parametrization.

References

  1. Jacobson (2009), p. 34, Ex. 14.