Spatial Transformation Matrix Math Methods

Affine spatial transformation matrices are used to represent the orientation and position of a coordinate system within a 3 dimensional "world" coordinate system. Spatial transformation matrices can be used to map from one to another. They are often used in Kinematics.

Any combination of translation, rotation, scaling, reflection and shear can be represented by a single 4 by 4 affine transformation matrix:

A =

a₁₁	a₁₂	a₁₃	x
a₂₁	a₂₂	a₂₃	y
a₃₁	a₃₂	a₃₃	z
0	0	0	1

This matrix A can be used to transform any point / vector p to a new point / vector p'.

DCM (Direction Cosine Matrix): The upper left 3x3 area is the rotation / orientation / scale / shear sub-section. It is made of 3 sections:

X unit vector: The first, leftmost, column [a₁₁, a₂₁, a₃₁] is a unit vector representing the change relative to the X axis.
Y unit vector: The 2nd, center, column [a₁₂, a₂₂, a₃₂] represents change from the Y axis. For example, if it is [0,0,-1] then we have rotated straight down, because a₃₂ is the z component of the y unit vector.
Z unit vector: The 3rd, rightmost, column [a₁₃, a₂₃, a₃₃] is the Z component.

Position: The upper right area translates position.

The 4th row is always [0, 0, 0, 1] to maintain transformation matrix format.

To apply the transformation, simply enter the desired values for the type of transformation (see below) and apply matrix math multiplication to the original point as a vector of x, y, z, and 1.

p' = Ap

Multiple translations can be applied sequentially, or can be merged into a single matrix operation by multiplying all the transformation matrixes together. Remember that unlike normal math, in matrix multiplication the order matters, and as applied here, this maintains the order of transformations.

Translation

A =

1	0	0	x
0	1	0	y
0	0	1	z
0	0	0	1

Where x, y, and z are the desired translation in each direction.

Scaling

A =

Sx	0	0	0
0	Sy	0	0
0	0	Sz	0
0	0	0	1

Where Sx, Sy, and Sz are the desired scaling factor in each dimension. Note that if you want to scale in all three axis by the same amount, you can just change the lower right "1" to factor and leave Sx, Sy, Sz as 1.

Rotation

Where theta O is the desired rotation in radians around the given axis. These can be combined, by multiplying all three matrix together, and separate the "thetas" to three different variables:

Rotation about X

A =

1	0	0	0
0	cos(O)	sin(O)	0
0	-sin(O)	cos(O)	0
0	0	0	1

Rotation about Y

A =

cos(O)	0	-sin(O)	0
0	1	0	0
sin(O)	0	cos(O)	0
0	0	0	1

Rotation about Z

A =

cos(O)	-sin(O)	0	0
sin(O)	cos(O)	0	0
0	0	1	0
0	0	0	1

Note that these are "right handed" rotations. Positive angles of O cause counter-clockwise rotation about an axis as seen from a point on the positive side of the axis looking in toward the origin.

These three transformations (in X, Y, and Z) can, theoretically, be applied in sequence to reach any orientation. This uses a set of three angles, which are referred to as Euler Angles.

Note that if you apply a rotation of 1/4 of a full turn around, for example, the X axis, then you have effectively aligned the other two axis, Y and Z, with each other. If you then apply any rotation around, for example, the Z axis, that would be exactly the same as rotation around the Y axes. This removes a degree of freedom from any further transformations, and is called "gimbal lock" which is one type of singularity. It is a limitation of Euler Angles.

To avoid this limitation, a better system for rotations in multiple axis can be developed by using Quaternions. These are 4 dimensional values where the extra dimensions are akin to complex numbers. Where a complex number has one real and one imaginary part, a quaternion has a vector component and 4 imaginary parts, called i, j, and k.

qw + i qx + j qy + k qz.

A quaternion holds the angles in the real and imaginary parts. Because they are all combined, any rotation around a combination of axis can be described. Quaternions can be converted into 3x3 rotational matrixes.

See also:

https://www.youtube.com/watch?v=htYh-Tq7ZBI An excellent explanation of quaternions via the multiplication of vectors and their byproducts.+
https://youtu.be/d4EgbgTm0Bg Visualizing quaternions
https://www.euclideanspace.com/maths/algebra/realNormedAlgebra/quaternions/transforms/index.htm

https://www.euclideanspace.com/maths/geometry/rotations/conversions/quaternionToMatrix/index.htm
If a quaternion is represented by qw + i qx + j qy + k qz , then the equivalent matrix, to represent the same rotation, is:

1 - 2qy² - 2qz²	2qxqy - 2qzqw	2qxqz + 2qyqw
2qxqy + 2qzqw	1 - 2qx² - 2qz²	2qyqz - 2qxqw
2qxqz - 2qyqw	2qyqz + 2qxqw	1 - 2qx² - 2qy²

file: /Techref/method/math/spatial-transformations.htm, 15KB, , updated: 2023/11/8 19:17, local time: 2024/11/21 13:44,

^{3.138.114.140:LOG IN}

©2024 These pages are served without commercial sponsorship. (No popup ads, etc...).Bandwidth abuse increases hosting cost forcing sponsorship or shutdown. This server aggressively defends against automated copying for any reason including offline viewing, duplication, etc... Please respect this requirement and DO NOT RIP THIS SITE. Questions?

Please DO link to this page! Digg it! / MAKE!

<A HREF="http://techref.massmind.org/techref/method/math/spatial-transformations.htm"> Spatial Transformation Matrix Math Methods</A>

Did you find what you needed?

1 - 2qy² - 2qz²	2qxqy - 2qzqw	2qxqz + 2qyqw
2qxqy + 2qzqw	1 - 2qx² - 2qz²	2qyqz - 2qxqw
2qxqz - 2qyqw	2qyqz + 2qxqw	1 - 2qx² - 2qy²