lwcompound_contains_point()

int lwcompound_contains_point	(	const LWCOMPOUND *	comp,
		const POINT2D *	pt
	)

Definition at line 145 of file lwcompound.c.

References LWCOMPOUND::geoms, LINETYPE, LW_BOUNDARY, LW_FALSE, LW_INSIDE, LW_OUTSIDE, lwerror(), lwgeom_as_lwcircstring(), lwgeom_as_lwline(), lwtype_name(), LWCOMPOUND::ngeoms, LWLINE::points, LWCIRCSTRING::points, ptarray_contains_point(), ptarray_contains_point_partial(), ptarrayarc_contains_point(), ptarrayarc_contains_point_partial(), and LWGEOM::type.

Referenced by lwgeom_contains_point().

{
        int i;
        LWLINE *lwline;
        LWCIRCSTRING *lwcirc;
        int wn = 0;
        int winding_number = 0;
        int result;

        for ( i = 0; i < comp->ngeoms; i++ )
        {
                LWGEOM *lwgeom = comp->geoms[i];
                if ( lwgeom->type == LINETYPE )
                {
                        lwline = lwgeom_as_lwline(lwgeom);
                        if ( comp->ngeoms == 1 )
                        {
                                return ptarray_contains_point(lwline->points, pt);
                        }
                        else
                        {
                                /* Don't check closure while doing p-i-p test */
                                result = ptarray_contains_point_partial(lwline->points, pt, LW_FALSE, &winding_number);
                        }
                }
                else
                {
                        lwcirc = lwgeom_as_lwcircstring(lwgeom);
                        if ( ! lwcirc ) {
                                lwerror("Unexpected component of type %s in compound curve", lwtype_name(lwgeom->type));
                                return 0;
                        }
                        if ( comp->ngeoms == 1 )
                        {
                                return ptarrayarc_contains_point(lwcirc->points, pt);
                        }
                        else
                        {
                                /* Don't check closure while doing p-i-p test */
                                result = ptarrayarc_contains_point_partial(lwcirc->points, pt, LW_FALSE, &winding_number);
                        }
                }

                /* Propogate boundary condition */
                if ( result == LW_BOUNDARY )
                        return LW_BOUNDARY;

                wn += winding_number;
        }

        /* Outside */
        if (wn == 0)
                return LW_OUTSIDE;

        /* Inside */
        return LW_INSIDE;
}

Here is the call graph for this function:

Here is the caller graph for this function: