rt_raster_geopoint_to_rasterpoint()

rt_errorstate rt_raster_geopoint_to_rasterpoint	(	rt_raster	raster,
		double	xw,
		double	yw,
		double *	xr,
		double *	yr,
		double *	igt
	)

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

Parameters

raster	: the raster to get info from
xw	: X ordinate of the geographical point
yw	: Y ordinate of the geographical point
xr	: output parameter, the x ordinate in raster space
yr	: output parameter, the y ordinate in raster space
igt	: input/output parameter, inverse geotransform matrix

Returns

ES_NONE if success, ES_ERROR if error

Definition at line 730 of file rt_raster.c.

  {
        double _igt[6] = {0};
 
        assert(NULL != raster);
        assert(NULL != xr && NULL != yr);
 
        if (igt != NULL)
                memcpy(_igt, igt, sizeof(double) * 6);
 
        /* matrix is not set */
        if (
                FLT_EQ(_igt[0], 0.) &&
                FLT_EQ(_igt[1], 0.) &&
                FLT_EQ(_igt[2], 0.) &&
                FLT_EQ(_igt[3], 0.) &&
                FLT_EQ(_igt[4], 0.) &&
                FLT_EQ(_igt[5], 0.)
        ) {
                if (rt_raster_get_inverse_geotransform_matrix(raster, NULL, _igt) != ES_NONE) {
                        rterror("rt_raster_geopoint_to_cell: Could not get inverse geotransform matrix");
                        return ES_ERROR;
                }
        }
 
        GDALApplyGeoTransform(_igt, xw, yw, xr, yr);
        RASTER_DEBUGF(4, "GDALApplyGeoTransform (g -> c) for (%f, %f) = (%f, %f)",
                xw, yw, *xr, *yr);
 
        return ES_NONE;
}

References ES_ERROR, ES_NONE, FLT_EQ, RASTER_DEBUGF, rt_raster_get_inverse_geotransform_matrix(), and rterror().

Referenced by RASTER_getPixelValueResample(), rt_raster_copy_to_geometry(), rt_raster_geopoint_to_cell(), and test_raster_get_pixel_bilinear().

Here is the call graph for this function:

Here is the caller graph for this function: