938{
942 GEOSGeometry *g1, *g3;
943
945
947 {
950 }
951
953
955
956 g3 = GEOSGetCentroid(g1);
957
959 GEOSSetSRID(g3, srid);
960
963
965
967}
char result[OUT_DOUBLE_BUFFER_SIZE]
GEOSGeometry * LWGEOM2GEOS(const LWGEOM *lwgeom, uint8_t autofix)
#define GEOS_FREE_AND_FAIL(...)
void lwgeom_geos_error(const char *fmt,...)
void(*) LWGEOM GEOS2LWGEOM)(const GEOSGeometry *geom, uint8_t want3d)
LWGEOM * lwpoint_as_lwgeom(const LWPOINT *obj)
#define FLAGS_GET_Z(flags)
LWPOINT * lwpoint_construct_empty(int32_t srid, char hasz, char hasm)
int lwgeom_has_m(const LWGEOM *geom)
Return LW_TRUE if geometry has M ordinates.
void lwnotice(const char *fmt,...) __attribute__((format(printf
Write a notice out to the notice handler.
static int lwgeom_is_empty(const LWGEOM *geom)
Return true or false depending on whether a geometry is an "empty" geometry (no vertices members)