PostGIS  2.5.0dev-r@@SVN_REVISION@@

◆ distance2d_pt_seg()

double distance2d_pt_seg ( const POINT2D p,
const POINT2D A,
const POINT2D B 
)

Definition at line 2330 of file measures.c.

References distance2d_pt_pt(), FP_ABS, r, s, POINT2D::x, and POINT2D::y.

Referenced by lwline_split_by_point_to(), and ptarray_locate_point().

2331 {
2332  double r,s;
2333 
2334  /*if start==end, then use pt distance */
2335  if ( ( A->x == B->x) && (A->y == B->y) )
2336  return distance2d_pt_pt(p,A);
2337 
2338  /*
2339  * otherwise, we use comp.graphics.algorithms
2340  * Frequently Asked Questions method
2341  *
2342  * (1) AC dot AB
2343  * r = ---------
2344  * ||AB||^2
2345  * r has the following meaning:
2346  * r=0 P = A
2347  * r=1 P = B
2348  * r<0 P is on the backward extension of AB
2349  * r>1 P is on the forward extension of AB
2350  * 0<r<1 P is interior to AB
2351  */
2352 
2353  r = ( (p->x-A->x) * (B->x-A->x) + (p->y-A->y) * (B->y-A->y) )/( (B->x-A->x)*(B->x-A->x) +(B->y-A->y)*(B->y-A->y) );
2354 
2355  if (r<0) return distance2d_pt_pt(p,A);
2356  if (r>1) return distance2d_pt_pt(p,B);
2357 
2358 
2359  /*
2360  * (2)
2361  * (Ay-Cy)(Bx-Ax)-(Ax-Cx)(By-Ay)
2362  * s = -----------------------------
2363  * L^2
2364  *
2365  * Then the distance from C to P = |s|*L.
2366  *
2367  */
2368 
2369  s = ( (A->y-p->y)*(B->x-A->x)- (A->x-p->x)*(B->y-A->y) ) /
2370  ( (B->x-A->x)*(B->x-A->x) +(B->y-A->y)*(B->y-A->y) );
2371 
2372  return FP_ABS(s) * sqrt(
2373  (B->x-A->x)*(B->x-A->x) + (B->y-A->y)*(B->y-A->y)
2374  );
2375 }
char * r
Definition: cu_in_wkt.c:24
double x
Definition: liblwgeom.h:330
double y
Definition: liblwgeom.h:330
char * s
Definition: cu_in_wkt.c:23
double distance2d_pt_pt(const POINT2D *p1, const POINT2D *p2)
Definition: measures.c:2312
#define FP_ABS(a)
Here is the call graph for this function:
Here is the caller graph for this function: