lwgeom_clip_to_ordinate_range()

LWCOLLECTION * lwgeom_clip_to_ordinate_range	(	const LWGEOM *	lwin,
		char	ordinate,
		double	from,
		double	to,
		double	offset
	)

Given a geometry clip based on the from/to range of one of its ordinates (x, y, z, m).

Use for m- and z- clipping.

Definition at line 823 of file lwlinearreferencing.c.

{
        LWCOLLECTION *out_col;
        LWCOLLECTION *out_offset;
        uint32_t i;
 
        /* Ensure 'from' is less than 'to'. */
        if (to < from)
        {
                double t = from;
                from = to;
                to = t;
        }
 
        if (!lwin)
                lwerror("lwgeom_clip_to_ordinate_range: null input geometry!");
 
        switch (lwin->type)
        {
        case LINETYPE:
                out_col = lwline_clip_to_ordinate_range((LWLINE *)lwin, ordinate, from, to);
                break;
        case MULTIPOINTTYPE:
                out_col = lwmpoint_clip_to_ordinate_range((LWMPOINT *)lwin, ordinate, from, to);
                break;
        case POINTTYPE:
                out_col = lwpoint_clip_to_ordinate_range((LWPOINT *)lwin, ordinate, from, to);
                break;
        case POLYGONTYPE:
                out_col = lwpoly_clip_to_ordinate_range((LWPOLY *)lwin, ordinate, from, to);
                break;
        case TRIANGLETYPE:
                out_col = lwtriangle_clip_to_ordinate_range((LWTRIANGLE *)lwin, ordinate, from, to);
                break;
        case TINTYPE:
        case MULTILINETYPE:
        case MULTIPOLYGONTYPE:
        case COLLECTIONTYPE:
        case POLYHEDRALSURFACETYPE:
                out_col = lwcollection_clip_to_ordinate_range((LWCOLLECTION *)lwin, ordinate, from, to);
                break;
        default:
                lwerror("This function does not accept %s geometries.", lwtype_name(lwin->type));
                return NULL;
        }
 
        /* Stop if result is NULL */
        if (!out_col)
                lwerror("lwgeom_clip_to_ordinate_range clipping routine returned NULL");
 
        /* Return if we aren't going to offset the result */
        if (FP_IS_ZERO(offset) || lwgeom_is_empty(lwcollection_as_lwgeom(out_col)))
                return out_col;
 
        /* Construct a collection to hold our outputs. */
        /* Things get ugly: GEOS offset drops Z's and M's so we have to drop ours */
        out_offset = lwcollection_construct_empty(MULTILINETYPE, lwin->srid, 0, 0);
 
        /* Try and offset the linear portions of the return value */
        for (i = 0; i < out_col->ngeoms; i++)
        {
                int type = out_col->geoms[i]->type;
                if (type == POINTTYPE)
                {
                        lwnotice("lwgeom_clip_to_ordinate_range cannot offset a clipped point");
                        continue;
                }
                else if (type == LINETYPE)
                {
                        /* lwgeom_offsetcurve(line, offset, quadsegs, joinstyle (round), mitrelimit) */
                        LWGEOM *lwoff = lwgeom_offsetcurve(out_col->geoms[i], offset, 8, 1, 5.0);
                        if (!lwoff)
                        {
                                lwerror("lwgeom_offsetcurve returned null");
                        }
                        lwcollection_add_lwgeom(out_offset, lwoff);
                }
                else
                {
                        lwerror("lwgeom_clip_to_ordinate_range found an unexpected type (%s) in the offset routine",
                                lwtype_name(type));
                }
        }
 
        return out_offset;
}

References COLLECTIONTYPE, FP_IS_ZERO, LWCOLLECTION::geoms, LINETYPE, lwcollection_add_lwgeom(), lwcollection_as_lwgeom(), lwcollection_clip_to_ordinate_range(), lwcollection_construct_empty(), lwerror(), lwgeom_is_empty(), lwgeom_offsetcurve(), lwline_clip_to_ordinate_range(), lwmpoint_clip_to_ordinate_range(), lwnotice(), lwpoint_clip_to_ordinate_range(), lwpoly_clip_to_ordinate_range(), lwtriangle_clip_to_ordinate_range(), lwtype_name(), MULTILINETYPE, MULTIPOINTTYPE, MULTIPOLYGONTYPE, LWCOLLECTION::ngeoms, POINTTYPE, POLYGONTYPE, POLYHEDRALSURFACETYPE, LWGEOM::srid, TINTYPE, TRIANGLETYPE, LWGEOM::type, and LWLINE::type.

Referenced by lwcollection_clip_to_ordinate_range(), lwgeom_locate_between(), lwgeom_solid_contains_lwgeom(), ST_LocateBetween(), ST_LocateBetweenElevations(), test_lwline_clip(), test_lwline_clip_big(), test_lwmline_clip(), test_lwpoly_clip(), and test_lwtriangle_clip().

Here is the call graph for this function:

Here is the caller graph for this function: