Datum geos_intersects

(

PG_FUNCTION_ARGS

)

Definition at line 2571 of file postgis/lwgeom_geos.c.

References PrepGeomCache::argnum, error_if_srid_mismatch(), errorIfGeometryCollection(), FALSE, gbox_overlaps_2d(), GetPrepGeomCache(), GetRtreeCache(), gserialized_get_gbox_p(), gserialized_get_srid(), gserialized_get_type(), gserialized_is_empty(), LW_FALSE, lwerror(), lwgeom_as_lwpoint(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_geos_errmsg, lwgeom_geos_error(), lwnotice(), lwpoint_free(), MULTIPOLYGONTYPE, point_in_multipolygon(), point_in_multipolygon_rtree(), point_in_polygon(), POINTTYPE, RTREE_POLY_CACHE::polyCount, POLYGONTYPE, POSTGIS2GEOS(), PrepGeomCache::prepared_geom, result, RTREE_POLY_CACHE::ringCounts, RTREE_POLY_CACHE::ringIndices, and TRUE.

{
        GSERIALIZED *geom1;
        GSERIALIZED *geom2;
        GSERIALIZED *serialized_poly;
        int result;
        GBOX box1, box2;
        int type1, type2, polytype;
        LWPOINT *point;
        LWGEOM *lwgeom;
        RTREE_POLY_CACHE *poly_cache;
        PrepGeomCache *prep_cache;

        geom1 = (GSERIALIZED *)PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
        geom2 = (GSERIALIZED *)PG_DETOAST_DATUM(PG_GETARG_DATUM(1));

        errorIfGeometryCollection(geom1,geom2);
        error_if_srid_mismatch(gserialized_get_srid(geom1), gserialized_get_srid(geom2));

        /* A.Intersects(Empty) == FALSE */
        if ( gserialized_is_empty(geom1) || gserialized_is_empty(geom2) )
                PG_RETURN_BOOL(false);

        /*
         * short-circuit 1: if geom2 bounding box does not overlap
         * geom1 bounding box we can prematurely return FALSE.
         * Do the test IFF BOUNDING BOX AVAILABLE.
         */
        if ( gserialized_get_gbox_p(geom1, &box1) &&
                gserialized_get_gbox_p(geom2, &box2) )
        {
                if ( gbox_overlaps_2d(&box1, &box2) == LW_FALSE )
                {
                        PG_RETURN_BOOL(FALSE);
                }
        }

        /*
         * short-circuit 2: if the geoms are a point and a polygon,
         * call the point_outside_polygon function.
         */
        type1 = gserialized_get_type(geom1);
        type2 = gserialized_get_type(geom2);
        if ( (type1 == POINTTYPE && (type2 == POLYGONTYPE || type2 == MULTIPOLYGONTYPE)) ||
             (type2 == POINTTYPE && (type1 == POLYGONTYPE || type1 == MULTIPOLYGONTYPE)))
        {
                POSTGIS_DEBUG(3, "Point in Polygon test requested...short-circuiting.");

                if ( type1 == POINTTYPE )
                {
                        point = lwgeom_as_lwpoint(lwgeom_from_gserialized(geom1));
                        lwgeom = lwgeom_from_gserialized(geom2);
                        serialized_poly = geom2;
                        polytype = type2;
                }
                else
                {
                        point = lwgeom_as_lwpoint(lwgeom_from_gserialized(geom2));
                        lwgeom = lwgeom_from_gserialized(geom1);
                        serialized_poly = geom1;
                        polytype = type1;
                }

                poly_cache = GetRtreeCache(fcinfo, serialized_poly);

                if ( poly_cache && poly_cache->ringIndices )
                {
                        result = point_in_multipolygon_rtree(poly_cache->ringIndices, poly_cache->polyCount, poly_cache->ringCounts, point);
                }
                else if ( polytype == POLYGONTYPE )
                {
                        result = point_in_polygon((LWPOLY*)lwgeom, point);
                }
                else if ( polytype == MULTIPOLYGONTYPE )
                {
                        result = point_in_multipolygon((LWMPOLY*)lwgeom, point);
                }
                else
                {
                        /* Gulp! Should not be here... */
                        elog(ERROR,"Type isn't poly or multipoly!");
                        PG_RETURN_NULL();
                }

                lwgeom_free(lwgeom);
                lwpoint_free(point);
                PG_FREE_IF_COPY(geom1, 0);
                PG_FREE_IF_COPY(geom2, 1);
                if ( result != -1 ) /* not outside */
                {
                        PG_RETURN_BOOL(TRUE);
                }
                else
                {
                        PG_RETURN_BOOL(FALSE);
                }
        }

        initGEOS(lwnotice, lwgeom_geos_error);
        prep_cache = GetPrepGeomCache( fcinfo, geom1, geom2 );

        if ( prep_cache && prep_cache->prepared_geom )
        {
                if ( prep_cache->argnum == 1 )
                {
                        GEOSGeometry *g = (GEOSGeometry *)POSTGIS2GEOS(geom2);
                        if ( 0 == g )   /* exception thrown at construction */
                        {
                                lwerror("Geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
                                PG_RETURN_NULL();
                        }
                        result = GEOSPreparedIntersects( prep_cache->prepared_geom, g);
                        GEOSGeom_destroy(g);
                }
                else
                {
                        GEOSGeometry *g = (GEOSGeometry *)POSTGIS2GEOS(geom1);
                        if ( 0 == g )   /* exception thrown at construction */
                        {
                                lwerror("Geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
                                PG_RETURN_NULL();
                        }
                        result = GEOSPreparedIntersects( prep_cache->prepared_geom, g);
                        GEOSGeom_destroy(g);
                }
        }
        else
        {
                GEOSGeometry *g1;
                GEOSGeometry *g2;
                g1 = (GEOSGeometry *)POSTGIS2GEOS(geom1);
                if ( 0 == g1 )   /* exception thrown at construction */
                {
                        lwerror("First argument geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
                        PG_RETURN_NULL();
                }
                g2 = (GEOSGeometry *)POSTGIS2GEOS(geom2);
                if ( 0 == g2 )   /* exception thrown at construction */
                {
                        lwerror("Second argument geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
                        GEOSGeom_destroy(g1);
                        PG_RETURN_NULL();
                }
                result = GEOSIntersects( g1, g2);
                GEOSGeom_destroy(g1);
                GEOSGeom_destroy(g2);
        }

        if (result == 2)
        {
                lwerror("GEOSIntersects: %s", lwgeom_geos_errmsg);
                PG_RETURN_NULL(); /* never get here */
        }

        PG_FREE_IF_COPY(geom1, 0);
        PG_FREE_IF_COPY(geom2, 1);

        PG_RETURN_BOOL(result);
}

Here is the call graph for this function: