529 ArrayType *arr_geoms = NULL;
530 ArrayType *arr_ids = NULL;
531 int num_geoms, num_ids, i = 0;
533 ArrayIterator iter_geoms, iter_ids;
534 bool null_geom, null_id;
535 Datum val_geom, val_id;
537 int is_homogeneous =
true;
538 uint32_t subtype = 0;
542 int64_t *idlist = NULL;
551 if ( PG_NARGS() < 2 || PG_ARGISNULL(0) || PG_ARGISNULL(1) )
554 arr_geoms = PG_GETARG_ARRAYTYPE_P(0);
555 arr_ids = PG_GETARG_ARRAYTYPE_P(1);
557 num_geoms = ArrayGetNItems(ARR_NDIM(arr_geoms), ARR_DIMS(arr_geoms));
558 num_ids = ArrayGetNItems(ARR_NDIM(arr_ids), ARR_DIMS(arr_ids));
560 if ( num_geoms != num_ids )
562 elog(ERROR,
"size of geometry[] and integer[] arrays must match");
569 iter_geoms = array_create_iterator(arr_geoms, 0, NULL);
570 iter_ids = array_create_iterator(arr_ids, 0, NULL);
572 while( array_iterate(iter_geoms, &val_geom, &null_geom) &&
573 array_iterate(iter_ids, &val_id, &null_id) )
578 if ( null_geom || null_id )
580 elog(NOTICE,
"ST_AsTWKB skipping NULL entry at position %d", i);
585 uid = DatumGetInt64(val_id);
595 idlist = palloc0(num_geoms *
sizeof(int64_t));
601 elog(ERROR,
"Geometries have different dimensionality");
602 PG_FREE_IF_COPY(arr_geoms, 0);
603 PG_FREE_IF_COPY(arr_ids, 1);
614 is_homogeneous =
false;
622 array_free_iterator(iter_geoms);
623 array_free_iterator(iter_ids);
627 elog(NOTICE,
"No valid geometry - id pairs found");
628 PG_FREE_IF_COPY(arr_geoms, 0);
629 PG_FREE_IF_COPY(arr_ids, 1);
632 if ( is_homogeneous )
641 if ( PG_NARGS() > 2 && ! PG_ARGISNULL(2) )
642 sp.precision_xy = PG_GETARG_INT32(2);
645 if ( PG_NARGS() > 3 && ! PG_ARGISNULL(3) )
646 sp.precision_z = PG_GETARG_INT32(3);
649 if ( PG_NARGS() > 4 && ! PG_ARGISNULL(4) )
650 sp.precision_m = PG_GETARG_INT32(4);
656 if ( PG_NARGS() > 5 && ! PG_ARGISNULL(5) && PG_GETARG_BOOL(5) )
660 if ( PG_NARGS() > 6 && ! PG_ARGISNULL(6) && PG_GETARG_BOOL(6) )
666 sp.precision_xy, sp.precision_z, sp.precision_m,
670 result = palloc(twkb_size + VARHDRSZ);
671 memcpy(VARDATA(result), twkb, twkb_size);
672 SET_VARSIZE(result, twkb_size + VARHDRSZ);
678 PG_FREE_IF_COPY(arr_geoms, 0);
679 PG_FREE_IF_COPY(arr_ids, 1);
681 PG_RETURN_BYTEA_P(result);
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
LWGEOM * lwcollection_as_lwgeom(const LWCOLLECTION *obj)
uint32_t lwtype_get_collectiontype(uint8_t type)
Given an lwtype number, what homogeneous collection can hold it?
int32_t lwgeom_get_srid(const LWGEOM *geom)
Return SRID number.
uint8_t * lwgeom_to_twkb_with_idlist(const LWGEOM *geom, int64_t *idlist, uint8_t variant, int8_t precision_xy, int8_t precision_z, int8_t precision_m, size_t *twkb_size)
Convert LWGEOM to a char* in TWKB format.
int lwgeom_has_z(const LWGEOM *geom)
Return LW_TRUE if geometry has Z ordinates.
LWCOLLECTION * lwcollection_construct_empty(uint8_t type, int32_t srid, char hasz, char hasm)
void lwcollection_free(LWCOLLECTION *col)
#define TWKB_DEFAULT_PRECISION
LWCOLLECTION * lwcollection_add_lwgeom(LWCOLLECTION *col, const LWGEOM *geom)
Appends geom to the collection managed by col.
int lwgeom_has_m(const LWGEOM *geom)
Return LW_TRUE if geometry has M ordinates.
static uint32_t lwgeom_get_type(const LWGEOM *geom)
Return LWTYPE number.