RASTER_rasterToWorldCoord()

Datum RASTER_rasterToWorldCoord

(

PG_FUNCTION_ARGS

)

Definition at line 693 of file rtpg_raster_properties.c.

{
        rt_pgraster *pgraster = NULL;
        rt_raster raster = NULL;
        int i = 0;
        int cr[2] = {0};
        bool skewed[2] = {false, false};
        double cw[2] = {0};
 
        TupleDesc tupdesc;
        Datum values[VALUES_LENGTH];
        bool nulls[VALUES_LENGTH];
        HeapTuple tuple;
        Datum result;
 
        POSTGIS_RT_DEBUG(3, "RASTER_rasterToWorldCoord: Starting");
 
        /* pgraster is null, return null */
        if (PG_ARGISNULL(0))
                PG_RETURN_NULL();
        pgraster = (rt_pgraster *) PG_DETOAST_DATUM_SLICE(PG_GETARG_DATUM(0), 0, sizeof(struct rt_raster_serialized_t));
 
        /* raster */
        raster = rt_raster_deserialize(pgraster, TRUE);
        if (!raster) {
                PG_FREE_IF_COPY(pgraster, 0);
                elog(ERROR, "RASTER_rasterToWorldCoord: Could not deserialize raster");
                PG_RETURN_NULL();
        }
 
        /* raster skewed? */
        skewed[0] = FLT_NEQ(rt_raster_get_x_skew(raster), 0.0) ? true : false;
        skewed[1] = FLT_NEQ(rt_raster_get_y_skew(raster), 0.0) ? true : false;
 
        /* column and row */
        for (i = 1; i <= 2; i++) {
                if (PG_ARGISNULL(i)) {
                        /* if skewed on same axis, parameter is required */
                        if (skewed[i - 1]) {
                                elog(NOTICE, "Pixel row and column required for computing longitude and latitude of a rotated raster");
                                rt_raster_destroy(raster);
                                PG_FREE_IF_COPY(pgraster, 0);
                                PG_RETURN_NULL();
                        }
 
                        continue;
                }
 
                cr[i - 1] = PG_GETARG_INT32(i);
        }
 
        /* user-provided value is 1-based but needs to be 0-based */
        if (rt_raster_cell_to_geopoint(
                raster,
                (double) cr[0] - 1, (double) cr[1] - 1,
                &(cw[0]), &(cw[1]),
                NULL
        ) != ES_NONE) {
                rt_raster_destroy(raster);
                PG_FREE_IF_COPY(pgraster, 0);
                elog(ERROR, "RASTER_rasterToWorldCoord: Could not compute longitude and latitude from pixel row and column");
                PG_RETURN_NULL();
        }
        rt_raster_destroy(raster);
        PG_FREE_IF_COPY(pgraster, 0);
 
        /* Build a tuple descriptor for our result type */
        if (get_call_result_type(fcinfo, NULL, &tupdesc) != TYPEFUNC_COMPOSITE) {
                ereport(ERROR, (
                        errcode(ERRCODE_FEATURE_NOT_SUPPORTED),
                        errmsg(
                                "function returning record called in context "
                                "that cannot accept type record"
                        )
                ));
        }
 
        BlessTupleDesc(tupdesc);
 
        values[0] = Float8GetDatum(cw[0]);
        values[1] = Float8GetDatum(cw[1]);
 
        memset(nulls, FALSE, sizeof(bool) * VALUES_LENGTH);
 
        /* build a tuple */
        tuple = heap_form_tuple(tupdesc, values, nulls);
 
        /* make the tuple into a datum */
        result = HeapTupleGetDatum(tuple);
 
        PG_RETURN_DATUM(result);
}

References ES_NONE, FALSE, FLT_NEQ, POSTGIS_RT_DEBUG, result, rt_raster_cell_to_geopoint(), rt_raster_deserialize(), rt_raster_destroy(), rt_raster_get_x_skew(), rt_raster_get_y_skew(), TRUE, and VALUES_LENGTH.

Here is the call graph for this function: