rt_band_set_pixel()

rt_errorstate rt_band_set_pixel	(	rt_band	band,
		int	x,
		int	y,
		double	val,
		int *	converted
	)

Set single pixel's value.

Parameters

band	: the band to set value to
x	: x ordinate (0-based)
y	: y ordinate (0-based)
val	: the pixel value
converted	: (optional) non-zero if value truncated/clamped/converted

Returns

ES_NONE on success, ES_ERROR on error

Definition at line 1140 of file rt_band.c.

  {
        rt_pixtype pixtype = PT_END;
        uint8_t *data = NULL;
        uint32_t offset = 0;
 
        int32_t checkvalint = 0;
        uint32_t checkvaluint = 0;
        float checkvalfloat = 0;
        double checkvaldouble = 0;
 
        assert(NULL != band);
 
        if (converted != NULL)
                *converted = 0;
 
        if (band->offline) {
                rterror("rt_band_set_pixel not implemented yet for OFFDB bands");
                return ES_ERROR;
        }
 
        pixtype = band->pixtype;
 
        if (
                x < 0 || x >= band->width ||
                y < 0 || y >= band->height
        ) {
                rterror("rt_band_set_pixel: Coordinates out of range");
                return ES_ERROR;
        }
 
        /* check that clamped value isn't clamped NODATA */
        if (band->hasnodata && pixtype != PT_64BF) {
                double newval;
                int corrected;
 
                rt_band_corrected_clamped_value(band, val, &newval, &corrected);
 
                if (corrected) {
#if POSTGIS_RASTER_WARN_ON_TRUNCATION > 0
                        rtwarn("Value for pixel %d x %d has been corrected as clamped value becomes NODATA", x, y);
#endif
                        val = newval;
 
                        if (converted != NULL)
                                *converted = 1;
                }
        }
 
        data = rt_band_get_data(band);
        offset = x + (y * band->width);
 
        switch (pixtype)
        {
        case PT_1BB: {
                data[offset] = rt_util_clamp_to_1BB(val);
                checkvalint = data[offset];
                break;
        }
                case PT_2BUI: {
                        data[offset] = rt_util_clamp_to_2BUI(val);
                        checkvalint = data[offset];
                        break;
                }
                case PT_4BUI: {
                        data[offset] = rt_util_clamp_to_4BUI(val);
                        checkvalint = data[offset];
                        break;
                }
                case PT_8BSI: {
                        data[offset] = (uint8_t)rt_util_clamp_to_8BSI(val);
                        checkvalint = (int8_t) data[offset];
                        break;
                }
                case PT_8BUI: {
                        data[offset] = rt_util_clamp_to_8BUI(val);
                        checkvalint = data[offset];
                        break;
                }
                case PT_16BSI: {
                        int16_t *ptr = (int16_t*) data; /* we assume correct alignment */
                        ptr[offset] = rt_util_clamp_to_16BSI(val);
                        checkvalint = (int16_t) ptr[offset];
                        break;
                }
                case PT_16BUI: {
                        uint16_t *ptr = (uint16_t*) data; /* we assume correct alignment */
                        ptr[offset] = rt_util_clamp_to_16BUI(val);
                        checkvalint = ptr[offset];
                        break;
                }
                case PT_32BSI: {
                        int32_t *ptr = (int32_t*) data; /* we assume correct alignment */
                        ptr[offset] = rt_util_clamp_to_32BSI(val);
                        checkvalint = (int32_t) ptr[offset];
                        break;
                }
                case PT_32BUI: {
                        uint32_t *ptr = (uint32_t *)data; /* we assume correct alignment */
                        ptr[offset] = rt_util_clamp_to_32BUI(val);
                        checkvaluint = ptr[offset];
                        break;
                }
                case PT_16BF: {
                        uint16_t *ptr = (uint16_t *)data; /* we assume correct alignment */
                        float clamped = rt_util_clamp_to_16F(val);
                        /* Pack via explicit converter to mirror GDAL's Float16 representation. */
                        ptr[offset] = rt_util_float_to_float16(clamped);
                        checkvalfloat = rt_util_float16_to_float(ptr[offset]);
                        break;
                }
                case PT_32BF: {
                        float *ptr = (float *)data; /* we assume correct alignment */
                        ptr[offset] = rt_util_clamp_to_32F(val);
                        checkvalfloat = ptr[offset];
                        break;
                }
                case PT_64BF: {
                        double *ptr = (double*) data; /* we assume correct alignment */
                        ptr[offset] = val;
                        checkvaldouble = ptr[offset];
                        break;
                }
                default: {
                        rterror("rt_band_set_pixel: Unknown pixeltype %d", pixtype);
                        return ES_ERROR;
                }
                }
 
        /* If the stored value is not NODATA, reset the isnodata flag */
        if (!rt_band_clamped_value_is_nodata(band, val)) {
                RASTER_DEBUG(3, "Band has a value that is not NODATA. Setting isnodata to FALSE");
                band->isnodata = FALSE;
        }
 
        /* Overflow checking */
        if (rt_util_dbl_trunc_warning(
                val,
                checkvalint, checkvaluint,
                checkvalfloat, checkvaldouble,
                pixtype
        ) && converted != NULL) {
                *converted = 1;
        }
 
        return ES_NONE;
}

References ES_ERROR, ES_NONE, FALSE, PT_16BF, PT_16BSI, PT_16BUI, PT_1BB, PT_2BUI, PT_32BF, PT_32BSI, PT_32BUI, PT_4BUI, PT_64BF, PT_8BSI, PT_8BUI, PT_END, RASTER_DEBUG, rt_band_clamped_value_is_nodata(), rt_band_corrected_clamped_value(), rt_band_get_data(), rt_util_clamp_to_16BSI(), rt_util_clamp_to_16BUI(), rt_util_clamp_to_16F(), rt_util_clamp_to_1BB(), 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(), rt_util_dbl_trunc_warning(), rt_util_float16_to_float(), rt_util_float_to_float16(), rterror(), and rtwarn().

Referenced by fillRasterToPolygonize(), RASTER_clip(), RASTER_mapAlgebra2(), RASTER_mapAlgebraExpr(), RASTER_mapAlgebraFct(), RASTER_mapAlgebraFctNgb(), RASTER_setPixelValue(), RASTER_setPixelValuesArray(), RASTER_setPixelValuesGeomval(), rt_band_reclass(), rt_band_reclass_exact(), rt_raster_gdal_rasterize(), rt_raster_iterator(), test_band_get_nearest_pixel(), test_band_get_pixel_line(), test_band_get_pixel_of_value(), test_band_pixtype_16BSI(), test_band_pixtype_16BUI(), test_band_pixtype_1BB(), test_band_pixtype_2BUI(), test_band_pixtype_32BF(), test_band_pixtype_32BSI(), test_band_pixtype_32BUI(), test_band_pixtype_4BUI(), test_band_pixtype_64BF(), test_band_pixtype_8BSI(), test_band_pixtype_8BUI(), test_band_reclass(), test_band_stats(), test_band_value_count(), test_gdal_to_raster(), test_gdal_warp(), test_pixel_set_to_array(), test_raster_colormap(), test_raster_fully_within_distance(), test_raster_geos_contains(), test_raster_geos_contains_properly(), test_raster_geos_covered_by(), test_raster_geos_covers(), test_raster_geos_overlaps(), test_raster_geos_touches(), test_raster_get_pixel_bilinear(), test_raster_intersects(), test_raster_iterator(), test_raster_perimeter(), test_raster_pixel_as_polygon(), test_raster_surface(), test_raster_to_gdal(), and test_raster_within_distance().

Here is the call graph for this function:

Here is the caller graph for this function: