PostGIS  3.7.0dev-r@@SVN_REVISION@@

◆ project_point_on_plane()

static double project_point_on_plane ( const POINT3DZ p,
PLANE3D pl,
POINT3DZ p0 
)
static

Finds a point on a plane from where the original point is perpendicular to the plane.

Definition at line 57 of file measures3d.c.

58 {
59  /*In our plane definition we have a point on the plane and a normal vector (pl.pv), perpendicular to the plane
60  this vector will be parallel to the line between our inputted point above the plane and the point we are
61  searching for on the plane. So, we already have a direction from p to find p0, but we don't know the distance.
62  */
63 
64  VECTOR3D v1;
65  double f;
66 
67  if (!get_3dvector_from_points(&(pl->pop), p, &v1))
68  return 0.0;
69 
70  f = DOT(pl->pv, v1);
71  if (FP_IS_ZERO(f))
72  {
73  /* Point is in the plane */
74  *p0 = *p;
75  return 0;
76  }
77 
78  f = -f / DOT(pl->pv, pl->pv);
79 
80  p0->x = p->x + pl->pv.x * f;
81  p0->y = p->y + pl->pv.y * f;
82  p0->z = p->z + pl->pv.z * f;
83 
84  return f;
85 }
#define FP_IS_ZERO(A)
static int get_3dvector_from_points(const POINT3DZ *p1, const POINT3DZ *p2, VECTOR3D *v)
Definition: measures3d.c:34
#define DOT(u, v)
Definition: measures3d.h:31
POINT3DZ pop
Definition: measures3d.h:57
VECTOR3D pv
Definition: measures3d.h:58
double z
Definition: liblwgeom.h:396
double x
Definition: liblwgeom.h:396
double y
Definition: liblwgeom.h:396
double z
Definition: measures3d.h:52
double x
Definition: measures3d.h:52
double y
Definition: measures3d.h:52

References DOT, FP_IS_ZERO, get_3dvector_from_points(), PLANE3D::pop, PLANE3D::pv, POINT3DZ::x, VECTOR3D::x, POINT3DZ::y, VECTOR3D::y, POINT3DZ::z, and VECTOR3D::z.

Referenced by lw_dist3d_point_poly(), lw_dist3d_point_tri(), lw_dist3d_ptarray_poly(), and lw_dist3d_ptarray_tri().

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