Vector rotation

Vector rotations are widely used not only in the sciences, such as physics, chemistry and mathematics, but are critical for graphics computations in video game programs and navigation in space. A typical example used in computer games would be calculating the graphics for a military tank rolling up a slanted hill, the relative rotation of the tank's turret, and the elevation of the tanks' barrel. Although a rotation matrix for each point of the tank could be calculated individually, a more economical method is to calculate a single rotation matrix for the entire tank and apply that solution to every current position of the tank as it rolls up the hill. Additional rotations are then used for the turret rotation by a second multiplication.

Methods

A variety of methods can be used to determine the rotation matrix (in 3D or 4D space) needed to convert vector V₁ into vector V₂, see rotation matrix and Euler angles for two alternative methods. Because they may lead to non-unique results, inverse trigonometry functions, such as arctan, should only be used with great caution. Quaternions, a 4-dimensional approach in 3D space, can also be used, and this method has devoted followers and critics. Although several 3D matrix rotation methods can be used, the method of Hughes (J. Graphics Tools, 2000) is particularly fast, because it avoids time-consuming inverse trigonometry and square root calculations, and avoids computational pitfalls of instability inherent to some of the previous methods.

Method of Möller and Hughes

The product of two reflections is a rotation, and the Möller-Hughes method relies on this well-known fact, using the reflections defined by the Householder matrix to construct a rotation matrix by multiplying two reflections about either the X, Y or Z-axis, using the axis most orthogonal to the starting vector. It is important to note that this method only works for unit vectors, so translation and scaling may be required before performing the vector rotation.

When the unit vector V₁ and desired final unit direction vector V₂ defined as:

$\,V_{1}=(V_{1x},V_{1y},V_{1z})$

$\,V_{2}=(V_{2x},V_{2y},V_{2z})$

are almost parallel then we choose the reflection vector P as the unit vector most orthogonal to V₁ in the following manner:

$P=(1,0,0)\quad {\hbox{if}}\quad |V_{1x}|<|V_{1y}|\quad {\hbox{and}}\quad |V_{1x}|<|V_{1z}|$

$P=(0,1,0)\quad {\hbox{if}}\quad |V_{1y}|<|V_{1x}|\quad {\hbox{and}}\quad |V_{1y}|<|V_{1z}|$

$P=(0,0,1)\quad {\hbox{if}}\quad |V_{1z}|<|V_{1y}|\quad {\hbox{and}}\quad |V_{1z}|<|V_{1x}|$

Define two new vectors u and v:

$\mathbf {u}$ = P - $V_{1}$

$\mathbf {v}$ = P - $V_{2}$

Then the elements $(r_{ij})$ of the rotation matrix R which rotates V₁ into V₂ (that is R(V₁) = V₂) can be calculated simply by multiplication and dot products, without using either inverse trigonometry or square roots, as:

$r_{ij}=\delta _{ij}-{\frac {2u_{i}u_{j}}{\mathbf {u} \cdot \mathbf {u} }}-{\frac {2v_{i}v_{j}}{\mathbf {v} \cdot \mathbf {v} }}+{\frac {4\mathbf {u} \cdot \mathbf {v} }{(\mathbf {u} \cdot \mathbf {u} )(\mathbf {v} \cdot \mathbf {v} )}}v_{i}\,u_{j}$

where $\delta _{ij}$ = 1 if i=j, otherwise $\delta _{ij}$ = 0

Product operator formalism

In NMR spectroscopy, magnetic resonance imaging and related sciences, a simplified form of vector rotation, the product operator formalism, can be used because all vectors are unit vectors starting from the origin (0,0,0). This formalism is a simplied form of matrix rotation.

References

Tomas Möller and John F. Hughes, J. Graphics Tools, 4(4):1-4, 1999.