lwgeom_covers_lwgeom_sphere()

int lwgeom_covers_lwgeom_sphere	(	const LWGEOM *	lwgeom1,
		const LWGEOM *	lwgeom2
	)

Calculate covers predicate for two lwgeoms on the sphere.

Currently only handles point-in-polygon.

Definition at line 2314 of file lwgeodetic.c.

References LWGEOM::bbox, COLLECTIONTYPE, GBOX::flags, LWCOLLECTION::geoms, getPoint2d_p(), LW_FALSE, LW_TRUE, lwerror(), lwgeom_calculate_gbox_geodetic(), lwpoly_covers_point2d(), lwtype_is_collection(), MULTIPOINTTYPE, MULTIPOLYGONTYPE, LWCOLLECTION::ngeoms, POINTTYPE, POLYGONTYPE, and LWGEOM::type.

Referenced by geography_covers().

{
        int type1, type2;
        GBOX gbox1, gbox2;
        gbox1.flags = gbox2.flags = 0;

        assert(lwgeom1);
        assert(lwgeom2);

        type1 = lwgeom1->type;
        type2 = lwgeom2->type;

        /* Currently a restricted implementation */
        if ( ! ( (type1 == POLYGONTYPE || type1 == MULTIPOLYGONTYPE || type1 == COLLECTIONTYPE) &&
                 (type2 == POINTTYPE || type2 == MULTIPOINTTYPE || type2 == COLLECTIONTYPE) ) )
        {
                lwerror("lwgeom_covers_lwgeom_sphere: only POLYGON covers POINT tests are currently supported");
                return LW_FALSE;
        }

        /* Make sure we have boxes */
        if ( lwgeom1->bbox )
                gbox1 = *(lwgeom1->bbox);
        else
                lwgeom_calculate_gbox_geodetic(lwgeom1, &gbox1);

        /* Make sure we have boxes */
        if ( lwgeom2->bbox )
                gbox2 = *(lwgeom2->bbox);
        else
                lwgeom_calculate_gbox_geodetic(lwgeom2, &gbox2);


        /* Handle the polygon/point case */
        if ( type1 == POLYGONTYPE && type2 == POINTTYPE )
        {
                POINT2D pt_to_test;
                getPoint2d_p(((LWPOINT*)lwgeom2)->point, 0, &pt_to_test);
                return lwpoly_covers_point2d((LWPOLY*)lwgeom1, &pt_to_test);
        }

        /* If any of the first argument parts covers the second argument, it's true */
        if ( lwtype_is_collection( type1 ) )
        {
                int i;
                LWCOLLECTION *col = (LWCOLLECTION*)lwgeom1;

                for ( i = 0; i < col->ngeoms; i++ )
                {
                        if ( lwgeom_covers_lwgeom_sphere(col->geoms[i], lwgeom2) )
                        {
                                return LW_TRUE;
                        }
                }
                return LW_FALSE;
        }

        /* Only if all of the second arguments are covered by the first argument is the condition true */
        if ( lwtype_is_collection( type2 ) )
        {
                int i;
                LWCOLLECTION *col = (LWCOLLECTION*)lwgeom2;

                for ( i = 0; i < col->ngeoms; i++ )
                {
                        if ( ! lwgeom_covers_lwgeom_sphere(lwgeom1, col->geoms[i]) )
                        {
                                return LW_FALSE;
                        }
                }
                return LW_TRUE;
        }

        /* Don't get here */
        lwerror("lwgeom_covers_lwgeom_sphere: reached end of function without resolution");
        return LW_FALSE;

}

Here is the call graph for this function:

Here is the caller graph for this function: