Laguerre formula

Short description: Mathematical formula

The Laguerre formula (named after Edmond Laguerre) provides the acute angle [math]\displaystyle{ \phi }[/math] between two proper real lines, as follows:

[math]\displaystyle{ \phi=|\frac {1}{2i}\operatorname{Log} \operatorname{Cr}(I_1,I_2,P_1,P_2)| }[/math]

where:

• [math]\displaystyle{ \operatorname{Log} }[/math] is the principal value of the complex logarithm
• [math]\displaystyle{ \operatorname{Cr} }[/math] is the cross-ratio of four collinear points
• [math]\displaystyle{ P_1 }[/math] and [math]\displaystyle{ P_2 }[/math] are the points at infinity of the lines
• [math]\displaystyle{ I_1 }[/math] and [math]\displaystyle{ I_2 }[/math] are the intersections of the absolute conic, having equations [math]\displaystyle{ x_0=x_1^2+x_2^2+x_3^2=0 }[/math], with the line joining [math]\displaystyle{ P_1 }[/math] and [math]\displaystyle{ P_2 }[/math].

The expression between vertical bars is a real number.

Laguerre formula can be useful in computer vision, since the absolute conic has an image on the retinal plane which is invariant under camera displacements, and the cross ratio of four collinear points is the same for their images on the retinal plane.

Derivation

It may be assumed that the lines go through the origin. Any isometry leaves the absolute conic invariant, this allows to take as the first line the x axis and the second line lying in the plane z=0. The homogeneous coordinates of the above four points are

[math]\displaystyle{ (0,1,i,0),\ (0,1,-i,0),\ (0,1,0,0),\ (0,\cos\phi,\pm\sin\phi,0), }[/math]

respectively. Their nonhomogeneous coordinates on the infinity line of the plane z=0 are [math]\displaystyle{ i }[/math], [math]\displaystyle{ -i }[/math], 0, [math]\displaystyle{ \pm\sin\phi/\cos\phi }[/math]. (Exchanging [math]\displaystyle{ I_1 }[/math] and [math]\displaystyle{ I_2 }[/math] changes the cross ratio into its inverse, so the formula for [math]\displaystyle{ \phi }[/math] gives the same result.) Now from the formula of the cross ratio we have [math]\displaystyle{ \operatorname {Cr}(I_1,I_2,P_1,P_2)=-\frac{-i\cos\phi\pm\sin\phi}{i\cos\phi\pm\sin\phi}=e^{\pm 2i\phi}. }[/math]

References

• O. Faugeras. Three-dimensional computer vision. MIT Press, Cambridge, London, 1999.