PostGIS  3.0.6dev-r@@SVN_REVISION@@

◆ test_lwgeom_tcpa()

static void test_lwgeom_tcpa ( void  )
static

Definition at line 1251 of file cu_measures.c.

1252 {
1253  LWGEOM *g1, *g2;
1254  double m, dist;
1255 
1256  /* Invalid input, lack of dimensions */
1257 
1258  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 0 0 1)", LW_PARSER_CHECK_NONE);
1259  g2 = lwgeom_from_wkt("LINESTRING (0 0 2, 0 0 5)", LW_PARSER_CHECK_NONE);
1260  m = lwgeom_tcpa(g1, g2, NULL);
1261  lwgeom_free(g1);
1262  lwgeom_free(g2);
1263  ASSERT_DOUBLE_EQUAL(m, -1.0);
1264  ASSERT_STRING_EQUAL(
1265  cu_error_msg,
1266  "Both input geometries must have a measure dimension");
1267 
1268  /* Invalid input, not linestrings */
1269 
1270  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 0 0 1)", LW_PARSER_CHECK_NONE);
1271  g2 = lwgeom_from_wkt("POINT M (0 0 2)", LW_PARSER_CHECK_NONE);
1272  m = lwgeom_tcpa(g1, g2, NULL);
1273  lwgeom_free(g1);
1274  lwgeom_free(g2);
1275  ASSERT_DOUBLE_EQUAL(m, -1.0);
1276  ASSERT_STRING_EQUAL(cu_error_msg,
1277  "Both input geometries must be linestrings");
1278 
1279  /* Invalid input, too short linestring */
1280 
1281  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 10 0 10)", LW_PARSER_CHECK_NONE);
1282  g2 = lwgeom_from_wkt("LINESTRING M(2 0 1)", LW_PARSER_CHECK_NONE);
1283  dist = -77;
1284  m = lwgeom_tcpa(g1, g2, &dist);
1285  lwgeom_free(g1);
1286  lwgeom_free(g2);
1287  ASSERT_DOUBLE_EQUAL(dist, -77.0); /* not touched */
1288  ASSERT_DOUBLE_EQUAL(m, -1.0);
1289  ASSERT_STRING_EQUAL(
1290  cu_error_msg,
1291  "Both input lines must have at least 2 points" /* should be accepted
1292  ? */
1293 
1294  );
1295 
1296  /* Invalid input, empty linestring */
1297 
1298  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 10 0 10)", LW_PARSER_CHECK_NONE);
1299  g2 = lwgeom_from_wkt("LINESTRING M EMPTY", LW_PARSER_CHECK_NONE);
1300  m = lwgeom_tcpa(g1, g2, NULL);
1301  lwgeom_free(g1);
1302  lwgeom_free(g2);
1303  ASSERT_DOUBLE_EQUAL(m, -1.0);
1304  ASSERT_STRING_EQUAL(cu_error_msg,
1305  "Both input lines must have at least 2 points"
1306 
1307  );
1308 
1309  /* Timeranges do not overlap */
1310 
1311  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 0 0 1)", LW_PARSER_CHECK_NONE);
1312  g2 = lwgeom_from_wkt("LINESTRING M(0 0 2, 0 0 5)", LW_PARSER_CHECK_NONE);
1313  m = lwgeom_tcpa(g1, g2, NULL);
1314  lwgeom_free(g1);
1315  lwgeom_free(g2);
1316  ASSERT_DOUBLE_EQUAL(m, -2.0); /* means timeranges do not overlap */
1317 
1318  /* One of the tracks is still, the other passes to that point */
1319 
1320  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 0 0 1)", LW_PARSER_CHECK_NONE);
1321  g2 = lwgeom_from_wkt("LINESTRING M(0 0 1, 0 0 5)", LW_PARSER_CHECK_NONE);
1322  dist = -1;
1323  m = lwgeom_tcpa(g1, g2, &dist);
1324  ASSERT_DOUBLE_EQUAL(m, 1.0);
1325  ASSERT_DOUBLE_EQUAL(dist, 0.0);
1326  CU_ASSERT( lwgeom_cpa_within(g1, g2, 0.0) == LW_TRUE );
1327  lwgeom_free(g1);
1328  lwgeom_free(g2);
1329 
1330  /* One of the tracks is still, the other passes at 10 meters from point */
1331 
1332  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 0 0 5)", LW_PARSER_CHECK_NONE);
1333  g2 = lwgeom_from_wkt("LINESTRING M(-10 10 1, 10 10 5)", LW_PARSER_CHECK_NONE);
1334  dist = -1;
1335  m = lwgeom_tcpa(g1, g2, &dist);
1336  ASSERT_DOUBLE_EQUAL(m, 3.0);
1337  ASSERT_DOUBLE_EQUAL(dist, 10.0);
1338  CU_ASSERT( lwgeom_cpa_within(g1, g2, 11.0) == LW_TRUE );
1339  CU_ASSERT( lwgeom_cpa_within(g1, g2, 10.0) == LW_TRUE );
1340  CU_ASSERT( lwgeom_cpa_within(g1, g2, 9.0) == LW_FALSE );
1341  lwgeom_free(g1);
1342  lwgeom_free(g2);
1343 
1344  /* Equal tracks, 2d */
1345 
1346  g1 = lwgeom_from_wkt("LINESTRING M(0 0 10, 10 0 20)", LW_PARSER_CHECK_NONE);
1347  g2 = lwgeom_from_wkt("LINESTRING M(0 0 10, 10 0 20)", LW_PARSER_CHECK_NONE);
1348  dist = -1;
1349  m = lwgeom_tcpa(g1, g2, &dist);
1350  lwgeom_free(g1);
1351  lwgeom_free(g2);
1352  ASSERT_DOUBLE_EQUAL(m, 10.0);
1353  ASSERT_DOUBLE_EQUAL(dist, 0.0);
1354 
1355  /* Reversed tracks, 2d */
1356 
1357  g1 = lwgeom_from_wkt("LINESTRING M(0 0 10, 10 0 20)", LW_PARSER_CHECK_NONE);
1358  g2 = lwgeom_from_wkt("LINESTRING M(10 0 10, 0 0 20)", LW_PARSER_CHECK_NONE);
1359  dist = -1;
1360  m = lwgeom_tcpa(g1, g2, &dist);
1361  lwgeom_free(g1);
1362  lwgeom_free(g2);
1363  ASSERT_DOUBLE_EQUAL(m, 15.0);
1364  ASSERT_DOUBLE_EQUAL(dist, 0.0);
1365 
1366  /* Parallel tracks, same speed, 2d */
1367 
1368  g1 = lwgeom_from_wkt("LINESTRING M(2 0 10, 12 0 20)", LW_PARSER_CHECK_NONE);
1369  g2 = lwgeom_from_wkt("LINESTRING M(13 0 10, 23 0 20)", LW_PARSER_CHECK_NONE);
1370  dist = -1;
1371  m = lwgeom_tcpa(g1, g2, &dist);
1372  lwgeom_free(g1);
1373  lwgeom_free(g2);
1374  ASSERT_DOUBLE_EQUAL(m, 10.0);
1375  ASSERT_DOUBLE_EQUAL(dist, 11.0);
1376 
1377  /* Parallel tracks, different speed (g2 gets closer as time passes), 2d */
1378 
1379  g1 = lwgeom_from_wkt("LINESTRING M(4 0 10, 10 0 20)", LW_PARSER_CHECK_NONE);
1380  g2 = lwgeom_from_wkt("LINESTRING M(2 0 10, 9 0 20)", LW_PARSER_CHECK_NONE);
1381  dist = -1;
1382  m = lwgeom_tcpa(g1, g2, &dist);
1383  lwgeom_free(g1);
1384  lwgeom_free(g2);
1385  ASSERT_DOUBLE_EQUAL(m, 20.0);
1386  ASSERT_DOUBLE_EQUAL(dist, 1.0);
1387 
1388  /* Parallel tracks, different speed (g2 left behind as time passes), 2d */
1389 
1390  g1 = lwgeom_from_wkt("LINESTRING M(4 0 10, 10 0 20)", LW_PARSER_CHECK_NONE);
1391  g2 = lwgeom_from_wkt("LINESTRING M(2 0 10, 6 0 20)", LW_PARSER_CHECK_NONE);
1392  dist = -1;
1393  m = lwgeom_tcpa(g1, g2, &dist);
1394  lwgeom_free(g1);
1395  lwgeom_free(g2);
1396  ASSERT_DOUBLE_EQUAL(m, 10.0);
1397  ASSERT_DOUBLE_EQUAL(dist, 2.0);
1398 
1399  /* Tracks, colliding, 2d */
1400 
1401  g1 = lwgeom_from_wkt("LINESTRING M(0 0 0, 10 0 10)", LW_PARSER_CHECK_NONE);
1402  g2 = lwgeom_from_wkt("LINESTRING M(5 -8 0, 5 8 10)", LW_PARSER_CHECK_NONE);
1403  dist = -1;
1404  m = lwgeom_tcpa(g1, g2, &dist);
1405  lwgeom_free(g1);
1406  lwgeom_free(g2);
1407  ASSERT_DOUBLE_EQUAL(m, 5.0);
1408  ASSERT_DOUBLE_EQUAL(dist, 0.0);
1409 
1410  /* Tracks crossing, NOT colliding, 2d */
1411 
1412  g1 = lwgeom_from_wkt("LINESTRING M(0 0 0, 10 0 10)", LW_PARSER_CHECK_NONE);
1413  g2 = lwgeom_from_wkt("LINESTRING M(8 -5 0, 8 5 10)", LW_PARSER_CHECK_NONE);
1414  dist = -1;
1415  m = lwgeom_tcpa(g1, g2, &dist);
1416  lwgeom_free(g1);
1417  lwgeom_free(g2);
1418  ASSERT_DOUBLE_EQUAL(m, 6.5);
1419  ASSERT_DOUBLE_EQUAL(rint(dist*100), 212.0);
1420 
1421  /* Same origin, different direction, 2d */
1422 
1423  g1 = lwgeom_from_wkt("LINESTRING M(0 0 1, 10 0 10)", LW_PARSER_CHECK_NONE);
1424  g2 = lwgeom_from_wkt("LINESTRING M(0 0 1, -100 0 10)", LW_PARSER_CHECK_NONE);
1425  dist = -1;
1426  m = lwgeom_tcpa(g1, g2, &dist);
1427  lwgeom_free(g1);
1428  lwgeom_free(g2);
1429  ASSERT_DOUBLE_EQUAL(m, 1.0);
1430  ASSERT_DOUBLE_EQUAL(dist, 0.0);
1431 
1432  /* Same ending, different direction, 2d */
1433 
1434  g1 = lwgeom_from_wkt("LINESTRING M(10 0 1, 0 0 10)", LW_PARSER_CHECK_NONE);
1435  g2 = lwgeom_from_wkt("LINESTRING M(0 -100 1, 0 0 10)", LW_PARSER_CHECK_NONE);
1436  dist = -1;
1437  m = lwgeom_tcpa(g1, g2, &dist);
1438  lwgeom_free(g1);
1439  lwgeom_free(g2);
1440  ASSERT_DOUBLE_EQUAL(m, 10.0);
1441  ASSERT_DOUBLE_EQUAL(dist, 0.0);
1442 
1443  /* Converging tracks, 3d */
1444 
1445  g1 = lwgeom_from_wkt("LINESTRING ZM(0 0 0 10, 10 0 0 20)", LW_PARSER_CHECK_NONE);
1446  g2 = lwgeom_from_wkt("LINESTRING ZM(0 0 8 10, 10 0 5 20)", LW_PARSER_CHECK_NONE);
1447  dist = -1;
1448  m = lwgeom_tcpa(g1, g2, &dist);
1449  lwgeom_free(g1);
1450  lwgeom_free(g2);
1451  ASSERT_DOUBLE_EQUAL(m, 20.0);
1452  ASSERT_DOUBLE_EQUAL(dist, 5.0);
1453 
1454  /* G1 stops at t=1 until t=4 to let G2 pass by, then continues */
1455  /* G2 passes at 1 meter from G1 t=3 */
1456 
1457  g1 = lwgeom_from_wkt("LINESTRING M(0 0 0, 0 1 1, 0 1 4, 0 10 13)", LW_PARSER_CHECK_NONE);
1458  g2 = lwgeom_from_wkt("LINESTRING M(-10 2 0, 0 2 3, 12 2 13)", LW_PARSER_CHECK_NONE);
1459  dist = -1;
1460  m = lwgeom_tcpa(g1, g2, &dist);
1461  lwgeom_free(g1);
1462  lwgeom_free(g2);
1463  ASSERT_DOUBLE_EQUAL(m, 3.0);
1464  ASSERT_DOUBLE_EQUAL(dist, 1.0);
1465 
1466  /* Test for https://trac.osgeo.org/postgis/ticket/3136 */
1467 
1468  g1 = lwgeom_from_wkt("LINESTRING M (0 0 1432291464,2 0 1432291536) ", LW_PARSER_CHECK_NONE);
1469  g2 = lwgeom_from_wkt("LINESTRING M (0 0 1432291464,1 0 1432291466.25,2 0 1432291500)", LW_PARSER_CHECK_NONE);
1470  dist = -1;
1471  m = lwgeom_tcpa(g1, g2, &dist);
1472  lwgeom_free(g1);
1473  lwgeom_free(g2);
1474  ASSERT_DOUBLE_EQUAL(m, 1432291464.0);
1475  ASSERT_DOUBLE_EQUAL(dist, 0.0);
1476 
1477  /* Tracks share a single point in time */
1478 
1479  g1 = lwgeom_from_wkt("LINESTRINGM(0 0 0, 1 0 2)", LW_PARSER_CHECK_NONE);
1480  g2 = lwgeom_from_wkt("LINESTRINGM(0 0 2, 1 0 3)", LW_PARSER_CHECK_NONE);
1481  dist = -1;
1482  m = lwgeom_tcpa(g1, g2, &dist);
1483  lwgeom_free(g1);
1484  lwgeom_free(g2);
1485  ASSERT_DOUBLE_EQUAL(m, 2.0);
1486  ASSERT_DOUBLE_EQUAL(dist, 1.0);
1487 
1488 }
cu_error_msg
char cu_error_msg[MAX_CUNIT_ERROR_LENGTH+1]
Definition: liblwgeom/cunit/cu_tester.c:19
ASSERT_DOUBLE_EQUAL
#define ASSERT_DOUBLE_EQUAL(o, e)
Definition: liblwgeom/cunit/cu_tester.h:31
ASSERT_STRING_EQUAL
#define ASSERT_STRING_EQUAL(o, e)
Definition: liblwgeom/cunit/cu_tester.h:43
LW_FALSE
#define LW_FALSE
Definition: liblwgeom.h:108
lwgeom_free
void lwgeom_free(LWGEOM *geom)
Definition: lwgeom.c:1138
LW_PARSER_CHECK_NONE
#define LW_PARSER_CHECK_NONE
Definition: liblwgeom.h:2060
lwgeom_tcpa
double lwgeom_tcpa(const LWGEOM *g1, const LWGEOM *g2, double *mindist)
Find the time of closest point of approach.
Definition: lwlinearreferencing.c:1144
lwgeom_from_wkt
LWGEOM * lwgeom_from_wkt(const char *wkt, const char check)
Definition: lwin_wkt.c:905
LW_TRUE
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:107
lwgeom_cpa_within
int lwgeom_cpa_within(const LWGEOM *g1, const LWGEOM *g2, double maxdist)
Is the closest point of approach within a distance ?
Definition: lwlinearreferencing.c:1311

References ASSERT_DOUBLE_EQUAL, ASSERT_STRING_EQUAL, cu_error_msg, LW_FALSE, LW_PARSER_CHECK_NONE, LW_TRUE, lwgeom_cpa_within(), lwgeom_free(), lwgeom_from_wkt(), and lwgeom_tcpa().

Referenced by measures_suite_setup().

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