1 #ifndef VOLUME_H_

2 #define VOLUME_H_

3

4 #include <vector>

5 #include "image.h"

6

7 class Volume {

8 public:

9 virtual ~Volume();

10

11 // returns 0 if the volume is continuously defined

12 virtual int num_samples(int dim) const;

13

14 // central differences offset will be delta / num_samples for discrete volumes

15 // and delta for continuous volumes

16 virtual void normalf(float *norm, float x, float y, float z, float delta = 1.0);

17 virtual void normali(float *norm, int x, int y, int z);

18

19 virtual float valuef(float x, float y, float z) const = 0;

20 virtual float valuei(int x, int y, int z) const = 0;

21 };

22

23 class VoxelVolume : public Volume {

24 protected:

25 int size[3];

26 std::vector<Image> slices;

27

28 public:

29 VoxelVolume();

30

31 bool load(const char *fname);

32

33 int num_samples(int dim) const;

34

35 float valuef(float x, float y, float z) const;

36 float valuei(int x, int y, int z) const;

37 };

38

39 #endif // VOLUME_H_