test_raster_envelope_geom()

static void test_raster_envelope_geom ( )

static

Definition at line 88 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_envelope_geom(), rt_raster_get_srid(), rt_raster_new(), rt_raster_set_offsets(), rt_raster_set_scale(), LWPOLY::srid, LWGEOM::type, POINT4D::x, and POINT4D::y.

Referenced by raster_geometry_suite_setup().

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

        /* NULL raster */
        CU_ASSERT_EQUAL(rt_raster_get_envelope_geom(NULL, &env), ES_NONE);
        CU_ASSERT(env == NULL);

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

        CU_ASSERT_EQUAL(rt_raster_get_envelope_geom(raster, &env), ES_NONE);
        CU_ASSERT_EQUAL(env->type, POINTTYPE);
        lwgeom_free(env);
        cu_free_raster(raster);

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

        CU_ASSERT_EQUAL(rt_raster_get_envelope_geom(raster, &env), ES_NONE);
        CU_ASSERT_EQUAL(env->type, LINETYPE);
        lwgeom_free(env);
        cu_free_raster(raster);

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

        CU_ASSERT_EQUAL(rt_raster_get_envelope_geom(raster, &env), ES_NONE);
        CU_ASSERT_EQUAL(env->type, LINETYPE);
        lwgeom_free(env);
        cu_free_raster(raster);

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

        rt_raster_set_offsets(raster, 0.5, 0.5);
        rt_raster_set_scale(raster, 1, -1);

        CU_ASSERT_EQUAL(rt_raster_get_envelope_geom(raster, &env), ES_NONE);
        poly = lwgeom_as_lwpoly(env);
        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, 5.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, 0.5, DBL_EPSILON);

        getPoint4d_p(ring, 2, &pt);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 5.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, -4.5, DBL_EPSILON);

        getPoint4d_p(ring, 3, &pt);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 0.5, DBL_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, -4.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(env);
        cu_free_raster(raster);
}

Here is the call graph for this function:

Here is the caller graph for this function: