Affine transformation

From HandWiki
Short description: Geometric transformation that preserves lines but not angles nor the origin
An image of a fern-like fractal (Barnsley's fern) 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 dark blue leaf and any of the light blue leaves by a combination of reflection, rotation, scaling, and translation.

In Euclidean geometry, an affine transformation or affinity (from the Latin, affinis, "connected with") is a geometric transformation that preserves lines and parallelism, but not necessarily Euclidean distances and angles.

More generally, an affine transformation is an automorphism of an affine space (Euclidean spaces are specific affine spaces), that is, a function which maps an affine space onto itself while preserving both the dimension of any affine subspaces (meaning that it sends points to points, lines to lines, planes to planes, and so on) and the ratios of the lengths of parallel line segments. Consequently, sets of parallel affine subspaces 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.

If X is the point set of an affine space, then every affine transformation on X can be represented as the composition of a linear transformation on X and a translation of X. Unlike a purely linear transformation, an affine transformation need not preserve the origin of the affine space. Thus, every linear transformation is affine, but not every affine transformation is linear.

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

Viewing an affine space as the complement of a hyperplane at infinity of a projective space, the affine transformations are the projective transformations of that projective space that leave the hyperplane at infinity invariant, restricted to the complement of that hyperplane.

A generalization of an affine transformation is an affine map[1] (or affine homomorphism or affine mapping) between two (potentially different) affine spaces over the same field k. Let (X, V, k) and (Z, W, k) be two affine spaces with X and Z the point sets and V and W the respective associated vector spaces over the field k. A map f: XZ is an affine map if there exists a linear map mf : VW such that mf (xy) = f (x) − f (y) for all x, y in X.[2]

Definition

Let X be an affine space over a field k, and V be its associated vector space. An affine transformation is a bijection f from X onto itself that is an affine map; this means that a linear map g from V to V is well defined by the equation [math]\displaystyle{ g(y-x) =f(y)-f(x); }[/math] here, as usual, the subtraction of two points denotes the free vector from the second point to the first one, and "well-defined" means that [math]\displaystyle{ y-x= y'-x' }[/math] implies that [math]\displaystyle{ f(y)-f(x)=f(y')-f(x'). }[/math]

If the dimension of X is at least two, a semiaffine transformation f of X is a bijection from X onto itself satisfying:[3]

  1. For every d-dimensional affine subspace S of X, then f (S) is also a d-dimensional affine subspace of X.
  2. If S and T are parallel affine subspaces of X, then f (S) and f (T) are parallel.

These two conditions are satisfied by affine transformations, and express what is precisely meant by the expression that "f preserves parallelism".

These conditions are not independent as the second follows from the first.[4] Furthermore, if the field k has at least three elements, the first condition can be simplified to: f is a collineation, that is, it maps lines to lines.[5]

Structure

By the definition of an affine space, V acts on X, so that, for every pair (x, v) in X × V there is associated a point y in X. We can denote this action by v(x) = y. Here we use the convention that v = v are two interchangeable notations for an element of V. By fixing a point c in X one can define a function mc : XV by mc(x) = cx. For any c, this function is one-to-one, and so, has an inverse function mc−1 : VX given by mc−1(v) = v(c). These functions can be used to turn X into a vector space (with respect to the point c) by defining:[6]

  • [math]\displaystyle{ x + y = m_c^{-1}\left(m_c(x) + m_c(y)\right),\text{ for all } x,y \text{ in } X, }[/math] and
  • [math]\displaystyle{ rx = m_c^{-1}\left(r m_c(x)\right), \text{ for all } r \text{ in } k \text{ and } x \text{ in } X. }[/math]

This vector space has origin c and formally needs to be distinguished from the affine space X, but common practice is to denote it by the same symbol and mention that it is a vector space after an origin has been specified. This identification permits points to be viewed as vectors and vice versa.

For any linear transformation λ of V, we can define the function L(c, λ) : XX by

[math]\displaystyle{ L(c, \lambda)(x) = m_c^{-1}\left(\lambda (m_c (x))\right) = c + \lambda (\vec{cx}). }[/math]

Then L(c, λ) is an affine transformation of X which leaves the point c fixed.[7] It is a linear transformation of X, viewed as a vector space with origin c.

Let σ be any affine transformation of X. Pick a point c in X and consider the translation of X by the vector [math]\displaystyle{ \bold{w} = \overrightarrow{c \sigma (c)} }[/math], denoted by Tw. Translations are affine transformations and the composition of affine transformations is an affine transformation. For this choice of c, there exists a unique linear transformation λ of V such that[8]

[math]\displaystyle{ \sigma (x) = T_{\bold{w}} \left( L(c, \lambda)(x) \right). }[/math]

That is, an arbitrary affine transformation of X is the composition of a linear transformation of X (viewed as a vector space) and a translation of X.

This representation of affine transformations is often taken as the definition of an affine transformation (with the choice of origin being implicit).[9][10][11]

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 an invertible matrix [math]\displaystyle{ A }[/math] and the translation as the addition of a vector [math]\displaystyle{ \mathbf{b} }[/math], an affine map [math]\displaystyle{ f }[/math] acting on a vector [math]\displaystyle{ \mathbf{x} }[/math] can be represented as

[math]\displaystyle{ \mathbf{y} = f(\mathbf{x}) = A \mathbf{x} + \mathbf{b}. }[/math]

Augmented matrix

File:Affine transformations.ogv

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 be augmented with a "1" at the end, and all matrices be 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 [math]\displaystyle{ A }[/math] is a matrix,

[math]\displaystyle{ \begin{bmatrix} \mathbf{y} \\ 1 \end{bmatrix} = \left[ \begin{array}{ccc|c} & A & & \mathbf{b} \\ 0 & \cdots & 0 & 1 \end{array} \right] \begin{bmatrix} \mathbf{x} \\ 1 \end{bmatrix} }[/math]

is equivalent to the following

[math]\displaystyle{ \mathbf{y} = A \mathbf{x} + \mathbf{b}. }[/math]

The above-mentioned augmented matrix is called an affine transformation matrix. In the general case, when the last row vector is not restricted to be [math]\displaystyle{ \left[ \begin{array}{ccc|c} 0 & \cdots & 0 & 1 \end{array} \right] }[/math], the matrix becomes a 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 [math]\displaystyle{ K^n }[/math] and [math]\displaystyle{ \operatorname{GL}(n, K) }[/math]. 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 [math]\displaystyle{ (0,0, \dotsc, 0, 1) }[/math]. 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

Suppose you have three points that define a non-degenerate triangle in a plane, or four points that define a non-degenerate tetrahedron in 3-dimensional space, or generally n + 1 points x1, ..., xn+1 that define a non-degenerate simplex in n-dimensional space. Suppose you have corresponding destination points y1, ..., yn+1, where these new points can lie in a space with any number of dimensions. (Furthermore, the new points need not be distinct from each other and need not form a non-degenerate simplex.) The unique augmented matrix M that achieves the affine transformation [math]\displaystyle{ \begin{bmatrix}\mathbf{y}\\1\end{bmatrix} = M \begin{bmatrix}\mathbf{x}\\1\end{bmatrix} }[/math] is [math]\displaystyle{ M = \begin{bmatrix}\mathbf{y}_1&\cdots&\mathbf{y}_{n+1}\\1&\cdots&1\end{bmatrix} \begin{bmatrix}\mathbf{x}_1&\cdots&\mathbf{x}_{n+1}\\1&\cdots&1\end{bmatrix}^{-1}. }[/math]

Properties

Properties preserved

An affine transformation preserves:

  1. collinearity between points: three or more points which lie on the same line (called collinear points) continue to be collinear after the transformation.
  2. parallelism: two or more lines which are parallel, continue to be parallel after the transformation.
  3. convexity of sets: a convex set continues to be convex after the transformation. Moreover, the extreme points of the original set are mapped to the extreme points of the transformed set.[12]
  4. ratios of lengths of parallel line segments: for distinct parallel segments defined by points [math]\displaystyle{ p_1 }[/math] and [math]\displaystyle{ p_2 }[/math], [math]\displaystyle{ p_3 }[/math] and [math]\displaystyle{ p_4 }[/math], the ratio of [math]\displaystyle{ \overrightarrow{p_1p_2} }[/math] and [math]\displaystyle{ \overrightarrow{p_3p_4} }[/math] is the same as that of [math]\displaystyle{ \overrightarrow{f(p_1)f(p_2)} }[/math] and [math]\displaystyle{ \overrightarrow{f(p_3)f(p_4)} }[/math].
  5. barycenters of weighted collections of points.

Groups

As an affine transformation is invertible, the square matrix [math]\displaystyle{ A }[/math] appearing in its matrix representation is invertible. The matrix representation of the inverse transformation is thus

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

The invertible affine transformations (of an affine space onto itself) form the affine group, which has the general linear group of degree [math]\displaystyle{ n }[/math] as subgroup and is itself a subgroup of the general linear group of degree [math]\displaystyle{ n + 1 }[/math].

The similarity transformations form the subgroup where [math]\displaystyle{ A }[/math] is a scalar times an orthogonal matrix. For example, if the affine transformation acts on the plane and if the determinant of [math]\displaystyle{ A }[/math] is 1 or −1 then the transformation is an equiareal mapping. Such transformations form a subgroup called the equi-affine group.[13] 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 orientation-preserving or positive affine transformations: those where the determinant of [math]\displaystyle{ A }[/math] is positive. In the last case this is in 3D the group of rigid transformations (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 maps

An affine map [math]\displaystyle{ f\colon\mathcal{A} \to \mathcal{B} }[/math] 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, [math]\displaystyle{ f }[/math] determines a linear transformation [math]\displaystyle{ \varphi }[/math] such that, for any pair of points [math]\displaystyle{ P, Q \in \mathcal{A} }[/math]:

[math]\displaystyle{ \overrightarrow{f(P)~f(Q)} = \varphi(\overrightarrow{PQ}) }[/math]

or

[math]\displaystyle{ f(Q)-f(P) = \varphi(Q-P) }[/math].

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

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

[math]\displaystyle{ f\colon (O+\vec{x}) \mapsto (B+\varphi(\vec{x})) }[/math].

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

[math]\displaystyle{ g\colon (O+\vec{x}) \mapsto (O'+\varphi(\vec{x})) }[/math],

followed by the translation by a vector [math]\displaystyle{ \vec{b} = \overrightarrow{O'B} }[/math].

The conclusion is that, intuitively, [math]\displaystyle{ f }[/math] consists of a translation and a linear map.

Alternative definition

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

[math]\displaystyle{ \sum_{i\in I}\lambda_i = 1 }[/math],

we have[14]

[math]\displaystyle{ f\left(\sum_{i\in I}\lambda_i a_i\right)=\sum_{i\in I}\lambda_i f(a_i) }[/math].

In other words, [math]\displaystyle{ f }[/math] preserves barycenters.

History

The word "affine" as a mathematical term is defined in connection with tangents to curves in Euler's 1748 Introductio in analysin infinitorum.[15] Felix Klein attributes the term "affine transformation" to Möbius and Gauss.[10]

Image transformation

In their applications to digital image processing, the affine transformations are analogous to printing on a sheet of rubber and stretching the sheet's edges parallel to the plane. This transform relocates pixels requiring intensity interpolation to approximate the value of moved pixels, bicubic interpolation is the standard for image transformations in image processing applications. Affine transformations scale, rotate, translate, mirror and shear images as shown in the following examples:[16]

Transformation name Affine matrix Example
Identity (transform to original image) [math]\displaystyle{ \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix} }[/math] Checkerboard identity.svg
Translation [math]\displaystyle{ \begin{bmatrix} 1 & 0 & v_x \gt 0\\ 0 & 1 & v_y = 0\\ 0 & 0 & 1 \end{bmatrix} }[/math] Checkerboard identity.svg
Reflection [math]\displaystyle{ \begin{bmatrix} -1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix} }[/math] Checkerboard reflection.svg
Scale [math]\displaystyle{ \begin{bmatrix} c_x=2 & 0 & 0 \\ 0 & c_y=1 & 0 \\ 0 & 0 & 1 \end{bmatrix} }[/math] Checkerboard scale.svg
Rotate [math]\displaystyle{ \begin{bmatrix} \cos(\theta) & -\sin(\theta) & 0 \\ \sin(\theta) & \cos(\theta) & 0 \\ 0 & 0 & 1 \end{bmatrix} }[/math] Checkerboard rotate.svg
where θ = π/6 =30°
Shear [math]\displaystyle{ \begin{bmatrix} 1 & c_x=0.5 & 0 \\ c_y=0 & 1 & 0 \\ 0 & 0 & 1 \end{bmatrix} }[/math] Checkerboard shear.svg

The affine transforms are applicable to the registration process where two or more images are aligned (registered). An example of image registration is the generation of panoramic images that are the product of multiple images stitched together.

Affine warping

The affine transform preserves parallel lines. However, the stretching and shearing transformations warp shapes, as the following example shows:

White on black circle image 256 by 256.png Affine transform sheared circle.png

This is an example of image warping. However, the affine transformations do not facilitate projection onto a curved surface or radial distortions.

In 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 do not 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

Over the real numbers

The functions [math]\displaystyle{ f\colon \R \to \R,\; f(x) = mx + c }[/math] with [math]\displaystyle{ m }[/math] and [math]\displaystyle{ c }[/math] in [math]\displaystyle{ \R }[/math] and [math]\displaystyle{ m\ne 0 }[/math], are precisely the affine transformations of the real line.

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 [math]\displaystyle{ \mathbb{R}^2 }[/math], the transformation shown at left is accomplished using the map given by:

[math]\displaystyle{ \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} }[/math]

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. Berger 1987, p. 38.
  2. Samuel 1988, p. 11.
  3. Snapper & Troyer 1989, p. 65.
  4. Snapper & Troyer 1989, p. 66.
  5. Snapper & Troyer 1989, p. 69.
  6. Snapper & Troyer 1989, p. 59.
  7. Snapper & Troyer 1989, p. 76,87.
  8. Snapper & Troyer 1989, p. 86.
  9. Wan 1993, pp. 19-20.
  10. 10.0 10.1 Klein 1948, p. 70.
  11. Brannan, Esplen & Gray 1999, p. 53.
  12. Reinhard Schultz. "Affine transformations and convexity". http://math.ucr.edu/~res/math145A-2014/affine+convex.pdf. 
  13. Oswald Veblen (1918) Projective Geometry, volume 2, pp. 105–7.
  14. Schneider, Philip K.; Eberly, David H. (2003). Geometric Tools for Computer Graphics. Morgan Kaufmann. p. 98. ISBN 978-1-55860-594-7. https://books.google.com/books?id=3Q7HGBx1uLIC&pg=PA98. 
  15. Euler, Leonhard (1748) (in la). Introductio in analysin infinitorum. II. https://gallica.bnf.fr/ark:/12148/bpt6k33529/f240.image.  Book II, sect. XVIII, art. 442
  16. Gonzalez, Rafael (2008). 'Digital Image Processing, 3rd'. Pearson Hall. ISBN 9780131687288. 

References

External links