lwgeom_solid_contains_lwgeom()

static int lwgeom_solid_contains_lwgeom ( const LWGEOM *  solid, const LWGEOM *  g  )

inlinestatic

Definition at line 382 of file measures3d.c.

{
        const GBOX *b1;
        const GBOX *b2;
 
        /* If solid isn't solid it can't contain anything */
        if (!FLAGS_GET_SOLID(solid->flags))
                return LW_FALSE;
 
        b1 = lwgeom_get_bbox(solid);
        b2 = lwgeom_get_bbox(g);
 
        /* If box won't contain box, shape won't contain shape */
        if (!gbox_contains_3d(b1, b2))
                return LW_FALSE;
        else /* Raycast upwards in Z */
        {
                POINT4D pt;
                /* We'll skew copies if we're not lucky */
                LWGEOM *solid_copy = lwgeom_clone_deep(solid);
                LWGEOM *g_copy = lwgeom_clone_deep(g);
 
                while (LW_TRUE)
                {
                        uint8_t is_boundary = LW_FALSE;
                        uint8_t is_inside = LW_FALSE;
 
                        /* take first point */
                        if (!lwgeom_startpoint(g_copy, &pt))
                                return LW_FALSE;
 
                        /* get part of solid that is upwards */
                        LWCOLLECTION *c = lwgeom_clip_to_ordinate_range(solid_copy, 'Z', pt.z, DBL_MAX, 0);
 
                        for (uint32_t i = 0; i < c->ngeoms; i++)
                        {
                                if (c->geoms[i]->type == POLYGONTYPE)
                                {
                                        /* 3D raycast along Z is 2D point in polygon */
                                        int t = lwpoly_contains_point((LWPOLY *)c->geoms[i], (POINT2D *)&pt);
 
                                        if (t == LW_INSIDE)
                                                is_inside = !is_inside;
                                        else if (t == LW_BOUNDARY)
                                        {
                                                is_boundary = LW_TRUE;
                                                break;
                                        }
                                }
                                else if (c->geoms[i]->type == TRIANGLETYPE)
                                {
                                        /* 3D raycast along Z is 2D point in polygon */
                                        LWTRIANGLE *tri = (LWTRIANGLE *)c->geoms[i];
                                        int t = ptarray_contains_point(tri->points, (POINT2D *)&pt);
 
                                        if (t == LW_INSIDE)
                                                is_inside = !is_inside;
                                        else if (t == LW_BOUNDARY)
                                        {
                                                is_boundary = LW_TRUE;
                                                break;
                                        }
                                }
                        }
 
                        lwcollection_free(c);
                        lwgeom_free(solid_copy);
                        lwgeom_free(g_copy);
 
                        if (is_boundary)
                        {
                                /* randomly skew a bit and restart*/
                                double cx = lwrandom_uniform() - 0.5;
                                double cy = lwrandom_uniform() - 0.5;
                                AFFINE aff = {1, 0, cx, 0, 1, cy, 0, 0, 1, 0, 0, 0};
 
                                solid_copy = lwgeom_clone_deep(solid);
                                lwgeom_affine(solid_copy, &aff);
 
                                g_copy = lwgeom_clone_deep(g);
                                lwgeom_affine(g_copy, &aff);
 
                                continue;
                        }
                        return is_inside;
                }
        }
}

References LWGEOM::flags, FLAGS_GET_SOLID, gbox_contains_3d(), LWCOLLECTION::geoms, LW_BOUNDARY, LW_FALSE, LW_INSIDE, LW_TRUE, lwcollection_free(), lwgeom_affine(), lwgeom_clip_to_ordinate_range(), lwgeom_clone_deep(), lwgeom_free(), lwgeom_get_bbox(), lwgeom_startpoint(), lwpoly_contains_point(), lwrandom_uniform(), LWCOLLECTION::ngeoms, LWTRIANGLE::points, POLYGONTYPE, ptarray_contains_point(), TRIANGLETYPE, LWGEOM::type, and POINT4D::z.

Referenced by lwgeom_mindistance3d_tolerance().

Here is the call graph for this function:

Here is the caller graph for this function: