PostGIS  3.4.0dev-r@@SVN_REVISION@@

◆ lw_dist2d_tri_poly()

int lw_dist2d_tri_poly ( LWTRIANGLE tri,
LWPOLY poly,
DISTPTS dl 
)

Definition at line 835 of file measures.c.

836 {
837  POINTARRAY *pa = tri->points;
838  const POINT2D *pt = getPoint2d_cp(pa, 0);
839 
840  /* If we are looking for maxdistance, just check the outer rings.*/
841  if (dl->mode == DIST_MAX)
842  return lw_dist2d_ptarray_ptarray(pa, poly->rings[0], dl);
843 
844  /* Triangle has a point outside poly. Check distance to outer ring only. */
845  if (ptarray_contains_point(poly->rings[0], pt) == LW_OUTSIDE)
846  {
847  if (!lw_dist2d_ptarray_ptarray(pa, poly->rings[0], dl))
848  return LW_FALSE;
849 
850  /* just a check if the answer is already given */
851  if (dl->distance <= dl->tolerance)
852  return LW_TRUE;
853 
854  /* Maybe poly is inside triangle? */
855  const POINT2D *pt2 = getPoint2d_cp(poly->rings[0], 0);
856  if (ptarray_contains_point(pa, pt2) != LW_OUTSIDE)
857  {
858  dl->distance = 0.0;
859  dl->p1.x = dl->p2.x = pt2->x;
860  dl->p1.y = dl->p2.y = pt2->y;
861  return LW_TRUE;
862  }
863  }
864 
865  for (uint32_t i = 1; i < poly->nrings; i++)
866  {
867  if (!lw_dist2d_ptarray_ptarray(pa, poly->rings[i], dl))
868  return LW_FALSE;
869 
870  /* just a check if the answer is already given */
871  if (dl->distance <= dl->tolerance && dl->mode == DIST_MIN)
872  return LW_TRUE;
873  }
874 
875  /* It's inside a hole, then the actual distance is the min ring distance */
876  for (uint32_t i = 1; i < poly->nrings; i++)
877  if (ptarray_contains_point(poly->rings[i], pt) != LW_OUTSIDE)
878  return LW_TRUE;
879 
880  /* Not in hole, so inside polygon */
881  dl->distance = 0.0;
882  dl->p1.x = dl->p2.x = pt->x;
883  dl->p1.y = dl->p2.y = pt->y;
884  return LW_TRUE;
885 }
#define LW_FALSE
Definition: liblwgeom.h:94
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:93
int ptarray_contains_point(const POINTARRAY *pa, const POINT2D *pt)
Return 1 if the point is inside the POINTARRAY, -1 if it is outside, and 0 if it is on the boundary.
Definition: ptarray.c:753
#define LW_OUTSIDE
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:101
int lw_dist2d_ptarray_ptarray(POINTARRAY *l1, POINTARRAY *l2, DISTPTS *dl)
test each segment of l1 against each segment of l2.
Definition: measures.c:1209
#define DIST_MIN
Definition: measures.h:44
#define DIST_MAX
Definition: measures.h:43
POINT2D p1
Definition: measures.h:52
POINT2D p2
Definition: measures.h:53
double tolerance
Definition: measures.h:56
int mode
Definition: measures.h:54
double distance
Definition: measures.h:51
POINTARRAY ** rings
Definition: liblwgeom.h:519
uint32_t nrings
Definition: liblwgeom.h:524
POINTARRAY * points
Definition: liblwgeom.h:495
double y
Definition: liblwgeom.h:390
double x
Definition: liblwgeom.h:390

References DIST_MAX, DIST_MIN, DISTPTS::distance, getPoint2d_cp(), lw_dist2d_ptarray_ptarray(), LW_FALSE, LW_OUTSIDE, LW_TRUE, DISTPTS::mode, LWPOLY::nrings, DISTPTS::p1, DISTPTS::p2, LWTRIANGLE::points, ptarray_contains_point(), LWPOLY::rings, DISTPTS::tolerance, POINT2D::x, and POINT2D::y.

Referenced by lw_dist2d_distribute_bruteforce().

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