Projective Geometry/Classic/Projective Transformations/Transformations of the projective line

Projective transformation

Let X be a point on the x-axis. A projective transformation can be defined geometrically for this line by picking a pair of points P, Q, and a line m, all within the same x-y plane which contains the x-axis upon which the transformation will be performed.

Points P and Q represent two different observers, or points of view. Point R is the position of some object they are observing. Line m is the objective world which they are observing, and the x-axis is the subjective perception of m.

Draw line l through points P and X. Line l crosses line m at point R. Then draw line t through points Q and R: line t will cross the x-axis at point T. Point T is the transform of point X [Paiva].

Analysis

The above is a synthetic description of a one-dimensional projective transformation. It is now desired to convert it to an analytical (Cartesian) description.

Let point X have coordinates (x₀,0). Let point P have coordinates $(P_x,P_y)$. Let point Q have coordinates $(Q_x,Q_y)$. Let line m have slope m (m is being overloaded in meaning).

The slope of line l is

$\frac{P_y}{P_x-x_0,}$

so an arbitrary point (x,y) on line l is given by the equation

$\frac{y}{x-x_0}=\frac{P_y}{P_x-x_0}$,

$y=\frac{P_y}{P_x-x_0}(x-x_0).(1)$

On the other hand, any point (x,y) on line m is described by

$y=mx+b.(2)$

The intersection of lines l and m is point R, and it is obtained by combining equations (1) and (2):

$mx+b=\frac{P_{y}x}{P_x-x_0}-\frac{P_y{x}_0}{P_x-x_0}.$

Joining the x terms yields

$(\frac{P_y}{P_x-x_0}-m)x=b+\frac{P_y{x}_0}{P_x-x_0}$

and solving for x we obtain

$x_1=\frac{b(P_x-x_0)+P_y{x}_0}{P_y-m(P_x-x_0)}.$

x₁ is the abscissa of R. The ordinate of R is

$y_1=m\left[\frac{b(P_x-x_0)+P_y{x}_0}{P_y-m(P_x-x_0)}\right]+b.$

Now, knowing both Q and R, the slope of line n is

$\frac{y_1-Q_y}{x_1-Q_x}.$

We want to find the intersection of line n and the x-axis, so let

$(Q_x,Q_y)+\lambda (x_1-Q_x,y_1-Q_y)=(x,0)(3)$

The value of λ must be adjusted so that both sides of vector equation (3) are equal. Equation (3) is actually two equations, one for abscissas and one for ordinates. The one for ordinates is

$Q_y+\lambda (y_1-Q_y)=0$

Solve for lambda,

$\lambda =\frac{-Q_y}{y_1-Q_y}(4)$

The equation for abscissas is

$x=Q_x+\lambda (x_1-Q_x)$

which together with equation (4) yields

$x=Q_x-Q_y\left(\frac{x_1-Q_x}{y_1-Q_y}\right)(5)$

which is the abscissa of T.

Substitute the values of x₁ and y₁ into equation (5),

$x=Q_x-Q_y\left[\frac{\frac{b(P_x-x_0)+P_y{x}_0}{P_y-m(P_x-x_0)}-Q_x}{\frac{mb(P_x-x_0)+m{P}_y{x}_0}{P_y-m(P_x-x_0)}+b-Q_y}\right].$

Dissolve the fractions in both numerator and denominator:

$x=Q_x-Q_y\left[\frac{b(P_x-x_0)+P_y{x}_0-Q_x{P}_y+m{Q}_x(P_x-x_0)}{mb(P_x-x_0)+m{P}_y{x}_0+b{P}_y-mb(P_x-x_0)-Q_y{P}_y+m{Q}_y(P_x-x_0)}\right].$

Simplify and relabel x as t(x):

$t(x)=Q_x-Q_y\left[\frac{(P_x-x_0)(b+m{Q}_x)+P_y(x_0-Q_x)}{(P_x-x_0)m{Q}_y+P_y(m{x}_0+b-Q_y)}\right].$

t(x) is the projective transformation.

Transformation t(x) can be simplified further. First, add its two terms to form a fraction:

$t(x)=\frac{(m{Q}_x{P}_y-Q_y{P}_y+b{Q}_y)x_0+(b{Q}_x{P}_y-b{Q}_y{P}_x)}{m(P_y-Q_y)x_0+(m{P}_x{Q}_y+P_y(b-Q_y))}(6)$

Then, define the coefficients α, β, γ and δ to be the following

$\alpha =m{Q}_x{P}_y-Q_y{P}_y+b{Q}_y,$

$\beta =b{Q}_x{P}_y-b{Q}_y{P}_x,$

$\gamma =m(P_y-Q_y),$

$\delta =m{P}_x{Q}_y+P_y(b-Q_y).$

Substitute these coefficients into equation (6), in order to produce

$t(x)=\frac{\alpha x+\beta }{\gamma x+\delta }$

This is the Möbius transformation or linear fractional transformation.

Inverse transformation

It is clear from the synthetic definition that the inverse transformation is obtained by exchanging points P and Q. This can also be shown analytically. If P ↔ Q, then α → α′, β → β′, γ → γ′, and δ → δ′, where

${\alpha }^{\prime }=m{P}_x{Q}_y-P_y{Q}_y+b{P}_y=\delta ,$

${\beta }^{\prime }=b{P}_x{Q}_y-b{P}_y{Q}_x=-\beta ,$

${\gamma }^{\prime }=m(Q_y-P_y)=-\gamma ,$

${\delta }^{\prime }=m{Q}_x{P}_y+b{Q}_y-Q_y{P}_y=\alpha .$

Therefore if the forwards transformation is

$t(x)=\frac{\alpha x+\beta }{\gamma x+\delta }$

then the transformation t′ obtained by exchanging P and Q (P ↔ Q) is:

$t^{\prime }(x)=\frac{\delta x-\beta }{-\gamma x+\alpha }.$

Then

$t^{\prime }(t(x))=\frac{\delta \left(\frac{\alpha x+\beta }{\gamma x+\delta }\right)-\beta }{-\gamma \left(\frac{\alpha x+\beta }{\gamma x+\delta }\right)+\alpha }$.

Dissolve the fractions in both numerator and denominator of the right side of this last equation:

$t^{\prime }(t(x))=\frac{\alpha \delta x+\beta \delta -\beta \gamma x-\beta \delta }{-\alpha \gamma x-\beta \gamma +\alpha \gamma x+\alpha \delta }$

$=\frac{\alpha \delta x-\beta \gamma x}{\alpha \delta -\beta \gamma }=x$.

Therefore t′(x) = t⁻¹(x): the inverse projective transformation is obtained by exchanging observers P and Q, or by letting α ↔ δ, β → −β, and γ → −γ. This is, by the way, analogous to the procedure for obtaining the inverse of a two-dimensional matrix:

$\left[\begin{array}{cc}\alpha & \beta \\ \gamma & \delta \end{array}\right]\left[\begin{array}{cc}\delta & -\beta \\ -\gamma & \alpha \end{array}\right]=\mathrm{\Delta }\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$

where Δ = α δ − β γ is the determinant.

Identity transformation

Also analogous with matrices is the identity transformation, which is obtained by letting α = 1, β = 0, γ = 0, and δ = 1, so that

$t_I(x)=x.$

Composition of transformations

It remains to show that there is closure in the composition of transformations. One transformation operating on another transformation produces a third transformation. Let the first transformation be t₁ and the second one be t₂:

$t_1(x)=\frac{{\alpha }_{1}x+{\beta }_1}{{\gamma }_{1}x+{\delta }_1},$

$t_2(x)=\frac{{\alpha }_{2}x+{\beta }_2}{{\gamma }_{2}x+{\delta }_2}.$

The composition of these two transformations is

$t_2(t_1(x))=\frac{{\alpha }_2\left(\frac{{\alpha }_{1}x+{\beta }_1}{{\gamma }_{1}x+{\delta }_1}\right)+{\beta }_2}{{\gamma }_2\left(\frac{{\alpha }_{1}x+{\beta }_1}{{\gamma }_{1}x+{\delta }_1}\right)+{\delta }_2}$

$=\frac{{\alpha }_2{\alpha }_{1}x+{\alpha }_2{\beta }_1+{\beta }_2{\gamma }_{1}x+{\beta }_2{\delta }_1}{{\gamma }_2{\alpha }_{1}x+{\gamma }_2{\beta }_1+{\delta }_2{\gamma }_{1}x+{\delta }_2{\delta }_1}$

$=\frac{({\alpha }_2{\alpha }_1+{\beta }_2{\gamma }_1)x+({\alpha }_2{\beta }_1+{\beta }_2{\delta }_1)}{({\gamma }_2{\alpha }_1+{\delta }_2{\gamma }_1)x+({\gamma }_2{\beta }_1+{\delta }_2{\delta }_1)}.$

Define the coefficients α₃, β₃, γ₃ and δ₃ to be equal to

${\alpha }_3={\alpha }_2{\alpha }_1+{\beta }_2{\gamma }_1,$

${\beta }_3={\alpha }_2{\beta }_1+{\beta }_2{\delta }_1,$

${\gamma }_3={\gamma }_2{\alpha }_1+{\delta }_2{\gamma }_1,$

${\delta }_3={\gamma }_2{\beta }_1+{\delta }_2{\delta }_1.$

Substitute these coefficients into $t_2(t_1(x))$ to obtain

$t_2(t_1(x))=\frac{{\alpha }_{3}x+{\beta }_3}{{\gamma }_{3}x+{\delta }_3}.$

Projections operate in a way analogous to matrices. In fact, the composition of transformations can be obtained by multiplying matrices:

$\left[\begin{array}{cc}{\alpha }_2& {\beta }_2\\ {\gamma }_2& {\delta }_2\end{array}\right]\left[\begin{array}{cc}{\alpha }_1& {\beta }_1\\ {\gamma }_1& {\delta }_1\end{array}\right]=\left[\begin{array}{cc}{\alpha }_2{\alpha }_1+{\beta }_2{\gamma }_1& {\alpha }_2{\beta }_1+{\beta }_2{\delta }_1\\ {\gamma }_2{\alpha }_1+{\delta }_2{\gamma }_1& {\gamma }_2{\beta }_1+{\delta }_2{\delta }_1\end{array}\right]=\left[\begin{array}{cc}{\alpha }_3& {\beta }_3\\ {\gamma }_3& {\delta }_3\end{array}\right].$

Since matrices multiply associatively, it follows that composition of projections is also associative.

Projections have: an operation (composition), associativity, an identity, an inverse and closure, so they form a group.

The cross-ratio defined by means of a projection

Let there be a transformation t_s such that t_s(A) = $\mathrm{\infty }$, t_s(B) = 0, t_s(C) = 1. Then the value of t_s(D) is called the cross-ratio of points A, B, C and D, and is denoted as [A, B, C, D]_s:

$[A,B,C,D{]}_s=t_s(D).$

Let

$t_s(x)=\frac{\alpha x+\beta }{\gamma x+\delta },$

then the three conditions for t_s(x) are met when

$t_s(A)=\frac{\alpha A+\beta }{\gamma A+\delta }=\mathrm{\infty },(7)$

$t_s(B)=\frac{\alpha B+\beta }{\gamma B+\delta }=0,(8)$

$t_s(C)=\frac{\alpha C+\beta }{\gamma C+\delta }=1.(9)$

Equation (7) implies that $\gamma A+\delta =0$, therefore $\delta =-\gamma A$. Equation (8) implies that $\alpha B+\beta =0$, so that $\beta =-\alpha B$. Equation (9) becomes

$\frac{\alpha C-\alpha B}{\gamma C-\gamma A}=1,$

which implies

$\gamma =\alpha \frac{C-B}{C-A}.$

Therefore

$t_s(D)=\frac{\alpha D-\alpha B}{\alpha \left(\frac{C-B}{C-A}\right)D-\gamma A}=\frac{\alpha (D-B)}{\alpha \left(\frac{C-B}{C-A}\right)D-\alpha \left(\frac{C-B}{C-A}\right)A}$

$=\frac{D-B}{C-B}\frac{C-A}{D-A}=\frac{A-C}{A-D}\frac{B-D}{B-C}.(10)$

In equation (10), it is seen that t_s(D) does not depend on the coefficients of the projection t_s. It only depends on the positions of the points on the "subjective" projective line. This means that the cross-ratio depends only on the relative distances among four collinear points, and not on the projective transformation which was used to obtain (or define) the cross-ratio. The cross ratio is therefore

$[A,B,C,D]=\frac{A-C}{A-D}\frac{B-D}{B-C}.(11)$

Conservation of cross-ratio

Transformations on the projective line preserve cross ratio. This will now be proven. Let there be four (collinear) points A, B, C, D. Their cross-ratio is given by equation (11). Let S(x) be a projective transformation:

$S(x)=\frac{\alpha x+\beta }{\gamma x+\delta }$

where $\alpha \delta \ne \beta \gamma$. Then

$[S(A)S(B)S(C)S(D)]=\frac{\frac{\alpha A+\beta }{\gamma A+\delta }-\frac{\alpha C+\beta }{\gamma C+\delta }}{\frac{\alpha A+\beta }{\gamma A+\delta }-\frac{\alpha D+\beta }{\gamma D+\delta }}\cdot \frac{\frac{\alpha B+\beta }{\gamma B+\delta }-\frac{\alpha D+\beta }{\gamma D+\delta }}{\frac{\alpha B+\beta }{\gamma B+\delta }-\frac{\alpha C+\beta }{\gamma C+\delta }}$

$=\frac{[(\alpha A+\beta )(\gamma C+\delta )-(\alpha C+\beta )(\gamma A+\delta )][(\alpha B+\beta )(\gamma D+\delta )-(\alpha D+\beta )(\gamma B+\delta )]}{[(\alpha A+\beta )(\gamma D+\delta )-(\alpha D+\beta )(\gamma A+\delta )][(\alpha B+\beta )(\gamma C+\delta )-(\alpha C+\beta )(\gamma B+\delta )]}$

$=\frac{[\alpha A\delta +\beta \gamma C-\alpha C\delta -\beta \gamma A][\alpha B\delta +\beta \gamma D-\alpha D\delta -\beta \gamma B]}{[\alpha A\delta +\beta \gamma D-\alpha D\delta -\beta \gamma A][\alpha B\delta +\beta \gamma C-\alpha C\delta -\beta \gamma B]}$

$=\frac{[\alpha \delta (A-C)+\beta \gamma (C-A)][\alpha \delta (B-D)+\beta \gamma (D-B)]}{[\alpha \delta (A-D)+\beta \gamma (D-A)][\alpha \delta (B-C)+\beta \gamma (C-B)]}$

$=\frac{(\alpha \delta -\beta \gamma )(A-C)(\alpha \delta -\beta \gamma )(B-D)}{(\alpha \delta -\beta \gamma )(A-D)(\alpha \delta -\beta \gamma )(B-C)}$

$=\frac{A-C}{A-D}\cdot \frac{B-D}{B-C}$

Therefore [S(A) S(B) S(C) S(D)] = [A B C D], Q.E.D.