LWGEOM*
lwmpoint_remove_repeated_points(LWMPOINT *mpoint)
{
uint32_t nnewgeoms;
uint32_t i, j;
LWGEOM **newgeoms;
newgeoms = lwalloc(sizeof(LWGEOM *)*mpoint->ngeoms);
nnewgeoms = 0;
for (i=0; i<mpoint->ngeoms; ++i)
{
/* Brute force, may be optimized by building an index */
int seen=0;
for (j=0; j<nnewgeoms; ++j)
{
if ( lwpoint_same((LWPOINT*)newgeoms[j],
(LWPOINT*)mpoint->geoms[i]) )
{
seen=1;
break;
}
}
if ( seen ) continue;
newgeoms[nnewgeoms++] = (LWGEOM*)lwpoint_clone(mpoint->geoms[i]);
}
return (LWGEOM*)lwcollection_construct(mpoint->type,
mpoint->srid, mpoint->bbox ? gbox_copy(mpoint->bbox) : NULL,
nnewgeoms, newgeoms);
}
/**
* @brief geom1 same as geom2
* iff
* + have same type
* + have same # objects
* + have same bvol
* + each object in geom1 has a corresponding object in geom2 (see above)
* @param lwgeom1
* @param lwgeom2
*/
char
lwgeom_same(const LWGEOM *lwgeom1, const LWGEOM *lwgeom2)
{
LWDEBUGF(2, "lwgeom_same(%s, %s) called",
lwtype_name(lwgeom1->type),
lwtype_name(lwgeom2->type));
if ( lwgeom1->type != lwgeom2->type )
{
LWDEBUG(3, " type differ");
return LW_FALSE;
}
if ( FLAGS_GET_ZM(lwgeom1->flags) != FLAGS_GET_ZM(lwgeom2->flags) )
{
LWDEBUG(3, " ZM flags differ");
return LW_FALSE;
}
/* Check boxes if both already computed */
if ( lwgeom1->bbox && lwgeom2->bbox )
{
/*lwnotice("bbox1:%p, bbox2:%p", lwgeom1->bbox, lwgeom2->bbox);*/
if ( ! gbox_same(lwgeom1->bbox, lwgeom2->bbox) )
{
LWDEBUG(3, " bounding boxes differ");
return LW_FALSE;
}
}
/* geoms have same type, invoke type-specific function */
switch (lwgeom1->type)
{
case POINTTYPE:
return lwpoint_same((LWPOINT *)lwgeom1,
(LWPOINT *)lwgeom2);
case LINETYPE:
return lwline_same((LWLINE *)lwgeom1,
(LWLINE *)lwgeom2);
case POLYGONTYPE:
return lwpoly_same((LWPOLY *)lwgeom1,
(LWPOLY *)lwgeom2);
case TRIANGLETYPE:
return lwtriangle_same((LWTRIANGLE *)lwgeom1,
(LWTRIANGLE *)lwgeom2);
case CIRCSTRINGTYPE:
return lwcircstring_same((LWCIRCSTRING *)lwgeom1,
(LWCIRCSTRING *)lwgeom2);
case MULTIPOINTTYPE:
case MULTILINETYPE:
case MULTIPOLYGONTYPE:
case MULTICURVETYPE:
case MULTISURFACETYPE:
case COMPOUNDTYPE:
case CURVEPOLYTYPE:
case POLYHEDRALSURFACETYPE:
case TINTYPE:
case COLLECTIONTYPE:
return lwcollection_same((LWCOLLECTION *)lwgeom1,
(LWCOLLECTION *)lwgeom2);
default:
lwerror("lwgeom_same: unsupported geometry type: %s",
lwtype_name(lwgeom1->type));
return LW_FALSE;
}
}
