Distance from a point to a line

Short description: Problem in coordinate geometry

In Euclidean geometry, the distance from a point to a line is the shortest distance from a given point to any point on an infinite straight line. It is the perpendicular distance of the point to the line, the length of the line segment which joins the point to nearest point on the line. The algebraic expression for calculating it can be derived and expressed in several ways.

Knowing the distance from a point to a line can be useful in various situations—for example, finding the shortest distance to reach a road, quantifying the scatter on a graph, etc. In Deming regression, a type of linear curve fitting, if the dependent and independent variables have equal variance this results in orthogonal regression in which the degree of imperfection of the fit is measured for each data point as the perpendicular distance of the point from the regression line.

Line defined by an equation

In the case of a line in the plane given by the equation ax + by + c = 0, where a, b and c are real constants with a and b not both zero, the distance from the line to a point (x₀, y₀) is

[math]\displaystyle{ \operatorname{distance}(ax+by+c=0, (x_0, y_0)) = \frac{|ax_0+by_0+c|}{\sqrt{a^2+b^2}}. }[/math]

The point on this line which is closest to (x₀, y₀) has coordinates:

[math]\displaystyle{ x = \frac{b(bx_0 - ay_0)-ac}{a^2 + b^2} \text{ and } y = \frac{a(-bx_0 + ay_0) - bc}{a^2+b^2}. }[/math]

Horizontal and vertical lines

In the general equation of a line, ax + by + c = 0, a and b cannot both be zero unless c is also zero, in which case the equation does not define a line. If a = 0 and b ≠ 0, the line is horizontal and has equation y = −c/b. The distance from (x₀, y₀) to this line is measured along a vertical line segment of length |y₀ − (−c/b)| = |by₀ + c|/|b| in accordance with the formula. Similarly, for vertical lines (b = 0) the distance between the same point and the line is |ax₀ + c|/|a|, as measured along a horizontal line segment.

Line defined by two points

If the line passes through two points P₁ = (x₁, y₁) and P₂ = (x₂, y₂) then the distance of (x₀, y₀) from the line is:

[math]\displaystyle{ \operatorname{distance}(P_1, P_2, (x_0, y_0)) = \frac{|(x_2-x_1)(y_1-y_0)-(x_1-x_0)(y_2-y_1)|}{\sqrt{(x_2-x_1)^2+(y_2-y_1)^2}}. }[/math]

The denominator of this expression is the distance between P₁ and P₂. The numerator is twice the area of the triangle with its vertices at the three points, (x₀, y₀), P₁ and P₂. See: Area of a triangle § Using coordinates. The expression is equivalent to h = 2A/b, which can be obtained by rearranging the standard formula for the area of a triangle: A = 1/2 bh, where b is the length of a side, and h is the perpendicular height from the opposite vertex.

Line defined by point and angle

If the line passes through the point P = (P_x, P_y) with angle θ, then the distance of some point (x₀, y₀) to the line is

[math]\displaystyle{ \operatorname{distance}(P, \theta, (x_0, y_0)) = |\cos(\theta)(P_y-y_0) -\sin(\theta)(P_x-x_0)| }[/math]

Proofs

An algebraic proof

This proof is valid only if the line is neither vertical nor horizontal, that is, we assume that neither a nor b in the equation of the line is zero.

The line with equation ax + by + c = 0 has slope −a/b, so any line perpendicular to it will have slope b/a (the negative reciprocal). Let (m, n) be the point of intersection of the line ax + by + c = 0 and the line perpendicular to it which passes through the point (x₀, y₀). The line through these two points is perpendicular to the original line, so

[math]\displaystyle{ \frac{y_0 - n}{x_0 - m}=\frac{b}{a}. }[/math]

Thus, [math]\displaystyle{ a(y_0 -n) - b(x_0 - m) = 0, }[/math] and by squaring this equation we obtain:

[math]\displaystyle{ a^2(y_0 - n)^2 + b^2(x_0 - m)^2 = 2ab(y_0 - n)(x_0 - m). }[/math]

Now consider,

[math]\displaystyle{ \begin{align} (a(x_0 - m) + b(y_0 - n))^2 & = a^2(x_0 - m)^2 + 2ab(y_0 -n)(x_0 - m) + b^2(y_0 - n)^2 \\ & = \left(a^2 + b^2\right) \left((x_0 - m)^2 + (y_0 - n)^2\right) \end{align} }[/math]

using the above squared equation. But we also have,

[math]\displaystyle{ (a(x_0 - m) + b(y_0 - n))^2 = (ax_0 + by_0 - am - bn)^2 = (ax_0 + by_0 + c)^2 }[/math]

since (m, n) is on ax + by + c = 0. Thus,

[math]\displaystyle{ \left(a^2 + b^2\right) \left((x_0 - m)^2 + (y_0 - n)^2\right) = (ax_0 + by_0 + c)^2 }[/math]

and we obtain the length of the line segment determined by these two points,

[math]\displaystyle{ d=\sqrt{(x_0 - m)^2+(y_0 - n)^2} = \frac{|ax_0+ by_0 +c|}{\sqrt{a^2+b^2}}. }[/math]

A geometric proof

Diagram for geometric proof

This proof is valid only if the line is not horizontal or vertical.

Drop a perpendicular from the point P with coordinates (x₀, y₀) to the line with equation Ax + By + C = 0. Label the foot of the perpendicular R. Draw the vertical line through P and label its intersection with the given line S. At any point T on the line, draw a right triangle TVU whose sides are horizontal and vertical line segments with hypotenuse TU on the given line and horizontal side of length |B| (see diagram). The vertical side of ∆TVU will have length |A| since the line has slope -A/B.

∆PRS and ∆TVU are similar triangles, since they are both right triangles and ∠PSR ≅ ∠TUV since they are corresponding angles of a transversal to the parallel lines PS and UV (both are vertical lines). Corresponding sides of these triangles are in the same ratio, so:

[math]\displaystyle{ \frac{|\overline{PR}|}{|\overline{PS}|} = \frac{|\overline{TV}|}{|\overline{TU}|}. }[/math]

If point S has coordinates (x₀,m) then |PS| = |y₀ - m| and the distance from P to the line is:

[math]\displaystyle{ |\overline{PR} | = \frac{|y_0 - m||B|}{\sqrt{A^2 + B^2}}. }[/math]

Since S is on the line, we can find the value of m,

[math]\displaystyle{ m = \frac{-Ax_0 - C}{B}, }[/math]

and finally obtain:

[math]\displaystyle{ |\overline{PR}| = \frac{|Ax_0 + By_0 + C|}{\sqrt{A^2 + B^2}}. }[/math]

A variation of this proof is to place V at P and compute the area of the triangle ∆UVT two ways to obtain that [math]\displaystyle{ D|\overline{TU}| = |\overline{VU}||\overline{VT}| }[/math] where D is the altitude of ∆UVT drawn to the hypotenuse of ∆UVT from P. The distance formula can then used to express [math]\displaystyle{ |\overline{TU}| }[/math], [math]\displaystyle{ |\overline{VU}| }[/math], and [math]\displaystyle{ |\overline{VT}| }[/math]in terms of the coordinates of P and the coefficients of the equation of the line to get the indicated formula.

A vector projection proof

Let P be the point with coordinates (x₀, y₀) and let the given line have equation ax + by + c = 0. Also, let Q = (x₁, y₁) be any point on this line and n the vector (a, b) starting at point Q. The vector n is perpendicular to the line, and the distance d from point P to the line is equal to the length of the orthogonal projection of [math]\displaystyle{ \overrightarrow{QP} }[/math] on n. The length of this projection is given by:

[math]\displaystyle{ d = \frac{|\overrightarrow{QP} \cdot \mathbf{n}|}{\| \mathbf{n}\|}. }[/math]

Now,

[math]\displaystyle{ \overrightarrow{QP} = (x_0 - x_1, y_0 - y_1), }[/math] so [math]\displaystyle{ \overrightarrow{QP} \cdot \mathbf{n} = a(x_0 - x_1) + b(y_0 - y_1) }[/math] and [math]\displaystyle{ \| \mathbf{n} \| = \sqrt{a^2 + b^2}, }[/math]

thus

[math]\displaystyle{ d = \frac{|a(x_0 - x_1) + b(y_0 - y_1)|}{\sqrt{a^2 + b^2}}. }[/math]

Since Q is a point on the line, [math]\displaystyle{ c = -ax_1 - by_1 }[/math], and so,

[math]\displaystyle{ d = \frac{|ax_0 + by_0 + c|}{\sqrt{a^2 + b^2}}. }[/math]

Although the distance is given as a modulus, the sign can be useful to determine which side of the line the point is on, in a sense determined by the direction of normal vector (a,b)

Another formula

It is possible to produce another expression to find the shortest distance of a point to a line. This derivation also requires that the line be not vertical or horizontal.

The point P is given with coordinates ([math]\displaystyle{ x_0, y_0 }[/math]). The equation of a line is given by [math]\displaystyle{ y=mx+k }[/math]. The equation of the normal of that line which passes through the point P is given [math]\displaystyle{ y=\frac{x_0-x}{m}+y_0 }[/math].

The point at which these two lines intersect is the closest point on the original line to the point P. Hence:

[math]\displaystyle{ mx+k=\frac{x_0-x}{m}+y_0. }[/math]

We can solve this equation for x,

[math]\displaystyle{ x=\frac{x_0+my_0-mk}{m^2+1}. }[/math]

The y coordinate of the point of intersection can be found by substituting this value of x into the equation of the original line,

[math]\displaystyle{ y=m\frac{(x_0+my_0-mk)}{m^2+1}+k. }[/math]

Using the equation for finding the distance between 2 points, [math]\displaystyle{ d=\sqrt{(X_2-X_1)^2+(Y_2-Y_1)^2} }[/math], we can deduce that the formula to find the shortest distance between a line and a point is the following:

[math]\displaystyle{ d=\sqrt{ \left( {\frac{x_0 + m y_0-mk}{m^2+1}-x_0 } \right) ^2 + \left( {m\frac{x_0+m y_0-mk}{m^2+1}+k-y_0 }\right) ^2 } = \frac{|k + m x_0 - y_0|}\sqrt{1 + m^2} . }[/math]

Recalling that m = -a/b and k = - c/b for the line with equation ax + by + c = 0, a little algebraic simplification reduces this to the standard expression.

Vector formulation

Illustration of the vector formulation.

The equation of a line can be given in vector form:

[math]\displaystyle{ \mathbf{x} = \mathbf{a} + t\mathbf{n} }[/math]

Here a is a point on the line, and n is a unit vector in the direction of the line. Then as scalar t varies, x gives the locus of the line.

The distance of an arbitrary point p to this line is given by

[math]\displaystyle{ \operatorname{distance}(\mathbf{x} = \mathbf{a} + t\mathbf{n}, \mathbf{p}) = \| (\mathbf{p}-\mathbf{a}) - ((\mathbf{p}-\mathbf{a}) \cdot \mathbf{n})\mathbf{n} \|. }[/math]

This formula can be derived as follows: [math]\displaystyle{ \mathbf{p}-\mathbf{a} }[/math] is a vector from a to the point p. Then [math]\displaystyle{ (\mathbf{p} - \mathbf{a}) \cdot \mathbf{n} }[/math] is the projected length onto the line and so

[math]\displaystyle{ \mathbf{a} + ((\mathbf{p} - \mathbf{a}) \cdot \mathbf{n})\mathbf{n} }[/math]

is a vector that is the projection of [math]\displaystyle{ \mathbf{p}-\mathbf{a} }[/math] onto the line and represents the point on the line closest to [math]\displaystyle{ \mathbf{p} }[/math]. Thus

[math]\displaystyle{ (\mathbf{p}-\mathbf{a}) - ((\mathbf{p}-\mathbf{a}) \cdot \mathbf{n})\mathbf{n} }[/math]

is the component of [math]\displaystyle{ \mathbf{p}-\mathbf{a} }[/math] perpendicular to the line. The distance from the point to the line is then just the norm of that vector. This more general formula is not restricted to two dimensions.

Another vector formulation

If the vector space is orthonormal and if the line goes through point a and has a direction vector n, the distance between point p and the line is

[math]\displaystyle{ \operatorname{distance}(\mathbf{x} = \mathbf{a} + t\mathbf{n}, \mathbf{p}) = \frac{\left\|(\mathbf{p}-\mathbf{a}) \times \mathbf{n}\right\|}{\|\mathbf{n}\|}. }[/math]

Note that cross products only exist in dimensions 3 and 7.

See also

• Hesse normal form
• Line-line intersection
• Distance between two lines
• Distance from a point to a plane
• Skew lines

References

• Anton, Howard (1994), Elementary Linear Algebra (7th ed.), John Wiley & Sons, ISBN 0-471-58742-7 
• Ballantine, J.P.; Jerbert, A.R. (1952), "Distance from a line or plane to a point", American Mathematical Monthly 59 (4): 242–243, doi:10.2307/2306514 
• Larson, Ron; Hostetler, Robert (2007), Precalculus: A Concise Course, Houghton Mifflin Co., ISBN 978-0-618-62719-6, https://archive.org/details/precalculusconci00lars 
• Spain, Barry (2007), Analytical Conics, Dover Publications, ISBN 978-0-486-45773-4 
• Weisstein, Eric W.. "Point-Line Distance--3-Dimensional". MathWorld. https://mathworld.wolfram.com/Point-LineDistance3-Dimensional.html. 

Further reading

• Deza, Michel Marie; Deza, Elena (2013), Encyclopedia of Distances (2nd ed.), Springer, p. 86, ISBN 9783642309588, https://books.google.com/books?id=QxX2CX5OVMsC&pg=PA86