test_raster_convex_hull()

static void test_raster_convex_hull ( )

static

Definition at line 166 of file cu_raster_geometry.c.

References cu_free_raster(), ES_NONE, getPoint4d_p(), LINETYPE, lwgeom_as_lwpoly(), lwgeom_free(), POINTARRAY::npoints, LWPOLY::nrings, POINTTYPE, rtrowdump::raster, LWPOLY::rings, rt_raster_get_convex_hull(), rt_raster_get_srid(), rt_raster_new(), rt_raster_set_offsets(), rt_raster_set_scale(), rt_raster_set_skews(), LWPOLY::srid, LWGEOM::type, POINT4D::x, and POINT4D::y.

Referenced by raster_geometry_suite_setup().

                                      {
        rt_raster raster = NULL;
        LWGEOM *hull = NULL;
        LWPOLY *poly = NULL;
        POINTARRAY *ring = NULL;
        POINT4D pt;

        /* NULL raster */
        CU_ASSERT_EQUAL(rt_raster_get_convex_hull(NULL, &hull), ES_NONE);
        CU_ASSERT(hull == NULL);

        /* width = 0, height = 0 */
        raster = rt_raster_new(0, 0);
        CU_ASSERT(raster != NULL);

        CU_ASSERT_EQUAL(rt_raster_get_convex_hull(raster, &hull), ES_NONE);
        CU_ASSERT_EQUAL(hull->type, POINTTYPE);
        lwgeom_free(hull);
        cu_free_raster(raster);

        /* width = 0 */
        raster = rt_raster_new(0, 256);
        CU_ASSERT(raster != NULL);

        CU_ASSERT_EQUAL(rt_raster_get_convex_hull(raster, &hull), ES_NONE);
        CU_ASSERT_EQUAL(hull->type, LINETYPE);
        lwgeom_free(hull);
        cu_free_raster(raster);

        /* height = 0 */
        raster = rt_raster_new(256, 0);
        CU_ASSERT(raster != NULL);

        CU_ASSERT_EQUAL(rt_raster_get_convex_hull(raster, &hull), ES_NONE);
        CU_ASSERT_EQUAL(hull->type, LINETYPE);
        lwgeom_free(hull);
        cu_free_raster(raster);

        /* normal raster */
        raster = rt_raster_new(256, 256);
        CU_ASSERT(raster != NULL);

        rt_raster_set_offsets(raster, 0.5, 0.5);
        rt_raster_set_scale(raster, 1, 1);
        rt_raster_set_skews(raster, 4, 5);

        CU_ASSERT_EQUAL(rt_raster_get_convex_hull(raster, &hull), ES_NONE);
        poly = lwgeom_as_lwpoly(hull);
        CU_ASSERT_EQUAL(poly->srid, rt_raster_get_srid(raster));
        CU_ASSERT_EQUAL(poly->nrings, 1);

        ring = poly->rings[0];
        CU_ASSERT(ring != NULL);
        CU_ASSERT_EQUAL(ring->npoints, 5);

        getPoint4d_p(ring, 0, &pt);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 0.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, 0.5, DBL_EPSILON);

        getPoint4d_p(ring, 1, &pt);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 256.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, 1280.5, DBL_EPSILON);

        getPoint4d_p(ring, 2, &pt);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 1280.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, 1536.5, DBL_EPSILON);

        getPoint4d_p(ring, 3, &pt);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 1024.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, 256.5, DBL_EPSILON);

        getPoint4d_p(ring, 4, &pt);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 0.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, 0.5, DBL_EPSILON);

        lwgeom_free(hull);
        cu_free_raster(raster);
}

Here is the call graph for this function:

Here is the caller graph for this function: