2268 GBOX clip, subbox1, subbox2;
2274 double pivot = DBL_MAX;
2275 double center = DBL_MAX;
2279 if (!geom)
return 0;
2281 if (!box_in)
return 0;
2289 if ( width == 0.0 && height == 0.0 )
2320 for ( i = 0; i < incol->
ngeoms; i++ )
2334 if ( depth > maxdepth )
2343 if (nvertices == 0)
return 0;
2346 if (nvertices <= maxvertices)
2352 split_ordinate = (width > height) ? 0 : 1;
2353 if (split_ordinate == 0)
2354 center = (clip.
xmin + clip.
xmax) / 2;
2356 center = (clip.
ymin + clip.
ymax) / 2;
2361 double ring_area = 0;
2362 double pivot_eps = DBL_MAX;
2363 double pt_eps = DBL_MAX;
2372 for (i = 1; i < lwpoly->
nrings; i++)
2375 if (current_ring_area >= ring_area)
2377 ring_area = current_ring_area;
2383 pa = lwpoly->
rings[ring_to_trim];
2386 for (i = 0; i < pa->
npoints; i++)
2389 if (split_ordinate == 0)
2393 pt_eps = fabs(pt - center);
2394 if (pivot_eps > pt_eps)
2407 if (split_ordinate == 0)
2412 subbox1.
xmax = subbox2.
xmin = center;
2419 subbox1.
ymax = subbox2.
ymin = center;
void gbox_duplicate(const GBOX *original, GBOX *duplicate)
Copy the values of original GBOX into duplicate.
LWGEOM * lwgeom_intersection(const LWGEOM *geom1, const LWGEOM *geom2)
#define POINTTYPE
LWTYPE numbers, used internally by PostGIS.
#define POLYHEDRALSURFACETYPE
const char * lwtype_name(uint8_t type)
Return the type name string associated with a type number (e.g.
LWCOLLECTION * lwcollection_add_lwgeom(LWCOLLECTION *col, const LWGEOM *geom)
Appends geom to the collection managed by col.
#define LW_TRUE
Return types for functions with status returns.
LWPOLY * lwpoly_construct_envelope(int srid, double x1, double y1, double x2, double y2)
const POINT2D * getPoint2d_cp(const POINTARRAY *pa, uint32_t n)
Returns a POINT2D pointer into the POINTARRAY serialized_ptlist, suitable for reading from.
double ptarray_signed_area(const POINTARRAY *pa)
Returns the area in cartesian units.
#define FP_TOLERANCE
Floating point comparators.
void lwgeom_simplify_in_place(LWGEOM *geom, double epsilon, int preserve_collapsed)
int lwgeom_is_collection(const LWGEOM *geom)
Determine whether a LWGEOM can contain sub-geometries or not.
LWGEOM * lwgeom_clone_deep(const LWGEOM *lwgeom)
Deep-clone an LWGEOM object.
uint32_t lwgeom_count_vertices(const LWGEOM *geom)
Count points in an LWGEOM.
int lwgeom_dimension(const LWGEOM *geom)
For an LWGEOM, returns 0 for points, 1 for lines, 2 for polygons, 3 for volume, and the max dimension...
const GBOX * lwgeom_get_bbox(const LWGEOM *lwg)
Get a non-empty geometry bounding box, computing and caching it if not already there.
int lwgeom_is_empty(const LWGEOM *geom)
Return true or false depending on whether a geometry is an "empty" geometry (no vertices members)
void lwgeom_free(LWGEOM *lwgeom)
static int lwgeom_subdivide_recursive(const LWGEOM *geom, uint8_t dimension, uint32_t maxvertices, uint32_t depth, LWCOLLECTION *col)
void lwerror(const char *fmt,...)
Write a notice out to the error handler.
static double pivot(double *left, double *right)