test_pixel_set_to_array()

static void test_pixel_set_to_array ( )

static

Definition at line 370 of file cu_band_misc.c.

                                     {
        rt_raster rast;
        rt_band band;
        rt_mask mask = NULL;
        int x, y;
        int rtn;
        const int maxX = 10;
        const int maxY = 10;
        int maskX = 3;
        int maskY = 3;
        rt_pixel npixels = NULL;
        int i;
        double **value;
        double val;
        int nod;
        int **nodata;
        int dimx;
        int dimy;
 
        rast = rt_raster_new(maxX, maxY);
        CU_ASSERT(rast != NULL);
 
        band = cu_add_band(rast, PT_32BF, 1, 0);
        CU_ASSERT(band != NULL);
 
        for (x = 0; x < maxX; x++) {
                for (y = 0; y < maxY; y++) {
                        rtn = rt_band_set_pixel(band, x, y, 1, NULL);
                }
        }
 
 
        rtn = rt_band_get_pixel(band,4,4,&val,&nod);
        CU_ASSERT_EQUAL(nod,0);
        CU_ASSERT_DOUBLE_EQUAL(val,1,.01);
 
        /* set up mask */
 
        mask = (rt_mask) rtalloc(sizeof(struct rt_mask_t) );
        CU_ASSERT(mask != NULL);
        mask->values = rtalloc(sizeof(double*)*maskY);
        mask->nodata = rtalloc(sizeof(int*)*maskY);
 
        for( i = 0; i < maskY;  i++) {
          mask->values[i] = rtalloc(sizeof(double) *maskX);
          mask->nodata[i] = rtalloc(sizeof(int) *maskX);
        }
 
        CU_ASSERT(mask->values != NULL);
        CU_ASSERT(mask->nodata != NULL);
 
        /* set mask to nodata */
 
        for(y = 0; y < maskY; y++) {
          for(x = 0; x < maskX; x++){
            mask->values[y][x]= 0;
            mask->nodata[y][x]= 1;
          }
        }
 
        mask->dimx = maskX;
        mask->dimy = maskY;
        mask->weighted = 0;
 
       rtn = rt_band_get_nearest_pixel(
                band,
                4,4,
                1, 1,
                1,
                &npixels
        );
 
       CU_ASSERT_EQUAL(rtn,8);
 
        rtn = rt_pixel_set_to_array(
                npixels, rtn, mask,
                4,4,
                1, 1,
                &value,
                &nodata,
                &dimx, &dimy
        );
 
 
       rtdealloc(npixels);
       CU_ASSERT_EQUAL(rtn, ES_NONE);
       CU_ASSERT_EQUAL(dimx, 3);
       CU_ASSERT_EQUAL(dimy, 3);
       CU_ASSERT_EQUAL(nodata[1][1],1);
 
        for (x = 0; x < dimx; x++) {
                rtdealloc(nodata[x]);
                rtdealloc(value[x]);
        }
 
        rtdealloc(nodata);
        rtdealloc(value);
 
        /* set mask to 1 */
 
        for(y = 0; y < maskY; y++) {
          for(x = 0; x < maskX; x++){
            mask->values[y][x]= 1;
            mask->nodata[y][x]= 0;
          }
        }
 
        mask->dimx = maskX;
        mask->dimy = maskY;
        mask->weighted = 0;
 
 
       rtn = rt_band_get_nearest_pixel(
                band,
                4,4,
                1, 1,
                1,
                &npixels
        );
 
       CU_ASSERT_EQUAL(rtn,8);
 
        rtn = rt_pixel_set_to_array(
                npixels, rtn, mask,
                4,4,
                1, 1,
                &value,
                &nodata,
                &dimx, &dimy
        );
 
 
       rtdealloc(npixels);
       CU_ASSERT_EQUAL(rtn, ES_NONE);
       CU_ASSERT_EQUAL(dimx, 3);
       CU_ASSERT_EQUAL(dimy, 3);
       CU_ASSERT_NOT_EQUAL(nodata[0][0],1);
       CU_ASSERT_DOUBLE_EQUAL(value[0][0],1,.01);
 
        for (x = 0; x < dimx; x++) {
                rtdealloc(nodata[x]);
                rtdealloc(value[x]);
        }
 
        rtdealloc(nodata);
        rtdealloc(value);
 
        /* set mask to 0.5 */
 
        for(y = 0; y < maskY; y++) {
          for(x = 0; x < maskX; x++){
            mask->values[y][x]= 0.5;
            mask->nodata[y][x]= 0;
          }
        }
 
        mask->dimx = maskX;
        mask->dimy = maskY;
        mask->weighted = 1;
 
       rtn = rt_band_get_nearest_pixel(
                band,
                4,4,
                1, 1,
                1,
                &npixels
        );
 
       CU_ASSERT_EQUAL(rtn,8);
 
        rtn = rt_pixel_set_to_array(
                npixels, rtn, mask,
                4,4,
                1, 1,
                &value,
                &nodata,
                &dimx, &dimy
        );
 
 
       rtdealloc(npixels);
       CU_ASSERT_EQUAL(rtn, ES_NONE);
       CU_ASSERT_EQUAL(dimx, 3);
       CU_ASSERT_EQUAL(dimy, 3);
       CU_ASSERT_NOT_EQUAL(nodata[0][0],1);
       CU_ASSERT_DOUBLE_EQUAL(value[0][0],0.5,0.1);
 
        for (x = 0; x < dimx; x++) {
                rtdealloc(nodata[x]);
                rtdealloc(value[x]);
        }
 
        rtdealloc(nodata);
        rtdealloc(value);
 
        for( i = 0; i < maskY;  i++) {
          rtdealloc(mask->values[i]);
          rtdealloc(mask->nodata[i]);
        }
 
        rtdealloc(mask->values);
        rtdealloc(mask->nodata);
        rtdealloc(mask);
 
        if (rtn)
          rtdealloc(npixels);
 
       cu_free_raster(rast);
 
}

References cu_add_band(), cu_free_raster(), rt_mask_t::dimx, rt_mask_t::dimy, ES_NONE, rt_mask_t::nodata, PT_32BF, rt_band_get_nearest_pixel(), rt_band_get_pixel(), rt_band_set_pixel(), rt_pixel_set_to_array(), rt_raster_new(), rtalloc(), rtdealloc(), rt_mask_t::values, and rt_mask_t::weighted.

Referenced by band_misc_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: