test_raster_surface()

static void test_raster_surface ( )

static

Definition at line 255 of file cu_raster_geometry.c.

                                  {
        rt_raster rast;
        rt_band band;
        const int maxX = 5;
        const int maxY = 5;
        int x, y;
        LWMPOLY *mpoly = NULL;
        int err;
        LWGEOM *gexpected, *gobserved;
 
        rast = rt_raster_new(maxX, maxY);
        CU_ASSERT(rast != NULL);
 
        rt_raster_set_offsets(rast, 0, 0);
        rt_raster_set_scale(rast, 1, -1);
 
        band = cu_add_band(rast, PT_32BUI, 1, 0);
        CU_ASSERT(band != NULL);
 
        for (y = 0; y < maxY; y++) {
                for (x = 0; x < maxX; x++) {
                        rt_band_set_pixel(band, x, y, 1, NULL);
                }
        }
 
        err = rt_raster_surface(rast, 0, &mpoly);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(mpoly != NULL);
        //wkt = lwgeom_to_text(lwmpoly_as_lwgeom(mpoly));
        gobserved = (LWGEOM *)lwmpoly_as_lwgeom(mpoly);
        gexpected = lwgeom_from_wkt("MULTIPOLYGON(((0 0,0 -5,5 -5,5 0,0 0)))", LW_PARSER_CHECK_NONE);
        CU_ASSERT_DOUBLE_EQUAL(lwgeom_area(gobserved), lwgeom_area(gexpected), FLT_EPSILON);
        lwmpoly_free(mpoly);
        mpoly = NULL;
 
        /* 0,0 NODATA */
        rt_band_set_pixel(band, 0, 0, 0, NULL);
 
        err = rt_raster_surface(rast, 0, &mpoly);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(mpoly != NULL);
        gobserved = (LWGEOM *)lwmpoly_as_lwgeom(mpoly);
        gexpected = lwgeom_from_wkt("MULTIPOLYGON(((1 0,1 -1,0 -1,0 -5,4 -5,5 -5,5 0,1 0)))", LW_PARSER_CHECK_NONE);
        CU_ASSERT_DOUBLE_EQUAL(lwgeom_area(gobserved), lwgeom_area(gexpected), FLT_EPSILON);
        lwmpoly_free(mpoly);
        mpoly = NULL;
 
        /* plus 1,1 NODATA */
        rt_band_set_pixel(band, 1, 1, 0, NULL);
 
        err = rt_raster_surface(rast, 0, &mpoly);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(mpoly != NULL);
        gobserved = (LWGEOM *)lwmpoly_as_lwgeom(mpoly);
        gexpected = lwgeom_from_wkt("MULTIPOLYGON(((1 0,1 -1,0 -1,0 -5,4 -5,5 -5,5 0,1 0),(1 -1,1 -2,2 -2,2 -1,1 -1)))",
                                    LW_PARSER_CHECK_NONE);
        CU_ASSERT_DOUBLE_EQUAL(lwgeom_area(gobserved), lwgeom_area(gexpected), FLT_EPSILON);
        lwmpoly_free(mpoly);
        mpoly = NULL;
 
        /* plus 2,2 NODATA */
        rt_band_set_pixel(band, 2, 2, 0, NULL);
 
        err = rt_raster_surface(rast, 0, &mpoly);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(mpoly != NULL);
        gobserved = (LWGEOM *)lwmpoly_as_lwgeom(mpoly);
        gexpected = lwgeom_from_wkt(
            "MULTIPOLYGON(((1 -1,1 0,5 0,5 -5,4 -5,0 -5,0 -1,1 -1),(1 -1,1 -2,2 -2,2 -1,1 -1),(2 -2,2 -3,3 -3,3 -2,2 -2)))",
            LW_PARSER_CHECK_NONE);
        CU_ASSERT_DOUBLE_EQUAL(lwgeom_perimeter(gobserved), lwgeom_perimeter(gexpected), FLT_EPSILON);
        lwmpoly_free(mpoly);
        mpoly = NULL;
 
        /* plus 3,3 NODATA */
        rt_band_set_pixel(band, 3, 3, 0, NULL);
 
        err = rt_raster_surface(rast, 0, &mpoly);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(mpoly != NULL);
        gobserved = (LWGEOM *)lwmpoly_as_lwgeom(mpoly);
        gexpected = lwgeom_from_wkt(
            "MULTIPOLYGON(((1 -1,1 0,5 0,5 -5,4 -5,0 -5,0 -1,1 -1),(1 -1,1 -2,2 -2,2 -1,1 -1),(2 -2,2 -3,3 -3,3 -2,2 -2),(3 -3,3 -4,4 -4,4 -3,3 -3)))",
            LW_PARSER_CHECK_NONE);
        CU_ASSERT_DOUBLE_EQUAL(lwgeom_perimeter(gobserved), lwgeom_perimeter(gexpected), FLT_EPSILON);
        lwmpoly_free(mpoly);
        mpoly = NULL;
 
        /* plus 4,4 NODATA */
        rt_band_set_pixel(band, 4, 4, 0, NULL);
 
        err = rt_raster_surface(rast, 0, &mpoly);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(mpoly != NULL);
        gobserved = (LWGEOM *)lwmpoly_as_lwgeom(mpoly);
        gexpected = lwgeom_from_wkt(
            "MULTIPOLYGON(((4 -4,4 -5,0 -5,0 -1,1 -1,1 -2,2 -2,2 -3,3 -3,3 -4,4 -4)),((1 -1,1 0,5 0,5 -4,4 -4,4 -3,3 -3,3 -2,2 -2,2 -1,1 -1)))",
            LW_PARSER_CHECK_NONE);
        CU_ASSERT_DOUBLE_EQUAL(lwgeom_area(gobserved), lwgeom_area(gexpected), FLT_EPSILON);
        lwmpoly_free(mpoly);
        mpoly = NULL;
 
        /* a whole lot more NODATA */
        rt_band_set_pixel(band, 4, 0, 0, NULL);
        rt_band_set_pixel(band, 3, 1, 0, NULL);
        rt_band_set_pixel(band, 1, 3, 0, NULL);
        rt_band_set_pixel(band, 0, 4, 0, NULL);
 
        err = rt_raster_surface(rast, 0, &mpoly);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(mpoly != NULL);
        gobserved = (LWGEOM *)lwmpoly_as_lwgeom(mpoly);
        gexpected = lwgeom_from_wkt(
            "MULTIPOLYGON(((1 -4,2 -4,2 -3,3 -3,3 -4,4 -4,4 -5,3 -5,1 -5,1 -4)),((1 -4,0 -4,0 -1,1 -1,1 -2,2 -2,2 -3,1 -3,1 -4)),((3 -2,4 -2,4 -1,5 -1,5 -4,4 -4,4 -3,3 -3,3 -2)),((3 -2,2 -2,2 -1,1 -1,1 0,4 0,4 -1,3 -1,3 -2)))",
            LW_PARSER_CHECK_NONE);
        CU_ASSERT_DOUBLE_EQUAL(lwgeom_area(gobserved), lwgeom_area(gexpected), FLT_EPSILON);
        lwmpoly_free(mpoly);
        mpoly = NULL;
 
        cu_free_raster(rast);
}

References cu_add_band(), cu_free_raster(), ES_NONE, LW_PARSER_CHECK_NONE, lwgeom_area(), lwgeom_from_wkt(), lwgeom_perimeter(), lwmpoly_as_lwgeom(), lwmpoly_free(), PT_32BUI, rt_band_set_pixel(), rt_raster_new(), rt_raster_set_offsets(), rt_raster_set_scale(), and rt_raster_surface().

Referenced by raster_geometry_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: