clusterintersecting_garray()

Datum clusterintersecting_garray

(

PG_FUNCTION_ARGS

)

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

{
        Datum* result_array_data;
        ArrayType *array, *result;
        int is3d = 0;
        uint32 nelems, nclusters, i;
        GEOSGeometry **geos_inputs, **geos_results;
        int32_t 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]);
        }
        lwfree(geos_results);
 
        get_typlenbyvalalign(array->elemtype, &elmlen, &elmbyval, &elmalign);
        result = 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);
}

References ARRAY2GEOS(), array_nelems_not_null(), cluster_intersecting(), GEOS2POSTGIS(), LW_SUCCESS, lwfree(), lwgeom_geos_error(), result, and SRID_UNKNOWN.

Referenced by pgis_geometry_clusterintersecting_finalfn().

Here is the call graph for this function:

Here is the caller graph for this function: