24#include <grass/ogsf.h>
25#include <grass/raster3d.h>
26#include <grass/glocale.h>
47 void *slice[MAX_VOL_SLICES];
50static geovol_file *Data[MAX_VOL_FILES];
51static geovol_file Df[MAX_VOL_FILES];
53static int Numfiles = 0;
55static int Cur_id =
LUCKY;
61void *
open_g3d_file(
const char *, IFLAG *,
double *,
double *);
69static int init_volfiles(
void)
74 for (i = 0; i < MAX_VOL_FILES; i++) {
79 Cur_max = MAX_VOL_FILES;
97static 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);
169 static char retstr[GPATH_MAX];
171 for (i = 0; i < Numfiles; i++) {
172 if (Data[i]->data_id ==
id) {
174 strcpy(retstr, fvf->file_name);
192 return (vf->file_type);
204 return (vf->data_type);
235 if (file_type == VOL_FTYPE_RASTER3D) {
253 if (type == VOL_FTYPE_RASTER3D) {
272 static int first = 1;
279 if (0 > init_volfiles()) {
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++;
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) {
384 fvf->file_name =
NULL;
388 for (j = i; j < (Numfiles - 1); j++) {
389 Data[j] = Data[j + 1];
426 G_warning(_(
"3D raster map <%s> not found"), filename);
431 map = Rast3d_open_cell_old(filename, mapset, RASTER3D_DEFAULT_WINDOW,
432 RASTER3D_TILE_SAME_AS_FILE,
433 RASTER3D_USE_CACHE_DEFAULT);
435 G_warning(_(
"Unable to open 3D raster map <%s>"), filename);
440 if (!Rast3d_range_load(map)) {
441 G_warning(_(
"Unable to read range of 3D raster map <%s>"), filename);
445 Rast3d_range_min_max(map,
min,
max);
448 itype = Rast3d_file_type_map(map);
449 if (itype == FCELL_TYPE)
450 *type = VOL_DTYPE_FLOAT;
451 if (itype == DCELL_TYPE)
452 *type = VOL_DTYPE_DOUBLE;
468 if (Rast3d_close((RASTER3D_Map *)map) != 1) {
469 G_warning(_(
"Unable to close 3D raster map <%s>"),
470 ((RASTER3D_Map *)map)->fileName);
492 case (VOL_DTYPE_FLOAT):
493 *((
float *)value) = Rast3d_get_float(map,
x, y, z);
497 case (VOL_DTYPE_DOUBLE):
498 *((
double *)value) = Rast3d_get_double(map,
x, y, z);
526 case (VOL_DTYPE_FLOAT):
528 for (y = 0; y <
Rows; y++) {
529 ((
float *)data)[
x + y *
Cols] =
530 Rast3d_get_float(map,
x, y, level);
537 case (VOL_DTYPE_DOUBLE):
539 for (y = 0; y <
Rows; y++) {
540 ((
double *)data)[
x + y *
Cols] =
541 Rast3d_get_double(map,
x, y, level);
571 case (VOL_DTYPE_FLOAT):
573 for (y = 0; y <
Rows; y++) {
574 for (z = 0; z <
Depths; z++) {
576 Rast3d_get_float(map,
x, y, z);
584 case (VOL_DTYPE_DOUBLE):
586 for (y = 0; y <
Rows; y++) {
587 for (z = 0; z <
Depths; z++) {
589 Rast3d_get_double(map,
x, y, z);
618 case (VOL_DTYPE_FLOAT):
619 return Rast3d_is_null_value_num(value, FCELL_TYPE);
623 case (VOL_DTYPE_DOUBLE):
624 return Rast3d_is_null_value_num(value, DCELL_TYPE);
655 case (VOL_DTYPE_FLOAT):
656 *((
float *)value) = ((
float *)data)[offset];
660 case (VOL_DTYPE_DOUBLE):
661 *((
double *)value) = ((
double *)data)[offset];
689 switch (vf->file_type) {
691 case (VOL_FTYPE_RASTER3D):
718 switch (vf->data_type) {
720 case (VOL_DTYPE_FLOAT):
721 if ((vf->buff = (
float *)G_malloc(
sizeof(
float) *
Cols *
Rows *
727 case (VOL_DTYPE_DOUBLE):
728 if ((vf->buff = (
double *)G_malloc(
sizeof(
double) *
Cols *
Rows *
765 switch (vf->file_type) {
767 case (VOL_FTYPE_RASTER3D):
812 slice_data *sd = (slice_data *)vf->buff;
814 switch (vf->data_type) {
816 case (VOL_DTYPE_FLOAT):
817 for (i = 0; i < sd->num; i++) {
819 (
float *)G_malloc(
sizeof(
float) *
Cols *
Rows)) ==
NULL)
825 case (VOL_DTYPE_DOUBLE):
826 for (i = 0; i < sd->num; i++) {
828 (
double *)G_malloc(
sizeof(
double) *
Cols *
Rows)) ==
NULL)
851 slice_data *sd = (slice_data *)vf->buff;
853 for (i = 0; i < sd->num; i++) {
873 slice_data *sd = (slice_data *)vf->buff;
875 switch (vf->file_type) {
877 case (VOL_FTYPE_RASTER3D):
899 slice_data *sd = (slice_data *)vf->buff;
903 for (i = 0; i < sd->num - 1; i++) {
904 sd->slice[i] = sd->slice[i + 1];
906 sd->slice[sd->num - 1] = tmp;
909 read_slice(vf, sd->num, sd->crnt + 1 + (sd->num - sd->base));
927 slice_data *sd = (slice_data *)vf->buff;
930 if ((z >= sd->crnt - (sd->base - 1)) &&
931 (z <= sd->crnt + sd->num - sd->base)) {
938 else if (z == sd->crnt - (sd->base - 1) + 1) {
992 sd = (slice_data *)vf->buff;
997 for (i = 0; i < (sd->num - sd->base + 1); i++)
1087 switch (vf->file_type) {
1089 case (VOL_FTYPE_RASTER3D):
1121 if (vf->mode == mode)
1131 if ((vf->buff = (slice_data *)G_malloc(
sizeof(slice_data))) ==
NULL)
1134 sd = (slice_data *)vf->buff;
1174 sd = (slice_data *)vf->buff;
void G_free(void *buf)
Free allocated memory.
int G_debug(int level, const char *msg,...)
Print debugging message.
const char * G_find_raster3d(const char *name, const char *mapset)
Search for a 3D raster map in current search path or in a specified mapset.
void G_fatal_error(const char *msg,...)
Print a fatal error message to stderr.
void G_warning(const char *msg,...)
Print a warning message to stderr.
void * GVL_get_window(void)
Get window.
int close_g3d_file(void *)
Close g3d 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_free_datah(int id)
Free geovol_file structure for given handle.
int gvl_file_get_data_type(geovol_file *vf)
Get data type for given handle.
int close_volfile(void *map, IFLAG type)
Close volume file.
int gvl_file_end_read(geovol_file *vf)
End read - free buffer memory.
void shift_slices(geovol_file *vf)
Read new slice into buffer.
int get_vol_value(geovol_file *vf, int x, int y, int z, void *value)
Get value from volume buffer.
int get_buff_value(IFLAG type, void *data, int offset, void *value)
Get value from buffer.
int gvl_file_newh(const char *name, IFLAG file_type)
Get handle for given file name and type.
int read_vol(geovol_file *vf)
Read all values from volume file.
void gvl_file_get_min_max(geovol_file *vf, double *min, double *max)
Get minimum and maximum value in volume file.
int get_direct_value(geovol_file *vf, int x, int y, int z, void *value)
Read value direct from volume file.
int gvl_file_set_mode(geovol_file *vf, IFLAG mode)
Set read mode.
int get_slice_value(geovol_file *vf, int x, int y, int z, void *value)
Get value from slice buffer.
int free_slice_buff(geovol_file *vf)
Free buffer for slice mode reading.
int alloc_slice_buff(geovol_file *vf)
Allocate 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.
int gvl_file_set_slices_param(geovol_file *vf, int n, int b)
Set parameters for slice reading.
int read_g3d_slice(IFLAG type, void *map, int level, void *data)
Read slice of values at level from g3d file.
int free_vol_buff(geovol_file *vf)
Free memory buffer memory.
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 gvl_file_start_read(geovol_file *vf)
Start read - allocate memory buffer a read first data into buffer.
char * gvl_file_get_name(int id)
Get file name for given handle.
void * open_g3d_file(const char *, IFLAG *, double *, double *)
Open 3d raster file.
int find_datah(const char *name, IFLAG type, int begin)
Find file with name and type in geovol_file array an return handle.
int is_null_g3d_value(IFLAG type, void *value)
Check for null value.
int gvl_file_get_file_type(geovol_file *vf)
Get file type for given handle.
int free_volfile_buffs(geovol_file *vf)
Free allocated buffers.
int gvl_file_is_null_value(geovol_file *vf, void *value)
Check for null value.
geovol_file * gvl_file_get_volfile(int id)
Get geovol_file structure for given handle.
int read_slice(geovol_file *vf, int s, int l)
Read slice of values at level from volume file.
int read_g3d_vol(IFLAG type, void *map, void *data)
Read all values from g3d file.
int alloc_vol_buff(geovol_file *vf)
Allocate buffer memory for full mode reading.
char * G_store(const char *s)
Copy string to allocated memory.