Grid3.h

#ifndef GRID_3_H
#define GRID_3_H

#include "Vector3.h"
#include "MetaDataContainer.h"
#include <algorithm>
#include <string>
#include <stdint.h>
#include <vector>
using namespace std;

template<class T>
class Grid3{

protected:
    std::vector<T> data;
    int _width, _height, _depth;
    double _spacingX, _spacingY, _spacingZ;
    MetaDataContainer metadata;

public:
    Grid3(){}

    Grid3(int width, int height, int depth=1, double x_spacing=1.0, double y_spacing=1.0, double z_spacing=1.0);

    Grid3(Vec3i dimensions, Vec3d spacing=Vec3d(1.0,1.0,1.0));

    void Init(int width, int height, int depth=1,  double x_spacing=1.0, double y_spacing=1.0, double z_spacing=1.0); 

    int GetWidth(){ return _width;}

    int GetHeight(){ return _height;}

    int GetDepth(){ return _depth; }

    Vec3i GetDimensions(){return Vec3i(_width, _height, _depth);}

    double GetSpacingX(){return _spacingX;}

    double GetSpacingY(){return _spacingY;}

    double GetSpacingZ(){return _spacingZ;}

    Vec3d GetSpacing(){return Vec3d(_spacingX, _spacingY, _spacingZ);}

    void SetSpacing(double x_spacing, double y_spacing, double z_spacing);

    void SetSpacing(Vec3d spacing);

    Vec3d WorldToGrid(Vec3d p);

    Vec3d GridToWorld(Vec3d p);


    MetaDataContainer GetMetaData(){return metadata;}

    void SetMetaData(MetaDataContainer md){metadata=md;}

    int NumberOfElements();

    T* Data();

    std::string GetDataType();
     
    int Offset(int x, int y, int z);

    int Offset(Vec3i p){ return Offset(p.x,p.y,p.z); }

    Vec3i Coordinate(long index);

    bool IsInside(int x, int y, int z);

    bool IsInside(Vec3i p);

    T& operator[](long index);

    T& operator()(int x,int y, int z=0, bool Safe=false);

    T& operator()(Vec3i p, bool Safe=false);

    T LinearAt(double x, double y, double z);

    T LinearAt(Vec3d p);

    void Get6Neighbors(long index, vector<long>& n);
    
    void Get6Neighbors(Vec3i p, vector<Vec3i>& n);

    void Fill(T value);

    T GetMaxValue();

    T GetMinValue();

    /*
    Grid3<T> operator*(double x);
        
    Grid3<T> operator/(double x);

    Grid3<T> operator+(T x);

    Grid3<T> operator-(T x);*/

    void CopyDataFrom(Grid3<T>& src);

    void SwapContentsWith(Grid3<T>& v);

    bool CreateFromFile(std::string filename); 

    bool ReadVTK(std::string filename);

    bool WriteVTK(std::string filename);
    
    /*
    bool ReadAmira(std::string filename);

    bool WriteAmira(std::string filename);*/


    bool ReadMetaImage(std::string filename);

    bool WriteMetaImage(std::string filename);

    bool ReadNIFTI(std::string filename);

    bool WriteNIFTI(std::string filename);

    bool ReadDICOM(std::string filename);

};

#include "Grid3_i.h"


#endif