Editing Plücker coordinates

{{short description|Method of assigning coordinates to every line in projective 3-space}}
{{About|the classic case of lines in projective 3-space|general Plücker coordinates|Plücker embedding}}

{{No footnotes|date=February 2011}}

In [[geometry]], '''Plücker coordinates''', introduced by [[Julius Plücker]] in the 19th century, are a way to assign six [[homogeneous coordinates]] to each [[line (mathematics)|line]] in [[projective space|projective 3-space]], {{tmath|\mathbb P^3}}. Because they satisfy a quadratic constraint, they establish a [[one-to-one correspondence]] between the 4-dimensional space of lines in {{tmath|\mathbb P^3}} and points on a [[quadric]] in {{tmath|\mathbb P^5}} (projective 5-space). A predecessor and special case of [[Grassmann coordinates]] (which describe {{mvar|k}}-dimensional linear subspaces, or [[Flat (geometry)|''flats'']], in an {{mvar|n}}-dimensional [[Euclidean space]]), Plücker coordinates arise naturally in [[geometric algebra]]. They have proved useful for [[computer graphics]], and also can be extended to coordinates for the [[screw theory|screws and wrenches]] in the theory of [[kinematics]] used for [[robot control]].

== Geometric intuition ==
[[File:Plücker line coordinate geometry.png|thumb|right|Displacement ''d'' (yellow arrow) and moment ''m'' (green arrow) of two points ''x'',''y'' on a line (in red)]]

A line {{mvar|L}} in 3-dimensional [[Euclidean space]] is determined by two distinct points that it contains, or by two distinct planes that contain it (a [[plane-plane intersection]]). Consider the first case, with points <math>x=(x_1,x_2,x_3)</math> and  <math>y=(y_1,y_2,y_3).</math> The [[vector displacement]] from {{mvar|x}} to {{mvar|y}} is nonzero because the points are distinct, and represents the ''[[direction (geometry)|direction]]'' of the line. That is, every displacement between points on the line {{mvar|L}} is a [[scalar multiplication|scalar multiple]] of {{math|1=''d'' = ''y'' − ''x''}}. If a physical particle of unit mass were to move from {{mvar|x}} to {{mvar|y}}, it would have a [[moment (physics)|moment]] about the origin of the coordinate system. The geometric equivalent to this moment is a vector whose direction is perpendicular to the plane containing the line {{mvar|L}} and the origin, and whose length equals twice the area of the triangle formed by the displacement and the origin. Treating the points as displacements from the origin, the moment is {{math|1='''m''' = '''x''' × '''y'''}}, where "×" denotes the vector [[cross product]]. For a fixed line, {{mvar|L}}, the area of the triangle is proportional to the length of the segment between {{mvar|x}} and {{mvar|y}}, considered as the base of the triangle; it is not changed by sliding the base along the line, parallel to itself. By definition the moment vector is perpendicular to every displacement along the line, so {{math|1='''d''' ⋅ '''m''' = 0}}, where "⋅" denotes the vector [[dot product]].

Although neither direction {{mvar|d}} nor moment {{mvar|m}} alone is sufficient to determine the line {{mvar|L}}, together the pair does so uniquely, up to a common (nonzero) scalar multiple which depends on the distance between {{mvar|x}} and {{mvar|y}}. That is, the coordinates

: <math>(\mathbf d : \mathbf m ) = (d_1:d_2:d_3\ :\ m_1:m_2:m_3)</math>

may be considered [[homogeneous coordinates]] for {{mvar|L}}, in the sense that all pairs {{math|(λ'''d''' : λ'''m''')}}, for {{math|λ ≠ 0}}, can be produced by points on {{mvar|L}} and only {{mvar|L}}, and any such pair determines a unique line so long as {{math|'''d'''}} is not zero and {{math|1='''d''' ⋅ '''m''' = 0}}. Furthermore, this approach extends to include [[point at infinity|points]], [[line at infinity|lines]], and a [[plane at infinity|plane "at infinity"]], in the sense of [[projective geometry]]. In addition a point <math>x</math> lies on the line {{mvar|L}} if and only if <math>x \times d = m</math>.

: '''Example.''' Let {{math|1='''x''' = (2, 3, 7)}} and {{math|1='''y''' = (2, 1, 0)}}. Then {{math|1=('''d''' : '''m''') = (0 : −2 : −7 : −7 : 14 : −4)}}.

Alternatively, let the equations for points {{math|'''x'''}} of two distinct planes containing {{mvar|L}} be

: <math>\begin{align}
0 &= a + \mathbf a \cdot \mathbf x, \\ 
0 &= b + \mathbf b \cdot \mathbf x.
\end{align}</math>

Then their respective planes are perpendicular to vectors {{math|'''a'''}} and {{math|'''b'''}}, and the direction of {{mvar|L}} must be perpendicular to both. Hence we may set {{math|1='''d''' = '''a''' × '''b'''}}, which is nonzero because {{math|'''a''', '''b'''}} are neither zero nor parallel (the planes being distinct and intersecting). If point {{math|'''x'''}} satisfies both plane equations, then it also satisfies the linear combination

:<math>\begin{align}
0 &= a (b + \mathbf b \cdot \mathbf x) - b(a+ \mathbf a \cdot \mathbf x) \\
  &= (a \mathbf b - b \mathbf a) \cdot \mathbf x
\end{align}</math>

That is,

:<math>\mathbf m = a \mathbf b - b \mathbf a</math>

is a vector perpendicular to displacements to points on {{mvar|L}} from the origin; it is, in fact, a moment consistent with the {{math|'''d'''}} previously defined from {{math|'''a'''}} and {{math|'''b'''}}.

{{Collapse top|title=Proof of geometric formulation}}

''Proof 1'': Need to show that 
:<math>\mathbf m = a \mathbf b - b \mathbf a = \mathbf r \times \mathbf d = \mathbf r \times (\mathbf a \times \mathbf b).</math><sup>''what is "'''r'''"?''</sup>

[[Without loss of generality]], let

:<math>\mathbf a \cdot \mathbf a = \mathbf b \cdot \mathbf b = 1.</math>

[[File:BaiduShurufa 2018-9-29 20-38-43.png|thumb|Plane orthogonal to line {{mvar|L}} and including the origin.]]
Point {{mvar|B}} is the origin. Line {{mvar|L}} passes through point {{mvar|D}} and is orthogonal to the plane of the picture. The two planes pass through {{mvar|CD}} and {{mvar|DE}} and are both orthogonal to the plane of the picture. Points {{mvar|C}} and {{mvar|E}} are the closest points on those planes to the origin {{mvar|B}}, therefore angles {{math|∠ ''BCD''}} and {{math| ∠ ''BED''}} are right angles and so the points {{mvar|B, C, D, E}} lie on a circle (due to a corollary of [[Thales's theorem]]). {{mvar|BD}} is the diameter of that circle.

: <math>\begin{align}
& \mathbf a := \frac{BE}{||BE||}, \quad  \mathbf b := \frac{BC}{||BC||}, \quad \mathbf r := BD; \\[4pt]
& - \! a = ||BE|| = ||BF||, \quad -b =  ||BC|| = ||BG||; \\[4pt]
& \mathbf m = a \mathbf b - b \mathbf a = FG \\[4pt]
& || \mathbf d || = || \mathbf a \times \mathbf b || = \sin\angle FBG
\end{align}</math>

Angle {{math|∠ ''BHF''}} is a right angle due to the following argument. Let {{math|1=ε := ∠ ''BEC''}}. Since {{math|△ ''BEC'' ≅ △ ''BFG''}} (by side-angle-side congruence), then {{math|1= ∠ ''BFG'' = ε}}. Since {{math|1=∠ ''BEC'' + ∠ ''CED'' = 90°}}, let {{math|1=ε' := 90° − ε = ∠ ''CED''}}. By the [[inscribed angle theorem]], {{math|1=∠ ''DEC'' = ∠ ''DBC''}}, so {{math|1= ∠ ''DBC'' = ε'}}. {{math|1=∠ ''HBF'' + ∠ ''BFH'' + ∠ ''FHB'' = 180°}}; {{math|1=ε' + ε + ∠ ''FHB'' = 180°}}, {{math|1=ε + ε' = 90°}}; therefore, {{math|1= ∠ ''FHB'' = 90°}}. Then {{math|∠ ''DHF''}} must be a right angle as well.

Angles {{math|∠ ''DCF'', ∠ ''DHF''}} are right angles, so the four points {{mvar|C, D, H, F}} lie on a circle, and (by the [[intersecting secants theorem]])

:<math>||BF|| \, ||BC|| = ||BH|| \, ||BD||</math>
that is, 
:<math>\begin{align}
&ab \sin\angle FBG = ||BH|| \, || \mathbf r || \sin\angle FBG , \\[4pt]
& 2 \, \text{Area}_{\triangle BFG} = ab \sin\angle FBG = ||BH|| \, ||FG|| = ||BH|| \, || \mathbf r || \sin\angle FBG, \\[4pt]
& || \mathbf m || = ||FG|| = || \mathbf r || \sin\angle FBG = || \mathbf r || \, || \mathbf d ||, \\[4pt]
& \mathbf m = \mathbf r \times \mathbf d. \blacksquare
\end{align}</math>

''Proof 2'':

Let

:<math>\mathbf a \cdot \mathbf a = \mathbf b \cdot \mathbf b = 1.</math>

This implies that

: <math>a = -||BE||, \quad b = -||BC||.</math>

According to the [[vector triple product]] formula,

: <math>\mathbf r \times (\mathbf a \times \mathbf b) = (\mathbf r \cdot \mathbf b) \mathbf a - (\mathbf r \cdot \mathbf a) \mathbf b.</math>

Then

:<math>\begin{align}
\mathbf r \times (\mathbf a \times \mathbf b) 
 &= \mathbf a \, || \mathbf r || \, || \mathbf b || \cos\angle DBC - \mathbf b \, ||\mathbf r || \, || \mathbf a || \cos\angle DBE \\[4pt]
 &= \mathbf a \, || \mathbf r || \cos\angle DBC - \mathbf b \, || \mathbf r || \cos\angle DBE \\[4pt]
 &= \mathbf a \, || BC || - \mathbf b \, || BE ||  \\[4pt]
 &= -b \mathbf a - (-a) \mathbf b \\[4pt]
 &= a \mathbf b - b \mathbf a\ \  \blacksquare
\end{align}</math>

When <math>|| \mathbf r || = 0,</math> the line {{mvar|L}} passes the origin with direction {{math|'''d'''}}. If <math>|| \mathbf r || > 0,</math>  the line has direction {{math|'''d'''}}; the plane that includes the origin and the line {{mvar|L}} has normal vector {{math|'''m'''}}; the line is tangent to a circle on that plane (normal to {{math|'''m'''}} and perpendicular to the plane of the picture) centered at the origin and with radius <math>|| \mathbf r ||.</math>

: '''Example.''' Let {{math|1=''a''<sub>0</sub> = 2}}, {{math|1='''a''' = (−1, 0, 0)}} and {{math|1=''b''<sub>0</sub> = −7}}, {{math|1='''b''' = (0, 7, −2)}}. Then {{math|1=('''d''' : '''m''') = (0 : −2 : −7 : −7 : 14 : −4)}}.

Although the usual algebraic definition tends to obscure the relationship, {{math|('''d''' : '''m''')}} are the Plücker coordinates of {{mvar|L}}.

{{Collapse bottom}}

== Algebraic definition ==

=== Primal coordinates ===
In a 3-dimensional projective space {{tmath|\mathbb P^3}}, let {{mvar|L}} be a line through distinct points {{mvar|x}} and {{mvar|y}} with [[homogeneous coordinates]] {{math|(''x''{{sub|0}} : ''x''{{sub|1}} : ''x''{{sub|2}} : ''x''{{sub|3}})}} and {{math|(''y''{{sub|0}} : ''y''{{sub|1}} : ''y''{{sub|2}} : ''y''{{sub|3}})}}.

The Plücker coordinates {{mvar|p{{sub|ij}}}} are defined as follows:
:<math>p_{ij} = \begin{vmatrix} x_{i} & y_{i} \\ x_{j} & y_{j}\end{vmatrix} = x_{i}y_{j}-x_{j}y_{i} . </math>

(the skew symmetric matrix whose elements are {{mvar|p<sub>ij</sub>}} is also called the [[Plücker matrix]] )<br>
This implies {{math|1=''p<sub>ii</sub>'' = 0}} and {{math|1=''p<sub>ij</sub>'' = −''p<sub>ji</sub>''}}, reducing the possibilities to only six (4 [[binomial coefficient|choose]] 2) independent quantities. The sextuple

: <math>(p_{01}:p_{02}:p_{03}:p_{23}:p_{31}:p_{12}) </math>

is uniquely determined by {{mvar|L}} up to a common nonzero scale factor. Furthermore, not all six components can be zero.
Thus the Plücker coordinates of {{mvar|L}} may be considered as homogeneous coordinates of a point in a 5-dimensional projective space, as suggested by the colon notation.

To see these facts, let {{mvar|M}} be the 4×2 matrix with the point coordinates as columns.

: <math> M = \begin{bmatrix} x_0 &  y_0 \\ x_1 & y_1 \\ x_2 & y_2 \\ x_3 & y_3 \end{bmatrix}</math>

The Plücker coordinate {{mvar|p<sub>ij</sub>}} is the determinant of rows {{mvar|i}} and {{mvar|j}} of {{mvar|M}}.
Because {{math|'''x'''}} and {{math|'''y'''}} are distinct points, the columns of {{mvar|M}} are [[linear independence|linearly independent]]; {{mvar|M}} has [[rank (linear algebra)|rank]] 2. Let {{mvar|M&prime;}} be a second matrix, with columns {{math|'''x&prime;''', '''y&prime;'''}} a different pair of distinct points on {{mvar|L}}. Then the columns of {{mvar|M&prime;}} are [[linear combination]]s of the columns of {{mvar|M}}; so for some 2×2 [[nonsingular matrix]] {{math|Λ}},

: <math> M' = M\Lambda . </math>

In particular, rows {{mvar|i}} and {{mvar|j}} of {{mvar|M&prime;}} and {{mvar|M}} are related by

: <math> \begin{bmatrix} x'_{i} & y'_{i}\\x'_{j}& y'_{j} \end{bmatrix} = \begin{bmatrix} x_{i} & y_{i}\\x_{j}& y_{j} \end{bmatrix} \begin{bmatrix} \lambda_{00} & \lambda_{01} \\ \lambda_{10} & \lambda_{11} \end{bmatrix} . </math>

Therefore, the determinant of the left side 2×2 matrix equals the product of the determinants of the right side 2×2 matrices, the latter of which is a fixed scalar, {{math|det Λ}}. Furthermore, all six 2×2 subdeterminants in {{mvar|M}} cannot be zero because the rank of {{mvar|M}} is 2.

=== Plücker map ===
Denote the set of all lines (linear images of {{tmath|\mathbb P^1}}) in {{tmath|\mathbb P^3}} by {{math|''G''<sub>1,3</sub>}}.  We thus have a map:
:<math>\begin{align}
\alpha \colon \mathrm{G}_{1,3} & \rightarrow \mathbb P^5 \\
L & \mapsto L^{\alpha},
\end{align}</math>
where
:<math> L^{\alpha}=(p_{01}:p_{02}:p_{03}:p_{23}:p_{31}:p_{12}) . </math>

=== Dual coordinates ===
Alternatively, a line can be described as the intersection of two planes. Let {{mvar|L}} 
be a line contained in distinct planes {{math|'''a'''}} and {{math|'''b'''}} with homogeneous coefficients {{math|(''a''{{sup|0}} : ''a''{{sup|1}} : ''a''{{sup|2}} : ''a''{{sup|3}})}} and {{math|(''b''{{sup|0}} : ''b''{{sup|1}} : ''b''{{sup|2}} : ''b''{{sup|3}})}}, respectively. (The first plane equation is <math display=inline>\sum_k a^k x_k =0,</math> for example.) The dual Plücker coordinate {{mvar|p<sup>ij</sup>}} is

:<math>p^{ij} = \begin{vmatrix} a^{i} & a^{j} \\ b^{i} & b^{j}\end{vmatrix}  = a^{i}b^{j}-a^{j}b^{i} . </math>

Dual coordinates are convenient in some computations, and they are equivalent to primary coordinates:
: <math>
(p_{01}:p_{02}:p_{03}:p_{23}:p_{31}:p_{12})=
(p^{23}:p^{31}:p^{12}:p^{01}:p^{02}:p^{03})
</math>
Here, equality between the two vectors in homogeneous coordinates means that the numbers on the right side are equal to the numbers on the left side up to some common scaling factor {{math|λ}}. Specifically, let {{math|(''i'', ''j'', ''k'', ''ℓ'')}} be an [[even permutation]] of {{math|(0, 1, 2, 3)}}; then

: <math>p_{ij} = \lambda p^{k\ell} .  </math>

=== Geometry ===
To relate back to the geometric intuition, take {{math|1=''x''<sub>0</sub> = 0}} as the plane at infinity; thus the coordinates of points ''not'' at infinity can be normalized so that {{math|1=''x''<sub>0</sub> = 1}}. Then {{mvar|M}} becomes

: <math> M = \begin{bmatrix} 1 &  1 \\ x_1 & y_1 \\ x_2& y_2 \\ x_3 & y_3 \end{bmatrix} , </math>

and setting <math>x=(x_1,x_2,x_3)</math> and <math>y=(y_1,y_2,y_3)</math>, we have <math>d=(p_{01},p_{02},p_{03})</math>and <math>m=(p_{23},p_{31},p_{12})</math>.

Dually, we have <math>d=(p^{23},p^{31},p^{12})</math> and <math>m=(p^{01},p^{02},p^{03}).</math>

==Bijection between lines and Klein quadric==
=== Plane equations ===
If the point <math>\mathbf z = (z_0:z_1:z_2:z_3)</math> lies on {{mvar|L}}, then the columns of

: <math> \begin{bmatrix} x_0 & y_0 & z_0 \\ x_1 & y_1 & z_1 \\ x_2 & y_2 & z_2 \\ x_3 & y_3 & z_3 \end{bmatrix} </math>

are [[linearly dependent]], so that the rank of this larger matrix is still 2. This implies that all 3×3 submatrices have determinant zero, generating four (4 choose 3) plane equations, such as

: <math>
\begin{align}
0 & = \begin{vmatrix} x_0 & y_0 & z_0 \\ x_1 & y_1 & z_1 \\ x_2 & y_2 & z_2 \end{vmatrix} \\[5pt]
& = \begin{vmatrix} x_1 & y_1 \\ x_2 & y_2 \end{vmatrix} z_0 - \begin{vmatrix} x_0 & y_0 \\ x_2 & y_2 \end{vmatrix} z_1 + \begin{vmatrix} x_0 & y_0 \\ x_1 & y_1 \end{vmatrix} z_2 \\[5pt]
& = p_{12} z_0 - p_{02} z_1 + p_{01} z_2 . \\[5pt]
& = p^{03} z_0 + p^{13} z_1 + p^{23} z_2 .
\end{align}
</math>

The four possible planes obtained are as follows.

: <math> \begin{matrix}
0 & = & {}+ p_{12} z_0 & {}- p_{02} z_1 & {}+ p_{01} z_2 & \\
0 & = & {}- p_{31} z_0 & {}- p_{03} z_1 & & {}+ p_{01} z_3 \\
0 & = & {}+p_{23} z_0 & & {}- p_{03} z_2 & {}+ p_{02} z_3 \\
0 & = & & {}+p_{23} z_1 & {}+ p_{31} z_2 & {}+ p_{12} z_3
\end{matrix} </math>

Using dual coordinates, and letting {{math|(''a''{{sup|0}} : ''a''{{sup|1}} : ''a''{{sup|2}} : ''a''{{sup|3}})}} be the line coefficients, each of these is simply {{math|1=''a<sup>i</sup>'' = ''p<sup>ij</sup>''}}, or

: <math> 0 = \sum_{i=0}^3 p^{ij} z_i , \qquad j = 0,\ldots,3 . </math>

Each Plücker coordinate appears in two of the four equations, each time multiplying a different variable; and as at least one of the coordinates is nonzero, we are guaranteed non-vacuous equations for two distinct planes intersecting in {{mvar|L}}. Thus the Plücker coordinates of a line determine that line uniquely, and the map α is an [[Injective function|injection]].

=== Quadratic relation ===
The image of {{math|α}} is not the complete set of points in {{tmath|\mathbb P^5}}; the Plücker coordinates of a line {{mvar|L}} satisfy the quadratic Plücker relation

: <math>
\begin{align}
0 & = p_{01}p^{01}+p_{02}p^{02}+p_{03}p^{03} \\
& = p_{01}p_{23}+p_{02}p_{31}+p_{03}p_{12}.
\end{align}
</math>

For proof, write this homogeneous polynomial as determinants and use [[Laplace expansion]] (in reverse).

: <math>
\begin{align}
0 & = \begin{vmatrix}x_0&y_0\\x_1&y_1\end{vmatrix}\begin{vmatrix}x_2&y_2\\x_3&y_3\end{vmatrix}+
\begin{vmatrix}x_0&y_0\\x_2&y_2\end{vmatrix}\begin{vmatrix}x_3&y_3\\x_1&y_1\end{vmatrix}+
\begin{vmatrix}x_0&y_0\\x_3&y_3\end{vmatrix}\begin{vmatrix}x_1&y_1\\x_2&y_2\end{vmatrix} \\[5pt]
& = (x_0 y_1-y_0 x_1)\begin{vmatrix}x_2&y_2\\x_3&y_3\end{vmatrix}-
(x_0 y_2-y_0 x_2)\begin{vmatrix}x_1&y_1\\x_3&y_3\end{vmatrix}+
(x_0 y_3-y_0 x_3)\begin{vmatrix}x_1&y_1\\x_2&y_2\end{vmatrix} \\[5pt]
& = x_0 \left(y_1\begin{vmatrix}x_2&y_2\\x_3&y_3\end{vmatrix}-
y_2\begin{vmatrix}x_1&y_1\\x_3&y_3\end{vmatrix}+
y_3\begin{vmatrix}x_1&y_1\\x_2&y_2\end{vmatrix}\right)
-y_0 \left(x_1\begin{vmatrix}x_2&y_2\\x_3&y_3\end{vmatrix}-
x_2\begin{vmatrix}x_1&y_1\\x_3&y_3\end{vmatrix}+
x_3\begin{vmatrix}x_1&y_1\\x_2&y_2\end{vmatrix}\right) \\[5pt]
& = x_0 \begin{vmatrix}x_1&y_1&y_1\\x_2&y_2&y_2\\x_3&y_3&y_3\end{vmatrix}
-y_0 \begin{vmatrix}x_1&x_1&y_1\\x_2&x_2&y_2\\x_3&x_3&y_3\end{vmatrix}
\end{align}
</math>

Since both 3×3 determinants have duplicate columns, the right hand side is identically zero.

Another proof may be done like this:
Since vector

: <math> d = \left( p_{01}, p_{02}, p_{03} \right) </math>

is perpendicular to vector

: <math> m = \left( p_{23}, p_{31}, p_{12} \right) </math>

(see above), the scalar product of {{mvar|d}} and {{mvar|m}} must be zero. q.e.d.

=== Point equations ===
Letting {{math|(''x''{{sub|0}} : ''x''{{sub|1}} : ''x''{{sub|2}} : ''x''{{sub|3}})}} be the point coordinates, four possible points on a line each have coordinates {{math|1=''x<sub>i</sub>'' = ''p<sub>ij</sub>''}}, for {{math|1=''j'' = 0, 1, 2, 3}}. Some of these possible points may be inadmissible because all coordinates are zero, but since at least one Plücker coordinate is nonzero, at least two distinct points are guaranteed.

=== Bijectivity ===
If <math>(q_{01}:q_{02}:q_{03}:q_{23}:q_{31}:q_{12})</math> are the homogeneous coordinates of a point in {{tmath|\mathbb P^5}}, without loss of generality assume that {{math|''q''<sub>01</sub>}} is nonzero. Then the matrix

: <math> M = \begin{bmatrix} q_{01} & 0 \\ 0 & q_{01} \\ -q_{12} & q_{02} \\ q_{31} & q_{03} \end{bmatrix} </math>

has rank 2, and so its columns are distinct points defining a line {{mvar|L}}. When the {{tmath|\mathbb P^5}} coordinates, {{mvar|q<sub>ij</sub>}}, satisfy the quadratic Plücker relation, they are the Plücker coordinates of {{mvar|L}}. To see this, first normalize {{math|''q''<sub>01</sub>}} to 1. Then we immediately have that for the Plücker coordinates computed from {{mvar|M}}, {{math|1=''p<sub>ij</sub>'' = ''q<sub>ij</sub>''}}, except for

: <math> p_{23} = - q_{03} q_{12} - q_{02} q_{31} . </math>

But if the {{mvar|q<sub>ij</sub>}} satisfy the Plücker relation 
:<math>q_{23} + q_{02}q_{31} + q_{03}q_{12} = 0,</math> 
then {{math|1=''p''<sub>23</sub> = ''q''<sub>23</sub>}}, 
completing the set of identities.

Consequently, {{math|α}} is a [[surjection]] onto the [[algebraic variety]] consisting of the set of zeros of the quadratic polynomial

: <math> p_{01}p_{23}+p_{02}p_{31}+p_{03}p_{12} . </math>

And since {{math|α}} is also an injection, the lines in {{tmath|\mathbb P^3}} are thus in [[bijection|bijective]] correspondence with the points of this [[quadric]] in {{tmath|\mathbb P^5}}, called the Plücker quadric or [[Klein quadric]].

==Uses==
Plücker coordinates allow concise solutions to problems of line geometry in 3-dimensional space, especially those involving [[incidence (geometry)|incidence]].

=== Line-line crossing ===
{{further|Line-line intersection}}

Two lines in {{tmath|\mathbb P^3}} are either [[skew lines|skew]] or [[coplanar]], and in the latter case they are either coincident or intersect in a unique point. If {{mvar|p<sub>ij</sub>}} and {{mvar|p&prime;<sub>ij</sub>}} are the Plücker coordinates of two lines, then they are coplanar precisely when

:<math>\mathbf d \cdot \mathbf m' + \mathbf m \cdot \mathbf d' = 0,</math>

as shown by

: <math>
\begin{align}
0 & = p_{01}p'_{23} + p_{02}p'_{31} + p_{03}p'_{12} + p_{23}p'_{01} + p_{31}p'_{02} + p_{12}p'_{03} \\[5pt]
& = \begin{vmatrix}x_0&y_0&x'_0&y'_0\\
x_1&y_1&x'_1&y'_1\\
x_2&y_2&x'_2&y'_2\\
x_3&y_3&x'_3&y'_3\end{vmatrix}.
\end{align}
</math>

When the lines are skew, the sign of the result indicates the sense of crossing: positive if a right-handed screw takes {{mvar|L}} into {{mvar|L&prime;}}, else negative.

The quadratic Plücker relation essentially states that a line is coplanar with itself.

=== Line-line join ===
In the event that two lines are coplanar but not parallel, their common plane has equation

:<math>0 = (\mathbf m \cdot \mathbf d')x_0 + (\mathbf d \times \mathbf d')\cdot \mathbf x,</math>

where <math>x=(x_1,x_2,x_3).</math>

The slightest perturbation will destroy the existence of a common plane, and near-parallelism of the lines will cause numeric difficulties in finding such a plane even if it does exist.

=== Line-line meet ===
Dually, two coplanar lines, neither of which contains the origin, have common point

: <math>(x_0:\mathbf x) = (\mathbf d \cdot \mathbf m': \mathbf m \times \mathbf m').</math>

To handle lines not meeting this restriction, see the references.

=== Plane-line meet ===
{{further|Plane-line intersection}}

Given a plane with equation

: <math> 0 =  a^0x_0 + a^1x_1 + a^2x_2 + a^3x_3 , </math>

or more concisely,
:<math>0 = a^0x_0 + \mathbf a \cdot \mathbf x;</math>
and given a line not in it with Plücker coordinates {{math|('''d''' : '''m''')}}, then their point of intersection is

: <math>(x_0 : \mathbf x) = (\mathbf a \cdot \mathbf d : \mathbf a \times \mathbf m - a_0\mathbf d) .</math>

The point coordinates, {{math|(''x''{{sub|0}} : ''x''{{sub|1}} : ''x''{{sub|2}} : ''x''{{sub|3}})}}, can also be expressed in terms of Plücker coordinates as

: <math> x_i = \sum_{j \ne i} a^j p_{ij} , \qquad i = 0 \ldots 3 . </math>

=== Point-line join ===
Dually, given a point {{math|(''y''<sub>0</sub> : '''y''')}} and a line not containing it, their common plane has equation

: <math>0 = (\mathbf y \cdot \mathbf m) x_0 + (\mathbf y \times \mathbf d - y_0 \mathbf m)\cdot \mathbf x.</math>

The plane coordinates, {{math|(''a''{{sup|0}} : ''a''{{sup|1}} : ''a''{{sup|2}} : ''a''{{sup|3}})}}, can also be expressed in terms of dual Plücker coordinates as

: <math> a^i = \sum_{j \ne i} y_j p^{ij} , \qquad i = 0 \ldots 3 . </math>

=== Line families ===
Because the [[Klein quadric]] is in {{tmath|\mathbb P^5}}, it contains linear subspaces of dimensions one and two (but no higher). These correspond to one- and two-parameter families of lines in {{tmath|\mathbb P^3}}.

For example, suppose {{mvar|L, L&prime;}} are distinct lines in {{tmath|\mathbb P^3}} determined by points {{math|'''x''', '''y'''}} and {{math|'''x'''&prime;, '''y'''&prime;}}, respectively. Linear combinations of their determining points give linear combinations of their Plücker coordinates, generating a one-parameter family of lines containing {{mvar|L}} and {{math|''L''&prime;}}. This corresponds to a one-dimensional linear subspace belonging to the Klein quadric.

==== Lines in plane ====
If three distinct and non-parallel lines are coplanar; their linear combinations generate a two-parameter family of lines, all the lines in the plane. This corresponds to a two-dimensional linear subspace belonging to the Klein quadric.

==== Lines through point ====
If three distinct and non-coplanar lines intersect in a point, their linear combinations generate a two-parameter family of lines, all the lines through the point. This also corresponds to a two-dimensional linear subspace belonging to the Klein quadric.

==== Ruled surface ====
A [[ruled surface]] is a family of lines that is not necessarily linear. It corresponds to a curve on the Klein quadric. For example, a [[hyperboloid of one sheet]] is a quadric surface in {{tmath|\mathbb P^3}} ruled by two different families of lines, one line of each passing through each point of the surface; each family corresponds under the Plücker map to a [[conic section]] within the Klein quadric in {{tmath|\mathbb P^5}}.

=== Line geometry ===
During the nineteenth century, ''line geometry'' was studied intensively. In terms of the bijection given above, this is a description of the intrinsic geometry of the Klein quadric.

==== Ray tracing ====
Line geometry is extensively used in [[Ray tracing (graphics)|ray tracing]] application where the geometry and intersections of rays need to be calculated in 3D. An implementation is described in
[http://www.flipcode.com/archives/Introduction_To_Plcker_Coordinates.shtml Introduction to Plücker Coordinates] written for the Ray Tracing forum by Thouis Jones.

== See also ==
* [[Real projective plane#The flat projective plane|Flat projective plane]]
* [[Plücker matrix]]

== References ==
* {{cite book
  | last = Hodge
  | first = W. V. D.
  | author-link = W. V. D. Hodge
  |author2=D. Pedoe 
  | title = Methods of Algebraic Geometry, Volume&nbsp;I (Book&nbsp;II)
  | publisher = [[Cambridge University Press]]
  | year = 1994
  | isbn = 978-0-521-46900-5
  | orig-year = 1947 | author2-link = D. Pedoe
 }}
* {{cite book
  |editor-last = Behnke
  |editor-first = H.
 |editor2=F. Bachmann |editor3=K. Fladt |editor4=H. Kunle
  | others = trans. S. H. Gould
  | title = Fundamentals of Mathematics, Volume II: Geometry
  | publisher = [[MIT Press]]
  | year = 1984
  | isbn = 978-0-262-52094-2 }}<br />From the German:  ''Grundzüge der Mathematik, Band II: Geometrie''. Vandenhoeck & Ruprecht.
* {{cite journal
  | last = Guilfoyle
  | first = B.
  |author2=W. Klingenberg
   | title = On the space of oriented affine lines in R^3
  | journal = [[Archiv der Mathematik]]
  | volume = 82
  | issue = 1
  | pages = 81–84
  | publisher = [[Birkhäuser]]
  | year = 2004
  | doi = 10.1007/s00013-003-4861-3
 | issn = 0003-889X
  | arxiv = math/0405189
  | s2cid = 118352042
 }}
*{{SpringerEOM|title=Plücker coordinates|first=L.P.|last= Kuptsov}}
* {{cite book
  | last = Mason
  | first = Matthew T.
  |author2=J. Kenneth Salisbury
   | title = Robot Hands and the Mechanics of Manipulation
  | publisher = [[MIT Press]]
  | year = 1985
  | isbn = 978-0-262-13205-3}}
* {{cite book
  | last = Hartley
  | first = R.~I.
  |author2=Zisserman A.
   | title = Multiple View Geometry in Computer Vision
  | publisher = [[Cambridge University Press]]
  | year = 2004
  | url = http://www.robots.ox.ac.uk/~vgg/hzbook/
  | isbn = 0521540518}}
* {{cite journal
  | last = Hohmeyer
  | first = M.
  |author2=S. Teller
   | title = Determining the Lines Through Four Lines
  | journal = [[Journal of Graphics Tools]]
  | volume = 4
  | issue = 3
  | pages = 11–22
  | publisher = [[A K Peters]]
  | year = 1999
  | doi = 10.1080/10867651.1999.10487506
 | url = http://people.csail.mit.edu/seth/pubs/TellerHohmeyerJGT2000.pdf
  | issn = 1086-7651}}
* {{cite book
  | last = Shafarevich
  | first = I. R. 
  | author-link = Igor Shafarevich
  |author2=A. O. Remizov
  | title = Linear Algebra and Geometry
  | publisher = [[Springer Science+Business Media|Springer]]
  | year = 2012
  | url = https://www.springer.com/mathematics/algebra/book/978-3-642-30993-9
  | isbn = 978-3-642-30993-9}}
*{{cite report
  | first = Yan-Bin
  | last = Jia
  | title = Plücker Coordinates for Lines in the Space
  | url = http://web.cs.iastate.edu/~cs577/handouts/plucker-coordinates.pdf
  | year = 2017}}
*{{cite web
 |url=http://www.realtimerendering.com/resources/RTNews/html/rtnv11n1.html#art3
 |archive-url=https://archive.today/20221018165206/https://www.realtimerendering.com/resources/RTNews/html/rtnv11n1.html%23art3
 |title=Plücker Coordinate Tutorial
 |last=Shoemake
 |first=Ken
 |date=1998
 |publisher=Ray Tracing News
 |access-date=4 July 2018
 |archive-date=18 October 2022
 }}
{{Authority control}}

{{DEFAULTSORT:Plucker Coordinates}}
[[Category:Projective geometry]]
[[Category:Multilinear algebra]]
[[Category:Geometric algebra]]
[[Category:Coordinate systems]]