PostGIS 3.0.6dev-r@@SVN_REVISION@@
Loading...
Searching...
No Matches

◆ lw_dist2d_line_poly()

int lw_dist2d_line_poly ( LWLINE line,
LWPOLY poly,
DISTPTS dl 
)

line to polygon calculation Brute force.

Test line-ring distance against each ring. If there's an intersection (distance==0) then return 0 (crosses boundary). Otherwise, test to see if any point is inside outer rings of polygon, but not in inner rings. If so, return 0 (line inside polygon), otherwise return min distance to a ring (could be outside polygon or inside a hole)

Definition at line 739 of file measures.c.

740{
741 POINTARRAY *pa = line->points;
742 const POINT2D *pt = getPoint2d_cp(pa, 0);
743
744 /* Line has a pount outside poly. Check distance to outer ring only. */
745 if (ptarray_contains_point(poly->rings[0], pt) == LW_OUTSIDE || dl->mode == DIST_MAX)
746 return lw_dist2d_ptarray_ptarray(pa, poly->rings[0], dl);
747
748 for (uint32_t i = 1; i < poly->nrings; i++)
749 {
750 if (!lw_dist2d_ptarray_ptarray(pa, poly->rings[i], dl))
751 return LW_FALSE;
752
753 /* just a check if the answer is already given */
754 if (dl->distance <= dl->tolerance && dl->mode == DIST_MIN)
755 return LW_TRUE;
756 }
757
758 /* It's inside a hole, then the actual distance is the min ring distance */
759 for (uint32_t i = 1; i < poly->nrings; i++)
760 if (ptarray_contains_point(poly->rings[i], pt) != LW_OUTSIDE)
761 return LW_TRUE;
762
763 /* Not in hole, so inside polygon */
764 if (dl->mode == DIST_MIN)
765 {
766 dl->distance = 0.0;
767 dl->p1.x = dl->p2.x = pt->x;
768 dl->p1.y = dl->p2.y = pt->y;
769 }
770 return LW_TRUE;
771}
#define LW_FALSE
Definition liblwgeom.h:108
#define LW_TRUE
Return types for functions with status returns.
Definition liblwgeom.h:107
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:740
#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:91
int lw_dist2d_ptarray_ptarray(POINTARRAY *l1, POINTARRAY *l2, DISTPTS *dl)
test each segment of l1 against each segment of l2.
Definition measures.c:1208
#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 * points
Definition liblwgeom.h:469
POINTARRAY ** rings
Definition liblwgeom.h:505
uint32_t nrings
Definition liblwgeom.h:510
double y
Definition liblwgeom.h:376
double x
Definition liblwgeom.h:376

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, LWLINE::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: