mitkPhantom.h

/*=========================================================================

  Program:   3DMed
  Date:      $Date: 2014-02-25 18:30:00 +0800 $
  Version:   $Version: 4.6.0 $
  Copyright: MIPG, Institute of Automation, Chinese Academy of Sciences

=========================================================================*/


#ifndef __mitkPhantom_h
#define __mitkPhantom_h

#include "mitkObject.h"
#include <list>
#include "mitkReconstructionIncludes.h"

typedef enum {
        PELEM_INVALID,
        PELEM_ELLIPSOID,
        PELEM_CUBOID
} PhmElemType;

typedef double (* pMat2D)[3];

class mitkPhantomElement : public mitkObject
{
public:
    
    MITK_TYPE(mitkPhantomElement, mitkObject)
    virtual void PrintSelf(ostream &os);
    
    mitkPhantomElement(const char* type, const double cx, const double cy, const double cz, 
        const double a, const double b, const double c, const double rotx, const double roty, 
        const double rotz, const double atten); 
    ~mitkPhantomElement();

    double GetAtten() const {return m_atten; }
    void SetAtten(double atten) { m_atten = atten;}

    double GetCentX() const {return m_centX;}
    double GetCentY() const {return m_centY;}
    double GetCentZ() const {return m_centZ;}

    void SetCentX(double centX) { m_centX = centX;}
    void SetCentY(double centY) { m_centY = centY;}
    void SetCentZ(double centZ) { m_centZ = centZ;}

    double GetSemiAxisLengthX() const {return m_semiAxisLengthX;}
    double GetSemiAxisLengthY() const {return m_semiAxisLengthY;}
    double GetSemiAxisLengthZ() const {return m_semiAxisLengthZ;}

    void SetSemiAxisLengthX(double x) { m_semiAxisLengthX = x;}
    void SetSemiAxisLengthY(double y) { m_semiAxisLengthY = y;}
    void SetSemiAxisLengthZ(double z) { m_semiAxisLengthZ = z;}

    double GetRotX() const {return m_rotX;}
    double GetRotY() const {return m_rotY;}
    double GetRotZ() const {return m_rotZ;}

    void SetRotX(double rotX) { m_rotX = rotX; }
    void SetRotY(double rotY) { m_rotY = rotY; }
    void SetRotZ(double rotZ) { m_rotZ = rotZ; }

    pMat2D GetRotMat() {return m_rotMat;}
    static PhmElemType convertNameToType(const char* const typeName);

protected:
    PhmElemType m_type;     // phantom element type (ellipsoid, cuboid)
    double m_atten;         // x-ray attenuation coefficient
    double m_centX, m_centY, m_centZ;   // center of phantom element
    double m_semiAxisLengthX, m_semiAxisLengthY, m_semiAxisLengthZ; 
    double m_rotX, m_rotY, m_rotZ;      // rotation angle (in rad)
    double m_rotMat[3][3];  // rotation matrix

    static const char* const convertTypeToName(PhmElemType iType);
    void makeRotateMatrix();

private:
    mitkPhantomElement (const mitkPhantomElement& rhs);        // copy constructor
    mitkPhantomElement& operator= (const mitkPhantomElement&); // assignment operator

};

typedef enum {
        PHM_INVALID,
        PHM_SHEPP_LOGAN,
        PHM_MODIFIED_SHEPP_LOGAN,
} PhmType;

class MITK_RECONSTRUCTION_API mitkPhantom : public mitkObject
{
public:
    MITK_TYPE(mitkPhantom, mitkObject)
    virtual void PrintSelf(ostream &os);

    mitkPhantom();
    mitkPhantom(const char* const phmName);
    
    ~mitkPhantom();
    bool createFromPhantom (const char* const phmName);
    bool createFromPhantom (PhmType phmtype);   
    bool createFromFile (const char* const fname);
    bool fileWrite (const char* const fname);
    
    void addPElem (const mitkPhantomElement& pelem);
    void addPElem (const char* type, const double cx, const double cy, const double cz, 
        const double a, const double b, const double c, const double rotx, const double roty, const double rotz, const double atten);

    PhmType convertPhmNameToType(const char* const phmName);

    std::list<mitkPhantomElement*> GetPElemList() {return m_listPElem;}
    int GetPElemNum() {return m_nPElem;}

protected:
//    PhmType m_type;                           // shepp-logan or user defined
//    std::string m_name;

    int m_nPElem;                           // number of pelems in phantom
    mutable std::list<mitkPhantomElement*> m_listPElem;      // pelem lists 

    void init();
    void addSheppLogan();
    void addModifiedSheppLogan();
    
    static const char* const convertPhmTypeToName(PhmType phmType);

private:
    mitkPhantom (const mitkPhantom& rhs);        // copy constructor
    mitkPhantom& operator= (const mitkPhantom&); // assignment operator
};

typedef std::list<mitkPhantomElement*>::iterator mitkPElemIterator;
typedef std::list<mitkPhantomElement*>::const_iterator mitkPElemConstIterator;



//#define DEFINED_mitkPhantom



#endif

---