32 #define MODE_PRELOAD 3 34 #define MODE_DEFAULT 0 36 #define STATUS_READY 0 53 static int Numfiles = 0;
55 static int Cur_id =
LUCKY;
69 static int init_volfiles(
void)
97 static int check_num_volfiles(
void)
99 if (Numfiles < Cur_max) {
121 for (i = 0; i < Numfiles; i++) {
122 if (Data[i]->data_id ==
id) {
144 start = begin ? 0 : i + 1;
146 for (i = start; i < Numfiles; i++) {
147 if (!strcmp(Data[i]->file_name, name)) {
148 if (Data[i]->file_type == type) {
149 return (Data[i]->data_id);
171 for (i = 0; i < Numfiles; i++) {
172 if (Data[i]->data_id ==
id) {
272 static int first = 1;
279 if (0 > init_volfiles()) {
286 if (0 <= (
id =
find_datah(name, file_type, 1))) {
287 for (i = 0; i < Numfiles; i++) {
288 if (Data[i]->data_id ==
id) {
297 if (0 > check_num_volfiles()) {
309 new = Data[Numfiles];
313 new->data_id = Cur_id++;
315 new->file_name =
G_store(name);
316 new->file_type = file_type;
321 new->data_type = data_type;
329 return (new->data_id);
366 int i, j, found = -1;
369 G_debug(5,
"gvl_file_free_datah(): id=%d",
id);
371 for (i = 0; i < Numfiles; i++) {
372 if (Data[i]->data_id ==
id) {
376 if (fvf->
count > 1) {
388 for (j = i; j < (Numfiles - 1); j++) {
389 Data[j] = Data[j + 1];
427 G_warning(
_(
"3D raster map <%s> not found"), filename);
436 G_warning(
_(
"Unable to open 3D raster map <%s>"), filename);
442 G_warning(
_(
"Unable to read range of 3D raster map <%s>"), filename);
470 G_warning(
_(
"Unable to close 3D raster map <%s>"),
528 for (x = 0; x <
Cols; x++) {
529 for (y = 0; y <
Rows; y++) {
530 ((
float *)data)[x + y *
Cols] =
539 for (x = 0; x <
Cols; x++) {
540 for (y = 0; y <
Rows; y++) {
541 ((
double *)data)[x + y *
Cols] =
573 for (x = 0; x <
Cols; x++) {
574 for (y = 0; y <
Rows; y++) {
575 for (z = 0; z <
Depths; z++) {
576 ((
float *)data)[x + y * Cols + z * Rows *
Cols] =
586 for (x = 0; x <
Cols; x++) {
587 for (y = 0; y <
Rows; y++) {
588 for (z = 0; z <
Depths; z++) {
589 ((
double *)data)[x + y * Cols + z * Rows *
Cols] =
657 *((
float *)value) = ((
float *)data)[offset];
662 *((
double *)value) = ((
double *)data)[offset];
814 slice_data *sd = (slice_data *) vf->
buff;
819 for (i = 0; i < sd->num; i++) {
828 for (i = 0; i < sd->num; i++) {
853 slice_data *sd = (slice_data *) vf->
buff;
855 for (i = 0; i < sd->num; i++) {
875 slice_data *sd = (slice_data *) vf->
buff;
901 slice_data *sd = (slice_data *) vf->
buff;
905 for (i = 0; i < sd->num - 1; i++) {
906 sd->slice[i] = sd->slice[i + 1];
908 sd->slice[sd->num - 1] = tmp;
911 read_slice(vf, sd->num, sd->crnt + 1 + (sd->num - sd->base));
929 slice_data *sd = (slice_data *) vf->
buff;
932 if ((z >= sd->crnt - (sd->base - 1)) &&
933 (z <= sd->crnt + sd->num - sd->base)) {
940 else if (z == sd->crnt - (sd->base - 1) + 1) {
994 sd = (slice_data *) vf->
buff;
999 for (i = 0; i < (sd->num - sd->base + 1); i++)
1123 if (vf->
mode == mode)
1136 sd = (slice_data *) vf->
buff;
1176 sd = (slice_data *) vf->
buff;
int is_null_g3d_value(IFLAG type, void *value)
Check for null value.
#define RASTER3D_USE_CACHE_DEFAULT
void * GVL_get_window()
Get window.
void void void void G_fatal_error(const char *,...) __attribute__((format(printf
void * Rast3d_open_cell_old(const char *, const char *, RASTER3D_Region *, int, int)
Opens existing g3d-file name in mapset. Tiles are stored in memory with type which must be any of FCE...
int close_g3d_file(void *)
Close g3d file.
int Rast3d_file_type_map(RASTER3D_Map *)
Returns the type with which tiles of map are stored on file.
int alloc_vol_buff(geovol_file *vf)
Allocate buffer memory for full mode reading.
geovol_file * gvl_file_get_volfile(int id)
Get geovol_file structure for given handle.
#define RASTER3D_DEFAULT_WINDOW
int Rast3d_range_load(RASTER3D_Map *)
Loads the range into the range structure of map.
void G_free(void *)
Free allocated memory.
int Rast3d_close(RASTER3D_Map *)
Close 3D raster map files.
int gvl_file_get_data_type(geovol_file *vf)
Get data type for given handle.
const char * G_find_raster3d(const char *, const char *)
Search for a 3D raster map in current search path or in a specified mapset.
int read_vol(geovol_file *vf)
Read all values from volume file.
char * gvl_file_get_name(int id)
Get file name for given handle.
int gvl_file_get_file_type(geovol_file *vf)
Get file type for given handle.
double Rast3d_get_double(RASTER3D_Map *, int, int, int)
Is equivalent to Rast3d_get_value (map, x, y, z, &value, DCELL_TYPE); return value.
int get_buff_value(IFLAG type, void *data, int offset, void *value)
Get value from buffer.
int close_volfile(void *map, IFLAG type)
Close volume file.
int read_g3d_value(IFLAG type, void *map, int x, int y, int z, void *value)
Eead value from g3d file.
int gvl_file_set_slices_param(geovol_file *vf, int n, int b)
Set parameters for slice reading.
int get_slice_value(geovol_file *vf, int x, int y, int z, void *value)
Get value from slice buffer.
void Rast3d_range_min_max(RASTER3D_Map *, double *, double *)
Returns in min and max the minimum and maximum values of the range.
int free_volfile_buffs(geovol_file *vf)
Free allocated buffers.
int gvl_file_newh(const char *name, IFLAG file_type)
Get handle for given file name and type.
int alloc_slice_buff(geovol_file *vf)
Allocate buffer for slice mode reading.
void shift_slices(geovol_file *vf)
Read new slice into buffer.
int gvl_file_start_read(geovol_file *vf)
Start read - allocate memory buffer a read first data into buffer.
int get_direct_value(geovol_file *vf, int x, int y, int z, void *value)
Read value direct from volume file.
int gvl_file_free_datah(int id)
Free geovol_file structure for given handle.
int get_vol_value(geovol_file *vf, int x, int y, int z, void *value)
Get value from volume buffer.
void * open_g3d_file(const char *, IFLAG *, double *, double *)
Open 3d raster file.
int gvl_file_end_read(geovol_file *vf)
End read - free buffer memory.
#define RASTER3D_TILE_SAME_AS_FILE
void G_warning(const char *,...) __attribute__((format(printf
int read_g3d_vol(IFLAG type, void *map, void *data)
Read all values from g3d file.
int gvl_file_is_null_value(geovol_file *vf, void *value)
Check for null value.
void gvl_file_get_min_max(geovol_file *vf, double *min, double *max)
Get minimum and maximum value in volume file.
char * G_store(const char *)
Copy string to allocated memory.
int Rast3d_is_null_value_num(const void *, int)
int free_slice_buff(geovol_file *vf)
Free buffer for slice mode reading.
void * open_volfile(const char *name, IFLAG file_type, IFLAG *data_type, double *min, double *max)
Open 3d raster file.
#define VOL_FTYPE_RASTER3D
int free_vol_buff(geovol_file *vf)
Free memory buffer memory.
int find_datah(const char *name, IFLAG type, int begin)
Find file with name and type in geovol_file array an return handle.
int gvl_file_set_mode(geovol_file *vf, IFLAG mode)
Set read mode.
int read_g3d_slice(IFLAG type, void *map, int level, void *data)
Read slice of values at level from g3d file.
float Rast3d_get_float(RASTER3D_Map *, int, int, int)
Is equivalent to Rast3d_get_value (map, x, y, z, &value, FCELL_TYPE); return value.
int G_debug(int, const char *,...) __attribute__((format(printf
int gvl_file_get_value(geovol_file *vf, int x, int y, int z, void *value)
Get value for volume file at x, y, z.
int read_slice(geovol_file *vf, int s, int l)
Read slice of values at level from volume file.