GRASS GIS 8 Programmer's Manual  8.2.2dev(2023)-3d2c704037
changetype.c
Go to the documentation of this file.
1 
2 #include <stdio.h>
3 #include <stdlib.h>
4 #include <sys/types.h>
5 #include <unistd.h>
6 #include <grass/raster3d.h>
7 
8 /*---------------------------------------------------------------------------*/
9 
10 
11 /*!
12  * \brief
13  *
14  * Makes a copy of <em>map</em> with name <em>nameOut</em> in which the
15  * cells are of type FCELL_TYPE if they are DCELL_TYPE in <em>map</em>,
16  * and in DCELL_TYPE otherwise.
17  * The source code can be found in <em>changetype.c</em>.
18  *
19  * \param map
20  * \param nameOut
21  * \return void
22  */
23 
24 void Rast3d_change_type(void *map, const char *nameOut)
25 {
26  void *map2;
27  int x, y, z, saveType;
28  void *data, *data2;
29  RASTER3D_Region region;
30  int tileSize;
31  int tileX, tileY, tileZ, typeIntern, typeIntern2;
32  int tileXsave, tileYsave, tileZsave, nx, ny, nz;
33 
34  saveType = Rast3d_get_file_type();
37  Rast3d_get_tile_dimension(&tileXsave, &tileYsave, &tileZsave);
38  Rast3d_get_tile_dimensions_map(map, &tileX, &tileY, &tileZ);
39  Rast3d_set_tile_dimension(tileX, tileY, tileZ);
40 
41  Rast3d_get_region_struct_map(map, &region);
42  map2 =
43  Rast3d_open_cell_new(nameOut, FCELL_TYPE, RASTER3D_USE_CACHE_DEFAULT, &region);
44 
45  if (map2 == NULL)
46  Rast3d_fatal_error("Rast3d_change_type: error in Rast3d_open_cell_new");
47 
48  Rast3d_set_file_type(saveType);
49  Rast3d_set_tile_dimension(tileXsave, tileYsave, tileZsave);
50 
51  data = Rast3d_alloc_tiles(map, 1);
52  if (data == NULL)
53  Rast3d_fatal_error("Rast3d_change_type: error in Rast3d_alloc_tiles");
54  data2 = Rast3d_alloc_tiles(map2, 1);
55  if (data2 == NULL)
56  Rast3d_fatal_error("Rast3d_change_type: error in Rast3d_alloc_tiles");
57 
58  Rast3d_get_nof_tiles_map(map2, &nx, &ny, &nz);
59  typeIntern = Rast3d_tile_type_map(map);
60  typeIntern2 = Rast3d_tile_type_map(map2);
61  tileSize = tileX * tileY * tileZ;
62 
63  for (z = 0; z < nz; z++)
64  for (y = 0; y < ny; y++)
65  for (x = 0; x < nx; x++) {
66  if (!Rast3d_read_tile(map, Rast3d_tile2tile_index(map, x, y, z), data,
67  typeIntern))
68  Rast3d_fatal_error("Rast3d_change_type: error in Rast3d_read_tile");
69  Rast3d_copy_values(data, 0, typeIntern, data2, 0, typeIntern2,
70  tileSize);
72  (map2, Rast3d_tile2tile_index(map2, x, y, z), data2,
73  typeIntern2))
74  Rast3d_fatal_error("Rast3d_change_type: error in Rast3d_write_tile");
75  }
76 
77  Rast3d_free_tiles(data);
78  Rast3d_free_tiles(data2);
79  if (!Rast3d_close(map2))
80  Rast3d_fatal_error("Rast3d_change_type: error in Rast3d_close");
81 }
#define RASTER3D_USE_CACHE_DEFAULT
Definition: raster3d.h:20
void Rast3d_get_region_struct_map(RASTER3D_Map *, RASTER3D_Region *)
Returns in region the region of map.
Definition: headerinfo.c:115
void Rast3d_set_tile_dimension(int, int, int)
set Tile Dimension
Definition: defaults.c:244
int Rast3d_file_type_map(RASTER3D_Map *)
Returns the type with which tiles of map are stored on file.
Definition: headerinfo.c:281
void Rast3d_get_tile_dimensions_map(RASTER3D_Map *, int *, int *, int *)
Returns the tile dimensions used for map.
Definition: headerinfo.c:142
int Rast3d_read_tile(RASTER3D_Map *, int, void *, int)
Reads tile with index tileIndex into the tile buffer. The cells are stored with type type which must ...
Definition: tileread.c:145
int Rast3d_close(RASTER3D_Map *)
Close 3D raster map files.
void * Rast3d_open_cell_new(const char *, int, int, RASTER3D_Region *)
Opens new g3d-file with name in the current mapset. Tiles are stored in memory with type which must b...
#define NULL
Definition: ccmath.h:32
#define x
void Rast3d_free_tiles(void *)
Is equivalent to Rast3d_free (tiles);
Definition: tilealloc.c:75
void Rast3d_fatal_error(const char *,...) __attribute__((format(printf
void Rast3d_change_type(void *map, const char *nameOut)
Makes a copy of map with name nameOut in which the cells are of type FCELL_TYPE if they are DCELL_TYP...
Definition: changetype.c:24
#define DCELL_TYPE
Definition: raster.h:13
void Rast3d_copy_values(const void *, int, int, void *, int, int, int)
Definition: raster3d/misc.c:53
int Rast3d_get_file_type(void)
get G3d file type
Definition: defaults.c:227
void * Rast3d_alloc_tiles(RASTER3D_Map *, int)
Is equivalent to Rast3d_alloc_tiles_type (map, nofTiles, Rast3d_file_type_map (map)).
Definition: tilealloc.c:50
int Rast3d_write_tile(RASTER3D_Map *, int, const void *, int)
Writes tile with index tileIndex to the file corresponding to map. It is assumed that the cells in ti...
Definition: tilewrite.c:128
void Rast3d_get_tile_dimension(int *, int *, int *)
get Tile Dimension
Definition: defaults.c:271
#define FCELL_TYPE
Definition: raster.h:12
void Rast3d_get_nof_tiles_map(RASTER3D_Map *, int *, int *, int *)
Returns the dimensions of the tile-cube used to tile the region of map. These numbers include partial...
Definition: headerinfo.c:51
int Rast3d_tile_type_map(RASTER3D_Map *)
Returns the type in which tiles of map are stored in memory.
Definition: headerinfo.c:161
void Rast3d_set_file_type(int)
set G3d file type
Definition: defaults.c:210
int Rast3d_tile2tile_index(RASTER3D_Map *, int, int, int)
Returns tile-index corresponding to tile-coordinates (xTile, yTile, zTile).
Definition: tilemath.c:52