1620 FuncCallContext *funcctx;
1631 if (SRF_IS_FIRSTCALL()) {
1632 MemoryContext oldcontext;
1639 double *search = NULL;
1642 bool exclude_nodata_value =
TRUE;
1653 funcctx = SRF_FIRSTCALL_INIT();
1656 oldcontext = MemoryContextSwitchTo(funcctx->multi_call_memory_ctx);
1658 if (PG_ARGISNULL(0)) {
1659 MemoryContextSwitchTo(oldcontext);
1660 SRF_RETURN_DONE(funcctx);
1662 pgraster = (
rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
1665 PG_FREE_IF_COPY(pgraster, 0);
1666 MemoryContextSwitchTo(oldcontext);
1667 elog(ERROR,
"RASTER_pixelOfValue: Could not deserialize raster");
1668 SRF_RETURN_DONE(funcctx);
1673 if (num_bands < 1) {
1674 elog(NOTICE,
"Raster provided has no bands");
1676 PG_FREE_IF_COPY(pgraster, 0);
1677 MemoryContextSwitchTo(oldcontext);
1678 SRF_RETURN_DONE(funcctx);
1682 if (!PG_ARGISNULL(1))
1683 nband = PG_GETARG_INT32(1);
1684 if (nband < 1 || nband > num_bands) {
1685 elog(NOTICE,
"Invalid band index (must use 1-based). Returning NULL");
1687 PG_FREE_IF_COPY(pgraster, 0);
1688 MemoryContextSwitchTo(oldcontext);
1689 SRF_RETURN_DONE(funcctx);
1693 array = PG_GETARG_ARRAYTYPE_P(2);
1694 etype = ARR_ELEMTYPE(array);
1695 get_typlenbyvalalign(etype, &typlen, &typbyval, &typalign);
1703 PG_FREE_IF_COPY(pgraster, 0);
1704 MemoryContextSwitchTo(oldcontext);
1705 elog(ERROR,
"RASTER_pixelOfValue: Invalid data type for pixel values");
1706 SRF_RETURN_DONE(funcctx);
1710 deconstruct_array(array, etype, typlen, typbyval, typalign, &e,
1713 search = palloc(
sizeof(
double) * n);
1714 for (i = 0, nsearch = 0; i < n; i++) {
1715 if (nulls[i])
continue;
1719 val = (double) DatumGetFloat4(e[i]);
1722 val = (double) DatumGetFloat8(e[i]);
1726 search[nsearch] = val;
1733 elog(NOTICE,
"No search values provided. Returning NULL");
1736 PG_FREE_IF_COPY(pgraster, 0);
1737 MemoryContextSwitchTo(oldcontext);
1738 SRF_RETURN_DONE(funcctx);
1740 else if (nsearch < n)
1741 search = repalloc(search,
sizeof(
double) * nsearch);
1744 if (!PG_ARGISNULL(3))
1745 exclude_nodata_value = PG_GETARG_BOOL(3);
1750 elog(NOTICE,
"Could not find band at index %d. Returning NULL",
nband);
1752 PG_FREE_IF_COPY(pgraster, 0);
1753 MemoryContextSwitchTo(oldcontext);
1754 SRF_RETURN_DONE(funcctx);
1759 band, exclude_nodata_value,
1766 PG_FREE_IF_COPY(pgraster, 0);
1770 elog(NOTICE,
"Could not get the pixels of search values for band at index %d",
nband);
1773 elog(NOTICE,
"No pixels of search values found for band at index %d",
nband);
1775 MemoryContextSwitchTo(oldcontext);
1776 SRF_RETURN_DONE(funcctx);
1780 funcctx->user_fctx = pixels;
1783 funcctx->max_calls =
count;
1786 if (get_call_result_type(fcinfo, NULL, &tupdesc) != TYPEFUNC_COMPOSITE) {
1788 errcode(ERRCODE_FEATURE_NOT_SUPPORTED),
1790 "function returning record called in context "
1791 "that cannot accept type record"
1796 BlessTupleDesc(tupdesc);
1797 funcctx->tuple_desc = tupdesc;
1799 MemoryContextSwitchTo(oldcontext);
1803 funcctx = SRF_PERCALL_SETUP();
1805 call_cntr = funcctx->call_cntr;
1806 max_calls = funcctx->max_calls;
1807 tupdesc = funcctx->tuple_desc;
1808 pixels2 = funcctx->user_fctx;
1811 if (call_cntr < max_calls) {
1820 pixels2[call_cntr].
x += 1;
1821 pixels2[call_cntr].
y += 1;
1823 values[0] = Float8GetDatum(pixels2[call_cntr].
value);
1824 values[1] = Int32GetDatum(pixels2[call_cntr].
x);
1825 values[2] = Int32GetDatum(pixels2[call_cntr].
y);
1828 tuple = heap_form_tuple(tupdesc, values, nulls);
1831 result = HeapTupleGetDatum(tuple);
1833 SRF_RETURN_NEXT(funcctx,
result);
1837 SRF_RETURN_DONE(funcctx);
char result[OUT_DOUBLE_BUFFER_SIZE]
int rt_band_get_pixel_of_value(rt_band band, int exclude_nodata_value, double *searchset, int searchcount, rt_pixel *pixels)
Search band for pixel(s) with search values.
void rt_raster_destroy(rt_raster raster)
Release memory associated to a raster.
void rt_band_destroy(rt_band band)
Destroy a raster band.
uint16_t rt_raster_get_num_bands(rt_raster raster)
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): ...
#define POSTGIS_RT_DEBUGF(level, msg,...)