lwgeom_clip_by_rect()

LWGEOM* lwgeom_clip_by_rect	(	const LWGEOM *	geom1,
		double	x1,
		double	y1,
		double	x2,
		double	y2
	)

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

{
        LWGEOM *result;
        GEOSGeometry *g1, *g3;
        int is3d;
 
        /* A.Intersection(Empty) == Empty */
        if ( lwgeom_is_empty(geom1) )
                return lwgeom_clone_deep(geom1);
 
        is3d = FLAGS_GET_Z(geom1->flags);
 
        initGEOS(lwnotice, lwgeom_geos_error);
 
        if (!(g1 = LWGEOM2GEOS(geom1, AUTOFIX)))
                GEOS_FAIL_DEBUG();
 
        if (!(g3 = GEOSClipByRect(g1, x1, y1, x2, y2)))
                GEOS_FREE_AND_FAIL_DEBUG(g1);
 
        GEOS_FREE(g1);
        result = GEOS2LWGEOM(g3, is3d);
        GEOS_FREE(g3);
 
        if (!result)
                GEOS_FAIL_DEBUG();
 
        result->srid = geom1->srid;
 
        return result;
}

References AUTOFIX, LWGEOM::flags, FLAGS_GET_Z, GEOS2LWGEOM(), GEOS_FAIL_DEBUG, GEOS_FREE, GEOS_FREE_AND_FAIL_DEBUG, LWGEOM2GEOS(), lwgeom_clone_deep(), lwgeom_geos_error(), lwgeom_is_empty(), lwnotice(), result, and LWGEOM::srid.

Referenced by mvt_unsafe_clip_by_box(), ST_ClipByBox2d(), and test_lwgeom_clip_by_rect().

Here is the call graph for this function:

Here is the caller graph for this function: