PostGIS  2.5.0dev-r@@SVN_REVISION@@
static int _lwt_InitEdgeEndByLine ( edgeend fee,
edgeend lee,
LWLINE edge,
POINT2D fp,
POINT2D lp 
)
static

Definition at line 1461 of file lwgeom_topo.c.

References _lwt_FirstDistinctVertex2D(), azimuth_pt_pt(), edgeend_t::ccwFace, edgeend_t::cwFace, LWDEBUG, LWDEBUGF, lwerror(), edgeend_t::myaz, edgeend_t::nextCCW, edgeend_t::nextCW, POINTARRAY::npoints, LWLINE::points, POINT2D::x, and POINT2D::y.

Referenced by lwt_ChangeEdgeGeom().

1463 {
1464  POINTARRAY *pa = edge->points;
1465  POINT2D pt;
1466 
1467  fee->nextCW = fee->nextCCW =
1468  lee->nextCW = lee->nextCCW = 0;
1469  fee->cwFace = fee->ccwFace =
1470  lee->cwFace = lee->ccwFace = -1;
1471 
1472  /* Compute azimuth of first edge end */
1473  LWDEBUG(1, "computing azimuth of first edge end");
1474  if ( ! _lwt_FirstDistinctVertex2D(pa, fp, 0, 1, &pt) )
1475  {
1476  lwerror("Invalid edge (no two distinct vertices exist)");
1477  return -1;
1478  }
1479  if ( ! azimuth_pt_pt(fp, &pt, &(fee->myaz)) ) {
1480  lwerror("error computing azimuth of first edgeend [%.15g %.15g,%.15g %.15g]",
1481  fp->x, fp->y, pt.x, pt.y);
1482  return -2;
1483  }
1484  LWDEBUGF(1, "azimuth of first edge end [%.15g %.15g,%.15g %.15g] is %g",
1485  fp->x, fp->y, pt.x, pt.y, fee->myaz);
1486 
1487  /* Compute azimuth of second edge end */
1488  LWDEBUG(1, "computing azimuth of second edge end");
1489  if ( ! _lwt_FirstDistinctVertex2D(pa, lp, pa->npoints-1, -1, &pt) )
1490  {
1491  lwerror("Invalid edge (no two distinct vertices exist)");
1492  return -1;
1493  }
1494  if ( ! azimuth_pt_pt(lp, &pt, &(lee->myaz)) ) {
1495  lwerror("error computing azimuth of last edgeend [%.15g %.15g,%.15g %.15g]",
1496  lp->x, lp->y, pt.x, pt.y);
1497  return -2;
1498  }
1499  LWDEBUGF(1, "azimuth of last edge end [%.15g %.15g,%.15g %.15g] is %g",
1500  lp->x, lp->y, pt.x, pt.y, lee->myaz);
1501 
1502  return 0;
1503 }
static int _lwt_FirstDistinctVertex2D(const POINTARRAY *pa, POINT2D *ref, int from, int dir, POINT2D *op)
Definition: lwgeom_topo.c:1422
LWT_ELEMID cwFace
Definition: lwgeom_topo.c:1403
#define LWDEBUG(level, msg)
Definition: lwgeom_log.h:83
LWT_ELEMID nextCCW
Definition: lwgeom_topo.c:1405
double x
Definition: liblwgeom.h:327
double myaz
Definition: lwgeom_topo.c:1409
double y
Definition: liblwgeom.h:327
LWT_ELEMID ccwFace
Definition: lwgeom_topo.c:1407
LWT_ELEMID nextCW
Definition: lwgeom_topo.c:1401
int azimuth_pt_pt(const POINT2D *p1, const POINT2D *p2, double *ret)
Compute the azimuth of segment AB in radians.
Definition: measures.c:2422
#define LWDEBUGF(level, msg,...)
Definition: lwgeom_log.h:88
void lwerror(const char *fmt,...)
Write a notice out to the error handler.
Definition: lwutil.c:190
POINTARRAY * points
Definition: liblwgeom.h:421
uint32_t npoints
Definition: liblwgeom.h:370

Here is the call graph for this function:

Here is the caller graph for this function: