Homogeneous coordinates

From Example Problems
Jump to: navigation, search

In mathematics, homogeneous coordinates, introduced by August Ferdinand Möbius, allow affine transformations to be easily represented by a matrix. Also they make calculations possible in projective space just as Cartesian coordinates do in Euclidean space. The homogeneous coordinates of a point of projective space of dimension n are usually written as (x : y : z : ... : w), a row vector of length n + 1, other than (0 : 0 : 0 : ... : 0). Two sets of coordinates that are proportional denote the same point of projective space: for any non-zero scalar c from the underlying field K, (cx : cy : cz : ... : cw) denotes the same point. Therefore this system of coordinates can be explained as follows: if the projective space is constructed from a vector space V of dimension n + 1, introduce coordinates in V by choosing a basis, and use these in P(V), the equivalence classes of proportional non-zero vectors in V.

Taking the example of projective space of dimension three, there will be homogeneous coordinates (x : y : z : w). The plane at infinity is usually identified with the set of points with w = 0. Away from this plane we can use (x/w, y/w, z/w) as an ordinary Cartesian system; therefore the affine space complementary to the plane at infinity is coordinatised in a familiar way, with a basis corresponding to (1 : 0 : 0 : 1), (0 : 1 : 0 : 1), (0 : 0 : 1 : 1).

If we try to intersect the two planes defined by equations x = w and x = 2w then we clearly will derive first w = 0 and then x = 0. That tells us that the intersection is contained in the plane at infinity, and consists of all points with coordinates (0 : y : z : 0). It is a line, and in fact the line joining (0 : 1 : 0 : 0) and (0 : 0 : 1 : 0). The line is given by the equation

(0:y:z:0)=\mu (1-\lambda )(0:1:0:0)+\mu \lambda (0:0:1:0)

where μ is a scaling factor. The scaling factor can be adjusted to normalize the coordinates (0 : y : z : 0), thereby eliminating one of the two degrees of freedom. The result is a set of points with only one degree of freedom, as is expected for a line.

Brackets versus parentheses

Consider projective 2-space: points in the projective plane are projections of points in 3-space ("3-D points"). Let the notation

(x:y:z)

refer to one of these 3-D points. Let

(u:v:w)

refer to another 3-D point. Then

(x:y:z)=(u:v:w)\leftrightarrow x=u\wedge y=v\wedge z=w.

On the other hand, let the notation

[x:y:z]

refer to the projection of 3-D point (x : y : z) onto the projective plane. The point [x : y : z] can be considered to be equal to an equivalence class of 3-D points which belong to the 3-D line passing through the points (x : y : z) and (0 : 0 : 0). If

[u:v:w]

is another projective point, then

[x:y:z]=[u:v:w]\leftrightarrow \exists \alpha (x=\alpha u\wedge y=\alpha v\wedge z=\alpha w).

Two 3-D points are equivalent if their projections onto the projective plane are equal:

(x:y:z)\equiv (u:v:w)\leftrightarrow \exists \alpha (x=\alpha u\wedge y=\alpha v\wedge z=\alpha w).

Thus,

(x:y:z)\equiv (u:v:w)\leftrightarrow [x:y:z]=[u:v:w].

Addition of homogeneous coordinates

This distinction between brackets and parentheses means that addition of points in homogeneous coordinates will be defined in two different ways, depending on whether the coordinates are enclosed with brackets or parentheses.

Consider once again the case of the projective plane. Addition of a pair of 3-D points is the same as for ordinary coordinates:

(a:b:c)+(x:y:z)=(a+x:b+y:c+z).

On the other hand, addition of a pair of projected points can be defined thus:

[a:b:c]+[x:y:z]=[za+xc:zb+yc:cz].

For projective 3-space, similar considerations apply. Addition of a pair of unprojected points is

(a:b:c:d)+(x:y:z:w)=(a+x:b+y:c+z:d+w)

whereas addition of a pair of projected points is

[a:b:c:d]+[x:y:z:w]=[wa+dx:wb+dy:wc+dz:dw].

Scalar multiplication of homogeneous coordinates

There are two kinds of scalar multiplication: one for unprojected points and another one for projected points.

Consider a scalar a and an unprojected 3-D point (x : y : z). Then

a(x:y:z)=(ax:ay:az).

Notice that

(x:y:z)\equiv a(x:y:z)

even though

(x:y:z)\neq a(x:y:z).

Now consider the scalar a and a projected point [x : y : z]. Then

a[x:y:z]=[ax:ay:z]

so that

[x:y:z]\neq a[x:y:z].

Linear combinations of points described with homogeneous coordinates

Let there be a pair of points A and B in projective 3-space, whose homogeneous coordinates are

{\mathbf  {A}}:[X_{A}:Y_{A}:Z_{A}:W_{A}],
{\mathbf  {B}}:[X_{B}:Y_{B}:Z_{B}:W_{B}].

It is desired to find their linear combination a{\mathbf  {A}}+b{\mathbf  {B}} where a and b are coefficients which can be adjusted at will. There are three cases to consider:

  • both points belong to affine 3-space,
  • both points belong to the plane at infinity,
  • one point is affine and the other one is at infinity.

The X, Y, and Z coordinates can be considered as numerators, whereas the W coordinate can be considered as a denominator. To add homogeneous coordinates it is necessary that the denominator be common. Otherwise it is necessary to rescale the coordinates until all the denominators are common. Homogeneous coordinates are equivalent up to any uniform rescaling.

Both points are affine

If both points are in affine 3-space, then W_{A}\neq 0 and W_{B}\neq 0. Their linear combination is

a[X_{A}:Y_{A}:Z_{A}:W_{A}]+b[X_{B}:Y_{B}:Z_{B}:W_{B}]\
=[aX_{A}:aY_{A}:aZ_{A}:W_{A}]+[bX_{B}:bY_{B}:bZ_{B}:W_{B}]\
=\left[a{X_{A} \over W_{A}}:a{Y_{A} \over W_{A}}:a{Z_{A} \over W_{A}}:1\right]+\left[b{X_{B} \over W_{B}}:b{Y_{B} \over W_{B}}:b{Z_{B} \over W_{B}}:1\right]
=\left[a{X_{A} \over W_{A}}+b{X_{B} \over W_{B}}:a{Y_{A} \over W_{A}}+b{Y_{B} \over W_{B}}:a{Z_{A} \over W_{A}}+b{Z_{B} \over W_{B}}:1\right].

Both points are at infinity

If both points are on the plane at infinity, then WA = 0 and WB = 0. Their linear combination is

a[X_{A}:Y_{A}:Z_{A}:W_{A}]+b[X_{B}:Y_{B}:Z_{B}:W_{B}]=[aX_{A}:aY_{A}:aZ_{A}:0]+[bX_{B}:bY_{B}:bZ_{B}:0]
=[aX_{A}+bX_{B}:aY_{A}+bY_{B}:aZ_{A}+bZ_{B}:0].

One point is affine and the other at infinity

Let the first point be affine, so that W_{A}\neq 0. Then

a[X_{A}:Y_{A}:Z_{A}:W_{A}]+b[X_{B}:Y_{B}:Z_{B}:0]
=a[0:0:0:0]+b[X_{B}:Y_{B}:Z_{B}:0],
=[bX_{B}:bY_{B}:bZ_{B}:0],

which means that the point at infinity is "dominant".

See also: barycentric coördinates.de:Homogene Koordinaten fr:Coordonnées homogènes