point_in_multipolygon_rtree()

int point_in_multipolygon_rtree	(	RTREE_NODE **	root,
		int	polyCount,
		int *	ringCounts,
		LWPOINT *	point
	)

Definition at line 931 of file lwgeom_functions_analytic.c.

{
        int i, p, r, in_ring;
        POINT2D pt;
        int result = -1;
 
        POSTGIS_DEBUGF(2, "point_in_multipolygon_rtree called for %p %d %p.", root, polyCount, point);
 
        /* empty is not within anything */
        if (lwpoint_is_empty(point)) return -1;
 
        getPoint2d_p(point->point, 0, &pt);
        /* assume bbox short-circuit has already been attempted */
 
        i = 0; /* the current index into the root array */
 
        /* is the point inside any of the sub-polygons? */
        for ( p = 0; p < polyCount; p++ )
        {
                /* Skip empty polygons */
                if( ringCounts[p] == 0 ) continue;
 
                in_ring = point_in_ring_rtree(root[i], &pt);
                POSTGIS_DEBUGF(4, "point_in_multipolygon_rtree: exterior ring (%d), point_in_ring returned %d", p, in_ring);
                if ( in_ring == -1 ) /* outside the exterior ring */
                {
                        POSTGIS_DEBUG(3, "point_in_multipolygon_rtree: outside exterior ring.");
                }
                else if ( in_ring == 0 ) /* on the boundary */
                {
                        POSTGIS_DEBUGF(3, "point_in_multipolygon_rtree: on edge of exterior ring %d", p);
                        return 0;
                } else {
                        result = in_ring;
 
                        for(r=1; r<ringCounts[p]; r++)
                        {
                                in_ring = point_in_ring_rtree(root[i+r], &pt);
                                POSTGIS_DEBUGF(4, "point_in_multipolygon_rtree: interior ring (%d), point_in_ring returned %d", r, in_ring);
                                if (in_ring == 1) /* inside a hole => outside the polygon */
                                {
                                        POSTGIS_DEBUGF(3, "point_in_multipolygon_rtree: within hole %d of exterior ring %d", r, p);
                                        result = -1;
                                        break;
                                }
                                if (in_ring == 0) /* on the edge of a hole */
                                {
                                        POSTGIS_DEBUGF(3, "point_in_multipolygon_rtree: on edge of hole %d of exterior ring %d", r, p);
                                        return 0;
                                }
                        }
                        /* if we have a positive result, we can short-circuit and return it */
                        if ( result != -1)
                        {
                                return result;
                        }
                }
                /* increment the index by the total number of rings in the sub-poly */
                /* we do this here in case we short-cutted out of the poly before looking at all the rings */
                i += ringCounts[p];
        }
 
        return result; /* -1 = outside, 0 = boundary, 1 = inside */
 
}

References getPoint2d_p(), lwpoint_is_empty(), LWPOINT::point, point_in_ring_rtree(), and r.

Referenced by pip_short_circuit().

Here is the call graph for this function:

Here is the caller graph for this function: