/**
* Deep-clone an #LWGEOM object. #POINTARRAY <em>are</em> copied.
*/
LWGEOM *
lwgeom_clone_deep(const LWGEOM *lwgeom)
{
LWDEBUGF(2, "lwgeom_clone called with %p, %s",
lwgeom, lwtype_name(lwgeom->type));
switch (lwgeom->type)
{
case POINTTYPE:
case LINETYPE:
case CIRCSTRINGTYPE:
case TRIANGLETYPE:
return (LWGEOM *)lwline_clone_deep((LWLINE *)lwgeom);
case POLYGONTYPE:
return (LWGEOM *)lwpoly_clone_deep((LWPOLY *)lwgeom);
case COMPOUNDTYPE:
case CURVEPOLYTYPE:
case MULTICURVETYPE:
case MULTISURFACETYPE:
case MULTIPOINTTYPE:
case MULTILINETYPE:
case MULTIPOLYGONTYPE:
case POLYHEDRALSURFACETYPE:
case TINTYPE:
case COLLECTIONTYPE:
return (LWGEOM *)lwcollection_clone_deep((LWCOLLECTION *)lwgeom);
default:
lwerror("lwgeom_clone_deep: Unknown geometry type: %s", lwtype_name(lwgeom->type));
return NULL;
}
}
static LWGEOM*
lwline_split_by_point(const LWLINE* lwline_in, const LWPOINT* blade_in)
{
LWMLINE* out;
out = lwmline_construct_empty(lwline_in->srid,
FLAGS_GET_Z(lwline_in->flags),
FLAGS_GET_M(lwline_in->flags));
if ( lwline_split_by_point_to(lwline_in, blade_in, out) < 2 )
{
lwmline_add_lwline(out, lwline_clone_deep(lwline_in));
}
/* Turn multiline into collection */
out->type = COLLECTIONTYPE;
return (LWGEOM*)out;
}
static LWGEOM*
lwline_split_by_mpoint(const LWLINE* lwline_in, const LWMPOINT* mp)
{
LWMLINE* out;
int i, j;
out = lwmline_construct_empty(lwline_in->srid,
FLAGS_GET_Z(lwline_in->flags),
FLAGS_GET_M(lwline_in->flags));
lwmline_add_lwline(out, lwline_clone_deep(lwline_in));
for (i=0; i<mp->ngeoms; ++i)
{
for (j=0; j<out->ngeoms; ++j)
{
lwline_in = out->geoms[j];
LWPOINT *blade_in = mp->geoms[i];
int ret = lwline_split_by_point_to(lwline_in, blade_in, out);
if ( 2 == ret )
{
/* the point splits this line,
* 2 splits were added to collection.
* We'll move the latest added into
* the slot of the current one.
*/
lwline_free(out->geoms[j]);
out->geoms[j] = out->geoms[--out->ngeoms];
}
}
}
/* Turn multiline into collection */
out->type = COLLECTIONTYPE;
return (LWGEOM*)out;
}
