rt_band_corrected_clamped_value()

rt_errorstate rt_band_corrected_clamped_value	(	rt_band	band,
		double	val,
		double *	newval,
		int *	corrected
	)

Correct value when clamped value is equal to clamped NODATA value.

Correction does NOT occur if unclamped value is exactly unclamped NODATA value.

Parameters

band	: the band whose NODATA value will be used for comparison
val	: the value to compare to the NODATA value and correct
*newval	: pointer to corrected value
*corrected	: (optional) non-zero if val was corrected

Returns

ES_NONE if success, ES_ERROR if error

Definition at line 1701 of file rt_band.c.

References ES_ERROR, ES_NONE, FLT_EQ, rt_band_t::nodataval, rt_band_t::pixtype, PT_16BSI, PT_16BUI, PT_1BB, PT_2BUI, PT_32BF, PT_32BSI, PT_32BUI, PT_4BUI, PT_64BF, PT_8BSI, PT_8BUI, rt_band_clamped_value_is_nodata(), rt_pixtype_get_min_value(), rt_util_clamp_to_16BSI(), rt_util_clamp_to_16BUI(), rt_util_clamp_to_2BUI(), rt_util_clamp_to_32BSI(), rt_util_clamp_to_32BUI(), rt_util_clamp_to_32F(), rt_util_clamp_to_4BUI(), rt_util_clamp_to_8BSI(), rt_util_clamp_to_8BUI(), and rterror().

Referenced by rt_band_set_pixel().

  {
        double minval = 0.;

        assert(NULL != band);
        assert(NULL != newval);

        if (corrected != NULL)
                *corrected = 0;

        /* no need to correct if clamped values IS NOT clamped NODATA */
        if (rt_band_clamped_value_is_nodata(band, val) != 1) {
                *newval = val;
                return ES_NONE;
        }

        minval = rt_pixtype_get_min_value(band->pixtype);
        *newval = val;

        switch (band->pixtype) {
                case PT_1BB:
                        *newval = !band->nodataval;
                        break;
                case PT_2BUI:
                        if (rt_util_clamp_to_2BUI(val) == rt_util_clamp_to_2BUI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_4BUI:
                        if (rt_util_clamp_to_4BUI(val) == rt_util_clamp_to_4BUI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_8BSI:
                        if (rt_util_clamp_to_8BSI(val) == rt_util_clamp_to_8BSI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_8BUI:
                        if (rt_util_clamp_to_8BUI(val) == rt_util_clamp_to_8BUI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_16BSI:
                        if (rt_util_clamp_to_16BSI(val) == rt_util_clamp_to_16BSI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_16BUI:
                        if (rt_util_clamp_to_16BUI(val) == rt_util_clamp_to_16BUI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_32BSI:
                        if (rt_util_clamp_to_32BSI(val) == rt_util_clamp_to_32BSI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_32BUI:
                        if (rt_util_clamp_to_32BUI(val) == rt_util_clamp_to_32BUI(minval))
                                (*newval)++;
                        else
                                (*newval)--;
                        break;
                case PT_32BF:
                        if (FLT_EQ(rt_util_clamp_to_32F(val), rt_util_clamp_to_32F(minval)))
                                *newval += FLT_EPSILON;
                        else
                                *newval -= FLT_EPSILON;
                        break;
                case PT_64BF:
                        break;
                default:
                        rterror("rt_band_corrected_clamped_value: Unknown pixeltype %d", band->pixtype);
                        return ES_ERROR;
        }

        if (corrected != NULL)
                *corrected = 1;

        return ES_NONE;
}

Here is the call graph for this function:

Here is the caller graph for this function: