rt_raster_same_alignment()

rt_errorstate rt_raster_same_alignment	(	rt_raster	rast1,
		rt_raster	rast2,
		int *	aligned,
		char **	reason
	)

Definition at line 45 of file rt_spatial_relationship.c.

References ES_ERROR, ES_NONE, FLT_EQ, FLT_NEQ, rt_raster_t::ipX, rt_raster_t::ipY, RASTER_DEBUGF, rt_raster_cell_to_geopoint(), rt_raster_geopoint_to_cell(), rt_raster_get_srid(), rterror(), rt_raster_t::scaleX, rt_raster_t::scaleY, rt_raster_t::skewX, and rt_raster_t::skewY.

Referenced by diff_rastinfo(), RASTER_addBandOutDB(), RASTER_mapAlgebra2(), RASTER_notSameAlignmentReason(), RASTER_sameAlignment(), rt_band_load_offline_data(), rt_raster_from_two_rasters(), rt_raster_iterator(), and test_raster_same_alignment().

  {
        double xr;
        double yr;
        double xw;
        double yw;
        int err = 0;

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

        err = 0;
        /* same srid */
        if (rt_raster_get_srid(rast1) != rt_raster_get_srid(rast2)) {
                if (reason != NULL) *reason = "The rasters have different SRIDs";
                err = 1;
        }
        /* scales must match */
        else if (FLT_NEQ(fabs(rast1->scaleX), fabs(rast2->scaleX))) {
                if (reason != NULL) *reason = "The rasters have different scales on the X axis";
                err = 1;
        }
        else if (FLT_NEQ(fabs(rast1->scaleY), fabs(rast2->scaleY))) {
                if (reason != NULL) *reason = "The rasters have different scales on the Y axis";
                err = 1;
        }
        /* skews must match */
        else if (FLT_NEQ(rast1->skewX, rast2->skewX)) {
                if (reason != NULL) *reason = "The rasters have different skews on the X axis";
                err = 1;
        }
        else if (FLT_NEQ(rast1->skewY, rast2->skewY)) {
                if (reason != NULL) *reason = "The rasters have different skews on the Y axis";
                err = 1;
        }

        if (err) {
                *aligned = 0;
                return ES_NONE;
        }

        /* raster coordinates in context of second raster of first raster's upper-left corner */
        if (rt_raster_geopoint_to_cell(
                        rast2,
                        rast1->ipX, rast1->ipY,
                        &xr, &yr,
                        NULL
        ) != ES_NONE) {
                rterror("rt_raster_same_alignment: Could not get raster coordinates of second raster from first raster's spatial coordinates");
                *aligned = 0;
                return ES_ERROR;
        }

        /* spatial coordinates of raster coordinates from above */
        if (rt_raster_cell_to_geopoint(
                rast2,
                xr, yr,
                &xw, &yw,
                NULL
        ) != ES_NONE) {
                rterror("rt_raster_same_alignment: Could not get spatial coordinates of second raster from raster coordinates");
                *aligned = 0;
                return ES_ERROR;
        }

        RASTER_DEBUGF(4, "rast1(ipX, ipxY) = (%f, %f)", rast1->ipX, rast1->ipY);
        RASTER_DEBUGF(4, "rast2(xr, yr) = (%f, %f)", xr, yr);
        RASTER_DEBUGF(4, "rast2(xw, yw) = (%f, %f)", xw, yw);

        /* spatial coordinates are identical to that of first raster's upper-left corner */
        if (FLT_EQ(xw, rast1->ipX) && FLT_EQ(yw, rast1->ipY)) {
                if (reason != NULL) *reason = "The rasters are aligned";
                *aligned = 1;
                return ES_NONE;
        }

        /* no alignment */
        if (reason != NULL) *reason = "The rasters (pixel corner coordinates) are not aligned";

        *aligned = 0;
        return ES_NONE;
}

Here is the call graph for this function:

Here is the caller graph for this function: