lwgeom_snap()

LWGEOM* lwgeom_snap	(	const LWGEOM *	geom1,
		const LWGEOM *	geom2,
		double	tolerance
	)

Snap vertices and segments of a geometry to another using a given tolerance.

Parameters

geom1	the geometry to snap
geom2	the geometry to snap to
tolerance	the distance under which vertices and segments are snapped

Definition at line 1514 of file liblwgeom/lwgeom_geos.c.

{
        LWGEOM* result;
        int32_t srid = RESULT_SRID(geom1, geom2);
        uint8_t is3d = (FLAGS_GET_Z(geom1->flags) || FLAGS_GET_Z(geom2->flags));
        GEOSGeometry *g1, *g2, *g3;
 
        if (srid == SRID_INVALID) return NULL;
 
        initGEOS(lwnotice, lwgeom_geos_error);
 
        if (!(g1 = LWGEOM2GEOS(geom1, AUTOFIX))) GEOS_FAIL();
        if (!(g2 = LWGEOM2GEOS(geom2, AUTOFIX))) GEOS_FREE_AND_FAIL(g1);
 
        g3 = GEOSSnap(g1, g2, tolerance);
 
        if (!g3) GEOS_FREE_AND_FAIL(g1, g2);
        GEOSSetSRID(g3, srid);
 
        if (!(result = GEOS2LWGEOM(g3, is3d)))
                GEOS_FREE_AND_FAIL(g1, g2, g3);
 
        GEOS_FREE(g1, g2, g3);
        return result;
}

References AUTOFIX, LWGEOM::flags, FLAGS_GET_Z, GEOS2LWGEOM(), GEOS_FAIL, GEOS_FREE, GEOS_FREE_AND_FAIL, LWGEOM2GEOS(), lwgeom_geos_error(), lwnotice(), result, RESULT_SRID, and SRID_INVALID.

Referenced by _lwt_toposnap(), and ST_Snap().

Here is the call graph for this function:

Here is the caller graph for this function: