RASTER_worldToRasterCoord()

Datum RASTER_worldToRasterCoord

(

PG_FUNCTION_ARGS

)

Definition at line 785 of file rtpg_raster_properties.c.

References ES_NONE, FALSE, FLT_NEQ, PG_FUNCTION_INFO_V1(), POSTGIS_RT_DEBUG, rtrowdump::raster, RASTER_setSRID(), rt_raster_deserialize(), rt_raster_destroy(), rt_raster_geopoint_to_cell(), rt_raster_get_x_skew(), rt_raster_get_y_skew(), TRUE, and VALUES_LENGTH.

Referenced by RASTER_rasterToWorldCoord().

{
        rt_pgraster *pgraster = NULL;
        rt_raster raster = NULL;
        int i = 0;
        double cw[2] = {0};
        double _cr[2] = {0};
        int cr[2] = {0};
        bool skewed = false;

        TupleDesc tupdesc;
        Datum values[VALUES_LENGTH];
        bool nulls[VALUES_LENGTH];
        HeapTuple tuple;
        Datum result;

        POSTGIS_RT_DEBUG(3, "RASTER_worldToRasterCoord: 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_worldToRasterCoord: Could not deserialize raster");
                PG_RETURN_NULL();
        }

        /* raster skewed? */
        skewed = FLT_NEQ(rt_raster_get_x_skew(raster), 0) ? true : false;
        if (!skewed)
                skewed = FLT_NEQ(rt_raster_get_y_skew(raster), 0) ? true : false;

        /* longitude and latitude */
        for (i = 1; i <= 2; i++) {
                if (PG_ARGISNULL(i)) {
                        /* if skewed, parameter is required */
                        if (skewed) {
                                elog(NOTICE, "Latitude and longitude required for computing pixel row and column of a rotated raster");
                                rt_raster_destroy(raster);
                                PG_FREE_IF_COPY(pgraster, 0);
                                PG_RETURN_NULL();
                        }

                        continue;
                }

                cw[i - 1] = PG_GETARG_FLOAT8(i);
        }

        /* return pixel row and column values are 0-based */
        if (rt_raster_geopoint_to_cell(
                raster,
                cw[0], cw[1],
                &(_cr[0]), &(_cr[1]),
                NULL
        ) != ES_NONE) {
                rt_raster_destroy(raster);
                PG_FREE_IF_COPY(pgraster, 0);
                elog(ERROR, "RASTER_worldToRasterCoord: Could not compute pixel row and column from longitude and latitude");
                PG_RETURN_NULL();
        }
        rt_raster_destroy(raster);
        PG_FREE_IF_COPY(pgraster, 0);

        /* force to integer and add one to make 1-based */
        cr[0] = ((int) _cr[0]) + 1;
        cr[1] = ((int) _cr[1]) + 1;

        /* 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] = Int32GetDatum(cr[0]);
        values[1] = Int32GetDatum(cr[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:

Here is the caller graph for this function: