---

VolumetricData.h

Go to the documentation of this file.

/***************************************************************************
 *cr                                                                       
 *cr            (C) Copyright 1995-2008 The Board of Trustees of the           
 *cr                        University of Illinois                       
 *cr                         All Rights Reserved                        
 *cr                                                                   
 ***************************************************************************/

/***************************************************************************
 * RCS INFORMATION:
 *
 *      $RCSfile: VolumetricData.h,v $
 *      $Author: johns $        $Locker:  $             $State: Exp $
 *      $Revision: 1.22 $       $Date: 2008/04/30 00:04:15 $
 *
 ***************************************************************************
 * DESCRIPTION:
 *
 * Base class for storing volumetric data and associated gradient data
 *
 ***************************************************************************/

class VolumetricData {
public:
  char *name;              
  double origin[3];        
  double xaxis[3];         
  double yaxis[3];         
  double zaxis[3];         
  int xsize, ysize, zsize; 
  float *data;             
  float *gradient;         
  float datamin, datamax;  

  VolumetricData(const char *name, const float *origin, 
                 const float *xaxis, const float *yaxis, const float *zaxis,
                 int xs, int ys, int zs, float *dataptr);

  ~VolumetricData();

  void cell_lengths(float *xl, float *yl, float *zl) const;

  void cell_axes(float *xax, float *yax, float *zax) const;

  void cell_dirs(float *xax, float *yax, float *zax) const;

  void voxel_coord_from_cartesian_coord(const float *carcoord, float *voxcoord, int shiftflag) const;

  int voxel_index_from_coord(float xpos, float ypos, float zpos) const;

  inline float voxel_value(int x, int y, int z) const {
    return data[z*xsize*ysize + y*xsize + x];
  }

  float voxel_value_safe(int x, int y, int z) const;

  float voxel_value_interpolate(float xv, float yv, float zv) const;

  float voxel_value_from_coord(float xpos, float ypos, float zpos) const;
  float voxel_value_interpolate_from_coord(float xpos, float ypos, float zpos) const;


  void compute_volume_gradient(void);

  void voxel_gradient_fast(int x, int y, int z, float *grad) const {
    int index = (z*xsize*ysize + y*xsize + x) * 3;
    grad[0] = gradient[index    ];
    grad[1] = gradient[index + 1];
    grad[2] = gradient[index + 2];
  }

  void voxel_gradient_safe(int x, int y, int z, float *grad) const;

  void voxel_gradient_interpolate(const float *voxcoord, float *gradient) const;

  void voxel_gradient_from_coord(const float *coord, float *gradient) const;
  void voxel_gradient_interpolate_from_coord(const float *coord, float *gradient) const;

};


//
// Fast and loose accessor macros, don't use unless you have to
// 

#define VOXEL_GRADIENT_FAST(v, x, y, z, grad) \
  { int index = ((z)*v->xsize*v->ysize + (y)*v->xsize + (x)) * 3; \
    (grad)[0] = v->gradient[index    ]; \
    (grad)[1] = v->gradient[index + 1]; \
    (grad)[2] = v->gradient[index + 2]; \
  }

---