rt_raster_geos_spatial_relationship()

static rt_errorstate rt_raster_geos_spatial_relationship ( rt_raster  rast1, int  nband1, rt_raster  rast2, int  nband2, rt_geos_spatial_test  testtype, int *  testresult  )

static

Definition at line 137 of file rt_spatial_relationship.c.

References ES_ERROR, ES_NONE, GSR_CONTAINS, GSR_CONTAINSPROPERLY, GSR_COVEREDBY, GSR_COVERS, GSR_OVERLAPS, GSR_TOUCHES, LWGEOM2GEOS(), lwgeom_geos_error(), lwmpoly_as_lwgeom(), lwmpoly_free(), RASTER_DEBUG, rt_raster_get_num_bands(), rt_raster_get_srid(), rt_raster_surface(), rterror(), and rtinfo().

Referenced by rt_raster_contains(), rt_raster_contains_properly(), rt_raster_coveredby(), rt_raster_covers(), rt_raster_overlaps(), and rt_raster_touches().

  {
        LWMPOLY *surface1 = NULL;
        LWMPOLY *surface2 = NULL;
        GEOSGeometry *geom1 = NULL;
        GEOSGeometry *geom2 = NULL;
        int rtn = 0;
        int flag = 0;

        RASTER_DEBUG(3, "Starting");

        assert(NULL != rast1);
        assert(NULL != rast2);
        assert(NULL != testresult);

        if (nband1 < 0 && nband2 < 0) {
                nband1 = -1;
                nband2 = -1;
        }
        else {
                assert(nband1 >= 0 && nband1 < rt_raster_get_num_bands(rast1));
                assert(nband2 >= 0 && nband2 < rt_raster_get_num_bands(rast2));
        }

        /* initialize to zero, false result of spatial relationship test */
        *testresult = 0;

        /* same srid */
        if (rt_raster_get_srid(rast1) != rt_raster_get_srid(rast2)) {
                rterror("rt_raster_geos_spatial_relationship: The two rasters provided have different SRIDs");
                return ES_ERROR;
        }

        initGEOS(rtinfo, lwgeom_geos_error);

        /* get LWMPOLY of each band */
        if (rt_raster_surface(rast1, nband1, &surface1) != ES_NONE) {
                rterror("rt_raster_geos_spatial_relationship: Could not get surface of the specified band from the first raster");
                return ES_ERROR;
        }
        if (rt_raster_surface(rast2, nband2, &surface2) != ES_NONE) {
                rterror("rt_raster_geos_spatial_relationship: Could not get surface of the specified band from the second raster");
                lwmpoly_free(surface1);
                return ES_ERROR;
        }

        /* either surface is NULL, spatial relationship test is false */
        if (surface1 == NULL || surface2 == NULL) {
                if (surface1 != NULL) lwmpoly_free(surface1);
                if (surface2 != NULL) lwmpoly_free(surface2);
                return ES_NONE;
        }

        /* convert LWMPOLY to GEOSGeometry */
        geom1 = LWGEOM2GEOS(lwmpoly_as_lwgeom(surface1), 0);
        lwmpoly_free(surface1);
        if (geom1 == NULL) {
                rterror("rt_raster_geos_spatial_relationship: Could not convert surface of the specified band from the first raster to a GEOSGeometry");
                lwmpoly_free(surface2);
                return ES_ERROR;
        }

        geom2 = LWGEOM2GEOS(lwmpoly_as_lwgeom(surface2), 0);
        lwmpoly_free(surface2);
        if (geom2 == NULL) {
                rterror("rt_raster_geos_spatial_relationship: Could not convert surface of the specified band from the second raster to a GEOSGeometry");
                return ES_ERROR;
        }

        flag = 0;
        switch (testtype) {
                case GSR_OVERLAPS:
                        rtn = GEOSOverlaps(geom1, geom2);
                        break;
                case GSR_TOUCHES:
                        rtn = GEOSTouches(geom1, geom2);
                        break;
                case GSR_CONTAINS:
                        rtn = GEOSContains(geom1, geom2);
                        break;
                case GSR_CONTAINSPROPERLY:
                        rtn = GEOSRelatePattern(geom1, geom2, "T**FF*FF*");
                        break;
                case GSR_COVERS:
                        rtn = GEOSRelatePattern(geom1, geom2, "******FF*");
                        break;
                case GSR_COVEREDBY:
                        rtn = GEOSRelatePattern(geom1, geom2, "**F**F***");
                        break;
                default:
                        rterror("rt_raster_geos_spatial_relationship: Unknown or unsupported GEOS spatial relationship test");
                        flag = -1;
                        break;
        }
        GEOSGeom_destroy(geom1);
        GEOSGeom_destroy(geom2);

        /* something happened in the spatial relationship test */
        if (rtn == 2) {
                rterror("rt_raster_geos_spatial_relationship: Could not run the appropriate GEOS spatial relationship test");
                flag = ES_ERROR;
        }
        /* spatial relationship test ran fine */
        else if (flag >= 0) {
                if (rtn != 0)
                        *testresult = 1;
                flag = ES_NONE;
        }
        /* flag < 0 for when testtype is unknown */
        else
                flag = ES_ERROR;

        return flag;
}

Here is the call graph for this function:

Here is the caller graph for this function: