test_band_reclass()

static void test_band_reclass ( )

static

Definition at line 796 of file cu_mapalgebra.c.

References ovdump::band, cu_add_band(), cu_free_raster(), rt_reclassexpr_t::dst, ES_NONE, rt_reclassexpr_t::rt_reclassrange::exc_max, rt_reclassexpr_t::rt_reclassrange::exc_min, rt_reclassexpr_t::rt_reclassrange::inc_max, rt_reclassexpr_t::rt_reclassrange::inc_min, rt_reclassexpr_t::rt_reclassrange::max, rt_reclassexpr_t::rt_reclassrange::min, PT_16BUI, PT_8BUI, rtrowdump::raster, rt_band_destroy(), rt_band_get_nodata(), rt_band_get_pixel(), rt_band_reclass(), rt_band_set_nodata(), rt_band_set_pixel(), rt_raster_new(), rtalloc(), rtdealloc(), rt_reclassexpr_t::src, pixval::x, and pixval::y.

Referenced by mapalgebra_suite_setup().

                                {
        rt_reclassexpr *exprset;

        rt_raster raster;
        rt_band band;
        uint16_t x;
        uint16_t y;
        double nodata;
        int cnt = 2;
        int i = 0;
        int rtn;
        rt_band newband;
        double val;

        raster = rt_raster_new(100, 10);
        CU_ASSERT(raster != NULL); /* or we're out of virtual memory */
        band = cu_add_band(raster, PT_16BUI, 0, 0);
        CU_ASSERT(band != NULL);
        rt_band_set_nodata(band, 0, NULL);

        for (x = 0; x < 100; x++) {
                for (y = 0; y < 10; y++) {
                        rtn = rt_band_set_pixel(band, x, y, x * y + (x + y), NULL);
                }
        }

        rt_band_get_nodata(band, &nodata);
        CU_ASSERT_DOUBLE_EQUAL(nodata, 0, DBL_EPSILON);

        exprset = rtalloc(cnt * sizeof(rt_reclassexpr));
        CU_ASSERT(exprset != NULL);

        for (i = 0; i < cnt; i++) {
                exprset[i] = rtalloc(sizeof(struct rt_reclassexpr_t));
                CU_ASSERT(exprset[i] != NULL);

                if (i == 0) {
                        /* nodata */
                        exprset[i]->src.min = 0;
                        exprset[i]->src.inc_min = 0;
                        exprset[i]->src.exc_min = 0;

                        exprset[i]->src.max = 0;
                        exprset[i]->src.inc_max = 0;
                        exprset[i]->src.exc_max = 0;

                        exprset[i]->dst.min = 0;
                        exprset[i]->dst.max = 0;
                }
                else {
                        /* range */
                        exprset[i]->src.min = 0;
                        exprset[i]->src.inc_min = 0;
                        exprset[i]->src.exc_min = 0;

                        exprset[i]->src.max = 1000;
                        exprset[i]->src.inc_max = 1;
                        exprset[i]->src.exc_max = 0;

                        exprset[i]->dst.min = 1;
                        exprset[i]->dst.max = 255;
                }
        }

        newband = rt_band_reclass(band, PT_8BUI, 0, 0, exprset, cnt);
        CU_ASSERT(newband != NULL);

        rtn = rt_band_get_pixel(newband, 0, 0, &val, NULL);
        CU_ASSERT_EQUAL(rtn, ES_NONE);
        CU_ASSERT_DOUBLE_EQUAL(val, 0, DBL_EPSILON);

        rtn = rt_band_get_pixel(newband, 49, 5, &val, NULL);
        CU_ASSERT_EQUAL(rtn, ES_NONE);
        CU_ASSERT_DOUBLE_EQUAL(val, 77, DBL_EPSILON);

        rtn = rt_band_get_pixel(newband, 99, 9, &val, NULL);
        CU_ASSERT_EQUAL(rtn, ES_NONE);
        CU_ASSERT_DOUBLE_EQUAL(val, 255, DBL_EPSILON);

        for (i = cnt - 1; i >= 0; i--) rtdealloc(exprset[i]);
        rtdealloc(exprset);
        cu_free_raster(raster);

        rt_band_destroy(newband);
}

Here is the call graph for this function:

Here is the caller graph for this function: