test_edge_intersection()

static void test_edge_intersection ( void  )

static

Definition at line 401 of file cu_geodetic.c.

References edge_intersection(), edge_set(), GEOGRAPHIC_EDGE::end, GEOGRAPHIC_POINT::lat, GEOGRAPHIC_POINT::lon, LW_FALSE, LW_TRUE, point_rad2deg(), and GEOGRAPHIC_EDGE::start.

Referenced by geodetic_suite_setup().

{
        GEOGRAPHIC_EDGE e1, e2;
        GEOGRAPHIC_POINT g;
        int rv;

        /* Covers case, end-to-end intersection */
        edge_set(50, -10.999999999999998224, -10.0, 50.0, &e1);
        edge_set(-10.0, 50.0, -10.272779983831613393, -16.937003313332997578, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /* Medford case, very short segment vs very long one */
        e1.start.lat = 0.74123572595649878103;
        e1.start.lon = -2.1496353191142714145;
        e1.end.lat = 0.74123631950116664058;
        e1.end.lon = -2.1496353248304860273;
        e2.start.lat = 0.73856343764436815924;
        e2.start.lon = -2.1461493501950630325;
        e2.end.lat = 0.70971354024834598651;
        e2.end.lon = 2.1082194552519770703;
        rv = edge_intersection(&e1, &e2, &g);
        CU_ASSERT_EQUAL(rv, LW_FALSE);

        /* Again, this time with a less exact input edge. */
        edge_set(-123.165031277506, 42.4696787216231, -123.165031605021, 42.4697127292275, &e1);
        rv = edge_intersection(&e1, &e2, &g);
        CU_ASSERT_EQUAL(rv, LW_FALSE);

        /* Second Medford case, very short segment vs very long one
        e1.start.lat = 0.73826546728290887156;
        e1.start.lon = -2.14426380171833042;
        e1.end.lat = 0.73826545883786642843;
        e1.end.lon = -2.1442638997530165668;
        e2.start.lat = 0.73775469118192538165;
        e2.start.lon = -2.1436035534281718817;
        e2.end.lat = 0.71021099548296817705;
        e2.end.lon = 2.1065275171200439353;
        rv = edge_intersection(e1, e2, &g);
        CU_ASSERT_EQUAL(rv, LW_FALSE);
        */

        /* Intersection at (0 0) */
        edge_set(-1.0, 0.0, 1.0, 0.0, &e1);
        edge_set(0.0, -1.0, 0.0, 1.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
        CU_ASSERT_DOUBLE_EQUAL(g.lat, 0.0, 0.00001);
        CU_ASSERT_DOUBLE_EQUAL(g.lon, 0.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /*  No intersection at (0 0) */
        edge_set(-1.0, 0.0, 1.0, 0.0, &e1);
        edge_set(0.0, -1.0, 0.0, -2.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        CU_ASSERT_EQUAL(rv, LW_FALSE);

        /*  End touches middle of segment at (0 0) */
        edge_set(-1.0, 0.0, 1.0, 0.0, &e1);
        edge_set(0.0, -1.0, 0.0, 0.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
#if 0
        printf("\n");
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e1.start.lon,  e1.start.lat, e1.end.lon,  e1.end.lat);
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e2.start.lon,  e2.start.lat, e2.end.lon,  e2.end.lat);
        printf("g = (%.15g %.15g)\n", g.lon, g.lat);
        printf("rv = %d\n", rv);
#endif
        CU_ASSERT_DOUBLE_EQUAL(g.lon, 0.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /*  End touches end of segment at (0 0) */
        edge_set(0.0, 0.0, 1.0, 0.0, &e1);
        edge_set(0.0, -1.0, 0.0, 0.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
#if 0
        printf("\n");
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e1.start.lon,  e1.start.lat, e1.end.lon,  e1.end.lat);
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e2.start.lon,  e2.start.lat, e2.end.lon,  e2.end.lat);
        printf("g = (%.15g %.15g)\n", g.lon, g.lat);
        printf("rv = %d\n", rv);
#endif
        CU_ASSERT_DOUBLE_EQUAL(g.lat, 0.0, 0.00001);
        CU_ASSERT_DOUBLE_EQUAL(g.lon, 0.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /* Intersection at (180 0) */
        edge_set(-179.0, -1.0, 179.0, 1.0, &e1);
        edge_set(-179.0, 1.0, 179.0, -1.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
        CU_ASSERT_DOUBLE_EQUAL(g.lat, 0.0, 0.00001);
        CU_ASSERT_DOUBLE_EQUAL(fabs(g.lon), 180.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /* Intersection at (180 0) */
        edge_set(-170.0, 0.0, 170.0, 0.0, &e1);
        edge_set(180.0, -10.0, 180.0, 10.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
        CU_ASSERT_DOUBLE_EQUAL(g.lat, 0.0, 0.00001);
        CU_ASSERT_DOUBLE_EQUAL(fabs(g.lon), 180.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /* Intersection at north pole */
        edge_set(-180.0, 80.0, 0.0, 80.0, &e1);
        edge_set(90.0, 80.0, -90.0, 80.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
        CU_ASSERT_DOUBLE_EQUAL(g.lat, 90.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /* Equal edges return true */
        edge_set(45.0, 10.0, 50.0, 20.0, &e1);
        edge_set(45.0, 10.0, 50.0, 20.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /* Parallel edges (same great circle, different end points) return true  */
        edge_set(40.0, 0.0, 70.0, 0.0, &e1);
        edge_set(60.0, 0.0, 50.0, 0.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
        CU_ASSERT_EQUAL(rv, 2); /* Hack, returning 2 as the 'co-linear' value */

        /* End touches arc at north pole */
        edge_set(-180.0, 80.0, 0.0, 80.0, &e1);
        edge_set(90.0, 80.0, -90.0, 90.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
#if 0
        printf("\n");
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e1.start.lon,  e1.start.lat, e1.end.lon,  e1.end.lat);
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e2.start.lon,  e2.start.lat, e2.end.lon,  e2.end.lat);
        printf("g = (%.15g %.15g)\n", g.lon, g.lat);
        printf("rv = %d\n", rv);
#endif
        CU_ASSERT_DOUBLE_EQUAL(g.lat, 90.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

        /* End touches end at north pole */
        edge_set(-180.0, 80.0, 0.0, 90.0, &e1);
        edge_set(90.0, 80.0, -90.0, 90.0, &e2);
        rv = edge_intersection(&e1, &e2, &g);
        point_rad2deg(&g);
#if 0
        printf("\n");
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e1.start.lon,  e1.start.lat, e1.end.lon,  e1.end.lat);
        printf("LINESTRING(%.15g %.15g, %.15g %.15g)\n", e2.start.lon,  e2.start.lat, e2.end.lon,  e2.end.lat);
        printf("g = (%.15g %.15g)\n", g.lon, g.lat);
        printf("rv = %d\n", rv);
#endif
        CU_ASSERT_DOUBLE_EQUAL(g.lat, 90.0, 0.00001);
        CU_ASSERT_EQUAL(rv, LW_TRUE);

}

Here is the call graph for this function:

Here is the caller graph for this function: