rt_raster_cell_to_geopoint()

rt_errorstate rt_raster_cell_to_geopoint	(	rt_raster	raster,
		double	xr,
		double	yr,
		double *	xw,
		double *	yw,
		double *	gt
	)

Convert an xr, yr raster point to an xw, yw point on map.

Parameters

raster	: the raster to get info from
xr	: the pixel's column
yr	: the pixel's row
xw	: output parameter, X ordinate of the geographical point
yw	: output parameter, Y ordinate of the geographical point
gt	: input/output parameter, 3x2 geotransform matrix

Returns

ES_NONE if success, ES_ERROR if error

Definition at line 755 of file rt_raster.c.

References ES_NONE, FLT_EQ, RASTER_DEBUGF, and rt_raster_get_geotransform_matrix().

Referenced by RASTER_rasterToWorldCoord(), RASTER_tile(), rt_raster_compute_skewed_raster(), rt_raster_from_two_rasters(), rt_raster_gdal_rasterize(), rt_raster_gdal_warp(), rt_raster_get_convex_hull(), rt_raster_get_envelope(), rt_raster_get_envelope_geom(), rt_raster_get_perimeter(), rt_raster_intersects(), rt_raster_intersects_algorithm(), rt_raster_same_alignment(), and test_raster_cell_to_geopoint().

  {
        double _gt[6] = {0};

        assert(NULL != raster);
        assert(NULL != xw && NULL != yw);

        if (NULL != gt)
                memcpy(_gt, gt, sizeof(double) * 6);

        /* scale of matrix is not set */
        if (
                FLT_EQ(_gt[1], 0) ||
                FLT_EQ(_gt[5], 0)
        ) {
                rt_raster_get_geotransform_matrix(raster, _gt);
        }

        RASTER_DEBUGF(4, "gt = (%f, %f, %f, %f, %f, %f)",
                _gt[0],
                _gt[1],
                _gt[2],
                _gt[3],
                _gt[4],
                _gt[5]
        );

        GDALApplyGeoTransform(_gt, xr, yr, xw, yw);
        RASTER_DEBUGF(4, "GDALApplyGeoTransform (c -> g) for (%f, %f) = (%f, %f)",
                xr, yr, *xw, *yw);

        return ES_NONE;
}

Here is the call graph for this function:

Here is the caller graph for this function: