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");
}
void
Rast3d_coord2tile_index(RASTER3D_Map * map, int x, int y, int z, int *tileIndex,
int *offset)
{
int xTile, yTile, zTile, xOffs, yOffs, zOffs;
Rast3d_coord2tile_coord(map, x, y, z,
&xTile, &yTile, &zTile, &xOffs, &yOffs, &zOffs);
*tileIndex = Rast3d_tile2tile_index(map, xTile, yTile, zTile);
*offset = zOffs * map->tileXY + yOffs * map->tileX + xOffs;
}
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);
}
int
Rast3d_flush_tile_cube(RASTER3D_Map * map, int xMin, int yMin, int zMin, int xMax,
int yMax, int zMax)
{
int x, y, z;
if (!map->useCache)
Rast3d_fatal_error
("Rast3d_flush_tile_cube: function invalid in non-cache mode");
for (x = xMin; x <= xMax; x++)
for (y = yMin; y <= yMax; y++)
for (z = zMin; z <= zMax; z++)
if (!Rast3d_flush_tile(map, Rast3d_tile2tile_index(map, x, y, z))) {
Rast3d_error("Rast3d_flush_tile_cube: error in Rast3d_flush_tile");
return 0;
}
return 1;
}
