PostGIS  2.3.8dev-r@@SVN_REVISION@@

◆ lw_dist3d_pt_poly()

int lw_dist3d_pt_poly ( POINT3DZ p,
LWPOLY poly,
PLANE3D plane,
POINT3DZ projp,
DISTPTS3D dl 
)

Checking if the point projected on the plane of the polygon actually is inside that polygon.

If so the mindistance is between that projected point and our original point. If not we check from original point to the bounadary. If the projected point is inside a hole of the polygon we check the distance to the boudary of that hole.

Definition at line 1020 of file measures3d.c.

References lw_dist3d_pt_pt(), lw_dist3d_pt_ptarray(), LW_TRUE, LWDEBUG, LWPOLY::nrings, pt_in_ring_3d(), and LWPOLY::rings.

Referenced by lw_dist3d_point_poly(), and lw_dist3d_ptarray_poly().

1021 {
1022  int i;
1023 
1024  LWDEBUG(2, "lw_dist3d_point_poly called");
1025 
1026 
1027  if(pt_in_ring_3d(projp, poly->rings[0], plane))
1028  {
1029  for (i=1; i<poly->nrings; i++)
1030  {
1031  /* Inside a hole. Distance = pt -> ring */
1032  if ( pt_in_ring_3d(projp, poly->rings[i], plane ))
1033  {
1034  LWDEBUG(3, " inside an hole");
1035  return lw_dist3d_pt_ptarray(p, poly->rings[i], dl);
1036  }
1037  }
1038 
1039  return lw_dist3d_pt_pt(p,projp,dl);/* If the projected point is inside the polygon the shortest distance is between that point and the inputed point*/
1040  }
1041  else
1042  {
1043  return lw_dist3d_pt_ptarray(p, poly->rings[0], dl); /*If the projected point is outside the polygon we search for the closest distance against the boundarry instead*/
1044  }
1045 
1046  return LW_TRUE;
1047 
1048 }
#define LWDEBUG(level, msg)
Definition: lwgeom_log.h:83
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:75
int lw_dist3d_pt_pt(POINT3DZ *thep1, POINT3DZ *thep2, DISTPTS3D *dl)
Compares incomming points and stores the points closest to each other or most far away from each othe...
Definition: measures3d.c:812
POINTARRAY ** rings
Definition: liblwgeom.h:456
int nrings
Definition: liblwgeom.h:454
int pt_in_ring_3d(const POINT3DZ *p, const POINTARRAY *ring, PLANE3D *plane)
pt_in_ring_3d(): crossing number test for a point in a polygon input: p = a point, pa = vertex points of a ring V[n+1] with V[n]=V[0] plane=the plane that the vertex points are lying on returns: 0 = outside, 1 = inside
Definition: measures3d.c:1241
int lw_dist3d_pt_ptarray(POINT3DZ *p, POINTARRAY *pa, DISTPTS3D *dl)
search all the segments of pointarray to see which one is closest to p Returns distance between point...
Definition: measures3d.c:717
Here is the call graph for this function:
Here is the caller graph for this function: