test_raster_perimeter()

static void test_raster_perimeter ( )

static

Definition at line 377 of file cu_raster_geometry.c.

                                    {
        rt_raster rast;
        rt_band band;
        const int maxX = 5;
        const int maxY = 5;
        int x, y;
        char *wkt = NULL;
        LWGEOM *geom = NULL;
        int err;
 
        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_get_perimeter(rast, -1, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((0 0,5 0,5 -5,0 -5,0 0))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        /* row 0 is NODATA */
        rt_band_set_pixel(band, 0, 0, 0, NULL);
        rt_band_set_pixel(band, 1, 0, 0, NULL);
        rt_band_set_pixel(band, 2, 0, 0, NULL);
        rt_band_set_pixel(band, 3, 0, 0, NULL);
        rt_band_set_pixel(band, 4, 0, 0, NULL);
 
        err = rt_raster_get_perimeter(rast, 0, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((0 -1,5 -1,5 -5,0 -5,0 -1))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        /* column 4 is NODATA */
        /* pixel 4, 0 already set to NODATA */
        rt_band_set_pixel(band, 4, 1, 0, NULL);
        rt_band_set_pixel(band, 4, 2, 0, NULL);
        rt_band_set_pixel(band, 4, 3, 0, NULL);
        rt_band_set_pixel(band, 4, 4, 0, NULL);
 
        err = rt_raster_get_perimeter(rast, 0, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((0 -1,4 -1,4 -5,0 -5,0 -1))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        /* row 4 is NODATA */
        rt_band_set_pixel(band, 0, 4, 0, NULL);
        rt_band_set_pixel(band, 1, 4, 0, NULL);
        rt_band_set_pixel(band, 2, 4, 0, NULL);
        rt_band_set_pixel(band, 3, 4, 0, NULL);
        /* pixel 4, 4 already set to NODATA */
 
        err = rt_raster_get_perimeter(rast, 0, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((0 -1,4 -1,4 -4,0 -4,0 -1))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        /* column 0 is NODATA */
        /* pixel 0, 0 already set to NODATA*/
        rt_band_set_pixel(band, 0, 1, 0, NULL);
        rt_band_set_pixel(band, 0, 2, 0, NULL);
        rt_band_set_pixel(band, 0, 3, 0, NULL);
        /* pixel 0, 4 already set to NODATA*/
 
        err = rt_raster_get_perimeter(rast, 0, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((1 -1,4 -1,4 -4,1 -4,1 -1))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        /* columns 1 and 3 are NODATA */
        /* pixel 1, 0 already set to NODATA */
        rt_band_set_pixel(band, 1, 1, 0, NULL);
        rt_band_set_pixel(band, 1, 2, 0, NULL);
        rt_band_set_pixel(band, 1, 3, 0, NULL);
        /* pixel 1, 4 already set to NODATA */
        /* pixel 3, 0 already set to NODATA */
        rt_band_set_pixel(band, 3, 1, 0, NULL);
        rt_band_set_pixel(band, 3, 2, 0, NULL);
        rt_band_set_pixel(band, 3, 3, 0, NULL);
        /* pixel 3, 4 already set to NODATA */
 
        err = rt_raster_get_perimeter(rast, 0, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((2 -1,3 -1,3 -4,2 -4,2 -1))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        /* more pixels are NODATA */
        rt_band_set_pixel(band, 2, 1, 0, NULL);
        rt_band_set_pixel(band, 2, 3, 0, NULL);
 
        err = rt_raster_get_perimeter(rast, 0, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((2 -2,3 -2,3 -3,2 -3,2 -2))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        /* second band */
        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_get_perimeter(rast, -1, &geom);
        CU_ASSERT_EQUAL(err, ES_NONE);
        CU_ASSERT(geom != NULL);
        wkt = lwgeom_to_text(geom);
        CU_ASSERT_STRING_EQUAL(wkt, "POLYGON((0 0,5 0,5 -5,0 -5,0 0))");
        rtdealloc(wkt);
        lwgeom_free(geom);
        geom = NULL;
 
        cu_free_raster(rast);
}

References cu_add_band(), cu_free_raster(), ES_NONE, lwgeom_free(), lwgeom_to_text(), PT_32BUI, rt_band_set_pixel(), rt_raster_get_perimeter(), rt_raster_new(), rt_raster_set_offsets(), rt_raster_set_scale(), and rtdealloc().

Referenced by raster_geometry_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: