_rti_iterator_arg_populate()

static int _rti_iterator_arg_populate ( _rti_iterator_arg  _param, rt_iterator  itrset, uint16_t  itrcount, uint16_t  distancex, uint16_t  distancey, int *  allnull, int *  allempty  )

static

Definition at line 531 of file rt_mapalgebra.c.

References _rti_iterator_arg_t::band, _rti_iterator_arg_t::columns, _rti_iterator_arg_t::count, _rti_iterator_arg_t::dimension, _rti_iterator_arg_t::distance, _rti_iterator_arg_t::hasnodata, _rti_iterator_arg_t::height, _rti_iterator_arg_t::isempty, _rti_iterator_arg_t::isnodata, _rti_iterator_arg_t::minval, pixval::nband, _rti_iterator_arg_t::nodataval, _rti_iterator_arg_t::offset, _rti_iterator_arg_t::raster, rt_iterator_t::raster, RASTER_DEBUGF, _rti_iterator_arg_t::rows, rt_band_get_hasnodata_flag(), rt_band_get_isnodata_flag(), rt_band_get_min_value(), rt_band_get_nodata(), rt_raster_get_band(), rt_raster_get_height(), rt_raster_get_width(), rt_raster_has_band(), rt_raster_is_empty(), rtalloc(), _rti_iterator_arg_t::rtband, rterror(), _rti_iterator_arg_t::width, _rti_iterator_arg_t::x, and _rti_iterator_arg_t::y.

Referenced by rt_raster_iterator().

  {
        int i = 0;
        int hasband = 0;

        _param->count = itrcount;
        _param->distance.x = distancex;
        _param->distance.y = distancey;
        _param->dimension.columns = distancex * 2 + 1;
        _param->dimension.rows = distancey * 2 + 1;

        /* allocate memory for children */
        _param->raster = rtalloc(sizeof(rt_raster) * itrcount);
        _param->isempty = rtalloc(sizeof(int) * itrcount);
        _param->width = rtalloc(sizeof(int) * itrcount);
        _param->height = rtalloc(sizeof(int) * itrcount);

        _param->offset = rtalloc(sizeof(double *) * itrcount);

        _param->band.rtband = rtalloc(sizeof(rt_band) * itrcount);
        _param->band.hasnodata = rtalloc(sizeof(int) * itrcount);
        _param->band.isnodata = rtalloc(sizeof(int) * itrcount);
        _param->band.nodataval = rtalloc(sizeof(double) * itrcount);
        _param->band.minval = rtalloc(sizeof(double) * itrcount);

        if (
                _param->raster == NULL ||
                _param->isempty == NULL ||
                _param->width == NULL ||
                _param->height == NULL ||
                _param->offset == NULL ||
                _param->band.rtband == NULL ||
                _param->band.hasnodata == NULL ||
                _param->band.isnodata == NULL ||
                _param->band.nodataval == NULL ||
                _param->band.minval == NULL
        ) {
                rterror("_rti_iterator_arg_populate: Could not allocate memory for children of _rti_iterator_arg");
                return 0;
        }

        *allnull = 0;
        *allempty = 0;

        /*
                check input rasters
                        not empty, band # is valid
                        copy raster pointers and set flags
        */
        for (i = 0; i < itrcount; i++) {
                /* initialize elements */
                _param->raster[i] = NULL;
                _param->isempty[i] = 0;
                _param->width[i] = 0;
                _param->height[i] = 0;

                _param->offset[i] = NULL;

                _param->band.rtband[i] = NULL;
                _param->band.hasnodata[i] = 0;
                _param->band.isnodata[i] = 0;
                _param->band.nodataval[i] = 0;
                _param->band.minval[i] = 0;

                /* set isempty */
                if (itrset[i].raster == NULL) {
                        _param->isempty[i] = 1;

                        (*allnull)++;
                        (*allempty)++;

                        continue;
                }
                else if (rt_raster_is_empty(itrset[i].raster)) {
                        _param->isempty[i] = 1;

                        (*allempty)++;

                        continue;
                }

                /* check band number */
                hasband = rt_raster_has_band(itrset[i].raster, itrset[i].nband);
                if (!hasband) {
                        if (!itrset[i].nbnodata) {
                                rterror("_rti_iterator_arg_populate: Band %d not found for raster %d", itrset[i].nband, i);
                                return 0;
                        }
                        else {
                                RASTER_DEBUGF(4, "Band %d not found for raster %d. Using NODATA", itrset[i].nband, i);
                        }
                }

                _param->raster[i] = itrset[i].raster;
                if (hasband) {
                        _param->band.rtband[i] = rt_raster_get_band(itrset[i].raster, itrset[i].nband);
                        if (_param->band.rtband[i] == NULL) {
                                rterror("_rti_iterator_arg_populate: Could not get band %d for raster %d", itrset[i].nband, i);
                                return 0;
                        }

                        /* hasnodata */
                        _param->band.hasnodata[i] = rt_band_get_hasnodata_flag(_param->band.rtband[i]);

                        /* hasnodata = TRUE */
                        if (_param->band.hasnodata[i]) {
                                /* nodataval */
                                rt_band_get_nodata(_param->band.rtband[i], &(_param->band.nodataval[i]));

                                /* isnodata */
                                _param->band.isnodata[i] = rt_band_get_isnodata_flag(_param->band.rtband[i]);
                        }
                        /* hasnodata = FALSE */
                        else {
                                /* minval */
                                _param->band.minval[i] = rt_band_get_min_value(_param->band.rtband[i]);
                        }
                }

                /* width, height */
                _param->width[i] = rt_raster_get_width(_param->raster[i]);
                _param->height[i] = rt_raster_get_height(_param->raster[i]);

                /* init offset */
                _param->offset[i] = rtalloc(sizeof(double) * 2);
                if (_param->offset[i] == NULL) {
                        rterror("_rti_iterator_arg_populate: Could not allocate memory for offsets");
                        return 0;
                }
        }

        return 1;
}

Here is the call graph for this function:

Here is the caller graph for this function: