Projective Geometry: Homography, Cross-ratio, Harmonic Conjugates, Polarity
[url=https://en.wikipedia.org/wiki/Projective_geometry]Projective geometry[/url] is the geometry of what we see in real life where paralel lines meet at the horizon and circles are seen as ellipses. Those effects are the result of cone-shaped projection: light rays from objects come to the eye, or in reverse, light rays from a flashlight shine objects and cast shadow on a screen.[br][br]Here we use circled squares as objects of projection between a "ground" plane and a slanted "screen": (see [url=https://www.geogebra.org/m/h3txa4nw]this 3D applet[/url])
We can see the object on the ground is projected into an image on the screen as well as the object on the screen is projected into an image on the ground. Thus each plane can be both [i]object plane[/i] and [i]image plane[/i]. As image planes, the screen contains a [i]vanishing line[/i] of the ground (ground horizon) and the ground contains a [i]vanishing line[/i] of the screen (screen horizon), which are yellow high lighted.[br][br]Let's turn aside to see clearer, reducing from 3D to 2D:[br]
Here on 2D space, we have object/image lines containing points (yellow high lighted) instead of object/image planes, and vanishing points (red marked) instead of vanishing lines. The object points and the corresponding image points are [i]in perspective[/i], and the mapping between them is called [b]perspectivity[/b], the simplest projective transformation. In reverse, for any two lines [i]h[/i] and [i]ℓ[/i] (see applet below), a triple of points ([i]N, P, E[/i]) on [i]h[/i] is [i]in perspective[/i] with a triple of points ([i]N', P', E'[/i]) on [i]ℓ[/i] only if the three projection lines [i]NN', PP', EE'[/i] are [i]concurrent[/i], i.e. they meet at a single point [i]V[/i]. In other words, any [i]two [b]couples [/b]of points ([i]N, P[/i]) and [i]([i]N', P'[/i])[/i] determine a perspectivity from V = NN' [i] ∩ PP'[/i][/i], which maps any other point, like [i]E ∈ h[/i] to [i]E'[/i] = [i]VE ∩ ℓ[/i] and [i]I[i] ∈ h[/i][/i] to [i]I'[/i] = [i]VI ∩ ℓ[/i]. [br][br]But how about two [i]arbitrary triples [/i]of colinear points?[br][br][i]Two arbitrary [b]triples [/b]of colinear points determine a [b]projectivity[/b][/i] or [url=https://en.wikipedia.org/wiki/Homography]homography[/url] which is the general form of projective transformation including both single perspectivity and [i]composition of perspectivities[/i]. The projectivity is not very much different from the perspectivity in that any projectivity can be decomposed into at most [i]two perspectivities[/i], and a self projectivity (identical mapping) can be decomposed into three distint perspectivities. That's the content of [url=https://youtu.be/3BRtaHlUItQ]the fundamental theorem of projective geometry[/url]:[br][br][math]\forall Obj,\forall Img:[/math][br][math]\left[\exists p\in Proj:Obj\mapsto_pImg\right]\Longleftrightarrow\left[\exists p_1,p_2\in Pers:Obj\mapsto_{p_1}mid\mapsto_{p_2}Img\right][/math]
Projective transformation does [i]not [/i]preserve distance [i]nor [/i]ratio between distances, but it preserves the [url=https://en.wikipedia.org/wiki/Cross-ratio]cross-ratio[/url] (the ratio of ratio) between four colinear points. So, cross-ratio is the crucial measure defining the structure of projective spaces. [br][br]The [i][b]cross-ratio[/b] of [b]quadruple [/b]of colinear points [/i]([i]N[/i], [i]P,[/i] [i]E[/i], [i]I[/i]) is defined using [i]signed [/i]segments (i.e. [i]NP[/i] = -[i]PN[/i]):[br][math]\left(N,P;E,I\right)=\left(E,I;N,P\right)=\frac{NE\cdot PI}{NI\cdot PE}=\frac{NE}{NI}:\frac{PE}{PI}=\frac{NE}{PE}:\frac{NI}{PI}[/math]
Here, when we perspectively project ([i]N, P, E, I[/i]) from [i]object line[/i] [i]h[/i] onto [i]image line[/i] [math]\ell[/math], we have the [i]perspective image[/i] ([i]N', P', E', I'[/i]) with the same cross-ratio ([i]N', P'; E', I'[/i]) = ([i]N, P; E, I[/i]), for [b]any [/b][i]center of perspectivity[/i] V (view point, eye or light source) and [b]any [/b]image line [math]\ell[/math]. That means if two images ([i]N, P, E, I[/i]) and ([i]N', P', E', I'[/i]) are [i]in perspective[/i], then their cross-ratios are the same. [br][br][math]\forall V,\forall\ell:\left[\left(N,P,E,I\right)\mapsto_{V,\ell}\left(N',P',E',I'\right)\Longrightarrow\left(N,P;E,I\right)=\left(N',P';E',I'\right)\right][/math][br][br]For general projectivity (homography), we use [url=https://en.wikipedia.org/wiki/Pappus%27s_hexagon_theorem]Pappus's hexagon theorem[/url] to construct the intermediate image line [i]g[/i] to decompose the projectivity [math]h\mapsto_{Proj}\ell[/math] into two perspectivities [math]h\mapsto_{N',g}g\mapsto_{N,g}\ell[/math].
Here, we have the equivalence between projectivity and cross-ratio equality:[br][br][math]\forall\left(N,P,E,I\right),\forall\left(N',P',E',I'\right):[/math][br][math]\left(N,P;E,I\right)=\left(N',P';E',I'\right)\Longleftrightarrow\exists p\in Proj:\left(N,P,E,I\right)\mapsto_p\left(N',P',E',I'\right)[/math][br][math]\Longleftrightarrow\exists g:\left(N,P,E,I\right)\mapsto_{N',g}\left(N'',P'',E'',I''\right)\mapsto_{N,g}\left(N',P',E',I'\right)[/math][br]
When the cross-ratio ([i]N[/i], [i]P[/i]; [i]E[/i], [i]I[/i]) = −1, they are called [b]harmonic conjugates[/b] ([i]harmonic points[/i]) where [i]E[/i] and [i]I[/i] [url=https://en.wikipedia.org/wiki/Projective_harmonic_conjugate]harmonically conjugate[/url] with each other with respect to [i]N[/i] and [i]P[/i]. And because ([i]E[/i], [i]I[/i]; [i]N[/i], [i]P[/i]) = ([i]N[/i], [i]P[/i]; [i]E[/i], [i]I[/i]), [i]N[/i] and [i]P[/i] also harmonically conjugate with each other with respect to [i]E[/i] and [i]I[/i]. Because cross-ratio is invariant under projectivity, for any view point V not lying on [i]h[/i], we have four [i]harmonic rays[/i] or [i]harmonic straights[/i] ([i]VN, VP, VE, VI[/i]) which project onto any line a quadruple of harmonic points like (N', P', E', I').
The quadruple of harmonic points ([i]N, P, E, I[/i]) determines the [i]basis of the [b]perspective metric[/b][/i]: [i]N[/i], [i]P[/i] are unit points (-1, +1) and [i]E, I[/i] are central points (0, ∞). This metric is reciprocally linear, i.e. the reciprocal of distance to the vanishing point [i]E[/i] increases linearly: [math]\frac{1}{EN},\frac{1}{EI},\frac{1}{EP},\frac{1}{EP_2},\frac{1}{EP_3},\frac{1}{EP_4},\frac{1}{EP_5},...,\infty[/math][br] [math]\frac{1}{EP_{k+1}}-\frac{1}{EP_k}=const\left(c\right)[/math], [b]([i]P[sub]k[/sub], P; E, [b]∞[/b][/i]) = c/k[/b], (?) where [i]N = P[sub]-1[/sub][/i], [i]I = P[sub]0[/sub][/i] and [i]P = P[sub]1[/sub][/i]. [br]Because ([i]N, P; E, I[/i]) = -1, the actual basis reduces to only [i]3 free points[/i] ([i]P, E, I[/i]), where the 4th point [i]N[/i] can be calculated from those 3 points. The basis triple ([i]P, E, I[/i]) of perspective metric corresponds to (1, ∞, 0) in linear metric, where any point [i]P[sub]-x[/sub][/i] has cross-ratio ([i]P[sub]-x[/sub], P; E, I[/i]) = (-x, 1; ∞, 0) = -x exactly ...?... [b]x = (x, 1; 0, [b]∞[/b])[/b] for arbitrary x ranging from -∞ to +∞.[br]Note that in this page, only odd points [i]P[sub]2k-1[/sub][/i] are shown, because [i]N[/i] is used in the place of [i]I[/i], i.e. (P[sub]2k-1[/sub], P; E, N) = k, and [i]I[/i] is ignored in the perspective metric and exponential metric ([i]I[/i] will takes different role in growth metric).[br]That [i]reciprocal linearity[/i] is how we see the array of light poles on the street. And that also means when we project it onto a screen [i]ℓ[/i] parallel to the ray [i]VE[/i] from view point [i]V[/i] to vanishing point [i]E[/i] (click "Make Para∥el", chose "ℓ ∥ VE" in the applet above), the metric becomes [i]linear [/i]as normal.
To see the whole image of perspective metric, we wrap the line [i]ℓ[/i]([i]N', P', E', I'[/i]) around to form a [url=https://en.wikipedia.org/wiki/Riemann_sphere]Riemann circle[/url] ([i]V, N[/i][i]', P', E[sub]c[/sub], I[sub]c[/sub][/i][sub][/sub]), via [url=https://en.wikipedia.org/wiki/Stereographic_projection]stereographic projection[/url]. There we see clearly the 4 points: "negative unit" [i]N[/i], "positive unit" [i]P[/i], "internal center" [i]I[/i], "external center" [i]E[/i]. [i]I[/i] & [i]E[/i] also have other names: absolutes, limits, boundaries, balance points, fixed points, attractors, extremities.
Perspective metric & linear metric (cross-ratio = -1)
In non-harmonic cases where cross-ratio ≠ -1, we have two other kinds of metrics: [url=http://www.budworkshop.co.uk/bw_algebra_of_measure.html]exponential metric and growth metric[/url] generalizing the linear metric and projective metric in the harmonic case. (In the applet above, uncheck "Harmonic" and choose "Exponential" or "Growth" in the Metric dropdown list). The growth metric is the most general one, which has [url=https://en.wikipedia.org/wiki/Cayley%E2%80%93Klein_metric#Special_relativity]application in Einstein's Relativity[/url].[br][br]The projective & linear metric is just the growth & exponential metric where the two limits (I, E) coincide
Growth metric & exponential metric (cross-ratio = -1.12)
In the following applet, mimicking the case of linear perspective metric, we define a linear metric m() on [i]h[/i] using a constant [i]c[/i] for the unit pair ([i]N, P[/i]) such that m([i]N[/i])=-[i]c[/i] and m([i]P[/i])=[i]c[/i]. Then the measurements of the central pair ([i]E, I[/i]), [i]d[/i]:=m([i]I[/i]) and [i]e[/i]:=m([i]E[/i]), are harmonic conjugates in the sense of [url=https://en.wikipedia.org/wiki/Inversive_geometry]inverse geometry[/url]: [math]d\cdot e=c^2[/math]. That's because:[br][math]\left(N,P;E,I\right)=\frac{NE\cdot PI}{NI\cdot PE}=\frac{\left(c+e\right)\left(d-c\right)}{\left(d+c\right)\left(e-c\right)}=\frac{de-ce+dc-c^2}{de+ce-dc-c^2}=-1[/math][br][math]\Longleftrightarrow2de-2c^2=0\Longleftrightarrow de=c^2[/math][br][br]Moreover, we can rewrite the harmonic condition ([i]N[/i], [i]P[/i]; [i]E[/i], [i]I[/i]) = −1 as a [i]zero sum[/i]:[br][math]\frac{NE}{PE}:\frac{NI}{PI}=-1\Longleftrightarrow\frac{NE}{PE}=-\frac{NI}{PI}\Longleftrightarrow\frac{NI}{PI}+\frac{NE}{PE}=0[/math]
Here, we construct a [b]quadrilateral[/b] [i]VTRS[/i] (in blue) using an arbitrary point [i]T[/i] on [i]VN[/i], and two intersection points [i]R[/i] = [i]VI[/i] ∩ [i]TP[/i], [i]S[/i] = [i]VP[/i] ∩ [i]TE[/i]. Then R, S are colinear with N. And in reverse, for [i]any quadrilateral RSTV[/i] (convex, concave, complex) with three pairs of opposite sides intersect at [i]N[/i] = [i]VT[/i] ∩ [i]SR[/i], [i]P[/i] = [i]VS[/i] ∩ [i]TR[/i], [i]Pole[/i] = [i]VR[/i] ∩ [i]TS[/i], and the line [i]NP[/i] intersects with the other two opposite sides at [i]E[/i] = [i]NP[/i] ∩ [i]VR[/i], [i]I[/i] = [i]NP[/i] ∩ [i]TS[/i], there are three harmonic quadruples ([i]N[/i], [i]P[/i]; [i]E[/i], [i]I[/i]), ([i]V[/i], [i]R[/i]; [i]I[/i], [i]Pole[/i]) and ([i]T[/i], [i]S[/i]; [i]E[/i], [i]Pole[/i]). Then using a variant of Pascal's theorem, we construct the [i][b]conic[/b] (RSTVU[sub]c[/sub]) circumscribing the quadrilateral RSTV[/i] (check "χExtras" to see the point [i]U[sub]c[/sub][/i]).[br][br]Then we can always construct a perspectivity to project that "conic circumscribed quadrilateral" into a "circled square", as shown in [url=https://www.geogebra.org/m/s6aqr24j]this 3D applet[/url]. That means [i]any quadrilateral is a projective image of a [b]square[/b][/i] as shown in [url=https://youtu.be/a-OYdVGNdJI]this video[/url], [url=https://people.southwestern.edu/~futamurf/Image-of-A-Square-2016-03-25.pdf]this article[/url], and [i]any [b]conic [/b]is a projective image of a [b]circle[/b][/i] as shown in the article "[url=https://www.geogebra.org/m/nknaryzz]Shine a light on Conic sections[/url]". So, the line [i]h[/i] here is the vanishing line (horizon), and N, P are two vanishing points of the sides of that "square" (square in the object plane).[br][br]Moreover, the intersect of two diagonals is the pole and the vanishing line is the [url=https://en.wikipedia.org/wiki/Pole_and_polar]polar[/url] of the circumscribed conic.[br]
Project an ellipse into a full circle
Project a hyperbola into a circle (only the upper part of the circle is visible infront of the eye)