point_in_ring_rtree()

static int point_in_ring_rtree ( RTREE_NODE *  root, const POINT2D *  point  )

static

Definition at line 722 of file lwgeom_functions_analytic.c.

References determineSide(), FP_CONTAINS_BOTTOM, LWMLINE::geoms, getPoint2d_cp(), isOnSegment(), LWMLINE::ngeoms, LWLINE::points, RTreeFindLineSegments(), POINT2D::x, and POINT2D::y.

Referenced by point_in_multipolygon_rtree(), and point_in_polygon_rtree().

{
        int wn = 0;
        int i;
        double side;
        const POINT2D *seg1;
        const POINT2D *seg2;
        LWMLINE *lines;

        POSTGIS_DEBUG(2, "point_in_ring called.");

        lines = RTreeFindLineSegments(root, point->y);
        if (!lines)
                return -1;

        for (i=0; i<lines->ngeoms; i++)
        {
                seg1 = getPoint2d_cp(lines->geoms[i]->points, 0);
                seg2 = getPoint2d_cp(lines->geoms[i]->points, 1);

                side = determineSide(seg1, seg2, point);

                POSTGIS_DEBUGF(3, "segment: (%.8f, %.8f),(%.8f, %.8f)", seg1->x, seg1->y, seg2->x, seg2->y);
                POSTGIS_DEBUGF(3, "side result: %.8f", side);
                POSTGIS_DEBUGF(3, "counterclockwise wrap %d, clockwise wrap %d", FP_CONTAINS_BOTTOM(seg1->y, point->y, seg2->y), FP_CONTAINS_BOTTOM(seg2->y, point->y, seg1->y));

                /* zero length segments are ignored. */
                if (((seg2->x - seg1->x)*(seg2->x - seg1->x) + (seg2->y - seg1->y)*(seg2->y - seg1->y)) < 1e-12*1e-12)
                {
                        POSTGIS_DEBUG(3, "segment is zero length... ignoring.");

                        continue;
                }

                /* a point on the boundary of a ring is not contained. */
                /* WAS: if (fabs(side) < 1e-12), see #852 */
                if (side == 0.0)
                {
                        if (isOnSegment(seg1, seg2, point) == 1)
                        {
                                POSTGIS_DEBUGF(3, "point on ring boundary between points %d, %d", i, i+1);

                                return 0;
                        }
                }

                /*
                 * If the point is to the left of the line, and it's rising,
                 * then the line is to the right of the point and
                 * circling counter-clockwise, so incremement.
                 */
                if ((seg1->y <= point->y) && (point->y < seg2->y) && (side > 0))
                {
                        POSTGIS_DEBUG(3, "incrementing winding number.");

                        ++wn;
                }
                /*
                 * If the point is to the right of the line, and it's falling,
                 * then the line is to the right of the point and circling
                 * clockwise, so decrement.
                 */
                else if ((seg2->y <= point->y) && (point->y < seg1->y) && (side < 0))
                {
                        POSTGIS_DEBUG(3, "decrementing winding number.");

                        --wn;
                }
        }

        POSTGIS_DEBUGF(3, "winding number %d", wn);

        if (wn == 0)
                return -1;
        return 1;
}

Here is the call graph for this function:

Here is the caller graph for this function: