rt_raster_geopoint_to_cell()

rt_errorstate rt_raster_geopoint_to_cell	(	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.

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 pixel's column
yr	: output parameter, the pixel's row
igt	: input/output parameter, inverse geotransform matrix

Returns

ES_NONE if success, ES_ERROR if error

Definition at line 806 of file rt_raster.c.

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

Referenced by RASTER_nearestValue(), RASTER_setPixelValuesGeomval(), RASTER_worldToRasterCoord(), rt_band_load_offline_data(), rt_raster_compute_skewed_raster(), rt_raster_from_two_rasters(), rt_raster_gdal_rasterize(), rt_raster_gdal_warp(), rt_raster_intersects(), rt_raster_intersects_algorithm(), rt_raster_same_alignment(), and test_raster_geopoint_to_cell().

  {
        double _igt[6] = {0};
        double rnd = 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);

        rnd = ROUND(*xr, 0);
        if (FLT_EQ(rnd, *xr))
                *xr = rnd;
        else
                *xr = floor(*xr);

        rnd = ROUND(*yr, 0);
        if (FLT_EQ(rnd, *yr))
                *yr = rnd;
        else
                *yr = floor(*yr);

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

        return ES_NONE;
}

Here is the call graph for this function:

Here is the caller graph for this function: