relate_pattern()

Datum relate_pattern

(

PG_FUNCTION_ARGS

)

Definition at line 900 of file lwgeom_geos_predicates.c.

{
        SHARED_GSERIALIZED *shared_geom1 = ToastCacheGetGeometry(fcinfo, 0);
        SHARED_GSERIALIZED *shared_geom2 = ToastCacheGetGeometry(fcinfo, 1);
        const GSERIALIZED *geom1 = shared_gserialized_get(shared_geom1);
        const GSERIALIZED *geom2 = shared_gserialized_get(shared_geom2);
 
        /* Ensure DE9IM pattern is no more than 9 chars */
        text *imPtr = DatumGetTextP(DirectFunctionCall2(text_left,
                PG_GETARG_DATUM(2), Int32GetDatum(9)));
        char *im = text_to_cstring(imPtr);
        char result;
        uint32_t i;
#if POSTGIS_GEOS_VERSION >= 31300
        PrepGeomCache *prep_cache;
#endif
 
        gserialized_error_if_srid_mismatch(geom1, geom2, __func__);
 
        /*
        ** Need to make sure 't' and 'f' are upper-case before handing to GEOS
        */
        for (i = 0; i < strlen(im); i++)
                im[i] = toupper(im[i]);
 
        initGEOS(lwpgnotice, lwgeom_geos_error);
 
#if POSTGIS_GEOS_VERSION >= 31300
        prep_cache = GetPrepGeomCache(fcinfo, shared_geom1, shared_geom2);
        if ( prep_cache && prep_cache->prepared_geom )
        {
                GEOSGeometry *g = NULL;
                if (prep_cache->gcache.argnum == 1)
                {
                        g = POSTGIS2GEOS(geom2);
                }
                else
                {
                        g = POSTGIS2GEOS(geom1);
                        imInvert(im);
                }
 
                if (!g) HANDLE_GEOS_ERROR("Geometry could not be converted to GEOS");
                result = GEOSPreparedRelatePattern(prep_cache->prepared_geom, g, im);
                GEOSGeom_destroy(g);
        }
        else
#endif
        {
                GEOSGeometry *g1, *g2;
                g1 = POSTGIS2GEOS(geom1);
                if (!g1) HANDLE_GEOS_ERROR("First argument geometry could not be converted to GEOS");
                g2 = POSTGIS2GEOS(geom2);
                if (!g2)
                {
                        HANDLE_GEOS_ERROR("Second argument geometry could not be converted to GEOS");
                        GEOSGeom_destroy(g1);
                }
                result = GEOSRelatePattern(g1, g2, im);
                GEOSGeom_destroy(g1);
                GEOSGeom_destroy(g2);
        }
 
        pfree(im);
        if (result == 2) HANDLE_GEOS_ERROR("GEOSRelatePattern");
 
        PG_RETURN_BOOL(result);
}

References PrepGeomCache::gcache, GetPrepGeomCache(), gserialized_error_if_srid_mismatch(), HANDLE_GEOS_ERROR, lwgeom_geos_error(), POSTGIS2GEOS(), PrepGeomCache::prepared_geom, and result.

Here is the call graph for this function: