View source for Kozai mechanism

{{Unreferenced|date=December 2009}}
[[Image:Plane-line intersection.svg|thumb|350px|right|The three possible plane-line intersections:<br /> 1. No intersection.<br />
2. Point intersection.<br />
3. Line intersection.]]
In analytic [[geometry]], the intersection of a [[line (mathematics)|line]] and a [[plane (mathematics)|plane]] can be the [[empty set]], 
a [[point (geometry)|point]], or 
a line. Distinguishing these cases, and determining equations for the point and line in the latter cases have use, for example, in [[computer graphics]], [[motion planning]], and [[collision detection]].

==Parametric form==
[[Image:Line plane.svg|thumb|300px|right|The intersection of line and plane.]]
A line is described by all points that are a given direction from a point. Thus a general point on a line can be represented as 

:<math>\mathbf{l}_a + (\mathbf{l}_b - \mathbf{l}_a)t, \quad t\in \mathbb{R}</math>

where <math>\mathbf{l}_a=(x_a, y_a, z_a)</math> and <math>\mathbf{l}_b=(x_b, y_b, z_b)</math> 
are two distinct points along the line.

Similarly a general point on a plane can be represented as 

:<math>\mathbf{p}_0 + (\mathbf{p}_1-\mathbf{p}_0)u + (\mathbf{p}_2-\mathbf{p}_0)v, \quad u,v\in\mathbb{R}</math>

where <math>\mathbf{p}_k=(x_k,y_k,z_k)</math>,
<math>k=0,1,2</math> are three points in the plane which are not co-linear.

The point at which the line intersects the plane is therefore described by setting the point on the line equal to the point on the plane, giving the parametric equation:
:<math>\mathbf{l}_a + (\mathbf{l}_b - \mathbf{l}_a)t = \mathbf{p}_0 + (\mathbf{p}_1-\mathbf{p}_0)u + (\mathbf{p}_2-\mathbf{p}_0)v</math>
This can be simplified to
:<math>\mathbf{l}_a - \mathbf{p}_0 =  (\mathbf{l}_a - \mathbf{l}_b)t + (\mathbf{p}_1-\mathbf{p}_0)u + (\mathbf{p}_2-\mathbf{p}_0)v,</math>
which can be expressed in matrix form as:
:<math> \begin{bmatrix} x_a - x_0 \\ y_a - y_0 \\ z_a - z_0 \end{bmatrix}  = \begin{bmatrix} x_a - x_b & x_1 - x_0 & x_2 - x_0 \\ y_a - y_b & y_1 - y_0 & y_2 - y_0 \\ z_a - z_b & z_1 - z_0 & z_2 - z_0 \end{bmatrix} \begin{bmatrix} t \\ u \\ v \end{bmatrix} </math>

The point of intersection is then equal to
:<math>\mathbf{l}_a + (\mathbf{l}_b - \mathbf{l}_a)t</math>

If the line is parallel to the plane then the vectors <math>\mathbf{l}_b - \mathbf{l}_a</math>, <math>\mathbf{p}_1-\mathbf{p}_0</math>, and <math>\mathbf{p}_2-\mathbf{p}_0</math> will be [[linearly dependent]] and the matrix will be singular. This situation will also occur when the line lies in the plane.

If the solution satisfies the condition <math>t \in [0,1],</math>, then the intersection point is on the line between <math>\mathbf{l}_a</math> and <math>\mathbf{l}_b</math>.

If the solution satisfies
:<math>u,v \in [0,1], \;\;\; (u+v) \leq 1,</math>
then the intersection point is in the plane inside the triangle spanned by the three points <math>\mathbf{p}_0</math>, <math>\mathbf{p}_1</math> and <math>\mathbf{p}_2</math>. 

This problem is typically solved by expressing it in matrix form, and [[Matrix inversion|inverting]] it:
:<math> \begin{bmatrix} t \\ u \\ v \end{bmatrix} = \begin{bmatrix} x_a - x_b & x_1 - x_0 & x_2 - x_0 \\ y_a - y_b & y_1 - y_0 & y_2 - y_0 \\ z_a - z_b & z_1 - z_0 & z_2 - z_0 \end{bmatrix}^{-1} \begin{bmatrix} x_a - x_0 \\ y_a - y_0 \\ z_a - z_0 \end{bmatrix}.</math>

==Algebraic form==
In [[vector notation]], a plane can be expressed as the set of points <math>\mathbf{p}</math> for which
:<math>(\mathbf{p}-\mathbf{p_0})\cdot\mathbf{n} = 0</math>
where <math>\mathbf{n}</math> is a [[normal vector]] to the plane and <math>\mathbf{p_0}</math> is a point on the plane.
The vector equation for a line is
:<math>\mathbf{p} = d\mathbf{l} + \mathbf{l_0} \quad    d\in\mathbb{R}</math>
where <math>\mathbf{l}</math> is a vector in the direction of the line, <math>\mathbf{l_0}</math> is a point on the line, and  <math>d</math> is a scalar in the real number domain. Substitute the line into the plane equation to get
:<math>(d \mathbf{l} + \mathbf{l_0} - \mathbf{p_0})\cdot\mathbf{n} = 0</math>
Distribute <math>\mathbf{n}</math> to get
:<math>d \mathbf{l}\cdot\mathbf{n} + (\mathbf{l_0}-\mathbf{p_0})\cdot\mathbf{n} = 0</math>
And solve for <math>d</math>
:<math>d = {(\mathbf{p_0}-\mathbf{l_0})\cdot\mathbf{n} \over \mathbf{l}\cdot\mathbf{n}}</math>

If the line starts outside the plane and is parallel to the plane, there is no intersection. In this case, the above denominator will be zero and the numerator will be non-zero.  If the line starts inside the plane and is parallel to the plane, the line intersects the plane everywhere. In this case, both the numerator and denominator above will be zero. In all other cases, the line intersects the plane once and <math>d</math> represents the intersection as the distance along the line from <math>\mathbf{l_0}</math>, i.e. <math>d\mathbf{l} + \mathbf{l_0}</math>

==Uses==
In the [[Ray tracing (graphics)|ray tracing]] method of [[computer graphics]] a surface can be represented as a set of pieces of planes. The intersection of a ray of light with each plane is used to produce an image of the surface. In vision-based 3D reconstruction, a subfield of computer vision, depth values are commonly measured by so-called triangulation method, which finds the intersection between light plane and ray reflected toward camera.

The algorithm can be generalised to cover intersection with other planar figures, in particular, the [[intersection of a polyhedron with a line]].

{{DEFAULTSORT:Line-Plane Intersection}}
[[Category:Euclidean geometry]]
[[Category:Computational physics]]
[[Category:Geometric algorithms]]

[[cs:Analytická geometrie#Vzájemná poloha dvou rovin v třírozměrném prostoru]]