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 2405 of file lwgeodetic.c.

References LWGEOM::bbox, GBOX::flags, LWCOLLECTION::geoms, getPoint2d_p(), LINETYPE, LW_FALSE, LW_TRUE, LWDEBUG, lwerror(), lwgeom_calculate_gbox_geodetic(), lwline_covers_lwline(), lwline_covers_lwpoint(), lwpoint_same(), lwpoly_covers_lwline(), lwpoly_covers_lwpoly(), lwpoly_covers_point2d(), lwtype_is_collection(), 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;

        /* dim(geom2) > dim(geom1) always returns false (because geom2 is bigger) */
        if ( (type1 == POINTTYPE && type2 == LINETYPE)
                || (type1 == POINTTYPE && type2 == POLYGONTYPE)
                || (type1 == LINETYPE && type2 == POLYGONTYPE) )
        {
                LWDEBUG(4, "dimension of geom2 is bigger than geom1");
                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);
        }
        else if ( type1 == POLYGONTYPE && type2 == LINETYPE)
        {
                return lwpoly_covers_lwline((LWPOLY*)lwgeom1, (LWLINE*)lwgeom2);
        }
        else if ( type1 == POLYGONTYPE && type2 == POLYGONTYPE)
        {
                return lwpoly_covers_lwpoly((LWPOLY*)lwgeom1, (LWPOLY*)lwgeom2);
        }
        else if ( type1 == LINETYPE && type2 == POINTTYPE)
        {
                return lwline_covers_lwpoint((LWLINE*)lwgeom1, (LWPOINT*)lwgeom2);
        }
        else if ( type1 == LINETYPE && type2 == LINETYPE)
        {
                return lwline_covers_lwline((LWLINE*)lwgeom1, (LWLINE*)lwgeom2);
        }
        else if ( type1 == POINTTYPE && type2 == POINTTYPE)
        {
                return lwpoint_same((LWPOINT*)lwgeom1, (LWPOINT*)lwgeom2);
        }

        /* 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: