rt_raster_generate_new_band()

int rt_raster_generate_new_band	(	rt_raster	raster,
		rt_pixtype	pixtype,
		double	initialvalue,
		uint32_t	hasnodata,
		double	nodatavalue,
		int	index
	)

Generate a new inline band and add it to a raster.

Memory is allocated in this function for band data.

Parameters

raster	: the raster to add a band to
pixtype	: the pixel type for the new band
initialvalue	: initial value for pixels
hasnodata	: indicates if the band has a nodata value
nodatavalue	: nodata value for the new band
index	: position to add the new band in the raster

Returns

identifier (position) for the just-added raster, or -1 on error

Definition at line 489 of file rt_raster.c.

  {
    rt_band band = NULL;
    int width = 0;
    int height = 0;
    int numval = 0;
    int datasize = 0;
    int oldnumbands = 0;
    int numbands = 0;
    void * mem = NULL;
 
    assert(NULL != raster);
 
    /* Make sure index is in a valid range */
    oldnumbands = rt_raster_get_num_bands(raster);
    if (index < 0)
        index = 0;
    else if (index > oldnumbands + 1)
        index = oldnumbands + 1;
 
    /* Determine size of memory block to allocate and allocate it */
    width = rt_raster_get_width(raster);
    height = rt_raster_get_height(raster);
    numval = width * height;
    datasize = rt_pixtype_size(pixtype) * numval;
 
    mem = (int *)rtalloc(datasize);
    if (!mem) {
        rterror("rt_raster_generate_new_band: Could not allocate memory for band");
        return -1;
    }
 
    band = rt_band_new_inline(width, height, pixtype, hasnodata, nodatavalue, mem);
    if (! band) {
        rterror("rt_raster_generate_new_band: Could not add band to raster. Aborting");
        rtdealloc(mem);
        return -1;
    }
 
        rt_band_init_value(band, initialvalue);
 
        rt_band_set_ownsdata_flag(band, 1); /* we DO own this data!!! */
        index = rt_raster_add_band(raster, band, index);
        numbands = rt_raster_get_num_bands(raster);
        if (numbands == oldnumbands || index == -1) {
                rterror("rt_raster_generate_new_band: Could not add band to raster. Aborting");
                rt_band_destroy(band);
                return -1;
        }
 
        /* set isnodata if hasnodata = TRUE and initial value = nodatavalue */
        if (hasnodata && FLT_EQ(initialvalue, nodatavalue))
                rt_band_set_isnodata_flag(band, 1);
 
    return index;
}

References ovdump::band, FLT_EQ, rtrowdump::raster, rt_band_destroy(), rt_band_init_value(), rt_band_new_inline(), rt_band_set_isnodata_flag(), rt_band_set_ownsdata_flag(), rt_pixtype_size(), rt_raster_add_band(), rt_raster_get_height(), rt_raster_get_num_bands(), rt_raster_get_width(), rtalloc(), rtdealloc(), and rterror().

Referenced by RASTER_addBand(), RASTER_clip(), RASTER_mapAlgebra2(), RASTER_mapAlgebraExpr(), RASTER_mapAlgebraFct(), RASTER_mapAlgebraFctNgb(), RASTER_tile(), RASTER_union_transfn(), rt_raster_from_gdal_dataset(), rt_raster_iterator(), and test_raster_get_pixel_bilinear().

Here is the call graph for this function:

Here is the caller graph for this function: