rtpg_union_callback()

static int rtpg_union_callback ( rt_iterator_arg  arg, void *  userarg, double *  value, int *  nodata  )

static

Definition at line 1787 of file rtpg_mapalgebra.c.

References rt_iterator_arg_t::columns, rt_iterator_arg_t::nodata, POSTGIS_RT_DEBUGF, rt_iterator_arg_t::rasters, rt_iterator_arg_t::rows, UT_COUNT, UT_FIRST, UT_LAST, UT_MAX, UT_MEAN, UT_MIN, UT_RANGE, UT_SUM, and rt_iterator_arg_t::values.

Referenced by RASTER_union_transfn().

  {
        rtpg_union_type *utype = (rtpg_union_type *) userarg;

        if (arg == NULL)
                return 0;

        if (
                arg->rasters != 2 ||
                arg->rows != 1 ||
                arg->columns != 1
        ) {
                elog(ERROR, "rtpg_union_callback: Invalid arguments passed to callback");
                return 0;
        }

        *value = 0;
        *nodata = 0;

        /* handle NODATA situations except for COUNT, which is a special case */
        if (*utype != UT_COUNT) {
                /* both NODATA */
                if (arg->nodata[0][0][0] && arg->nodata[1][0][0]) {
                        *nodata = 1;
                        POSTGIS_RT_DEBUGF(4, "value, nodata = %f, %d", *value, *nodata);
                        return 1;
                }
                /* second NODATA */
                else if (!arg->nodata[0][0][0] && arg->nodata[1][0][0]) {
                        *value = arg->values[0][0][0];
                        POSTGIS_RT_DEBUGF(4, "value, nodata = %f, %d", *value, *nodata);
                        return 1;
                }
                /* first NODATA */
                else if (arg->nodata[0][0][0] && !arg->nodata[1][0][0]) {
                        *value = arg->values[1][0][0];
                        POSTGIS_RT_DEBUGF(4, "value, nodata = %f, %d", *value, *nodata);
                        return 1;
                }
        }

        switch (*utype) {
                case UT_FIRST:
                        *value = arg->values[0][0][0];
                        break;
                case UT_MIN:
                        if (arg->values[0][0][0] < arg->values[1][0][0])
                                *value = arg->values[0][0][0];
                        else
                                *value = arg->values[1][0][0];
                        break;
                case UT_MAX:
                        if (arg->values[0][0][0] > arg->values[1][0][0])
                                *value = arg->values[0][0][0];
                        else
                                *value = arg->values[1][0][0];
                        break;
                case UT_COUNT:
                        /* both NODATA */
                        if (arg->nodata[0][0][0] && arg->nodata[1][0][0])
                                *value = 0;
                        /* second NODATA */
                        else if (!arg->nodata[0][0][0] && arg->nodata[1][0][0])
                                *value = arg->values[0][0][0];
                        /* first NODATA */
                        else if (arg->nodata[0][0][0] && !arg->nodata[1][0][0])
                                *value = 1;
                        /* has value, increment */
                        else
                                *value = arg->values[0][0][0] + 1;
                        break;
                case UT_SUM:
                        *value = arg->values[0][0][0] + arg->values[1][0][0];
                        break;
                case UT_MEAN:
                case UT_RANGE:
                        break;
                case UT_LAST:
                default:
                        *value = arg->values[1][0][0];
                        break;
        }

        POSTGIS_RT_DEBUGF(4, "value, nodata = %f, %d", *value, *nodata);


        return 1;
}

Here is the caller graph for this function: