View source for Intersection (geometry)

{{short description|Shape formed from points common to other shapes}}
[[File:Schnittpunkt-2g.svg|200px|thumb|The red dot represents the point at which the two lines intersect.]]

In [[Geometry|geometry]], an '''intersection''' is a point, line, or curve common to two or more objects (such as lines, curves, planes, and surfaces). The simplest case in [[Euclidean geometry]] is the [[Line–line intersection|line–line intersection]] between two distinct [[Line (geometry)|lines]], which either is one [[Point (geometry)|point]] (sometimes called a ''[[Vertex (geometry)|vertex]]'') or does not exist (if the lines are parallel). Other types of geometric intersection include:
* [[Line–plane intersection]]
* [[Line–sphere intersection]]
* [[Intersection of a polyhedron with a line]]
* [[Line segment intersection]]
* [[Intersection curve]]

Determination of the intersection of [[Flat (geometry)|flats]] – linear geometric objects embedded in a higher-[[Physics:Dimension|dimension]]al space – is a simple task of [[Linear algebra|linear algebra]], namely the solution of a [[System of linear equations|system of linear equations]]. In general the determination of an intersection leads to non-linear equations, which can be solved numerically, for example using Newton iteration. Intersection problems between a line and a [[Conic section|conic section]] (circle, ellipse, parabola, etc.) or a [[Quadric|quadric]] (sphere, cylinder, hyperboloid, etc.) lead to [[Quadratic equation|quadratic equation]]s that can be easily solved. Intersections between quadrics lead to [[Quartic equation|quartic equation]]s that can be solved [[Algebraic equation|algebraically]].

== On a plane ==
{{further|Plane (geometry)|Two-dimensional space}}

=== Two lines ===
{{main|Line–line intersection}}
For the determination of the intersection point of two non-parallel lines

<math>a_1x+b_1y=c_1, \ a_2x+b_2y=c_2 </math> 

one gets, from [[Cramer's rule]] or by substituting out a variable, the coordinates of the intersection point <math>(x_s,y_s)</math> :
:<math> x_s=\frac{c_1b_2-c_2b_1}{a_1b_2-a_2b_1} , \quad y_s=\frac{a_1c_2-a_2c_1}{a_1b_2-a_2b_1}. \ </math> 
(If <math> a_1b_2-a_2b_1=0</math> the lines are parallel and these formulas cannot be used because they involve dividing by 0.)

=== Two line segments ===<!-- caution: an internal link -->

[[File:Is-linesegm.svg|300px|thumb|Intersection of two line segments]]
For two non-parallel [[Line segment|line segment]]s <math>(x_1,y_1),(x_2,y_2)</math> and <math>(x_3,y_3),(x_4,y_4)</math> there is not necessarily an intersection point (see diagram), because the intersection point <math>(x_0,y_0)</math> of the corresponding lines need not to be contained in the line segments. In order to check the situation one uses parametric representations of the lines:
: <math> (x(s),y(s))=(x_1+s(x_2-x_1),y_1+s(y_2-y_1)),</math>
: <math> (x(t),y(t))=(x_3+t(x_4-x_3),y_3+t(y_4-y_3)). </math>
The line segments intersect only in a common point <math>(x_0,y_0)</math> of the corresponding lines if the corresponding parameters <math>s_0,t_0</math> fulfill the condition <math> 0\le s_0,t_0 \le 1 </math>. 
The parameters <math>s_0,t_0 </math> are the solution of the linear system 
: <math>s(x_2-x_1)-t(x_4-x_3)=x_3-x_1,</math>
:<math> s(y_2-y_1)-t(y_4-y_3)=y_3-y_1 \ .</math>
It can be solved for ''s'' and ''t'' using Cramer's rule (see above). If the condition <math> 0\le s_0,t_0 \le 1 </math> is fulfilled one inserts <math>s_0</math> or <math>t_0</math> into the corresponding parametric representation and gets the intersection point <math>(x_0,y_0)</math>.

''Example:''  For the line segments <math>(1,1),(3,2)</math> and <math>(1,4),(2,-1)</math> one gets the linear system
:<math> 2s-t=0</math>
:<math>s+5t=3</math>
and <math>s_0=\tfrac{3}{11}, t_0=\tfrac{6}{11}</math>. That means: the lines intersect at point <math>(\tfrac{17}{11},\tfrac{14}{11})</math>.

''Remark:'' Considering lines, instead of segments, determined by pairs of points, each condition <math> 0\le s_0,t_0 \le 1 </math> can be dropped and the method yields the intersection point of the lines (see above).

[[File:Is-circle-line-s.svg|thumb|Line–circle intersection]]

=== A line and a circle ===
{{further|Line–sphere intersection}}
For the intersection of
*line <math>ax+by=c</math> and [[Circle|circle]] <math>x^2+y^2=r^2</math> 
one solves the line equation for {{mvar|x}} or {{mvar|y}} and [[Substitution (algebra)|substitutes]] it into the equation of the circle and gets for the solution (using the formula of a quadratic equation) <math>(x_1,y_1),(x_2,y_2)</math> with
:<math>x_{1/2}= \frac{ac\pm b\sqrt{r^2(a^2+b^2)-c^2}}{a^2+b^2} \ ,</math> 
:<math>y_{1/2}= \frac{bc\mp a\sqrt{r^2(a^2+b^2)-c^2}}{a^2+b^2} \ , </math>
if <math> r^2(a^2+b^2)-c^2 > 0 \ .</math> If this condition holds with strict inequality, there are two intersection points; in this case the line is called a [[Secant line|secant line]] of the circle, and the line segment connecting the intersection points is called a [[Chord (geometry)|chord]] of the circle.

If <math> r^2(a^2+b^2)-c^2=0 </math> holds,  there exists only one intersection point and the line is tangent to the circle. If the weak inequality does not hold, the line does not intersect the circle.
<br /> 

If the circle's midpoint is not the origin, see.<ref>Erich Hartmann: [http://www.mathematik.tu-darmstadt.de/~ehartmann/cdgen0104.pdf ''Geometry and Algorithms for COMPUTER AIDED DESIGN'']. Lecture notes, Technische Universität Darmstadt, October 2003, p. 17</ref> The intersection of a line and a parabola or hyperbola may be treated analogously.

=== Two circles ===

The determination of the intersection points of two circles
* <math>(x-x_1)^2+(y-y_1)^2=r_1^2 ,\ \quad  (x-x_2)^2+(y-y_2)^2=r_2^2</math>
can be reduced to the previous case of intersecting a line and a circle. By subtraction of the two given equations one gets the line equation: 
:<math>2(x_2-x_1)x+2(y_2-y_1)y=r_1^2-x_1^2-y_1^2-r_2^2+x_2^2+y_2^2. </math>
This special line is the radical line of the two circles.
[[File:Is-circle-circle-2x.svg|thumb|upright=1.2|Intersection of two circles with centers on the x-axis, their radical line is dark red]]
'''Special case <math>\;x_1=y_1=y_2=0 </math> :'''<br> 
In this case the origin is the center of the first circle and the second center lies on the x-axis (s. diagram). The equation of the radical line simplifies to <math>\;2x_2x=r_1^2-r_2^2+x_2^2\;</math> and the points of intersection can be written as <math>(x_0,\pm y_0)</math> with
:<math>x_0=\frac{r_1^2-r_2^2+x_2^2}{2x_2},\quad y_0 =\sqrt{r_1^2-x_0^2}\ .</math>
In case of <math>r_1^2<x_0^2</math> the circles have no points in common. <br>
In case of <math>r_1^2=x_0^2</math> the circles have one point in common and the radical line is a common tangent.

Any general case as written above can be transformed by a shift and a rotation into the special case.

The intersection of two [[Disk (mathematics)|disks]] (the interiors of the two circles) forms a shape called a [[Lens (geometry)|lens]].
[[File:Is-circle-ellipse-s.svg|thumb|circle–ellipse intersection]]

=== Two conic sections ===
The problem of intersection of an ellipse/hyperbola/parabola with another [[Conic section|conic section]] leads to a [[System of polynomial equations|system of quadratic equations]], which can be solved in special cases easily by elimination of one coordinate. Special properties of conic sections may be used to obtain a [[Conic section#Intersecting two conics|solution]]. In general the intersection points can be determined by solving the equation by a Newton iteration. If a) both conics are given implicitly (by an equation) a 2-dimensional Newton iteration b) one implicitly and the other parametrically given a 1-dimensional Newton iteration is necessary. See next section.

=== Two smooth curves ===
[[File:Schnittp2d-transv-s.svg|250px|thumb|A transversal intersection of two curves]]
[[File:Beruehr-schnitt-s.svg|350px|thumb|touching intersection (left),  touching (right)]]
Two curves in <math>\R^2</math> (two-dimensional space), which are continuously differentiable (i.e. there is no sharp bend),
have an intersection point, if they have a point of the plane in common and have at this point (see diagram):
: a) different tangent lines ('''''transversal intersection''''', after ''[[Transversality (mathematics)|transversality]]''), or
: b) the tangent line in common and they are crossing each other ('''''touching intersection''''', after ''tangency'').

If both the curves have a point {{mvar|S}} and the tangent line there in common but do not cross each other, they are just ''touching'' at point {{mvar|S}}.

Because touching intersections appear rarely and are difficult to deal with, the following considerations omit this case. In any case below all necessary differential conditions are presupposed. The determination of intersection points always leads to one or two non-linear equations which can be solved by Newton iteration. A list of the appearing cases follows:

[[File:Schnittp2d-pi-s.svg|thumb|intersection of a parametric curve and an implicit curve]]
[[File:Schnittp2d-ii-s.svg|thumb|intersection of two implicit curves]]
*If ''both curves are explicitly'' given: <math> y=f_1(x), \ y=f_2(x)</math>, equating them yields the equation
:: <math>f_1(x)=f_2(x) \ .</math>
*If ''both curves are parametrically'' given: <math>C_1: (x_1(t),y_1(t)), \ C_2: (x_2(s),y_2(s)).</math>
: Equating them yields two equations in two variables:
::  <math>x_1(t)=x_2(s), \ y_1(t)=y_2(s) \ .</math>
*If ''one curve is parametrically and the other implicitly'' given: <math>C_1: (x_1(t),y_1(t)), \ C_2: f(x,y)=0.</math>
:This is the simplest case besides the explicit case. One has to insert the parametric representation of <math>C_1</math> into the equation <math>f(x,y)=0</math> of curve <math>C_2</math> and one gets the equation:
::<math>f(x_1(t),y_2(t))=0 \ .</math>
*If ''both curves are implicitly'' given: <math>C_1: f_1(x,y)=0, \ C_2: f_2(x,y)=0.</math>
: Here, an intersection point is a solution of the system
::<math>f_1(x,y)=0, \ f_2(x,y)=0 \ .</math>
Any Newton iteration needs convenient starting values, which can be derived by a visualization of both the curves. A parametrically or explicitly given curve can easily be visualized, because to any parameter {{mvar|t}} or {{mvar|x}} respectively it is easy to calculate the corresponding point. For implicitly given curves this task is not as easy. In this case one has to determine a curve point with help of starting values and an iteration. See
.<ref>Erich Hartmann: [http://www.mathematik.tu-darmstadt.de/~ehartmann/cdgen0104.pdf ''Geometry and Algorithms for COMPUTER AIDED DESIGN'']. Lecture notes, Technische Universität Darmstadt, October 2003, p. 33</ref>

''Examples:''
:1: <math>C_1: (t,t^3)</math> and circle <math>C_2: (x-1)^2+(y-1)^2-10=0</math> (see diagram). 
:: The Newton iteration <math>t_{n+1}:=t_n-\frac{f(t_n)}{f'(t_n)}</math> for function 
:::<math>f(t)=(t-1)^2+(t^3-1)^2-10</math> has to be done. As start values one can choose −1 and 1.5.
::The intersection points are: (−1.1073, −1.3578), (1.6011, 4.1046)
:2:<math>C_1: f_1(x,y)=x^4+y^4-1=0,</math>
:: <math>C_2: f_2(x,y)=(x-0.5)^2+(y-0.5)^2-1=0 </math> (see diagram). 
:: The Newton iteration 
:::<math>{x_{n+1}\choose y_{n+1}}={x_{n}+\delta_x\choose y_n+\delta_y}</math> has to be performed, where <math>{\delta_x \choose \delta_y}</math> is the solution of the linear system
:::<math>\begin{pmatrix}
  \frac{\partial f_1}{\partial x} & \frac{\partial f_1}{\partial y} \\
  \frac{\partial f_2}{\partial x} & \frac{\partial f_2}{\partial y} 
 \end{pmatrix}{\delta_x \choose \delta_y}={-f_1\choose -f_2}
</math> at point <math>(x_n,y_n)</math>. As starting values one can choose(−0.5, 1) and (1, −0.5).
:: The linear system can be solved by Cramer's rule.
::The intersection points are (−0.3686, 0.9953) and (0.9953, −0.3686).

=== Two polygons ===
[[File:Is-polygpolyg-s.svg|300px|thumb|intersection of two polygons: window test]]
If one wants to determine the intersection points of two [[Polygon|polygon]]s, one can check the intersection of any pair of line segments of the polygons (see above). For polygons with many segments this method is rather time-consuming. In practice one accelerates the intersection algorithm by using ''window tests''. In this case one divides the polygons into small sub-polygons and determines the smallest window (rectangle with sides parallel to the coordinate axes) for any sub-polygon. Before starting the time-consuming determination of the intersection point of two line segments any pair of windows is tested for common points. See.<ref>Erich Hartmann: [http://www.mathematik.tu-darmstadt.de/~ehartmann/cdg-skript-1998.pdf ''CDKG: Computerunterstützte Darstellende und Konstruktive Geometrie'']. Lecture notes, TU Darmstadt, 1997, p. 79 (PDF; 3,4&nbsp;MB)</ref>

== In space (three dimensions) ==
{{more information|three-dimensional space}}
In 3-dimensional space there are intersection points (common points) between curves and surfaces. In the following sections we consider ''[[Transversality (mathematics)|transversal]] intersection'' only.

=== A line and a plane ===
{{main|Line–plane intersection}}
[[File:Schnittp-ger-eb-s.svg|thumb|Line–plane intersection]]
The intersection of a line and a plane ''in [[General position|general position]]'' in three dimensions is a point.

Commonly a line in space is represented parametrically  <math> (x(t),y(t),z(t)) </math> and a plane by an equation <math>ax+by+cz=d</math>. Inserting the parameter representation into the equation yields the linear equation 
:<math>ax(t)+by(t)+cz(t)=d\ ,</math>
for parameter <math>t_0</math> of the intersection point <math>(x(t_0),y(t_0),z(t_0))</math>.

If the linear equation has no solution, the line either lies on the plane or is parallel to it.

=== Three planes ===
If a line is defined by two intersecting planes <math>\varepsilon_i: \ \vec n_i\cdot\vec x=d_i, \ i=1,2</math> and should be intersected by a third plane <math>\varepsilon_3: \ \vec n_3\cdot\vec x=d_3 </math>, the common intersection point of the three planes has to be evaluated.

Three planes <math>\varepsilon_i: \ \vec n_i\cdot\vec x=d_i, \ i=1,2,3 </math> with linear independent normal vectors <math> \vec n_1,\vec n_2, \vec n_3</math> have the intersection point 
:<math> \vec p_0=\frac{d_1(\vec n_2\times \vec n_3) +d_2(\vec n_3\times \vec n_1) + d_3(\vec n_1\times \vec n_2)}{\vec n_1\cdot(\vec n_2\times \vec n_3)} \ .</math> 
For the proof one should establish <math>\vec n_i\cdot\vec p_0=d_i, \ i=1,2,3 , </math> using the rules of a [[Scalar triple product|scalar triple product]]. If the scalar triple product equals to 0, then planes either do not have the triple intersection or it is a line (or a plane, if all three planes are the same).

=== A curve and a surface ===
[[File:Is-pcurve-isurface-s.svg|250px|thumb|intersection of curve <math>(t,t^2,t^3)</math> with surface <math>x^4+y^4+z^4=1 </math>]]
Analogously to the plane case the following cases lead to non-linear systems, which can be solved using a 1- or 3-dimensional Newton iteration.<ref>Erich Hartmann: [http://www.mathematik.tu-darmstadt.de/~ehartmann/cdgen0104.pdf ''Geometry and Algorithms for COMPUTER AIDED DESIGN'']. Lecture notes, Technische Universität Darmstadt, October 2003, p. 93</ref>
*parametric curve <math>C: (x(t),y(t),z(t)) </math> and 
:parametric surface <math>S: (x(u,v),y(u,v),z(u,v))\ ,</math>
*parametric curve <math>C: (x(t),y(t),z(t)) </math> and 
:implicit surface <math> S: f(x,y,z)=0\ .</math>

'''Example:'''
:parametric curve <math>C: (t,t^2,t^3)</math>  and 
:implicit surface <math>S: x^4+y^4+z^4-1=0</math> (s. picture).
:The intersection points are: (−0.8587, 0.7374, −0.6332), (0.8587, 0.7374, 0.6332).

A [[Line–sphere intersection|line–sphere intersection]] is a simple special case.

Like the case of a line and a plane, the intersection of a curve and a surface ''in [[General position|general position]]'' consists of discrete points, but a curve may be partly or totally contained in a surface.

=== A line and a polyhedron ===
{{main|Intersection of a polyhedron with a line}}

=== Two surfaces ===
{{main|Intersection curve}}
Two transversally intersecting surfaces give an [[Intersection curve|intersection curve]]. The most simple case is the intersection line of two non-parallel planes.

=== A sphere and a plane ===
{{also|Spherical circle}}
When the intersection of a sphere and a plane is not empty or a single point, it is a circle. This can be seen as follows:

Let ''S'' be a sphere with center ''O'', ''P'' a plane which intersects ''S''. Draw {{overline|''OE''}} perpendicular to ''P'' and meeting ''P'' at ''E''. Let ''A'' and ''B'' be any two different points in the intersection. Then ''AOE'' and ''BOE'' are right triangles with a common side, ''OE'', and hypotenuses ''AO'' and ''BO'' equal. Therefore, the remaining sides ''AE'' and ''BE'' are equal. This proves that all points in the intersection are the same distance from the point ''E'' in the plane ''P'', in other words all points in the intersection lie on a circle ''C'' with center ''E''.<ref>Proof follows Hobbs, Prop. 304</ref> This proves that the intersection of ''P'' and ''S'' is contained in ''C''. Note that ''OE'' is the axis of the circle.

Now consider a point ''D'' of the circle ''C''. Since ''C'' lies in ''P'', so does ''D''. On the other hand, the triangles ''AOE'' and ''DOE'' are right triangles with a common side, ''OE'', and legs ''EA'' and ''ED'' equal. Therefore, the hypotenuses ''AO'' and ''DO'' are equal, and equal to the radius of ''S'', so that ''D'' lies in ''S''. This proves that ''C'' is contained in the intersection of ''P'' and ''S''.

As a corollary, on a sphere there is exactly one circle that can be drawn through three given points.<ref>Hobbs, Prop. 308</ref>

The proof can be extended to show that the points on a circle are all a common angular distance from one of its poles.<ref>Hobbs, Prop. 310</ref>

Compare also [[Conic section|conic section]]s, which can produce [[Oval|oval]]s.

=== Two spheres ===
To show that a non-trivial intersection of two spheres is a circle, assume (without loss of generality) that one sphere (with radius <math>R</math>) is centered at the origin. Points on this sphere satisfy
:<math>x^2 + y^2 + z^2 = R^2.</math>
Also without loss of generality, assume that the second sphere, with radius <math>r</math>, is centered at a point on the positive x-axis, at distance <math>a</math> from the origin. Its points satisfy
:<math>(x-a)^2 + y^2 + z^2 = r^2.</math>
The intersection of the spheres is the set of points satisfying both equations. Subtracting the equations gives
:<math>\begin{align}
(x-a)^2 - x^2 & = r^2 - R^2 \\
    a^2 - 2ax & = r^2 - R^2 \\
            x & = \frac{a^2 + R^2 - r^2}{2a}.
\end{align}</math>
In the singular case <math>a = 0</math>, the spheres are concentric. There are two possibilities: if <math>R = r</math>, the spheres coincide, and the intersection is the entire sphere; if <math>R \not= r</math>, the spheres are disjoint and the intersection is empty.
When ''a'' is nonzero, the intersection lies in a vertical plane with this x-coordinate, which may intersect both of the spheres, be tangent to both spheres, or external to both spheres.
The result follows from the previous proof for sphere-plane intersections.

==See also==
* Line-plane intersection
* [[Line–sphere intersection]]
* [[Line-cylinder intersection]]
* [[Analytic geometry]]
* [[Computational geometry]]
* Equation of a line
* [[Intersection (set theory)]]
* [[Intersection theory]]

== Notes ==
{{reflist}}

== References ==
*{{cite book |title=Solid Geometry|url=https://archive.org/details/solidgeometry00hobbgoog|first=C.A.|last=Hobbs|publisher=G.H. Kent|year=1921|pages=[https://archive.org/details/solidgeometry00hobbgoog/page/n151 397] ff}}

==Further reading==
* {{cite web |url= https://www.realtimerendering.com/intersections.html |title=Intersections (Ray Tracing Resources Page) |last= Haines |first= Eric |date= June 6, 2021 |website=Real-Time Rendering |access-date=December 14, 2023 |quote=a grid of intersection routines for various popular objects, pointing to resources in books and on the web.}}
* Nicholas M. Patrikalakis and Takashi Maekawa, ''Shape Interrogation for Computer Aided Design and Manufacturing'', Springer, 2002, {{ISBN|3540424547}}, 9783540424543, pp. 408. [https://books.google.com/books?id=qoxBeti5ySoC&q=intersection]
*{{cite book |title=Solid Geometry
|url=https://archive.org/details/solidgeometry01comsgoog
|first1=M.|last1=Sykes|first2=C.E.|last2=Comstock|publisher=Rand McNally
|year=1922|pages=[https://archive.org/details/solidgeometry01comsgoog/page/n93 81] ff}}


[[Category:Geometric intersection|*]]

{{Sourceattribution|Intersection (geometry)}}