RASTER_getGeometryValues()

Datum RASTER_getGeometryValues

(

PG_FUNCTION_ARGS

)

Definition at line 258 of file rtpg_pixel.c.

{
        rt_raster raster = NULL;
        rt_pgraster *pgraster = NULL;
        GSERIALIZED *gser;
        LWGEOM *lwgeom_in, *lwgeom_out;
        rt_resample_type resample_type = RT_NEAREST;
        rt_errorstate err;
        char dimension;
        const char *func_name;
        uint16_t num_bands;
        int32_t band;
 
        text *resample_text = PG_GETARG_TEXT_P(2);
 
        /* Dimension depends on the name of calling SQL function */
        /* ST_SetZ()? or ST_SetM()? */
        func_name = get_func_name(fcinfo->flinfo->fn_oid);
        if (strcmp(func_name, "st_setz") == 0)
                dimension = 'z';
        else if (strcmp(func_name, "st_setm") == 0)
                dimension = 'm';
        else
                elog(ERROR, "%s called from unexpected SQL signature", __func__);
 
        /* Geometry */
        gser = (GSERIALIZED*)PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
        if (gserialized_is_empty(gser)) {
                elog(ERROR, "Cannot copy value into an empty geometry");
                PG_RETURN_NULL();
        }
 
        /* Raster */
        pgraster = (rt_pgraster *)PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
        raster = rt_raster_deserialize(pgraster, FALSE);
        num_bands = rt_raster_get_num_bands(raster);
        if (!raster) {
                elog(ERROR, "Could not deserialize raster");
                PG_RETURN_NULL();
        }
 
        /* Bandnidex is 1-based */
        band = PG_GETARG_INT32(3);
        if (band < 1 || band > num_bands) {
                elog(NOTICE, "Invalid band index %d. Must be between 1 and %u", band, num_bands);
                PG_RETURN_NULL();
        }
 
        /* SRID consistency */
        if (gserialized_get_srid(gser) != rt_raster_get_srid(raster)) {
                elog(ERROR, "Raster and geometry do not have the same SRID");
                PG_RETURN_NULL();
        }
 
        /* Process arguments */
        resample_type = resample_text_to_type(resample_text);
 
        /* Get the geometry */
        lwgeom_in = lwgeom_from_gserialized(gser);
 
        /* Run the sample */
        err = rt_raster_copy_to_geometry(
                raster,
                band - 1,      /* rtcore uses 0-based band number */
                dimension,     /* 'z' or 'm' */
                resample_type, /* bilinear or nearest */
                lwgeom_in,
                &lwgeom_out
        );
 
        rt_raster_destroy(raster);
        lwgeom_free(lwgeom_in);
        if (err != ES_NONE || !lwgeom_out) {
                PG_RETURN_NULL();
        }
        PG_RETURN_POINTER(gserialized_from_lwgeom(lwgeom_out, NULL));
}

References ES_NONE, FALSE, gserialized_from_lwgeom(), gserialized_get_srid(), gserialized_is_empty(), lwgeom_free(), lwgeom_from_gserialized(), resample_text_to_type(), RT_NEAREST, rt_raster_copy_to_geometry(), rt_raster_deserialize(), rt_raster_destroy(), rt_raster_get_num_bands(), and rt_raster_get_srid().

Here is the call graph for this function: