rt_raster_get_envelope()

rt_errorstate rt_raster_get_envelope	(	rt_raster	raster,
		rt_envelope *	env
	)

Get raster's envelope.

The envelope is the minimum bounding rectangle of the raster

Parameters

raster	: the raster to get envelope of
env	: pointer to rt_envelope

Returns

ES_NONE if success, ES_ERROR if error

Definition at line 871 of file rt_raster.c.

  {
        int i;
        int rtn;
        int set = 0;
        double _r[2] = {0.};
        double _w[2] = {0.};
        double _gt[6] = {0.};
 
        assert(raster != NULL);
        assert(env != NULL);
 
        rt_raster_get_geotransform_matrix(raster, _gt);
 
        for (i = 0; i < 4; i++) {
                switch (i) {
                        case 0:
                                _r[0] = 0;
                                _r[1] = 0;
                                break;
                        case 1:
                                _r[0] = 0;
                                _r[1] = raster->height;
                                break;
                        case 2:
                                _r[0] = raster->width;
                                _r[1] = raster->height;
                                break;
                        case 3:
                                _r[0] = raster->width;
                                _r[1] = 0;
                                break;
                }
 
                rtn = rt_raster_cell_to_geopoint(
                        raster,
                        _r[0], _r[1],
                        &(_w[0]), &(_w[1]),
                        _gt
                );
                if (rtn != ES_NONE) {
                        rterror("rt_raster_get_envelope: Could not compute spatial coordinates for raster pixel");
                        return ES_ERROR;
                }
 
                if (!set) {
                        set = 1;
                        env->MinX = _w[0];
                        env->MaxX = _w[0];
                        env->MinY = _w[1];
                        env->MaxY = _w[1];
                }
                else {
                        if (_w[0] < env->MinX)
                                env->MinX = _w[0];
                        else if (_w[0] > env->MaxX)
                                env->MaxX = _w[0];
 
                        if (_w[1] < env->MinY)
                                env->MinY = _w[1];
                        else if (_w[1] > env->MaxY)
                                env->MaxY = _w[1];
                }
        }
 
        return ES_NONE;
}

References ES_ERROR, ES_NONE, rt_envelope::MaxX, rt_envelope::MaxY, rt_envelope::MinX, rt_envelope::MinY, rtrowdump::raster, rt_raster_cell_to_geopoint(), rt_raster_get_geotransform_matrix(), and rterror().

Referenced by rt_raster_get_envelope_geom(), and test_raster_envelope().

Here is the call graph for this function:

Here is the caller graph for this function: