3053 Datum* result_array_data;
3054 ArrayType *array, *
result;
3056 uint32 nelems, nclusters, i;
3057 GEOSGeometry **geos_inputs, **geos_results;
3066 if (PG_ARGISNULL(0))
3069 array = PG_GETARG_ARRAYTYPE_P(0);
3072 POSTGIS_DEBUGF(3,
"clusterintersecting_garray: number of non-null elements: %d", nelems);
3074 if ( nelems == 0 ) PG_RETURN_NULL();
3081 geos_inputs =
ARRAY2GEOS(array, nelems, &is3d, &srid);
3089 elog(ERROR,
"clusterintersecting: Error performing clustering");
3094 if (!geos_results) PG_RETURN_NULL();
3096 result_array_data = palloc(nclusters *
sizeof(Datum));
3097 for (i=0; i<nclusters; ++i)
3099 result_array_data[i] = PointerGetDatum(
GEOS2POSTGIS(geos_results[i], is3d));
3100 GEOSGeom_destroy(geos_results[i]);
3104 get_typlenbyvalalign(array->elemtype, &elmlen, &elmbyval, &elmalign);
3105 result = construct_array(result_array_data, nclusters, array->elemtype, elmlen, elmbyval, elmalign);
3109 elog(ERROR,
"clusterintersecting: Error constructing return-array");
3113 PG_RETURN_POINTER(
result);
char result[OUT_DOUBLE_BUFFER_SIZE]
void lwgeom_geos_error(const char *fmt,...)
int cluster_intersecting(GEOSGeometry **geoms, uint32_t num_geoms, GEOSGeometry ***clusterGeoms, uint32_t *num_clusters)
Takes an array of GEOSGeometry* and constructs an array of GEOSGeometry*, where each element in the c...
#define SRID_UNKNOWN
Unknown SRID value.
uint32_t array_nelems_not_null(ArrayType *array)
GSERIALIZED * GEOS2POSTGIS(GEOSGeom geom, char want3d)
GEOSGeometry ** ARRAY2GEOS(ArrayType *array, uint32_t nelems, int *is3d, int *srid)