PostGIS  3.7.0dev-r@@SVN_REVISION@@

◆ lw_dist2d_ptarray_ptarrayarc()

int lw_dist2d_ptarray_ptarrayarc ( const POINTARRAY pa,
const POINTARRAY pb,
DISTPTS dl 
)

Test each segment of pa against each arc of pb for distance.

Definition at line 1268 of file measures.c.

1269 {
1270  uint32_t t, u;
1271  const POINT2D *A1;
1272  const POINT2D *A2;
1273  const POINT2D *B1;
1274  const POINT2D *B2;
1275  const POINT2D *B3;
1276  int twist = dl->twisted;
1277 
1278  LWDEBUGF(2, "lw_dist2d_ptarray_ptarrayarc called (points: %d-%d)", pa->npoints, pb->npoints);
1279 
1280  if (pb->npoints % 2 == 0 || pb->npoints < 3)
1281  {
1282  lwerror("lw_dist2d_ptarray_ptarrayarc called with non-arc input");
1283  return LW_FALSE;
1284  }
1285 
1286  if (dl->mode == DIST_MAX)
1287  {
1288  lwerror("lw_dist2d_ptarray_ptarrayarc does not currently support DIST_MAX mode");
1289  return LW_FALSE;
1290  }
1291  else
1292  {
1293  A1 = getPoint2d_cp(pa, 0);
1294  for (t = 1; t < pa->npoints; t++) /* For each segment in pa */
1295  {
1296  A2 = getPoint2d_cp(pa, t);
1297  B1 = getPoint2d_cp(pb, 0);
1298  for (u = 1; u < pb->npoints; u += 2) /* For each arc in pb */
1299  {
1300  B2 = getPoint2d_cp(pb, u);
1301  B3 = getPoint2d_cp(pb, u + 1);
1302  dl->twisted = twist;
1303 
1304  lw_dist2d_seg_arc(A1, A2, B1, B2, B3, dl);
1305 
1306  /* If we've found a distance within tolerance, we're done */
1307  if (dl->distance <= dl->tolerance && dl->mode == DIST_MIN)
1308  return LW_TRUE;
1309 
1310  B1 = B3;
1311  }
1312  A1 = A2;
1313  }
1314  }
1315  return LW_TRUE;
1316 }
#define LW_FALSE
Definition: liblwgeom.h:94
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:93
#define LWDEBUGF(level, msg,...)
Definition: lwgeom_log.h:106
void void lwerror(const char *fmt,...) __attribute__((format(printf
Write a notice out to the error handler.
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:97
int lw_dist2d_seg_arc(const POINT2D *A1, const POINT2D *A2, const POINT2D *B1, const POINT2D *B2, const POINT2D *B3, DISTPTS *dl)
Calculate the shortest distance between an arc and an edge.
Definition: measures.c:1373
#define DIST_MIN
Definition: measures.h:44
#define DIST_MAX
Definition: measures.h:43
int twisted
Definition: measures.h:55
double tolerance
Definition: measures.h:56
int mode
Definition: measures.h:54
double distance
Definition: measures.h:51
uint32_t npoints
Definition: liblwgeom.h:427

References DIST_MAX, DIST_MIN, DISTPTS::distance, getPoint2d_cp(), lw_dist2d_seg_arc(), LW_FALSE, LW_TRUE, LWDEBUGF, lwerror(), DISTPTS::mode, POINTARRAY::npoints, DISTPTS::tolerance, and DISTPTS::twisted.

Referenced by lw_dist2d_line_circstring(), lw_dist2d_tri_circstring(), and test_lw_dist2d_ptarray_ptarrayarc().

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