RASTER_convex_hull()

Datum RASTER_convex_hull

(

PG_FUNCTION_ARGS

)

Definition at line 119 of file rtpg_geometry.c.

References ES_NONE, FALSE, gserialized_from_lwgeom(), lwgeom_free(), pixval::nband, PG_FUNCTION_INFO_V1(), rtrowdump::raster, RASTER_dumpAsPolygons(), rt_raster_deserialize(), rt_raster_destroy(), rt_raster_get_convex_hull(), rt_raster_get_perimeter(), rt_raster_has_band(), and TRUE.

Referenced by RASTER_envelope().

{
        rt_pgraster *pgraster;
        rt_raster raster;
        LWGEOM *geom = NULL;
        GSERIALIZED* gser = NULL;
        size_t gser_size;
        int err = ES_NONE;

        bool minhull = FALSE;

        if (PG_ARGISNULL(0))
                PG_RETURN_NULL();

        /* # of args */
        if (PG_NARGS() > 1)
                minhull = TRUE;

        if (!minhull) {
                pgraster = (rt_pgraster *) PG_DETOAST_DATUM_SLICE(PG_GETARG_DATUM(0), 0, sizeof(struct rt_raster_serialized_t));
                raster = rt_raster_deserialize(pgraster, TRUE);
        }
        else {
                pgraster = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
                raster = rt_raster_deserialize(pgraster, FALSE);
        }

        if (!raster) {
                PG_FREE_IF_COPY(pgraster, 0);
                elog(ERROR, "RASTER_convex_hull: Could not deserialize raster");
                PG_RETURN_NULL();
        }

        if (!minhull)
                err = rt_raster_get_convex_hull(raster, &geom);
        else {
                int nband = -1;

                /* get arg 1 */
                if (!PG_ARGISNULL(1)) {
                        nband = PG_GETARG_INT32(1);
                        if (!rt_raster_has_band(raster, nband - 1)) {
                                elog(NOTICE, "Invalid band index (must use 1-based). Returning NULL");
                                rt_raster_destroy(raster);
                                PG_FREE_IF_COPY(pgraster, 0);
                                PG_RETURN_NULL();
                        }
                        nband = nband - 1;
                }

                err = rt_raster_get_perimeter(raster, nband, &geom);
        }

        rt_raster_destroy(raster);
        PG_FREE_IF_COPY(pgraster, 0);

        if (err != ES_NONE) {
                elog(ERROR, "RASTER_convex_hull: Could not get raster's convex hull");
                PG_RETURN_NULL();
        }
        else if (geom == NULL) {
                elog(NOTICE, "Raster's convex hull is NULL");
                PG_RETURN_NULL();
        }

        gser = gserialized_from_lwgeom(geom, &gser_size);
        lwgeom_free(geom);

        SET_VARSIZE(gser, gser_size);
        PG_RETURN_POINTER(gser);
}

Here is the call graph for this function:

Here is the caller graph for this function: