Affine transformation

From Infogalactic: the planetary knowledge core
Jump to: navigation, search
An image of a fern-like fractal that exhibits affine self-similarity. Each of the leaves of the fern is related to each other leaf by an affine transformation. For instance, the red leaf can be transformed into both the small dark blue leaf and the large light blue leaf by a combination of reflection, rotation, scaling, and translation.

In geometry, an affine transformation, affine map[1] or an affinity (from the Latin, affinis, "connected with") is a function between affine spaces which preserves points, straight lines and planes. Also, sets of parallel lines remain parallel after an affine transformation. An affine transformation does not necessarily preserve angles between lines or distances between points, though it does preserve ratios of distances between points lying on a straight line.

Examples of affine transformations include translation, scaling, homothety, similarity transformation, reflection, rotation, shear mapping, and compositions of them in any combination and sequence.

If X and Y are affine spaces, then every affine transformation f : X \to Y is of the form x \mapsto Mx + b, where M is a linear transformation on X and b is a vector in Y. Unlike a purely linear transformation, an affine map need not preserve the zero point in a linear space. Thus, every linear transformation is affine, but not every affine transformation is linear.

All Euclidean spaces are affine, but there are affine spaces that are non-Euclidean. In affine coordinates, which include Cartesian coordinates in Euclidean spaces, each output coordinate of an affine map is a linear function (in the sense of calculus) of all input coordinates. Another way to deal with affine transformations systematically is to select a point as the origin; then, any affine transformation is equivalent to a linear transformation (of position vectors) followed by a translation.

Mathematical definition

An affine map[1] f:\mathcal{A} \to \mathcal{B} between two affine spaces is a map on the points that acts linearly on the vectors (that is, the vectors between points of the space). In symbols, f determines a linear transformation \varphi such that, for any pair of points P, Q \in \mathcal{A}:

\overrightarrow{f(P)~f(Q)} = \varphi(\overrightarrow{PQ})

or

f(Q)-f(P) = \varphi(Q-P).

We can interpret this definition in a few other ways, as follows.

If an origin O \in \mathcal{A} is chosen, and B denotes its image f(O) \in \mathcal{B}, then this means that for any vector \vec{x}:

f: (O+\vec{x}) \mapsto (B+\varphi(\vec{x})) .

If an origin O' \in \mathcal{B} is also chosen, this can be decomposed as an affine transformation g : \mathcal{A} \to \mathcal{B} that sends O \mapsto O', namely

g: (O+\vec{x}) \mapsto (O'+\varphi(\vec{x})) ,

followed by the translation by a vector \vec{b} = \overrightarrow{O'B}.

The conclusion is that, intuitively, f consists of a translation and a linear map.

Alternative definition

Given two affine spaces \mathcal{A} and \mathcal{B}, over the same field, a function f: \mathcal{A} \to \mathcal{B} is an affine map if and only if for every family \{(a_i, \lambda_i)\}_{i\in I} of weighted points in \mathcal{A} such that

\sum_{i\in I}\lambda_i = 1 ,

we have[2]

f\left(\sum_{i\in I}\lambda_i a_i\right)=\sum_{i\in I}\lambda_i f(a_i)\, .

In other words, f preserves barycenters.

Representation

As shown above, an affine map is the composition of two functions: a translation and a linear map. Ordinary vector algebra uses matrix multiplication to represent linear maps, and vector addition to represent translations. Formally, in the finite-dimensional case, if the linear map is represented as a multiplication by a matrix A and the translation as the addition of a vector \vec{b}, an affine map f acting on a vector \vec{x} can be represented as


\vec{y} = f(\vec{x}) = A \vec{x} + \vec{b}.

Augmented matrix

Affine transformations on the 2D plane can be performed in three dimensions. Translation is done by shearing along over the z axis, and rotation is performed around the z axis.

Using an augmented matrix and an augmented vector, it is possible to represent both the translation and the linear map using a single matrix multiplication. The technique requires that all vectors are augmented with a "1" at the end, and all matrices are augmented with an extra row of zeros at the bottom, an extra column—the translation vector—to the right, and a "1" in the lower right corner. If A is a matrix,


\begin{bmatrix} \vec{y} \\ 1 \end{bmatrix} = \left[ \begin{array}{ccc|c} \, & A & & \vec{b} \ \\ 0 & \ldots & 0 & 1 \end{array} \right] \begin{bmatrix} \vec{x} \\ 1 \end{bmatrix}

is equivalent to the following


\vec{y} = A \vec{x} + \vec{b}.

The above-mentioned augmented matrix is called affine transformation matrix, or projective transformation matrix (as it can also be used to perform projective transformations).

This representation exhibits the set of all invertible affine transformations as the semidirect product of Kn and GL(n, k). This is a group under the operation of composition of functions, called the affine group.

Ordinary matrix-vector multiplication always maps the origin to the origin, and could therefore never represent a translation, in which the origin must necessarily be mapped to some other point. By appending the additional coordinate "1" to every vector, one essentially considers the space to be mapped as a subset of a space with an additional dimension. In that space, the original space occupies the subset in which the additional coordinate is 1. Thus the origin of the original space can be found at (0,0, ... 0, 1). A translation within the original space by means of a linear transformation of the higher-dimensional space is then possible (specifically, a shear transformation). The coordinates in the higher-dimensional space are an example of homogeneous coordinates. If the original space is Euclidean, the higher dimensional space is a real projective space.

The advantage of using homogeneous coordinates is that one can combine any number of affine transformations into one by multiplying the respective matrices. This property is used extensively in computer graphics, computer vision and robotics.

Example augmented matrix

If the vectors \vec{x}_1, \ldots, \vec{x}_{n+1} are a basis of the domain's projective vector space and if \vec{y}_1, \ldots, \vec{y}_{n+1} are the corresponding vectors in the codomain vector space then the augmented matrix M that achieves this affine transformation

\left[\begin{array}{c}\vec{y}\\1\end{array}\right] = M \left[\begin{array}{c}\vec{x}\\1\end{array}\right]

is

M = \left[\begin{array}{ccc}\vec{y}_1&\ldots&\vec{y}_{n+1}\\1&\ldots&1\end{array}\right] \left[\begin{array}{ccc}\vec{x}_1&\ldots&\vec{x}_{n+1}\\1&\ldots&1\end{array}\right]^{-1}.

This formulation works irrespective of whether any of the domain, codomain and image vector spaces have the same number of dimensions.

For example, the affine transformation of a vector plane is uniquely determined from the knowledge of where the three vertices of a non-degenerate triangle are mapped to.

Properties

An affine transformation preserves:

  1. The collinearity relation between points; i.e., points which lie on the same line (called collinear points) continue to be collinear after the transformation.
  2. Ratios of vectors along a line; i.e., for distinct collinear points p_1,\, p_2,\, p_3, the ratio of \overrightarrow{p_1p_2} and \overrightarrow{p_2p_3} is the same as that of \overrightarrow{f(p_1)f(p_2)} and \overrightarrow{f(p_2)f(p_3)}.
  3. More generally barycenters of weighted collections of points.

An affine transformation is invertible if and only if A is invertible. In the matrix representation, the inverse is:


\left[ \begin{array}{ccc|c} & A^{-1} & & -A^{-1}\vec{b} \ \\ 0 & \ldots & 0 & 1 \end{array} \right]

The invertible affine transformations (of an affine space onto itself) form the affine group, which has the general linear group of degree n as subgroup and is itself a subgroup of the general linear group of degree n + 1.

The similarity transformations form the subgroup where A is a scalar times an orthogonal matrix. For example, if the affine transformation acts on the plane and if the determinant of A is 1 or −1 then the transformation is an equiareal mapping. Such transformations form a subgroup called the equi-affine group[3] A transformation that is both equi-affine and a similarity is an isometry of the plane taken with Euclidean distance.

Each of these groups has a subgroup of transformations which preserve orientation: those where the determinant of A is positive. In the last case this is in 3D the group of rigid body motions (proper rotations and pure translations).

If there is a fixed point, we can take that as the origin, and the affine transformation reduces to a linear transformation. This may make it easier to classify and understand the transformation. For example, describing a transformation as a rotation by a certain angle with respect to a certain axis may give a clearer idea of the overall behavior of the transformation than describing it as a combination of a translation and a rotation. However, this depends on application and context.

Affine transformation of the plane

A central dilation. The triangles A1B1Z, A1C1Z, and B1C1Z get mapped to A2B2Z, A2C2Z, and B2C2Z, respectively.

Affine transformations in two real dimensions include:

  • pure translations,
  • scaling in a given direction, with respect to a line in another direction (not necessarily perpendicular), combined with translation that is not purely in the direction of scaling; taking "scaling" in a generalized sense it includes the cases that the scale factor is zero (projection) or negative; the latter includes reflection, and combined with translation it includes glide reflection,
  • rotation combined with a homothety and a translation,
  • shear mapping combined with a homothety and a translation, or
  • squeeze mapping combined with a homothety and a translation.

To visualise the general affine transformation of the Euclidean plane, take labelled parallelograms ABCD and A′B′C′D′. Whatever the choices of points, there is an affine transformation T of the plane taking A to A′, and each vertex similarly. Supposing we exclude the degenerate case where ABCD has zero area, there is a unique such affine transformation T. Drawing out a whole grid of parallelograms based on ABCD, the image T(P) of any point P is determined by noting that T(A) = A′, T applied to the line segment AB is A′B′, T applied to the line segment AC is A′C′, and T respects scalar multiples of vectors based at A. [If A, E, F are collinear then the ratio length(AF)/length(AE) is equal to length(AF′)/length(AE′).] Geometrically T transforms the grid based on ABCD to that based in A′B′C′D′.

Affine transformations don't respect lengths or angles; they multiply area by a constant factor

area of A′B′C′D′ / area of ABCD.

A given T may either be direct (respect orientation), or indirect (reverse orientation), and this may be determined by its effect on signed areas (as defined, for example, by the cross product of vectors).

Examples of affine transformations

Affine transformations over the real numbers

Functions f : RR, f(x) = mx + c with m and c constant, are commonplace affine transformations.

Affine transformation over a finite field

The following equation expresses an affine transformation in GF(28):


\{\,a'\,\} = M\{\,a\,\} \oplus \{\,v\,\},
where [M] is the matrix and {v} is the vector
M\{\,a\,\}=
\begin{bmatrix}
1&0&0&0&1&1&1&1 \\
1&1&0&0&0&1&1&1 \\
1&1&1&0&0&0&1&1 \\
1&1&1&1&0&0&0&1 \\
1&1&1&1&1&0&0&0 \\
0&1&1&1&1&1&0&0 \\
0&0&1&1&1&1&1&0 \\
0&0&0&1&1&1&1&1
\end{bmatrix}
 :\{\,v\,\}= \begin{bmatrix} 1 \\ 1 \\ 0 \\ 0 \\ 0 \\ 1 \\ 1 \\ 0 \end{bmatrix}.

For instance, the affine transformation of the element {a} = y7 + y6 + y3 + y = {11001010} in big-endian binary notation = {CA} in big-endian hexadecimal notation, is calculated as follows:

a_0' = a_0 \oplus a_4 \oplus a_5 \oplus a_6 \oplus a_7 \oplus 1 = 0 \oplus 0 \oplus 0 \oplus 1 \oplus 1 \oplus 1 = 1
a_1' = a_0 \oplus a_1 \oplus a_5 \oplus a_6 \oplus a_7 \oplus 1 = 0 \oplus 1 \oplus 0 \oplus 1 \oplus 1 \oplus 1 = 0
a_2' = a_0 \oplus a_1 \oplus a_2 \oplus a_6 \oplus a_7 \oplus 0 = 0 \oplus 1 \oplus 0 \oplus 1 \oplus 1 \oplus 0 = 1
a_3' = a_0 \oplus a_1 \oplus a_2 \oplus a_3 \oplus a_7 \oplus 0 = 0 \oplus 1 \oplus 0 \oplus 1 \oplus 1 \oplus 0 = 1
a_4' = a_0 \oplus a_1 \oplus a_2 \oplus a_3 \oplus a_4 \oplus 0 = 0 \oplus 1 \oplus 0 \oplus 1 \oplus 0 \oplus 0 = 0
a_5' = a_1 \oplus a_2 \oplus a_3 \oplus a_4 \oplus a_5 \oplus 1 = 1 \oplus 0 \oplus 1 \oplus 0 \oplus 0 \oplus 1 = 1
a_6' = a_2 \oplus a_3 \oplus a_4 \oplus a_5 \oplus a_6 \oplus 1 = 0 \oplus 1 \oplus 0 \oplus 0 \oplus 1 \oplus 1 = 1
a_7' = a_3 \oplus a_4 \oplus a_5 \oplus a_6 \oplus a_7 \oplus 0 = 1 \oplus 0 \oplus 0 \oplus 1 \oplus 1 \oplus 0 = 1.

Thus, {a′} = y7 + y6 + y5 + y3 + y2 + 1 = {11101101} = {ED}.

Affine transformation in plane geometry

A simple affine transformation on the real plane
Effect of applying various 2D affine transformation matrices on a unit square. Note that the reflection matrices are special cases of the scaling matrix.

In ℝ2, the transformation shown at left is accomplished using the map given by:

\begin{bmatrix} x \\ y\end{bmatrix} \mapsto \begin{bmatrix} 0&1\\ 2&1 \end{bmatrix}\begin{bmatrix} x \\ y\end{bmatrix} + \begin{bmatrix} -100 \\ -100\end{bmatrix}

Transforming the three corner points of the original triangle (in red) gives three new points which form the new triangle (in blue). This transformation skews and translates the original triangle.

In fact, all triangles are related to one another by affine transformations. This is also true for all parallelograms, but not for all quadrilaterals.

See also

Notes

  1. 1.0 1.1 Berger, Marcel (1987), p. 38.
  2. Lua error in package.lua at line 80: module 'strict' not found.
  3. Oswald Veblen (1918) Projective Geometry, volume 2, pp. 105–7.

References

  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.

External links