4140 if (PG_ARGISNULL(0))
4146 elog(ERROR,
"RASTER_colorMap: Could not initialize argument structure");
4151 pgraster = (
rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
4157 PG_FREE_IF_COPY(pgraster, 0);
4158 elog(ERROR,
"RASTER_colorMap: Could not deserialize raster");
4163 if (!PG_ARGISNULL(1))
4164 arg->
nband = PG_GETARG_INT32(1);
4169 elog(NOTICE,
"Raster does not have band at index %d. Returning empty raster", arg->
nband);
4174 PG_FREE_IF_COPY(pgraster, 0);
4175 elog(ERROR,
"RASTER_colorMap: Could not create empty raster");
4180 PG_FREE_IF_COPY(pgraster, 0);
4184 if (pgraster == NULL)
4187 SET_VARSIZE(pgraster, ((
rt_pgraster*) pgraster)->size);
4188 PG_RETURN_POINTER(pgraster);
4193 if (arg->
band == NULL) {
4196 PG_FREE_IF_COPY(pgraster, 0);
4197 elog(ERROR,
"RASTER_colorMap: Could not get band at index %d",
nband);
4202 if (!PG_ARGISNULL(3)) {
4203 char *method = NULL;
4204 char *tmp = text_to_cstring(PG_GETARG_TEXT_P(3));
4211 if (strcmp(method,
"INTERPOLATE") == 0)
4213 else if (strcmp(method,
"EXACT") == 0)
4215 else if (strcmp(method,
"NEAREST") == 0)
4218 elog(NOTICE,
"Unknown value provided for method. Defaulting to INTERPOLATE");
4228 if (PG_ARGISNULL(2)) {
4230 PG_FREE_IF_COPY(pgraster, 0);
4231 elog(ERROR,
"RASTER_colorMap: Value must be provided for colormap");
4236 char *colormap = text_to_cstring(PG_GETARG_TEXT_P(2));
4245 if (!strlen(colormap)) {
4247 PG_FREE_IF_COPY(pgraster, 0);
4248 elog(ERROR,
"RASTER_colorMap: Value must be provided for colormap");
4256 PG_FREE_IF_COPY(pgraster, 0);
4257 elog(ERROR,
"RASTER_colorMap: Could not process the value provided for colormap");
4265 PG_FREE_IF_COPY(pgraster, 0);
4266 elog(ERROR,
"RASTER_colorMap: Could not allocate memory for colormap entries");
4272 for (i = 0; i < arg->
nentry; i++) {
4286 if (!strlen(_entry)) {
4296 PG_FREE_IF_COPY(pgraster, 0);
4297 elog(ERROR,
"RASTER_colorMap: Could not process colormap entry %d", i + 1);
4301 elog(NOTICE,
"More than five elements in colormap entry %d. Using at most five elements", i + 1);
4310 for (j = 0; j < arg->
nelement; j++) {
4319 char *percent = NULL;
4323 strcmp(_element,
"NV") == 0 ||
4324 strcmp(_element,
"NULL") == 0 ||
4325 strcmp(_element,
"NODATA") == 0
4330 elog(NOTICE,
"More than one NODATA entry found. Using only the first one");
4339 else if ((percent = strchr(_element,
'%')) != NULL) {
4351 PG_FREE_IF_COPY(pgraster, 0);
4352 elog(ERROR,
"RASTER_colorMap: Could not get band's summary stats to process percentages");
4358 tmp = palloc(
sizeof(
char) * (percent - _element + 1));
4362 PG_FREE_IF_COPY(pgraster, 0);
4363 elog(ERROR,
"RASTER_colorMap: Could not allocate memory for value of percentage");
4367 memcpy(tmp, _element, percent - _element);
4368 tmp[percent - _element] =
'\0';
4373 value = strtod(tmp, NULL);
4375 if (errno != 0 || _element == junk) {
4378 PG_FREE_IF_COPY(pgraster, 0);
4379 elog(ERROR,
"RASTER_colorMap: Could not process percent string to value");
4385 elog(NOTICE,
"Percentage values cannot be less than zero. Defaulting to zero");
4388 else if (
value > 100.) {
4389 elog(NOTICE,
"Percentage values cannot be greater than 100. Defaulting to 100");
4401 if (errno != 0 || _element == junk) {
4404 PG_FREE_IF_COPY(pgraster, 0);
4405 elog(ERROR,
"RASTER_colorMap: Could not process string to value");
4416 value = (int) strtod(_element, &junk);
4417 if (errno != 0 || _element == junk) {
4420 PG_FREE_IF_COPY(pgraster, 0);
4421 elog(ERROR,
"RASTER_colorMap: Could not process string to value");
4426 elog(NOTICE,
"RGBA value cannot be greater than 255. Defaulting to 255");
4429 else if (
value < 0) {
4430 elog(NOTICE,
"RGBA value cannot be less than zero. Defaulting to zero");
4439 POSTGIS_RT_DEBUGF(4,
"colormap->entry[%d] (isnodata, value, R, G, B, A) = (%d, %f, %d, %d, %d, %d)",
4460 PG_FREE_IF_COPY(pgraster, 0);
4461 elog(ERROR,
"RASTER_colorMap: Could not create new raster with applied colormap");
4466 PG_FREE_IF_COPY(pgraster, 0);
4472 if (pgraster == NULL)
4475 SET_VARSIZE(pgraster, ((
rt_pgraster*) pgraster)->size);
4476 PG_RETURN_POINTER(pgraster);
rt_raster rt_raster_colormap(rt_raster raster, int nband, rt_colormap colormap)
Returns a new raster with up to four 8BUI bands (RGBA) from applying a colormap to the user-specified...
void rt_raster_destroy(rt_raster raster)
Release memory associated to a raster.
int rt_raster_has_band(rt_raster raster, int nband)
Return TRUE if the raster has a band of this number.
void * rt_raster_serialize(rt_raster raster)
Return this raster in serialized form.
rt_raster rt_raster_clone(rt_raster raster, uint8_t deep)
Clone an existing raster.
rt_bandstats rt_band_get_summary_stats(rt_band band, int exclude_nodata_value, double sample, int inc_vals, uint64_t *cK, double *cM, double *cQ)
Compute summary statistics for a band.
rt_raster rt_raster_deserialize(void *serialized, int header_only)
Return a raster from a serialized form.
rt_band rt_raster_get_band(rt_raster raster, int bandNum)
Return Nth band, or NULL if unavailable.
raster
Be careful!! Zeros function's input parameter can be a (height x width) array, not (width x height): ...
char * rtpg_strtoupper(char *str)
char * rtpg_trim(const char *input)
char ** rtpg_strsplit(const char *str, const char *delimiter, uint32_t *n)
char * rtpg_strreplace(const char *str, const char *oldstr, const char *newstr, int *count)
static void rtpg_colormap_arg_destroy(rtpg_colormap_arg arg)
static rtpg_colormap_arg rtpg_colormap_arg_init()
#define POSTGIS_RT_DEBUG(level, msg)
#define POSTGIS_RT_DEBUGF(level, msg,...)
enum rt_colormap_t::@10 method