test_gdal_polygonize()

static void test_gdal_polygonize ( )

static

Definition at line 154 of file cu_gdal.c.

                                   {
        int i;
        rt_raster rt;
        int nPols = 0;
        double total_area = 0;
        double total_val = 0;
        rt_geomval gv = NULL;
        LWGEOM *gobserved;
        //char *wkt = NULL;
 
        rt = fillRasterToPolygonize(1, -1.0);
        CU_ASSERT(rt_raster_has_band(rt, 0));
 
        nPols = 0;
        gv = rt_raster_gdal_polygonize(rt, 0, TRUE, &nPols);
        CU_ASSERT_DOUBLE_EQUAL(nPols, 4, FLT_EPSILON);
                total_area = 0; total_val = 0;
        for (i = 0; i < nPols; i++) {
                total_val += gv[i].val;
                gobserved = (LWGEOM *) gv[i].geom;
                total_area += lwgeom_area(gobserved);
                lwgeom_free((LWGEOM *) gv[i].geom);
        }
        printf("total area, total val, nPols  = %f, %f, %i\n", total_area, total_val, nPols);
        CU_ASSERT_DOUBLE_EQUAL(total_val, 1.8 + 0.0 + 2.8 + 0, FLT_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(total_area, 81, FLT_EPSILON);
 
        rtdealloc(gv);
        cu_free_raster(rt);
 
        /* Second test: NODATA value = 1.8 */
        rt = fillRasterToPolygonize(1, 1.8);
 
        /* We can check rt_raster_has_band here too */
        CU_ASSERT(rt_raster_has_band(rt, 0));
 
        nPols = 0;
        gv = rt_raster_gdal_polygonize(rt, 0, TRUE, &nPols);
        CU_ASSERT_DOUBLE_EQUAL(nPols, 4, FLT_EPSILON);
        total_area = 0; total_val = 0;
        for (i = 0; i < nPols; i++) {
                total_val += gv[i].val;
                gobserved = (LWGEOM *) gv[i].geom;
                total_area += lwgeom_area(gobserved);
                lwgeom_free((LWGEOM *) gv[i].geom);
        }
        printf("total area, total_val, polys = %f, %f, %i\n", total_area, total_val, nPols);
        CU_ASSERT_DOUBLE_EQUAL(total_val, 4.6, FLT_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(total_area, 81, FLT_EPSILON);
 
 
        rtdealloc(gv);
        cu_free_raster(rt);
 
        /* Third test: NODATA value = 2.8 */
        rt = fillRasterToPolygonize(1, 2.8);
 
        /* We can check rt_raster_has_band here too */
        CU_ASSERT(rt_raster_has_band(rt, 0));
 
        nPols = 0;
        gv = rt_raster_gdal_polygonize(rt, 0, TRUE, &nPols);
        CU_ASSERT_DOUBLE_EQUAL(nPols, 4, FLT_EPSILON);
        total_area = 0; total_val = 0;
        for (i = 0; i < nPols; i++) {
                total_val += gv[i].val;
                gobserved = (LWGEOM *) gv[i].geom;
                total_area += lwgeom_area(gobserved);
                lwgeom_free((LWGEOM *) gv[i].geom);
        }
 
        printf("total area, total_val, polys = %f, %f, %i\n", total_area, total_val, nPols);
        CU_ASSERT_DOUBLE_EQUAL(total_val, 4.6, FLT_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(total_area, 81, FLT_EPSILON);
 
        rtdealloc(gv);
        cu_free_raster(rt);
 
        /* Fourth test: NODATA value = 0 */
        rt = fillRasterToPolygonize(1, 0.0);
        /* We can check rt_raster_has_band here too */
        CU_ASSERT(rt_raster_has_band(rt, 0));
 
        nPols = 0;
        gv = rt_raster_gdal_polygonize(rt, 0, TRUE, &nPols);
 
        CU_ASSERT_DOUBLE_EQUAL(nPols, 2, FLT_EPSILON);
        total_area = 0; total_val = 0;
        for (i = 0; i < nPols; i++) {
                total_val += gv[i].val;
                gobserved = (LWGEOM *) gv[i].geom;
                total_area += lwgeom_area(gobserved);
                lwgeom_free((LWGEOM *) gv[i].geom);
        }
 
        printf("total area, total_val, polys = %f, %f, %i\n", total_area, total_val, nPols);
        CU_ASSERT_DOUBLE_EQUAL(total_val, 4.6, FLT_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(total_area, 28, FLT_EPSILON);
 
        rtdealloc(gv);
        cu_free_raster(rt);
 
        /* Last test: There is no NODATA value (all values are valid) */
        rt = fillRasterToPolygonize(0, 0.0);
        /* We can check rt_raster_has_band here too */
        CU_ASSERT(rt_raster_has_band(rt, 0));
 
        nPols = 0;
        gv = rt_raster_gdal_polygonize(rt, 0, TRUE, &nPols);
 
        CU_ASSERT_DOUBLE_EQUAL(nPols, 4, FLT_EPSILON);
        total_area = 0; total_val = 0;
        for (i = 0; i < nPols; i++) {
                total_val += gv[i].val;
                gobserved = (LWGEOM *) gv[i].geom;
                total_area += lwgeom_area(gobserved);
                lwgeom_free((LWGEOM *) gv[i].geom);
        }
 
        printf("total area, total_val, polys = %f, %f, %i\n", total_area, total_val, nPols);
        CU_ASSERT_DOUBLE_EQUAL(total_val, 1.8 + 0.0 + 2.8 + 0.0, FLT_EPSILON);
        CU_ASSERT_DOUBLE_EQUAL(total_area, 81, FLT_EPSILON);
        rtdealloc(gv);
        cu_free_raster(rt);
}

References cu_free_raster(), fillRasterToPolygonize(), lwgeom_area(), lwgeom_free(), rt_raster_gdal_polygonize(), rt_raster_has_band(), rtdealloc(), TRUE, and rt_geomval_t::val.

Referenced by gdal_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: