RASTER_rasterToWorldCoord()

Datum RASTER_rasterToWorldCoord

(

PG_FUNCTION_ARGS

)

Definition at line 691 of file rtpg_raster_properties.c.

References ES_NONE, FALSE, FLT_NEQ, PG_FUNCTION_INFO_V1(), POSTGIS_RT_DEBUG, rtrowdump::raster, RASTER_worldToRasterCoord(), 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.

{
        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) ? true : false;
        skewed[1] = FLT_NEQ(rt_raster_get_y_skew(raster), 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);
}

Here is the call graph for this function: