PostGIS  3.6.1dev-r@@SVN_REVISION@@

◆ ptarray_closest_segment_2d()

int ptarray_closest_segment_2d ( const POINTARRAY pa,
const POINT2D qp,
double *  dist 
)
Parameters
pathe subject pointarray
qpthe query point
distoptional output for actual distance from segment
Returns
0-based segment index for the closest segment (earliest segment in case of same distance)

Definition at line 1447 of file ptarray.c.

1448 {
1449  const POINT2D *start = getPoint2d_cp(pa, 0), *end = NULL;
1450  uint32_t t, seg=0;
1451  double mindist=DBL_MAX;
1452 
1453  /* Loop through pointarray looking for nearest segment */
1454  for (t=1; t<pa->npoints; t++)
1455  {
1456  double dist_sqr;
1457  end = getPoint2d_cp(pa, t);
1458  dist_sqr = distance2d_sqr_pt_seg(qp, start, end);
1459 
1460  if (dist_sqr < mindist)
1461  {
1462  mindist = dist_sqr;
1463  seg=t-1;
1464  if ( mindist == 0 )
1465  {
1466  LWDEBUG(3, "Breaking on mindist=0");
1467  break;
1468  }
1469  }
1470 
1471  start = end;
1472  }
1473 
1474  if ( dist ) *dist = sqrt(mindist);
1475  return seg;
1476 }
double distance2d_sqr_pt_seg(const POINT2D *p, const POINT2D *A, const POINT2D *B)
Definition: measures.c:2354
#define LWDEBUG(level, msg)
Definition: lwgeom_log.h:101
static const POINT2D * getPoint2d_cp(const POINTARRAY *pa, uint32_t n)
Returns a POINT2D pointer into the POINTARRAY serialized_ptlist, suitable for reading from.
Definition: lwinline.h:97
uint32_t npoints
Definition: liblwgeom.h:427

References distance2d_sqr_pt_seg(), getPoint2d_cp(), LWDEBUG, and POINTARRAY::npoints.

Referenced by lwt_GetFaceContainingPoint(), and test_ptarray_closest_segment_2d().

Here is the call graph for this function:
Here is the caller graph for this function: