int Rast3d_get_standard3d_params(int *useTypeDefault, int *type,
int *useCompressionDefault, int *doCompression,
int *usePrecisionDefault, int *precision,
int *useDimensionDefault, int *tileX, int *tileY,
int *tileZ)
{
*useTypeDefault = *useCompressionDefault = 0;
*usePrecisionDefault = *useDimensionDefault = 0;
Rast3d_init_defaults();
if (strcmp(param->type->answer, "double") == 0)
*type = DCELL_TYPE;
else if (strcmp(param->type->answer, "float") == 0)
*type = FCELL_TYPE;
else {
*type = Rast3d_get_file_type();
*useTypeDefault = 1;
}
Rast3d_get_compression_mode(doCompression, precision);
if (strcmp(param->precision->answer, "default") != 0) {
if (strcmp(param->precision->answer, "max") == 0)
*precision = -1;
else if ((sscanf(param->precision->answer, "%d", precision) != 1) ||
(*precision < 0)) {
Rast3d_error(_("Rast3d_get_standard3d_params: precision value invalid"));
return 0;
}
}
else
*usePrecisionDefault = 1;
if (strcmp(param->compression->answer, "default") != 0) {
if (strcmp(param->compression->answer, "zip") == 0)
*doCompression = RASTER3D_COMPRESSION;
else
*doCompression = RASTER3D_NO_COMPRESSION;
} else {
*useCompressionDefault = 1;
}
Rast3d_get_tile_dimension(tileX, tileY, tileZ);
if (strcmp(param->dimension->answer, "default") != 0) {
if (sscanf(param->dimension->answer, "%dx%dx%d",
tileX, tileY, tileZ) != 3) {
Rast3d_error(_("Rast3d_get_standard3d_params: tile dimension value invalid"));
return 0;
}
}
else
*useDimensionDefault = 1;
Rast3d_free(param);
return 1;
}
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
modifyNull(char *name, d_Mask * maskRules, int changeNull, double newNullVal)
{
void *map, *mapOut;
RASTER3D_Region region;
int tileX, tileY, tileZ, x, y, z;
double value;
int doCompress, doLzw, doRle, precision;
int cacheSize;
cacheSize = Rast3d_cache_size_encode(RASTER3D_USE_CACHE_XY, 1);
if (NULL == G_find_raster3d(name, ""))
Rast3d_fatal_error(_("3D raster map <%s> not found"), name);
fprintf(stderr, "name %s Mapset %s \n", name, G_mapset());
map = Rast3d_open_cell_old(name, G_mapset(), RASTER3D_DEFAULT_WINDOW,
DCELL_TYPE, cacheSize);
if (map == NULL)
Rast3d_fatal_error(_("Unable to open 3D raster map <%s>"), name);
Rast3d_get_region_struct_map(map, ®ion);
Rast3d_get_tile_dimensions_map(map, &tileX, &tileY, &tileZ);
Rast3d_get_compression_mode(&doCompress, &doLzw, &doRle, &precision);
mapOut = Rast3d_open_new_param(name, DCELL_TYPE, RASTER3D_USE_CACHE_XY,
®ion, Rast3d_file_type_map(map),
doLzw, doRle, Rast3d_tile_precision_map(map), tileX,
tileY, tileZ);
if (mapOut == NULL)
Rast3d_fatal_error(_("modifyNull: error opening tmp file"));
Rast3d_min_unlocked(map, RASTER3D_USE_CACHE_X);
Rast3d_autolock_on(map);
Rast3d_unlock_all(map);
Rast3d_min_unlocked(mapOut, RASTER3D_USE_CACHE_X);
Rast3d_autolock_on(mapOut);
Rast3d_unlock_all(mapOut);
for (z = 0; z < region.depths; z++) {
if ((z % tileZ) == 0) {
Rast3d_unlock_all(map);
Rast3d_unlock_all(mapOut);
}
for (y = 0; y < region.rows; y++)
for (x = 0; x < region.cols; x++) {
value = Rast3d_get_double_region(map, x, y, z);
if (Rast3d_is_null_value_num(&value, DCELL_TYPE)) {
if (changeNull) {
value = newNullVal;
}
}
else if (Rast3d_mask_d_select((DCELL *) & value, maskRules)) {
Rast3d_set_null_value(&value, 1, DCELL_TYPE);
}
Rast3d_put_double(mapOut, x, y, z, value);
}
if ((z % tileZ) == 0) {
if (!Rast3d_flush_tiles_in_cube
(mapOut, 0, 0, MAX(0, z - tileZ), region.rows - 1,
region.cols - 1, z))
Rast3d_fatal_error(_("modifyNull: error flushing tiles in cube"));
}
}
if (!Rast3d_flush_all_tiles(mapOut))
Rast3d_fatal_error(_("modifyNull: error flushing all tiles"));
Rast3d_autolock_off(map);
Rast3d_unlock_all(map);
Rast3d_autolock_off(mapOut);
Rast3d_unlock_all(mapOut);
if (!Rast3d_close(map))
Rast3d_fatal_error(_("Unable to close raster map"));
if (!Rast3d_close(mapOut))
Rast3d_fatal_error(_("modifyNull: Unable to close tmp file"));
}