1007{
1011 GEOSGeometry *g1, *g3;
1012
1014
1016 {
1019 }
1020
1022
1024
1025 g3 = GEOSPointOnSurface(g1);
1026
1028 GEOSSetSRID(g3, srid);
1029
1032
1034
1036}
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)