Plücker matrix explained

The Plücker matrix is a special skew-symmetric 4 × 4 matrix, which characterizes a straight line in projective space. The matrix is defined by 6 Plücker coordinates with 4 degrees of freedom. It is named after the German mathematician Julius Plücker.

Definition

A straight line in space is defined by two distinct points

A = \left(A _{0,} A _{1,} A _{2,} A ^{\top} _{3 \right)} \in R l{P} ^{3} and B = \left(B _{0,} B _{1,} B _{2,} B ^{\top} _{3 \right)} \in R l{P} ^{3} in homogeneous coordinates of the projective space. Its Plücker matrix is:[ L ] _{x} \propto A B ^{\top} - B A ^{\top} = \left(\begin{array}{cccc} 0 & - L _{01} & - L _{02} & - L _{03} \\ L _{01} & 0 & - L _{12} & - L _{13} \\ L _{02} & L _{12} & 0 & - L _{23} \\ L _{03} & L _{13} & L _{23} & 0 \end{array}\right) L \propto(L _{01} , L _{02} , L _{03} , L _{12} , L _{13} , L _{23} ) ^{\top}

with

L _{ij} = A _{iB} _{j} - B _{iA} _{j.}

Plücker coordinates fulfill the Grassmann–Plücker relations

L _{01} L _{23} - L _{02} L _{13} + L _{03} L _{12} = 0

and are defined up to scale. A Plücker matrix has only rank 2 and four degrees of freedom (just like lines in

R ^{3} ). They are independent of a particular choice of the points A and B and can be seen as a generalization of the line equation i.e. of the cross product for both the intersection (meet) of two lines, as well as the joining line of two points in the projective plane.

Properties

The Plücker matrix allows us to express the following geometric operations as matrix-vector product:

• Plane contains line:

0 = [ L ] _{x} E X = [ L ] _{x} E is the point of intersection of the line L and the plane E ('Meet')

• Point lies on line:

0 = [\tilde{ L }]_\mathbfE = [\tilde{ L }]_\mathbf is the common plane E , which contains both the point X and the line L ('Join').

• Direction of a line:

[ L ] _{x} \pi ^{infty} = [ L ] _{x} ( 0 , 0 , 0 , 1 ) ^{\top} = \left( - L _{03} , - L _{13} , - L _{23} , 0 \right) ^{\top} (Note: The latter can be interpreted as a plane orthogonal to the line passing through the coordinate origin)

• Closest point to the origin

X _{0} \cong [ L ] _{x} [ L ] _{x} \pi ^{infty} .

Uniqueness

Two arbitrary distinct points on the line can be written as a linear combination of

A and B :A ^{\prime} \propto A \alpha + B \beta and B ^{\prime} \propto A \gamma + B \delta.

Their Plücker matrix is thus:

\begin{align} {[} L ^{\prime{]}} _{x} & = A ^{\primeB} ^{\prime} - B ^{\prime} A ^{\prime} \\[ 6 pt] & = ( A \alpha + B \beta)( A \gamma + B \delta) ^{\top} - ( A \gamma + B \delta)( A \alpha + B \beta) ^{\top} \\[ 6 pt] & = \underbrace{(\alpha\delta - \beta\gamma)} _{λ[L]} _{x , \end{align}}

up to scale identical to

[ L ] _{x} .

Intersection with a plane

Let

E = \left(E _{0} , E _{1} , E _{2} , E _{3} \right) ^{\top} \in R l{P} ^{3} denote the plane with the equationE _{0} x + E _{1} y + E _{2} z + E _{3} = 0 .

which does not contain the line

L . Then, the matrix-vector product with the Plücker matrix describes a pointX = [ L ] _{x} E = A \underset{\alpha}{\underbrace{ B ^{\top} E }} - \mathbf\underset = \mathbf\alpha + \mathbf\beta,

which lies on the line

L because it is a linear combination of A and B . X is also contained in the plane E E ^{\top} X = E ^{\top} [ L ] _{x} E = \underset{\alpha}{\underbrace{ E ^{\top} A }}\underset - \underset\underset = 0,

and must therefore be their point of intersection.

In addition, the product of the Plücker matrix with a plane is the zero-vector, exactly if the line

L is contained entirely in the plane:\alpha = \beta = 0 \iff E contains L .

Dual Plücker matrix

In projective three-space, both points and planes have the same representation as 4-vectors and the algebraic description of their geometric relationship (point lies on plane) is symmetric. By interchanging the terms plane and point in a theorem, one obtains a dual theorem which is also true.

In case of the Plücker matrix, there exists a dual representation of the line in space as the intersection of two planes:

E = \left(E _{0,} E _{1,} E _{2,} E ^{\top} _{3 \right)} \in R l{P} ^{3}

and

F = \left(F _{0,} F _{1,} F _{2,} F ^{\top} _{3 \right)} \in R l{P} ^{3}

in homogeneous coordinates of projective space. Their Plücker matrix is:

\left[\tilde{ L }\right]_ = \mathbf\mathbf^ - \mathbf\mathbf^

and

G = \left[\tilde{ L }\right]_\mathbf

describes the plane

G which contains both the point X and the line L .

Relationship between primal and dual Plücker matrices

As the vector

X = [ L ] _{x} E , with an arbitrary plane E , is either the zero-vector or a point on the line, it follows:\forall E \in R l{P} ^{3} : X = [ L ] _{x} E lieson L \iff \left[\tilde{ L }\right]_\mathbf = \mathbf.

Thus:

\left([\tilde{ L }]_[\mathbf{L}]_\right)^ = [\mathbf{L}]_\left[\tilde{\mathbf{L}}\right]_ = \mathbf \in \mathbb^.

The following product fulfills these properties:

\begin{align} &\left(\begin{array}{cccc} 0 & L _{23} & - L _{13} & L _{12} \\ - L _{23} & 0 & L _{03} & - L _{02} \\ L _{13} & - L _{03} & 0 & L _{01} \\ - L _{12} & L _{02} & - L _{01} & 0 \end{array}\right) \left(\begin{array}{cccc} 0 & - L _{01} & - L _{02} & - L _{03} \\ L _{01} & 0 & - L _{12} & - L _{13} \\ L _{02} & L _{12} & 0 & - L _{23} \\ L _{03} & L _{13} & L _{23} & 0 \end{array}\right) \\[ 10 pt] = {} &\left(L _{01} L _{23} - L _{02} L _{13} + L _{03} L _{12} \right) ⋅ \left(\begin{array}{cccc} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{array}\right) = 0 , \end{align}

due to the Grassmann–Plücker relation. With the uniqueness of Plücker matrices up to scalar multiples, for the primal Plücker coordinates

L = \left(L _{01} ,L _{02} ,L _{03} ,L _{12} ,L _{13} ,L _{23} \right) ^{\top}

we obtain the following dual Plücker coordinates:

\tilde{ L } = \left(L_,\,-L_,\,L_,\,L_,\,-L_,\,L_\right)^.

In the projective plane

The 'join' of two points in the projective plane is the operation of connecting two points with a straight line. Its line equation can be computed using the cross product:

l \propto a x b = \left(\begin{array}{c} a _{1} b _{2} - b _{1} a _{2} \\ b _{0} a _{2} - a _{0} b _{2} \\ a _{0} b _{1} - a _{1} b _{0} \end{array}\right) = \left(\begin{array}{c} l _{0} \\ l _{1} \\ l _{2} \end{array}\right).

Dually, one can express the 'meet', or intersection of two straight lines by the cross-product:

x \propto l x m

The relationship to Plücker matrices becomes evident, if one writes the cross product as a matrix-vector product with a skew-symmetric matrix:

[ l ] _{x} = a b ^{\top} - b a ^{\top} = \left(\begin{array}{ccc} 0 & l _{2} & - l _{1} \\ - l _{2} & 0 & l _{0} \\ l _{1} & - l _{0} & 0 \end{array}\right)

and analogously

[ x ] _{x} = l m ^{\top} - m l ^{\top}

Geometric interpretation

Let

d = \left( - L _{03} , - L _{13} , - L _{23} \right) ^{\top} and m = \left(L _{12} , - L _{02} , L _{01} \right) ^{\top} , then we can write[ L ] _{x} = \left(\begin{array}{cc} [ m ] _{x} & d \\ - d & 0 \end{array}\right)

and

[\tilde{ L }]_ = \left(\begin [-\mathbf{d}]_ & \mathbf\\ -\mathbf & 0 \end\right),

where

d is the displacement and m is the moment of the line, compare the geometric intuition of Plücker coordinates.

References

• Book: Richter-Gebert, Jürgen . Springer Science & Business Media . Perspectives on Projective Geometry: A Guided Tour Through Real and Complex Projective Geometry . 2011 . 978-3-642-17286-1 .
• Book: Jorge Stolfi . Jorge Stolfi . Oriented Projective Geometry: A Framework for Geometric Computations . . 1991 . 978-1483247045 .
  From original Stanford University 1988 Ph.D. dissertation, Primitives for Computational Geometry, available as http://www.hpl.hp.com/techreports/Compaq-DEC/SRC-RR-36.pdf.
• Blinn. James F.. Jim Blinn . ACM SIGGRAPH Computer Graphics . 11 . 2 . A homogeneous formulation for lines in 3 space . 0097-8930 . 237–241 . Aug 1977 . 10.1145/965141.563900 .