rt_raster_from_band()

rt_raster rt_raster_from_band	(	rt_raster	raster,
		uint32_t *	bandNums,
		int	count
	)

Construct a new rt_raster from an existing rt_raster and an array of band numbers.

Parameters

raster	: the source raster
bandNums	: array of band numbers to extract from source raster and add to the new raster (0 based)
count	: number of elements in bandNums

Returns

a new rt_raster or NULL on error

Definition at line 1439 of file rt_raster.c.

References rt_raster_t::bands, genraster::count, window::gt, rt_raster_t::height, rtpixdump::rast, RASTER_DEBUGF, rt_band_destroy(), rt_raster_copy_band(), rt_raster_destroy(), rt_raster_get_geotransform_matrix(), rt_raster_get_num_bands(), rt_raster_new(), rt_raster_set_geotransform_matrix(), rt_raster_set_srid(), rterror(), rt_raster_t::srid, and rt_raster_t::width.

Referenced by RASTER_addBandRasterArray(), RASTER_band(), RASTER_union_finalfn(), rt_raster_clone(), and test_raster_from_band().

                                                                     {
        rt_raster rast = NULL;
        int i = 0;
        int j = 0;
        int idx;
        int32_t flag;
        double gt[6] = {0.};

        assert(NULL != raster);
        assert(NULL != bandNums);

        RASTER_DEBUGF(3, "rt_raster_from_band: source raster has %d bands",
                rt_raster_get_num_bands(raster));

        /* create new raster */
        rast = rt_raster_new(raster->width, raster->height);
        if (NULL == rast) {
                rterror("rt_raster_from_band: Out of memory allocating new raster");
                return NULL;
        }

        /* copy raster attributes */
        rt_raster_get_geotransform_matrix(raster, gt);
        rt_raster_set_geotransform_matrix(rast, gt);

        /* srid */
        rt_raster_set_srid(rast, raster->srid);

        /* copy bands */
        for (i = 0; i < count; i++) {
                idx = bandNums[i];
                flag = rt_raster_copy_band(rast, raster, idx, i);

                if (flag < 0) {
                        rterror("rt_raster_from_band: Could not copy band");
                        for (j = 0; j < i; j++) rt_band_destroy(rast->bands[j]);
                        rt_raster_destroy(rast);
                        return NULL;
                }

                RASTER_DEBUGF(3, "rt_raster_from_band: band created at index %d",
                        flag);
        }

        RASTER_DEBUGF(3, "rt_raster_from_band: new raster has %d bands",
                rt_raster_get_num_bands(rast));
        return rast;
}

Here is the call graph for this function:

Here is the caller graph for this function: