void Rast3d_change_precision(void *map, int precision, const char *nameOut)
{
void *map2;
int x, y, z, savePrecision, saveCompression;
char *data;
RASTER3D_Region region;
int typeIntern;
int nx, ny, nz;
int tileXsave, tileYsave, tileZsave, tileX, tileY, tileZ, saveType;
saveType = Rast3d_get_file_type();
/* Rast3d_set_file_type (Rast3d_file_type_map (map)); */
Rast3d_get_compression_mode(&saveCompression, &savePrecision);
Rast3d_set_compression_mode(RASTER3D_COMPRESSION, precision);
Rast3d_get_tile_dimension(&tileXsave, &tileYsave, &tileZsave);
Rast3d_get_tile_dimensions_map(map, &tileX, &tileY, &tileZ);
Rast3d_set_tile_dimension(tileX, tileY, tileZ);
typeIntern = Rast3d_tile_type_map(map);
Rast3d_get_region_struct_map(map, ®ion);
map2 =
Rast3d_open_cell_new(nameOut, typeIntern, RASTER3D_USE_CACHE_DEFAULT, ®ion);
if (map2 == NULL)
Rast3d_fatal_error("Rast3d_change_precision: error in Rast3d_open_cell_new");
Rast3d_set_file_type(saveType);
Rast3d_set_compression_mode(saveCompression, savePrecision);
Rast3d_set_tile_dimension(tileXsave, tileYsave, tileZsave);
data = Rast3d_alloc_tiles(map, 1);
if (data == NULL)
Rast3d_fatal_error("Rast3d_change_precision: error in Rast3d_alloc_tiles");
Rast3d_get_nof_tiles_map(map2, &nx, &ny, &nz);
for (z = 0; z < nz; z++)
for (y = 0; y < ny; y++)
for (x = 0; x < nx; x++) {
if (!Rast3d_read_tile(map, Rast3d_tile2tile_index(map, x, y, z), data,
typeIntern))
Rast3d_fatal_error
("Rast3d_change_precision: error in Rast3d_read_tile");
if (!Rast3d_write_tile
(map2, Rast3d_tile2tile_index(map2, x, y, z), data,
typeIntern))
Rast3d_fatal_error
("Rast3d_change_precision: error in Rast3d_write_tile");
}
Rast3d_free_tiles(data);
if (!Rast3d_close(map2))
Rast3d_fatal_error("Rast3d_change_precision: error in Rast3d_close");
}
static void
retileNocache(void *map, const char *nameOut, int tileX, int tileY, int tileZ)
{
void *map2;
int x, y, z, saveType, nx, ny, nz;
int typeIntern;
void *data;
int tileXsave, tileYsave, tileZsave;
RASTER3D_Region region;
saveType = Rast3d_get_file_type();
Rast3d_set_file_type(Rast3d_file_type_map(map));
Rast3d_get_tile_dimension(&tileXsave, &tileYsave, &tileZsave);
Rast3d_set_tile_dimension(tileX, tileY, tileZ);
typeIntern = Rast3d_tile_type_map(map);
Rast3d_get_region_struct_map(map, ®ion);
map2 = Rast3d_open_cell_new(nameOut, typeIntern, RASTER3D_NO_CACHE, ®ion);
if (map2 == NULL)
Rast3d_fatal_error("Rast3d_retile: error in Rast3d_open_cell_new");
Rast3d_set_file_type(saveType);
Rast3d_set_tile_dimension(tileXsave, tileYsave, tileZsave);
data = Rast3d_alloc_tiles(map2, 1);
if (data == NULL)
Rast3d_fatal_error("Rast3d_retile: error in Rast3d_alloc_tiles");
Rast3d_get_nof_tiles_map(map2, &nx, &ny, &nz);
for (z = 0; z < nz; z++) {
G_percent(z, nz, 1);
for (y = 0; y < ny; y++)
for (x = 0; x < nx; x++) {
Rast3d_get_block(map, x * tileX, y * tileY, z * tileZ,
tileX, tileY, tileZ, data, typeIntern);
if (!Rast3d_write_tile
(map2, Rast3d_tile2tile_index(map2, x, y, z), data,
typeIntern))
Rast3d_fatal_error
("Rast3d_retileNocache: error in Rast3d_write_tile");
}
}
G_percent(1, 1, 1);
Rast3d_free_tiles(data);
Rast3d_close(map2);
}
void
Rast3d_retile(void *map, const char *nameOut, int tileX, int tileY, int tileZ)
{
void *map2;
double value;
int x, y, z, saveType;
int rows, cols, depths, typeIntern;
int xTile, yTile, zTile;
int xOffs, yOffs, zOffs, prev;
int tileXsave, tileYsave, tileZsave;
RASTER3D_Region region;
if (!Rast3d_tile_use_cache_map(map)) {
retileNocache(map, nameOut, tileX, tileY, tileZ);
return;
}
saveType = Rast3d_get_file_type();
Rast3d_set_file_type(Rast3d_file_type_map(map));
Rast3d_get_tile_dimension(&tileXsave, &tileYsave, &tileZsave);
Rast3d_set_tile_dimension(tileX, tileY, tileZ);
typeIntern = Rast3d_tile_type_map(map);
Rast3d_get_region_struct_map(map, ®ion);
map2 =
Rast3d_open_cell_new(nameOut, typeIntern, RASTER3D_USE_CACHE_DEFAULT, ®ion);
if (map2 == NULL)
Rast3d_fatal_error("Rast3d_retile: error in Rast3d_open_cell_new");
Rast3d_set_file_type(saveType);
Rast3d_set_tile_dimension(tileXsave, tileYsave, tileZsave);
Rast3d_coord2tile_coord(map2, 0, 0, 0,
&xTile, &yTile, &zTile, &xOffs, &yOffs, &zOffs);
prev = zTile;
x = 0;
y = 0;
Rast3d_get_coords_map(map, &rows, &cols, &depths);
for (z = 0; z < depths; z++) {
G_percent(z, depths, 1);
Rast3d_coord2tile_coord(map2, x, y, z, &xTile, &yTile, &zTile,
&xOffs, &yOffs, &zOffs);
if (zTile > prev) {
if (!Rast3d_flush_all_tiles(map2))
Rast3d_fatal_error("Rast3d_retile: error in Rast3d_flush_all_tiles");
prev++;
}
for (y = 0; y < rows; y++)
for (x = 0; x < cols; x++) {
Rast3d_get_value_region(map, x, y, z, &value, typeIntern);
if (!Rast3d_put_value(map2, x, y, z, &value, typeIntern))
Rast3d_fatal_error("Rast3d_retile: error in Rast3d_put_value");
}
}
G_percent(1, 1, 1);
if (!Rast3d_flush_all_tiles(map2))
Rast3d_fatal_error("Rast3d_retile: error in Rast3d_flush_all_tiles");
if (!Rast3d_close(map2))
Rast3d_fatal_error("Rast3d_retile: error in Rast3d_close");
}