n-raster iterator.
840 double offset[4] = {0.};
861 assert(itrset != NULL && itrcount > 0);
862 assert(rtnraster != NULL);
868 if (callback == NULL) {
869 rterror(
"rt_raster_iterator: Callback function not provided");
875 rterror(
"rt_raster_iterator: Custom extent cannot be empty if extent type is ET_CUSTOM");
881 rterror(
"rt_raster_iterator: Pixel type cannot be PT_END");
887 rterror(
"rt_raster_iterator: Could not initialize internal variables");
893 rterror(
"rt_raster_iterator: Could not populate for internal variables");
899 if (allnull == itrcount) {
900 RASTER_DEBUG(3,
"all rasters are NULL, returning NULL");
907 else if (allempty == itrcount) {
908 RASTER_DEBUG(3,
"all rasters are empty, returning empty raster");
913 if (rtnrast == NULL) {
914 rterror(
"rt_raster_iterator: Could not create empty raster");
919 *rtnraster = rtnrast;
930 RASTER_DEBUG(4,
"using custom extent as reference raster");
935 for (i = 0; i < itrcount; i++) {
937 RASTER_DEBUGF(4,
"using raster at index %d as reference raster", i);
946 rterror(
"rt_raster_iterator: Could not find reference raster to use for alignment tests");
959 rterror(
"rt_raster_iterator: Could not test for alignment between reference raster and custom extent");
971 for (i = 0; i < itrcount; i++) {
977 rterror(
"rt_raster_iterator: Could not test for alignment between reference raster and raster %d", i);
983 RASTER_DEBUGF(5,
"raster at index %d alignment: %d", i, aligned);
994 rterror(
"rt_raster_iterator: The set of rasters provided (custom extent included, if appropriate) do not have the same alignment");
1003 switch (extenttype) {
1008 if (rtnrast == NULL) {
1009 rterror(
"rt_raster_iterator: Could not allocate memory for output raster");
1016 for (i = 0; i < itrcount; i++) {
1023 rtnrast->
bands = NULL;
1027 for (i = i + 1; i < itrcount; i++) {
1035 rterror(
"rt_raster_iterator: Could not compute %s extent of rasters",
1036 extenttype ==
ET_UNION ?
"union" :
"intersection"
1044 rtinfo(
"rt_raster_iterator: Computed raster for %s extent is empty",
1045 extenttype ==
ET_UNION ?
"union" :
"intersection"
1075 if (i < 0) i = itrcount - 1;
1078 if (_param->
raster[i] == NULL) {
1079 RASTER_DEBUGF(3,
"returning NULL as %s raster is NULL and extent type is ET_%s",
1080 (i == 0 ?
"first" : (i == 1 ?
"second" :
"last")),
1081 (i == 0 ?
"FIRST" : (i == 1 ?
"SECOND" :
"LAST"))
1089 else if (_param->
isempty[i]) {
1090 RASTER_DEBUGF(3,
"returning empty raster as %s raster is empty and extent type is ET_%s",
1091 (i == 0 ?
"first" : (i == 1 ?
"second" :
"last")),
1092 (i == 0 ?
"FIRST" : (i == 1 ?
"SECOND" :
"LAST"))
1098 if (rtnrast == NULL) {
1099 rterror(
"rt_raster_iterator: Could not create empty raster");
1104 *rtnraster = rtnrast;
1111 if (rtnrast == NULL) {
1112 rterror(
"rt_raster_iterator: Could not allocate memory for output raster");
1119 switch (extenttype) {
1121 memcpy(rtnrast, customextent,
sizeof(
struct rt_raster_t));
1129 rtnrast->
bands = NULL;
1136 RASTER_DEBUGF(4,
"rtnrast (width, height, ulx, uly, scalex, scaley, skewx, skewy, srid) = (%d, %d, %f, %f, %f, %f, %f, %f, %d)",
1150 rterror(
"rt_raster_iterator: Could not initialize empty values and NODATA");
1163 hasnodata, nodataval,
1166 rterror(
"rt_raster_iterator: Could not add new band to output raster");
1176 if (rtnband == NULL) {
1177 rterror(
"rt_raster_iterator: Could not get new band from output raster");
1190 rterror(
"rt_raster_iterator: Could not initialize callback function argument");
1200 for (i = 0; i < itrcount; i++) {
1207 rterror(
"rt_raster_iterator: Could not compute raster offsets");
1216 _param->
offset[i][0] = offset[2];
1217 _param->
offset[i][1] = offset[3];
1218 RASTER_DEBUGF(4,
"rast %d offset: %f %f", i, offset[2], offset[3]);
1224 for (_y = 0; _y < _height; _y++) {
1225 for (_x = 0; _x < _width; _x++) {
1226 RASTER_DEBUGF(4,
"iterating output pixel (x, y) = (%d, %d)", _x, _y);
1231 for (i = 0; i < itrcount; i++) {
1241 (_param->
band.
rtband[i] == NULL && itrset[i].nbnodata) ||
1244 RASTER_DEBUG(4,
"empty raster, band does not exist or band is NODATA. using empty values and NODATA");
1256 x = _x - (int) _param->
offset[i][0];
1257 y = _y - (
int) _param->
offset[i][1];
1266 if (distancex > 0 && distancey > 0) {
1272 distancex, distancey,
1277 rterror(
"rt_raster_iterator: Could not get pixel neighborhood");
1290 (
x >= 0 && x < _param->width[i]) &&
1291 (
y >= 0 && y < _param->height[i])
1300 rterror(
"rt_raster_iterator: Could not get the pixel value of band");
1312 RASTER_DEBUG(4,
"Outside band extent, setting value to NODATA");
1331 if (npixels == NULL) {
1332 rterror(
"rt_raster_iterator: Could not reallocate memory for neighborhood");
1341 npixels[status - 1].
x =
x;
1342 npixels[status - 1].
y =
y;
1343 npixels[status - 1].
nodata = 1;
1347 if ((!_param->
band.
hasnodata[i] && inextent) || !isnodata) {
1348 npixels[status - 1].
nodata = 0;
1354 npixels, status, mask,
1356 distancex, distancey,
1363 rterror(
"rt_raster_iterator: Could not create 2D array of neighborhood");
1377 status = callback(_param->
arg, userarg, &
value, &nodata);
1384 rterror(
"rt_raster_iterator: Callback function returned an error");
1399 else if (!hasnodata) {
1401 RASTER_DEBUGF(4,
"burning pixel (%d, %d) with minval: %f", _x, _y, minval);
1407 rterror(
"rt_raster_iterator: Could not set pixel value");
1421 *rtnraster = rtnrast;
void rterror(const char *fmt,...)
Wrappers used for reporting errors and info.
void * rtalloc(size_t size)
Wrappers used for managing memory.
#define RASTER_DEBUG(level, msg)
int32_t rt_raster_get_srid(rt_raster raster)
Get raster's SRID.
#define RASTER_DEBUGF(level, msg,...)
double rt_raster_get_x_skew(rt_raster raster)
Get skew about the X axis.
double rt_raster_get_x_offset(rt_raster raster)
Get raster x offset, in projection units.
int rt_raster_generate_new_band(rt_raster raster, rt_pixtype pixtype, double initialvalue, uint32_t hasnodata, double nodatavalue, int index)
Generate a new inline band and add it to a raster.
void rtinfo(const char *fmt,...)
void rt_raster_set_scale(rt_raster raster, double scaleX, double scaleY)
Set scale in projection units.
rt_errorstate rt_band_get_pixel(rt_band band, int x, int y, double *value, int *nodata)
Get pixel value.
void rt_raster_destroy(rt_raster raster)
Release memory associated to a raster.
rt_raster rt_raster_new(uint32_t width, uint32_t height)
Construct a raster with given dimensions.
struct rt_pixel_t * rt_pixel
double rt_raster_get_x_scale(rt_raster raster)
Get scale X in projection units.
double rt_band_get_min_value(rt_band band)
Returns the minimal possible value for the band according to the pixel type.
rt_errorstate rt_band_set_pixel(rt_band band, int x, int y, double val, int *converted)
Set single pixel's value.
void rt_band_destroy(rt_band band)
Destroy a raster band.
void * rtrealloc(void *mem, size_t size)
uint16_t rt_raster_get_height(rt_raster raster)
rt_errorstate rt_pixel_set_to_array(rt_pixel npixel, uint32_t count, rt_mask mask, int x, int y, uint16_t distancex, uint16_t distancey, double ***value, int ***nodata, int *dimx, int *dimy)
uint16_t rt_raster_get_width(rt_raster raster)
uint32_t rt_band_get_nearest_pixel(rt_band band, int x, int y, uint16_t distancex, uint16_t distancey, int exclude_nodata_value, rt_pixel *npixels)
Get nearest pixel(s) with value (not NODATA) to specified pixel.
void rtdealloc(void *mem)
rt_errorstate rt_raster_from_two_rasters(rt_raster rast1, rt_raster rast2, rt_extenttype extenttype, rt_raster *rtnraster, double *offset)
double rt_raster_get_y_scale(rt_raster raster)
Get scale Y in projection units.
rt_errorstate rt_raster_same_alignment(rt_raster rast1, rt_raster rast2, int *aligned, char **reason)
double rt_raster_get_y_skew(rt_raster raster)
Get skew about the Y axis.
double rt_raster_get_y_offset(rt_raster raster)
Get raster y offset, in projection units.
int rt_raster_is_empty(rt_raster raster)
Return TRUE if the raster is empty.
rt_band rt_raster_get_band(rt_raster raster, int bandNum)
Return Nth band, or NULL if unavailable.
static int _rti_iterator_arg_callback_init(_rti_iterator_arg _param)
static int _rti_iterator_arg_empty_init(_rti_iterator_arg _param)
static _rti_iterator_arg _rti_iterator_arg_init()
static int _rti_iterator_arg_populate(_rti_iterator_arg _param, rt_iterator itrset, uint16_t itrcount, uint16_t distancex, uint16_t distancey, int *allnull, int *allempty)
static void _rti_iterator_arg_callback_clean(_rti_iterator_arg _param)
static void _rti_iterator_arg_destroy(_rti_iterator_arg _param)
struct _rti_iterator_arg_t::@10 band
struct _rti_iterator_arg_t::@13 empty