478 ArrayIterator iterator;
483 int nelems = 0, geoms_size = 0, curgeom = 0,
count = 0;
487 GEOSGeometry *g = NULL;
488 GEOSGeometry *g_union = NULL;
489 GEOSGeometry **geoms = NULL;
496 if ( PG_ARGISNULL(0) )
499 array = PG_GETARG_ARRAYTYPE_P(0);
500 nelems = ArrayGetNItems(ARR_NDIM(array), ARR_DIMS(array));
503 if ( nelems == 0 ) PG_RETURN_NULL();
506 iterator = array_create_iterator(array, 0, NULL);
507 while (array_iterate(iterator, &
value, &isnull))
510 if (isnull)
continue;
513 array_free_iterator(iterator);
521 if (
count == 1 && nelems == 1 )
523 #if __GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 6)
524 #pragma GCC diagnostic push
525 #pragma GCC diagnostic ignored "-Wsign-compare"
527 PG_RETURN_POINTER((
GSERIALIZED *)(ARR_DATA_PTR(array)));
528 #if __GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 6)
529 #pragma GCC diagnostic pop
540 geoms = palloc(
sizeof(GEOSGeometry*) * geoms_size);
546 iterator = array_create_iterator(array, 0, NULL);
547 while (array_iterate(iterator, &
value, &isnull))
552 if (isnull)
continue;
570 if (gser_type > empty_type)
572 empty_type = gser_type;
573 POSTGIS_DEBUGF(4,
"empty_type = %d gser_type = %d", empty_type, gser_type);
584 "One of the geometries in the set "
585 "could not be converted to GEOS");
589 if ( curgeom == geoms_size )
592 geoms = repalloc( geoms,
sizeof(GEOSGeometry*) * geoms_size );
600 array_free_iterator(iterator);
608 g = GEOSGeom_createCollection(GEOS_GEOMETRYCOLLECTION, geoms, curgeom);
609 if (!g)
HANDLE_GEOS_ERROR(
"Could not create GEOS COLLECTION from geometry array");
611 g_union = GEOSUnaryUnion(g);
615 GEOSSetSRID(g_union, srid);
617 GEOSGeom_destroy(g_union);
623 if ( empty_type > 0 )
640 PG_RETURN_POINTER(gser_out);
void gserialized_error_if_srid_mismatch_reference(const GSERIALIZED *g1, const int32_t srid2, const char *funcname)
int32_t gserialized_get_srid(const GSERIALIZED *g)
Extract the SRID from the serialized form (it is packed into three bytes so this is a handy function)...
int gserialized_is_empty(const GSERIALIZED *g)
Check if a GSERIALIZED is empty without deserializing first.
int gserialized_has_z(const GSERIALIZED *g)
Check if a GSERIALIZED has a Z ordinate.
uint32_t gserialized_get_type(const GSERIALIZED *g)
Extract the geometry type from the serialized form (it hides in the anonymous data area,...
void lwgeom_geos_error(const char *fmt,...)
#define SRID_UNKNOWN
Unknown SRID value.
LWGEOM * lwgeom_construct_empty(uint8_t type, int32_t srid, char hasz, char hasm)
#define HANDLE_GEOS_ERROR(label)
GSERIALIZED * GEOS2POSTGIS(GEOSGeom geom, char want3d)
GEOSGeometry * POSTGIS2GEOS(const GSERIALIZED *pglwgeom)