Axis–angle representation

The angle axis vector

θ = θ e

is a unit vector

e

multiplied by an angle

θ

In mathematics, the axis–angle representation of a rotation parameterizes a rotation in a three-dimensional Euclidean space by three 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 the unit vector $e$ because its magnitude is constrained. The angle $θ$ scalar multiplied by the unit vector $e$ is the axis-angle vector

\boldsymbol{\theta} = \theta\mathbf{e} \,.

The vector itself does not perform rotations, but is used to construct transformations on vectors that correspond to rotations. The rotation occurs in the sense prescribed by the right-hand rule. The rotation axis is sometimes called the Euler axis.

It is one of many rotation formalisms in three dimensions. The axis–angle representation is predicated on Euler's rotation theorem, which dictates that any rotation or sequence of rotations of a rigid body in a three-dimensional space is equivalent to a pure rotation about a single fixed axis.

Rotation vector

The axis–angle representation is equivalent to the more concise rotation vector, also called the Euler vector. In this case, both the rotation axis and the angle are represented by a vector codirectional with the rotation axis whose length is the rotation angle $θ$ ,

\boldsymbol{\theta} = \theta \mathbf{e} ~.

It is used for the exponential and logarithm maps involving this representation.

Example

Say you are standing on the ground and you pick the direction of gravity to be the negative z direction. Then if you turn to your left, you will rotate π/2 radians (or 90°) about the z axis. Viewing the axis-angle representation as an ordered pair, this would be

( \mathrm{axis}, \mathrm{angle} ) = \left( \begin{bmatrix} e_x \\ e_y \\ e_z \end{bmatrix},\theta \right) = \left( \begin{bmatrix} 0 \\ 0 \\ 1 \end{bmatrix},\frac{\pi}{2}\right).

The above example can be represented as a rotation vector with a magnitude of $π ⁄ 2$ pointing in the $z$ direction,

\begin{bmatrix} 0 \\ 0 \\ \frac{\pi}{2} \end{bmatrix}.

Uses

The axis–angle representation is convenient when dealing with rigid body dynamics. It is useful to both characterize rotations, and also for converting between different representations of rigid body motion, such as homogeneous transformations and twists.

When a rigid body rotates around a fixed axis, its axis–angle data are a constant rotation axis and the rotation angle continuously dependent on time.

Rotating a vector

Main article: Rodrigues' rotation formula

The Rodrigues' rotation formula, named after Olinde Rodrigues, is an efficient algorithm for rotating a Euclidean vector, given a rotation axis and an angle of rotation. In other words, the Rodrigues' formula provides an algorithm to compute the exponential map from $so (3)$ to $SO(3)$ without computing the full matrix exponential.

If $v$ is a vector in $ℝ 3$ and $e$ is a unit vector describing an axis of rotation about which $v$ is rotated by an angle $θ$ the Rodrigues' rotation formula to obtain the rotated vector is

\mathbf{v}_\mathrm{rot} = (\cos\theta) \mathbf{v} + (\sin\theta) (\mathbf{e} \times \mathbf{v}) + (1 - \cos\theta) (\mathbf{e} \cdot \mathbf{v}) \mathbf{e}~.

For the rotation of a single vector it may be more efficient than converting $e$ and $θ$ into a rotation matrix to rotate the vector.

Relationship to other representations

Further information: Charts on SO(3)

There are several ways to represent a rotation. It is useful to understand how different representations relate to one another, and how to convert between them. Here the unit vector is denoted $ω$ instead of $e$ .

Exponential map from so(3) to SO(3)

Further information: Matrix exponential, Orthogonal matrix, Lie algebra § Relation to Lie groups and Rotation group SO(3) § Exponential map

The exponential map effects a transformation from the axis-angle representation of rotations to rotation matrices,

\exp\colon \mathfrak{so}(3) \to \mathrm{SO}(3) ~.

Essentially, by using a Taylor expansion one derives a closed form relation between these two representations. Given a unit vector ω ∈ $\mathfrak{so}$ (3) = ℝ³ representing the unit rotation axis, and an angle, θ ∈ ℝ, an equivalent rotation matrix R is given as follows, where K is the cross product matrix of ω.

That is, K v = ω × v for all vectors v ∈ ℝ³,

R = \exp(\theta \mathbf{K}) = \sum_{k=0}^\infty\frac{(\theta \mathbf{K})^k}{k!} = I + \mathbf{K} \theta + \frac{1}{2!}(\theta \mathbf{K})^2 + \frac{1}{3!}(\theta \mathbf{K})^3 + \cdots

Because K is skew-symmetric, and the sum of the squares of its above-diagonal entries is 1, the characteristic polynomial $P (t)$ of K is $P (t) = det(K - t I) = -(t 3 + t)$ . Since, by the Hamilton-Cayley theorem, $P (K)$ = 0, this implies that

K 3 = - K

As a result, K⁴ = –K², K⁵ = K, K⁶ = K², K⁷ = –K .

This cyclic pattern continues indefinitely, and so all higher powers of K can be expressed in terms of K and K². Thus, from the above equation, it follows that

R = I + \left(\theta - \frac{\theta^3}{3!} + \frac{\theta^5}{5!} - \cdots\right) \mathbf{K} + \left(\frac{\theta^2}{2!} - \frac{\theta^4}{4!} + \frac{\theta^6}{6!} - \cdots\right) \mathbf{K}^2 ~,

that is,

R = I + \sin(\theta) \mathbf{K} + (1-\cos(\theta)) \mathbf{K}^2.

This is a Lie-algebraic derivation, in contrast to the geometric one in the article Rodrigues' rotation formula.^[1]

Due to the existence of the above-mentioned exponential map, the unit vector ω representing the rotation axis, and the angle θ are sometimes called the exponential coordinates of the rotation matrix R.

Log map from SO(3) to so(3)

Further information: Rotation group SO(3) and Infinitesimal transformation

Let K continue to denote the 3×3 matrix that effects the cross product with the rotation axis ω: K(v) = ω × v for all vectors v in what follows.

To retrieve the axis–angle representation of a rotation matrix, calculate the angle of rotation from the trace of the rotation matrix

\theta = \arccos\left( \frac{\mathrm{trace}(R) - 1}{2} \right)

and then use that to find the normalized axis,

\mathbf{\omega} = \frac{1}{2 \sin(\theta)} \begin{bmatrix} R(3,2)-R(2,3) \\ R(1,3)-R(3,1) \\ R(2,1)-R(1,2) \end{bmatrix} ~.

Note also that the Matrix logarithm of the rotation matrix R is

\log R = \left\{ \begin{matrix} 0 & \mathrm{if} \; \theta = 0 \\ \frac{\theta}{2 \sin(\theta)} (R - R^\mathsf{T}) & \mathrm{if} \; \theta \ne 0 \; \mathrm{and} \; \theta \in (-\pi, \pi) \end{matrix}\right.

An exception occurs when R has eigenvalues equal to −1. In this case, the log is not unique. However, even in the case where θ=π the Frobenius norm of the log is

\| \log(R) \|_F = \sqrt{2} | \theta | ~.

Given rotation matrices A and B,

d_g(A,B) := \| \log(A^\mathsf{T} B)\|_F

is the geodesic distance on the 3D manifold of rotation matrices.

For small rotations, the above computation of θ may be numerically imprecise as the derivative of arccos goes to infinity as θ→0. In that case, the off-axis terms will actually provide better information about θ since, for small angles, $R \approx I + θ K$ . (This is because these are the first two terms of the Taylor series for exp(θ K).)

This formulation also has numerical problems at θ=π, where the off-axis terms don't give information about the rotation axis (which is still defined up to a sign ambiguity). In that case, we must reconsider the above formula.

R = I + \mathbf{K} \sin(\theta) + \mathbf{K}^2 (1-\cos(\theta))

At θ=π, we have

R = I + 2 \mathbf{K}^2 = I + 2(\mathbf{\omega} \otimes \mathbf{\omega} - I) = 2 \mathbf{\omega} \otimes \mathbf{\omega} - I

and so let

B:=\mathbf{\omega} \otimes \mathbf{\omega} = \frac{1}{2}(R+I) ~,

so the diagonal terms of B are the squares of the elements of ω and the signs (up to sign ambiguity) can be determined from the signs of the off-axis terms of B.

Unit quaternions

Main article: Quaternions and spatial rotation

the following expression transforms axis–angle coordinates to versors (unit quaternions):

Q = \left(\cos\left(\frac{\theta}{2}\right), \mathbf{\omega} \sin\left(\frac{\theta}{2}\right)\right)

Given a versor $q = s + \mathbf{x}$ represented with its scalar s and vector x, the axis–angle coordinates can be extracted using the following:

\theta = 2\,\arccos(s)\,

\mathbf{\omega} = \left\{ \begin{matrix} \frac{x}{ \sin( \theta/2 ) } , & \mathrm{if} \; \theta \neq 0 \\ 0, & \mathrm{otherwise} \end{matrix} \right.

A more numerically stable expression of the rotation angle uses the atan2 function:

\theta = 2\,\operatorname{atan2}(|x|,s)\,,

where |x| is the Euclidean norm of the 3-vector x.

References

↑ This holds for the triplet representation of the rotation group, i.e., spin 1. For higher dimensional representations/spins, see ref Curtright, T L; Fairlie, D B; Zachos, C K (2014). "A compact formula for rotations as spin matrix polynomials". SIGMA 10: 084. doi:10.3842/SIGMA.2014.084.

This article is issued from Wikipedia - version of the Tuesday, January 12, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.