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 637 of file rt_raster.c.

  {
        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.0) || FLT_EQ(_gt[5], 0.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;
}

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().

Here is the call graph for this function:

Here is the caller graph for this function: