test_ptarrayarc_contains_point()

static void test_ptarrayarc_contains_point ( )

static

Definition at line 432 of file cu_ptarray.c.

{
        /* int ptarrayarc_contains_point(const POINTARRAY *pa, const POINT2D *pt) */
 
        LWCIRCSTRING *lwcirc;
        POINTARRAY *pa;
        POINT2D pt;
        int rv;
 
        /*** Collection of semi-circles surrounding unit square ***/
        lwcirc = lwgeom_as_lwcircstring(lwgeom_from_text("CIRCULARSTRING(-1 -1, -2 0, -1 1, 0 2, 1 1, 2 0, 1 -1, 0 -2, -1 -1)"));
        pa = lwcirc->points;
 
        /* Point in middle of square */
        pt.x = 0;
        pt.y = 0;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_INSIDE);
 
        /* Point in left lobe */
        pt.x = -1.1;
        pt.y = 0.1;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_INSIDE);
 
        /* Point on boundary of left lobe */
        pt.x = -1;
        pt.y = 0;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_INSIDE);
 
        /* Point on boundary vertex */
        pt.x = -1;
        pt.y = 1;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_BOUNDARY);
 
        /* Point outside */
        pt.x = -1.5;
        pt.y = 1.5;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_OUTSIDE);
 
        /*** Two-edge ring made up of semi-circles (really, a circle) ***/
        lwcircstring_free(lwcirc);
        lwcirc = lwgeom_as_lwcircstring(lwgeom_from_text("CIRCULARSTRING(-1 0, 0 1, 1 0, 0 -1, -1 0)"));
        pa = lwcirc->points;
 
        /* Point outside */
        pt.x = -1.5;
        pt.y = 1.5;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_OUTSIDE);
 
        /* Point more outside */
        pt.x = 2.5;
        pt.y = 1.5;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_OUTSIDE);
 
        /* Point more outside */
        pt.x = 2.5;
        pt.y = 2.5;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_OUTSIDE);
 
        /* Point inside at middle */
        pt.x = 0;
        pt.y = 0;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_INSIDE);
 
        /* Point inside offset from middle */
        pt.x = 0.01;
        pt.y = 0.01;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_INSIDE);
 
        /* Point on edge vertex */
        pt.x = 0;
        pt.y = 1;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_BOUNDARY);
 
        /*** Two-edge ring, closed ***/
        lwcircstring_free(lwcirc);
        lwcirc = lwgeom_as_lwcircstring(lwgeom_from_text("CIRCULARSTRING(1 6, 6 1, 9 7, 6 10, 1 6)"));
        pa = lwcirc->points;
 
        /* Point to left of ring */
        pt.x = 20;
        pt.y = 4;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_OUTSIDE);
 
        /*** One-edge ring, closed circle ***/
        lwcircstring_free(lwcirc);
        lwcirc = lwgeom_as_lwcircstring(lwgeom_from_text("CIRCULARSTRING(-1 0, 1 0, -1 0)"));
        pa = lwcirc->points;
 
        /* Point inside */
        pt.x = 0;
        pt.y = 0;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_INSIDE);
 
        /* Point outside */
        pt.x = 0;
        pt.y = 2;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_OUTSIDE);
 
        /* Point on boundary */
        pt.x = 0;
        pt.y = 1;
        rv = ptarrayarc_contains_point(pa, &pt);
        CU_ASSERT_EQUAL(rv, LW_BOUNDARY);
 
        /*** Overshort ring ***/
        lwcircstring_free(lwcirc);
        lwcirc = lwgeom_as_lwcircstring(lwgeom_from_text("CIRCULARSTRING(-1 0, 1 0)"));
        pa = lwcirc->points;
        cu_error_msg_reset();
        rv = ptarrayarc_contains_point(pa, &pt);
        //printf("%s\n", cu_error_msg);
        ASSERT_STRING_EQUAL(cu_error_msg, "ptarrayarc_raycast_intersections called with even number of points");
 
        /*** Unclosed ring ***/
        lwcircstring_free(lwcirc);
        lwcirc = lwgeom_as_lwcircstring(lwgeom_from_text("CIRCULARSTRING(-1 0, 1 0, 2 0)"));
        pa = lwcirc->points;
        cu_error_msg_reset();
        rv = ptarrayarc_contains_point(pa, &pt);
        ASSERT_STRING_EQUAL(cu_error_msg, "ptarrayarc_contains_point called on unclosed ring");
 
        lwcircstring_free(lwcirc);
}

References ASSERT_STRING_EQUAL, cu_error_msg, cu_error_msg_reset(), LW_BOUNDARY, LW_INSIDE, LW_OUTSIDE, lwcircstring_free(), lwgeom_as_lwcircstring(), lwgeom_from_text(), LWCIRCSTRING::points, ptarrayarc_contains_point(), POINT2D::x, and POINT2D::y.

Referenced by ptarray_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: