polygonize_garray()

Datum polygonize_garray

(

PG_FUNCTION_ARGS

)

Definition at line 3132 of file postgis/lwgeom_geos.c.

References ARRAY2GEOS(), array_nelems_not_null(), clusterintersecting_garray(), GEOS2POSTGIS(), lwgeom_geos_error(), PG_FUNCTION_INFO_V1(), and SRID_UNKNOWN.

Referenced by GEOSnoop(), and pgis_geometry_polygonize_finalfn().

{
        ArrayType *array;
        int is3d = 0;
        uint32 nelems, i;
        GSERIALIZED *result;
        GEOSGeometry *geos_result;
        const GEOSGeometry **vgeoms;
        int srid=SRID_UNKNOWN;
#if POSTGIS_DEBUG_LEVEL >= 3
        static int call=1;
#endif

#if POSTGIS_DEBUG_LEVEL >= 3
        call++;
#endif

        if (PG_ARGISNULL(0))
                PG_RETURN_NULL();

        array = PG_GETARG_ARRAYTYPE_P(0);
        nelems = array_nelems_not_null(array);

        if (nelems == 0)
                PG_RETURN_NULL();

        POSTGIS_DEBUGF(3, "polygonize_garray: number of non-null elements: %d", nelems);

        /* Ok, we really need geos now ;) */
        initGEOS(lwpgnotice, lwgeom_geos_error);

        vgeoms = (const GEOSGeometry**) ARRAY2GEOS(array, nelems, &is3d, &srid);

        POSTGIS_DEBUG(3, "polygonize_garray: invoking GEOSpolygonize");

        geos_result = GEOSPolygonize(vgeoms, nelems);

        POSTGIS_DEBUG(3, "polygonize_garray: GEOSpolygonize returned");

        for (i=0; i<nelems; ++i) GEOSGeom_destroy((GEOSGeometry *)vgeoms[i]);
        pfree(vgeoms);

        if ( ! geos_result ) PG_RETURN_NULL();

        GEOSSetSRID(geos_result, srid);
        result = GEOS2POSTGIS(geos_result, is3d);
        GEOSGeom_destroy(geos_result);
        if ( result == NULL )
        {
                elog(ERROR, "%s returned an error", __func__);
                PG_RETURN_NULL(); /*never get here */
        }

        PG_RETURN_POINTER(result);
}

Here is the call graph for this function:

Here is the caller graph for this function: