lwmpoint_remove_repeated_points()

LWGEOM* lwmpoint_remove_repeated_points	(	const LWMPOINT *	in,
		double	tolerance
	)

Definition at line 92 of file lwmpoint.c.

References LWMPOINT::bbox, gbox_copy(), LWMPOINT::geoms, lwalloc(), lwcollection_construct(), lwgeom_clone_deep(), lwgeom_mindistance2d(), LWMPOINT::ngeoms, LWMPOINT::srid, and LWMPOINT::type.

Referenced by lwgeom_remove_repeated_points().

{
        uint32_t nnewgeoms;
        uint32_t i, j;
        LWGEOM **newgeoms;
        LWGEOM *lwpt1, *lwpt2;

        newgeoms = lwalloc(sizeof(LWGEOM *)*mpoint->ngeoms);
        nnewgeoms = 0;
        for (i=0; i<mpoint->ngeoms; ++i)
        {
                lwpt1 = (LWGEOM*)mpoint->geoms[i];
                /* Brute force, may be optimized by building an index */
                int seen=0;
                for (j=0; j<nnewgeoms; ++j)
                {
                        lwpt2 = (LWGEOM*)newgeoms[j];
                        if ( lwgeom_mindistance2d(lwpt1, lwpt2) <= tolerance )
                        {
                                seen=1;
                                break;
                        }
                }
                if ( seen ) continue;
                newgeoms[nnewgeoms++] = lwgeom_clone_deep(lwpt1);
        }

        return (LWGEOM*)lwcollection_construct(mpoint->type,
                                               mpoint->srid,
                                                                                   mpoint->bbox ? gbox_copy(mpoint->bbox) : NULL,
                                               nnewgeoms, newgeoms);

}

Here is the call graph for this function:

Here is the caller graph for this function: