clusterintersecting_garray()

Datum clusterintersecting_garray

(

PG_FUNCTION_ARGS

)

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

References ARRAY2GEOS(), array_nelems_not_null(), cluster_intersecting(), cluster_within_distance_garray(), GEOS2POSTGIS(), LW_SUCCESS, lwgeom_geos_error(), PG_FUNCTION_INFO_V1(), and SRID_UNKNOWN.

Referenced by pgis_geometry_clusterintersecting_finalfn(), and polygonize_garray().

{
        Datum* result_array_data;
        ArrayType *array, *result;
        int is3d = 0;
        uint32 nelems, nclusters, i;
        GEOSGeometry **geos_inputs, **geos_results;
        int srid=SRID_UNKNOWN;

        /* Parameters used to construct a result array */
        int16 elmlen;
        bool elmbyval;
        char elmalign;

        /* Null array, null geometry (should be empty?) */
    if (PG_ARGISNULL(0))
        PG_RETURN_NULL();

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

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

        if ( nelems == 0 ) PG_RETURN_NULL();

    /* TODO short-circuit for one element? */

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

        geos_inputs = ARRAY2GEOS(array, nelems, &is3d, &srid);
        if(!geos_inputs)
        {
                PG_RETURN_NULL();
        }

        if (cluster_intersecting(geos_inputs, nelems, &geos_results, &nclusters) != LW_SUCCESS)
        {
                elog(ERROR, "clusterintersecting: Error performing clustering");
                PG_RETURN_NULL();
        }
        pfree(geos_inputs); /* don't need to destroy items because GeometryCollections have taken ownership */

        if (!geos_results) PG_RETURN_NULL();

        result_array_data = palloc(nclusters * sizeof(Datum));
        for (i=0; i<nclusters; ++i)
        {
                result_array_data[i] = PointerGetDatum(GEOS2POSTGIS(geos_results[i], is3d));
                GEOSGeom_destroy(geos_results[i]);
        }
        pfree(geos_results);

        get_typlenbyvalalign(array->elemtype, &elmlen, &elmbyval, &elmalign);
        result = (ArrayType*) construct_array(result_array_data, nclusters, array->elemtype, elmlen, elmbyval, elmalign);

        if (!result)
        {
                elog(ERROR, "clusterintersecting: Error constructing return-array");
                PG_RETURN_NULL();
        }

        PG_RETURN_POINTER(result);
}

Here is the call graph for this function:

Here is the caller graph for this function: