test_ptarrayarc_contains_point()

static void test_ptarrayarc_contains_point ( )

static

Definition at line 552 of file cu_ptarray.c.

References ASSERT_STRING_EQUAL, cu_error_msg, cu_error_msg_reset(), LW_BOUNDARY, LW_INSIDE, LW_OUTSIDE, lwgeom_as_lwline(), lwgeom_from_text(), lwline_free(), LWLINE::points, ptarrayarc_contains_point(), POINT2D::x, and POINT2D::y.

Referenced by ptarray_suite_setup().

{
        /* int ptarrayarc_contains_point(const POINTARRAY *pa, const POINT2D *pt) */

        LWLINE *lwline;
        POINTARRAY *pa;
        POINT2D pt;
        int rv;

        /*** Collection of semi-circles surrounding unit square ***/
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(-1 -1, -2 0, -1 1, 0 2, 1 1, 2 0, 1 -1, 0 -2, -1 -1)"));
        pa = lwline->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) ***/
        lwline_free(lwline);
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(-1 0, 0 1, 1 0, 0 -1, -1 0)"));
        pa = lwline->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 ***/
        lwline_free(lwline);
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(1 6, 6 1, 9 7, 6 10, 1 6)"));
        pa = lwline->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 ***/
        lwline_free(lwline);
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(-1 0, 1 0, -1 0)"));
        pa = lwline->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 ***/
        lwline_free(lwline);
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(-1 0, 1 0)"));
        pa = lwline->points;
        cu_error_msg_reset();
        rv = ptarrayarc_contains_point(pa, &pt);
        //printf("%s\n", cu_error_msg);
        ASSERT_STRING_EQUAL("ptarrayarc_contains_point called with even number of points", cu_error_msg);

        /*** Unclosed ring ***/
        lwline_free(lwline);
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(-1 0, 1 0, 2 0)"));
        pa = lwline->points;
        cu_error_msg_reset();
        rv = ptarrayarc_contains_point(pa, &pt);
        ASSERT_STRING_EQUAL("ptarrayarc_contains_point called on unclosed ring", cu_error_msg);

        lwline_free(lwline);
}

Here is the call graph for this function:

Here is the caller graph for this function: